0% encontró este documento útil (0 votos)
57 vistas5 páginas

Ts 3

Descargar como pdf o txt
Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1/ 5

Tema # 3

MÉTODOS IMPLÍCITOS

MÉTODO DE EULER MODIFICADO O CRANK-NICOLSON

Otra forma de usar las técnicas de extrapolación para la estimación


de es incluir el punto entre los puntos por los que debe pasar
el polinomio de extrapolación.

La diferencia esencial con el método de Euler es que debemos resolver una


ecuación generalmente no lineal para encontrar .

Figura 1 Método de Euler Modificado

El método de Euler modificado tiene dos motivaciones. La primera es que


es más preciso que el anterior. La segunda es que es más estable.

Este método se obtiene al aplicar la regla del trapecio para integrar


.

(1)

Si es lineal en y, la ecuación (1) se puede resolver fácilmente en


términos de . Por ejemplo, si la EDO está dada por:

La ecuación (1) queda:


Por lo tanto, al despejar se obtiene:

Si es una función no lineal de y, la ecuación (1) es una función no lineal


de , por lo que se requiere un algoritmo para resolver ecuaciones no lineales.
Uno de los métodos ampliamente utilizados para resolver ecuaciones no lineales
es el de la sustitución sucesiva.

(2)

Donde es la k-ésima aproximación iterativa de , y es una


estimación inicial de . Esta iteración se detiene si es menor
que una cierta tolerancia preestablecida. La estimación inicial es igual
a .Entonces, el primer paso de iteración es idéntico al método de Euler hacia
adelante. En el caso de que sólo se utilice un paso más de iteración, el esquema
se convierte en el método de Runge-Kutta de segundo orden o, en forma
equivalente, en el método predictor-corrector de Euler. Pero en el método de
Euler modificado, la iteración sigue hasta satisfacer la tolerancia.

El siguiente ejemplo muestra una aplicación del método de Euler


modificado al caso de una EDO no lineal de primer orden.

Ejemplo

Resuelva la siguiente ecuación mediante el método de Euler modificado,


con h = 0.1:
Con . Imprima los resultados hasta t=1

Solución:

El método de Euler modificado es:

Para n=0:

es la mejor estimación del lado derecho para . Sustituimos


del lado derecho y obtenemos:

Incorporamos este valor de en la ecuación (A) y tenemos:

Repetimos la sustitución unas veces más y obtenemos:


Los resultados calculados para diez pasos son:

T y
0 10
0,1 7,4932
0,2 5,8586
0,3 4,7345
0,4 3,9298
0,5 3,3357
0,6 2,8859
0,7 2,5386
0,8 2,2658
0,9 2,0487
1 1,8738

¿Por qué es mejor la precisión del método modificado con respecto al método de
Euler hacia adelante? Para dar una explicación analítica a esta pregunta, conside-
remos una ecuación de prueba, Entonces, la ecuación para este problema
es:

O bien, al despejar ,

Desarrollamos los coeficientes de esa ecuación y tenemos:

Al hacer la comparación de este desarrollo con la serie de Taylor del valor


exacto , tenemos que la ecuación es exacta hasta el
término de segundo orden. Así, el método de Euler modificado es un método
(preciso) de segundo orden. Por otro lado, un análisis similar del método de
Euler hacia adelante muestra que tiene una precisión de primer orden.

El error local (generado en cada paso) del método de Euler hacia adelante
es proporcional a , mientras que el error global es proporcional a h; en tanto
que el error local del método de Euler modificado es proporcional a y su error
global es proporcional a . El orden del error del método de Euler hacia atrás es
igual al del método de Euler hacia adelante.

Al aplicar el método de Euler modificado a un conjunto de EDO, la


ecuación debe resolverse en forma simultánea o "implícita". Sin embargo, la
ventaja de la solución implícita consiste en el hecho de que el método es más
estable que el método de Euler hacia adelante, por lo que permite un mayor
intervalo de tiempo.

MÉTODOS DE PREDICCIÓN – CORRECCIÓN

Combinación de fórmulas explícitas e implícitas para mejorar la


extrapolación y evitar los problemas numéricos de los métodos implícitos.

Ejemplo:

Método predictor-corrector de Euler

Predicción:

Corrección:

Nota: Puede iterarse con el corrector.

También podría gustarte