Решение дифференцированного уравнения

@Vaha 1995 · Регистрация: 13.11.2012

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

Здраствуйте!!!Вы бы не могли мене помочь в решении еще одного диф.уравнения.....у меня всё работает, но чуть чуть некоректно!!!
Вот диф.ур.:7*у(4)+6*у(3)-5*у(2)+4*у(1)+3*у+21=tg(t) Начальные условия [2 2 2 0], время [1 22]
y(4)-это порядок диф.ура
И вот, то что я сделал, но оно не работает на всьом интервале!!!
Файл-функція

Matlab M

1
2
3

function F=dif(t,y) 
 
 F=[y(2); y(3); y(4); -(6/7)*y(4)+(5/7)*y(3)-(4/7)*y(2)-(3/7)*y(1)-3+(1/7)*tan(t)];

Файл-розв’язок

Matlab M

options=odeset('RelTol', 1.0e-03);
 Time=[1 22]; 
 Y0=[2 2 2 0]; 
 [T,Y]=ode45('dif', Time, Y0, options);
 plot(T,Y(:,1),'yd--','LineWidth',2);
 hold 
 grid 
 xlabel('\bf\fontname{Times New Roman}\fontsize{14}Час t, с'); 
 ylabel('\bf\fontname{Times New Roman}\fontsize{14}y1,y2,y3,y4');
 title('\bf\fontname{Times New Roman}\fontsize{14}Рішення диференційного рівняння четвертого порядку ');
 options=odeset('RelTol',1.0e-03); 
 [T,Y]=ode45('dif',Time,Y0,options);
 plot(T,Y(:,2),'yd--','LineWidth',2);
 
 grid;
 options=odeset('RelTol',1.0e-03);
 [T,Y]=ode45('dif',Time,Y0,options);
 plot(T,Y(:,3),'yd--','LineWidth',2);
 
 grid;
 options=odeset('RelTol',1.0e-03);
 [T,Y]=ode45('dif',Time,Y0,options);
 plot(T,Y(:,4),'yd--','LineWidth',2);
 legend('y1','y2','y3','y4',3)

Зосима · 21.04.2013, 19:54

Не знаю, у меня что-то нарисовалось

Зосима 5243 / 3571 / 379 Регистрация: 02.04.2012 Сообщений: 6,473 Записей в блоге: 17
	21.04.2013, 19:54	2
	Не знаю, у меня что-то нарисовалось 0