数学建模社区-数学中国

标题: 为什么用ode45求解如下方程不对？ [打印本页]

作者: Reader_Founder 时间: 2013-8-26 23:45
标题: 为什么用ode45求解如下方程不对？
function dy=cb(t,y);
dy=[-0.01*y(1)*y(4); 0.01*y(1)*y(4)-0.5*y(2); 0.5*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:24;
x0=[1, 0, 0, 0, 0];
[T,Y]=ode45('cb',ts,x0);
plot(T,Y(:,1),'-');

什么原因，在线等~谢谢

作者: Reader_Founder 时间: 2013-8-26 23:46
纠结好长时间了

作者: 逝水lcc无痕 时间: 2014-12-7 20:17
你定义的函数中多了一个后中括号！初值改一下会有不同的结果！

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