Die Aussage von van der Houwen ist richtig, aber keine Aussage über alle Runge-Kutta-Methoden fünfter Ordnung. Die "Taylor-Polynome", auf die er sich bezieht, sind (wie Sie zu wissen scheinen) nur die Polynome des Grades , die sich annähern , um zu ordnen :exp ( z ) ppexp( z)p
Pp( z) = ∑j = 1pzjj !
Für das Polynom fünfter Ordnung zeigt sich, dass für kleines , so dass der Stabilitätsbereich einer Methode mit als Stabilitätspolynom keine Umgebung des Ursprungs enthält auf der imaginären Achse . Genau das sagt van der Houwen.ϵ P 5 ( z )| P5( i ϵ ) | > 1ϵP5( z)
Die wahrscheinlichste Ursache für Ihre Verwirrung ist, was mit "der Runge-Kutta-Methode fünfter Ordnung" gemeint ist. Es gibt (unendlich) viele Runge-Kutta-Methoden fünfter Ordnung, aber die bekanntesten haben nicht als Stabilitätspolynom. Warum? Wie John Butcher bekanntermaßen bewiesen hat , muss eine Runge-Kutta-Methode fünfter Ordnung mindestens sechs Stufen haben . Normalerweise hat das Stabilitätspolynom einer Methode mit sechs (oder mehr) Stufen den Grad sechs (oder mehr). Beispielsweise verwendet jede der auf dieser Wikipedia-Seite aufgelisteten Methoden fünfter Ordnung genau sechs Stufen und weist ein Stabilitätspolynom vom Grad sechs auf.P5( z)
Kann eine Methode fünfter Ordnung als Stabilitätspolynom haben? Ja; Dies würde eine explizite Extrapolationsmethode fünfter Ordnung (wie die in diesem Artikel beschriebenen bekannten ) tun. Es ist auch zu beachten, dass eine stufige Runge-Kutta-Methode mit dem Stabilitätspolynom für lineare ODEs eine Genauigkeit von 5 aufweist, jedoch nicht für nichtlineare ODEs.P5( z)pP5( z)
Schließlich ist es leicht, Fehler bei der Bestimmung des Umfangs des imaginären Stabilitätsintervalls für Runge-Kutta-Methoden höherer Ordnung zu machen. Dies liegt daran, dass die Grenze des Stabilitätsbereichs für solche Verfahren extrem nahe an der imaginären Achse liegt . Daher können Rundungsfehler zu falschen Schlussfolgerungen führen. Es sollten nur exakte Berechnungen verwendet werden (natürlich könnte die Relevanz der Stabilitätsbereichsgrenze für praktische Zwecke unter diesen Umständen mit Sicherheit diskutiert werden).
Zum Beispiel ist hier eine Darstellung des Stabilitätsbereichs der Methode fünfter Ordnung aus dem Fehlberg 5 (4) -Paar:
Das imaginäre Stabilitätsintervall ist leer, aber Sie können es bei dieser Auflösung nicht aus dem Bild erkennen! Es ist zu beachten, dass der Bereich eindeutig einen Teil der imaginären Achse enthält, jedoch keinen Abstand zum Ursprung.
In der Zwischenzeit ist hier der Plot für die Methode fünfter Ordnung aus dem Dormand-Prince 5 (4) -Paar:
Es hat ein imaginäres Stabilitätsintervall von ungefähr .[ - 1 , 1 ]
Für eine genaue Charakterisierung der Stabilitätsbereichsgrenze in der Nähe der imaginären Achse für (was ziemlich faszinierend ist!) Siehe meinen kürzlich erschienenen Aufsatz .Pp(z)
Sie könnten auch an dem NodePy-Paket interessiert sein , mit dem die obigen Darstellungen erstellt wurden und mit dem beispielsweise das imaginäre Stabilitätsintervall einer Methode genau bestimmt werden kann (Haftungsausschluss: Ich habe NodePy erstellt).