אנליזה נומרית/פתרון משוואות דיפרנציאליות: הבדלים בין גרסאות בדף

תוכן שנמחק תוכן שנוסף
Mintz l (שיחה | תרומות)
מאין תקציר עריכה
Mintz l (שיחה | תרומות)
שורה 67:
{{גזור ושמור|שיטה רב-צעדית:<br /><math>\ y_{i+1}= y_i+ h\left[ 1+{\nabla\over 2}+ {5 \over 12}\nabla^2+ {3\over 8}\nabla^3+... \right]f(x_i,y_i)</math>}}
 
*שימו לב כי אם ניקח איבר בודד, נקבל את שיטת אוילר.
===דוגמאות===
 
* אם ניקח איבר בודד, נקבל את שיטת אוילר.
* '''===שיטת Adams-Bashforth''': ===
נפתח את 4 האיברים הראשונים שציינו לעיל:
:<math>\ y_{i+1}= y_i+{h\over 24} \left[ 55f(x_i,y_i)-59f(x_{i-1},y_{i-1})+37f(x_{i-2},y_{i-2})-9f(x_{i-3},y_{i-3}) \right]+ O(h^5)\ ,\quad i\ge 3</math>
:*שימו לב כי לשיטה זו יש לספק 3 תנאי התחלה, אשר בדרך כלל מוצאים אותם בשיטת רונגה-קוטה או בשיטת אוילר. כמו כן, על אותה שיטה להשתמש באותו מרווח h.
:*נשווה שיטה זו לשיטת רונגה-קוטה מסדר 4, מבחינת מספר פעולות החישוב: <math>\ \tfrac{R-K}{A-B}= \tfrac{4N}{4\cdot 3+(N-3)}=\tfrac{4N}{N+9}\approx 4</math>, כלומר שיטת רונגה-קוטה איטית פי 4.
:*על מנת להשתמש בשיטה זו עבור '''מערכת מד"ר''', נחליף את הכתיב הסקלרי בכתיב הוקטורי:
::<math>\ \vec y_{i+1}= \vec y_i+{h\over 24} \left[ 55\vec f(x_i,\vec y_i)-59\vec f(x_{i-1},\vec y_{i-1})+37\vec f(x_{i-2},\vec y_{i-2})-9\vec f(x_{i-3},\vec y_{i-3}) \right]+ O(h^5)\ ,\quad i\ge 3</math>
 
* שיטות נוספות:
===שיטת הטרפז===
:* <math>\ y_{i+1}= y_{i-1}+2hf(x_i,y_i)+ O(h^3)</math>
(להשלים)
:* <math>\ y_{i+1}= y_{i-3}+ \tfrac{4h}{3}[2f(x_i,y_i)-f(x_{i-1},y_{i-1})+ 2f(x_{i-2},y_{i-2})] +O(h^5)</math>
 
* ===שיטות נוספות:===
:* <math>\ y_{i+1}= y_{i-1}+2hf(x_i,y_i)+ O(h^3)</math>
:* <math>\ y_{i+1}= y_{i-3}+ \tfrac{4h}{3}[2f(x_i,y_i)-f(x_{i-1},y_{i-1})+ 2f(x_{i-2},y_{i-2})] +O(h^5)</math>
 
===קישורים חיצוניים===