Решение системы ДУ 4-го порядка методом Рунге-Кутта

@Ksenia1009 · Регистрация: 19.05.2018

Студворк — интернет-сервис помощи студентам

Необходимо решить данную систему методом Рунге-Кутта. Шаг по времени 0.001.
Помогите, пожалуйста, реализовать в matlab. Уже который день сижу, ничего не выходит..

Заранее спасибо.

@SSC · 08.11.2018, 07:55

Вот вариант в одном файле

Matlab M

function mathieu
    clear all; clc; close all;
    global Eps Psi k Betta
    Eps=1.6;
    Psi=0.2;
    k=1.8;
    Betta=0.6;
  
    dt=0.001;
    t=[0:dt:1.9];
    options=odeset('RelTol',1e-6);
    [T,Y]=ode45(@syst, t, [0 0.1 0 0.2],options);
 
    plot(T,Y(:,1),'r',T,Y(:,2),'g',T,Y(:,3),'b',T,Y(:,4),'y');
    grid on
    legend('fi1','y1','fi2','y2');
end
 
function Dy=syst(t,y)
    global Eps Psi k Betta
    
    Dy=[y(2);
        -Eps*y(2)-sin(y(1))+Psi+k*sin(y(3)-y(1)-Betta);
        y(4);
        -Eps*y(4)-sin(y(3))+Psi+k*sin(y(1)-y(3)-Betta)];
end

@Ksenia1009 0 / 0 / 0 Регистрация: 19.05.2018 Сообщений: 2
		1
	Решение системы ДУ 4-го порядка методом Рунге-Кутта 07.11.2018, 21:16. Показов 2228. Ответов 1 Метки нет (Все метки) Необходимо решить данную систему методом Рунге-Кутта. Шаг по времени 0.001. Помогите, пожалуйста, реализовать в matlab. Уже который день сижу, ничего не выходит.. Заранее спасибо. Миниатюры 0

Опции темы