Explizites Euler-Verfahren

Verfahren zur numerischen Lösung eines Anfangswertproblems

Das eulersche Polygonzugverfahren oder explizite Euler-Verfahren (auch Euler-Cauchy-Verfahren oder Euler-vorwärts-Verfahren) ist das einfachste Verfahren zur numerischen Lösung eines Anfangswertproblems.

Es wurde von Leonhard Euler 1768 in seinem Buch Institutiones Calculi Integralis vorgestellt. Cauchy benutzte es, um einige Eindeutigkeitsresultate für gewöhnliche Differentialgleichungen zu beweisen.

Das Verfahren wird manchmal in der Physik als Methode der kleinen Schritte bezeichnet.

Das Verfahren

Bearbeiten
 
Zwei Schritte des expliziten Euler-Verfahrens

Zur numerischen Lösung des Anfangswertproblems:

 

für eine gewöhnliche Differentialgleichung wähle man eine Diskretisierungs-Schrittweite  , betrachte die diskreten Zeitpunkte

 

und berechne die Werte

 

Die berechneten Werte   stellen Approximationen an die tatsächlichen Werte   der exakten Lösung des Anfangswertproblems dar. Je kleiner die Schrittweite   gewählt ist, desto mehr Rechenarbeit ist nötig, aber desto genauer werden auch die approximierten Werte.

Eine Modifikation des Verfahrens besteht hier darin, dass man die Schrittweite variabel wählt. Eine sinnvolle Veränderung der Schrittweite setzt einen Algorithmus zur Schrittweitensteuerung voraus, der den Fehler im aktuellen Schritt abschätzt und dann die Schrittweite für den nächsten Schritt dementsprechend wählt.

Wird ein Verfahren über   definiert, erhält man das implizite Euler-Verfahren. Dieses ist A-stabil und daher für steife Anfangswertprobleme besser geeignet.

Herleitung

Bearbeiten

Für die Herleitung von Einschrittverfahren wird das Anfangswertproblem meist in die dazu äquivalente Integralgleichung umgeformt[1]

 

Nun besteht die Idee, beim expliziten Euler-Verfahren eine simple Quadraturformel für das Integral zu benutzen: die linksseitige Boxregel. Man wählt in jedem  -ten Schritt den Integranden als konstanten Wert an der linken Intervallgrenze[2]

 

Eigenschaften

Bearbeiten
 
Stabilitätsgebiet des expliziten Euler-Verfahrens

Das explizite Euler-Verfahren hat Konsistenz- und Konvergenzordnung 1. Die Stabilitätsfunktion ist   und sein Stabilitätsgebiet daher das Innere des Kreises um −1 mit Radius 1 in der komplexen Zahlenebene.

Verbessertes explizites Euler-Verfahren

Bearbeiten

Anstatt die Boxregel für die numerische Integration zu verwenden, kann man auch die Mittelpunktsregel anwenden.

 

Nun wendet man wieder das explizite Euler-Verfahren zur Approximation von   an

 

Zusammen führt dies auf das verbesserte explizite Euler-Verfahren (oder Euler-Verfahren mit kleinerer Schrittweite)[3]

 

Verallgemeinerungen

Bearbeiten

Es lässt sich im Wesentlichen durch zwei verschiedene Ideen auf effizientere Verfahren verallgemeinern.

  • Die erste Idee ist, bei der Berechnung des nächsten Schrittes mehr als nur einen der zuvor berechneten Werte mit einzubeziehen. Auf diese Weise erhält man Verfahren höherer Ordnung in der Klasse der linearen Mehrschrittverfahren.
  • Die zweite Idee ist, bei der Berechnung des nächsten Schrittes die Funktion   auf dem Intervall   an mehreren Stellen auszuwerten. Auf diese Weise erhält man die Klasse der Runge-Kutta-Verfahren.

Die Klasse der allgemeinen linearen Verfahren bezieht beide Ideen der Verallgemeinerung mit ein und enthält die Klasse der linearen Mehrschrittverfahren sowie die Klasse der Runge-Kutta-Verfahren als Spezialfall.

Literatur

Bearbeiten
  • 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

Einzelnachweise

Bearbeiten
  1. Arnold Reusken: Numerik für Ingenieure und Naturwissenschaftler. Springer, Berlin 2006, ISBN 3-540-25544-3, S. 378.
  2. Arnold Reusken: Numerik für Ingenieure und Naturwissenschaftler. Springer, Berlin 2006, ISBN 3-540-25544-3, S. 381.
  3. Arnold Reusken: Numerik für Ingenieure und Naturwissenschaftler. Springer, Berlin 2006, ISBN 3-540-25544-3, S. 382.