24小时热门版块排行榜    

查看: 279  |  回复: 1

xwndf250

银虫 (小有名气)

[求助] 用ode45求解常微分方程出现错误

错误如下:Warning: Failure at t=1.400456e+001.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.842171e-014) at time t.
m文件:
function dy=system(t,y)
a=0.1;b=0.5;r=0.1;ka=0.5;kir=0.01;kil=0.01;kls=0.1;
dy=[r*b*y(1)^2-r*(2*a-b)*y(2)*y(1)+r*b*(y(3)+y(4))*y(1)+kls*y(4);
-r*b*y(1)^2-r*(-2*a+b)*y(2)*y(1)-r*b*(y(3)+y(4))*y(1)-(r-ka)*y(1);
r*y(1)+kir*y(3);
ka*y(1)-kil*y(3)-kir*y(3);
kil*y(3)-kls*y(4)];
end
命令:
[t,x]=ode45('system',[0,50],[1/6,5/6,0,0,0]);

是不是不能用ode45?我用的是matlab。。。求帮助~

sirm1.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

柴柴湖

银虫 (小有名气)

【答案】应助回帖

★ ★
xwndf250: 金币+1, 有帮助 2013-02-21 23:50:59
fegg7502: 金币+1, 3ks 2013-02-22 12:02:44
一般不会是ode45的事情,多半是方程或者参数设置不合理
2楼2013-02-21 19:47:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xwndf250 的主题更新
信息提示
请填处理意见