וריאציית הפרמטר

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

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

==היסטוריה==
שיטת וריאציית הפרמטר הוצגה לראשונה במאה ה-18 על ידי המתמטיקאי השווייצרי לאונרד אוילר והשלים אותה המתמטיקאי האיטלקי ז'וזף לואי לגרנז'.

==מד"ר מסדר ראשון==
נתונה משוואה דיפרנציאלית מסדר ראשון. נסדר אותה באופן הבא, כשהיא "מנורמלת" (מקדם <math>y'</math> ערכו 1):
:<math> y' + p(x)y = q(x) </math>
הפתרון הכללי של המשוואה מורכב מסכום הפתרון הכללי של המשוואה ההומוגנית המתאימה (הרשומה מטה) ומפתרון פרטי של המשוואה הלא-הומוגנית.

ראשית נמצא פתרון למשוואה ההומוגנית המתאימה:
: <math> y' + p(x)y = 0</math>.
ניתן לפתור את המשוואה ההומוגנית במגוון דרכים, ביניהן הפרדת משתנים:

<math>\frac{\mathrm{d}}{\mathrm{d}x} y + p(x)y = 0 </math>

:<math>\frac{\mathrm{d}y}{\mathrm{d}x}=-p(x)y </math>
<math>\frac{\mathrm{d}y}{y} = -p(x)\,\mathrm{d}x</math>
:<math>\int \frac{1}{y} \, \mathrm{d}y = -\int p(x) \, \mathrm{d}x </math>
:<math>\ln |y| = -\int p(x) \, \mathrm{d}x + C_0 </math>

<math>|y|=C_1 e^{-\int p(x)dx}, C_1>0</math>

:<math>y = \pm C_1 e^{-\int p(x) \, \mathrm{d}x} = C_2 e^{-\int p(x) \, \mathrm{d}x} , C_2 \neq 0</math>
:אך מכיוון שגם <math>y=0</math> הוא פתרון ניתן לכלול גם אותו, ולהרשות כל <math>C_2 \in \mathbb{R}</math>.

כלומר, פתרון המשוואה ההומוגנית הוא:
:<math>y_H(x) = C e^{-\int p(x) \, dx}</math>.

כעת נחזור למשוואה המקורית הלא-הומגנית:

: <math> y' + p(x)y = q(x)</math>

נשתמש בשיטת וריאציית הפרמטר - באמצעות החלפת הפרמטר <math>C</math> של הפתרון הפרטי בפונקציה <math>C(x)</math>:

:<math>y_p(x) = C(x) e^{-\int p(x) \, \mathrm{d}x}</math>

באמצעות הצבת הפתרון הפרטי לתוך המשוואה המקורית הלא-הומוגנית ניתן למצוא <math>C(x)</math>:

:<math> C' (x) e^{-\int p(x) \, \mathrm{d}x} - C(x) p(x) e^{-\int p(x) \, \mathrm{d}x} + p(x) C(x) e^{-\int p(x) \, \mathrm{d}x} = q(x)</math>

: <math> C' (x) e^{-\int p(x) \, dx} = q(x)</math>

: <math> C' (x) = q(x) e^{\int p(x) \, dx} </math>

: <math> C(x) =\int q(x) e^{\int p(x) \, dx} \, dx + C_1 </math>

אנו נדרשים רק לפתרון פרטי, לכן נבחר <math>C_1=0</math> לצורך פשטות, כאשר הפתרון הפרטי הוא: <math>y_p =e^{-\int p(x) \, dx} \int q(x) e^{\int p(x) \, dx} \, \mathrm{d}x</math>

לבסוף נסכום את שני הפתרונות הפרטיים שמצאנו לכדי פתרון כללי
:<math>y(x) = y_H(x) + y_p(x)</math>
:<math> y(x) = C e^{-\int p(x) \, \mathrm{d}x} + e^{-\int p(x) \, \mathrm{d}x} \int q(x) e^{\int p(x) \, \mathrm{d}x} \, \mathrm{d}x</math>

== מד"ר מסדר שני ==
בהינתן משוואה דיפרנציאלית רגילה מסדר שני: <math>y+p(x)y'+q(x)y=g(x)</math>, נניח שבידינו פתרונות המשוואה ההומוגנית <math>y_1(x),y_2(x)</math> וננסה פתרון מהצורה:

<math>y(x)=u_1(x)y_1(x)+u_2(x)y_2(x)</math>.

נגזור לפי <math>x</math> ונקבל: <math>y'=u_1' y_1+u_1 y_1'+u_2' y_2+u_2 y_2'</math> 

בנוסף נדרוש שהפונקציות <math>u_1,u_2</math> תקיימנה את התנאי: <math>y_1u_1'+y_2u_2'=0</math>

ולכן <math>y'=u_1 y_1'+u_2 y_2'</math>.

נגזור פעם שנייה לפי <math>x</math> ונקבל:

<math>y=u_1' y_1'+u_1 y_1+u_2' y_2'+u_2 y_2</math>.

נציב את הנגזרות של הפתרון לתוך המשוואה הדיפרנציאלית: <math>u_1' y_1'+u_1 y_1+u_2' y_2'+u_2 y_2+p(u_1 y_1'+u_2 y_2')+q(u_1y_1+u_2y_2)=g(x)</math>.

נסדר את הביטוי מחדש: <math>u_1(y_1+py_1'+qy_1)+u_2(y_2+py_2'+qy_2)+u_1' y_1'+u_2' y_2'=g(x)</math>.

<math>y_1, y_2</math> הם פתרונות של המשוואה ההומוגנית, ולכן <math>y_{1/2}+py_{1/2}'+qy_{1/2}=0
</math>.

בינתיים יש לנו שתי משוואות שבעזרתן נוכל למצוא את <math>u_1,u_2</math>: 

<math>\begin{cases} y_1u_1'+y_2u_2'=0 \\ y_1'u_1'+ y_2'u_2'=g(x) \end{cases}</math>

נוכל לחלץ את הערכים של <math>u_1',u_2'</math>, ולרשום אותם דרך הורונסקיאן:

<math>u_1'=-{y_2(x)g(x) \over W(y_1,y_2)}</math>,<math>u_2'={y_1(x)g(x) \over W(y_1,y_2)}</math>.

ומכאן: <math>u_1=-\int {y_2(x)g(x) \over W(y_1,y_2)} \mathrm{d}x</math>,<math>u_2=\int {y_1(x)g(x) \over W(y_1,y_2)}\mathrm{d}x</math>.

נרשום את הפתרון הפרטי למשוואה: <math>y_p=-y_1\int {y_2(x)g(x) \over W(y_1,y_2)} \mathrm{d}x+y_2\int {y_1(x)g(x) \over W(y_1,y_2)}\mathrm{d}x</math>.

למעשה קיבלנו את הפתרון הכללי למשוואה הדיפרנציאלית: <math>y(x)=Ay_1+By_2+y_p</math> כאשר הפרמטרים <math>A</math> ו-<math>B</math> מתקבלים מתנאי ההתחלה.

==קישורים חיצוניים==
* 
* 

==הערות שוליים==

קטגוריה:משוואות דיפרנציאליות
קטגוריה:לאונרד אוילר