שיטת רונגה-קוטה

באנליזה נומרית, רוּנְגֶה-קוּטַה היא משפחת שיטות איטרטיביות, ישירות ועקיפות (Implicit and explicit), הכוללות את השיטה הפשוטה והידועה, שיטת אוילר, המבוססת על דיסקרטיזציה טמפורלית להערכת פתרונות מקורבים של משוואות דיפרנציאליות רגילות.תבנית:הערה שיטות אלו פותחו בסביבות שנת 1900 על ידי המתמטיקאים הגרמנים קרל דייוויד רונגה ומרטין וילהלם קוטה. לשיטות שימושים רבים בפיזיקה, בכימיה ובביולוגיה.

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

תיאור

שיטות רונגה-קוטה ממוספרות על פי הסדר הגבוה ביותר של פונקציונלי העזר המשמשים במסגרתן. כך למשל, שיטת רונגה-קוטה מסדר שני עושה שימוש בשני פונקציונלים על מנת לקבל כלל נסיגה המגדיר את וקטור הפתרונות. בנוסף, מתוך מספור השיטה ניתן לקבל מידע על הסדר של שגיאת הקיטוע המקומית.

השיטה המוכרת ביותר במשפחה זו, היא שיטת רונגה-קוטה מסדר רביעי, המוכרת גם בתור "שיטת רונגה-קוטה הקלאסית" או בתור הקיצור RK4. להלן המרשם לפתרון משוואה דיפרנציאלית רגילה מסדר ראשון.

תהי בעיית תנאי שפה (מד"ר מסדר ראשון) מהצורה הבאה:

$\frac{dy}{dt}=f(t,y),y(t_0)=y_0$

כאשר $y$ היא פונקציה לא-ידועה של הזמן אותה נרצה להעריך במסגרת הקירוב. הפסוק הנ"ל אומר כי הנגזרת של $y$ היא פונקציה של הזמן ושל $y$ עצמה. עוד נתון תנאי ההתחלה על $y$, אשר יהווה מרכיב מרכזי בפתרון בשיטת רונגה-קוטה להלן.

כעת נבחר צעד $h$ אי-שלילי כלשהו ונגדיר ארבעה פונקציונלים $k_i$:

$\begin{array}{cc}{k}_1& =f(t_n,y_n),\\ k_2& =f(t_n+\frac{h}{2},y_n+h\frac{k_1}{2}),\\ k_3& =f(t_n+\frac{h}{2},y_n+h\frac{k_2}{2}),\\ k_4& =f(t_n+h,y_n+h{k}_3).\end{array}$

כאשר ניתן לתת פרשנות גאומטרית לכל אחד מהם:

• $k_1$ הוא השיפוע בתחילת המקטע (התחום בין הנקודה הידועה $y_n$ וע"י הנקודה הבאה אחריה $y_{n+1}$).
• $k_2$ הוא השיפוע במחצית המקטע, כמחושב בעזרת $y$ ו-$k_1$.
• $k_3$ הוא שוב, השיפוע במחצית המקטע, אך הפעם כמחושב בעזרת $y$ ו-$k_2$.
• $k_4$ הוא השיפוע בסוף המקטע.

על בסיס ארבעת הפונקציונלים הללו, נבנה את כלל הנסיגה הבא:

$\begin{array}{cc}{y}_{n+1}& =y_n+\frac{1}{6}h(k_1+2k_2+2k_3+k_4),\\ t_{n+1}& =t_n+h\\ \end{array}$

לעיתים, נהוג הסימון $T_4$ לממוצע הפונקציונלים המשוקלל המופיע בכלל הנסיגה לעיל. ניתן לראות כי משקל רב יותר ניתן לשיפוע במחצית המקטע.

אם $f$ לא תלויה ב-$y$ הרי שהאינטגרציה טריוויאלית ושיטת רונגה-קוטה מסדר רביעי שקולה להפעלתו של תבנית:קישור שפה.

שגיאת הקיטוע המקומית בשיטה הנ"ל היא מסדר גודל של $h^5$, בעוד השגיאה הנצברת הכוללת היא מסדר גודל של $h^4$.

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

• תבנית:MathWorld