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 der Stabilitätsfunktion eines Verfahrens

mit

ist das Stabilitätsgebiet daher gegeben durch:

Besonders interessant ist der Fall, wenn das Stabilitätsgebiet die komplette linke Halbebene enthält, dann heißt das 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