Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 2 vistas

    9999_merged

    Cargado por

    amelgarbaltazar2017
    asdfa

    Copyright:

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

    9999_merged

    Cargado por

    amelgarbaltazar2017
    2 vistas9 páginas
    asdfa

    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?

    asdfa

    Copyright:

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

    9999_merged

    Cargado por

    amelgarbaltazar2017
    asdfa

    Copyright:

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

    “Año del Bicentenario, de la consolidación de nuestra Independencia, y de la conmemoración de las

    Heroicas Batallas de Junín y Ayacucho”

    UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

    “Universidad del Perú. Decana de América”

    FACULTAD DE INGENIERÍA ELECTRÓNICA Y ELÉCTRICA

    Trabajo N° 9

    Asignatura:
    Métodos Numéricos

    Docente:
    Elmer Lévano

    Alumno:

    Melgar Baltazar Jesus Antonio - 23190065

    2024 - II
    LIMA - PERÚ
    Solución al problema 1
    Queremos resolver el siguiente problema de valor inicial (PVI):
    ( √
    y ′ (t) = (1 + 2t) y,
    y(0) = 1,

    usando el método de Runge-Kutta de orden 2 (RK2) para diferentes tamaños


    de paso h = {0.5, 0.25, 0.125} en el intervalo [0, 1].

    Método de Runge-Kutta de Orden 2


    La fórmula del método RK2 es:

    yn+1 = yn + h · Φ(tn , yn ),

    donde  
    1 h h
    Φ(tn , yn ) = f (tn , yn ) + hf tn + , yn + f (tn , yn ) .
    2 2 2
    En este caso, la función f (t, y) está dada por:
    √
    f (t, y) = (1 + 2t) y.

    Cálculos para h = 0.5


    1−0
    Dividimos el intervalo [0, 1] en 0.5 = 2 pasos.

    Paso 0 (t0 = 0, y0 = 1)
    √
    k1 = h · f (t0 , y0 ) = 0.5 · (1 + 2 · 0) 1 = 0.5 · 1 = 0.5,
    √ √
     
    h k1
    k2 = h · f t0 + , y0 + = 0.5 · (1 + 2 · 0.25) 1 + 0.25 = 0.5 · 1.5 · 1.25.
    2 2
    Calculamos:
    1
    y1 = y0 + (k1 + k2 ).
    2

    Paso 1 (t1 = 0.5, y1 )


    Repetimos el proceso usando el nuevo valor de y1 :

    k1 = h · f (t1 , y1 ),
     
    h k1
    k2 = h · f t1 + , y1 + ,
    2 2
    1
    y2 = y1 + (k1 + k2 ).
    2

    1
    Cálculos para h = 0.25
    1−0
    De manera similar, dividimos el intervalo en 0.25 = 4 pasos. Para cada paso,
    calculamos:
    k1 = h · f (tn , yn ),
     
    h k1
    k2 = h · f tn + , yn + ,
    2 2
    1
    yn+1 = yn + (k1 + k2 ).
    2

    Resultados
    t h = 0.5 h = 0.25 h = 0.125
    0.0 1.0000 1.0000 1.0000
    0.125 −−− −−− 1.1450
    0.25 −−− 1.3315 1.3355
    0.375 −−− −−− 1.5795
    0.5 1.8385 1.8758 1.8866
    0.625 −−− −−− 2.2677
    0.75 −−− 2.7135 2.7352
    0.875 −−− −−− 3.3029
    1.0 3.8215 3.9488 3.9862

    2
    Resolución del PVI usando el Método de Runge-
    Kutta de Orden 2
    Dado el problema:
    (
    y ′ (t) = y · sin3 (t),
    y(0) = 1,
    aplicamos el método de Runge-Kutta de orden 2 para aproximar los valores
    de y(t) en el intervalo [0, 1].

    Pasos del Método de Runge-Kutta de Orden 2


    El método se define por las siguientes fórmulas:

    k1 = h · f (tn , yn ),
     
    h k1
    k2 = h · f tn + , yn + ,
    2 2
    yn+1 = yn + k2 ,
    donde f (t, y) = y · sin3 (t). A continuación, calculamos los valores para cada
    paso.

    Cálculos para h = 0.5


    1. Paso inicial: t0 = 0, y0 = 1.
    k1 = 0.5 · (1 · sin3 (0)) = 0,
      
    0 3 0.5
    = 0.5 · 1 · sin3 (0.25) = 0.0154,
    
    k2 = 0.5 · 1 + · sin 0 +
    2 2
    y1 = 1 + 0.0154 = 1.0154.
    2. Paso siguiente: t1 = 0.5, y1 = 1.0154.
    k1 = 0.5 · (1.0154 · sin3 (0.5)) = 0.1227,
     
    0.1227
    k2 = 0.5 · 1.0154 + · sin3 (0.75) = 0.1735,
    2
    y2 = 1.0154 + 0.1735 = 1.1889.
    3. Paso siguiente: t2 = 1.0, y2 = 1.1889.
    k1 = 0.5 · (1.1889 · sin3 (1.0)) = 0.2493,
     
    0.2493
    k2 = 0.5 · 1.1889 + · sin3 (1.25) = 0.3678,
    2
    y3 = 1.1889 + 0.3678 = 1.5567.

    1
    [width=0.5]e6699441-fc30-48a0-b315-2479636e10e9.png

    Cálculos para h = 0.25


    Realizando cálculos similares, se obtienen los siguientes valores:

    t = [0.0, 0.25, 0.5, 0.75, 1.0],

    y = [1.0, 1.0039, 1.0311, 1.1193, 1.2781].

    Cálculos para h = 0.125


    Con un paso más pequeño, los resultados son:

    t = [0.0, 0.125, 0.25, 0.375, 0.5, 0.625, 0.75, 0.875, 1.0],

    y = [1.0, 1.0005, 1.0020, 1.0046, 1.0089, 1.0153, 1.0237, 1.0348, 1.0486].

    Resultados en Tabla
    Organizamos los resultados calculados para tamaños de paso h = 0.5, 0.25, 0.125
    en la siguiente tabla:

    t h = 0.5 h = 0.25 h = 0.125


    0.0 1.0000 1.0000 1.0000
    0.125 −−− −−− 1.0005
    0.25 −−− 1.0039 1.0020
    0.375 −−− −−− 1.0046
    0.5 1.0154 1.0311 1.0089
    0.625 −−− −−− 1.0153
    0.75 −−− 1.1193 1.0237
    0.875 −−− −−− 1.0348
    1.0 1.5567 1.2781 1.0486

    2
    Solución del Problema 3
    Resolveremos el problema de valor inicial (PVI):
    (
    y ′ (t) = t2 − y,
    y(0) = 1,

    utilizando el método de Runge-Kutta de orden 2, con tamaños de paso h =


    0.5, 0.25, 0.125, y para el intervalo t ∈ [0, 1]. A continuación, detallaremos paso
    a paso el proceso.

    Método de Runge-Kutta de orden 2


    El método de Runge-Kutta de orden 2 para un PVI y ′ (t) = f (t, y), y(t0 ) = y0
    se define como:
    k1 = f (tn , yn ),
    k2 = f (tn + h, yn + h · k1 ),
    h
    yn+1 = yn + · (k1 + k2 ).
    2
    Donde tn+1 = tn + h. Este proceso se repite para calcular los valores aprox-
    imados de y en cada paso.

    Cálculos para cada tamaño de paso


    Caso h = 0.5
    Para h = 0.5, calcularemos los valores de y en t = 0, 0.5, 1.0:
    • Paso 1: t0 = 0, y0 = 1
    k1 = f (0, 1) = 02 − 1 = −1,
    k2 = f (0.5, 1 + 0.5 · (−1)) = f (0.5, 0.5) = 0.52 − 0.5 = −0.25,
    0.5
    y1 = 1 + · (−1 + (−0.25)) = 1 + 0.25 · (−1.25) = 1 − 0.3125 = 0.6875.
    2
    • Paso 2: t1 = 0.5, y1 = 0.6875
    k1 = f (0.5, 0.6875) = 0.52 − 0.6875 = −0.4375,
    k2 = f (1.0, 0.6875+0.5·(−0.4375)) = f (1.0, 0.46875) = 1.02 −0.46875 = 0.53125,
    0.5
    y2 = 0.6875+ ·(−0.4375+0.53125) = 0.6875+0.25·0.09375 = 0.6875+0.0234375 = 0.7109375.
    2
    Por lo tanto, los valores aproximados para h = 0.5 son:
    t y
    0 1
    0.5 0.6875
    1.0 0.7109375

    1
    Caso h = 0.25
    Realizamos cálculos similares, pero ahora con un paso h = 0.25. Los puntos de
    evaluación serán t = 0, 0.25, 0.5, 0.75, 1.0. Calcularemos y1 , y2 , y3 , y4 :
    • Paso 1: t0 = 0, y0 = 1

    k1 = f (0, 1) = 02 − 1 = −1,

    k2 = f (0.25, 1 + 0.25 · (−1)) = f (0.25, 0.75) = 0.252 − 0.75 = −0.6875,


    0.25
    y1 = 1+ ·(−1+(−0.6875)) = 1+0.125·(−1.6875) = 1−0.2109375 = 0.7890625.
    2
    • Paso 2: t1 = 0.25, y1 = 0.7890625

    k1 = f (0.25, 0.7890625) = 0.252 − 0.7890625 = −0.7265625,

    k2 = f (0.5, 0.7890625+0.25·(−0.7265625)) = f (0.5, 0.607421875) = 0.52 −0.607421875 = −0.107421875,


    0.25
    y2 = 0.7890625+ ·(−0.7265625+(−0.107421875)) = 0.7890625+0.125·(−0.833984375) = 0.7890625−0
    2
    El procedimiento se repite para calcular y3 y y4 .

    Caso h = 0.125
    Procedemos de manera similar para h = 0.125, realizando cálculos en t =
    0, 0.125, 0.25, . . . , 1.0.

    Tabla Comparativa de Resultados


    t y (h = 0.5) y (h = 0.25) y (h = 0.125)
    0.0 1.0000 1.0000 1.0000
    0.5 0.6875 0.6848 0.6836
    1.0 0.7109 0.7031 0.7005

    También podría gustarte