משוואת קושי-אוילר

מתוך המכלול, האנציקלופדיה היהודית
קפיצה לניווט קפיצה לחיפוש

משוואת קושי-אוילר (לעיתים משוואת אוילר) היא משוואה דיפרנציאלית רגילה, אשר לה דרך פתרון ייחודית שקשורה לפתרון משוואות דיפרנציאליות לינאריות עם מקדמים קבועים. קרויה על שמותיהם של המתמטיקאים אוגוסטן לואי קושי ולאונרד אוילר.

הגדרה פורמלית

משוואות אוילר ההומוגנית נתונה בצורה הבאה: , כאשר המקדמים הם מספרים ממשיים, והמקדם של הנגזרת מהסדר הגבוה ביותר שונה מאפס.

זוהי משוואה לינארית הומוגנית, ולכן יש לה n פתרונות בלתי תלויים המהווים מרחב וקטורי.

משוואות אוילר הלא הומוגנית היא מהצורה: , כאשר f היא פונקציה כלשהי שלא שווה זהותית לאפס.

לעיתים דנים במשוואה המנורמלת, כלומר כאשר המקדם של הנגזרת מהסדר הגבוה ביותר הוא 1. תמיד ניתן להגיע לצורה זו על ידי חלוקה במקדם, ולכן מעתה נתייחס למשוואה כבצורה:

דוגמאות

נציג מספר משוואות אוילר.

1)

אפשר להיווכח בכך ש- וגם הם פתרונות בלתי תלויים למשוואה, ולכן צירוף לינארי שלהם הוא הפתרון הכללי.

2) הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle { x }^{ 2 }{ y }''-x{ y }'+y=0}

כאן הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle y(x)=x} הוא פתרון למשוואה. גם הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle y(x)=xln(x)} הוא פתרון למשוואה. היות שהם בלתי תלויים לינארית, שניהם מהווים בסיס ולכן הפתרון הכללי הוא צירוף לינארי של שני אלו.

3) הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle { x }^{ 2 }{ y }''+xy'+y=0}

מכאן שהפתרונות הבלתי תלויים למשוואה הם הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle { y }_{ 1 }(x)=sin(ln(x)), { y }_{ 2 }(x)=cos(ln(x))} .

בהמשך נראה את השיטה לפתרון המשוואה.

מעבר למשוואה עם מקדמים קבועים

אחד הדברים החשובים בנוגע למשוואת אוילר, הוא שהיא שקולה למשוואה לינארית עם מקדמים קבועים. כלומר, ניתן לקחת כל משוואה כזו ולהפוך אותה למשוואה לינארית עם מקדמים קבועים. עם זאת, לא זוהי השיטה לפתור את המשוואה בפועל - שיטה זו תוצג בהמשך.

לצורך הפשטות, נניח שלפנינו משוואת אוילר הבאה: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle { x }^{ 2 }y''+{a}_{1}xy'+{a}_{0}{y}=0}

אם נגביל את הדיון לx>0, נשתמש בהצבה (עבור x<0 אפשר בדומה להציב ).

מתוך ההצבה נקבל הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle t=ln(x)} .

הנגזרות נתונות, לפי כלל השרשרת, על ידי:

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle y'=\frac { dy }{ dx } =\frac { dy }{ dt } \frac { dt }{ dx } =\frac { 1 }{ x } \frac { dy }{ dt } ={ e }^{ -t }\frac { dy }{ dt } }

הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle y''=\frac { d }{ dx } \frac { dy }{ dx } =\frac { d }{ dt } ({ e }^{ -t }\frac { dy }{ dt } )\frac { dt }{ dx } =(-{ e }^{ -t }\frac { dy }{ dt } +{ e }^{ -t }\frac { { d }^{ 2 }y }{ d{ t }^{ 2 } } ){ e }^{ -t }}

אם נציב חזרה על המשוואה המקורית נקבל: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \frac { { d }^{ 2 }y }{ d{ t }^{ 2 } } +({ a }_{ 1 }-1)\frac { dy }{ dt } +{ a }_{ 0 }y=0}

