שינויים

קפיצה אל: ניווט, חיפוש

מדר קיץ תשעב/סיכומים/תקציר

נוספו 753 בתים, 18:39, 13 באוקטובר 2012
* '''וריאציית הפרמטרים:''' נתונים <math>y_1,\dots,y_n</math> פתרונות בת״ל של המד״ר הלינארית־הומוגנית המתאימה. אזי הפתרון הכללי של המד״ר הוא <math>\sum_{k=1}^n y_k(x)\int c_k'(x)\mathrm dx</math> כאשר <math>\begin{pmatrix}y_1(x)&\cdots&y_n(x)\\y_1'(x)&\cdots&y_n'(x)\\\vdots&\ddots&\vdots\\y_1^{(n-1)}(x)&\cdots&y_n^{(n-1)}(x)\end{pmatrix}\begin{pmatrix}c_1'\\c_2'\\\vdots\\c_n'\end{pmatrix}=\begin{pmatrix}0\\0\\\vdots\\f(x)\end{pmatrix}</math>. באופן שקול: <math>c_k'(x)=\frac{W_k(x)}{W(x)}</math>, כאשר <math>W_k(x)=\begin{vmatrix}y_1(x)&\cdots&y_{k-1}(x)&0&y_{k+1}(x)&\cdots&y_n(x)\\y_1'(x)&\cdots&y_{k-1}'(x)&0&y_{k+1}'(x)&\cdots&y_n'(x)\\\vdots&\ddots&\vdots&\vdots&\vdots&\ddots&\vdots\\y_1^{(n-1)}(x)&\cdots&y_{k-1}^{(n-1)}(x)&f(x)&y_{k+1}^{(n-1)}(x)&\cdots&y_n^{(n-1)}(x)\end{vmatrix}</math>.
* נניח שהמד״ר לינארית־הומוגנית עם מקדמים קבועים. אזי נציב <math>y=\mathrm e^{rx}</math>, ולכן <math>y^{(k)}=r^k\mathrm e^{rx}</math> וגם <math>r^n+\sum_{k=0}^{n-1} a_k r^k</math> (זה הפולינום האופייני של המשוואה) שווה ל־0. אם השורשים השונים זה מזה הם <math>r_1,\dots,r_m</math> והריבויים שלהם <math>d_1,\dots,d_m</math> בהתאמה אזי הפתרון הכללי הוא <math>y=\sum_{k=1}^m\mathrm e^{r_kx}\sum_{i=0}^{d_k-1}c_{k,i}x^i</math>. אם <math>r_k</math> אינו ממשי ניתן לכתוב <math>r_k=\alpha+\beta\mathrm i</math> ואז, כיוון ש־<math>\overline{r_k}</math> שורש עם אותו ריבוי, נציב <math>C_1\mathrm e^{r_kx}+C_2\mathrm e^{\overline{r_k}x}=\mathrm e^{\alpha x}\Big(c_1\cos(\beta x)+c_2\sin(\beta x)\Big)</math>.
:* '''שיטת הניחוש/הבחירה/המקדמים הנעלמים:''' נניח שהמד״ר לינארית עם מקדמים קבועים וכן <math>f(x)=\mathrm e^{\lambda x}\sum_{k=0}^m b_k x^k</math>, כאשר <math>\lambda</math> קבועה (יכולה להיות גם 0), והריבוי של <math>\lambda</math> בפולינום האופייני הוא <math>d</math> (במידה ו־<math>\lambda</math> לא שורש נאמר <math>d=0</math>). אזי קיים פתרון פרטי מהצורה <math>\mathrm e^{\lambda x}x^d\sum_{k=0}^m B_k x^k</math> כאשר <math>b_m,B_m\ne0</math>. ''הערה:'' אם <math>f(x)=g(x)+h(x)</math> נוכל לפתור עבור <math>g(x),h(x)</math> בנפרד ולסכום את הפתרונות הפרטיים.* '''משוואת אוילר(־לגראנג׳–לגראנג׳)''' היא מד״ר לינארית מהצורה <math>(x-x_0)^ny^{(n)}+\sum_{k=0}^{n-1}a_k (x-x_0)^k y^{(k)}=f(x)</math> עם <math>\forall k:\ a_k=\text{const.}</math>. מציבים <math>x-x_0=\begin{cases}\mathrm e^t,&x>x_0\\-\mathrm e^t,&x<x_0\end{cases}</math> במד״ר ההומוגנית ואז <math>y'=\frac{\mathrm dy}{\mathrm dt}\frac{\mathrm dt}{\mathrm dx}=\frac{\mathrm dy}{\mathrm dt}\frac1{x-x_0},\ y''=\mathrm e^{-2t}\left(\frac{\mathrm d^2y}{\mathrm dt^2}-\frac{\mathrm dy}{\mathrm dt}\right),\ \dots</math>. נקבל משוואה לינארית־הומוגנית עם מקדמים קבועים, וניתן להמשיך לפתור אותה באופן זה. לחלופין, אנו לומדים מכך שאפשר להציב <math>y=(x-x_0)^r</math> במד״ר ההומוגנית ולקבל <math>r^n+\sum_{k=0}^{n-1} b_k r^k=0</math> (משוואה אינדיציאלית). אם השורשים השונים זה מזה הם <math>r_1,\dots,r_m</math> והריבויים שלהם <math>d_1,\dots,d_m</math> בהתאמה אזי הפתרון ההומוגני הכללי הוא <math>y=\sum_{k=1}^m (x-x_0)^{r_k}\sum_{i=0}^{d_k-1}c_{k,i}\ln^i(x-x_0)</math>. אם <math>r_k</math> אינו ממשי ניתן לכתוב <math>r_k=\alpha+\beta\mathrm i</math> ואז, כיוון ש־<math>\overline{r_k}</math> שורש עם אותו ריבוי, נציב <math>C_1(x-x_0)^{r_k}+C_2(x-x_0)^{\overline{r_k}}=(x-x_0)^\alpha\Big(c_1\cos(\beta\ln(x-x_0))+c_2\sin(\beta\ln(x-x_0))\Big)</math>.:* אם <math>f(x)=(x-x_0)^\lambda\sum_{k=0}^m b_k \ln^k(x-x_0)</math> כאשר <math>\lambda</math> קבועה (יכולה להיות גם 0), והריבוי של <math>\lambda</math> במשוואה האינדיציאלית הוא <math>d</math> (אם לא שורש <math>d=0</math>). אזי קיים פתרון פרטי מהצורה <math>(x-x_0)^\lambda\ln^d(x-x_0)\sum_{k=0}^m B_k \ln^k(x-x_0)</math> כאשר <math>b_m,B_m\ne0</math>.
===== פתרון מד״ר באמצעות טורי חזקות =====
* נתונה מד״ר מהצורה <math>y^{(n)}+\sum_{k=0}^{n-1} a_k(x)y^{(k)}=f(x)</math> כאשר <math>\forall k:\ f(x),a_k(x)\in C(a,b)</math> ותהי <math>x_0\in(a,b)</math>. אם <math>f</math> וכל המקדמים <math>a_k</math> אנליטיים סביב <math>x_0</math> עם רדיוס התכנסות <math>R</math> או יותר אזי קיים פתרון אנליטי סביב <math>x_0</math> של המד״ר עם רדיוס התכנסות <math>R</math> או יותר.
* '''טור פרוביניוס''' הוא טור מהצורה <math>(x-x_0)^r\sum_{k=0}^\infty a_k(x-x_0)^k</math>.
* '''שיטת פרוביניוס:''' בהנתן <math>a_2(x)y''+a_1(x)y'+a_0(x)y=0</math> נחלק ב־<math>a_2(x)</math>. תהי <math>x_0</math> נקודה סינגולרית של <math>\frac1{a_2(x)}</math>. אם קיימים הגבולות <math>L_kL_1=\lim_{x\to x_0}(x-x_0)^{2-k}\frac{a_ka_1(x)}{a_2(x)}</math> הנקודה נקראת סינגולרית־רגולרית. בקרבת <math>x_0</math> נקבל <math>0,L_2=(x-x_0)^2y''+\fraclim_{a_1(x)}{a_2(x)\to x_0}(x-x_0)y'+^2\frac{a_0(x)}{a_2(x)}y=</math> הנקודה נקראת סינגולרית־רגולרית. בקרבת <math>x_0</math> נקבל <math>(x-x_0)^2y''+(L_1+o(1))(x-x_0)y'+(L_0L_2+o(1))y=0</math>. לפי משפט, אם <math>x_0</math> נקודה סינגולרית־רגולרית אזי קיים פתרון אנליטי למד״ר סביב <math>x_0</math> בצורת בצורת טור פרוביניוס. לכן נפתור עבור נחשב את הפולינום האופייני של המד״ר עם <math>o(1)=0</math>, נציב ע״י הצבת <math>y=(x-x_0)^r</math> ונקבל את הפתרונות בצורת טורים של המד״ר עם , ואם פתרונות הפולינום הם <math>o(1)=0r_1,r_2</math> (אם אזי יש שני פתרונות הפולינום האופייני של המד״ר עם פרטיים בת״ל מהצורה <math>oy_k=(1x-x_0)^{r_k}\sum_{i=0}^\infty b_{k,i}(x-x_0)^i,\quad k\in\{1,2\}</math> הם כאשר <math>r_1,-r_2\not\in\mathbb Z</math> אז ו־<math>y_ky_1=(x-x_0)^{r_kr_1}\sum_{i=0}^\infty b_{k1,i}(x-x_0)^i,y_2=y_1\ln(x-x_0)+(x-x_0)^{r_2}\sum_{i=0}^\infty b_{2,i}(x-x_0)^i</math> פתרון פרטי)כאשר <math>r_1-r_2\in\mathbb Z</math> ומתקיים בה״כ <math>r_1\ge r_2</math>. נציב פתרונות אלו אותם במד״ר המקורית ונקבל את מקדמי הטורים.<br>''הערה:'' נאמר ש־<math>f\in o(g)</math> אם <math>\lim_{x\to x_0}\frac{f(x)}{g(x)}=0</math>. לעתים כותבים "<math>o(1)</math>" לציון איבר הנמצא בקבוצה זו, ולא הקבוצה עצמה.:* '''פונציית גמא:''' <math>\forall x>0:\ \Gamma(x):=\int\limits_0^\infty t^{x-1}\mathrm e^{-t}\mathrm dt</math>. היא מקיימת <math>\Gamma(xּּ+1)=x\Gamma(x)</math> וגם <math>\forall n\in\mathbb N:\ \Gamma(n)=(n-1)!</math>. ערך חשוב: <math>\Gamma\!\left(\frac12\right)=\sqrt\pi</math>.:* '''פונקציית בסל:''' <math>J_m(x)=\sum_{k=0}^\infty \frac{(-1)^k \Gamma(m+1)}{2^{2k}k!\Gamma(m+k+1)}x^{2k+m}</math> כאשר <math>m</math> נקראת דרגת הפונקציה.:* '''משוואת בסל:''' <math>x^2y''+xy'+(x^2-m^2)y=0</math>. מתקיים <math>y''+\frac1xy'+\left(1-\frac{m^2}{x^2}\right)=0</math> ולכן <math>\lim x\frac1x=1,\ \lim x^2\left(1-\frac{m^2}{x^2}\right)=-m^2</math>, כלומר 0 סיגולריות־רגולרית. השורשים של הפולינום האופייני הם <math>0\pm m</math> סיגולריות־רגולרית.:* '''פונציית גמא:''' ולכן, כאשר <math>m\Gamma(x):=not\intin\limits_0^frac12\infty t^{x-1}\mathrm e^{-t}\mathrm dtmathbb Z</math>. היא מקיימת , הפתרון הכללי הוא <math>\Gammac_1 J_m(xּּx)+1)=x\Gammac_2 J_{-m}(x)</math> וגם . אם <math>\forall nm\in\frac12\mathbb N:\ \GammaZ</math> אז הפתרון הכללי הוא <math>c_1 J_m(nx)=+c_2 N_m(n-1x)!</math>.:* '''משוואת אוילר:''' כאשר <math>N_m(x^2y''+xy')=\frac{J_m\cos(\pi m)-J_{-m^2y=0}(x)}{\sin(\pi m)}</math>.