Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 4 vistas

    Modelacion y Optimización

    Cargado por

    Bryan Suarez Mgg
    Metodo iterartivo de RK4 aplicado en el SOFTWARE MATLAB, si quieren usarlo solo copian el codigo en un nuevo script

    Copyright:

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

    Modelacion y Optimización

    Cargado por

    Bryan Suarez Mgg
    4 vistas5 páginas
    Metodo iterartivo de RK4 aplicado en el SOFTWARE MATLAB, si quieren usarlo solo copian el codigo en un nuevo script

    Título original

    modelacion y optimización

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

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

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Metodo iterartivo de RK4 aplicado en el SOFTWARE MATLAB, si quieren usarlo solo copian el codigo en un nuevo script

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    4 vistas5 páginas

    Modelacion y Optimización

    Cargado por

    Bryan Suarez Mgg
    Metodo iterartivo de RK4 aplicado en el SOFTWARE MATLAB, si quieren usarlo solo copian el codigo en un nuevo script

    Copyright:

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

    FACULTAD DE CIENCIAS DE LA INGENIERIA Y APLICADAS

    INGENIERÍA HIDRÁULICA
    MATERIA
    MODELACIÓN Y OPTIMIZACIÓN
    TRABAJO
    TEMA:
    APLICACIÓN DE RK4 PARA UN SISTEMA DE ECUACIONES

    AUTOR:
    BRYAN ISRAEL SUÁREZ SANDOVAL
    Un sistema resonante de muelles sobre el que se ejerce una fuerza externa periódica se
    modela
    mediante la ecuación
    x" (t) = 4 sin(5t) - 25x(t), x(0) = x'(0) = 0.
    1) Transformar el PVI en un sistema de ecuaciones diferenciales de primer orden
    2) Ingresar la ecuación y resolver el sistema usando RK4 para sistemas en el intervalo
    [0,2] con N=40.
    3) Represente gráficamente la solución x(t) en el intervalo mostrado
    4) Indicar en una tabla los valores de x(t) para t {0, 0.25, 0.5, 0.75, 1, 1.25, 1.50, 1.75,
    2}
    1) Para transformar a un sistema de ecuaciones consideramos que:
    ' ' ''
    y=x y =x

    Reemplazamos
    '
    y (t )=4 sin ( 5 t )−25 x (t )

    Y para MATLAB x= y ' 1 y= y '2

    Reemplazamos
    ' '
    x =y 2
    ' '
    y (t )=4 sin ( 5 t )−25∗( y 1)(t )

    Por lo tanto, tendremos nuestro sistema de ecuaciones

    2)
    function [t, y] = RK4SE(f, a, b, y0, N)
    %Método de Runge-Kutta: RK4
    h = (b - a) / N;
    t = a:h:b;
    t = t(:);
    y = zeros(length(t), length(y0));
    y(1, :) = y0;
    for k = 1:N
    k1 = f(t(k), y(k, :))';
    k2 = f(t(k) + h / 2, y(k, :) + h * k1 / 2)';
    k3 = f(t(k) + h / 2, y(k, :) + h * k2 / 2)';
    k4 = f(t(k + 1), y(k, :) + h * k3)';
    y(k + 1, :) = y(k, :) + (h / 6) * (k1 + 2 * k2 + 2 * k3 + k4);
    end
    end

    function dydt = Ejemplo4(t, y)


    dydt = [y(2); 4 * sin(5 * t) - 25 * y(1)];
    end

    % Aplicación del método


    [t2, y2] = RK4SE(@Ejemplo4, 0, 2, [0 2], 40);

    % Crear la tabla de resultados en los puntos específicos


    step = 0.25;
    evaluation_points = 0:step:2;
    evaluated_y = interp1(t2, y2, evaluation_points);

    % Mostrar los resultados en una tabla


    results_table = table(evaluation_points', evaluated_y(:,1), evaluated_y(:,2),
    'VariableNames', {'t', 'y1', 'y2'});
    disp(results_table)

    % Graficar los resultados


    plot(t2, y2)
    legend('y_1', 'y_2')

    3) GRÁFICA DE LA SOLUCIÓN DE x(t)

    4) TABLA DE RESULTADOS PARA UN t con los intervalos que se muestran en la


    tabla
    t y1 y2
    0 0 2
    0.25 0.42398 1.1052
    0.50 0.44751 -1.0037
    0.75 -0.028144 -2.4984
    1 -0.57372 -1.3508
    1.25 -0.5157 1.9156
    1.50 0.24216 3.5074
    1.75 0.84644 0.62547
    2 0.41026 -3.8538

    TABLA CREADA EN MATLAB

    También podría gustarte