Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 51 vistas

    Runge Kutta de Cuarto Orden

    Cargado por

    Jorge Perez
    Este documento presenta la solución numérica de un PVI utilizando el método de Runge-Kutta de orden 4. Se calcula la solución en el intervalo [1,2] con N=10 y se compara con la solución exacta, obteniendo un acuerdo exacto entre ambas soluciones.

    Copyright:

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

    Runge Kutta de Cuarto Orden

    Cargado por

    Jorge Perez
    51 vistas4 páginas
    Este documento presenta la solución numérica de un PVI utilizando el método de Runge-Kutta de orden 4. Se calcula la solución en el intervalo [1,2] con N=10 y se compara con la solución exacta, obteniendo un acuerdo exacto entre ambas soluciones.

    Título original

    Runge kutta de cuarto orden

    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 presenta la solución numérica de un PVI utilizando el método de Runge-Kutta de orden 4. Se calcula la solución en el intervalo [1,2] con N=10 y se compara con la solución exacta, obteniendo un acuerdo exacto entre ambas soluciones.

    Copyright:

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

    Runge Kutta de Cuarto Orden

    Cargado por

    Jorge Perez
    Este documento presenta la solución numérica de un PVI utilizando el método de Runge-Kutta de orden 4. Se calcula la solución en el intervalo [1,2] con N=10 y se compara con la solución exacta, obteniendo un acuerdo exacto entre ambas soluciones.

    Copyright:

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

    𝑑𝑦

    Calcule la solución numérica para el PVI: 𝑑𝑥 = 𝑦𝑥 2 − 1,1𝑦, con 𝑦(1)=1, en 𝐼=[1,2], con el
    método de Runge-Kutta 4, con 𝑁=10 y redondeo a 4 cifras decimales, compare el resultado con
    la solución exacta:
    1 3 2,3
    𝑦(𝑡) = e3𝑡 −1,1𝑡+
    3

    Solución por el método de Runge kutta 4 orden

    ecuaciones a utilizar para la implementación del método:

    ℎ
    𝑦(𝑡𝑖+1 ) = 𝑦(𝑡0 + (𝑖 + 1)ℎ) ≈ 𝑦𝑖+1 = 𝑦𝑖 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4)
    6
    𝑘1 = 𝑓(𝑡𝑖 , 𝑦𝑖 )
    1 1
    𝑘2 = 𝑓 (𝑡𝑖 + ℎ, 𝑦𝑖 + ℎ𝑘1)
    2 2
    1 1
    𝑘3 = 𝑓 (𝑡𝑖 + ℎ, 𝑦𝑖 + ℎ𝑘2)
    2 2
    𝑘4 = 𝑓(𝑡𝑖 + ℎ, 𝑦𝑖 + ℎ𝑘3)
    𝑡𝑁 − 𝑡0
    ℎ=
    𝑁
    inicialmente tenemos:
    𝑡0 = 1, 𝑦0 = 1, 𝑡𝑛 = 2
    Pasando la ecuación diferencial inicial a términos de t nos queda:
    𝑓(𝑡𝑖 , 𝑦𝑖 ) = 𝑦𝑖 𝑡𝑖2 − 1,1𝑦𝑖

    Calculamos el valor de h:
    𝑡𝑁 − 𝑡0
    ℎ=
    𝑁
    2−1
    ℎ= = 0,1
    10
    𝑡𝑖 = 𝑡0 + 𝑖ℎ, 𝑖 = 0,1,2,3,4 … 9
    i=0
    𝑘1 = 𝑓(𝑡0 , 𝑦0 ) = 𝑓(1,1) = −0,1
    1 1 1 1
    𝑘2 = 𝑓 (𝑡0 + ℎ, 𝑦0 + ℎ𝑘1) = 𝑓 (1 + 0.1 , 1 + 0.1𝑘1) = 𝑓(1.0500 , 0.9950) = 0.0025
    2 2 2 2
    1 1 1 1
    𝑘3 = 𝑓 (𝑡0 + ℎ, 𝑦0 + ℎ𝑘2) = 𝑓 (1 + 0.1 , 1 + 0.1𝑘2) = 𝑓(1.0500 , 1.0001) = 0.0025
    2 2 2 2
    𝑘4 = 𝑓(𝑡0 + ℎ, 𝑦0 + ℎ𝑘3) = 𝑓(1 + 0.1 , 1 + 0.1𝑘3) = 𝑓(1.1000 , 1.0003) = 0.1100

    ℎ
    𝑦(𝑡1) = 𝑦(𝑡0 + 1ℎ) ≈ 𝑦1 = 𝑦0 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4)
    6
    0.1
    𝑦1 = 1 + (−0.1 + 2 ∗ 0.0025 + 2 ∗ 0.0025 + 0.1100)
    6
    𝑦1 = 1,0003

    i=1
    𝑘1 = 𝑓(𝑡1 , 𝑦1 ) = 𝑓(1.1 , 1.0003) = 0.1100
    1 1 1 1
    𝑘2 = 𝑓 (𝑡1 + ℎ, 𝑦1 + ℎ𝑘1) = 𝑓 (1,1 + 0.1 , 1.0003 + 0.1𝑘1) = 𝑓(1.1500 , 1.0058) = 0.2238
    2 2 2 2
    1 1 1 1
    𝑘3 = 𝑓 (𝑡1 + ℎ, 𝑦1 + ℎ𝑘1) = 𝑓 (1,1 + 0.1 , 1.0003 + 0.1𝑘2) = 𝑓(1.1500 , 1.0115) = 0.2251
    2 2 2 2
    𝑘4 = 𝑓(𝑡1 + ℎ, 𝑦1 + ℎ𝑘3) = 𝑓(1.1 + 0.1 , 1.0003 + 0.1𝑘3) = 𝑓(1.2000 , 1.0228) = 0.3478

    ℎ
    𝑦(𝑡2) = 𝑦(𝑡1 + 1ℎ) ≈ 𝑦2 = 𝑦1 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4)
    6
    0.1
    𝑦2 = 1.0003 + (0.1100 + 2 ∗ 0.2238 + 2 ∗ 0.2251 + 0.3478)
    6
    𝑦2 = 1,0229

    Seguimos haciendo las iteraciones siguiendo los mismos pasos hasta llegar a i=9 donde
    obtendremos el valor de 𝑦10 , los datos de las demás iteraciones se resumen en la siguiente tabla:
    i t y k1 k2 k3 k4
    -
    0 1,0 1,0000 0,1000 0,0025 0,0025 0,1100
    1 1,1 1,0003 0,1100 0,2238 0,2251 0,3478
    2 1,2 1,0229 0,3478 0,4811 0,4842 0,6321
    3 1,3 1,0714 0,6321 0,7969 0,8029 0,9905
    4 1,4 1,1518 0,9906 1,2043 1,2151 1,4643
    5 1,5 1,2734 1,4644 1,7539 1,7728 2,1179
    6 1,6 1,4506 2,1179 2,5255 2,5585 3,0546
    7 1,7 1,7063 3,0543 3,6483 3,7066 4,4447
    8 1,8 2,0765 4,4436 5,3386 5,4425 6,5780
    9 1,9 2,6195 6,5750 7,9677 8,1559 9,9618
    10 2,0 3,4326 9,9545 12,1938 12,5412 15,5130

    Los datos correspondientes a la columna de las y son la solución numérica de la ecuación


    diferencial.

    A continuación, se presenta la comparación entre la solución numérica y la solución exacta, para


    obtener los datos de la solución exacta evaluamos los valores de x en la siguiente función:

    1 3 2,3
    𝑦(𝑡) = e3𝑡 −1,1𝑡+
    3

    Obteniendo los siguientes datos.

    x 1 1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1,9 2


    Solución numérica 1 1 1,023 1,071 1,152 1,273 1,451 1,706 2,076 2,620 3,433
    Solución exacta 1 1 1,023 1,071 1,152 1,273 1,451 1,706 2,076 2,620 3,433
    solucion numerica y solucion exacta
    4

    3,5

    2,5

    1,5

    0,5

    0
    0 0,5 1 1,5 2 2,5

    Numerica Exacta

    También podría gustarte