Método de Runge-Kutta

El método de Runge-Kutta es un método genérico de resolución numérica de ecuaciones diferenciales. Este conjunto de métodos fue inicialmente desarrollado alrededor del año 1900 por los matemáticos C. Runge y M. W. Kutta.

Contenido

1 Descripción
- 1.1 Ejemplo
- 1.2 Variantes
2 Métodos de Runge-Kutta
- 2.1 Métodos de Runge-Kutta de cuarto orden
3 Véase también
4 Referencias
5 Enlaces externos

Descripción

Los métodos de Runge-Kutta (RK) son un conjuntos de métodos iterativos (implícitos y explícitos) para la aproximación de soluciones de ecuaciones diferenciales ordinarias, concretamente, del problema de valor inicial.

Sea

$y'(t) = f(t, y(t)) \,$

una ecuación diferencial ordinaria, con $f: \Omega \subset \mathbb{R} \times \mathbb{R}^n \to \mathbb{R}^n$ donde $\Omega \,$ es un conjunto abierto, junto con la condición de que el valor inicial de ƒ sea

$(t_0, y_0) \in \Omega.$

Entonces el método RK (de orden s) tiene la siguiente expresión, en su forma más general:

$y_{n+1} = y_n + h\,\sum_{i=1}^s b_ik_i$ ,

donde h es el paso por iteración, o lo que es lo mismo, el incremento $Δ t n$ entre los sucesivos puntos $t n$ y $t n + 1$ . Los coeficientes $k i$ son términos de aproximación intermedios, evaluados en ƒ de manera local

$k_i = f \left( t_n + h\, c_i\, , y_n + h\,\sum_{j=1}^s a_{ij}k_j \right ) \quad i=1,...,s.$

con $a i j, b i, c i$ coeficientes propios del esquema numérico elegido, dependiente de la regla de cuadratura utilizada. Los esquemas Runge-Kutta pueden ser explícitos o implícitos dependiendo de las constantes $a i j$ del esquema. Si esta matriz es triangular inferior con todos los elementos de la diagonal principal iguales a cero; es decir, $a i j = 0$ para $j = i,..., s$ , los esquemas son explícitos.

Ejemplo

Esquema Runge-Kutta de dos etapas, una en $t = t n$ y otra en $t = t n + Δ t n$ . ƒ(t,y(t)) en la primera etapa es:

$f_n=k_1 = f(t_n, y_n) \,$

Para estimar ƒ(t,y) en $t = t n + Δ t n$ se usa un esquema Euler

$f_{n+1}=k_2=f( t_n + \Delta t_n\, , y_n+\Delta t_n k_1 ). \,$

Con estos valores de ƒ, se sustituyen en la ecuación

$y_{n+1} = y_n + \int_{t_n}^{t_{n+1}} f(t,y(t))\,dt,$

de manera que se obtiene la expresión:

$y_{n+1}=y_{n} + {{\Delta t_n}\over 2} (k_1 + k_2).$

Los coeficientes propios de este esquema son: $b 1 = b 2 = 1 / 2; a 21 = 1; c 2 = 1.$

Variantes

Existen variantes del método de Runge-Kutta clásico, también llamado Runge-Kutta explícito, tales como la versión implícita del procedimiento o las parejas de métodos Runge-Kutta (o métodos Runge-Kutta-Fehlberg).

Este último consiste en ir aproximando la solución de la ecuación mediante dos algoritmos Runge-Kutta de órdenes diferentes, para así mantener el error acotado y hacer una buena elección de paso.

Métodos de Runge-Kutta

El método de Runge-Kutta no es sólo un único método, sino una importante familia de métodos iterativos, tanto implícitos como explícitos, para aproximar las soluciones de ecuaciones diferenciales ordinarias (E.D.O´s); estas técnicas fueron desarrolladas alrededor de 1900 por los matematicos alemanes Carl David Tolmé Runge y Martin Wilhelm Kutta.

Métodos de Runge-Kutta de cuarto orden

Un miembro de la familia de los métodos Runge-Kutta es usado tan comúnmente que a menudo es referenciado como “RK4” o como “el método Runge-Kutta”.

Definamos un problema de valor inicial como:

$y' = f(x, y), \quad y(x_0) = y_0$

Entonces el método RK4 para este problema está dado por la siguiente ecuación:

$y_{i+1} = y_i + {1 \over 6}\left ( k_1 + 2k_2 + 2k_3 + k_4 \right ) h$

Donde

