| 查看: 1258 | 回复: 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 请大家指点一下！ |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有9人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
课题组会帮你通过考核么？
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】求大神，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
