数学建模社区-数学中国

标题: 求解一个五阶微分方程组，用ode命令失效怎么办？ [打印本页]

作者: 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
方程贴出来。
可能是初始条件给的不好，会求出常数也很正常。。

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=[100000000, 0, 0, 100, 0];
[t,y]=ode15s('cb',ts,x0);

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)