Runge-Kutta Rechner (RK2, RK3, RK4, RK45 Fehlberg)

Runge-Kutta-Rechner für ODE: RK2 Heun und Mittelpunkt, RK3 Kutta, RK4 klassisch (Gold-Standard), RK45 Fehlberg adaptiv mit eingebetteter Fehlerschätzung.

Berechnungsparameter
War dieser Rechner nützlich?

Runge-Kutta-Verfahren für ODE

Die Runge-Kutta-Verfahren sind eine Familie von Einschritt-Algorithmen zur numerischen Integration gewöhnlicher Differentialgleichungen y' = f(t,y). Entwickelt von Carl Runge (1895) und Wilhelm Kutta (1901), verallgemeinern sie das Euler-Verfahren auf mehrere Stufen für höhere Konvergenzordnungen.

Allgemeine Idee

Ein s-stufiges RK-Verfahren berechnet y_{n+1} = y_n + h * (b_1 k_1 + ... + b_s k_s) wobei k_i = f(t_n + c_i h, y_n + h sum_{j} a_{ij} k_j). Koeffizienten in "Butcher-Tableau" organisiert.

RK2 Heun (Ordnung 2)

Prädiktor-Korrektor: k1 = f(t_n, y_n), k2 = f(t_n + h, y_n + h k1), y_{n+1} = y_n + h (k1 + k2) / 2. Globaler Fehler O(h^2). 2 Auswertungen pro Schritt.

Klassisches RK4 (Ordnung 4) - Gold-Standard

4 Stufen: k1, k2 (Mitte mit k1), k3 (Mitte mit k2), k4 (Ende mit k3). y_{n+1} = y_n + h (k1 + 2 k2 + 2 k3 + k4) / 6. Fehler O(h^4). Standard in Computerphysik, Planetensimulationen, Molekulardynamik.

RK45 Fehlberg adaptiv

Berechnet simultan RK4 und RK5 mit gleichen 6 k-Werten. Fehlerschätzung: |y_5 - y_4|. Passt h automatisch an. Standard in scipy.integrate.RK45.

Konvergenz und Ordnung

Globaler Fehler eines Verfahrens der Ordnung p ist O(h^p). Halbieren von h reduziert Fehler um 2^p. RK4: Fehler /16 bei jeder Halbierung. Kompromiss: höhere Ordnung = mehr Auswertungen aber größeres h möglich.

Stabilität

RK4 endliches Stabilitätsgebiet in komplexer Ebene, nicht A-stabil. Für steife Probleme: implizite Verfahren (BDF, Radau, Rosenbrock).

Praktischer Vergleich

Für y' = -y, T = 10 (exakt = e^-10):

  • Explizites Euler h=0.1: Fehler 4%
  • RK2: Fehler 0.06%
  • RK3: Fehler 5e-5%
  • RK4: Fehler 4e-7%
  • RK45 adaptiv: ~5 Schritte

Butcher-Tableau

Standardnotation eingeführt von John Butcher (1963): Matrix A, Vektor b (Gewichte), Vektor c (Knoten). Für klassisches RK4: c=[0,1/2,1/2,1], b=[1/6,1/3,1/3,1/6].

Fortgeschrittene Verfahren

  • Cash-Karp, Dormand-Prince: eingebettete Varianten.
  • RK Gauss-Legendre: A-stabil, max. Ordnung.
  • Rosenbrock: linear implizit für steif.
  • Verlet: symplektisch für Hamiltonianer.

Geschichte

Runge 1895, Kutta 1901, Fehlberg 1969 (NASA), Dormand-Prince 1980. Monographie: Hairer-Norsett-Wanner "Solving ODEs" (1987-1996).

Die Ordnungsbedingungen und die Butcher-Schranke

Damit ein Runge-Kutta-Verfahren eine gegebene Genauigkeitsordnung erreicht, mussen seine Koeffizienten einen genauen Satz algebraischer Gleichungen erfullen, die Ordnungsbedingungen. Die einfachste, die Konsistenz, verlangt, dass die Summe der Gewichte b_i gleich 1 ist; fur Ordnung zwei kommt die Bedingung an die Knotenmomente hinzu, und fur hohere Ordnungen wachst die Zahl der Bedingungen rasch. John Butcher formalisierte diese Theorie in den 1960er Jahren, indem er jede Bedingung mit einem Wurzelbaum verknupfte, einem eleganten kombinatorischen Werkzeug, das die Einschrankungen zahlt und ordnet.

Bis Ordnung vier gibt es ein Verfahren mit genau so vielen Stufen wie die Ordnung, weshalb RK4 vier Auswertungen verwendet. Uber Ordnung vier hinaus greift jedoch die Butcher-Schranke: kein explizites Verfahren funfter Ordnung kann nur funf Stufen haben, mindestens sechs sind notig, und fur noch hohere Ordnungen wachst die Lucke. Diese theoretische Begrenzung erklart, warum RK4 ein so gunstiges Gleichgewicht zwischen Genauigkeit und Rechenaufwand darstellt.

