Implizites Euler-Verfahren
aus Wikipedia, der freien Enzyklopädie
Das implizite Euler-Verfahren (nach Leonhard Euler) ist das einfachste implizite Verfahren zur numerischen Lösung eines Anfangswert-Problems.
[Bearbeiten] Das Verfahren
Zur numerischen Lösung des Anfangswert-Problems:
für eine gewöhnliche Differentialgleichung wähle man eine Diskretisierungs-Schrittweite h > 0, betrachte die diskreten Zeitpunkte
und berechne die iterierten Werte
Der Wert f(tk + 1,xk + 1) ist hierbei nicht explizit gegeben, sondern nur implizit. Zu seiner Berechnung muss also in jedem Iterationsschritt ein (je nach Art der rechten Seite f) lineares oder nichtlineares Gleichungssystem gelöst werden.
Die berechneten Werte xk stellen dann Approximationen an die tatsächlichen Werte x(tk) der exakten Lösung des Anfangswert-Problems dar. Je kleiner man die Schrittweite h wählt, desto mehr Rechenarbeit hat man, aber desto besser werden auch die approximierten Werte.
Wird ein Verfahren über xk + 1 = xk + hf(tk,xk) definiert, erhält man das explizite Euler-Verfahren.
[Bearbeiten] Eigenschaften
Das implizite Euler-Verfahren hat Konsistenz- und Konvergenzordnung 1. Es ist A-stabil, sein Stabilitätsgebiet enthält also die komplette linke Halbebene der komplexen Zahlenebene. Es gibt damit für das implizite Euler-Verfahren keine Einschränkungen an die Zeitschritte aufgrund von Stabilitätseinschränkungen, was den Zwang des Lösens von Gleichungssystemen in jedem Schritt wettmacht. Aufgrund der geringen Ordnung ist es damit besonders für Probleme interessant, bei denen die Iteration in einen stabilen Endzustand hineinläuft und die Genauigkeit der Zwischenergebnisse nicht interessant ist.
[Bearbeiten] Literatur
- E. Hairer, S.P. Norsett, G. Wanner: Solving Ordinary Differential Equations I, Springer Verlag
- M. Hermann: Numerik gewöhnlicher Differentialgleichungen, Anfangs- und Randwertprobleme, Oldenbourg Verlag, München und Wien, 2004, ISBN 3-486-27606-9