Método de Euler y Runge Kutta

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 6

UNIVERSIDAD CENTRAL DEL ECUADOR

FACULTAD DE INGENIERÍA QUÍMICA

ANÁLISIS NUMÉRICO

NOMBRES: BAUTISTA MELANY, LARA CARLA, SALAZAR ABIGAIL

CURSO: P2

FECHA: 2020-02-04

Soluciones numéricas para Ecuaciones Diferenciales Ordinarias

 Método de Euler

Para utilizar el Método de Euler con el fin de generar una solución numérica a un problema de
valor inicial de la forma:
y ' =f ( x , y )
y ( x 0 )= y 0

Primero se decide en qué intervalo, a partir de la condición inicial, deseamos encontrar la


solución. Cortamos este intervalo en pequeñas subdivisiones de longitud h.
Luego, usando la condición inicial como punto de partida, generamos el resto de la solución
usando las fórmulas iterativas:
x n+1=x n +h
y n+ 1= y n +hf ( x n , y n)

para así encontrar las coordenadas de los puntos de nuestra solución numérica. Terminamos
este proceso cuando hemos llegado al final del intervalo deseado.
Desde el punto de vista geométrico, tenemos en definitiva que el Método de Euler aproxima a
la función solución por medio de una línea poligonal, la aproximación será tanto peor cuanto
mayor sea en número de pasos, es decir, cuanto más “lejos” nos encontremos del punto inicial
(x0, y0). Por otro lado, el error será evidentemente tanto mayor cuanto más grande sea el
“paso” del método, h.

 Método de 4º orden Runge-Kutta para resolver la ecuación diferencial

Dada la siguiente información;


dy
1. Una ecuación diferencial ordinaria que define el valor de en la forma x y y.
dx
2. El valor inicial de y, es decir, y(0)
Por lo tanto, se nos da a continuación.

dx
=f ( x , y ) , y ( 0 )= y 0
dy

La tarea es encontrar el valor de la función desconocida y en un punto dado x.


El método de Runge-Kutta encuentra el valor aproximado de y para una x dada.
Sólo las ecuaciones diferenciales ordinarias de primer orden pueden ser resueltas usando el
método de Runge Kutta de 4º orden.
A continuación, se muestra la fórmula utilizada para calcular el siguiente valor yn+1 a partir
del valor anterior yn. Los valores de n son 0, 1, 2, 3, ....(x - x0)/h.
Aquí h es la altura del paso y x n+1=x 0 +h
Un rango de intervalo más bajo significa más precisión.
K 1=hf ( x n , y n)
h k1
K 2=hf ( xn + , y n + )
2 2
h k2
K 3=hf ( x n+ , y n+ )
2 2
K 4=hf ( xn + h , y n + k 3)
k1 k2 k3 k 4 5
y n+ 1= y n + + + + +O(h )
6 3 3 6

La fórmula básicamente calcula el siguiente valor yn+1 usando el yn actual más el promedio
ponderado de cuatro incrementos.
 k1 es el incremento basado en la pendiente al principio del intervalo, utilizando y
 k2 es el incremento basado en la pendiente en el punto medio del intervalo, usando y +
hk1/2.
 k3 es de nuevo el incremento basado en la pendiente en el punto medio, usando y +
hk2/2.
 k4 es el incremento basado en la pendiente al final del intervalo, usando y + hk3.
El método es un método de cuarto orden, lo que significa que el error de truncamiento local es
del orden de O(h5), mientras que el error acumulado total es del orden de O(h4).

Resolución del Ejercicio planteado

ECUACIÓN:
dy
=0,1 √ y +0,4 x2
dx
con y ( 2 )=4
GRÁFICA:

Bibliografía
Barker, C. (2017). Numerical Methods for Solving Differential Equations EULER'S METHOD. San
Joaquin DELTA COLLEGE.

Butcher, J. (2005). Runge–Kutta methods for Differential Equations. The University of Auckland.

Ejercicio en clase
 Método de Runge- Kutta de cuarto orden
Resolver la siguiente ecuación diferencial:

f ( x , y )=0.1 √ y + 0.4 x2

Con las siguientes condiciones: y(2)= 4 y h= 0.1. Encontrar el valor para y(2.5).

Algoritmo en Matlab para el Método de Runge- Kutta de cuarto orden

clear all
clc
fprintf('\n \t ECUACÓN DIFERENCIAL POR EL MÉTODO DE RUNGE-KUTTA DE CUARTO
ORDEN\n')
f=input('\n Ingrese la ecuacion diferencial\n','s');
x0=2
x1=2.5
y0=4
n=6;
h=0.1;
xs=x0:h:x1;
disp('-------------------------');
fprintf('\n Nº x0 y(x1)\n');
disp('-------------------------');
for i=1:n
it=i-1;
x0=xs(i);
x=x0;
y=y0;
k1=h*eval(f);
x=x0+h/2;
y=y0+k1/2;
k2=h*eval(f);
x=x0+h/2;
y=y0+k2/2;
k3=h*eval(f);
x=x0+h;
y=y0+k3;
k4=h*eval(f);
y0=y0+(k1+2*k2+2*k3+k4)/6;
fprintf('\n%2.0f%10.6f%10.6f\n',it,x0,y0);
end
disp('-------------------------');
fprintf('\n El punto aproximado y(x1) es = %8.6f\n',y0);

Resultados
El valor aproximado para y(2.5)= 5.40600

Conclusiones
 El método de Runge-Kutta de cuarto orden, en análisis numérico es uno de los
métodos más usados para la resolución de ecuaciones diferenciales, mediante
la iteración se aproxima al valor del problema, haciendo uso del promedio de
las cuatro pendientes que posee la ecuación diferencial.

También podría gustarte