Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 74 vistas

    Euler Modificado

    Cargado por

    Adolfo Velásquez
    Este documento describe el método de Euler modificado (también conocido como método de Heun) para resolver ecuaciones diferenciales numéricamente. El método de Euler modificado mejora la precisión del método de Euler al aproximar el área bajo la curva como el promedio de las pendientes en los puntos extremos, en lugar de solo usar la pendiente en el punto inicial. A continuación, se presentan las fórmulas de recurrencia para este método y ejemplos de su aplicación para resolver sistemas de ecuaciones diferenciales.

    Copyright:

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

    Euler Modificado

    Cargado por

    Adolfo Velásquez
    74 vistas18 páginas
    Este documento describe el método de Euler modificado (también conocido como método de Heun) para resolver ecuaciones diferenciales numéricamente. El método de Euler modificado mejora la precisión del método de Euler al aproximar el área bajo la curva como el promedio de las pendientes en los puntos extremos, en lugar de solo usar la pendiente en el punto inicial. A continuación, se presentan las fórmulas de recurrencia para este método y ejemplos de su aplicación para resolver sistemas de ecuaciones diferenciales.

    Título original

    Euler_modificado

    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?

    Este documento describe el método de Euler modificado (también conocido como método de Heun) para resolver ecuaciones diferenciales numéricamente. El método de Euler modificado mejora la precisión del método de Euler al aproximar el área bajo la curva como el promedio de las pendientes en los puntos extremos, en lugar de solo usar la pendiente en el punto inicial. A continuación, se presentan las fórmulas de recurrencia para este método y ejemplos de su aplicación para resolver sistemas de ecuaciones diferenciales.

    Copyright:

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

    Euler Modificado

    Cargado por

    Adolfo Velásquez
    Este documento describe el método de Euler modificado (también conocido como método de Heun) para resolver ecuaciones diferenciales numéricamente. El método de Euler modificado mejora la precisión del método de Euler al aproximar el área bajo la curva como el promedio de las pendientes en los puntos extremos, en lugar de solo usar la pendiente en el punto inicial. A continuación, se presentan las fórmulas de recurrencia para este método y ejemplos de su aplicación para resolver sistemas de ecuaciones diferenciales.

    Copyright:

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

    Método de Euler Modificado

    Dr. Mario González Cardel


    Euler Modificado Análisis Numérico

    Como su nombre o indica este método es una


    modificación al método de Euler para mejorar su precisión
    y se conoce también como Método de Heun.
    Euler Modificado Análisis Numérico

    Cuando aplicamos el método de Euler,


    xn+1 = xn + h
    yn+1 = f ( xn , yn )h + yn
    el producto:
    f (xn , yn )h
    Representa una aproximación al área bajo la curva:
    Euler Modificado Análisis Numérico

    Aproximando la pendiente por su valor en el punto extremo


    izquierdo.
    Se obtiene una mejor aproximación si tomamos el promedio
    de sus valores en los puntos extremos, entonces el área bajo
    la curva es aproximada por el área del rectángulo
    sombreado.

    xn+1 = xn + h
    yn+1 = yn +  f ( xn , yn ) + f ( xn+1, yn+1 )h
    Euler Modificado Análisis Numérico

    La expresión de recurrencia en la diapositiva anterior para


    calcular yn+1 depende de yn+1, pero podemos aproximar la
    yn+1 por la de Euler, es decir que la expresión de
    recurrencia queda:
    xn+1 = xn + h

    yn+1 = yn
    H n
    +
    ) +
    f (x , y
    f ( xH
    n+1
    n , y E
    n+1 )h
    2
    El superíndice H significa aproximación de Heun, y el
    superíndice E significa una aproximación de Euler, es decir:
    xn+1 = xn + h
    ynE+1 = yn + f ( xn , ynH )h

    n+1 
    y = yn + f ( xn , y ) + f ( xn+1, y )
    H H
    n
    E h
    n+1
    2

    Ésta es la fórmula de Euler modificado o fórmula de Heun.
    Euler Modificado Análisis Numérico

    Ejemplo.
    Resolver la siguiente ecuación diferencial, con condiciones
    iníciales, usando el método de Taylor.
    con la condición inicial y(0) = 1

    Solución:
    Sabemos que:
    x0 = 0, es el centro de convergencia.
    y(0) = 1
    Evaluando en la ecuación diferencial, tenemos:
    Euler Modificado Análisis Numérico

    Solución:
    Sabemos que:
    x0 = 0, es el centro de convergencia.
    y(0) = 1
    La iteración número 1 es:
    x1 = x0 + h = 0 + 0.1 = 0.1
    y1E = y0 + f ( x0 , y0H )h = 1 + (1 − 0 + 4(1))(0.1) = 1.5
     
    E h
    y1 = y0 + f ( x0 , y0 ) + f ( x1, y1 )
    H H
    2
    = 1 + (1 − 0 + 4(1)) + (1 − 0.1 + 4(1.5)) = 1.5950
    0.1
    2
    Euler Modificado Análisis Numérico

    La iteración número 2 es:

    x2 = x1 + h = 0.1 + 0.1 = 0.2


    y2E = y1H + f ( x1, y1H )h = 1.595+ (1 − 0.1 + 4(1.595))(0.1) = 2.323
    H H
     H

    E h
    y2 = y1 + f ( x1, y1 ) + f ( x2 , y2 )
    2
    = 1.595+ (1 − 0.1 + 4(1.595)) + (1 − 0.2 + 4(2.323)) = 2.4636
    0.1
    2
    Euler Modificado Análisis Numérico

    La iteración
    . número 3 es:

    x3 = x2 + h = 0.2 + 0.1 = 0.3


    y3E = y2H + f ( x2 , y2H )h = 2.4636+ (1 − 0.2 + 4(2.4636))(0.1) = 3.52904
    H H
     H
    E h
    y3 = y2 + f ( x2 , y2 ) + f ( x3 , y3 )
    2
    = 2.4636+ (1 − 0.2 + 4(2.4636)) + (1 − 0.3 + 4(3.52904)) = 3.7371
    0.1
    2
    Euler Modificado Análisis Numérico

    Las primeras 10 iteraciones se muestran en la siguiente tabla:

    n xn f(xn, ynh) y nE f(xn, yn+1h) y nH


    0 0 5 1 1
    1 0.1 7.28 1.5 6.9 1.595
    2 0.2 10.6544 2.323 10.092 2.4636
    3 0.3 15.64851 3.52904 14.81616 3.73712
    4 0.4 23.03979 5.30197 21.80791 5.60994
    5 0.5 33.9789 7.91392 32.15571 8.36972
    6 0.6 50.16877 11.76761 47.47046 12.44219
    7 0.7 74.12978 17.45907 70.13628 18.45744
    8 0.8 109.59208 25.87042 103.68169 27.34802
    9 0.9 162.07627 38.30722 153.32891 40.49406
    10 1 56.70169 226.80679 59.93822
    Euler Modificado Análisis Numérico

    Comparación de la solución analítica y la


    numérica
    80000
    y(x)
    f(x)

    60000

    40000

    20000

    0
    0 0.5 1 1.5 2 2.5 3
    Euler Modificado Análisis Numérico

    Ejemplo.
    Determinar la solución del sistema de ecuaciones
    diferenciales por el método de Euler modificado.
    x' (t ) = y + cos(2t )
    y' (t ) = −x
    Euler Modificado Análisis Numérico
    ,

    Solución:
    Sabemos que: x(0) = 1 y(0) = 4
    Las expresiones de recurrencia se modifican, quedando:
    tn+1 = tn + h
    x E
    n+1 = x + f ( x , y )h
    H
    n
    H
    n
    H
    n

    H
    n+1
    H
    n 
    x = x + f (x , y ) + f (x , y )
    H
    n
    H
    n
    E
    n+1
    E
    n+1 
    h
    2
    ynE+1 = ynH + g ( xnH , ynH )h
    H

    yn+1 = yn + g ( xn , yn ) + g ( xn+1, yn+1 )
    H H H H E

    h
    2
    Euler Modificado Análisis Numérico

    La iteración número 1 es:


    t1 = t0 + h = 0 + 0.1 = 0.1
    x1E = x0H + f ( x0H , ynH )h = 1 + (4 + cos(2 * (0)))(0.1) = 1.5
    H
     H H E

    x1 = x0 + f ( x0 , y0 ) + f ( x1 , y1 )
    H E h
    2
    = 1 + (4 + cos(2 * (0))) + (4 + cos(2(0.1)) = 1.49900
    0.1
    2
    y1E = y0H + g ( x1H , y0H )h = 4 + (−1.5)(0.1) = 3.85
    H
     E H H

    y1 = y0 + g ( x1 , y0 ) + g ( x1 , y1 )
    H E h
    2
    = 4 + (−1.500) + (−1.4900) = 3.85045
    0.1
    2
    Euler Modificado Análisis Numérico

    La iteración número 2 es:


    t2 = t1 + h = 0 + 0.2 = 0.3
    x2E = x1H + f ( x1H , y1H )h = 1.49900+ (3.85005+ cos(2 * (0.1)))(0.1) = 1.98201
    H
     H H

    x2 = x1 + f ( x1 , y1 ) + f ( x2 , y2 )
    H E E h
    2
    = 1 + (1.49900+ cos(2 * (0.1))) + (1.98201+ cos(2(0.2)) = 1.96915
    0.1
    2
    y2E = y1H + g ( x2H , y1H )h = 3.85045+ (−1.96915)(0.1) = 3.65180
     
    y2H = y1H + g ( x2E , y1H ) + g ( x2H , y2E )
    h
    2
    = 3.85045+ (−1.49900) + (−1.96915) = 3.65249
    0.1
    2
    Euler Modificado Análisis Numérico

    Las primeras 10 iteraciones se muestran en la siguiente tabla:


    n t xE xH yE yH
    . 0 0 1 1 4 4
    1 0.1 1.5 1.499 3.85 3.85004
    2 0.2 1.98201 1.96915 3.65179 3.65249
    3 0.3 2.4265 2.40955 3.40914 3.41068
    4 0.4 2.83315 2.81248 3.12583 3.1284
    5 0.5 3.19499 3.17106 2.80633 2.8101
    6 0.6 3.50611 3.47949 2.45572 2.46082
    7 0.7 3.76181 3.73312 2.07954 2.08607
    8 0.8 3.95873 3.92865 1.68366 1.6917
    9 0.9 4.0949 4.06412 1.27417 1.28375
    10 1 4.16978 4.139 0.85719 0.86831
    Euler Modificado Análisis Numérico

    Comparación de la solución analítica y la numérica.


    Euler Modificado Análisis Numérico

    Conclusión.
    • El método de Heun brinda mejor precisión que los métodos
    vistos anteriormente, es aplicable a ecuaciones
    diferenciales ordinarias y a sistemas de ecuaciones de
    ecuaciones acopladas.
    • Es un método fácil de programar.

    • No podemos estimar el error.

    También podría gustarte