זו משוואה לינארית עם מקדמים קבועים.

גם המעבר בכיוון הנגדי אפשרי, על ידי ההצבה ההפוכה.

(במקרה הכללי פועלים באופן דומה - מוצאים את כל הנגזרות לפי כלל השרשרת, מציבים ומגיעים לאותה התוצאה למקרה של משוואה מסדר n).

שיטת פתרון המשוואה ההומוגנית

ראינו לעיל שמשוואת אוילר שקולה למשוואה לינארית עם מקדמים קבועים. כידוע, הפתרון של משוואה לינארית עם מקדמים קבועים הוא אקספוננט. לכן, גם פתרון המשוואה הנ"ל לפי t הוא אקספוננט: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle y(t)={ e }^{ rt }} . אם נחזור למשתנה x, נקבל: .

לכן, הפתרון של משוואת אוילר נתון על ידי חזקה של x.

לאחר הצגת עובדה זו, אין שימוש בפועל בשיטה של מעבר למשוואה לינארית. אלא, מנחשים פתרון מהצורה הנ"ל, ומוצאים את r שיקיים אותו.

אם שוב נעסוק במקרה של משוואה מסדר 2, נקבל:

נניח שניחשנו את הפתרון הבא: . אז נציב:

כלומר .

זה גורר שבהכרח מתאפס הפולינום לפי r: הפענוח נכשל (MathML עם גיבוי SVG או PNG (מומלץ לדפדפנים מודרניים ולכלי נגישות): תגובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle r(r-1)+{ a }_{ 1 }r+{ a }_{ 0 }=0}

פולינום זה נקרא המשוואה האינדיציאלית של המד"ר.

כעת נראה דוגמה לפתרון לפי השיטה הזו. למשל, במשוואה שהובאה בדוגמה 1 לעיל: , המשוואה האינדיציאלית היא: ששורשיה הם . לכן הפתרונות הם וגם , כנאמר לעיל.

השורשים של המשוואה האינדיציאלית יכולים להתחלק לכמה סוגים:

  • שורש ממשי עם ריבוי - לכל שורש המופיע בריבוי במשוואה האינדיציאלית, הפתרונות הם: .
  • שורש מרוכב - אם הוא שורש, אז גם הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle a-bi} הוא שורש, והפתרון עבורם נתון על ידי: הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle {x}^{a}\cos(b\cdot ln(x)),\quad {x}^{a}\sin(b\cdot \ln(x))}
  • שורש מרוכב עם ריבוי - אם הוא שורש מרוכב עם ריבוי m, אזי הפתרונות הם: הפענוח נכשל (שגיאת המרה. השרת ("https://en.wikipedia.org/api/rest_") השיב: "Cannot get mml. Server problem."): {\displaystyle {x}^{a}\cos(b\cdot ln(x)),{x}^{a}\sin(b\cdot \ln(x)),...,{\ln }^{m-1}(x)\cdot {x}^{a}\cos(b\cdot \ln(x)),{ln}^{m-1}(x)\cdot {x}^{a}\sin(b\cdot \ln(x))}

מסקנות אלו נכונות בגלל היותן נכונות ל משוואות דיפרנציאליות לינאריות במקדמים קבועים.

שיטת פתרון המשוואה הלא הומוגנית

כמו בכל משוואה לינארית, גם כאן הפתרון למשוואה הלא הומוגנית הוא פתרון ההומוגנית ועוד פתרון פרטי של הלא הומוגנית.

לכן, פתרון בעיה זו מתחלק לשני שלבים. ראשית, יש למצוא את הפתרון למשוואה ההומוגנית, בשיטה שתוארה לעיל. שנית, יש למצוא פתרון פרטי אחד למשוואה הלא הומוגנית. במציאת פתרון כזה נעסוק כעת.

תהי המשוואה:

ונניח שכבר מצאנו את הפתרון להומוגנית. מה שנכון עבור משוואות דיפרנציאליות לינאריות, נכון אנלוגית גם כאן.

בצורה הכללית ביותר, אם

ואם הוא שורש של המשוואה האינדיציאלית מסדר m, אז הפתרונות הם:

ראו גם