Hay muchos métodos de Runge–Kutta. La que tú has descrito es probablemente la más popular y ampliamente utilizado. No voy a mostrar cómo derivar este método en particular – en lugar de eso voy a derivar la fórmula general para el explícito de segundo orden de Runge–Kutta y usted puede generalizar las ideas.
En lo que sigue vamos a necesitar los siguientes series de Taylor de las expansiones.
x(t+Δt)=x(t)+(Δt)x′(t)+(Δt)22!x′′(t)+higher order terms.f(t+Δt,x+Δx)=f(t,x)+(Δt)ft(t,x)+(Δx)fx(t,x)+higher order terms.(1)(2)
Estamos interesados en la siguiente ODA:
x′(t)=f(t,x(t)).
El valor de x(t) es conocido y x(t+h) es deseado. Podemos resolver esta ODA mediante la integración de:
x(t+h)=x(t)+∫t+htf(τ,x(τ))dτ.
Por desgracia, en realidad, haciendo que la integración exactamente a menudo es bastante difícil (o imposible), así que un aproximado usando cuadratura:
x(t+h)≈x(t)+hN∑i=1ωif(t+νih,x(t+νih)).
La exactitud de la cuadratura depende del número de términos en la suma (el orden de los métodos de Runge–Kutta), los pesos, ωi, y la posición de los nodos, νi.
Incluso esta cuadratura puede ser bastante difícil de calcular, ya que en el lado derecho tenemos a x(t+νih), lo que no sabemos. Hemos de conseguir alrededor de este problema de la siguiente manera:
Deje ν1=0, lo que hace que el primer término de la cuadratura K1:=hf(t,x(t)). Esto nos lo hacen saber y nosotros también lo puede utilizar para aproximar x(t+ν2h) escribiendo el segundo término de la cuadratura como K2:=hf(t+αh,x(t)+βK1).
Con esto, la cuadratura fórmula es:
x(t+h)=x(t)+ω1K1+ω2K2.(3)
Algunas notas:
- Si queríamos encontrar una de tercer orden, método, se introduciría K3=hf(t+~αh,x+~βK1+~γK2). Si queríamos un cuarto método de orden, se introduciría K4 en una manera similar.
- Este es un método explícito, es decir, hemos escogido K2 a depender de K1 pero K1 no depende de K2. Del mismo modo, K3 (si nos presentamos a) depende de la K1 K2 pero no dependen de K3. Se podría permitir que la dependencia a ejecutar ambas maneras, pero el método es implícita y se vuelve mucho más difícil de resolver.
- Todavía tenemos que elegir α, β y el ωi. Vamos a hacer eso ahora, usando la serie de Taylor de las expansiones me introdujo en el principio.
En la ecuación (3), sustituimos en la expansión en series de Taylor (1) en el lado izquierdo:
x(t)+hx′(t)+h22!x′′(t)+O(h3)=x(t)+ω1K1+ω2K2.
Desde x′=f e lo x′′=ft+ffx (supresión de los argumentos para la facilidad de notación), tenemos:
hf+h22(ft+ffx)+O(h3)=ω1K1+ω2K2.
Ahora sustituye por K1K2:
hf+h22(ft+ffx)+O(h3)=ω1hf+ω2hf(t+αh,x+βK1).
Taylor-ampliar en el lado derecho de usar (2):
hf+h22(ft+ffx)+O(h3)=ω1hf+ω2(hf+αh2ft+βh2ffx)+O(h3).
Así, el método de Runge–Kutta de acuerdo con la serie de Taylor aproximación a O(h3) si elegimos:
ω1+ω2=1,
αω2=12,
βω2=12.
La elección canónica para el segundo orden de Runge–Kutta es α=β=1 ω1=ω2=1/2.
El mismo procedimiento puede ser utilizado para encontrar las restricciones sobre los parámetros del cuarto orden de Runge–Kutta. La canónica de elección en este caso es el método que se describe en la pregunta.