Regla trapezoidal (ecuaciones diferenciales)

En análisis numérico y computación científica, la regla trapezoidal es un método numérico de resolución de ecuaciones diferenciales ordinarias derivado de la regla del trapecio para calcular integrales. La regla trapezoidal es un método implícito de segundo orden, que puede considerarse tanto como un método de Runge-Kutta como un método lineal multipaso.

Método

Supongase que se quiere resolver la ecuación diferencial:

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

La regla trapezoidal viene dada por la fórmula

$y_{n+1}=y_{n}+{\tfrac {1}{2}}h{\Big (}f(t_{n},y_{n})+f(t_{n+1},y_{n+1}){\Big )},$

donde $h=t_{n+1}-t_{n}$ es el tamaño del paso.​

Este es un método implícito: el valor $y_{n+1}$ aparece en ambos lados de la ecuación, y para calcularlo realmente, hay que resolver una ecuación que será generalmente no lineal. Un método posible para resolver esta ecuación es el método de Newton. Se puede utilizar el método de Euler para obtener una estimación bastante buena para la solución, que puede utilizarse a su vez como valor inicial del método de Newton.​

Justificación

Integrando la ecuación diferencial de $t_{n}$ a $t_{n+1}$, se tiene que

$y(t_{n+1})-y(t_{n})=\int _{t_{n}}^{t_{n+1}}f(t,y(t))\,\mathrm {d} t.$

La regla del trapecio indica que la integral en el lado derecho puede ser aproximada como

$\int _{t_{n}}^{t_{n+1}}f(t,y(t))\,\mathrm {d} t\approx {\tfrac {1}{2}}h{\Big (}f(t_{n},y(t_{n}))+f(t_{n+1},y(t_{n+1})){\Big )}.$

Ahora, combinando ambas fórmulas y utilizando $y_{n}\approx y(t_{n})$ y $y_{n+1}\approx y(t_{n+1})$ se obtiene la regla trapezoidal para resolver ecuaciones diferenciales ordinarias.​

Análisis de errores

Se deduce del análisis de error de la regla trapezoidal para cuadratura que el error local de truncado $\tau _{n}$ de la regla trapezoidal para resolver ecuaciones diferenciales puede ser delimitado como:

$|\tau _{n}|\leq {\tfrac {1}{12}}h^{3}\max _{t}|y'''(t)|.$

Por lo tanto, la regla trapezoidal es un método de segundo orden. Este resultado se puede utilizar para demostrar que el error global es $O(h^{2})$ cuando el tamaño del paso $h$ tiende a cero (véase cota superior asintótica para el significado de esta afirmación).​

Estabilidad

En color rosa, la región de estabilidad del método trapezoidal.

La región de absoluta estabilidad para la regla trapezoidal viene dada por:

$\{z\in \mathbb {C} \mid \operatorname {Re} (z)<0\}.$

Esto incluye el semiplano izquierdo, por lo que la regla trapezoidal es A-estable. La segunda barrera de Dahlquist establece que la regla trapezoidal es la más precisa entre los métodos lineales multipaso estable. Más precisamente, un método lineal de varios pasos que es A-estable tiene como máximo el orden dos, y la constante de error de un método lineal multipaso estable de segundo orden no puede ser mejor que la constante de error de la regla trapezoidal.​

De hecho, la región de estabilidad absoluta para la regla trapezoidal es precisamente el semiplano izquierdo. Esto significa que si la regla trapezoidal se aplica a la ecuación de prueba lineal y = λyy, la solución numérica tiende a cero si y solo si la solución exacta también lo hace.

Referencias

• Iserles, Arieh (1996), A First Course in the Numerical Analysis of Differential Equations, Cambridge University Press, ISBN 978-0-521-55655-2 ..
• Süli, Endre; Mayers, David (2003), An Introduction to Numerical Analysis, Cambridge University Press, ISBN 0521007941 ..

Véase también

• Método de Crank-Nicolson
• Método de Runge-Kutta

• Datos: Q7835583