| 查看: 939 | 回复: 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 请大家指点一下！ |
回复此楼» 猜你喜欢
一志愿北交大材料工程总分358求调剂
已经有6人回复
-
312求调剂
已经有7人回复
-
材料334求调剂
已经有19人回复
-
生物与医药273求调剂
已经有6人回复
-
301求调剂
已经有10人回复
-
专硕0854初试考材科基，求调剂
已经有5人回复
-
0703化学
已经有15人回复
-
0703调剂，一志愿天津大学319分
已经有9人回复
-
085600材料与化工301分求调剂院校
已经有10人回复
-
求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】求大神，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
