| 查看: 2890 | 回復(fù): 13 | ||
hua99dushu鐵蟲 (小有名氣)
|
[求助]
用1stopt軟件做微分方程組的參數(shù)估計(jì) 已有1人參與
|
|
|
近一個(gè)月里才知道1stopt這個(gè)軟件,可以做微分方程組的參數(shù)估計(jì)問題。 之前一直想用matlab的神經(jīng)網(wǎng)絡(luò)來做,由于水平有限,也做不出來。 網(wǎng)上下載的1stopt軟件大多是無法寫代碼或者未注冊(cè),根據(jù)一些實(shí)例寫了我的代碼,可能有錯(cuò)誤。 請(qǐng)幫忙改正后,跑一下軟件。 Title "test 1"; Parameters a1,a2,a3,b1,b2,b3,c1,c2,c3,M,C,E,N[0,]; Sharedmodel; Variable t, x,y,z; Function x'=a1*x*(y/M-1)-a2*y+a3*z; y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E); z'=c1*x*(x/N-1)-c2*y-c3*z; RowData; 0,1,2,3,4,5,6,7,8,9,10; 265,291,325,378,453,557,647,770,911,1006,1201; 7217,7344,7800,8796,10607,13266,14893,16219,17303,17677,19786; 8554,9457,10607,12443,15004,18599,21742,26018,30981,34457,41425; |
鐵桿木蟲 (職業(yè)作家)
|
感覺你的公式有些問題,似乎可以簡(jiǎn)化,下面結(jié)果參考下: 均方差(RMSE):666.824975358025 殘差平方和(SSE):13339666.4328369 相關(guān)系數(shù)(R): 0.99281958012936 相關(guān)系數(shù)之平方(R^2): 0.985690718688239 確定系數(shù)(DC): 0.752302595710559 F統(tǒng)計(jì)(F-Statistic): -0.355256071121647 參數(shù) 最佳估算 -------------------- ------------- a1 0.063244391232729 a2 26.6537700533949 a3 425.865501005945 b1 3256273550.99063 b2 1385.92982355728 b3 190853579.985187 c1 1967.65266459696 c2 7557.60987652871 c3 174107.388897955 m 24291.4733546288 c 0.10182112015089 e 1.27923766610631E22 n 0.0972369125973638 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) x 計(jì)算 x 目標(biāo) y 計(jì)算 y 目標(biāo) z 計(jì)算 z 1 1 291 295.919384351006 7344 7409.69953342258 9457 8964.92474439791 2 2 325 329.82579568465 7800 8490.90903548783 10607 10769.8722203566 3 3 378 365.520171155061 8796 9774.32316004032 12443 13050.8936628585 4 4 453 403.287476232365 10607 11128.0651594216 15004 15713.5200315771 5 5 557 443.400579346621 13266 12544.2720271355 18599 18789.0601075391 6 6 647 486.105229396296 14893 14028.6589510091 21742 22331.0481756248 7 7 770 531.633932176283 16219 15588.082085985 26018 26401.9444748202 8 8 911 580.212516031404 17303 17228.8956295004 30981 31071.3625784165 9 9 1006 632.066439410749 17677 18957.0942374736 34457 36416.4169412114 10 10 1201 687.422431267362 19786 20778.230707146 41425 42522.166618056 |
鐵蟲 (小有名氣)
鐵蟲 (小有名氣)
鐵蟲 (小有名氣)
|
Title "test 1"; Parameters a1[0,],a2[0,],a3[0,],b1[0,],b2[0,],b3[0,],c1[0,],c2[0,],c3[0,],M[8000,],C[y,],E[z,],N[300,]; Sharedmodel; Variable t, x,y,z; Function x'=a1*x*(y/M-1)-a2*y+a3*z; y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E); z'=c1*x*(x/N-1)-c2*y-c3*z; RowData; 0,1,2,3,4,5,6,7,8,9,10; 265,291,325,378,453,557,647,770,911,1006,1201; 7217,7344,7800,8796,10607,13266,14893,16219,17303,17677,19786; 8554,9457,10607,12443,15004,18599,21742,26018,30981,34457,41425; |
鐵蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
| 內(nèi)容已刪除 |
鐵蟲 (小有名氣)
鐵蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工 +4 | 安全上岸! 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|