Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 157 vistas

    Método de Heum para Ecuaciones Diferenciales

    Cargado por

    Paul Mariatte Blue
    Este documento presenta un informe de laboratorio sobre el método de Heun para resolver ecuaciones diferenciales ordinarias de valor inicial. Explica que el método de Heun mejora el método de Euler al aproximar la pendiente mediante el promedio de las derivadas en el punto inicial y final del intervalo. Luego, describe las ecuaciones y variables utilizadas por el algoritmo del método de Heun para obtener una aproximación numérica a la solución del problema de valor inicial planteado.

    Copyright:

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

    Método de Heum para Ecuaciones Diferenciales

    Cargado por

    Paul Mariatte Blue
    157 vistas4 páginas
    Este documento presenta un informe de laboratorio sobre el método de Heun para resolver ecuaciones diferenciales ordinarias de valor inicial. Explica que el método de Heun mejora el método de Euler al aproximar la pendiente mediante el promedio de las derivadas en el punto inicial y final del intervalo. Luego, describe las ecuaciones y variables utilizadas por el algoritmo del método de Heun para obtener una aproximación numérica a la solución del problema de valor inicial planteado.

    Descripción original:

    Método de Heum para ecuaciones diferenciales

    Título original

    Método de Heum para ecuaciones diferenciales

    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 un informe de laboratorio sobre el método de Heun para resolver ecuaciones diferenciales ordinarias de valor inicial. Explica que el método de Heun mejora el método de Euler al aproximar la pendiente mediante el promedio de las derivadas en el punto inicial y final del intervalo. Luego, describe las ecuaciones y variables utilizadas por el algoritmo del método de Heun para obtener una aproximación numérica a la solución del problema de valor inicial planteado.

    Copyright:

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

    Método de Heum para Ecuaciones Diferenciales

    Cargado por

    Paul Mariatte Blue
    Este documento presenta un informe de laboratorio sobre el método de Heun para resolver ecuaciones diferenciales ordinarias de valor inicial. Explica que el método de Heun mejora el método de Euler al aproximar la pendiente mediante el promedio de las derivadas en el punto inicial y final del intervalo. Luego, describe las ecuaciones y variables utilizadas por el algoritmo del método de Heun para obtener una aproximación numérica a la solución del problema de valor inicial planteado.

    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 4

    Informe de Laboratorio N01-02

    Del alumno

    Helbert Justo Luque Zevallos

    Al profesor

    Angel Sangiacomo

    Asunto

    Mtodo de Heun

    Fecha

    02-09-2015

    Tengo a bien informar a usted sobre la prctica en el laboratorio sobre el


    mtodo de Heun
    Qu es?
    Es un algoritmo que implementa el mtodo de Heun para la solucin de E.D.O.
    con valor inicial.
    Cmo es?
    La mejora del mtodo de Heun consiste en la aproximacin a la pendiente
    mediante la aplicacin de dos derivadas del intervalo, una en el punto inicial y
    otra en punto final. La aproximacin mejorada de la pendiente ser el promedio
    de las dos derivadas.
    Recordando el mtodo de Euler, la pendiente al principio del intervalo es
    y ' =f (xi , y i)

    que se emplea en extrapolar linealmente el VALOR de y al final del intervalo:


    y 0i+1= yi + f ( xi , y i )h

    Este valor no es la solucin final sino una prediccin intermedia, por lo que se
    ha distinguido a sta con el superndice 0. Esta ecuacin se denomina
    ecuacin predictora y proporciona una aproximacin del valor de y al final del
    intervalo. Este valor nos permite a su vez calcular la pendiente aproximada en
    dicho punto.
    y 'i+1= yi + f ( x i+1 , y 0i +1 )h

    Combinando las dos pendientes obtenemos el promedio del intervalo

    f ( xi , y i ) +f ( x i +1 , y i+1 )
    y +y
    y = i i +1 =
    2
    2
    '

    '

    '

    Esta pendiente promedio se usa para extrapolar linealmente el valor de la


    funcin en el siguiente punto usando el mtodo de Euler:

    f ( x i , yi ) + f ( x i+1 , y 0i+1 )
    y i+1= yi +
    h
    2

    Esta ecuacin se denomina ecuacin correctora

    Qu hace?
    Obtiene una aproximacin al problema de valor inicial bien planteado.
    at b

    y
    =f ( t , y )
    t
    y 0=a
    k 1=hf ( t i , y i )
    2
    2
    k 2=hf t i + h , y i+ k 1
    3
    3

    y i+1= yi +

    1
    [ k +3k 2 ]
    4 1

    Variables:
    fu: Funcin f(t,y) en formato texto debe contener variables t y y
    vi: Valor de Valor Inicial
    a: Valor de a
    b: Valor de b
    n: Numero de iteraciones
    h: Tamao de paso
    ti: Valor Inicial de la variable t
    k1: contiene la primera parte del mtodo de heun
    k2: contiene la segunda parte del mtodo de heun
    yii: Valores de la variable y obtenido por el mtodo de heun

    Algoritmo
    function f=f(z,t,y)
    f=eval(z);
    z : funcin en formato texto debe contener variables t y y
    mheun.m
    clc;
    disp(' METODO DE EULER ');
    fu=input('ingrese funcion f(t,y) = ');
    vi=input('ingrese Valor de Valor Inicial = ');
    a=input('ingrese Valor de A = ');
    b=input('ingrese Valor de B = ');
    n=input('ingrese iteraciones = ');
    h=(b-a)/n;
    disp('h=');
    disp(h);
    ti=a;
    yii=vi;
    %--------------------hold off;
    xlabel('Pasos')
    ylabel('Valor aproximado')
    %--------------------plot(ti,yii,'+');
    disp(' ti
    y(ti)
    ');
    for i=1:n+1
    fprintf ('\n%5.5f
    %5.5f
    ',ti,yii);
    yii=yii+f(fu,ti,yii)*h;
    ti=ti+h;
    %--------------------hold on;
    plot(ti,yii,'+');
    %--------------------end
    %--------------------hold on;
    legend('Metodo de euler');
    %--------------------disp('');

    METODO DE HEUM
    ingrese funcion f(t,y) = 'y+t'
    ingrese Valor de Valor Inicial = 0
    ingrese Valor de A = 0
    ingrese Valor de B = 1
    ingrese iteraciones = 10
    h=
    0.100000000000000
    ti
    y(ti)
    0.00000
    0.10000
    0.20000
    0.30000
    0.40000
    0.50000
    0.60000
    0.70000
    0.80000
    0.90000
    1.00000

    0.00000
    0.00500
    0.02103
    0.04923
    0.09090
    0.14745
    0.22043
    0.31157
    0.42279
    0.55618
    0.71408

    También podría gustarte