| 查看: 336 | 回复: 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变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
听说今天filecode变了
已经有45人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
-
filecode
已经有19人回复
-
基金中了
已经有11人回复
-
固定端突然变了，今天
已经有8人回复
-
filecode
已经有11人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
2楼2010-05-30 13:37:14
求助参数估计求解
|
你好，最近一直在专研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
100