| 查看: 931 | 回复: 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 请大家指点一下！ |
回复此楼» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有10人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有9人回复
-
今年春晚有几个节目很不错，点赞！
已经有10人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有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