Mtodo Heun
Mtodo Heun
Mtodo Heun
Se utiliza para
extrapolar linealmente
así, se combinan las dos pendientes para obtener una pendiente promedio en el
intervalo:
Este método consiste en una mejora del método de Euler para resolver
ecuaciones diferenciales de primer orden y conocido el valor inicial. En este caso,
lo que se realiza es un promedio entre el valor obtenido por Euler y otro obtenido a
partir de la aproximación del valor de la función en el punto siguiente, también por
Euler.
Yk+1 = Yk + h/2 [f (x k, Y k) + f (x k+1, Y k+1)]
Donde Y es solución de la ecuación diferencial, f es la ecuación diferencial en
función de las variables independientes y la solución de Y ˙ k+1 es una
aproximación de Euler.
Ejemplo Recordando la ecuación diferencial del ejercicio anterior:
i´(t) = sen (2πt) − i(t) = f(t, i)
Obtenemos la fórmula de iteración según el método de Heun:
I k+1 = I k + h /2 [ f (tk, Ik) + f (tk+1, ˙ I k+1)]
Nuevamente, los valores iniciales son t0 = 0, I0 = 0. En este caso, definiremos un
valor ˙ Ik+1 que corresponde a la aproximación de Euler en el paso siguiente.
Realizando los cálculos:
k=0
I1 = I0 + h · f (0; 0) = 0
I1 = I0 + h / 2 [f (0; 0) + f (0, 25; ˙I1)]= 0, 125 = i (0, 25)
k=1
I2 = I1 + h · f (0, 25; 0, 125) = 0, 3437
I2 = I1 + h / 2 [f (0, 25; 0, 125) + f (0, 5; ˙I2)] = 0, 1914 = i (0, 5)
k=2
I3 = I2 + h · f (0, 5; 0, 1914) = 0, 1435
I3 = I2 + h / 2 [ f (0, 5; 0, 1914) + f (0, 75; ˙I3)]= 0, 0245 = i (0, 75)
k=3
I4 = I3 + h · f (0, 75; 0, 0245) = −0, 2316
I4 = I3 + h / 2 [f (0, 75; 0, 245) + f (0, 75; ˙I 4)] = −0, 0745 = i(1)
Comparando con los valores teóricos:
2022.
2022.
2022.
2022.