Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 334 vistas

    6.3.-Metodo de Runge Kutta

    Cargado por

    Julio César Millán Díaz
    El documento describe el método de Runge-Kutta para resolver ecuaciones diferenciales ordinarias. Explica que existen varias versiones del método dependiendo del número de términos utilizados y presenta tres métodos de segundo orden específicos: el método de Heun, el método del punto medio y el método de Ralston.

    Copyright:

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

    6.3.-Metodo de Runge Kutta

    Cargado por

    Julio César Millán Díaz
    334 vistas11 páginas
    El documento describe el método de Runge-Kutta para resolver ecuaciones diferenciales ordinarias. Explica que existen varias versiones del método dependiendo del número de términos utilizados y presenta tres métodos de segundo orden específicos: el método de Heun, el método del punto medio y el método de Ralston.

    Descripción original:

    Apuntes de la Academia de Ciencias Básicas d ela ESIME U. P. Ticomán

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

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

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento describe el método de Runge-Kutta para resolver ecuaciones diferenciales ordinarias. Explica que existen varias versiones del método dependiendo del número de términos utilizados y presenta tres métodos de segundo orden específicos: el método de Heun, el método del punto medio y el método de Ralston.

    Copyright:

    © All Rights Reserved
    Descargue como PPTX, PDF, TXT o lea en línea desde Scribd
    Descargar como pptx, pdf o txt
    334 vistas11 páginas

    6.3.-Metodo de Runge Kutta

    Cargado por

    Julio César Millán Díaz
    El documento describe el método de Runge-Kutta para resolver ecuaciones diferenciales ordinarias. Explica que existen varias versiones del método dependiendo del número de términos utilizados y presenta tres métodos de segundo orden específicos: el método de Heun, el método del punto medio y el método de Ralston.

    Copyright:

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

    INSTITUTO POLITECNICO NACIONAL

    Escuela Superior de Ingeniera


    Mecnica y Elctrica.

    Unidad Profesional Ticomn.


    Academia de Ciencias Bsicas.
    UNIDAD VI: Ecuaciones Diferenciales Ordinarias.

    La Tcnica al Servicio de la Patria

    Ing. Julio Csar Milln Daz.

    6.3.- Mtodo de Runge Kutta.

    Existen muchas variantes del mtodo de Euler, pero todas tienen la forma
    generalizada de la ecuacin:

    y i 1 y i ( xi , y i , h)h
    Donde ( xi , y i , h) se conoce como funcin incremento, la cual puede
    interpretarse como una pendiente representativa en el intervalo. La
    funcin incremento se escribe en forma general como:

    a1k1 a2 k 2 ... an k n

    Donde las a son constantes y las k son:

    k1 f ( x i , y i )
    k 2 f ( xi p1 h, y i q11 k1 h)
    k3 f ( xi p2 h, yi q21k1h q22 k 2 h)

    k n f ( xi p n 1 h, y i q n 1,1 k1 h q n 1, 2 k 2 h ... q n 1,n 1 k n 1 h)

    Donde las p y las q son constantes. Observe que las k son relaciones de
    recurrencia. Es decir, k1 aparece en la ecuacin k2, la cual aparece en la
    ecuacin k3, etctera. Como cada k es una evaluacin final, esta
    recurrencia vuelve eficientes a los mtodos RK parra clculos en
    computadora.
    Es posible tener varios tipos de mtodos de Runge Kutta empleando
    diferentes nmeros de trminos en la funcin incremento especificada por
    n. Observe que el mtodo de Runge Kutta (RK) de primer orden con n =
    1 es, de hecho, el mtodo de Euler.
    Una vez que se elige n, se evalan las a, p y q igualando yi 1 yi ( xi , yi , h)h a
    los trminos en la expansin de Taylor. As, al menos para las versiones
    de orden inferior, el numero de trminos, n, por lo comn representa el
    orden de la aproximacin.

    k1 = h f (xi , yi )

    k4 = h f ( xi + h, yi + k3 )

    h
    k1
    k2 = h f xi + , yi +

    2
    2

    h
    k2
    k3 = h f xi + , yi +

    2
    2

    1
    k prom = ( k1 + 2k2 + 2k3 + k4 )
    6
    yi+1 = yi + k prom

    6.3.1.- Mtodo de Runge Kutta de segundo orden.


    La versin de segundo orden de la ecuacin y i 1 y i ( xi , y i , h)h , es:

    y i 1 y i (a1 k1 a 2 k 2 )h
    Donde:

    k1 f ( x i , y i )
    k 2 f ( xi p1 h, y i q11 k1 h)

    Los valores de a1, a2, p1 y q11 se evalan al igualar yi 1 yi (a1k1 a2 k 2 )h con


    la expansin de Taylor hasta el trmino de segundo orden. Al hacerlo,
    desarrollamos tres ecuaciones para evaluar las cuatro constantes
    desconocidas. Las tres ecuaciones son:

    a1 a 2 1
    1
    a 2 p1
    2
    1
    a 2 q11
    2

    Como tenemos tres trminos con cuatro incgnitas, debemos dar el valor
    de una de estas incgnitas para determinar las otras tres. Suponga que
    damos un valor para a2. Entonces se resuelven de manera simultnea las
    ecuaciones para a1 a 2 1 y a2 q11 1 / 2 , obteniendo:

    a1 1 a 2

    1
    p1 q11
    2a 2

    Debido a que podemos elegir un nmero de valores para a2, hay un


    nmero infinito de mtodos de RK de segundo orden. Cada versin dara
    exactamente los mismos resultados si la solucin de la EDO fuera
    cuadrtica, lineal o una constante. Sin embargo se obtienen diferentes
    resultados cuando (como tpicamente es el caso) la solucin es ms
    complicada.

    6.3.1.1.-Mtodo de Heun con un solo corrector (a = 1/2).


    Si suponemos que a2 es de las ecuaciones a1 1 a 2 y p1 q11 1 2a2
    puede obtenerse a1 = y p1 = q11 = 1. Estos parmetros, al sustituirse
    en la ecuacin yi 1 yi (a1k1 a 2 k 2 )h , dan como resultado:

    1
    1
    y i 1 y i k1 k 2 h
    2
    2
    Donde:

    k1 f ( x i , y i )
    k 2 f ( xi h, y i k1 h)
    Observe que k1 es la pendiente al inicio del intervalo y que k2 es la
    pendiente al final del intervalo. En consecuencia, este mtodo de Runge
    Kutta de segundo orden es, de hecho, la tcnica de Heun sin iteracin.

    6.3.1.2.-El mtodo del punto medio (a2 = 1).


    Si suponemos que a2 es 1, entonces a1 = 0, p1 = q11 = , y la ecuacin
    , se convierte en:
    y i 1 y i (a1 k1 a 2 k 2 )h

    yi 1 yi k 2 h
    Donde:

    k1 f ( x i , y i )
    1
    1

    k 2 f x i h, y i k 1 h
    2
    2

    Este es el mtodo del punto medio.

    6.3.1.3.-Mtodo de Ralston (a = 2/3).


    Ralston (1962) y Ralston & Rabinowitz (1978) determinaron que al
    seleccionar a2 = 2/3 se obtiene un mnimo en el error de truncamiento
    para los algoritmos RK de segundo orden. Con esta versin, a1 = 1/3 y p1
    = q11 = y se obtiene:

    y i 1

    2
    1
    y i k1 k 2 h
    3
    3

    Donde:

    k1 f ( x i , y i )

    3
    3

    k 2 f x i h, y i k 1 h
    4
    4

    También podría gustarte