Das Runge-Kutta-Verfahren

Dem Polygonzugverfahren liegt die Voraussetzung zugrunde, dass die Änderung der Funktion zwischen x = x i und x + h = x i + 1 allein durch den Anstieg G ( x i ;   y ¯ i ) an der Stelle x i bestimmt ist. Da die exakte Lösung f diese Voraussetzung im Allgemeinen nicht erfüllt – ihr Anstieg kann sich zwischen x i  und x i + 1 ändern – entstehen Abweichungen zwischen y i + 1 = f ( x i + 1 ) und y ¯ i + 1 = f ¯ ( x i + 1 ) . Diese können durch ein kleineres h und damit größeren Rechenaufwand reduziert werden.

Eine andere Möglichkeit zur Reduzierung der Abweichungen besteht darin, den Anstieg der Funktion f zwischen x i  und x i + 1 besser als nur durch G ( x i ;  y i ) (Anstieg am Anfang des Intervalls) abzuschätzen. Hierfür lieferten die deutschen Mathematiker CARL RUNGE (1856 bis1927) und MARTIN WILHELM KUTTA (1867 bis 1944) eine weit verbreitete Variante.

Der Term G ( x i ; y ¯ i ) für den Anstieg m i im Punkt P ( x i ; y ¯ i ) des Richtungsfeldes wird dabei ersetzt durch einen gemittelten Anstiegswert über vier Punkte P 1  bis P 4 .

Bild

Die Berechnung der numerischen Lösung folgt dann wieder der Grundgleichung y ¯ i + l = y ¯ i + m i ( r u -ka ) h (Berechnung von m i ( r u -ka ) siehe Tabelle).

Die Berechnung einer numerischen Lösung nach RUNGE und KUTTA sollte aus Effektivitätsgründen am PC durchgeführt werden. Die Programmierung des Algorithmus kann in einer Programmiersprache oder mittels Tabellenkalkulation erfolgen.

Anwendung des Polygonzugverfahrens und des RUNGE-KUTTA-Verfahrens

Beispiel: Das Anfangswertproblem f ( x ) + f ( x ) = 2  mit  f ( 0 ) = 0 soll nach dem Polygonzugverfahren und dem RUNGE-KUTTA-Verfahren mit h = 0,8 im Intervall [ 0 ; 4 ] näherungsweise gelöst und mit der exakten Lösung y = f ( x ) = 2 e x + 2 verglichen werden.

Die Differenzialgleichung f ( x ) + f ( x ) = 2 wird umgewandelt in f ( x ) = f ( x ) + 2 .
Für die Funktion G gilt G ( x ;   f ( x ) ) = f ( x ) + 2  bzw G ( x i ;   y ¯ i ) = y ¯ i + 2 .

  • Polygonzugverfahren:
    y ¯ i + l ( p o l y ) = y ¯ i ( p o l y ) + m i ( p o l y ) h
  • RUNGE- KUTTA -Verfahren:
    y ¯ i + l ( r u -ka ) = y ¯ i ( r u -ka ) + m i ( r u -ka ) h

Die Berechnungsergebnisse für beide Verfahren sowie die Werte y i der exakten Lösung sind in folgender Tabelle enthalten.

Bild

Sowohl die Tabelle als auch die grafische Darstellung deuten an, dass im Allgemeinen numerische Lösungen nach dem RUNGE-KUTTA-Verfahren kleinere Abweichungen gegenüber der exakten Lösung aufweisen als numerische Lösungen nach dem Polygonzugverfahren.

Vergleich des Polygonzugverfahrens mit dem RUNGE-KUTTA-Verfahren

Vergleich des Polygonzugverfahrens mit dem RUNGE-KUTTA-Verfahren

Stand: 2010
Dieser Text befindet sich in redaktioneller Bearbeitung.

Lexikon Share
Mathe Note verbessern?
 

Kostenlos bei Duden Learnattack registrieren und ALLES 48 Stunden testen.

Kein Vertrag. Keine Kosten.

  • 40.000 Lern-Inhalte in Mathe, Deutsch und 7 weiteren Fächern
  • Hausaufgabenhilfe per WhatsApp
  • Original Klassenarbeiten mit Lösungen
  • Deine eigene Lern-Statistik
  • Kostenfreie Basismitgliedschaft

Einloggen