24小时热门版块排行榜    

查看: 312  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 函评 +3 布布和一二 2026-06-13 3/150 2026-06-13 23:35 by 芋圆香飘飘
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[公派出国] T^T +5 q1zg6x4lf3 2026-06-12 5/250 2026-06-13 18:53 by lkmi35jdzt
[文学芳草园] 咸菜 +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[考研] o>_<o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[考博] T^T +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 14:45 by 80j16ujcu8
[考研] O(∩_∩)O +3 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 14:42 by 80j16ujcu8
[考研] ╭(╯ε╰)╮ +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 14:22 by 80j16ujcu8
[考博] T^T +4 jk31rb6vyz 2026-06-12 4/200 2026-06-13 12:09 by 80j16ujcu8
[考博] T^T +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[考博] →_→ +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 08:49 by t4hz8376wy
[基金申请] 面上项目没有好文章就没希望了吗? +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[考博] [急招]吉林大学招收26年9月入学考核制博士生1名,截止今晚24点前 +3 hehejin 2026-06-10 4/200 2026-06-12 18:24 by 冻雨天
[教师之家] 重庆理工大学副校长遇刺身亡 传涉案副教授疑因积怨行凶 +6 yexuqing 2026-06-07 7/350 2026-06-12 16:46 by wwwkkk83
[基金申请] 前几天时间戳更新了 +8 vito刘 2026-06-11 13/650 2026-06-12 14:17 by BlakeReary
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[考博] 还有课题组有博士名额吗 +5 冻雨天 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
请填处理意见