| 查看: 337 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】这是个什么问题？【已完成】 已有2人参与
|
|||
|
在用matlab进行参数估计的时候出现如下错误 ??? Error using ==> optim\private\snls lsqnonlin cannot continue: user function is returning Inf or NaN values. [ Last edited by adu886886 on 2010-7-1 at 08:29 ] |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有3人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
求助参数估计求解
|
你好，最近一直在专研matlab求解动力学程序，但苦于基础不好，已多方求助无果，现特向您求救。 k0=[0.5 0.5 0.5 0.5 0.5] x0=[0.0510 0.309 0.56 0.064 0.151] theta=[0.0485 0.0459 0.0720 0.0686 0.1504]; V31=[1.13 1.01 1.07 1.04 1.17]; V32=[1.01 0.94 0.89 0.92 1.05] z=[0 1]; 初值 k0 = [1 1 1 1 1];x0=[0.0565 0.3017 0.5357 0.0564 0.1355]； 需要先求解动力学方程，然后采用最小二乘lsqnonlin估计。 dxdz = ... (theta*V31*k(1)*x(3)) (theta*V32*k(2)*x(4)) (theta*(k(1)+k(2)+k(3)+k(5))*x(3)) (theta*V31*k(3)*x(3)-theta*k(4)*x(4) (theta*V32*k(5)*x(3)+theta*V31*k(4)*x(4))] 结果为在反应器出口处测得的结果，即在Z=1处测得的结果 结果(五组)如下：结果与theta(i),V31(i),V32(i)向对应 exp= [ 1 0.134 0.377 0.162 0.09 0.217 2 0.132 0.390 0.281 0.076 0.168 3 0.142 0.406 0.201 0.086 0.163 4 0.152 0.422 0.243 0.08 0.152 5 0.102 0.329 0.113 0.121 0.286] 采用ODE45求解 lsqnonlin估计参数 ncdxluoqiang@163.com [ Last edited by johnl on 2010-5-30 at 14:26 ] |
3楼2010-05-30 13:56:00
