Métodos Multipaso
Dr. César Gutiérrez Cuba
Dr. César Gutiérrez Cuba
Métodos Predictor-Corrector
Los métodos de Euler, RK se llaman métodos de paso simple o
métodos de un solo paso, porque en el cálculo de cada punto sólo
se utiliza la información del punto previo, esto es, sólo se utiliza el
punto (t0,y0) para calcular (t1,y1). En general sólo se necesita
conocer yn para calcular yn+1.
Una vez calculados varios puntos, éstos pueden ser tomados en
cuenta para calcular yn+1, esto permite construir mejores métodos
de aproximación que se llaman métodos de varios pasos o
métodos multipaso.
Métodos de un solo paso: Sólo tienen en cuenta lo que
ocurrió en el paso anterior
Método Multipaso: Tienen en cuenta lo ocurrido en varios
pasos anteriores
Dr. César Gutiérrez Cuba
Métodos predictores-correctores
• En los métodos predictor-corrector, el valor inicial de y se estima
utilizando una fórmula predictiva; entonces la fórmula correctora se
usa para mejorar este valor. Por ejemplo, el método de Heun:
La aproximación inicial yi+1(0) se deriva por el método de Euler (fórmula predictiva), y este
valor se mejora mediante el uso del método de Heun (fórmula correctora). En esta
sección, derivaremos dos métodos predictores-correctores importantes conocidos como
métodos de Milne y Adams.
Los métodos de Milne y Adams requieren de valores en puntos anteriores para calcular
los valores en puntos sucesivos.
Dr. César Gutiérrez Cuba
Método de Milne (Método de Milne-Simpson)
• Este método requiere tres valores y1, y2, y3:
que deben calcularse utilizando cualquier otro método como RK.
El algoritmo de Milne está dado por:
Donde:
Dr. César Gutiérrez Cuba
Procedimientos de cálculo método de Milne:
1) Para resolver:
Primero calcular los valores por ejemplo con un método RK4.
2) Usar la fórmula predictora para calcular el valor de
El valor de puede ser calculado en el punto , entonces:
3) Emplear la fórmula correctora:
Para corregir el valor de y4.
Dr. César Gutiérrez Cuba
4) Calcular valores empleando la fórmula generalizada:
Donde:
Dr. César Gutiérrez Cuba
Ejercicio:
• Calcule los valores de y(0.1), y(0.2) e y(0.3) con el método de cuarto orden
Runge – Kutta para el PVI:
• Use el método de Milne para calcular y(0.5) con un tamaño de paso h = 0.1
Dr. César Gutiérrez Cuba
Sol:
Empleando el método RK4 al PVI de valor inicial, se obtienen los siguientes resultados:
Ahora, calcularemos y4(0.4) e y5(0.5) usando el método de Milne.
Dr. César Gutiérrez Cuba
Dr. César Gutiérrez Cuba
Dr. César Gutiérrez Cuba
Dr. César Gutiérrez Cuba
Dr. César Gutiérrez Cuba
Ejercicio:
• Utilice el método predictor-corrector de Milne para calcular y(0.4) a
partir de la ecuación diferencial:
y los siguientes valores:
Dr. César Gutiérrez Cuba
Sol:
Dr. César Gutiérrez Cuba
Tarea 2:
• Resolver:
𝑑𝑦
= 𝑦 − 2𝑥 2 + 1, y(0)=0,5
𝑑𝑥
• Calcular y(1,5) con tamaño de paso h=0,2 empleando el método de
Milne.
Dr. César Gutiérrez Cuba
Método de Adams (Adams-Bashfort predictor
y Adams-Moulton corrector)
• El algoritmo de Adams está dado por las ecuaciones:
Dónde:
Dr. César Gutiérrez Cuba
Procedimientos de cálculo método de Adams:
1) Para resolver el PVI:
Calcular los valores de y empleando método RK4.
2) Use la fórmula predictora de Adams-Bashforth para predecir la aproximación
inicial .
El valor de la función se puede calcular en el punto :
Dr. César Gutiérrez Cuba
3) Utilice la fórmula correctora Adams-Moulton:
Para corregir el valor de y1 a la exactitud deseada
4) La siguiente fórmula generalizada puede ser utilizada para calcular
los valores de y2, y3, …..:
Donde:
Dr. César Gutiérrez Cuba
Ejercicio:
• Para el PVI:
Use el método Adams para calcular y(0.5) con el tamaño del paso, h = 0.1.
Calcule los valores iniciales y(0.1), y(0.2) e y(0.3) con el método de cuarto
orden Runge – Kutta.
Dr. César Gutiérrez Cuba
Sol:
Empleando RK4:
Ahora, calcularemos y1(0.4) e y2(0.5) usando el método Adams.
Dr. César Gutiérrez Cuba
Dr. César Gutiérrez Cuba
Dr. César Gutiérrez Cuba
Dr. César Gutiérrez Cuba
• Usando la fórmula correctora Adams-Moulton, obtenemos:
Dr. César Gutiérrez Cuba
Ejercicio:
• Utilice el método predictor-corrector de Adams para calcular y(0.4) a
partir de la ecuación diferencial:
y los siguientes valores:
Dr. César Gutiérrez Cuba
Sol:
Ing. César Gutiérrez Cuba
Tarea
Resolver el siguiente PVI:
𝑑𝑦 2
= −2𝑦 + 𝑒 −2 𝑡−6
𝑑𝑡
𝑦 0 =1
Empleando método multipaso de Adams-Moulton para el rango 0≤
𝑡 ≤10, con tamaño de paso h=0,1.
Grafique sus resultados.
Dr. César Gutiérrez Cuba