| 查看: 297 | 回复: 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 ] |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
求推荐英文EI期刊
已经有5人回复
-
最失望的一年
已经有8人回复
-
存款400万可以在学校里躺平吗
已经有27人回复
-
请教限项目规定
已经有4人回复
-
推荐一本书
已经有16人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 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