Dormand-Prince und moderne Loser

Fehlbergs RK45-Verfahren war das erste eingebettete Schema, das populaer wurde, doch heute ist der De-facto-Standard das Dormand-Prince-Paar (DP54), das vom ode45-Loser von MATLAB und von der RK45-Methode von SciPy verwendet wird. Wie RKF45 berechnet es gleichzeitig eine Schatzung vierter und funfter Ordnung zur Fehlerkontrolle, doch seine Koeffizienten sind darauf optimiert, den Fehler der funften-Ordnung-Losung zu minimieren, die tatsachlich fortgepflanzt wird.

Dormand-Prince besitzt zudem die FSAL-Eigenschaft (First Same As Last): die letzte Auswertung eines Schritts stimmt mit der ersten des nachsten uberein und spart so eine Funktionsauswertung pro Schritt. Viele moderne Loser bieten auch dichte Ausgabe, eine Interpolation, die die Losung an jedem Punkt zwischen den Knoten liefert, ohne den Schritt zu verkleinern. Diese Verfeinerungen, allesamt Nachkommen der ursprunglichen Idee von Runge und Kutta, machen adaptive Loser zum Standardwerkzeug fur ODEs in der computergestutzten Wissenschaft.

Aktualisiert im Juli 2026

Häufig gestellte Fragen

Was sind Runge-Kutta-Verfahren?

Es ist eine Familie numerischer Verfahren zur Losung gewohnlicher Differentialgleichungen y' = f(t,y). In jedem Schritt werten sie die Steigung an mehreren Zwischenpunkten aus und bilden einen gewichteten Mittelwert, was bei gleichem Schritt eine weit hohere Genauigkeit als das einfache Euler-Verfahren ergibt.

Was ist das klassische RK4-Verfahren?

Es ist das Runge-Kutta vierter Ordnung, das in der Praxis meistverwendete: es kombiniert vier Steigungsauswertungen (k1 am Anfang, k2 und k3 zur Schrittmitte, k4 am Ende) im gewichteten Mittel (k1+2k2+2k3+k4)/6. Es bietet einen ausgezeichneten Kompromiss zwischen Genauigkeit und Kosten und gilt als Goldstandard.

Was bedeutet die Ordnung eines Verfahrens?

Die Ordnung gibt an, wie schnell der Fehler bei Verkleinerung des Schritts h abnimmt: ein Verfahren der Ordnung p hat globalen Fehler proportional zu h^p. Bei Halbierung von h reduziert RK2 den Fehler auf ein Viertel, RK4 auf ein Sechzehntel. Hohere Ordnungen geben viel mehr Genauigkeit mit weniger Schritten.

Was ist ein adaptives Verfahren wie RK45?

Es ist ein Verfahren, das den Schritt h wahrend der Rechnung automatisch anpasst. RK45 (Fehlberg oder Dormand-Prince) berechnet zwei Schatzungen verschiedener Ordnung (4 und 5) und leitet aus ihrer Differenz den lokalen Fehler ab, verkleinert den Schritt, wo die Losung sich rasch andert, und vergrossert ihn, wo sie glatt ist.

Was ist das Butcher-Tableau?

Es ist die kompakte Tabelle, die die Koeffizienten eines Runge-Kutta-Verfahrens zusammenfasst: die Zeitknoten, die Gewichte des Endmittels und die Koeffizienten, die die Zwischenstufen kombinieren. Von John Butcher eingefuhrt, erlaubt es, jedes RK-Verfahren zu beschreiben und seine Ordnungsbedingungen zu prufen.

Warum verwendet RK4 genau vier Stufen?

Weil bis Ordnung vier ein explizites Verfahren mit so vielen Stufen wie die Ordnung konstruierbar ist. Uber Ordnung vier hinaus greift die Butcher-Schranke: ein Verfahren funfter Ordnung braucht mindestens sechs Stufen, was RK4 zu einem optimalen Gleichgewicht zwischen Genauigkeit und Kosten macht.

Was ist der Unterschied zum Euler-Verfahren?

Das Euler-Verfahren ist das einfachste Runge-Kutta, erster Ordnung, und nutzt eine einzige Steigungsauswertung pro Schritt. RK-Verfahren hoherer Ordnung nutzen mehr an Zwischenpunkten und reduzieren den Fehler drastisch: RK4 ist bei gleichem Schritt typischerweise millionenfach genauer als Euler.

Welcher Loser wird in der Praxis verwendet?

In wissenschaftlicher Software ist der Standard das Dormand-Prince-Paar (DP54): der ode45-Loser von MATLAB und die RK45-Methode von SciPy verwenden es. Es ist ein adaptives Verfahren der Ordnung 4-5 mit automatischer Fehlerkontrolle und der FSAL-Eigenschaft, die eine Auswertung pro Schritt spart.

Kommentare (0)

Noch keine Kommentare. Sei der Erste!