משוואה דיפרנציאלית ליניארית
במתמטיקה, משוואה דיפרנציאלית ליניארית היא משוואה דיפרנציאלית רגילה בפונקציה הנעלמת , שאפשר להציג בצורה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y^{(n)}+p_{n-1}(x)y^{(n-1)}+\dots+p_1(x)y'+p_0(x)y=g(x)} , כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ p_{n-1},\dots,p_0,g} הן פונקציות של המשתנה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ x} בלבד. בניגוד למרבית טיפוסי המשוואות הדיפרנציאליות, למשוואות הליניאריות יש תאוריה מפותחת וטכניקות פתרון שיטתיות, והן מופיעות בתחומי מדע רבים.
רישום בצורה אופרטורית
הסדר של משוואה דיפרנציאלית הוא הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ n} הגבוה ביותר שעבורו מופיעה במשוואה הנגזרת הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y^{(n)}} (אם קיים כזה).
משוואה ליניארית מסדר n אפשר לרשום בצורה מקוצרת אם מגדירים אופרטור כך: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ L=\frac{d^n}{dx^n}+p_{n-1}\frac{d^{n-1}}{dx^{n-1}}+\dots+p_0} , ואז המשוואה נרשמת כך: . המשוואה נקראת "ליניארית" שכן אופרטור זה הוא ליניארי: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ L\left[\lambda_1 y_1+\lambda_2 y_2\right]=\lambda_1 L\left[y_1\right]+\lambda_2 L\left[y_2\right]} .
אם הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ g(x)\equiv 0} המשוואה נקראת הומוגנית.
פתרון משוואות ליניאריות
באופן כללי, פונקציות רבות יכולות להיות פתרון של אותה משוואה. לפונקציה מסוימת שהיא פתרון של המשוואה נקרא "פתרון פרטי" של המשוואה. עבור משוואה הומוגנית הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ L\left[y\right]=0} סכום של כל שני פתרונות הוא פתרון וכפל בסקלר של פתרון הוא פתרון. לכן אוסף הפתרונות של משוואה ליניארית הומוגנית הוא מרחב וקטורי ויש לו בסיס, כלומר קבוצת פונקציות הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1,\dots,y_n} כך שכל פתרון של המשוואה ההומוגנית יכול להיכתב כצירוף ליניארי שלהן: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y=c_1 y_1+\dots+c_n y_n} . עבור בסיס של מרחב הפתרונות, נקרא לצירוף הליניארי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ c_1 y_1+\dots+c_n y_n} , כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ c_1,\dots,c_n} הם קבועים, "פתרון כללי של המשוואה ההומוגנית".
למשוואה ליניארית התכונה שהפרש של כל שני פתרונות של המשוואה הוא פתרון של המשוואה ההומוגנית הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ L\left[y\right]=0} . אכן, אם הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ L\left[y_1\right]=L\left[y_2\right]=g(x)} אז הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ L\left[y_1-y_2\right]=L\left[y_1\right]-L\left[y_2\right]=g(x)-g(x)=0} .
מכך נובעת תכונה חשובה של משוואות ליניאריות: כל פתרון של משוואה ליניארית ניתן לכתיבה כסכום של פתרון פרטי של המשוואה הליניארית, ופתרון כללי של המשוואה ההומוגנית המתאימה לה.
תכונה זו ניתן לקבל בצורה ישירה: אם הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y} הוא פתרון כלשהו של המשוואה, ו- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_p} הוא הפתרון הפרטי שאנו משתמשים בו, אז ההפרש הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y-y_p} , כפי שראינו, הוא פתרון של המשוואה ההומוגנית, וכל פתרון של המשוואה ההומוגנית ניתן לביטוי באמצעות הפתרון הכללי של המשוואה ההומוגנית.
מכך נובע שכדי לפתור בצורה כללית משוואה ליניארית לא הומוגנית יש לעשות שני דברים:
- לפתור את המשוואה ההומוגנית המתאימה לה.
- למצוא פתרון פרטי אחד למשוואה הלא הומוגנית.
קיימת שיטה בשם "וריאציית הפרמטר", המאפשרת למצוא בצורה שיטתית פתרון פרטי למשוואה הלא הומוגנית בהינתן הפתרון הכללי של המשוואה ההומוגנית. עם זאת, השיטה עלולה לכלול עבודה טכנית רבה, ונדרשת בה אינטגרציה שבה לא בהכרח ניתן לבטא את התוצאה בצורה מפורשת. שיטה אחרת, הנקראת "שיטת המקדמים הלא ידועים", מבוססת על ניחוש מושכל של צורת הפתרון הפרטי והשימוש בה הוא פשוט ונוח, אך היא טובה רק למקרים מסוימים.
באופן כללי אין פתרון שיטתי ופשוט למשוואה ליניארית הומוגנית. עם זאת, אם ידוע פתרון אחד של המשוואה, ניתן לקבל ממנו פתרון נוסף, שאינו תלוי בו. שיטה זו ידועה בשם "שיטת ד'אלמבר להורדת סדר המשוואה". כמו כן, במקרה הפרטי המיוחד שבו המקדמים הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ p_{n-1},\dots,p_0} הם כולם קבועים קיים פתרון שיטתי המתבסס על מה שמכונה "המשוואה האופיינית" (או: "הפולינום האופייני") של המשוואה הליניארית.
משוואה הומוגנית במקדמים קבועים
תהא משוואה ליניארית הומוגנית כאשר כל המקדמים הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ A_{n-1},\dots,A_0} הם מספרים קבועים. הפתרון הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y=0} הוא טריוויאלי, ואנו רוצים למצוא פתרונות נוספים. מכיוון שכל המקדמים קבועים, יהיה נוח לחפש פתרון שהוא פונקציה שהשינוי היחיד שהיא עוברת במהלך גזירתה הוא כפל בקבוע. זוהי בדיוק פונקציית האקספוננט: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y=e^{\lambda x}} כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \lambda } הוא קבוע שאנו רוצים למצוא. נשים לב כי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y^{(k)}=\lambda^k e^{\lambda x}} . לכן, לאחר הצבת הפתרון המשוער, נקבל:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \lambda^n e^{\lambda x}+A_{n-1}\lambda^{n-1} e^{\lambda x}+ \dots +A_0\ e^{\lambda x}=0} .
ניתן לצמצם ב-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ e^{\lambda x}} כי פונקציה זו תמיד שונה מאפס. נקבל:
- .
זוהי משוואה בנעלם הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \lambda} הנקראת המשוואה האופיינית של המשוואה הדיפרנציאלית שלנו. המשפט היסודי של האלגברה מבטיח עבור משוואה ממעלה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ n} שקיימים לה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ n} שורשים (לא בהכרח שונים זה מזה, אך בעלי ריבוי כולל של הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ n} ).
לעיתים בוחרים להסתכל על אגף שמאל של המשוואה בתור פולינום ב-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \lambda} , ולהציג את הבעיה כבעיה של מציאת שורשי הפולינום, הנקרא הפולינום האופייני של המשוואה הדיפרנציאלית. שתי הבעיות זהות לגמרי.
יהא הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \alpha} שורש של הפולינום (כלומר, פתרון של המשוואה). אז הוא פתרון של המשוואה, כפי שנובע מהניחוש שלנו.
ייתכן כי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \alpha} הוא שורש מריבוי גדול מאחד, ואז כל פונקציה מהצורה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ x^k e^{\alpha x}} , כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ k} הוא מספר טבעי הקטן מריבוי השורש, היא פתרון של המשוואה שאינו תלוי בשאר הפתרונות שהתקבלו בדרך זו. כך עבור משוואה מסדר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ n} ניתן לקבל הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ n} פתרונות בלתי תלויים.
אם כל מקדמי המשוואה הם ממשיים וחלק מהשורשים שיתקבלו הם מרוכבים, ניתן לקבל מכל פתרון מרוכב פתרון ממשי, שאינו כולל מספרים מרוכבים, בצורה זו: מכיוון שכל מקדמי המשוואה ממשיים, מספר מרוכב הוא פתרון שלה רק אם גם הצמוד שלו הוא פתרון שלה. נניח כי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \alpha+i\beta,\alpha-i\beta} הם שני פתרונות שכאלו. אז הפונקציות המתאימות להם הן .
כעת, בהסתמך על העובדה שסכום וכפל בקבוע של פתרונות גם הוא פתרון, נקבל שני פתרונות ממשיים על ידי שימוש בנוסחת אוילר:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ e^{\alpha x} \frac{e^{i\beta x}+e^{-i\beta x}}{2}=e^{\alpha x}\cos(\beta x), e^{\alpha x} \frac{e^{i\beta x}-e^{-i\beta x}}{2i}=e^{\alpha x}\sin(\beta x) } .
ניתן להוכיח על ידי בדיקה (למשל באמצעות הורונסקיאן) כי אלו פתרונות בלתי תלויים.
שיטת ד'אלמבר להורדת סדר המשוואה
חלק ניכר מן התאוריה הבסיסית של משוואות דיפרנציאליות נשען על האנלוגיה בין משוואות דיפרנציאליות ליניאריות למשוואות פולינומיות. בין היתר, כאשר ידוע שורש אחד הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \alpha} של משוואה פולינומית ממעלה n, אז ניתן לפרק לגורמים את הפולינום באמצעותו (כלומר, לחלק את הפולינום בגורם ליניארי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle x-\alpha} , המערב את השורש הידוע) ולקבל פולינום חדש ממעלה n-1. באופן דומה, כאשר ידוע פתרון אחד למשוואה דיפרנציאלית ליניארית הומוגנית כללית, ניתן להוריד בעזרתו את סדר המשוואה ב-1. יותר מכך, אם ידועים k פתרונות בלתי תלויים למשוואה, אז ניתן להוריד באמצעותם את סדר המשוואה עד ל-n-k.
כאשר נתונה משוואה ליניארית הומוגנית כללית (כלומר, המקדמים הם פונקציות של הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ x} ולא בהכרח קבועים), וידוע לנו פתרון לא טריוויאלי אחד של המשוואה, נרצה למצוא פתרון נוסף, בלתי תלוי בו. הגיוני לחפש פתרון שיהיה דומה בצורתו לפתרון הקיים, וההבדל ביניהם מתבטא בכפל בפונקציה לא ידועה כלשהי. לכן נסמן הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_2(x)=y_1(x)v(x)} כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ v(x)} היא פונקציה בלתי ידועה ו-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1} הוא הפתרון שידוע לנו. כעת נציב את הפתרון החדש למשוואה ונחלץ את הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ v(x)} . אם נצליח, קיבלנו פתרון נוסף, הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1(x)v(x)} .
נדגים את התהליך עבור משוואה ליניארית ממעלה שנייה: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y''+p(x)y'+q(x)y=0} . נניח שנתון פתרון לא טריוויאלי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1} ואנו מנחשים פתרון הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_2=y_1\cdot v} . נחשב את הנגזרת הראשונה והשנייה של הפתרון:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_2'=y_1'v+y_1v'}
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_2''=y_1''v+2y_1'v'+y_1v''} .
נציב במשוואה המקורית ונקבל:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1''v+2y_1'v'+y_1v''+p(x)(y_1'v+y_1v')+q(x)y_1v} .
לאחר פתיחת סוגריים והוצאת גורם משותף נקבל:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ v(y_1''+p(x)y_1'+q(x)y_1)+v'(2y_1'+p(x)y_1)+v''y_1=0} .
נשים לב כי הביטוי בסוגריים השמאליים מתאפס, כי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1} . הוא פתרון של המשוואה הליניארית (ולכן הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1''+p(x)y_1'+q(x)y_1=0} ) לכן בסך הכול קיבלנו:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ v'(2y_1'+p(x)y_1)+v''y_1=0} .
כעת ניתן להציב הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ t=v'} ולקבל אחרי חילוק ב-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1} (מותר לנו כי על פי משפט הקיום והיחידות, אם הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1} מתאפס בנקודה כלשהי הוא חייב להיות הפתרון הטריוויאלי הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y=0} ) את המשוואה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ t'+(2\frac{y_1'}{y_1}+p(x))t=0} . זוהי משוואה ליניארית מסדר ראשון, כלומר הורדנו את סדר המשוואה.
פתרון של משוואה זו נתון על ידי
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ t=\exp(-\int(2\frac{y_1'}{y_1}+p(x))dx)=\exp(-\int p(x)dx-2\ln(y_1))= e^{-\int p(x)dx}\cdot e^{\ln(y_1^{-2})}=\frac{e^{-\int p(x)dx}}{y_1^2}}
כדי למצוא את הפתרון השני של המשוואה יש לעבור מ-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ t} ל- על ידי אינטגרציה נוספת, ולבסוף לכפול ב-הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1} . נקבל את התוצאה הסופית:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_2=y_1\cdot\int\frac{e^{-\int p(x)dx}}{y_1^2}} .
ישנן דרכים נוספות להגיע לנוסחה זו, איך היתרון בדרך שהוצגה כאן הוא שאין היא מצריכה זכירת נוסחאות או משפטים, אלא רק שבהינתן פתרון הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1} כלשהו ניתן למצוא פתרון נוסף הנתון על ידי הצורה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1v} , גזירתו והצבתו במשוואה.
נשים לב כי שיטה זו אינה מבטיחה שהפתרון הנוסף יהיה ניתן לכתיבה על ידי פונקציות אלמנטריות: אין זה מובטח שניתן יהיה למצוא פתרון לשני האינטגרלים שבנוסחא. עם זאת, כאשר פונקציה מוצגת על ידי אינטגרלים לרוב קל יותר לעבוד איתה מאשר במצב שבו הצגתה היחידה היא על ידי משוואה דיפרנציאלית.
המקרה הכללי
ניתן להכליל את השיטות וההגדרות למשוואה דיפרנציאלית בפונקציות וקטוריות כלומר המשוואה מהצורה:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{d}{dx} \vec y (x) = A(x) \cdot \vec y (x) + \vec b (x)} כאשר הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec b , \ \vec y : \mathbb{R} \to \mathbb{R}^n} , ו- A היא מטריצה של פונקציות ממשיות.
אפשר להציג את y במפורש כוקטור שרכיביו הן פונקציות ממשיות גזירות, ולקבל מערכת של משוואות דיפרנציאליות עם n פונקציות נעלמות- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1 , . . . y_n} . בהמשך החלק נוותר על סימון החץ מעל הפונקציה y, כאשר ההקשר יהיה ברור.
שיטות כלליות
כמו בפתרון הבעיות הרגילות, גם בממדים גבוהים אוסף כל הפתרונות של משוואה דיפרנציאלית ליניארית הוא מרחב וקטורי, וממדו הוא ממד הווקטור y. כדי לייצג את מרחב הפתרונות באופן קומפקטי משתמשים במטריצה של פונקציות הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \Phi (x)} , שהיא הפיכה לכל t בתחום הפתרון של המשוואה, ומקיימת את השוויון הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \Phi '(x) = A(x) \Phi (x)} . מטריצה זו נקראת המטריצה היסודית של המשוואה, או המטריצה הפונדמנטלית. באמצעות מטריצה זו ניתן להציג כל פתרון של המשוואה ההומוגנית, הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y' = Ay} , כיוון שכל פתרון כזה מתקבל על ידי כפל של וקטור קבוע במטריצה היסודית. פתרונות של המשוואה הלא הומוגנית הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y' = Ay + b} מתקבלים על ידי כפל של המטריצה היסודית בוקטור של פונקציות גזירות שמקיים את המשוואה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \Phi (x) c'(x) = b(x)} .
נפתור קודם כל את המקרה הפשוט ביותר: הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y' = ay , \ \ y(0)=b } כאשר y פונקציה רגילה מהממשיים לעצמם. נשתמש באיטרציות של פיקאר כדי לפתור את המשוואה הזו:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \varphi_0 (x) = b}
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \varphi _1 (x) = b + \int_{0}^x a\varphi_0 (t) \, dt = (1+ax)b}
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \varphi _{n+1} (x) = b + \int_{0}^x a\varphi _n (t) \, dt = b + \int_{0}^x a\cdot \left( \sum_{k=0}^n \frac{a^k t^k}{k!} \right) b \, dt = \left( \sum_{k=0}^{n+1} \frac{a^k x^k}{k!} \right) b}
כאשר את המעבר השני אפשר להוכיח באינדוקציה. קל לראות שהאיטרציות מתכנסות לפונקציה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y = e^{ax} b} , שפותרת את המשוואה. בצורה דומה, ניתן לראות שהפתרון הכללי עבור תנאי התחלה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y(\tau ) = b} הוא:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle y = e^{(t-\tau ) a} \cdot b}
למעשה, אותה שיטה בדיוק יכולה לפתור גם את המשוואה הדיפרנציאלית הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y' = A y} כאשר A מטריצה ריבועית. על ידי האיטרציות של פיקאר מתקבל, כמו קודם, שהפתרון למשוואה עם תנאי ההתחלה הוא הגבול:
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y = \lim_{n \to \infty } \sum_{k=0}^n \frac{(x - \tau )^n A^n}{n!} \cdot b = \sum_{k=0}^\infty \frac{ \left( (x - \tau)A \right) ^n}{n!} \cdot b = \ e^{(x - \tau ) A} \cdot b }
האיבר האחרון בשוויונות הוא האקספוננט המטריציאלי של המטריצה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ (x - \tau ) A} מוכפל בוקטור b. המטריצה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ \Phi (x) = e^{ x A }} היא מטריצה יסודית של המשוואה.
באופן דומה למקרה הכללי, אפשר להראות שהמטריצה היסודית שמתאימה למשוואה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y' = A(x) y} היא המטריצה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ e^{ \int A(x)dx }} .
רדוקציה למשוואה מסדר n
כל משוואה ליניארית מסדר n ניתנת להצגה כמערכת של n משוואות ליניאריות על ידי המעבר הבא: אם הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y^{(n)}+p_{n-1}(x)y^{(n-1)}+\dots+p_1(x)y'+p_0(x)y=g(x)}
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y=y_1}
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_1 ' = y_2}
- .
- .
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_{n-1} ' = y_n}
- הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y_n ' = -p_{n-1}(x)y_n - \dots - p_1(x) y_2 - p_0 (x) y_1 + g(x) }
המשוואה האחרונה מתקבלת מהמשוואה המקורית הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y^{ (n) } = -(\ p_{n-1}(x)y^{(n-1)}+\dots+p_1(x)y'+p_0(x)y\ ) + g(x)} על ידי ההחלפה הפענוח נכשל (SVG (אפשר להפעיל MathML בעזרת הרחבת דפדפן): תשובה בלתי־תקינה ("Math extension cannot connect to Restbase.") מהשרת "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ y^{ (k) } = y_{k+1}} .
באופן הזה מתקבלת מערכת משוואות ליניאריות שפתרונה שקול לפתרון המשוואה המקורית.
קישורים חיצוניים
מיזמי קרן ויקימדיה |
---|
ערך מילוני בוויקימילון: מדלי"ה |
35676339משוואה דיפרנציאלית ליניארית