| 查看: 349 | 回复: 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 ] |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有24人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有18人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
-
重要来源：本周末出结果
已经有10人回复
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
