11.08.2013 Views

Apunte

Apunte

Apunte

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

3. MÉTODOS MULTIPASO LINEALES 185<br />

3. Métodos multipaso lineales<br />

Hasta ahora para aproximar las soluciones de x ′ (t) = f(t, x(t)) nos basamos en el punto inmediato<br />

anterior para calcular el siguiente, es decir, empezando en x0 las iteraciones de un paso son<br />

de la forma<br />

xn+1 = xn + hΦ(tn, xn, h).<br />

La filosofía de los métodos multipaso es usar la información obtenida hasta el momento (varias<br />

aproximaciones anteriores a la que se quiere calcular). Es decir, se quiere aproximar x en el paso<br />

(n + 1)-ésimo, x(tn+1), utilizando algunos de los valores ya calculados como xn, xn−1, xn−2, . . .<br />

Si se considera un registro de k datos anteriores el método tiene la forma:<br />

k−1<br />

αkxn+k = − αjxn+j + h<br />

Un método como en (8.21) se conoce como método multipaso de k pasos.<br />

j=0<br />

Por ejemplo, si aproximamos la derivada por<br />

x ′ (t) ∼<br />

considerando puntos equiespaciados nos queda<br />

y como<br />

podemos poner<br />

es decir<br />

k<br />

j=0<br />

x(t + h) − x(t − h)<br />

2h<br />

x ′ (tn) ∼ xn+1 − xn−1<br />

2h<br />

x ′ (tn) = f(tn, x(tn)) ∼ f(tn, xn)<br />

xn+1 − xn−1<br />

2h<br />

= f(tn, xn)<br />

xn+1 = xn−1 + 2hf(tn, xn).<br />

βjf(tn+j, xn+j) (8.21)<br />

Éste, resulta ser un método de dos pasos puesto que para calcular xn+1 se usan xn y xn−1.

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

Saved successfully!

Ooh no, something went wrong!