Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 85 vistas

    Euler - Rodrigo

    Cargado por

    selvino azorza
    El documento describe dos métodos para resolver ecuaciones diferenciales de primer orden: el método de Euler y el método de Heun. El método de Euler encuentra iterativamente la solución avanzando con un paso h. El método de Heun es una mejora que toma el promedio entre el valor de Euler y una aproximación en el punto siguiente, obteniendo resultados más precisos que Euler con la misma cantidad de puntos.

    Copyright:

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

    Euler - Rodrigo

    Cargado por

    selvino azorza
    85 vistas4 páginas
    El documento describe dos métodos para resolver ecuaciones diferenciales de primer orden: el método de Euler y el método de Heun. El método de Euler encuentra iterativamente la solución avanzando con un paso h. El método de Heun es una mejora que toma el promedio entre el valor de Euler y una aproximación en el punto siguiente, obteniendo resultados más precisos que Euler con la misma cantidad de puntos.

    Descripción original:

    euler

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento describe dos métodos para resolver ecuaciones diferenciales de primer orden: el método de Euler y el método de Heun. El método de Euler encuentra iterativamente la solución avanzando con un paso h. El método de Heun es una mejora que toma el promedio entre el valor de Euler y una aproximación en el punto siguiente, obteniendo resultados más precisos que Euler con la misma cantidad de puntos.

    Copyright:

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

    Euler - Rodrigo

    Cargado por

    selvino azorza
    El documento describe dos métodos para resolver ecuaciones diferenciales de primer orden: el método de Euler y el método de Heun. El método de Euler encuentra iterativamente la solución avanzando con un paso h. El método de Heun es una mejora que toma el promedio entre el valor de Euler y una aproximación en el punto siguiente, obteniendo resultados más precisos que Euler con la misma cantidad de puntos.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 4

    Metodo de Euler

    El metodo de Euler consiste en encontrar iterativamente la solucion de una ecuacion diferencial de


    primer orden y valores iniciales conocidos para un rango de valores. Partiendo de un valor inicial x0
    y avanzando con un paso h, se pueden obtener los valores de la solucion de la siguiente manera:

    Yk+1 = Yk + h f (xk , Yk )

    Donde Y es solucion de la ecuacion diferencial y f es la ecuacion diferencial en funcion de las


    variables independientes.

    Ejemplo

    Se quiere obtener el valor de la corriente en el siguiente circuito RL hasta un segundo con un paso
    de un cuarto de segundo:

    R = 1

    +
    v = sen(2t) L = 1Hy

    Planteando la malla correspondiente se puede obtener la ecuacion del circuito de la figura:

    di(t)
    v = R i(t) + L
    dt
    sen(t) = i(t) + i0 (t)
    i0 (t) = sen(2t) i(t) = f (t, i)

    Reescribiendo la formula de iteracion con nuestras variables del problema llegamos a:

    Ik+1 = Ik + h f (tk , Ik )
    En primer lugar, reconocemos nuestros datos iniciales como t0 = 0 e I0 = 0. Ademas, la funcion
    sera, en este caso, la derivada de la corriente. Respecto a tk , dado que hay un paso constante h, se
    puede observar que en general tk = t0 +hk; por lo tanto, se puede observar que tomara tres iteraciones
    llegar a t = 1 debido a que los valores de salida siempre estan un paso adelantado. Segun nuestros
    ndices, la primer iteracion sera la numero 0:

    k=0

    I1 = I0 + h f (0; 0) = 0

    k=1

    I2 = I1 + h f (0, 25; 0) = 0, 25

    k=2

    1
    I3 = I2 + h f (0, 5; 1) = 0, 1875

    k=3

    I4 = I3 + h f (0, 75; 1) = 0, 1093

    La solucion de la ecuacion diferencial, teoricamente, es:


    1
    i(t) = (2et + sen(2t) 2 cos(2t))
    1 + 4 2
    Comparando con los valores teoricos:

    t Valor teorico Valor aproximado Error


    0 0 0 0
    0,25 0,1455 0 0,1455
    0,5 0,2493 0,25 0,0007
    0,75 0,0486 0,1875 0,1389
    1 -0,0981 -0,1093 0,0112

    Para mejorar la aproximacion, se puede aumentar el tamano de puntos (es decir, reducir el tamano
    del paso h). Por otro lado, se puede utilizar un metodo de mayor orden para obtener una mejor apro-
    ximacion usando la misma cantidad de puntos.

    Metodo de Heun

    Este metodo consiste en una mejora del metodo 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 aproximacion del valor de la funcion en el punto
    siguiente, tambien por Euler.

    hh )
    i
    Yk+1 = Yk + f (xk , Yk ) + f (xk+1 , Yk+1
    2
    Donde Y es solucion de la ecuacion diferencial, f es la ecuacion diferencial en funcion de las va-
    es una aproximacion de Euler.
    riables independientes y la solucion de Yk+1

    Ejemplo

    Recordando la ecuacion diferencial del ejercicio anterior:

    i0 (t) = sen(2t) i(t) = f (t, i)


    Obtenemos la formula de iteracion segun el metodo de Heun:

    hh )
    i
    Ik+1 = Ik + f (tk , Ik ) + f (tk+1 , Ik+1
    2
    que
    Nuevamente, los valores iniciales son t0 = 0, I0 = 0. En este caso, definiremos un valor Ik+1
    corresponde a la aproximacion de Euler en el paso siguiente. Realizando los calculos:

    2
    k=0

    I1 = I0 + h f (0; 0) = 0
    hh i
    I1 = I0 + f (0; 0) + f (0, 25; I1 ) = 0, 125 = i(0, 25)
    2
    k=1

    I2 = I1 + h f (0, 25; 0, 125) = 0, 3437


    hh i
    I2 = I1 + f (0, 25; 0, 125) + f (0, 5; I2 ) = 0, 1914 = i(0, 5)
    2
    k=2

    I3 = I2 + h f (0, 5; 0, 1914) = 0, 1435


    hh i
    I3 = I2 + f (0, 5; 0, 1914) + f (0, 75; I3 ) = 0, 0245 = i(0, 75)
    2
    k=3

    I4 = I3 + h f (0, 75; 0, 0245) = 0, 2316


    hh i
    I4 = I3 + f (0, 75; 0, 245) + f (0, 75; I4 ) = 0, 0745 = i(1)
    2
    Comparando con los valores teoricos:

    t Valor teorico Valor aproximado Error


    0 0 0 0
    0,25 0,1455 0,125 0,00205
    0,5 0,2493 0,1914 0,0579
    0,75 0,0486 0,0245 0,00241
    1 -0,0981 -0,0745 0,00242

    A continuacion, se ilustran la solucion en un grafico que relaciona la curva esperada con los puntos
    obtenidos por el metodo de Heun y de Euler:

    3
    0.3
    Solucion de la ecuacion diferencial
    Heun
    Euler
    Solucion
    0.2

    0.1
    i(t)

    0.0

    0.1

    0.0 0.2 0.4 0.6 0.8 1.0


    t

    Figura 1: Graficos de aproximaciones y valores teoricos

    Se puede observar que, si bien para ciertos puntos el comportamiento es mejor para Euler, para el
    metodo de Heun, los puntos se comportan mejor segun la forma de la solucion respecto a la solucion
    de Euler.

    También podría gustarte