求解一个五阶微分方程组,用ode命令失效怎么办?
画的图显示出各因变量随时间不变,怎么处理?在线等,谢谢~
本帖最后由 甘乐桑 于 2013-8-27 15:45 编辑
方程贴出来。
可能是初始条件给的不好,会求出常数也很正常。。 甘乐桑 发表于 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]