Euler-Verfahren ODE Rechner (Explizit, Implizit, Heun)
Euler-Verfahren-Rechner für ODE y' = f(t,y): explizites Euler, implizites, Heun und Mittelpunkt. Konvergenzordnung 1-2, logistische und Newton-Beispiele.
Berechnungsfehler melden
Hilf uns zu verbessern. Beschreibe, was an der Berechnung nicht stimmt.
Ergebnis
| Voce | Valore |
|---|
Alle Details anzeigen
Ergebnis teilen
Haftungsausschluss: Diese Berechnung ist rein indikativ. Für wichtige Entscheidungen konsultieren Sie einen Fachmann.
Euler-Verfahren für gewöhnliche Differentialgleichungen
Die Euler-Verfahren sind die einfachste Klasse von Algorithmen zum numerischen Lösen gewöhnlicher Differentialgleichungen (ODE) der Form y' = f(t,y), y(t_0) = y_0. Entwickelt von Leonhard Euler 1768 (Institutiones Calculi Integralis), sind sie der Ausgangspunkt der numerischen ODE-Analyse.
Explizites Euler (Forward Euler)
Formel: y_{n+1} = y_n + h f(t_n, y_n) wobei h = (T - t_0)/N der Schritt ist. Globale Konvergenzordnung 1: Fehler = O(h). Schnell aber ungenau. Instabil für steife Probleme.
Implizites Euler (Backward Euler)
Formel: y_{n+1} = y_n + h f(t_{n+1}, y_{n+1}). Erfordert Lösung impliziter Gleichung in jedem Schritt. Doppelte Kosten aber A-stabil: keine h-Beschränkung. Ideal für steife ODEs.
Heun-Verfahren
Prädiktor-Korrektor: expliziter Euler-Prädiktor, Trapezkorrektor. Ordnung 2: Fehler = O(h^2). Kosten: 2 Auswertungen pro Schritt. Optimaler Genauigkeit/Kosten-Kompromiss.
Mittelpunkt (Modified Euler)
k1 = f(t_n, y_n), y_{n+1} = y_n + h f(t_n + h/2, y_n + h/2 k1). Ordnung 2 wie Heun. Stabil für oszillierende ODEs. Variante des zweistufigen Runge-Kutta.
Konvergenz und Ordnung
Lokaler Trunkierungsfehler explizites Euler: -h^2/2 y''(xi). Summierung von N = (T-t0)/h Schritten: globaler Fehler O(h). Heun und Mittelpunkt: O(h^2). Klassisches RK4: O(h^4). Kompromiss: höhere Ordnung = mehr Auswertungen pro Schritt, größeres h möglich.
Stabilität
Für y' = lambda y, explizites Euler stabil wenn |1 + h lambda| < 1. Implizites A-stabil: immer. Für steife ODEs implizit essentiell.
Steife Probleme
Steife ODEs haben sehr unterschiedliche Zeitskalen. Explizit benötigt winziges h. Implizit hat diese Beschränkung nicht. Beispiel: Robertson-Gleichung in chemischer Kinetik.
Bemerkenswerte ODEs
- Exponentieller Zerfall: y' = -k y.
- Logistisches Wachstum: y' = r y (1 - y/K).
- Newton-Abkühlung: T' = -k (T - T_amb).
- Riccati: y' = y^2 - t.
- Lotka-Volterra: Räuber-Beute.
- Harmonischer Oszillator: y'' + y = 0.
Wann welches Verfahren
- Explizites Euler: didaktisch, nicht steif, kleines h akzeptabel.
- Implizites Euler: steif (Chemie, Strömungen, Gleichgewichte).
- Heun: Lehrcode, mittlere Genauigkeit.
- RK4: Praxisstandard für nicht steife ODEs.
- Adaptive (RK45): variable Geschwindigkeiten.
Geschichte
Euler 1768 führt Verfahren ein. Cauchy 1840 beweist rigoros Konvergenz. Heun 1900 und Kutta 1901 erweitern auf höhere Ordnungen. Dahlquist 1956 führt A-Stabilität ein. Moderne Software (LSODE, scipy.integrate.solve_ivp) enthält Dutzende spezialisierter Verfahren.
Lokaler und globaler Fehler: woher die Ordnung kommt
Die Unterscheidung zwischen lokalem und globalem Fehler erklart, warum das explizite Euler-Verfahren ein Verfahren erster Ordnung heisst. Der lokale Abbruchfehler ist der in einem einzelnen Schritt begangene Fehler, wobei der Startwert als exakt angenommen wird: entwickelt man die Losung in eine Taylor-Reihe, sieht man, dass Euler den Term (h^2/2) y'' vernachlassigt, sein lokaler Fehler also proportional zu h^2 ist. Der globale Fehler, akkumuliert uber das gesamte in etwa 1/h Schritte geteilte Intervall, ist hingegen proportional zu h: eine Halbierung des Schritts halbiert den Endfehler.
Verfahren zweiter Ordnung wie Heun und der Mittelpunkt heben auch den quadratischen Term auf und haben lokalen Fehler h^3 und globalen Fehler h^2: eine Halbierung von h reduziert den Fehler auf ein Viertel. Ein praktischer Weg, den Fehler ohne Kenntnis der exakten Losung zu schatzen, ist die Richardson-Extrapolation, die die mit zwei verschiedenen Schritten (h und h/2) erhaltenen Ergebnisse vergleicht, um die effektive Ordnung und eine quantitative Fehlerschatzung abzuleiten.
Auf dem Weg zu Runge-Kutta: die Leiter der Einschrittverfahren
Die Euler-Verfahren sind die ersten Sprossen einer breiteren Hierarchie, der der Runge-Kutta-Einschrittverfahren. Das explizite Euler ist das Runge-Kutta erster Ordnung; der Mittelpunkt und das Heun-Verfahren sind Varianten zweiter Ordnung, die die Steigung an Zwischenpunkten auswerten, um die Genauigkeit zu verdoppeln. Das beruhmte RK4 vierter Ordnung kombiniert vier Auswertungen pro Schritt und bietet einen ausgezeichneten Kompromiss zwischen Genauigkeit und Kosten und ist das Arbeitspferd der numerischen ODE-Losung.
Der nachste Schritt ist die adaptive Schrittweitensteuerung: Verfahren wie Runge-Kutta-Fehlberg (RKF45) berechnen gleichzeitig eine Schatzung vierter und funfter Ordnung und leiten aus ihrer Differenz den lokalen Fehler ab, wobei sie h automatisch anpassen, um ihn unter einer Schwelle zu halten. So verwendet der Loser grosse Schritte, wo die Losung glatt ist, und kleine Schritte, wo sie sich rasch andert, und optimiert die Zahl der Auswertungen: eine Idee, die direkt von den einfachen Euler-Verfahren abstammt.
Aktualisiert im Juli 2026
Kommentare (0)
Anmelden um zu kommentieren.
Noch keine Kommentare. Sei der Erste!