您所在位置：主题：求帮忙龙格库塔法解微分方程组
fxlzzy
[专家分：30] 发布于 2009-02-20 17:24:00
大家帮帮忙啊，帮我看看 这方程组能用ODE45解吗
function dy=xitong(t,y)
clc;clear;
dy=zeros(8,1);
y=zeros(8,1);
dy(1)=0.7*y(1)-0.45*y(2)-3.96*y(5)-1.38*y(8)+0.20;
dy(2)=-0.15*y(1)-0.18*y(2)+0.1*y(5)-0.30;
dy(3)=0.60*y(3)-0.43*y(4)-4.08*y(6)-1.46*y(8)+0.40;
dy(4)=-0.16*y(3)-0.18*y(4)+0.1*y(6)-0.30;
dy(5)=0.7*y(1)-0.13*y(2)-1.03*y(5)+0.57*y(8);
dy(6)=0.73*y(3)-0.14*y(4)-1.08*y(6)+0.61*y(8);
dy(7)=-0.003*y(7)-0.009*y(8);
dy(8)=0.004*y(7)-0.001*y(5)-0.001*y(6);
我用ODE45去解，结果怎么会单调变化的呢，他应该是波动的啊，初值都是0，有问题请大家帮忙指出来，谢谢大家了
function dy=xitong(t,y)
clc;clear;
dy=zeros(8,1);
y=zeros(8,1);
dy(1)=0.7*y(1)-0.45*y(2)-3.96*y(5)-1.38*y(8)+0.20;
dy(2)=-0.15*y(1)-0.18*y(2)+0.1*y(5)-0.30;
dy(3)=0.60*y(3)-0.43*y(4)-4.08*y(6)-1.46*y(8)+0.40;
dy(4)=-0.16*y(3)-0.18*y(4)+0.1*y(6)-0.30;
dy(5)=0.7*y(1)-0.13*y(2)-1.03*y(5)+0.57*y(8);
dy(6)=0.73*y(3)-0.14*y(4)-1.08*y(6)+0.61*y(8);
dy(7)=-0.003*y(7)-0.009*y(8);
dy(8)=0.004*y(7)-0.001*y(5)-0.001*y(6);
我用ODE45去解，结果怎么会单调变化的呢，他应该是波动的啊，初值都是0，有问题请大家帮忙指出来，谢谢大家了