| 10 | 1/1 | 返回列表 |
| 查看: 1326 | 回復(fù): 9 | |||
ziliaohan銅蟲 (初入文壇)
|
[求助]
求助大神幫忙用1stopt擬合下模型參數(shù) 已有1人參與
|
|
在下小碩 想請(qǐng)大神幫忙用1stopt跑下模型參數(shù) 數(shù)據(jù) X=3091 3281 3508 3904 4108 Y=2874 3054 3271 3634 3832 Z=0.336 0.360 0.382 0.385 0.397 模型 dx/dt=-a1*x-a2*y+a3*z+a4*x*z+a5*x*y dy/dt=b1*x+b2*y-b3*z-b4*y*z dz/dt=-c1*x+c2*y+c3*z-c4*z*z 參數(shù) a1 a2 a3 a4 a5 b1 b2 b3 b4 c1 c2 c3 c4 如果可能的化,希望誤差可以<10^(-5) 感謝大神的幫助! |
鐵桿木蟲 (職業(yè)作家)
|
給你一組參考: 均方差(RMSE):0.898138086625654 殘差平方和(SSE):9.67982427177108 相關(guān)系數(shù)(R): 0.936001558277226 相關(guān)系數(shù)之平方(R^2): 0.876098917097396 確定系數(shù)(DC): -32.6666746162752 F統(tǒng)計(jì)(F-Statistic): -78192.6674901789 參數(shù) 最佳估算 -------------------- ------------- a1 0.41660768767696 a2 -0.596899303711546 a3 -1672.37950977148 a4 0.057929293839314 a5 -2.54346768480167E-6 b1 0.391376424549718 b2 -0.295297695034362 b3 1025.02219900667 b4 0.0520438936287644 c1 -0.0915214045518645 c2 -0.098248199763037 c3 -2.47400188768122 c4 -10.6412526763723 數(shù)據(jù)應(yīng)該有些問題:x、y與z大小存在數(shù)量級(jí)的差別,會(huì)嚴(yán)重影響數(shù)量級(jí)小(比如z)的擬合精度。 |
銅蟲 (初入文壇)
銅蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
銅蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
|
計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:30:03:00 均方差(RMSE):0.000238749203430912 殘差平方和(SSE):6.84014185666743E-7 相關(guān)系數(shù)(R): 0.999852515745996 相關(guān)系數(shù)之平方(R^2): 0.999705053243598 確定系數(shù)(DC): 0.999659468900762 F統(tǒng)計(jì)(F-Statistic): -6383.89529824175 參數(shù) 最佳估算 -------------------- ------------- a1 29.5516082698601 a2 -37.0355063680216 a3 -4.33534796194097 a4 3.37837345497585 a5 -18.2156770829613 b1 -31.1552076553908 b2 32.0949018446089 b3 -0.525118791751513 b4 6.63207723836531 c1 -23.5149166671026 c2 -23.5224081410208 c3 -0.976591458659144 c4 -4.870173216137 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) x 計(jì)算 x 目標(biāo) y 計(jì)算 y 目標(biāo) z 計(jì)算 z 1 1 0.1834 0.183057385327924 0.1833 0.183837598389866 0.1935 0.193615350874645 2 2 0.1961 0.196269323838988 0.1963 0.195960460145724 0.2054 0.205140005956856 3 3 0.2182 0.218285645083849 0.2181 0.218090044028013 0.207 0.207201021981877 4 4 0.2296 0.229627009142356 0.2299 0.229855468869106 0.2134 0.21335309384729 |
銅蟲 (初入文壇)
銅蟲 (初入文壇)
|
剛才把結(jié)果給老板看了下 比較滿意,和預(yù)想的有些出入 建議把模型中的最后一項(xiàng)修改成a5*y*y,會(huì)更好 即,模型為 dx/dt=-a1*x-a2*y+a3*z+a4*x*z+a5*y*y dy/dt=b1*x+b2*y-b3*z-b4*y*z dz/dt=-c1*x+c2*y+c3*z-c4*z*z 參數(shù)還是 X = 0.1728 0.1834 0.1961 0.2182 0.2296 Y = 0.1725 0.1833 0.1963 0.2181 0.2299 Z = 0.1806 0.1935 0.2054 0.2070 0.2134 能不能還請(qǐng)跑一下 要求有點(diǎn)過分 還請(qǐng)大神原諒 先謝了 |
銅蟲 (初入文壇)
|
還想請(qǐng)大神再幫忙跑一下 主要是模型中的a5x*y項(xiàng)需要修改成a5y*y, 這樣模型才符合理論推導(dǎo) 原先那個(gè)項(xiàng)推導(dǎo)錯(cuò)了 老板要求重來一次 dx/dt=-a1*x-a2*y+a3*z+a4*x*z+a5*y*y dy/dt=b1*x+b2*y-b3*z-b4*y*z dz/dt=-c1*x+c2*y+c3*z-c4*z*z 參數(shù)還是 X = 0.1728 .1834 .1961 .2182 .2296 Y = 0.1725 .1833 .1963 .2181 .2299 Z = 0.1806 .1935 .2054 .2070 .2134 感謝大神 |
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩英一數(shù)二306 +4 | z1z2z3879 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 12/600 |
|
|
[考研] 314求調(diào)劑 +8 | 無懈可擊的巨人 2026-03-12 | 8/400 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[基金申請(qǐng)] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|