11.08.2013 Views

Apunte

Apunte

Apunte

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

186 8. RESOLUCIÓN DE ECUACIONES DIFERENCIALES ORDINARIAS<br />

Los métodos de integración también nos proporcionan ejemplos de métodos multipaso. Por<br />

ejemplo, partimos del cálculo de la integral<br />

x(tn+2) − x(tn) =<br />

tn+2<br />

tn<br />

x ′ (s) ds,<br />

y aproximamos el valor de la integral por la regla de Simpson, obteniendo<br />

x(tn+2) − x(tn) ∼ h<br />

3 (x′ (tn) + 4x ′ (tn+1) + x ′ (tn+2)).<br />

Como estamos resolviendo la ecuación x ′ (t) = f(t, x(t)) podemos proponer el siguiente método<br />

xn+2 − xn = h<br />

3 (f(tn, xn) + 4f(tn+1, xn+1) + f(tn+2, xn+2))<br />

que es un método de dos pasos.<br />

En adelante, pondremos fn por f(tn, xn). Con esta notación, los métodos multipaso de k pasos<br />

admiten la escritura<br />

k<br />

αjxn+j = h<br />

j=0<br />

k<br />

βjfn+j, (8.22)<br />

donde h es la distancia entre dos nodos, αj y βj son coeficientes conocidos y αk = 0.<br />

Si βk = 0, el término βkfn+k aparece en la suma de la derecha con lo cuál el valor a calcular<br />

xn+k aparece involucrado en la ecuación en la evaluación de f. Por este motivo se dice que el<br />

método es implícito.<br />

Cuando βk = 0 el valor a encontrar xn+k se puede despejar de datos ya conocidos y el método<br />

se dice explícito.<br />

Con los ejemplos vistos hasta ahora tenemos que,<br />

xn+2 − xn = h<br />

3 (fn + 4fn+1 + fn+2) (8.23)<br />

es un método implícito, que se conoce como el método de Milne, mientras que el que resulta de<br />

aproximar la derivada<br />

es un método de dos pasos explícito.<br />

j=0<br />

xn+2 = xn + 2hfn+1<br />

Si usamos una fórmula de cuadratura basada en interpolación de Lagange para aproximar la<br />

integral<br />

tn+k<br />

tn+k−1<br />

x ′ (s) ds,

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!