| 5 | 1/1 | 返回列表 |
| 查看: 2699 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
zhl717銅蟲 (小有名氣)
|
[求助]
1st opt擬合方程——————求助! 已有1人參與
|
|
|
有高版本1st opt的大神幫忙跑下,謝謝!!! Test1; Parameters a,b,c,d,e,f,g,h; SharedModel; Variable t,X,P,S; ODEFunction X'=a*X*(1-X/b)*(1+S/c)^(-1); P'=d*X'+e*X; S'=f*X'+g*P'+h*X; Data; t,X,P,S; 0 0.35 0 99.8 12 2.75 0.55 88.96 24 4.6 3.72 73.49 36 6.27 8.19 61.57 48 8.07 14.12 47.38 60 9.73 19.37 33.63 72 10.41 22.54 24.65 84 10.7 27.61 11.19 96 10.53 32.49 3.2 108 10.56 35.62 0 120 10.59 37.58 0 132 10.4 38.44 0 144 10.77 39.88 0 |

銅蟲 (小有名氣)

木蟲之王 (職業(yè)作家)
|
首先這個(gè)不是sharemodel,而且你代碼寫的也有錯(cuò)。你的ode方程可以簡(jiǎn)化形式的,把第一個(gè)X’代入第二個(gè)ode方程,把第二個(gè)新的方程和第一個(gè)方程代入第三個(gè)方程,那么整個(gè)問題就簡(jiǎn)化成一個(gè)微分方程了。方程如下: S'=f*(a*X*(1-X/b)*(1+S/c)^(-1))+g*(d*(a*X*(1-X/b)*(1+S/c)^(-1))+e*X)+h*X 代碼改成如下: Parameters a,b,c,d,e,f,g,h; Variable t,X,P,S; ODEFunction S'=f*(a*X*(1-X/b)*(1+S/c)^(-1))+g*(d*(a*X*(1-X/b)*(1+S/c)^(-1))+e*X)+h*X; Data; 0 0.35 0 99.8 12 2.75 0.55 88.96 24 4.6 3.72 73.49 36 6.27 8.19 61.57 48 8.07 14.12 47.38 60 9.73 19.37 33.63 72 10.41 22.54 24.65 84 10.7 27.61 11.19 96 10.53 32.49 3.2 108 10.56 35.62 0 120 10.59 37.58 0 132 10.4 38.44 0 144 10.77 39.88 0 結(jié)果如下: 均方差(RMSE):1.70125017832388 殘差平方和(SSE):34.7310260309643 相關(guān)系數(shù)(R): 0.998920481960877 相關(guān)系數(shù)之平方(R^2): 0.997842129280951 確定系數(shù)(DC): 0.996942448980196 F統(tǒng)計(jì)(F-Statistic): 186.89093582654 參數(shù) 最佳估算 -------------------- ------------- a -0.307553362555249 b 11.7009599853419 c -0.0108960231582921 d 17354.8847310689 e -0.526630322610445 f 0.0237395701129598 g -0.824893382148287 h -0.326710215811974 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) s 計(jì)算 s 1 12 88.96 90.6377657133101 2 24 73.49 77.4840744651627 3 36 61.57 61.5302531902025 4 48 47.38 45.0285350315543 5 60 33.63 33.8481939424802 6 72 24.65 25.1096267702726 7 84 11.19 11.1348638405025 8 96 3.2 0.0109328172964804 9 108 0 0.0109000671739397 10 120 0 0.0109000671739397 11 132 0 0.0109000671739397 12 144 0 0.0109000671739397 |

銅蟲 (小有名氣)
|
參數(shù)值也不太合理,可能是我方程有問題,你在幫我擬合一下這個(gè)吧,感激不盡啊 ![]() http://www.gaoyang168.com/bbs/viewthread.php?tid=9291567 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +4 | Zhao anqi 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗(yàn)) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|