Reader_Founder 发表于 2013-8-27 13:11

求解一个五阶微分方程组,用ode命令失效怎么办?

画的图显示出各因变量随时间不变,怎么处理?
在线等,谢谢~

甘乐桑 发表于 2013-8-27 15:44

本帖最后由 甘乐桑 于 2013-8-27 15:45 编辑

方程贴出来。
可能是初始条件给的不好,会求出常数也很正常。。

Reader_Founder 发表于 2013-8-27 18:47

甘乐桑 发表于 2013-8-27 15:44 static/image/common/back.gif
方程贴出来。
可能是初始条件给的不好,会求出常数也很正常。。

function dy=cb(t,y);

dy=[-0.01*y(1)*y(4); 0.01*y(1)*y(4)-0.021*y(2); 0.021*y(2)-0.55*y(3); 0.5*y(3)-0.1*y(4); 0.55*y(3)+0.1*y(4)];


end

clear;
ts=0:3:72;
x0=;
=ode15s('cb',ts,x0);
页: [1]
查看完整版本: 求解一个五阶微分方程组,用ode命令失效怎么办?