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.

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

Ahora,<br />

tn+3<br />

tn+1<br />

x ′ (t) dt =<br />

2<br />

donde a, b corresponden al corrimiento t = as + b:<br />

Es decir,<br />

0<br />

x ′ (as + b)a ds,<br />

<br />

s = 0 → t = tn+1 = tn + h<br />

s = 2 → t = tn+3 = tn + 3h ⇒<br />

<br />

a = h<br />

b = tn + h .<br />

tn+3<br />

tn+1<br />

x ′ (t) dt ∼ h[ 1<br />

3 x′ (tn) − 2<br />

3 x′ (tn+1) + 7<br />

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

Ahora, como x ′ (tn) = f(tn, x(tn)) se aproxima por fn = f(tn, xn) y además<br />

proponemos el método<br />

x(tn+3) − x(tn+1) ∼ h[ 1<br />

3 x′ (tn) − 2<br />

3 x′ (tn+1) + 7<br />

3 x′ (tn+2)],<br />

xn+3 − xn+1 = h[ 1<br />

3 fn − 2<br />

3 fn+1 + 7<br />

3 fn+2],<br />

que es un método de multipaso explícito de 3 pasos.<br />

3.1. Convergencia de los métodos multipaso. Empecemos estudiando el error de<br />

truncamiento local para un método de k pasos. Para esto, recordemos que en los métodos de 1<br />

paso de la forma<br />

xi+1 = xi + hΦ(ti, xi, h),<br />

el error local es τ definido por medio de la expresión<br />

x(t + h) − x(t) − hΦ(t, x(t), h) = hτ,<br />

donde x es la solución exacta del problema x ′ (t) = f(t, x(t)). Sabemos que el método resulta<br />

convergente cuando τ tiende a cero. Aplicando la misma idea para un método de k pasos de la<br />

forma<br />

k<br />

k<br />

αjxn+j = h<br />

j=0<br />

j=0<br />

βjfn+j<br />

teniendo en cuenta que x ′ = f, definimos τ, el error de truncamiento local, de la siguiente manera<br />

k<br />

αjx(t + jh) − h<br />

j=0<br />

k<br />

βjx ′ (t + jh) = hτ. (8.24)<br />

j=0

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

Saved successfully!

Ooh no, something went wrong!