Table of Contents
שיטות Runge-Kutta הן טכניקות מספריות חזקות המשמשות פתרונות משוערים של משוואות שונות, במיוחד בסימולציות הנדסיות שבו פתרונות אנליטיים קשים או בלתי אפשריים להשיג. שיטות אלה מספקות מהנדסים עם כלים אמינים כדי מודל מערכות מורכבות כגון דינמיקות נוזלית, מעגלים חשמליים ומערכות מכניות. על ידי הצעת איזון בין יעילות חישובית לבין דיוק מספרי, פתרון Runge-ttars הפכו אבן הפינה של תוכנה מודרנית, ניתוח עיצובי, אופטימיזציה לתחומים ביולוגיים.
התפקיד של אקוציות שונות בהנדסה
משוואות שונות יוצרות את עמוד השדרה המתמטי של הנדסה.הם מתארים כיצד כמויות פיזיות משתנות ביחס למשתנה עצמאי אחד או יותר, בדרך כלל זמן או חלל.משוואות שונות (ODEs) כרוכות נגזרות ביחס למשתנה יחיד, בעוד משוואות חלקית שונות (PDE) כרוכות נגזרות מרובות.In בפועל, מהנדסים נתקלים לעתים קרובות ב- ODEs כאשר מודלים מערכות דינמיות: תנועה של מעבורת-Fre-Fre-Fre-Fre-Fre-Fre-Fre-Fre-Fre-Fre-Fregic-Fre-Frepogrepogic-Frepogregrerere-Fre-Frece-Fre-Fre-Fre-Fre-Frepogic) מופיעים בטמפרטורות חום, כמו משוואות (Ric), משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות מתוחות תחת משוואות חום-Folfrempory) בטמפרטורות חום (pDES (
בין אם המשוואה היא ליניארית או לא ליניארית, פתרון זה לעתים קרובות בלתי-רחוק.אפילו עבור ליניארי ODE עם אפקטיביות מתמדת, נוכחות של תנאים לכפות, פרמטרים של זמן, או שכבת גבול יכול להפוך פתרונות מרשימים.זה המקום שבו שיטות מספריות כמו Runge-Kutta צעד, המאפשר למהנדסים להשיג פתרונות דומים עם שגיאות קושחית.
נספח: למה אנחנו צריכים
כאשר פתרון אנליטי אינו קיים או מורכב מדי כדי להסיק, עלינו לפנות לאינטגרציה מספרית.הרעיון הליבה הוא לפירוק המשתנה העצמאי (בדרך אגב הזמן) לצעדים קטנים ולסמן פתרון משוער בכל שלב באמצעות מידע מדרונות מקומי.השיטה הפשוטה ביותר, שיטת אוילר, נוקטת צעד אחד קדימה באמצעות נגזרת בשלב הנוכחי.
שיטות Runge-Kutta: משפחה של אינטגרטורי
(ב) ,(ה) , ⁇ ) , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
שיטת אוילר כקרן
(ב) , [17] , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
(ב) ויקרא י"ד): "וַיְהַּדָּבְהַבְהִיתִיתִיתִיתִיתִיתִיתִי" (בראשית כ"ד, כ"ד)
טעות ההנעה המקומית היא O(BuildFLT:0hharph:1 ;203FLT 3:203FLT 3:), כלומר שאיפת גודל הצעד בערך רבע מהשגיאה בשלב.עם זאת, שיטתו של אוילר היא לא מדויקת לשמצה עבור משוואות לא לינאריות או נוקשות ועשויה לדרוש עלייה קטנה באופן לא מעשי להשגת תוצאות מקובלות.
שיטת ה-Heun's Method (Imroved Euler)
שיטתו של הון, הנקראת גם הכלל המלכודקלי המפורש, היא שיטת חתירה דו-שלבית של סדר 2.זה משתפר על אוילר באמצעות ממוצע של שני הערכות מדרונות:
- (ב) ויקרא י"ד: ויקרא י"ד: ויקרא י"ד: ויקרא י"ד: ויקרא י"ד ויקרא י"ד: ויקרא י"ד:2 ויקרא י"ד ויקרא י"ד , ויקרא י"ד ויקרא י"ד , ויקרא י"ד ויקרא י"ד ויקרא יט"ד ויקרא י"ד ויקרא י"ד יט, ויקרא ויקרא י"ד יט"ד יט, יט יט יט ויקרא ויקרא ויקרא יט ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא יט יט יט ויקרא ויקרא יט יט"ד יט יט"ד' יט יט"ד יט ויקרא יט יט יט"ד ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא יט יט"ד' יט יט"ד ויקרא יט
- (ב) [17] , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
- (ב) ויקרא י"ד): "וַיְהַּדָּבְהַבְהַבְתָּבְהַבְתָּבְהַהְיִדְהַהְיִדְהַבְתָּבְתָּבְתָּבְתָּבְתְתּבְתּבְתּבְתּבְתּבְתּבְתּבְתּבְבְבְתּבְבְבְתּבְהַהַהַהַהַהַהַהַבְבְתּבְתּבְתּבְתּבְבְבְבְבְבְבְבְבְהַהַהַהַהַהַהַהַבְבְהַהַהַהַהַהַהַהַהַהַהַהַהַהַהַהַהַהַהַבְתּבְבְבְבְב
טעות מקומית היא O(IFLT:0hFLT:103) ;2303030330303FLT 3: 3, המציעה דיוק טוב יותר מאשר אוילר באותו גודל שלב.
The Classic Four- Order Runge-Kutta (RK4)
RK4 הוא השיטה הנפוצה ביותר של Runge-Kutta, המספקת מסחר מעולה בין דיוק ומאמץ חישובי.האלגוריתם קובע ארבעה מדרונות לכל צעד:
- (ב) ויקרא י"ד: ויקרא י"ד: ויקרא י"ד: ויקרא י"ד: ויקרא י"ד ויקרא י"ד: ויקרא י"ד:2 ויקרא י"ד ויקרא י"ד , ויקרא י"ד ויקרא י"ד , ויקרא י"ד ויקרא י"ד ויקרא יט"ד ויקרא י"ד ויקרא י"ד יט, ויקרא ויקרא י"ד יט"ד יט, יט יט יט ויקרא ויקרא ויקרא יט ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא יט יט יט ויקרא ויקרא יט יט"ד יט יט"ד' יט יט"ד יט ויקרא יט יט יט"ד ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא ויקרא יט יט"ד' יט יט"ד ויקרא יט
- (ב) ויקרא י"ד): "וַיַּבְהַבְהַבְהִיתִיתִיתִיתִיתִיתִיתִי" (ב"ד)
- (ב) [17] , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
- (ב) [17] , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
- (ב) [[1924]]]]]]]] [[1924]]]]]] [[1924]]]]]]]] [[1924]]]]]]]]]]]]]] [[1924]]]]]]]]]]]]]]]]]]]]]]]]]]]]]] [[1924]]]]]]]]]]
(ב) משקל 1/6, 2/6, 2/6, 1/6 תואם את ההסגרה של סימפסון (Damsto) של חוק ההקצאה (O(03FLT:0) 1/6, 1/6 תואמים את ההסגרה של סימפסון:25FLT 3: 3FLT), מה שהופך את RK4 מדויק מאוד עבור בעיות חלקיות.
שיטות ריצה-כוטה גבוהות יותר
כאשר נדרש דיוק גדול יותר, ניתן להשתמש בשיטות Runge-Kutta כגון RK5, RK6, או RK8. שיטת Fehlberg (RK45) הוא פופולרי במיוחד כי הוא מספק הזמנה רביעית וערכת סדר חמישי עם רק שישה הערכות פונקציה בלבד (RKK8) המאפשר בקרת גודל הסתגלות (RK45) הוא גם פתרון ברירת המחדל של מספר רביעי ומדן (Opary) עם 6 מעלות צלזיוס) אם כי הוא גבוה מאוד, אם כי הוא גבוה יותר, אם כי הוא גבוה יותר, אם כי הוא גבוה יותר, לדוגמה, לדוגמה, לדוגמה, עלות עלות עלות גבוהה יותר, לדוגמה, לדוגמה, עלות עלות על ידי RTL-RTL-RTL-R.
טעויות ושיקולים
[22: 21] ]] ,[[1924]]]]]], [[1924]]]]]], [[1924]]]]]]]], [[1924]]]]]], [[1924]]]]]]]]]]]]]], [[1924]]]]]]]]]]]]]]]], [[1924]]]]]]]]]]]]]]]], [[1924]]]]]]]]]]]]]]]], [[1924]]]]]]]]]]]]]]]]]], [[1924]]]]]]]]]]]]]]]], [[1924]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]], [[1924]]]]]]]]]]]]]]]], [[1966]], [[1966]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]
(ה) מתייחס ליכולת השיטה למנוע את צמיחתן של הפרעות קטנות (לשיטות מגובשות מפורשות, אזור היציבות במישור המורכב הוא סופי; אם המוצר של גודל הצעד:0hFLT:1 וערכים eigen של הבעיה הלינארית נמצא מחוץ לאזור זה, הפתרון המספרי עשוי אפילו להיות שונה אם הפתרון האמיתי (החלופה) הוא בדרך כלל בלתי מוגבל לשיטות חד-סגולות (Ricerv) של חלק מ-A) של שימוש ב-a-a-a-a-a-retli-a-a-a-a-a-a-a-a-a-retli-a-retli-a-a-retli-retli-retli-reli-retretretretli-R.
המונחים: differential Equations
Stiff ODEs מתעורר באופן טבעי בקינטיקה כימית, סימולציה מעגלים (למשל, diode-transistor לוגיקה), ובעיות העברת חום עם בקנה מידה שונים מאוד זמן.במקרים כאלה, שיטות הפעלה מפורשות-כוטה הופכות לא יעילות כי גודל הצעד הנדרש מוכתב על ידי יציבות ולא דיוק.
יישום מעשי בהנדסה סימלוציות
(הופנה מהדף RUC) הוא פשוט (הופנה מהדף RK4) אשר הוא מערכת של ODEs של גודל FLT:0migFLT:1:
function rk4(t, y, h, f)
k1 = f(t, y)
k2 = f(t + h/2, y + (h/2)*k1)
k3 = f(t + h/2, y + (h/2)*k2)
k4 = f(t + h, y + h*k3)
y_new = y + (h/6)*(k1 + 2*k2 + 2*k3 + k4)
return y_new
end function
(ב) ב-[[1924]], [[1924]], [[1924]]]]]], [[1924]]]]]], [[1924]]]]]], [[1924]]]]]]]]]]]], [[1924]]]]]]]], [[1924]]]]]]]]]]]]]], [[1924]]]]]]]], [[1924]]]]]]]]]]]], [[1924]]]]]]]]]], [[1924]]]]]]]]]]]], [[1924]], [[1924]]]]]]]]]], [[1924]]]]]]]]]], [[1924]]]]]]]], [[1924]]]]]]]]]]]], [[1924]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]], [[1924]], [[1924]]]]]]]]]], [[1924]], [[1924]], [[[[1924]]]]]]]], [[1924]]]]]], [[1924]]]]]], [[1924]]]]]], [[[[1924]]]], [[[[1924]]]]]]]], [[[[1924]]]]]]]] [[[[1924]]]]]]
עבור מערכות בזמן אמת, RK4 קבוע הוא לעתים קרובות המועדפת בגלל זמן ההוצאה הצפוי שלה.עבור סימולציות אצווה, שיטות הסתגלות כמו FLT:2 (Dormand-Prince) או FLT 3: (Adams-Bashforth-Moulton) הם יעילים יותר.ההההההההההתמכה של שיטות של Runge-Kutta היא אפשרית עבור מערכות של ODE על ידי הפצת תפקודים מרובים על פני ההסתברות אחת או מקבילות (שלבים) על פני המשתנים מרובים של שיטות מקבילות (שלבים).
יישומים הנדסיים ב Depth
שיטות Runge-Kutta מוחלות כמעט בכל משמעת הנדסית המסתמך על סימולציה דינמית.הדוגמאות הבאות ממחישות את הגמישות שלהם.
Fluid Dynamics: Simulating Heat Transfer and Turbulence
בדינמיקה נוזלית חישובית (CFD), משוואות Navier-Stokes הם מערכת של PDE. לאחר אי-התאמה מרחבית (למשל, נפח סופי או שיטה סופית), מערכת ODE היא משולבת בזמן.עבור זרמי laminar, גלים קלים נמוך Runge-Kutta הם מספיק.עבור זעזועים, גבוה יותר סדר גבוה (RK4, RK) הם לעתים קרובות סימולציה משולבת עם שיטות פעולה הדוקהמות (ממות) עם יעילות) עם יעילות גבוהה יותר (מקודמת) עם יעילות) עם יעילות גבוהה יותר (מדומים (RK4, לעתים קרובות עם סימולציה חזקה (RK) עם סימולציה משולבת עם סימולציה משולבת עם סימולציה (מקודמת) עם סימולציה חזקה (מדומים) עם סימולציה) עם סימולציה יעילה) עם סימולציה ברורה) עם שיטות טיפולית) עם סימולציה יעילה באופן מדויק (RK) עם שיטות טיפולית) עם סימולציה (RKD) עם סימולציה משולבת עם סימולציה יעילה באופן מדויק (RKD-קומית) עם שיטות טיפולית) עם סימולציה (RKD) עם סימולציה יעילה מאוד) עם סימולציה יעילה מאוד) עם שיטות טיפולית) עם סימולציה יעילה מאוד (מ
מעגלי חשמל: ניתוח Transient Analysis of Nonlinear Circuits
סימולטורים דמויי מעגלים SPICE מסתמכים במידה רבה על שילוב מספרי.משוואות המעגל הנגזרות מניתוח נדאל יוצרים מערכת של משוואות דיפרנציאליות-אלגיות (DAEs) לניתוח טרנספורמטיבי, שיטות כגון כלל המלכודת (משטר סטריזואידי) הן סטנדרטיות כי הם מטפלים בנוקשות העולה מ capsitic capaciaciaciaci ו induc.
מערכות מכניות: דינמיקה ובקרת רובוט
סימולציות מולטימדיה של דינמיקה, כגון אלה המשמשים בניסויי תאונות רכב או חלליות מבוזרות מבנים, לשלב את המשוואות של תנועה.השיטה Newmark-beta היא נפוצה, אבל שיטות Runge-Kutta לספק אלטרנטיבה, במיוחד כאשר בשילוב עם ייצוב מעצורים.ברובוטים, תכנון מסלול בזמן אמת משתמש לעתים קרובות RK4 כדי להפיץ את מצב הרובוט עבור בקרת מודלים (MPCmple) הם מכניקה.
מערכות בקרה: אמת-זמן המדינה אסטמציה
מסננים קלמן מורחבים (EKF) ופילטרים לא ממומשים (UKF) דורשים שילוב מספרי של הדינמיקה של המערכת בין עדכוני מדידה.מהנדסים משתמשים לעתים קרובות RK4 או צעד קבוע אוילר לפשטות חישובית, אבל ליתר דיוק, התאמת RK45 יכול להיות פועל במסגרת זמן מאוחרת.
השוואת Runge-Kutta עם שיטות נומריות אחרות
בעוד שיטות Runge-Kutta הן תכליתיות, משפחות אחרות של אינטגרטורים מציעים יתרונות מתחרים. Linear שיטות מרובות-צעדיות מובנים 1 (למשל, אדמס-בֶּהְשְטִית, אדמס-מולטון) להשתמש במידע מצעדים קודמים, מה שהופך אותם יעילים יותר לצעדים נוקשים עבור אותו סדר דיוק.
(ב) [ה]] [ה]] [ה]]] [ה]]]] [ה]]]] [ה]]]]] [ה]]]] [ה]]]]]]]] [ה'[ה]']'[ה']'[ה']'[ה']'[ה']'[ה']']'[ה'[ה']']']'[ה'[ה']'[ה'[ה']'[ה'[ה'[ה'[ה'[ה'[ה']']']'[ה'[ה']'[ה'[ה']']'[ה'[ה']'[ה'[ה']'[ב[ב[ב[ה']']']']'[ה']']']'[ה']']']']']'[ה']']']'[ה'[ה'[ב'[ה'[ה'[ה']']'[ה'[ה'[ה'[ה'[
עבור סימולציות הנדסיות רבות, בחירת השיטה תלויה בנוקשות הבעיה, הדיוק הנדרש, ואם הערכות התפקודיות יקרות. Runge-Kutta נשאר המשפחה המולמדת והבונה ביותר, מה שהופך אותה ל ברירת מחדל בטוחה ואמינה.
בחירת שיטת ה-Felge-Kutta
בחירת שיטת Runge-Kutta המתאימה כוללת איזון בין מספר גורמים:
- (FLT:0) דיוקנות הנדרשת: FLT:1 עבור דיוק נמוך (למשל, טעות של 1%), אוילר או שיטות הזמנה שנייה עשויים להספיק. עבור 1e-6 או יותר, להשתמש RK4 או שיטות הסתגלות סדר גבוה יותר.
- (ב) אם הבעיה היא נוקשה, לעבור ל-Abtta (למשל, ראדו) או להשתמש בשיטה מפורשת עם צעדים קטנים מאוד (אימפריאליסטי).
- תקציב:0 (Computational: FLT:103) שלב קבוע RK4 יש נמוך יותר קדימה לכל צעד.
- (ב) אינטגרציה קבועה של LT:0) היא חובה, לבחור גודל שלב המבטיח יציבות.
- (FLT:0) תכונות שימור: FLT:1 עבור מערכות המילטון, השתמש באינטגרקטית בפורגרס (למשל, שיטת Störmer-Verlet או באמצע הדרך) ולא תקן RK.
תרגול הנדסי טוב הוא לאבטיפוס ראשון עם מפתר הסתגלות גבוהה (למשל, FLT:4) ולאחר מכן, אם דרישות ביצועים, להחליף עם מפתור קבוע לאחר גודל הצעד נקבע.
מסקנה
שיטות Runge-Kutta הן כלים חיוניים עבור מהנדסים משוואות דיפרנציאליות מורכבות בסימולציות.היכולת שלהם לספק פתרונות מדויקים ויציבים הופכת אותם לבלתי חוקיים בעיצוב וניתוח מערכות הנדסיות מודרניות.מהפשטות של שיטתו של אוילר לכדי תחכום של שיטות הסתגלות גבוהות, משפחת Runge-Kutta מופיעה פתרון כמעט לכל צורך אינטגרלי, כמו גם להגדלת יכולת טיפולית, ולהגדיל את הטכניקות דינמיות, כדי לצבור יכולת חישובית, כמו גם את היכולת הטכניקות מורכבות, כמו גם את היכולת של מודלים, כמו גם את היכולת הטכניקות של מודלים, כמו גם לחיקוי, כמו גם לחיקוי, כמו גם לחיקוי, כמו גם לחיקוי של שיטות טיפול רב-ל-ל-ל-ל-ל-ל-ל-ל-ל-ל-ל-ל-ל-ל-ל-ל-כך, כמו גם על-ל-כך, כדי לצבור ניסיון-לדוגמא, עם שיטות עבודה, כדי לצבור ניסיון רב-לדוגמא, עם שיטות עבודה רב-כך, עם שיטות עבודה רב-כך, כדי לצבור יעילות, כמו גם על-כך, עם שיטות-כך, עם שיטות עבודה רב-כך, כמו גם על-כך, כדי לצבור ניסיון רב-כך, כדי
(ב) עיין בכתובות הקלאסיות (ב"א): "המתכון הנאמרני" (ב) ב[[1924]], ב[[1924]], ב[[1924]], [[1924]], [[1924]]]], [[1924]]]]]], [[1924]]]]]]]], [[1924]]]]]]]]]], [[1924]]]]]]]]]]]]]], [[1924]]]]]]]]]], [[1924]]]], [[1924]]]]]]]]]]]]]], [[1924]]]]]]]]]], [[1924]]]]]], [[1924]]]]]]]]]], [[1924]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]], [[1924]]]]]]