| 5 | 1/1 | 返回列表 |
| 查看: 2697 | 回復: 3 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
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 |

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

木蟲之王 (職業(yè)作家)
|
首先這個不是sharemodel,而且你代碼寫的也有錯。你的ode方程可以簡化形式的,把第一個X’代入第二個ode方程,把第二個新的方程和第一個方程代入第三個方程,那么整個問題就簡化成一個微分方程了。方程如下: 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 結果如下: 均方差(RMSE):1.70125017832388 殘差平方和(SSE):34.7310260309643 相關系數(shù)(R): 0.998920481960877 相關系數(shù)之平方(R^2): 0.997842129280951 確定系數(shù)(DC): 0.996942448980196 F統(tǒng)計(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 ====== 結果輸出 ====== 文件: 數(shù)據文件-1 No t 目標 s 計算 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 |

銅蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 269專碩求調劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 求調劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 一志愿華中科技大學071000,求調劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 286求調劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研] 286分人工智能專業(yè)請求調劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 0805 316求調劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 297求調劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 279分求調劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 332求調劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 265求調劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 一志愿中南化學(0703)總分337求調劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 材料學碩297已過四六級求調劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 353求調劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 085601材料工程專碩求調劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 286求調劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 收復試調劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 考研求調劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 277調劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 308求調劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|