Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 15 vistas

    2-Ec. Dif 1er Orden

    Cargado por

    José Ignacio Salgado
    Este documento presenta varios métodos numéricos para resolver ecuaciones diferenciales ordinarias de primer orden, incluyendo el método de Euler, método de Heun, método del punto medio, método de Ralston y métodos de Runge-Kutta de tercer, cuarto y quinto orden. Explica las fórmulas matemáticas utilizadas en cada método y provee ejemplos de código para implementarlos numéricamente.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd

    2-Ec. Dif 1er Orden

    Cargado por

    José Ignacio Salgado
    15 vistas5 páginas
    Este documento presenta varios métodos numéricos para resolver ecuaciones diferenciales ordinarias de primer orden, incluyendo el método de Euler, método de Heun, método del punto medio, método de Ralston y métodos de Runge-Kutta de tercer, cuarto y quinto orden. Explica las fórmulas matemáticas utilizadas en cada método y provee ejemplos de código para implementarlos numéricamente.

    Título original

    2-Ec. Dif 1er orden

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta varios métodos numéricos para resolver ecuaciones diferenciales ordinarias de primer orden, incluyendo el método de Euler, método de Heun, método del punto medio, método de Ralston y métodos de Runge-Kutta de tercer, cuarto y quinto orden. Explica las fórmulas matemáticas utilizadas en cada método y provee ejemplos de código para implementarlos numéricamente.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    15 vistas5 páginas

    2-Ec. Dif 1er Orden

    Cargado por

    José Ignacio Salgado
    Este documento presenta varios métodos numéricos para resolver ecuaciones diferenciales ordinarias de primer orden, incluyendo el método de Euler, método de Heun, método del punto medio, método de Ralston y métodos de Runge-Kutta de tercer, cuarto y quinto orden. Explica las fórmulas matemáticas utilizadas en cada método y provee ejemplos de código para implementarlos numéricamente.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    Está en la página 1de 5

    UNIVERSIDAD SANTIAGO DE CHILE

    FACULTAD DE INGENIER´IA - DEPARTAMENTO DE MECA


    ´ NICA ME´TODOS NUME´RICOS PARA INGENIER´IA MECA
    ´ NICA

    Ayudante: Ren´e Torres ([email protected])


    Profesor: Francisco Sepu´lveda ([email protected])

    APUNTE
    ECUACIONES DIFERENCIALEES ORDINARIAS (O DE PRIMER ORDEN)

    1. M´etodo de Euler
    Se tiene la siguiente ecuacion diferencial de primer orden de la forma
    dy
    = f (x, y)
    dx
    Es importante ver que dy

    representa la pendiente o la tasa de cambio de la variable y con respecto de x en cada


    punto de la gr´afica de fd(x, y).

    Recordar que queremos obtener la funci´on y=f(x), se puede predecir cada valor yi+1 a partir del valor
    presente yi
    y de la pendiente ϕ (ver figura).
    La ecuaci´on para obtener la funci´on mediante el m´etodo de Euler es la siguiente donde ϕ = f (xi , yi ) :

    yi+1 = yi + f (xi, yi) h


    Notar que depende de h que es el taman˜o de paso. Entre menor sea el taman˜o del paso es mejor la extrapolaci
    ´on.

    1
    C´odigo de programaci´on ejemplo
    %Se quiere resolver la siguiente ecuaci´on diferencial:
    %dy/dx = -2x^3 +12x^2 -20x + 8.5 desde x=0 a x=4 con un paso de h=0.5.
    %la soluci´on o condici´on inicial es: y(x=0)=1
    clc;clear;close all;
    %M´etodo de euler
    dydx=@(x,y) -2*x.^3 +12*x.^2 -20*x + 8.5
    h=0.5;
    x=0:h:4;
    n=length(x);
    y=zeros(1,n)
    %cond inicial
    y(1)=1

    for i=1:n-1
    y(i+1)=y(i)+dydx(x(i),y(i))*h
    end
    %sol real
    x1=0:0.1:4
    f=@(x) -0.5*x.^4 +4*x.^3-10*x.^2+8.5*x+1

    hold on
    plot(x,y,'-*r');plot(x1,f(x1),'-b')
    legend('Sol Num´erica','Sol real')

    (a) h=0.5 (b) h=0.1

    Figura 1: Resultados empleando el m´etodo de Euler

    Se puede observar de acuerdo con la figura (a) que el taman˜o de paso 0.5 no es muy adecuado para estimar
    la funci´on, se requiere un paso m´as pequen˜o para capturar mejor las curvas, por ejemplo si se emplea un paso
    de 0.1 la soluci´on mejora.
    Por supuesto hay mejores formas de obtener la soluci´on. Vayamos a explorarlas!

    2
    2. M´etodo de Heum
    El m´etodo de Heum emplea 2 derivadas promediadas para mejorar la estimaci´on de la pendiente. La
    formulaci´on del m´etodo de Heum es la siguiente:
    ∗
    yi+1 = yi + f (xi , yi ) h
    ∗
    f (xi , yi ) + f xi+1 , yi +1
    yi+1 = + h
    yi

    y lo2usamos en yi+1 para la aproximaci´on de la funci


    ∗
    Es decir estimamos primero por el m´etodo de Euler yi +1
    ´on.

    C´odigo de programaci´on - M´etodo de Heum


    clc;clear;close all;
    %M´etodo de euler
    dydx=@(x,y) -2*x.^3 +12*x.^2 -20*x + 8.5
    h=0.5;
    x=0:h:4;
    n=length(x);
    y0=zeros(1,n)
    y=zeros(1,n);
    y(1)=1;
    %cond inicial

    for i=1:n-1
    yaux(i+1)=y(i)+dydx(x(i),y(i))*h y(i+1)=y(i)+
    (dydx(x(i),y(i))+ dydx(x(i+1),yaux(i+1)))*(h/2)
    end
    %sol real
    x1=0:0.1:4
    f=@(x) -0.5*x.^4 +4*x.^3-10*x.^2+8.5*x+1

    hold on
    plot(x,y,'-*r');plot(x1,f(x1),'-b')
    legend('Sol Num´erica','Sol real')

    Figura 2: M´etodo de Heum para un paso de h=0.5

    3
    3. M´etodo del punto medio
    Los m´etodos siguientes se resuelven de la misma manera, lo que va a cambiar es la forma en que debemos
    iterar para llegar a la soluci´on. La formula del m´etodo del punto medio es:

    yi+1 = yi + k2h

    k1 = f (xi, yi)
    1 1
    k =f
    2 xi + h, yi +
    2 2
    4. M´etodo de
    Ralston

    y =y 1 2
    i
    + h
    i k1 +
    3 3
    k1 = f (xi, yi)
    3 3
    k =f
    2 xi + h, yi +
    4 4
    5. M´etodos de Runge-Kutta
    5.1. Runge-Kutta de tercer orden (RK3)
    1
    yi+1 = yi + (k1 + 4k2 + k3) h
    6
    k1 = f (xi, yi)
    1 1
    k =f
    2 xi + h, yi +
    2 2
    k3 = f (xi + h, yi − k1h + 2k2h)

    5.2. Runge-Kutta de cuarto orden (RK4)


    1
    yi+1 = yi + (k1 + 2k2 + 2k3 + k4) h
    6
    k1 = f (xi, yi)
    1 1
    k =f
    2 xi + h, yi +
    12 12
    k =f
    3 xi + h, yi +
    2 2
    k4 = f (xi + h, yi + k3h)

    4
    5.3. Runge-Kutta de quinto orden (RK5)
    1
    yi+1 = yi + (7k1 + 32k3 + 12k4 + 32k5 + 7k6) h
    90
    k1 = f (xi, yi)
    1 1
    k =f
    2 xi + h, yi +
    14 14 1
    k =f
    3 xi + h, yi + k1h +
    4 8 8
    1 1
    k 4= f − k 2h + k 3h
    xi + 2
    2
    3 3 9
    k 5= f k h1 + k h4
    xi + h, yi + 1
    4
    3 12 12 8
    k 6= f x k h−
    3 k4h +
    2 7
    + h, y − 7 7
    i i k1h +
    7 7

    También podría gustarte