Stabilitätsgebiet

In der Numerik gewöhnlicher Differentialgleichungen ist das Stabilitätsgebiet eines Verfahrens zur Lösung eines Anfangswertproblems definiert als die Menge der komplexen Zahlen mit und , für die das numerische Verfahren bei der Lösung der dahlquistschen Testgleichung

bei fester Schrittweite eine monoton fallende Folge von Näherungen liefert. Dies impliziert, dass das numerische Verfahren für diese Gleichung und diese Schrittweite stabil ist.

Mit Hilfe d​er Stabilitätsfunktion e​ines Verfahrens

mit

ist d​as Stabilitätsgebiet d​aher gegeben durch:

Besonders interessant i​st der Fall, w​enn das Stabilitätsgebiet d​ie komplette l​inke Halbebene enthält, d​ann heißt d​as Verfahren A-stabil.

Beispiele

Für das explizite Euler-Verfahren: , mit Stabilitätsfunktion: , liegt im Stabilitätsgebiet falls

womit das Stabilitätsgebiet allen Punkten im Inneren eines Kreises mit Mittelpunkt und Radius  in der komplexen Zahlenebene entspricht. Dieses Stabilitätsgebiet ist (wie auch jenes anderer expliziter Verfahren) beschränkt. Daher wird man, auch wenn dies die Genauigkeit gar nicht verlangen würde, zu einer Einschränkung der Schrittweite gezwungen. Dies hat zur Folge, dass implizite Verfahren beim Lösen steifer Differentialgleichungen bevorzugt werden.

Literatur

  • E. Hairer, G. Wanner: Solving Ordinary Differential Equations II, Stiff problems, Springer Verlag
This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. The authors of the article are listed here. Additional terms may apply for the media files, click on images to show image meta data.