$\begin{cases} k_1 & = f \left( x_i, y_i \right) \\ k_2 & = f \left( x_i + {1 \over 2}h, y_i + {1 \over 2}k_1 h \right) \\ k_3 & = f \left( x_i + {1 \over 2}h, y_i + {1 \over 2}k_2 h \right) \\ k_4 & = f \left( x_i + h, y_i + k_3h \right) \\ \end{cases}$

Así, el siguiente valor (y_n+1) es determinado por el presente valor(y_n) más el producto del tamaño del intervalo (h) por una pendiente estimada. La pendiente es un promedio ponderado de pendientes, donde $k 1$ es la pendiente al principio del intervalo, $k 2$ es la pendiente en el punto medio del intervalo, usando $k 1$ para determinar el valor de y en el punto $\scriptstyle x_n + \frac{h}{2}$ usando el método de Euler. $k 3$ es otra vez la pendiente del punto medio, pero ahora usando $k 2$ para determinar el valor de y $k 4$ es la pendiente al final del intervalo, con el valor de y determinado por $k 3$ . Promediando las cuatro pendientes, se le asigna mayor peso a las pendientes en el punto medio:

$\mbox{pendiente} = \frac{k_1 + 2k_2 + 2k_3 + k_4}{6}.$

Esta forma del método de Runge-Kutta, es un método de cuarto orden lo cual significa que el error por paso es del orden de $O (h 5)$ , mientras que el error total acumulado tiene el orden $O (h 4)$ .

Véase también

Referencias

Ascher, Uri M.; Petzold, Linda Ruth (1998) (en inglés). Computer methods for ordinary differential equations and differential-algebraic equations (1ª edición). Philadelphia (USA): SIAM. ISBN 0898714125.

Enlaces externos

Weisstein, Eric W. «Runge-Kutta Method» (en inglés). MathWorld. Wolfram Research.
TEMA 1 Problemas de Valor Inicial para Ecuaciones Diferenciales Ordinarias: Métodos Numéricos de un paso.
Explicación gráfica del Método Runge Kutta

Categoría:

Ecuaciones diferenciales numéricas

Wikimedia foundation. 2010.

Игры ⚽ Поможем сделать НИР

Mira otros diccionarios:

Método de Runge-Kutta — El método de Runge Kutta es un método genérico de resolución numérica de ecuaciones diferenciales. Este conjunto de métodos fueron desarrollados alrededor del año 1900 por los matemáticos C. Runge y M. W. Kutta. Se trata de un método por etapas… … Enciclopedia Universal
Método de Dormand-Prince — En análisis numérico, Dormand Prince es un método para la resolución de ecuaciones diferenciales ordinarias. Pertenece a la familia de métodos Runge Kutta. Evalúa seis veces la función para calcular las soluciones de cuarto y quinto orden. La… … Wikipedia Español
Martin Wilhelm Kutta — Martin W. Kutta Martin Wilhelm Kutta Nacimiento 3 de noviembre de 1867 Pitschen, Alta Silesia … Wikipedia Español
Carle David Tolmé Runge — Saltar a navegación, búsqueda Carle David Tolmé Runge Carle David Tolmé Runge (30 de agosto de 1856 – 3 de enero de 1927) fue un matemático, físico y espectroscopista alemán. Su primer nombre suele escr … Wikipedia Español
Integración numérica — En análisis numérico, la integración numérica constituye una amplia gama de algoritmos para calcular el valor numérico de una integral definida y, por extensión, el término se usa a veces para describir algoritmos numéricos para resolver… … Wikipedia Español
Aproximación — es una representación inexacta que, sin embargo, es suficientemente fiel como para ser útil. Aunque en matemáticas la aproximación típicamente se aplica a números, también puede aplicarse a objetos tales como las funciones matemáticas, figuras… … Wikipedia Español
Computación Científica — La Computación Científica (o Ciencia Computacional) es el campo de estudio relacionado con la construcción de modelos matemáticos y técnicas numéricas para resolver problemas científicos, de ciencias sociales y problemas de ingeniería.… … Wikipedia Español
Ecuación diferencial ordinaria — Saltar a navegación, búsqueda En matemáticas, una ecuación diferencial ordinaria (comúnmente abreviada EDO ) es una relación que contiene funciones de una sola variable independiente, y una o más de sus derivadas con respecto a esa variable. Las… … Wikipedia Español
Análisis numérico — El análisis numérico o cálculo numérico es la rama de las matemáticas que se encarga de diseñar algoritmos para, a través de números y reglas matemáticas simples, simular procesos matemáticos más complejos aplicados a procesos del mundo real. El… … Wikipedia Español

Los diccionarios y las enciclopedias sobre el Académico

Método de Runge-Kutta

Contenido