| 查看: 935 | 回复: 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 请大家指点一下！ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】求大神，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
