| 查看: 1117 | 回复: 4 | |||
jibing_alex铜虫 (小有名气)
|
[交流]
【求助】ode23求助 已有3人参与
|
|
大家好，请教大家一个问题： 我现在使用matlab中的ode23函数建一个微分方程，但是在运算的过程中，出现这样的warning: Warning: Failure at t=3.073417e+003. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (7.275958e-012) at time t. > In ode23 at 362 In fitness at 64 In gads\private\validate>@(x)fitness(x,FitnessFcnArgs{:}) at 135 In gads\private\fcnvectorizer at 14 In gads\private\makeState at 48 In gads\private\galincon at 18 In ga at 289 In Genetic_algoritnm at 16 In op_mian at 14 Warning: Failure at t=2.436106e+003. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (7.275958e-012) at time t. > In ode23 at 362 In fitness at 64 In gads\private\validate>@(x)fitness(x,FitnessFcnArgs{:}) at 135 In gads\private\fcnvectorizer at 14 In gads\private\makeState at 48 In gads\private\galincon at 18 In ga at 289 In Genetic_algoritnm at 16 In op_mian at 14 请大家指点一下！ |
回复此楼» 猜你喜欢
微信指数没变化，科研之友没阅读
已经有4人回复
-
这种情况还有戏吗
已经有9人回复
-
同事接到电话了，我却没有
已经有4人回复
-
准备明年的基金了
已经有5人回复
-
今年的WR进展到哪一步了？
已经有6人回复
-
时间戳他又来了
已经有11人回复
-
小木虫看见有人已经知道结果了
已经有15人回复
-
这个自发加氧反应的机理是什么？
已经有6人回复
-
J-J-W不为人知的一面
已经有18人回复
-
没消息就是被刷了呗
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】求大神，matlab程序运行出错
已经有6人回复
- 【求助】解刚性微分方程组
已经有1人回复
2楼2010-05-05 21:41:35
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
3楼2010-05-05 22:03:30
jibing_alex
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 49
- 散金: 20
- 帖子: 90
- 在线: 72.7小时
- 虫号: 939289
- 注册: 2010-01-08
- 性别: GG
- 专业: 人工智能与知识工程
|
源程序如下： 这是阶微分方程的 options = odeset('RelTol',1e-4,'AbsTol',1e-4); [t,x]=ode23(@Trabecular_Normal,[tstart tfinal],x0,options); --------------函数的定义 function xdot = Trabecular_Normal(t,x); %%----- System Variables defination % x(1) = Osteoclast % x(2) = Osteoblast % x(3) = Bone thickness global a; global b; global c; global d; global e; global f; global Koc; global Kob; Xb=0; %three equtions for the model xdot(1) = a*x(1)-b*sqrt(x(2))*x(1);%+Fob /2.1 xdot(2) = c*(x(1)^2)/(Koc+x(1)^2)*x(2)-d*x(2);%8 xdot(3) = (x(2)/(Kob+x(2)))*(e+Xb-x(3))-f*x(1)^2;%+ g*x(1); %xdot(3) = e*x(1)- f*x(2)*x(3); xdot=xdot'; % transpose to make a column vector |
4楼2010-05-05 22:33:25
jibing_alex
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 49
- 散金: 20
- 帖子: 90
- 在线: 72.7小时
- 虫号: 939289
- 注册: 2010-01-08
- 性别: GG
- 专业: 人工智能与知识工程
5楼2010-05-05 22:35:05