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.

2. ANÁLISIS DEL ERROR Y CONVERGENCIA 177<br />

Ejemplo 8.5. Aplicar el método de Euler Modificado y el método de Runge-Kutta de orden 4<br />

dado en (8.15) para aproximar √ 2 como solución del problema<br />

x ′ (t) = 1<br />

2x ,<br />

x(1) = 1.<br />

considerando el intervalo [1, 2].<br />

En cualquier caso, la función a considerar es f(t, x) = 1<br />

2x . Para el método de Euler Modificado<br />

la sucesión de recurrencias correspondiente a (8.10) es<br />

xi+1 = xi + h<br />

2<br />

1<br />

2xi + h<br />

2<br />

1<br />

2xi<br />

= xi + h<br />

Para el método de orden 4 vamos a escribir K1, K2, K3 y K4 sin desarrollar completamente las<br />

expresiones.<br />

4x 2 i<br />

xi<br />

+ h.<br />

K1 = 1<br />

2x , K3<br />

1<br />

= ,<br />

2x + hK2<br />

1<br />

1<br />

K2 = , K4 = .<br />

2x + hK1<br />

2x + hK3<br />

En la Tabla 8.3 damos los valores obtenidos con ambos métodos para un paso h = 0,25. El gáfico<br />

correspondiente se tiene en la Figura 8.3.<br />

t x = √ t Runge-Kutta 2 Runge-Kutta 4<br />

1.000 1.000000 1.000000 1.000000<br />

1.250 1.118034 1.058824 1.119073<br />

1.500 1.224745 1.114734 1.226467<br />

1.750 1.322876 1.168117 1.325074<br />

2.000 1.414214 1.219278 1.416760<br />

Cuadro 8.3. Métodos de Runge-Kutta, x ′ = 1<br />

2x ; x(1) = 1; paso h = 0,25.<br />

2. Análisis del error y convergencia<br />

El objetivo de esta sección es analizar cuál es el error que se comete al utilizar los métodos<br />

numéricos para aproximar la solución de una ecuación de la forma (8.2). Para esto introducimos<br />

dos conceptos, el error de truncamiento local y el error global. El primero corresponde al error<br />

que se cometería al aproximar el valor de x(t+h) si se conociera el valor exacto de x(t), mientras<br />

que el segundo corresponde al error que se acumula al aplicar n pasos de un método, es decir<br />

x(tn) − xn.

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

Saved successfully!

Ooh no, something went wrong!