שינויים

קפיצה אל: ניווט, חיפוש
/* מבחנים לדוגמא */
*[[מדיה:22ODETestA.pdf|מבחן מועד א' תשפ"ב]]
*[[מדיה:22ODETestB.pdf|מבחן מועד ב' תשפ"ב]]
*[[מדיה:23ODETestA.pdf|מבחן מועד א' תשפ"ג]]
=הרצאות=
 *טענה - עבור פונקציות <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> הוא פתרון פרטי של המד"ר.**הוכחה:
מתקיים כי <math>y_p=c_1(x)y_1+...+c_n(x)y_n</math> הוא פתרון פרטי של המד"ר.
 
 
*הוכחה:
**<math>y_p'=c_1'y_1+\cdots+c_n'y_n+c_1y_1'+\cdots+c_ny_n'=c_1y_1'+\cdots+c_ny_n'</math>. (לפי המשוואה הראשונה.)
**באופן דומה <math>y_p''=c_1y_1''+\cdots+c_ny_n''</math>. (לפי המשוואה השנייה.)
***<math>y_p^{(n)}+a_{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>y_1,...,y_n</math> פתרונות למד"ר ההומוגנית הביטויים בסוגריים מתאפסים וסה"כ קיבלנו כי אכן <math>y_p'''+a_2(x)y_p''+a_1(x)y_p'+a_0(x)y_p=f(x)</math>.
 
 
*נכתוב '''שוב''' את ההוכחה, בעזרת סימן הסכימה (עשוי להיות נוח יותר או פחות):
**ראשית, ניתן להוכיח באינדוקציה כי לכל <math>0\leq m\leq n-1</math> מתקיים כי
**<math>D^m y_p = D^m \sum_{k=1}^n c_k(x)y_k = \sum_{k=1}^n c_k(x)D^m y_k</math>
**כעת בעזרת המשוואה האחרונה נקבל כי
**<math>D^n y_p = D D^{n-1}y_p = D\sum_{k=1}^nc_k(x)D^{n-1}y_k=\sum_{k=1}^n c'_k(x)D^{n-1}y_k + \sum_{k=1}^nc_k(x)D^ny_k=f(x)+\sum_{k=1}^nc_k(x)D^ny_k</math>
**נציב במד"ר ונקבל
**<math>Ty_p=D^ny_p +\sum_{t=0}^{n-1}a_t(x)D^ty_p=f(x)+\sum_{k=1}^nc_k(x)D^ny_k + \sum_{t=0}^{n-1}a_t(x)\left(\sum_{k=1}^n c_k(x)D^t y_k\right)=</math>
**<math>=f(x)+\sum_{k=1}^n c_k(x)\left(D^ny_k + \sum_{t=0}^{n-1}a_t(x)D^t y_k\right) = f(x)+0</math>
====שתי מסות על קפיץ - מערכת מד"ר מסדר שני====
*נביט בשתי מסות המחוברות לשני צידי קפיץ.
*נניח כי <math>y_1,y_2<y_1</math> מודדות את מיקום המסות ביחס לנקודת האפס שלהן, וצד ימין הוא הכיוון החיובי בשתיהן.
*נניח כי כאשר כל אחת מהמסות במקום אפס, אזי הקפיץ במנוחה.
*נניח כי המסות זהות בגודלן, ושוות אחד.
*לכן מתקבלת מערכת המד"ר <math>\begin{cases}y_1''=-k(y_1-y_2) \\ y_2''=-k(y_2-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>.
==הרצאה 10 התמרת לפלס==
*התמרת לפלס היא העתקה לינארית בין מרחבי פונקציות.
*עבור הפונקציה <math>y(t)</math> המוגדרת בקטע <math>[0,\infty)</math> נגדיר את התמרת הלפלס <math>F(s)=\mathcal{L}(y)=\int_0^\infty e^{-st}fy(t)dt</math>.
*שימו לב שנהוג לסמן את הפונקציה לפני ההתמרה עם המשתנים x או t, ולאחר ההתמרה נהוג להתמש במשתנה s.
*אם מתקיים כי <math>|y(t)|\leq Me^{at}</math> אזי ההתמרה מתכנסת לכל <math>s>a</math>.
**הביטוי האחרון מתכנס לכל <math>s>a</math>.
*נניח כי כל הפונקציות שאנו עוסקים בהן חסומות על ידי אקספוננט באופן דומה.
 
 
*דוגמא - נחשב את ההתמרה של הפונקציה <math>e^{at}</math>.
**<math>F(s)=\mathcal{L}(e^{at})=\int_0^\infty e^{-st}e^{at}dt = \int_0^\infty e^{(a-s)t}dt = \left[\frac{e^{(a-s)t}}{a-s}\right]_0^\infty</math>
**לכל <math>s\geq a</math> האינטגרל הלא אמיתי מתכנס ונקבל כי <math>F(s)=\frac{1}{s-a}</math>
**במילים פשוטות התמרת לפלס של הפונקציה <math>e^{at}</math> הינה הפונקציה <math>\frac{1}{s-a}</math>.
 
 
*דוגמא - נחשב את ההתמרה של הפונקציה <math>sin(at)</math>.
**<math>F(s)=\mathcal{L}(sin(at)) = \int_0^\infty e^{-st}sin(at)dt</math>
**נבצע אינטגרציה בחלקים
**<math>\int_0^\infty e^{-st}sin(at)dt = \left[\frac{e^{-st}}{-s}sin(at)\right]_0^\infty + \frac{a}{s}\int_0^\infty e^{-st}cos(at)dt = \frac{a}{s}\int_0^\infty e^{-st}cos(at)dt </math>
**נבצע אינטגרציה בחלקים על האינטגרל החדש
**<math>\mathcal{L}(cos(at))=\int_0^\infty e^{-st}cos(at)dt = \left[\frac{e^{-st}}{-s}cos(at)\right]_0^\infty - \frac{a}{s}\int_0^\infty e^{-st}sin(at)dt = \frac{1}{s} - \frac{a}{s}F(s)</math>
**ביחד נקבל כי
**<math>F(s) = \frac{a}{s} \left[\frac{1}{s} - \frac{a}{s}F(s)\right]</math>
**נבודד את <math>F(s)</math> ונקבל כי
**<math>\mathcal{L}(sin(at)) = F(s) = \frac{a}{s^2+a^2}</math>
 
 
 
*דוגמא - נחשב את ההתמרה של הפונקציה <math>cos(at)</math>.
**במהלך הדוגמא הקודמת קיבלו את השיוויון
**<math>\mathcal{L}(sin(at)) = \frac{a}{s} \mathcal{L}(cos(at))</math>.
**ולכן <math>\mathcal{L}(cos(at)) = \frac{s}{a}\mathcal{L}(sin(at)) = \frac{s}{a}\cdot\frac{a}{s^2+a^2}=\frac{s}{s^2+a^2}</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>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 - המשך התמרת לפלס==