שינויים
/* מבחנים לדוגמא */
[[88-240 משוואות דיפרנציאליות רגילות]]
=מבחנים לדוגמא=
*[[מדיה:18EngODEExmpTest1.pdf|מבחן לדוגמא 1]], [[מדיה:18EngODEExmpTest1Sol.pdf|פתרון]]
*[[מדיה:18EngODEExmpTest2.pdf|מבחן לדוגמא 2]], [[מדיה:18EngODEExmpTest2Sol.pdf|פתרון]]
*[[מדיה:18EngODETestA.pdf|מבחן מועד א' תשע"ח]], [[מדיה:18EngODETestASol.pdf|פתרון]]
*[[מדיה:18EngODETestB.pdf|מבחן מועד ב' הנדסה תשע"ח]], [[מדיה:18EngODETestBSol.pdf|פתרון מבחן מועד ב' הנדסה תשע"ח]]
*[[מדיה:19ODETestA.pdf|מבחן מועד א' תשע"ט]], [[מדיה:19ODETestASol.pdf|פתרון]]
*[[מדיה:19ODETestB.pdf|מבחן מועד ב' תשע"ט]], [[מדיה:19ODETestBSol.pdf|פתרון מבחן מועד ב' תשע"ט]]
*[[מדיה:21ODETestA.pdf|מבחן מועד א' תשפ"א]], [[מדיה:21ODETestASol.pdf|פתרון מבחן מועד א' תשפ"א]]
*[[מדיה:21ODETestB.pdf|מבחן מועד ב' תשפ"א]], [[מדיה:21ODETestBSol.pdf|פתרון מבחן מועד ב' תשפ"א]]
*[[מדיה:22ODETestA.pdf|מבחן מועד א' תשפ"ב]], [[מדיה:22ODETestASol.pdf|פתרון מבחן מועד א' תשפ"ב]]
*[[מדיה:22ODETestB.pdf|מבחן מועד ב' תשפ"ב]], [[מדיה:22ODETestBSol.pdf|פתרון מבחן מועד ב' תשפ"ב]]
*[[מדיה:23ODEQuiz.pdf|בוחן תשפ"ג]], [[מדיה:23ODEQuizSol.pdf|פתרון בוחן תשפ"ג]]
*[[מדיה:23ODETestA.pdf|מבחן מועד א' תשפ"ג]], [[מדיה:23ODETestASol.pdf|פתרון מבחן מועד א' תשפ"ג]]
*[[מדיה:23ODETestB.pdf|מבחן מועד ב' תשפ"ג]], [[מדיה:23ODETestBSol.pdf|פתרון מבחן מועד ב' תשפ"ג]]
*[[מדיה:23EngODEQuiz.pdf|בוחן הנדסה תשפ"ג]], [[מדיה:23EngODEQuizSol.pdf|פתרון בוחן הנדסה תשפ"ג]]
*[[מדיה:23EngODETestA.pdf|מבחן מועד א' הנדסה תשפ"ג]], [[מדיה:23EngODETestASol.pdf|פתרון]]
*[[מדיה:23EngODETestB.pdf|מבחן מועד ב' הנדסה תשפ"ג]], [[מדיה:23EngODETestBSol.pdf|פתרון]]
===מבחנים של מד"ר למדעי המוח===
*[[מדיה:23BSODETestA.pdf|מבחן מועד א' תשפ"ג]], [[מדיה:23BSODETestAPartialSol.pdf|פתרון חלקי מבחן מועד א' תשפ"ג]]
*[[מדיה:23BSODETestB.pdf|מבחן מועד ב' תשפ"ג]], [[מדיה:23BSODETestBPartialSol.pdf|פתרון חלקי מבחן מועד ב' תשפ"ג]]
=הרצאות=
[https://www.youtube.com/playlist?list=PLzSjdxrZD_hlMJrtihLjrl0d55Zk4Ggy6 פלייליסט של ההרצאות למחלקת מתמטיקה שנת תשפ"א]
==הרצאה 1 הקדמה ומשוואה פרידה==
*כאן הפונקציה נתונה, ואנו מחפשים ערך של המשתנה שמקיים את המשוואה.
*המלצה: ניתן להעזר בספר המצויין על מד"ר של סמי זעפרני ב[https://samyzaf.com/technion/ode/ode.pdf קישור הבא].
===נפילה חופשית===
===סדר ומעלההמד"ר===
*משוואה דיפרנציאלית נקראת '''מסדר''' n אם הנגזרת הגבוהה ביותר היא מסדר n.
**המשוואה <math>y''=g</math> היא משוואה מסדר שני.
**המשוואה <math>y'=ry</math> היא משוואה מסדר ראשון.
===משוואות פרידות===
*<math>f(y)y'=g(x)</math>
*הקדומות של שני הצדדים שוות עד כדי קבוע.
*<math>\int f(y)y'dx=\{t=y(x),dt=y'dx\}=\int f(t)dt=F(y)</math>*במקום t נשאר עם המשתנה לכן ביחד נקבל <math>F(y ובעצם )=G(x)+c</math>*בעצם אנו מחשבים אינטגרלים לשני הצדדים <math>f(y)dy=g(x)dx</math>, כל אחד לפי המשתנה שלו!
*בתרגיל זה איננו צריכים, כי מצאנו את הפתרון הכללי בדרך פשוטה יותר למעלה.
*בהמשך, משפט הקיום והיחידות יעזור לנו להתמודד עם השאלה הזו, אך באופן כללי לא נעסוק הרבה במקרי קצה בקורס זה.
====המרדף====
*דוגמא יפה וחשובה מ[http://people.uncw.edu/hermanr/mat361/ODEBook/ODE1.pdf הספר הזה] עמוד 19 של הספר (33 של הPDF)
*מרצה צועד במהירות קבועה <math>b</math> בקו ישר בשדרה שמוביל אל בניין 507.
*סטודנט שרוצה עוד שתי נקודות לעובר רואה את המרצה, ונע לכיוון המרצה במהירות קבועה <math>c</math>.
*המרצה מתחיל בנקודה <math>(0,0)</math> ונע בכיוון החיובי של ציר y, הסטודנט מתחיל בנקודה <math>(a,0)</math> עבור <math>a>0</math>.
*באיזה מסלול ינוע הסטודנט? באילו תנאים הוא יתפוס את המרצה?
*נסמן את פונקצית המסלול של הסטודנט ב<math>y(x)</math>
*כיוון שהסטודנט תמיד נע בכיוון המרצה, המשיק של הפונקציה בכל נקודה במסלול הסטודנט צריך לפגוש את המרצה באותו הזמן.
*בזמן <math>t</math> המרצה נמצא בנקודה <math>(0,b\cdot t)</math> והסטודנט נמצא בנקודה <math>(x,y)</math>.
*השיפוע בין המרצה לסטודנט הוא הנגזרת של פונקצית המסלול, כלומר <math>y'=\frac{y-bt}{x}</math>
*כעת יש לנו שלושה משתנים <math>t,x,y</math>, כיצד נפטר מאחד מהם? לא השתמשנו במהירות הסטודנט!
*המסלול שהסטודנט עבר צריך להיות שווה ל<math>c\cdot t</math>, כלומר <math>\int_x^a \sqrt{y'^2+1}=ct</math>
*מהמשוואה לעיל אנו יודעים כי <math>t=\frac{y-xy'}{b}</math>
*ביחד נקבל כי <math>\int_x^a \sqrt{y'^2+1}=c\cdot \frac{y-xy'}{b}</math>
*נגזור את שני הצדדים ונקבל כי:
*<math>-\sqrt{y'^2+1}=\frac{c}{b}\cdot (-xy'')</math>
*<math>\frac{c}{b}xy''=\sqrt{y'^2+1}</math>
*נסמן <math>y'=z</math> ונקבל <math>\frac{c}{b}xz'=\sqrt{z^2+1}</math>
*זו מד"ר פרידה
*<math>\frac{c}{b\sqrt{z^2+1}}dz=\frac{1}{x}dx</math>
*באמצעות [[מדיה:09Infi2Universal.pdf|ההצבה האוניברסאלית המתאימה]] <math>z=tan(t)</math> נפתור את האינטגרל של הצד השמאלי ונקבל כי
*<math>\frac{c}{b}ln(\sqrt{z^2+1}+z)=ln(x)+D</math>
*ברגע הראשון התקיים כי <math>x=a</math> והתלמיד כיוון לראשית הצירים כלומר <math>y'(a)=0</math> כלומר <math>z(a)=0</math>
*לכן <math>\frac{c}{b}ln(\sqrt{z^2+1}+z)=ln(x)-ln(a)</math>
*<math>ln(\sqrt{z^2+1}+z)=\frac{b}{c}ln(\frac{x}{a})</math>
*<math>\sqrt{z^2+1}+z=\left(\frac{x}{a}\right)^{\frac{b}{c}}</math>
*כעת קצת אלגברה:
*<math>z+\sqrt{z^2+1}=A</math>
*<math>\frac{-1}{z-\sqrt{z^2+1}}=A</math>
*<math>z-\sqrt{z^2+1}=-\frac{1}{A}</math>
*נחבר למשוואה הראשונה
*<math>z=\frac{1}{2}\left(A-\frac{1}{A}\right)</math>
*הרי <math>z=y'</math>, ולכן ביחד:
*<math>y'=\frac{1}{2}\left(\left(\frac{x}{a}\right)^{\frac{b}{c}}-\left(\frac{x}{a}\right)^{-\frac{b}{c}}\right)</math>
*ולכן אחרי אינטגרציה נקבל כי:
*<math>y=\frac{a}{2}\left(\frac{1}{\frac{b}{c}+1}\left(\frac{x}{a}\right)^{\frac{b}{c}+1} -
\frac{1}{1-\frac{b}{c}}\left(\frac{x}{a}\right)^{1-\frac{b}{c}}\right) + K</math>
*כאשר אנחנו מקבלים את הקבוע <math>K</math> מהנתון <math>y(a)=0</math>
*באופן טבעי, אם מהירות המרצה גדולה ממהירות הסטודנט <math>b>c</math> נקבל שאיפה לאינסוף כאשר <math>x\to 0</math> והסטודנט לא יגיע למרצה.
*אם <math>b<c</math> הסטודנט יגיע לשדירה ויתפוס את המרצה.
*אם <math>b=c</math> האינטגרציה שלנו שגוייה, וכאשר נחשב אותה נכון שוב נקבל שאיפה לאינסוף (באופן טבעי)
===הפיכת משוואה לפרידה===
*על מנת לדעת אם לא פספסנו פתרונות אחרים, נעזר בהמשך במשפט הקיום והיחידות.
*אבל כאמור - אנחנו לא נתייחס באופן כזה לכל מקרה קצה בהמשך הקורס.
==הרצאה 2 מד"ר הומוגנית, מד"ר לינאריות מסדר ראשון ומשוואת ברנולי==
===מד"ר הומוגנית===
*מד"ר הומוגנית (בניגוד למד"ר לינארית הומוגנית שנראה בהמשך) היא משוואה מהצורה <math>y'=g(\frac{y}{x})</math>.
*נפתור מד"ר הומוגנית באמצעות ההצבה <math>z=\frac{y}{x}</math> באופן הבא:
**ראשית נסמן <math>y'=g(\frac{y}{x})</math>.
**כעת נגזור את שני צידי המשוואה <math>zx=y</math>, ונקבל כי <math>z'x+z=y'</math>.
**לכן לאחר החלפת המשתנה קיבלנו משוואה '''פרידה''' <math>z'x+z=g(z)</math>.
**נפריד את המשתנים <math>\frac{1}{g(z)-z}dz=\frac{1}{x}dx</math>.
**ולכן <math>\int \frac{1}{g(z)-z}dz=\ln|x|+C</math>.
**נמצא את <math>z</math> ונציב בחזרה <math>y=zx</math>.
*פונקציה <math>f(x,y)</math> נקראת הומוגנית מסדר k אם לכל <math>\lambda\neq 0</math> מתקיים כי <math>f(\lambda x,\lambda y)=\lambda^k f(x,y)</math>.
*לדוגמא <math>f(x,y)=\frac{x^2+xy}{x+y}</math> הומוגנית מסדר 1.
*טענה: פונקציה <math>f(x,y)</math> היא מהצורה <math>\varphig(\frac{y}{x})</math> לכל <math>x\neq 0</math> אם"ם היא הומוגנית מסדר <math>0</math> לכל <math>x\neq 0</math>.
*הוכחה:
**אם <math>f(x,y)=\varphig(\frac{y}{x})</math> אזי לכל <math>x\neq 0</math> מתקיים <math>f(\lambda x,\lambda y)=\varphig(\frac{\lambda y}{\lambda x})=\varphig(\frac{y}{x})=\lambda^0 f(x,y)</math>.**אם <math>f(\lambda x,\lambda y)=f(x,y)</math>, נציב <math>\lambda=\frac{1}{x}</math> ונקבל כי <math>f(x,y)=f(1,\frac{y}{x})=\varphig(\frac{y}{x})</math>.
*דוגמא - נפתור את המשוואה <math>y'=\frac{x^2+y^2}{xy}</math>
**<math>\varphig(\frac{y}{x})=f(1,\frac{y}{x})=\frac{1+(\frac{y}{x})^2}{\frac{y}{x}}</math>**<math>\int \frac{1}{\varphig(z)-z}dz=\int \frac{1}{\frac{1+z^2}{z}-z}dz=\int z dz=\frac{z^2}{2} </math>
**<math>\frac{z^2}{2}=ln|x|+C</math>
**<math>z=\pm\sqrt{ln(x^2)+C}</math>
*דוגמא - נפתור את המשוואה <math>xdy-\left(x\cdot\cos^2(\frac{y}{x})+y\right)dx=0</math>
**<math>y'=\frac{x\cdot\cos^2(\frac{y}{x})+y}{x}</math>
**<math>\varphig(\frac{y}{x})=f(1,\frac{y}{x})=\cos^2(\frac{y}{x})+\frac{y}{x}</math>**<math>\int \frac{1}{\varphig(z)-z}dz=\int \frac{1}{\cos^2(z)}dz=\tan(z)</math>
**<math>\tan(z)=\ln|x|+c</math>
**<math>z=\arctan(ln|x|+C)</math>
**<math>y=x\cdot \arctan(ln|x|+C)</math>
===מד"ר לינארית מסדר ראשון===
*הגדרה: משוואה מסדר ראשון נקראת לינארית אם היא מהצורה <math>y'+pa(x)\cdot y=qb(x)</math>.*מד"ר לינארית הומוגנית (בניגוד למד"ר הומוגנית שראינו לעיל) היא מהצורה <math>y'+pa(x)\cdot y=0</math>.
*נחשב נוסחא לפתרון מד"ר לינארית כללית ע"י מציאת פתרון למשוואה לינארית הומוגנית ובאמצעות שיטת וריאצית המקדמים.
*נשים לב כי המשוואה הלינארית ההומוגנית <math>y'+pa(x)\cdot y=0</math> היא '''פרידה'''.*נפריד את המשתנים ונקבל <math>\frac{1}{y}dy=-pa(x)dx</math>.*נבצע אינטגרציה ונקבל כי <math>ln|y|=-\int pA(x)dx +C</math>.*ולכן <math>y=C\cdot e^{-\int pA(x)dx}</math>
*כלומר, נציב <math>y=C(x)\cdot e^{-\int pA(x)dx}</math> במשוואה <math>y'+pa(x)y=qb(x)</math>.*נקבל <math>C'(x)\cdot e^{-\int pA(x)dx}-pa(x)\cdot C(x)\cdot e^{-\int pA(x)dx} + pa(x)\cdot C(x) \cdot e^{-\int pA(x)dx}=qb(x)</math>*משוואה זו מתקיימת אם"ם <math>C'(x)\cdot e^{-\int pA(x)dx}=qb(x)</math>.*כלומר <math>C'(x)=qb(x)\cdot e^{\int pA(x)dx}</math>.*לכן נבחר <math>C(x)=\int \left[qb(x)\cdot e^{\int pA(x)dx}\right]dx+C</math>
*סה"כ הפתרון הכללי למד"ר הלינארית <math>y'+pa(x)\cdot y=qb(x)</math> הוא: <math>e^{-\int pA(x)dx}\cdot\left(C+\int\left(qb(x)\cdot e^{\int pA(x)dx}\right)+Cdx\right)</math>
*דוגמא - המשוואה החביבה עלינו <math>y'=ry</math>:
**ראשית, נשים לב כי <math>pa(x)=-r</math> ו<math>qb(x)=0</math>.
**כלומר זו מד"ר לינארית הומוגנית, והפתרון הכללי הוא <math>y=C\cdot e^{-\int (-r)dx}=C\cdot e^{rx}</math>
*נבצע פירוק לשברים חלקיים:
*<math>\frac{1}{g-\frac{b}{m}v^2}=\frac{1}{(\sqrt{g}-\sqrt{\frac{b}{m}}\cdot v)(\sqrt{g}+\sqrt{\frac{b}{m}}\cdot v)}=\frac{1}{2\sqrt{g}}\left(\frac{1}{\sqrt{g}+\sqrt{\frac{b}{m}}\cdot v}+\frac{1}{\sqrt{g}-\sqrt{\frac{b}{m}}\cdot v}\right)</math>
*ולכן <math>\int \frac{1}{g-\frac{b}{m}v^2}dv=\frac{\sqrt{bm}}{2\sqrt{g\cdot mb}}\ln\left|\frac{\sqrt{g}+\sqrt{\frac{b}{m}}\cdot v}{\sqrt{g}-\sqrt{\frac{b}{m}}\cdot v}\right|</math>
*מצד שני <math>\int dt=t+c</math>
*לכן <math>\frac{\sqrt{g}+\sqrt{\frac{b}{m}}\cdot v}{\sqrt{g}-\sqrt{\frac{b}{m}}\cdot v}=Ce^{\left(2\sqrt{\frac{g\cdot mb}{bm}}t\right)}</math>*נסדר קצת <math>v=\sqrt{\frac{g\cdot m}{b}}\cdot \left(1-\frac{2}{1+Ce^{\left(2\sqrt{\frac{g\cdot mb}{bm}}t\right)}}\right)</math>
*נשים לב שכאשר <math>t\to\infty</math> אנו מתכנסים ל[https://en.wikipedia.org/wiki/Terminal_velocity מהירות הסופית] <math>\sqrt{\frac{g\cdot m}{b}}</math>.
*אם זו הייתה המהירות ההתחלתית היינו מקבלים פונקצית מהירות קבועה.
=====במהירות נמוכה=====
*ולכן הפתרון הוא <math>v=e^{-\frac{b}{m}t}\cdot\left(\int ge^{\frac{b}{m}t}dt+C\right)=\frac{g\cdot m}{b}+Ce^{-\frac{b}{m}t}</math>.
*וכאשר <math>t\to\infty</math> המהירות שואפת למהירות הסופית <math>\frac{g\cdot m}{b}</math>.
===משוואת ברנולי===
**לכן <math>z'-bz=d</math>
**נפתור את המשוואה הדיפרנציאלית:
***<math>z=e^{bt}\cdot (de\frac{d}{-b}e^{-bt}+C)=d+Ce^{bt}-\frac{d}{b}</math>**ולכן <math>v=\frac{1}{d+Ce^{bt}-\frac{d}{b}}</math>
**כמובן שכאשר <math>t\to\infty</math> המהירות מתכנסת מהר מאד לאפס.
*דוגמא - [https://en.wikipedia.org/wiki/Logistic_function#Applications המשוואה הלוגיסטית]**קצב הגדילה של אוכלוסיה פרופורציונלית לגודל האוכלוסיה כפול כמות המשאבים הפנויים.**המשאבים קטנים באופן פרופורציונלי לגודל האוכלוסיה.**<math>y'=a\cdot y\cdot (1-by)</math> ==הרצאה 3 משוואות מדוייקות ומשפט הקיום והיחידות==
===הקדמה - פונקציות בשני משתנים===
*מד"ר מסדר ראשון נקראת מדוייקת אם היא מהצורה <math>U_x(x,y)dx+U_y(x,y)dy=0</math>, עבור <math>U(x,y)</math> דיפרנציאבילית.
*פתרון המד"ר ניתן בצורה סתומה על ידי המשוואה <math>U(x,uy)=C</math>, כאשר C קבוע כלשהו.*תהי מד"ר מהצורה <math>Pdx+Qdy=0</math> היא מדוייקת אם"ם <math>P_y=Q_x</math> וכאשר <math>P,Q</math> בעלות נגזרות רציפות.אזי המד"ר מדוייקת אם"ם <math>P_y=Q_x</math>
*הוכחה לפתרון המד"ר המדויקת:
*דוגמא: נפתור מצאו משוואה המתארת את המשוואה הפתרון למד"ר הבאה באופן סתום <math>(2x+6y)dx+(6x+3y^2)dy=0</math>.
**ראשית נוודא שמדובר במשוואה מדוייקת: <math>P_y=Q_x=6</math>.
**נבצע אינטגרציה <math>U=\int Pdx +c(y)= x^2+6xy +c(y)</math>.
**לכן <math>c'(y)=Q-6x=3y^2</math>.
**לכן <math>c(y)=y^3</math> וסה"כ <math>U(x,y)=x^2+6xy+y^3</math>.
**לכן הפתרון למד"ר הוא נתון באופן סתום ע"י <math>x^2+6xy+y^3=C</math>.
***<math>c'(y)=y-x+x=y</math>.
***<math>c(y)=\frac{y^2}{2}</math>.
***סה"כ הפתרון למד"ר הוא נתון באופן סתום ע"י <math>U(x,y)=-\frac{1}{x}-yx+\frac{y^2}{2}=C</math>.***אפשר באמצעות השלמה לריבוע לבודד את y.
==הרצאה 4 משפט הקיום והיחידות==
*מציאת פתרון למד"ר <math>y'=f(x,y)</math> המקיימת <math>y(x_0)=y_0</math>
===המשוואה האינטגרלית===
*בעיית הקושי <math>y'=f(x,y)</math> עם <math>y(x_0)=y_0</math> שקולה למשוואה <math>y(x)=y_0+\int_{x_0}^xf(t,y(t))dt</math>.
**בכיוון אחד - נניח כי המשוואה הדיפרנציאלית ותנאי ההתחלה נתונים.
***אזי <math>\int_{x_0}^x y'(t)dt=\int_{x_0}^xf(t,y(t))dt</math>.
***לכן <math>y(x)-y(x_0)=\int_{x_0}^xf(t,y(t))dt</math>.
***ולפי תנאי ההתחלה נקבל כי <math>y(x)-y_0=\int_{x_0}^xf(t,y(t))dt</math>.
**בכיוון שני, נניח כי המשוואה האינטגרלית נתונה.
***נגזור את שני הצדדים ונקבל את המשוואה הדיפרנציאלית (נגזרת של פונקצית שטח של פונקציה רציפה).
***נציב במשוואה האינטגרלית את <math>x_0</math> ונקבל <math>y(x_0)=y_0+\int_{x_0}^{x_0}f(t,y(t))dt=y_0</math>.
*מאוחר יותר נוכיח כי סדרת הפונקציות מתכנסת לפתרון של המד"ר.
**<math>\varphi_1=y_0+\int_{x_0}^x(-ry_0)dt=y_0+y_0(-r(x-x_0))</math>
**<math>\varphi_2=y_0+\int_{x_0}^x\left(-r)\cdot(y_0-r\cdot y_0(t-x_0)\right)dt=y_0+y_0(-r(x-x_0))+y_0\frac{(-r(x-x_0))^2}{2}</math>
**<math>\varphi_3=y_0+\int_{x_0}^x\varphi_2dt=y_0+y_0(-r(x-x_0))+y_0\frac{(-r(x-x_0))^2}{2}+y_0\frac{(-r(x-x_0))^23}{3!}</math>
**נמשיך כך, ונקבל סדרת פונקציות המתכנסת ל<math>\varphi_n(x)\to y(x)=y_0e^{-r(x-x_0)}</math>
**אם נתון תנאי ההתחלה <math>y(0)=C</math> נקבל בדיוק את הפתרון <math>y=Ce^{-rx}</math>.
===ניסוח משפט הקיום והיחידות===
*תהי <math>f(x,y)</math> רציפה ובעלת נגזרת <math>f_y</math> רציפה במלבן הסגור <math>|x-x_0|\leq a, |y-y_0|\leq b</math>.
*נביט בבעיית הקושי <math>y'=f(x,y)</math>, עם תנאי ההתחלה <math>y(x_0)=y_0</math>
*נבחר <math>M</math> חסם כך ש <math>|f(x,y)|<M</math> במלבן הנתון, ונסמן <math>a'=\min\{a,\frac{b}{M}\}</math>.
*הערות:
*שימו לב שהמשפט מבטיח פתרון בתחום מצומצם.
**אכן ראינו מד"ר שהייתה מוגדרת ורציפה בכל הממשיים, אך לא היה פתרון שמוגדר בכל הממשיים(<math>y'=(x+y)^2</math>).
**לכל נקודה יש פתרון מסביבה, גם אם אין פתרון שמוגדר בכל מקום.
*שימו לב שאם מצאנו פתרון בצורה כלשהי, אנחנו יודעים שהוא יחיד בזכות המשפט (לפחות בסביבה מסויימת).
==הרצאה 4 =הוכחת משפט הקיום והיחידות=====המשוואה האינטגרלית===*בעיית נוכיח שסדרת הפונקציות בשיטת פיקרד מתכנסת לפתרון לבעיית הקושי <math>y'=f(x,y)</math> עם <math>y(x_0)=y_0</math> שקולה למשוואה <math>y(x)=y_0+\int_{x_0}^xf(t,y(t))dt</math>.**בכיוון אחד - נניח כי המשוואה הדיפרנציאלית ותנאי ההתחלה נתונים.***אזי הערה: נוכיח עבור <math>\int_{x_0}^x y'(t)dt=\int_{geq x_0}^xf(t,y(t))dt</math>.***לכן ההוכחות עבור <math>y(x)-y(x_0)=\int_{x_0}^xf(t,y(t))dt</math>.***ולפי תנאי ההתחלה נקבל כי <math>y(x)-y_0=\int_{x_0}^xf(t,y(t))dt</math>.**בכיוון שני, נניח כי המשוואה האינטגרלית נתונה.***נגזור את שני הצדדים ונקבל את המשוואה הדיפרנציאלית (נגזרת של פונקצית שטח של פונקציה רציפה).***נציב במשוואה האינטגרלית את <math>x_0</math> ונקבל <math>y(x_0)=y_0+\int_{x_0}^{x_0}f(t,y(t))dt=y_0</math>דומות.
*ראשית כעת, נשים לב לתכונה הבאה:
**כיוון ש<math>f_y</math> רציפה במלבן סגור היא חסומה נניח ע"י K.
**לפי משפט לגראנז' נקבל כי <math>|f(x,y_1)-f(x,y_2)|\leq K|y_1-y_2|</math>
*כעת נוכיח שסדרת הפונקציות מתכנסת (במ"ש):
**ראשית, נשים לב כי <math>\varphi_n-y_0=\varphi_n-\varphi_0=\varphi_n-\varphi_{n-1}+\varphi_{n-1}-\varphi_{n-2}+...+\varphi_1-\varphi_0</math>.
**לכן עלינו להוכיח כי הטור <math>\sum_{i=1}^n\infty\left(\varphi_i-\varphi_{i-1}\right)</math> מתכנס כאשר במ"ש (כי הסס"ח שלו היא <math>n\to\inftyvarphi_n</math>פחות קבוע).
**ראשית, <math>|\varphi_1-\varphi_0|=|y_0+\int_{x_0}^xf(t,y_0)dt-y_0|\leq M(x-x_0)</math>
**כעת <math>|\varphi_2-\varphi_1|\leq\int_{x_0}^x|f(t,\varphi_1)-f(t,\varphi_0)|dt\leq \int_{x_0}^xK|\varphi_1-\varphi_0|dt\leq KM\frac{(x-x_0)^2}{2}</math>
\sum_{i=1}^nK^{n-1}M\frac{(a')^n}{n!}
</math>
**זה טור מתכנס לפי מבחן המנה, ולפי וכן לפי מבחן הM של קושי, הטור המקורי מתכנס במידה שווה.
**הערה: כיוון ש<math>\left|f(x,\varphi_n(x))-f(x,\varphi_{n-1}(x))\right|\leq K|\varphi_n(x)-\varphi_{n-1}(x)|</math> אזי גם הסדרה <math>f(x,\varphi_n(x))</math> מתכנסת במ"ש באופן דומה.
**הערה: האינטגרל של הסדרה שואף לאינטגרל של פונקצית הגבול בזכות ההתכנסות במ"ש.
===הוכחת היחידות===
*טענת עזר - תהי <math>g</math> חסומה כך שלכל <math>x\geq x_0</math> בקטע <math>|x-x_0|\leq a</math> מתקיים כי <math>|g|\leq K\int_{x_0}^x|g(t)|dt</math> אזי <math>g=0</math> לכל <math>x\geq x_0</math> בקטע.
**<math>|g|\leq M</math>.
**<math>|g|\leq K\int_{x_0}^x|g|dt\leq KM(x-x_0)</math>.
**<math>|g|\leq K\int_{x_0}^x|g|dt\leq K\int_{x_0}^x KM(t-x_0)dt=K^2M\frac{(x-x_0)^2}{2}</math>.
**נמשיך כך ונקבל שלכל n מתקיים כי <math>|g|\leq K^nM\frac{(x-x_0)^n}{n!}</math>.
**לכן <math>|g|\leq K^n M\frac{a^n}{n!}\to 0</math>.
*יהיו שני פתרונות <math>y_1,y_2</math> לבעיית הקושי, נוכיח כי <math>y_1=y_2</math>:
**<math>|y_2-y_1|=\left|\int_{x_0}^x(f(t,y_1y_2)-f(t,y_2y_1))dt\right|\leq \int_{x_0}^x|f(t,y_1y_2)-f(t,y_2y_1)|dt\leq K\int_{x_0}^x|y_2-y_1|dt</math>.
**לכן לפי טענת העזר, <math>y_1=y_2</math>.
==הרצאה 5 מד"ר מסדר גבוה (ובפרט סדר שני), מד"ר לינארית מסדר גבוה==
**הכוח הפועל על המסה הוא <math>-kX</math>.
**לכן לפי החוק השני של ניוטון <math>mX''=-kX</math>.
*דוגמא:
**נביט בסירה במים המחוברת בקפיץ למזח.
**מלבד הכוח שהקפיץ מפעיל, המים מתנגדים לסירה באופן פרופורציוני למהירות שלה.
**<math>mX''=-kX-dX'</math>
**היחס בין קבוע הקפיץ לקבוע התנגדות המים ישפיע על התנועה - האם הסירה תתקדם בכיוון אחד, או תעשה תנועה מחזורית (בכל מקרה היא תאט).
*דוגמא:
**מסה מחוברת לקפיץ עם חיכוך
*דוגמא:
**מסה תלוייה על קפיץ במאונך עם או בלי התנגדות אוויר ועם השפעת כוח המשיכה (לא הומוגני)
**נביט בפונקצית המהירות <math>V=X'</math> ונקבל את המשוואה <math>mV'=C</math> מסדר ראשון.
====מדהורדת סדר למד"ר מסדר גבוה שני ללא x====*אם x אינו מופיע במשוואה תהי מד"ר מהצורה <math>y''=f(y',y)</math>.*ראשית נחפש פונקציה של <math>p</math> המקיימת את המד"ר מסדר ראשון <math>p'(t)p(t)=f(p(t),t)</math> **נהוג לרשום את שם המשתנה כאן y כך שיתקיים ולא t, אך אני לא עושה את זה כעת על מנת למנוע בלבול מיותר.*כעת נחפש פונקציה y המקיימת את המד"ר עבור p שמצאנו <math>y'=p(y)</math>*פונקציה כזו תקיים כי <math>y''=p'(y)y'=p'(y)p(y)=f(p(y),y)=f(y',y)</math>*כלומר היא מהווה פתרון למד"ר.
**נחזור לדוגמא של מסה המחוברת לקפיץ, ולצורך הנוחות נחליף את פונקצית המיקום X בפונקציה y (המשתנה ישאר t).
**נניח כי המסה היא חלק מקבוע הקפיץ ונביט במשוואה <math>y''=-ky</math>.
**נחפש פונקציה p של y המקיימת <math>y'=p(y)</math>.**לכן <math>y''=p'(y)y'=p'\cdot p</math>.**לכן אנחנו רוצים למצוא p פונקציה של y המקיימת את המשוואה <math>pp'=-ky</math>.
***זו משוואה פרידה <math>pdp=-kydy</math> ולכן <math>\frac{p^2}{2}=-\frac{ky^2}{2}+C</math>.
***לכן <math>p(y)=\pm\sqrt{C-ky^2}</math>.
**שימו לב שישנם שני קבועים בפתרון. זה הגיוני, כי אנו צריכים שני תנאי התחלה - מיקום המסה, והמהירות שלה.
=====דוגמא - מהירות מילוט=====
*גוף בעל מסה <math>m</math> נזרק מכדור הארץ כלפי מעלה במהירות <math>v_0</math>, נסמן את מרחק הגוף ממרכז כדור הארץ ב<math>r</math>.
**מצאו את פונקצית מהירות הגוף ביחס לגובה שלו <math>v(r)</math>.
**מהי מהירות המילוט של הגוף? כלומר עבור איזו מהירות התחלתית מתקיים כי <math>r(t)\to \infty</math> כאשר <math>t\to \infty</math>?
*נסמן את מסת כדור הארץ ב<math>m_e</math>, את רדיוס כדור הארץ ב<math>R_e</math>, את קבוע הכבידה האוניברסאלי ב<math>G</math> ואת תאוצת הנפילה בכדור הארץ ב<math>g</math>
**ראשית נשים לב כי כוח המשיכה של כדור הארץ המופעל על מסה <math>m</math> הוא בקירוב <math>mg=\frac{Gm_e m}{R_e^2}</math> כלומר <math>g=\frac{Gm_e}{R_e^2}</math> ולכן <math>gR_e^2 = Gm_e</math>
**המשוואה המתארת את תנועת הגוף היא:
***<math>mr''=-\frac{Gm_e m}{r^2}</math> כלומר <math>r''=-\frac{Gm_e}{r^2}=-\frac{gR_e^2}{r^2}</math>
**זו משוואה מסדר שני שחסר בה המשתנה <math>t</math>
**נחפש <math>p</math> עבורה <math>p(r)=r'</math> ולכן <math>pp'=r''</math>
***<math>pp'=-\frac{gR_e^2}{r^2}</math>
***נעשה אינטגרציה למד"ר הפרידה שקיבלנו ונקבל
***<math>\frac{p^2}{2}=\frac{gR_e^2}{r}+C</math>
***לכן <math>p(r)=\pm\sqrt{C+\frac{2gR_e^2}{r}}</math>
*כיוון שהמהירות ההתחלתית היא חיובית נקבל כי
**<math>r'=\sqrt{C+\frac{2gR_e^2}{r}}</math>
*על מנת למצוא את הקבוע, נציב את תנאי ההתחלה:
**הגובה הראשוני הוא <math>r=R_e</math> ובו המהירות היא <math>v_0</math>
**<math>v_0=\sqrt{C+2gR_e}</math>
**<math>C=v_0^2-2gR_e</math>
*הערה: ניתן לפתור את המד"ר הזו על מנת למצוא את הגובה כפונקציה של הזמן, אך לא התבקשנו לעשות כן.
*סה"כ נקבל כי <math>v(r)=\sqrt{v_0^2-2gR_e+\frac{2gR_e^2}{r}}</math>
*מהירות המילוט היא המהירות ההתחלתית הנמוכה ביותר המבטיחה כי הגוף לא יגיע למהירות אפס.
*לכן מהירות המילוט מקיימת כי <math>v_0^2 = 2gR_e</math> ולכן <math>v_0 =\sqrt{2gR_e}</math>
**לכל מהירות נמוכה יותר הביטוי בתוך השורש מתחיל מ<math>v_0</math> ושואף למספר שלילי (בהנחת השלילה ש <math>r\to \infty</math>), ולכן יגיע לאפס. במהירות אפס החפץ לא ימשיך לנוע.
**לכל מהירות התחלתית גבוהה יותר, המהירות גדולה יותר מערך חיובי קבוע, ולכן <math>r\to\infty</math>.
**אם המהירות ההתחלתית היא בדיוק מהירות המילוט, ניתן לפתור את המד"ר בקלות ולראות כי <math>r\to\infty</math>.
===מד"ר לינארית===
*משפט קיום ויחידות: אם <math>a_i(x),f(x)</math> רציפות בקטע <math>I</math> ויהי <math>x_0\in I</math>, אזי קיים פתרון יחיד בקטע <math>I</math> לבעיית הקושי.
*נגדיר את אופרטור הגזירה <math>D</math> על מרחב הפונקציות הגזירות אינסוף פעמים.
*<math>a(x)D</math> גם הוא אופרטור לינארי
*לכן ניתן לכתוב מד"ר לינארית כ <math>Ty=f(x)</math> כאשר <math>T=D^n+\sum_{k=1}^{n-1} a_k(x)\cdot D^k + I </math> אופרטור לינארי.
====מד"ר לינארית הומוגנית====
*אוסף הפתרונות של מד"ר לינארית הומוגנית הוא תת מרחב וקטורי.
**פונקצית האפס מקיימת את המשוואה.**אם זה הרי הגרעין של האופרטור <math>y_1,y_2T</math> פתרונות, ו<math>c\in\mathbb{R}</math> קבוע אזי קל לראות על ידי הצבה ישירה שגם <math>y_1+cy_2</math> הוא פתרון.המתואר לעיל
*אם <math>W(x_0)=0</math> עבור <math>x_0\in I</math> כלשהו עבור <math>y_1,...,y_n</math> '''פתרונות של מד"ר לינארית הומוגנית'''עם מקדמים רציפים בקטע <math>I</math>, אזי הפתרונות ת"ל ו<math>W(x)\equiv 0</math>.**כיוון ש<math>W(x_0)=0</math> קיים פתרון לא טריוויאלי למערכת כך שלכל <math>0\leq k\leq n-1</math> מתקיים כי <math>c_1y_1^{(k)}(x_0)+...+c_ny_n^{(n-1k)}(x_0)=0</math>.
**נביט בפונקציה <math>g(x)=c_1y_1(x)+...+c_ny_n(x)</math>, לפי לינאריות גם <math>g(x)</math> פתרון של המד"ר.
**כיוון שלכל <math>0\leq k\leq n-1</math> מתקיים כי <math>g^{(k)}(x_0)=0</math> ולפי יחידות הפתרון, נובע כי <math>g(x)\equiv 0</math> (הרי פונקצית האפס היא פתרון שמקיים את אותם תנאיי ההתחלה).
\lambda_1^{n-1}& \cdots & \lambda_n^{n-1}
\end{pmatrix}\right|</math>
**זו מטריצת ונדרמונדה ונדרמונד ולכן <math>W(x)=e^{(\lambda_1+...+\lambda_n)x}\prod_{i<j}(\lambda_j-\lambda_i)</math>
**לכן הפונקציות בת"ל אם ורק אם כל הקבועים שונים זה מזה <math>\lambda_i\neq\lambda_j</math>
*הוכחה לחישוב הדטרמיננטה של מטריצת ונדרמונד:
:<math>
\left|\begin{pmatrix}
1 & 1 &\cdots & 1 \\
\lambda_1 & \lambda_2 &\cdots & \lambda_n\\
\vdots & && \vdots \\
\lambda_1^{n-2}&\lambda_2^{n-2}&\cdots&\lambda_n^{n-2}\\
\lambda_1^{n-1}& \lambda_2^{n-1}&\cdots & \lambda_n^{n-1}
\end{pmatrix}\right|=
</math>
:נבצע את פעולות השורה<math>R_n-\lambda_1 R_{n-1}\\R_{n-1}-\lambda_1 R_{n-2}\\\vdots\\R_2-\lambda_1 R_1</math>
:<math>=\left|\begin{pmatrix}
1 & 1 & \cdots & 1 \\
0&\lambda_2-\lambda_1&\cdots&\lambda_n-\lambda_1\\
\vdots & && \vdots \\
0&\lambda_2^{n-3}(\lambda_2-\lambda_1)&\cdots&\lambda_n^{n-3}(\lambda_n-\lambda_1)\\
0&\lambda_2^{n-2}(\lambda_2-\lambda_1)& \cdots & \lambda_n^{n-2}(\lambda_n-\lambda_1)
\end{pmatrix}\right|=
(\lambda_2-\lambda_1)\cdots(\lambda_n-\lambda_1)\cdot
\left|\begin{pmatrix}
1 & 1 &\cdots & 1 \\
\lambda_2 & \lambda_3 &\cdots & \lambda_n\\
\vdots & && \vdots \\
\lambda_2^{n-2}&\lambda_3^{n-2}&\cdots&\lambda_n^{n-2}\\
\lambda_2^{n-1}& \lambda_3^{n-1}&\cdots & \lambda_n^{n-1}
\end{pmatrix}\right|
</math>
:כאשר המעבר הוא חישוב דטרמיננטה לפי העמודה הראשונה
:ומכאן סיימנו באינדוקציה
==הרצאה 6 מד"ר לינארית עם מקדמים קבועים==
ראשית נציג גישה אחת לנושא, ומאוחר יותר נציג גרסא מעודכנת (2022) המבוססות יותר על אופרטורים.
===פולינום אופייני===
*כעת נטפל במקרה בו שורש חוזר על עצמו:
**ראשית, נביט באופרטור הלינארי <math>D=\frac{d}{dx}</math> ששולח פונקציה לנגזרת שלה, ונסמן ב<math>I</math> את אופרטור הזהות.**למשל המד"ר <math>y''-2y+y=0</math> ניתנת להצגה כ<math>\left(\frac{d}{dx}D^2-2\frac{d}{dx}2D+I\right)y=0</math>.**לכן <math>\left(\frac{d}{dx}D-I\right)\left(\frac{d}{dx}D-I\right)y=0</math>.
**הפולינום האופייני של המד"ר הוא <math>(x-1)^2=0</math> ולכן <math>y=e^x</math> הוא פתרון.
**כעת, נראה כי גם <math>xe^x</math> הוא פתרון של המד"ר.
***<math>\left(\frac{d}{dx}D-I\right)\left(\frac{d}{dx}D-I\right)xe^x=\left(\frac{d}{dx}D-I\right)(e^x+xe^x-xe^x)=0</math>
**באופן דומה אפשר להוכיח שאם ריבוי השורש הוא <math>n</math> אזי לכל <math>0\leq k \leq n-1</math> הביטוי <math>x^ke^{\lambda x}</math> הוא פתרון.
===סיכום מציאת פתרון כללי למד"ר הומוגנית עם מקדמים קבועים===
*לכל שורש מרוכב <math>a+bi</math> מריבוי <math>n</math> (ידוע שגם הצמוד שלו שורש מאותו ריבוי) מתאימים הפתרונות <math>e^{ax}\cos(bx),e^{ax}\sin(bx),xe^{ax}\cos(bx),xe^{ax}\sin(bx),...,x^{n-1}e^{ax}\cos(bx),x^{n-1}e^{ax}\sin(bx)</math>
*סה"כ מצאנו למד"ר מסדר n בדיוק n פתרונות.
*הפתרונות הללו בת"ל (ללא הוכחה), ולכן הפתרון הכללי הוא צירוף לינארית לינארי שלהם.**נוכיח שהפתרונות בת"ל (מעל המרוכבים).**<math>P_1e^{\lambda_1 x}+...+P_ne^{\lambda_n x} \equiv 0</math>.**נניח ש<math>|\lambda_i|\leq|\lambda_n|</math>, נחלק ב<math>e^{\lambda_n x}</math>.**נציב <math>x=t\overline{\lambda_n}</math> ונשאיף את <math>t\to\infty</math>.**נקבל כי הפולינום המקדם של האקספוננט הגדול ביותר חייב להיות אפס.**לכן באינדוקציה כל הפולינומים חייבים להיות אפס, ולכן כל אחד מהקבועים חייב להיות אפס.**כיוון שהפתרונות בת"ל מעל המרוכבים, אפשר ליצור איתם כל תנאי התחלה, ולקבל פונקציות ממשיות שפותרות אותו.
***<math>y''(0)=-2+2c_3=0</math> ולכן <math>c_3=1</math>.
**סה"כ הפתרון הוא <math>y=e^{-x}(x+x^2)</math>.
===גישה מבוססת אופרטורים===
*נציג את המד"ר הלינארית עם מקדמים קבועים באמצעות אופרטור הגזירה:
*<math>y^{(n)}+a_{n-1}y^{(n-1)}+\cdots+a_0y = (D^n+a_{n-1}D^{n-1}+\cdots+a_0 I)y=Ty</math>
*נגדיר את הפולינום האופייני <math>p(x)=x^n+a_{n-1}x^{n-1}+\cdots+a_0</math>
*סה"כ האופרטור של המד"ר הוא <math>T=p(D)</math>
*נפרק את הפולינום האופייני לגורמים לינאריים מעל המרוכבים
*<math>p(x)=(x-\lambda_1)(x-\lambda_2)\cdots(x-\lambda_n)</math>
*ונקבל כי <math>T=p(D)=(D-\lambda_1 I)\cdots (D-\lambda_n I)</math>
**שימו לב כי מותר לפתוח סוגריים באופן טבעי ואפשר להחליף בין סדר הגורמים כיוון ש <math>D,\lambda I</math> אופרטורים מתחלפים.
*כיוון שמותר להחליף את סדר הגורמים נובע כי אם <math>\lambda</math> שורש של הפולינום האופייני מריבוי <math>k</math> אזי
**<math>\ker\left((D-\lambda I)^k\right)\subseteq \ker T</math>
בטקסט לעיל, למדנו איך למצוא בסיס לגרעין הזה.
==הרצאה 7 מציאת פתרון פרטי למד"ר לינארית לא הומוגנית==
*אם <math>f(x)=P_m(x)</math> פולינום מדרגה m:
**<math>0</math> '''אינו''' שורש של הפולינום האופייני, ננחש <math>y_p=Q_m(x)</math> פולינום מדרגה m.
**אם <math>0</math> שורש של הפולינום האופייני מדרגה מריבוי k ננחש <math>y_p=x^kQ_m(x)</math>.
*אם <math>f(x)=e^{ax}P_m(x)</math>:
**אם <math>a</math> '''אינו''' שורש של הפולינום האופייני ננחש <math>y_p=e^{ax}Q_m(x)</math>.
**אם <math>a</math> שורש של הפולינום האופייני מדרגה מריבוי k ננחש <math>y_p=x^ke^{ax}Q_m(x)</math>.
**לכן הפתרון הפרטי הוא <math>y_p=x^2-4x+6</math>.
**סה"כ הפתרון הכללי הוא <math>c_1e^{-x}+c_2xe^{-x}+x^2-4x+6</math>.
===וריאצית מקדמים יחד עם שיטת קרמר למד"ר לינארית===
*טענה - עבור פונקציות <math>c_1(x),...,c_n(x)</math> המקיימות את מערכת המשוואות <math>\begin{cases}
c_1'y_1+...+c_n'y_n=0 \\
c_1'y_1'+...+c_n'y_n'=0 \\
c_1'y_1^{(n-2)} +...+c_n'y_n^{(n-2)}=0\\
c_1'y_1^{(n-1)}+...+c_n'y_n^{(n-1)}=f(x)
\end{cases}</math> מתקיים כי <math>y_p=c_1(x)y_1+...+c_n(x)y_n</math> הוא פתרון פרטי של המד"ר.* *הוכחה:**עבור פשטות הרישום, נבצע את ההוכחה עבור n=3, אמנם ההוכחה הכללית דומה לחלוטין.**<math>y_p'=c_1'y_1+c_2'y_2\cdots+c_3c_n'y_3y_n+c_1y_1'+c_2y_2'\cdots+c_3y_3c_ny_n'=c_1y_1'+c_2y_2'\cdots+c_3y_3c_ny_n'</math>. (לפי המשוואה הראשונה.)**באופן דומה <math>y_p''=c_1y_1''+c_2y_2''\cdots+c_3y_3c_ny_n''</math>. (לפי המשוואה השנייה.)**נמשיך כך עד שנקבל <math>y_p^{(n-1)} = c_1y_1^{(n-1)}+\cdots +c_ny_n^{(n-1)}</math>**כעת נגזור ונקבל <math>y_p'''^{(n)}=f(x)+c_1y_1'''^{(n)}+c_2y_2'''\cdots+c_3y_3'''c_ny_n^{(n)}</math>, לפי המשוואה האחרונה.
**נציב במד"ר המקורית:
***<math>y_p'''^{(n)}+a_2a_{n-1}(x)y_p''^{(n-1)}+\cdots +a_1(x)y_p'+a_0(x)y_p=f(x)+c_1(y_1^{(n)}+\cdots+a_0(x)y_1)+\cdots+c_n(y_n^{(n)}+\cdots+a_0(x)y_n)</math>***כיוון ש<math>+c_1(y_1,...,y_n</math> פתרונות למד"ר ההומוגנית הביטויים בסוגריים מתאפסים וסה"כ קיבלנו כי אכן <math>y_p'''+a_2(x)y_1y_p''+a_1(x)y_1y_p'+a_0(x)y_1y_p=f(x)+</math>. *נכתוב '''שוב''' את ההוכחה, בעזרת סימן הסכימה (עשוי להיות נוח יותר או פחות):**ראשית, ניתן להוכיח באינדוקציה כי לכל <math>+c_20\leq m\leq n-1</math> מתקיים כי **<math>D^m y_p = D^m \sum_{k=1}^n c_k(y_2'''+a_2x)y_k = \sum_{k=1}^n c_k(x)y_2''+a_1D^m y_k</math>**כעת בעזרת המשוואה האחרונה נקבל כי**<math>D^n y_p = D D^{n-1}y_p = D\sum_{k=1}^nc_k(x)y_2D^{n-1}y_k=\sum_{k=1}^n c'_k(x)D^{n-1}y_k +a_0\sum_{k=1}^nc_k(x)y_2D^ny_k=f(x)+\sum_{k=1}^nc_k(x)D^ny_k</math>**נציב במד"ר ונקבל**<math>Ty_p=D^ny_p +c_3\sum_{t=0}^{n-1}a_t(y_3'''+a_2x)D^ty_p=f(x)y_3''+a_1\sum_{k=1}^nc_k(x)y_3'D^ny_k +a_0\sum_{t=0}^{n-1}a_t(x)y_3\left(\sum_{k=1}^n c_k(x)D^t y_k\right)=</math>**כיוון ש<math>y_1,y_2,y_3</math> פתרונות למד"ר ההומוגנית הביטויים בסוגריים מתאפסים וסה"כ קיבלנו כי אכן <math>y_p'''+a_2=f(x)y_p''+a_1\sum_{k=1}^n c_k(x)y_p'\left(D^ny_k +a_0\sum_{t=0}^{n-1}a_t(x)y_pD^t y_k\right) =f(x)+0</math>.
\end{pmatrix}
\right|
}=sin^2(x)cos(x)
</math>
***לכן <math>c_1(x)=\int (-sin^3(x))dx = \int (1-cos^2(x))(-sin(x))dx=\{t=cos(x)\}=\int (1-t^2)dt=t-\frac{t^3}{3}=cos(x)-\frac{cos^3(x)}{3}</math>
*לכן <math>y_p=y_{p_1}+y_{p_2}</math> הוא פתרון פרטי למד"ר <math>y''+y=sin^2(x)</math> מתוך לינאריות.
==הרצאה 8 פתרון מד"ר באמצעות טורי טיילור ומערכות מד"ר==
===שימוש בטורי טיילור===
*ננחש שהפתרון הוא טור חזקות, ואם אכן יש פתרון כזה, נמצא את המקדמים.
*דוגמא מצאו את הפתרון הכללי למד"ר ההומוגנית <math>xy''-(x+2)y'+2y=0</math>.
*עבור <math>x\neq 0</math> מדובר במד"ר לינארית הומוגנית בעלת שני פתרונות בת"ל.
*<math>\sum_{n=2}^\infty n(n-1)a_nx^{n-1} -\sum_{n=1}^\infty na_nx^{n} - \sum_{n=1}^\infty 2na_nx^{n-1}+\sum_{n=0}^\infty 2a_nx^n=0</math>
*<math>\sum_{k=1}^\infty (k+1)ka_{k+1}x^{k} -\sum_{k=1}^\infty ka_kx^{k} - \sum_{k=0}^\infty 2(k+1)a_{k+1}x^{k}+\sum_{k=0}^\infty 2a_kx^k=0</math>
*<math>-2a_1+2a_0+\sum_{k=01}^\infty \left((k^2-k-2)a_{k+1}-(k-2)a_k\right)x^k=0</math>
*לכן:
**<math>a_0=a_1</math>
**לכל <math>k\geq 1</math> מתקיים <math>(k^2-k-2)a_{k+1}-(k-2)a_k=0</math>.
***עבור <math>k=2</math> מקבלים <math>0=0</math>.***עבור <math>k\neq 2</math> נחלק ב<math>k-2</math> ונקבל <math>(k+1)a_{k+1}=a_k</math>.
**<math>a_2=\frac{1}{2}a_1</math>
**<math>a_4=\frac{1}{4}a_3</math>
**<math>a_5=\frac{1}{5}a_4</math>
**וכן הלאה.
====מציאת פתרון פרטי====*דוגמא - מצאו את הפתרון הכללי למד"ר <math>xy''-(x+2)y'+2y=x^2e3e^x</math>.*ראשית נעביר את המד"ר לצורה סטנדרטית <math>y''-\frac{x+2}{x}y'+\frac{2}{x}y=x^2e^x</math>*נשתמש בשיטת וריאצית המקדמים על הפתרון למד"ר ההומוגנית יחד עם כלל קרמר.
**נחפש פתרון מהצורה <math>y_p=c_1(x)e^x+c_2(x)\left(\frac{1}{2}x^2+x+1\right)</math>.
**כעת <math>c_1'=\frac{\left|\begin{pmatrix} 0 & \frac{1}{2}x^2+x+1 \\ x^2e^x & x+1\end{pmatrix}\right|}{W(x)}=x^2+2x+2</math>
**לכן <math>c_2(x)=-2e^x</math>.
*סה"כ הפתרון הפרטי הינו <math>y_p=\left(\frac{1}{3}x^3+x^2+2x\right)e^x-2e^x\left(\frac{1}{2}x^2+x+1\right) = e^x\left(\frac{1}{3}x^3-2\right)</math>
*לכן הפתרון הכללי הינו <math>y=c_1e^x+c_2\left(\frac{1}{2}x^2+x+1\right)+e^x\left(\frac{1}{3}x^3-2\right)</math>
==הרצאה 9 מערכות מד"ר==
===מערכת מד"ר לינארית מסדר ראשון עם מקדמים קבועים===
*לעיתים יש לנו מד"ר העוסקות במספר פונקציות שונות.
*נניח שיש לנו סיר מים מתבשל על הגז.
*A היא מסת המים בסיר, וB היא מסת המים שהתאדו אל המכסה.
*נניח שקצב התאדות המים מהסיר אל המכסה הוא <math>\alpha\cdot A</math> וקצב התעבות המים מהמכסה בחזרה לסיר הוא <math>\beta\cdot B</math>.
*לכן <math>\begin{cases}A'=\beta B - \alpha A \\ B' = \alpha A - \beta B\end{cases}</math>
*נסמן את שתי הפונקציות ב<math>y_1,y_2</math> ונניח כי <math>\alpha =1, \beta=2</math>.
*נקבל את המערכת <math>\vec{y}'=A\vec{y}</math> כלומר <math>\begin{pmatrix}y_1'\\y_2'\end{pmatrix}=\begin{pmatrix}-1 & 2 \\ 1 &-2\end{pmatrix}\begin{pmatrix}y_1\\y_2\end{pmatrix} </math>
*נראה כיצד לכסון המטריצה A יעזור לנו לפתור את המערכת.
*במקרה בו A אינה לכסינה לא נטפל, אך אפשר לפתור אותו באופן כללי.
*עבור ו"ע מתקיים כי <math>A\vec{v}=\lambda \vec{v}</math>.
*כיוון שהוקטור <math>\vec{v}</math> הוא וקטור קבועים, <math>\left(\vec{v}e^{\lambda x}\right)'=\lambda\vec{v}e^{\lambda x} = A\left(\vec{v}e^{\lambda x}\right)</math>.
*כלומר, <math>\vec{y}=\vec{v}e^{\lambda x}</math> הוא פתרון למערכת.
*בחזרה לדוגמא:
**הע"ע של <math>\begin{pmatrix}-1 & 2 \\ 1 &-2\end{pmatrix}</math> הם <math>0,-3</math>.
**הו"ע המתאימים הם <math>\begin{pmatrix}2\\1\end{pmatrix},\begin{pmatrix}1\\-1\end{pmatrix}</math>
**הפתרון הכללי הוא <math>\vec{y}=c_1\begin{pmatrix}2\\1\end{pmatrix}e^0+c_2\begin{pmatrix}1\\-1\end{pmatrix}e^{-3x}</math>
**כלומר <math>y_1=2c_1+c_2e^{-3x}</math> ו<math>y_2=c_1-c_2e^{-3x}</math>
*שימו לב שככל שעובר הזמן היחס בין המים בסיר למים על המכסה שואף להיות קבוע.
*שימו לב ש<math>c_1=\frac{y_1(0)+y_2(0)}{3}</math>, זה הגיוני כיוון שמסת המים אינה משתנה בתהליך.
====שתי מסות על קפיץ - מערכת מד"ר מסדר שני====
*נביט בשתי מסות המחוברות לשני צידי קפיץ.
*נניח כי <math>y_2<y_1</math> מודדות את מיקום המסות ביחס לנקודת האפס שלהן, וצד ימין הוא הכיוון החיובי בשתיהן.
*נניח כי כאשר כל אחת מהמסות במקום אפס, אזי הקפיץ במנוחה.
*נניח כי המסות זהות בגודלן, ושוות אחד.
*לכן מתקבלת מערכת המד"ר <math>\begin{cases}y_1''=-k(y_1-y_2) \\ y_2''=k(y_1-y_2)\end{cases}</math>
*שימו לב שכאשר הקפיץ מתוח הוא מושך את שתי המסות למרכז, כלומר את המסה הראשונה (הימנית) הוא מושך שמאלה (בכיוון השלילי), ואת המסה השנייה (השמאלית) הוא מושך ימינה (בכיוון החיובי)
*נסמן <math>A=\begin{pmatrix}-k & k \\ k & -k\end{pmatrix}</math>, ולכן <math>\vec{y}''=A\vec{y}</math>.
*הע"ע של A הינם <math>0,-2k</math>.
*עבור הו"ע <math>\vec{v}=\begin{pmatrix}1\\1\end{pmatrix}</math> המתאים לע"ע <math>0</math> מתקיים כי <math>A\vec{v}=0</math>.
**לכן אם נבחר <math>f(t)</math> כך ש<math>f''=0</math>, ונבחר <math>\vec{y}=\vec{v}f(t)</math> אזי נקבל <math>\vec{y}''=0=A\vec{v}f(t)=A\vec{y}</math>.
**כלומר <math>\vec{y}=\vec{v}(c_1t+c_2)</math> הוא פתרון למערכת.
*עבור הו"ע <math>\vec{v}=\begin{pmatrix}1\\-1\end{pmatrix}</math> המתאים לע"ע <math>-2k</math> מתקיים כי <math>A\vec{v}=-2k\vec{v}</math>.
**לכן אם נבחר <math>f(t)</math> כך ש<math>f''=-2kf</math> ונבחר <math>\vec{y}=\vec{v}f(t)</math> אזי נקבל <math>\vec{y}''=-2k\vec{v}f(t)=A\vec{v}f(t)=A\vec{y}</math>.
**לכן <math>\vec{y}=\left(c_3cos\left(\sqrt{2k}t\right)+c_4sin\left(\sqrt{2k}t\right)\right)\vec{v}</math> הוא פתרון למשוואה.
*ביחד קיבלנו פתרון כללי <math>\vec{y}=(c_1t+c_2)\begin{pmatrix}1\\1\end{pmatrix}+\left(c_3cos\left(\sqrt{2k}t\right)+c_4sin\left(\sqrt{2k}t\right)\right)\begin{pmatrix}1\\-1\end{pmatrix} </math>
*תנאי ההתחלה הם המיקומים והמהירויות של כל אחת מהמסות.
====קשר בין מד"ר מסדר גבוה למערכת מד"ר מסדר ראשון====
*נביט במד"ר <math>f(x,y,y',...,y^{(n)})=0</math>.
*נסמן <math>y_1=y,y_2=y',...,y_n=y^{(n-1)}</math>.
*לכן המד"ר שקולה למערכת מסדר ראשון <math>\begin{cases}y_1'=y_2 \\ \vdots \\ y_{n-1}'=y_n \\ f(x,y_1,...,y_n,y_n')=0\end{cases}</math>.
*בפרט, המד"ר הלינארית <math>y^{(n)}+a_{n-1}y^{(n-1)}+...+a_0y=0</math> שקולה למערכת <math>\begin{cases}y_1'=y_2 \\ \vdots \\ y_{n-1}'=y_n \\ y_n'=-a_{n-1}y_{n}-...-a_0y_1\end{cases}</math>
*בכתיב מטריצות קיבלנו את המערכת <math>\vec{y}'=A\vec{y}</math> כאשר:
**<math>\vec{y}=\begin{pmatrix}y_1\\ \vdots \\ y_n\end{pmatrix}</math>
**<math>A=\begin{pmatrix}
& 1 \\
& & 1 \\
& & & \ddots \\
& & & & 1\\
-a_0 & -a_1 & -a_2 & \cdots & -a_{n-1}
\end{pmatrix}</math>
*הפולינום האופייני של <math>A</math> הוא:
**<math>p_A(x)=\left|\begin{pmatrix}
x & -1 \\
& x & -1 \\
& & \ddots & \ddots \\
& & & x& -1\\
a_0 & a_1 & \cdots & a_{n-2} & x+a_{n-1}
\end{pmatrix}\right|</math>
**ניתן להוכיח באינדוקציה כי <math>p_A(x)=x^n+a_{n-1}x^{n-1}+...+a_0</math>, בדיוק הפולינום האופייני של המד"ר המקורית, לא במפתיע.
==הרצאה 10 התמרת לפלס==
*התמרת לפלס היא העתקה לינארית בין מרחבי פונקציות.
*עבור הפונקציה <math>y(t)</math> המוגדרת בקטע <math>[0,\infty)</math> נגדיר את התמרת הלפלס <math>F(s)=\mathcal{L}(y)=\int_0^\infty e^{-st}y(t)dt</math>.
*שימו לב שנהוג לסמן את הפונקציה לפני ההתמרה עם המשתנים x או t, ולאחר ההתמרה נהוג להתמש במשתנה s.
*אם מתקיים כי <math>|y(t)|\leq Me^{at}</math> אזי ההתמרה מתכנסת לכל <math>s>a</math>.
**<math>\left|\int_0^\infty e^{-st}y(t)dt\right|\leq \int_0^\infty\left|e^{-st}y(t)\right|dt\leq \int_0^\infty Me^{(a-s)t}dt=\left[M\frac{e^{(a-s)t}}{a-s}\right]_0^\infty</math>
**הביטוי האחרון מתכנס לכל <math>s>a</math>.
*נניח כי כל הפונקציות שאנו עוסקים בהן חסומות על ידי אקספוננט באופן דומה.
*דוגמא - נחשב את ההתמרה של הפונקציה <math>1</math>.
**<math>\mathcal{L}(1)=\int_0^\infty e^{-st}dt = \left[\frac{e^{-st}}{-s}\right]_0^\infty = \frac{1}{s}</math>
*בויקיפדיה ניתן למצוא [https://he.wikipedia.org/wiki/%D7%94%D7%AA%D7%9E%D7%A8%D7%AA_%D7%9C%D7%A4%D7%9C%D7%A1#%D7%98%D7%91%D7%9C%D7%AA_%D7%94%D7%AA%D7%9E%D7%A8%D7%95%D7%AA_%D7%9C%D7%A4%D7%9C%D7%A1 טבלה של התמרות לפלס שימושיות].
*שימו לב לשימוש בפונקצית המדרגה <math>u(t)=\begin{cases}1 & t\geq 0\\ 0 & t<0\end{cases}</math> שמאפסת את כל החלק השלילי של ציר הx.
**הפונקציה <math>u(t-a)</math> מאפסת את ציר הx בקטע <math>(-\infty,a)</math>.
===תכונות התמרת לפלס===
*יחידות:
**אם <math>y_1,y_2</math> רציפות, ו<math>\mathcal{L}(y_1)=\mathcal{L}(y_2)</math> אזי <math>y_1=y_2</math>.
**[http://ctr.maths.lu.se/media/MATC12/2013ht2013/uniqueness.pdf הוכחה]
*לינאריות:
**<math>\mathcal{L}(y_1+ay_2) = \mathcal{L}(y_1)+a\mathcal{L}(y_2)</math>
*התמרת הנגזרת הראשונה:
**<math>\mathcal{L}(y')=s\mathcal{L}(y)-y(0)</math>
*התמרת נגזרת כללית:
**<math>\mathcal{L}(y^{(n)})=s^n\mathcal{L}(y)-s^{n-1}y(0)-s^{n-2}y'(0)-...-y^{(n-1)}(0)</math>
*הזזה של המשתנה s:
**אם <math>F(s)=\mathcal{L}(y)</math> אזי <math>F(s-a)=\mathcal{L}(e^{at}y)</math>
*הזזה של המשתנה t:
**אם <math>F(s)=\mathcal{L}(y)</math> אזי <math>e^{-as}F(s)=\mathcal{L}(u(t-a)y(t-a))</math>
*תכונות נוספות:
**אם <math>F(s)=\mathcal{L}(y)</math> אזי <math>\mathcal{L}(ty)=-F'(s)</math>
**אם <math>F(s)=\mathcal{L}(y)</math> אזי <math>\mathcal{L}(ty')=-F(s)-sF'(s)</math>
**אם <math>F(s)=\mathcal{L}(y)</math> אזי <math>\mathcal{L}(ty'')=-2sF(s)-s^2F'(s)+y(0)</math>
*נוכיח חלק מהתכונות לעיל כעת ובהרצאה הבאה.
*נוכיח עבור y החסומה ע"י אקספוננט כי <math>\mathcal{L}(y')=sF(s)-y(0)</math>
**<math>\mathcal{L}(y')=\int_0^\infty e^{-st}y'(t)dt</math>
**נבצע אינטגרציה בחלקים
**<math>\int_0^\infty e^{-st}y'(t)dt=\left[e^{-st}y(t)\right]_0^\infty+s\int_0^\infty e^{-st}y(t)dt = -y(0)+sF(s)</math>
*כעת <math>\mathcal{L}(y'')=s\mathcal{L}(y')-y'(0) = s^2F(s)-sy(0)-y'(0)</math>.
*וכן הלאה, עבור נגזרות מסדר גבוה.
===דוגמאות===
*דוגמא - נמצא את ההתמרה של האקספוננט
*נציב בנוסחא <math>\mathcal{L}(y')=s\mathcal{L}(y)-y(0)</math> את <math>y=e^{ax}</math>
*<math>\mathcal{L}(ae^{ax})=s\mathcal{L}(e^{ax})-1</math>
*סה"כ נקבל כי <math>\mathcal{L}(e^{ax})=\frac{1}{s-a}</math>
*דוגמא - נמצא פתרון למד"ר <math>y'=ry</math>.
*נבצע התמרת לפלס:
**<math>0=\mathcal{L}(y'-ry)=sF(s)-y(0)-rF(s)</math>
**<math>F(s)=\frac{y(0)}{s-r}</math>
**לכן <math>y=y(0)e^{rt}</math>
*דוגמא - נמצא את ההתמרה של סינוס וקוסינוס
*נסמן <math>F(s)=\mathcal{L}(\sin(ax))</math>, <math>G(s)=\mathcal{L}(\cos(ax))</math>
*נציב בנוסחא <math>\mathcal{L}(y')=s\mathcal{L}(y)-y(0)</math>:
**נציב <math>y=\sin(ax)</math> ונקבל <math>\mathcal{L}(a\cos(ax))=s\mathcal{L}(\sin(ax))-0</math> כלומר <math>aG(s)=sF(s)</math>
**נציב <math>y=\cos(ax)</math> ונקבל <math>\mathcal{L}(-a\sin(ax))=s\mathcal{L}(\cos(ax))-1</math> כלומר <math>-aF(s)=sG(s)-1</math>
*נקבל סה"כ כי
**<math>\mathcal{L}(sin(ax))=F(s)=\frac{a}{s^2+a^2}</math>
**<math>\mathcal{L}(cos(ax))=G(s)=\frac{s}{s^2+a^2}</math>
==הרצאה 11 - המשך התמרת לפלס==
*נוכיח כי <math>\mathcal{L}(e^{at}y(t)) = F(s-a)</math>
**<math>\mathcal{L}(e^{at}y(t))=\int_0^\infty e^{-st}e^{at}y(t)dt = \int_0^\infty e^{-(s-a)t}y(t)dt=F(s-a)</math>
*נפתור את המד"ר <math>y''-2y'+2y=0</math> עם תנאי ההתחלה <math>y(0)=0,y'(0)=1</math>.
*שימו לב שכבר למדנו איך לפתור מד"ר זו - למצוא פתרון כללי ולהציב תנאי ההתחלה.
*התמרת לפלס עשוייה לחסוך לנו קצת זמן.
*נבצע התמרת לפלס:
**<math>s^2F(s)-sy(0)-y'(0)-2(sF(s)-y(0))+F(s)=0</math>
**<math>F(s)=\frac{1}{s^2-2s+2} = \frac{1}{(s-1)^2+1}</math>
*ידוע ש<math>G(s)=\frac{1}{s^2+1}</math> הינה ההתמרה של <math>sin(t)</math>.
*לכן <math>F(s)=G(s-1)</math> הינה ההתמרה של <math>e^tsin(t)</math>, וזהו פתרון המד"ר.
*נוכיח כי אם <math>F(s)=\mathcal{L}(y)</math> אזי <math>\mathcal{L}(ty)=-F'(s)</math>
*<math>F(s)=\int_0^\infty e^{-st}y(t)dt</math>.
*נגזור את שני הצדדים לפי <math>s</math> ונקבל כי
**<math>F'(s)=\frac{\partial}{\partial s} \int_0^\infty e^{-st}y(t)dt=\int_0^\infty -te^{-st}y(t)dt=-\mathcal{L}(ty)</math>
**את העובדה שגזרנו בתוך האינטגרל לא נצדיק כאן, היא נכונה עבור פונקציות שחסומות על ידי אקספוננט.
*לכן, <math>\mathcal{L}(ty') = -\frac{\partial}{\partial s}\mathcal{L}(y') = -\frac{\partial}{\partial s}(sF(s)-y(0)) = -F(s)-sF'(s)</math>
*כמו כן, <math>\mathcal{L}(ty'') = -\frac{\partial}{\partial s}\mathcal{L}(y'') = -\frac{\partial}{\partial s}(s^2F(s)-sy(0)-y'(0)) = -(2sF(s)+s^2F'(s)-y(0))</math>
*דוגמא - נחשב את <math>\mathcal{L}(t^n)</math>.
**ידוע כי <math>\mathcal{L}(1)=\frac{1}{s}</math>
**לכן <math>\mathcal{L}(t)=-\frac{\partial}{\partial s}\mathcal{L}(1)= \frac{1}{s^2}</math>
**לכן <math>\mathcal{L}(t^2)=-\frac{\partial}{\partial s}\mathcal{L}(t)= \frac{2}{s^3}</math>
**לכן <math>\mathcal{L}(t^3)=-\frac{\partial}{\partial s}\mathcal{L}(t^2)= \frac{3!}{s^4}</math>
**ובאופן כללי <math>\mathcal{L}(t^n)=\frac{n!}{s^{n+1}}</math>
===דוגמא===
*נפתור את המד"ר <math>xy''-(x+2)y'+2y=0</math>.
*נבצע התמרת לפלס:
**<math>\mathcal{L}(xy''-(x+2)y'+2y)=\mathcal{L}(xy'')-\mathcal{L}(xy')-2\mathcal{L}(y')+2\mathcal{L}(y)=</math>
**<math>=-2sF(s)-s^2F'(s)+y(0)+F(s)+sF'(s)-2sF(s)+2y(0)+2F(s)</math>
**לכן קבלנו את המשוואה <math>(s-s^2)F'(s)+(3-4s)F(s)=-3y(0)</math>
*קיבלנו מד"ר לינארית.
*לצורך הנוחות, נחליף זמנית את הסימון ונפתור את <math>y'+\frac{3-4x}{x-x^2}y=\frac{-3y_0}{x-x^2}</math>
**נסמן <math>P(x)=\frac{3-4x}{x-x^2}=\frac{3}{x}+\frac{1}{x-1}</math>, ו<math>Q(x)=\frac{-3y_0}{x-x^2}</math>
**לכן <math>e^{-\int P(x)}=\frac{1}{x^3(x-1)}</math>.
**כמו כן <math>\int Q(x)e^{\int P(x)} = \int \frac{-3y_0}{x-x^2}x^3(x-1) = \int 3y_0x^2=y_0x^3</math>
**סה"כ הפתרון למד"ר הלינארית הוא <math>y=\frac{1}{x^3(x-1)}\left(y_0x^3+C\right)=\frac{y_0}{x-1}+\frac{C}{x^3(x-1)}</math>
*נחזור לסימון התמרת הלפלס:
**<math>F(s)=\frac{y(0)}{s-1}+\frac{C}{s^3(s-1)}=\frac{y(0)+C}{s-1} - C\left(\frac{1}{s}+\frac{1}{s^2}+\frac{1}{s^3}\right)</math>
*נבצע התמרה הפוכה על מנת לקבל את הפתרון למשוואה המקורית:
**<math>y=\mathcal{L}^{-1}(F(s))=(y(0)+C)e^x - C(1+x+\frac{1}{2}x^2)</math>
===דוגמא===
*נמצא פתרון למד"ר <math>ty''+2y'+ty=0</math> המקיים <math>y(0)=1</math>.
**נבצע התמרת לפלס <math>-2sF(s)-s^2F'(s)+1+2sF(s)-2-F'(s)=0</math>.
**לכן <math>F'(s)=-\frac{1}{1+s^2}</math>
**לכן <math>\mathcal{L}(ty)=\frac{1}{1+s^2}</math>
**לכן <math>ty=sin(t)</math>
**לכן <math>y=\frac{sin(t)}{t}</math>
*הערות:
**הפונקציה שקיבלנו רציפה אם נגדיר אותה ב0 להיות 1, ואכן מקיימת את תנאי ההתחלה.
**מצאנו רק פתרון אחד, כיוון שלפתרון השני <math>\frac{cos(t)}{t}</math> אין התמרת לפלס (האינטגרל לא מתכנס באיזור 0).
==הרצאה 12 - הדלתא של דירק==
===הדלתא של דירק===
*נתחיל ונאמר כי ישנן מספר גישות אל הדלתא של דירק, אנחנו נציג גישה אחת שרלוונטית אלינו.
*הדלתא של דירק '''אינה פונקציה''', אלא מייצגת תהליך.
*למרות האמור, אנחנו נתייחס לתוצאה הסופית של התהליך, כאילו היה מדובר בפונקציה ממש.
*מטרה עיקרית: 'פונקצית הדלתא' מקיימת את התכונה <math>\int_{-\infty}^\infty f(x)\delta(x)dx=f(0)</math> לכל פונקציה <math>f(x)</math> הרציפה ב<math>0</math>.
*כמו כן, <math>\int_{-\infty}^\infty f(x)\delta(x-a)dx=\{t=x-a\}=\int_{-\infty}^\infty f(t+a)\delta(t)dt=f(a)</math> לכל פונקציה הרציפה בa.
*בצורה מדוייקת יותר, נביט בסדרת הפונקציות <math>\delta_n(x)=\begin{cases}n & 0\leq x \leq \frac{1}{n}\\ 0 & x< 0 \vee x>\frac{1}{n}\end{cases}</math>
*כאשר <math>n\to\infty</math> לכל <math>x\neq 0</math> מתקיים כי <math>\delta_n(x)\to 0</math> ועבור <math>x=0</math> מקבלים כי <math>\delta_n(x)\to \infty</math>.
*לכל <math>n</math> מתקיים כי <math>\int_{-\infty}^\infty \delta_n(x)dx=1</math>.
*עקרונית הסדרה מייצגת פונקציות בעלות שטח אחד, ההולך ומתרכז בנקודה אפס.
*עבור <math>f(x)</math> הרציפה בסביבה של <math>0</math> מתקיים כי:
**<math>\int_{-\infty}^\infty f(x)\delta_n(x)dx=\int_0^{\frac{1}{n}}nf(x)dx</math>
**לפי משפט ערך הממוצע האינטגרלי <math>\int_0^{\frac{1}{n}}nf(x)dx=nf(c_n)\cdot \frac{1}{n}=f(c_n)\to f(0)</math>
*נגדיר את <math>\int_{-\infty}^\infty f(x)\delta(x)dx=\lim_{n\to \infty}\int_{-\infty}^\infty f(x)\delta_n(x)dx</math>
*נשים לב כי לפי גישה זו <math>\int_{-\infty}^0f(x)\delta(x)dx=0</math> ו<math>\int_0^\infty f(x)\delta(x)dx =f(0)</math>.
*נחשב את התמרת הלפלס של הדלתא של דירק:
*לכל <math>a\geq 0</math> מתקיים <math>\mathcal{L}(\delta(t-a))=\int_0^\infty e^{-st}\delta(t-a)dt=e^{-sa}</math>
*בפרט <math>\mathcal{L}(\delta(t))=1</math>
===תגובת הלם===
*נביט במערכת של מסה המחוברת לקפיץ, המתחילה במנוחה.
*נניח שברגע <math>t=a</math> מישהו נתן 'פליק' למסה.
*הדרך שלנו לבטא כוח נקודתי שכזה היא הדלתא של דירק, המכונה גם 'פונקצית הלם'.
*כלומר הכוח החיצוני על המערכת הוא <math>\delta(t-a)</math>, בנוסף לכוח המופעל על ידי הקפיץ.
*למעשה אנו מעוניינים בפתרון למד"ר <math>y''+ky=\delta(t-a)</math>
*באופן דומה להגדרת האינטגרל, ניתן לחשוב על הפתרון כגבול הפתרונות למערכות המקורבות <math>y''+ky=\delta_n(t-a)</math>.
*על מנת שיהיה פתרון למד"ר עלינו לבחור הפעם סדרה של פונקציות גזירות ב<math>[0,\infty)</math> כמו <math>\delta_n(x)=\begin{cases}ne^{-nx} & x\geq 0 \\ 0 & x<0\end{cases}</math>
*נוכיח כעת את הנוסחא <math>e^{-sa}\mathcal{L}(y(t))=\mathcal{L}(u(t-a)y(t-a))</math> עבור <math>a>0</math>:
**<math>\mathcal{L}(u(t-a)y(t-a))=\int_0^\infty e^{-st}u(t-a)y(t-a)dt = \int_a^\infty e^{-st}y(t-a)dt=</math>
**נבצע את ההצבה <math>x=t-a</math> ונקבל:
**<math>=\int_0^\infty e^{-s(x+a)}y(x)dx =e^{-sa}\int_0^\infty e^{-sx}y(x)dx=e^{-sa}\mathcal{L}(y(t))</math>.
*נפתור את המערכת עם התמרת לפלס:
**<math>\mathcal{L}(y''+ky)=s^2F(s)-sy(0)-y'(0)+kF(s)=e^{-sa}</math>.
**כיוון שהמערכת התחילה במנוחה, <math>y(0)=y'(0)=0</math>.
**לכן <math>F(s)=\frac{e^{-sa}}{s^2+k}</math>.
**ולכן <math>y=u(t-a)\frac{sin(\sqrt{k}(t-a))}{\sqrt{k}}</math>.
**(הרי <math>\mathcal{L}(sin(\sqrt{k}t))=\frac{\sqrt{k}}{s^2+k}</math>).
*אכן, עד רגע <math>t=a</math> המערכת במנוחה <math>y=0</math>.
*לאחר מכן, אנו מקבלים את הפתרון המקיים <math>y(a)=0,y'(a)=1</math>.
*כלומר ה'הלם' תפקד במקרה זה כמו תנאי התחלה על המהירות - זה בדיוק ה'פליק' שהכנו במסה.
*נפתור את המערכת <math>y''+ky=\delta(x-\frac{2\pi}{\sqrt{k}})</math> עם תנאי ההתחלה <math>y(0)=0,y'(0)=-1</math>.
**נפעיל התמרת לפלס <math>s^2F(s)-sy(0)-y'(0)+kF(s)=e^{-\frac{2\pi}{\sqrt{k}}s}</math>
**לכן <math>F(s)=\frac{e^{-\frac{2\pi}{\sqrt{k}}s}-1}{s^2+k}</math>
**לכן <math>y(t)=\frac{1}{\sqrt{k}}\left(u(t-\frac{2\pi}{\sqrt{k}})sin(\sqrt{k}(t-\frac{2\pi}{\sqrt{k}}))-sin(\sqrt{k}t)\right)</math>
**לכן <math>y(t)=\frac{u(t-\frac{2\pi}{\sqrt{k}})-1}{\sqrt{k}}sin(\sqrt{k}t)</math>
**כלומר בזמן <math>t=\frac{2\pi}{\sqrt{k}}</math> ההלם עוצר את התנועה במערכת, והפתרון מתאפס.
*דוגמא - נפתור את המד"ר <math>y'''-y=\delta(t-1)</math> עבור תנאי ההתחלה <math>y(0)=y'(0)=y''(0)=0</math>.
**נבצע התמרת לפלס ונקבל כי <math>s^3F(s)-F(s)=e^{-s}</math>.
**לכן <math>F(s)=\frac{e^{-s}}{s^3-1}=e^{-s}\frac{1}{3}\left(\frac{1}{s-1}-\frac{s+2}{s^2+s+1}\right)</math>
**ראשית נמצא את ההתמרה ההפוכה <math>\mathcal{L}^{-1}\left(\frac{1}{s-1}-\frac{s+2}{s^2+s+1}\right)</math>:
***<math>\mathcal{L}^{-1}\left(\frac{1}{s-1}\right)=e^t</math>
***<math>\frac{s+2}{s^2+s+1}=\frac{s+2}{\left(s+\frac{1}{2}\right)^2+\frac{3}{4}}=\frac{s+\frac{1}{2}+\frac{3}{2}}{\left(s+\frac{1}{2}\right)^2+\frac{3}{4}}</math>
***<math>\mathcal{L}^{-1}\left(\frac{s+\frac{1}{2}}{\left(s+\frac{1}{2}\right)^2+\frac{3}{4}}\right)=e^{-\frac{t}{2}}cos\left(\frac{\sqrt{3}}{2}t\right)</math>
***<math>\mathcal{L}^{-1}\left(\frac{3}{2}\frac{1}{\left(s+\frac{1}{2}\right)^2+\frac{3}{4}}\right)
=\mathcal{L}^{-1}\left(
\sqrt{3}\frac{\frac{\sqrt{3}}{2}}
{\left(s+\frac{1}{2}\right)^2+\frac{3}{4}}
\right)
=\sqrt{3}e^{-\frac{t}{2}}sin\left(\frac{\sqrt{3}}{2}t\right)
</math>
***לכן <math>\mathcal{L}^{-1}\left(\frac{1}{s-1}-\frac{s+2}{s^2+s+1}\right)=e^t-e^{-\frac{t}{2}}cos\left(\frac{\sqrt{3}}{2}t\right)-\sqrt{3}e^{-\frac{t}{2}}sin\left(\frac{\sqrt{3}}{2}t\right)</math>
**ולכן סה"כ הפתרון למד"ר הינו <math>
y=\frac{u(t-1)}{3}\left[
e^{t-1}-e^{-\frac{t-1}{2}}cos\left(\frac{\sqrt{3}}{2}(t-1)\right)-
\sqrt{3}e^{-\frac{t-1}{2}}sin\left(\frac{\sqrt{3}}{2}(t-1)\right)
\right]</math>
==הרצאה 13 - משוואת אוילר==
*משוואת אוילר הומוגנית היא משוואה מהצורה:
**<math>a_nx^ny^{(n)}+a_{n-1}x^{n-1}y^{(n-1)}+...+a_0y=0</math>
*נסמן את פונקצית האקפוננט <math>\exp(t)=e^t</math>
*נפתור את המד"ר ל<math>x>0</math>
*נגדיר <math>u=y\circ \exp</math> כלומר <math>u(t)=y(e^t)</math>.
*נקבל כי
**<math>u'(t)=e^ty'(e^t)</math>
**<math>u''(t)=e^{2t}y''(e^t)+e^ty'(e^t) = e^{2t}y''(e^t)+u'(t)</math>
**<math>u'''(t)=e^{3t}y'''(e^t) + 2e^{2t}y''(e^t)+u''(t) = e^{3t}y'''(e^t)+2(u''(t)-u'(t))+u''(t)</math>
**באופן כללי ניתן להוכיח באינדוקציה כי <math>u^{(m)}(t)=e^{mt}y^{(m)}(e^t)+\sum_{k=1}^{m-1} b_ku^{(k)}(t)</math> עבור קבועים כלשהם.
*נסמן את האופרטור המתאים למד"ר:
**<math>H=a_n x^n D^n +...+a_0 I</math>
**לכן <math>Hy\circ\exp (t)=a_n e^{nt}y^{(n)}(e^t)+...+a_0y(e^t)</math>
**לפי הפיתוח לעיל, זה שווה ל<math>c_nu^{(n)}(t)+...+c_0u(t)=0</math> עבור קבועים כלשהם.
*נסמן את האופרטור המתאים למד"ר זו ב<math>K=c_nD^n+...+c_0I</math>
**סה"כ הוכחנו כי <math>Hy\circ\exp=K(y\circ\exp)</math>
*את הגרעין של <math>K</math> אנחנו יודעים למצוא כיוון שזו מד"ר לינארית הומוגנית עם מקדמים קבועים.
*אם <math>u</math> פתרון למד"ר המתאים ל<math>K</math> אז עבור <math>y=u\circ \ln</math> מתקיים כי <math>K(y\circ\exp)=0</math>
*לכן <math>Hy\circ \exp =0</math> ולכן <math>Hy=0</math> בחיוביים, שהרי זו התמונה של <math>\exp</math>.
*אבל איך נמצא את הפתרונות ל<math>Ku=0</math>? צריך למצוא את הפולינום האופייני.
*עבור <math>y=x^r</math> נקבל כי <math>Hy\circ\exp=K(y\circ\exp)=K(e^{rt})=c_nr^n e^{rt}+...+c_0 e^{rt}</math>
*אם נחלק ב<math>e^{rt}</math> נקבל את הפולינום האופייני של המד"ר <math>Ku=0</math>, זו נקראת '''המשוואה האינדנציאלית''' של משוואת האוילר המקורית.
*במילים פשוטות, על מנת לחשב את המשוואה האינדנציאלית:
**נציב <math>x^r</math> במשוואת האוילר
**נציב <math>x=e^t</math> ונחלק ב<math>e^{rt}</math> (או בעצם נחלק מראש ב<math>x^r</math> שזה שקול)
*השורשים של המשוואה האינדנציאלית נותנים לנו את הפתרונות לגרעין של <math>K</math>, נרכיב אותם על <math>ln(x)</math> ונקבל את הפתרונות למשוואת האוילר.
*סה"כ אם r שורש ממשי מריבוי k של המשוואה האינדנציאלית אזי:
**<math>u(t)=t^me^{rt}</math> פתרון של המד"ר <math>c_nu^{(n)}(t)+...+c_0u(t)=0</math> לכל <math>0\leq m\leq k-1</math>.
**ולכן <math>y(x)=u(ln(x))=ln^m(x)x^r</math> פתרון של משוואת אוילר המקורית, לכל <math>0\leq m\leq k-1</math>.
*אם <math>r=a\pm bi</math> זוג שורשים מרוכבים צמודים מריבוי k כל אחד אזי:
**<math>u(t)=t^me^{at}cos(bt),t^me^{at}sin(bt)</math> פתרונות של המד"ר <math>c_nu^{(n)}(t)+...+c_0u(t)=0</math>, לכל <math>0\leq m\leq k-1</math>.
**לכן <math>y(x)=ln^m(x)x^acos(bln(x)),ln^m(x)x^asin(bln(x))</math> פתרונות של משוואת אוילר המקורית, לכל <math>0\leq m\leq k-1</math>.
*דוגמא:
**<math>x^3y'''-x^2y''+2xy'-2y=0</math>
**נציב <math>y=x^r</math> ונקבל את המשוואה האינדנציאלית <math>r(r-1)(r-2)-r(r-1)+2r-2=0</math>.
**לכן <math>r(r-1)(r-2)-(r-2)(r-1)=0</math>.
**כלומר <math>(r-2)(r-1)(r-1)=0</math>.
**לכן הפתרון הכללי הינו <math>y=c_1x^2+c_2x+x_3xln(x)</math>
*דוגמא:
**<math>xy''+y'+\frac{y}{x}=0</math>
**נעביר לצורה של משוואת אוילר <math>x^2y''+xy'+y=0</math>.
**המשוואה האינדנציאלית היא <math>r(r-1)+r+1=0</math>.
**כלומר <math>r^2+1=0</math>.
**לכן הפתרון הכללי הינו <math>y=c_1sin(ln(x))+c_2cos(ln(x))</math>
*דוגמא:
**מצאו פתרון כלשהו למד"ר <math>x^2y''-2xy'+2y=x^3e^x</math>
**ראשית נמצא את הפתרונות למד"ר ההומוגנית, שהיא משוואת אוילר.
**לאחר מכן נמצא פתרון פרטי באמצעות וריאצית המקדמים.
