| 5 | 1/1 | 返回列表 |
| 查看: 309 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】這是個(gè)什么問(wèn)題?【已完成】 已有2人參與
|
|||
|
在用matlab進(jìn)行參數(shù)估計(jì)的時(shí)候出現(xiàn)如下錯(cuò)誤 ??? 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 ] |
|
你好,最近一直在專研matlab求解動(dòng)力學(xué)程序,但苦于基礎(chǔ)不好,已多方求助無(wú)果,現(xiàn)特向您求救。 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]; 需要先求解動(dòng)力學(xué)方程,然后采用最小二乘lsqnonlin估計(jì)。 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))] 結(jié)果為在反應(yīng)器出口處測(cè)得的結(jié)果,即在Z=1處測(cè)得的結(jié)果 結(jié)果(五組)如下:結(jié)果與theta(i),V31(i),V32(i)向?qū)?yīng) 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估計(jì)參數(shù) ncdxluoqiang@163.com [ Last edited by johnl on 2010-5-30 at 14:26 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 289求調(diào)劑 +13 | 新時(shí)代材料 2026-03-27 | 13/650 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|