37 votos

Explicación y Prueba de la cuarta orden de Runge-Kutta

Método de Runge-Kutte de 4º orden, el método es una técnica numérica utilizada para resolver diferenciales ordinarias ecuación de la forma dy/dx=f(x,y),y(0)=y0

Da yi+1 en forma yi+1=yi+(a1k1+a2k2+a3k3+a4k4)h, donde las constantes se encuentran: yi+1=yi+16(k1+2k2+2k3+k4)h y: k1=f(xi,yi) k2=f(xi+12,yi+12k1h) k3=f(xi+12,yi+12k2h) k4=f(xi+h,yi+k3h)

También he aprendido que estos se derivan de los primeros cuatro términos de la serie de Taylor: yi+1=yi+f(xi,yi)h+12!f(xi,yi)h2+13!f(xi,yi)h3+14!f(xi,yi)h4

Estas son las cosas que yo entendía. Sin embargo, he aquí lo que no puedo:

  • No puedo ver por qué es este método de "trabajo"?
  • ¿Cómo son las aiki términos de derivada y calcula a partir de la serie de Taylor?
  • Lo que justifica o prueba de este resultado?

Por favor, tenga en cuenta que he visto muchos artículos en estas cosas cosas, pero ninguno de ellos contenía una plena prueba alguna... también me gustaría mucho como para obtener algunos ejemplos en los que la función f es en realidad una función de x, y no sólo de x. Por favor me podrian enviar algunos enlaces con los más avanzados y concisa de los enlaces de esta pregunta?

57voto

snicker Puntos 123

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.

(1)x(t+Δt)=x(t)+(Δt)x(t)+(Δt)22!x(t)+higher order terms.(2)f(t+Δt,x+Δx)=f(t,x)+(Δt)ft(t,x)+(Δx)fx(t,x)+higher order terms.

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)+tt+hf(τ,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)+hi=1Nω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: (3)x(t+h)=x(t)+ω1K1+ω2K2.

Algunas notas:

  1. 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.
  2. 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.
  3. 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.

3voto

ValdaR Puntos 11

Le sugiero que revise los enlaces de abajo, usted puede encontrar documentos, videos, matlab (mathematica, maple) códigos y ejemplos sobre el método de rk

De 2º orden

4 de la orden

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X