| 5 | 1/1 | 返回列表 |
| 查看: 916 | 回復: 4 | ||
m512990444新蟲 (初入文壇)
|
[求助]
求助1 stopt 求出參數(shù),已知編程代碼
|
|
Parameters k1[0,],k2[0,],k3[0,],k4[0,],k5[0,],k6[0,]; Constant E0=1.729; Variable t,F,S,M,GOS3,GOS4,GOS5; ODEFunction S'=-(k1*k2*(k3*M+k4*GOS3+k5*GOS4)*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); M'=-(k1*k2*k3*S*M*E0)/((k1+k2+k6 )*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); F'=(k1*k2*k8*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); GOS3'=((k3*M-k4*GOS3)*k1*k2*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); GOS4'=((k4*GOS3-k5*GOS4)*k1*k2*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); GOS5'=(k5*GOS4*k1*k2*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S); Data; 0 0 146.64 78.45 0 0 0 60 16.85 129.79 65 10.06 1.53 0 240 53.69 92.95 46.97 22.65 12.39 0 480 83.3 63.34 34.79 25.78 23.75 3.02 660 107.47 39.17 29 25.78 32.36 5.04 Data; 0 0 147.99 148.86 0 0 0 60 22.79 125.2 128.16 16.25 2.88 0 240 68.22 79.77 99.47 41 12.53 0.65 480 99.91 48.08 81.77 51.47 25.02 2.53 660 129.82 18.17 66.9 57.22 35.48 4.85 Data; 0 0.00 153.17 314.62 0.00 0.00 0.00 60 21.54 122.86 292.67 20.46 2.15 0.00 240 74.42 69.98 235.98 58.22 15.92 0.00 480 111.48 32.92 202.04 82.67 25.45 1.02 660 132.25 9.23 177.58 91.32 35.57 2.90 |
鐵桿木蟲 (職業(yè)作家)
|
供參考: 均方差(RMSE): 13.9366245195253 殘差平方和(SSE): 13984.5242158729 相關(guān)系數(shù)(R): 0.981929994660137 相關(guān)系數(shù)之平方(R^2): 0.964186514413256 確定系數(shù)(DC): 0.730266057946131 F統(tǒng)計(F-Statistic): -0.656005330812249 參數(shù) 最佳估算 -------------------- ------------- k8 0.136195429337862 k1 0.0202551743597175 k2 0.0302696407861668 k3 0.000850831443242199 k4 0.00106431484597736 k5 0.000284975127225702 k6 0.000333626838395123 |
新蟲 (初入文壇)
|
F'=(k1*k2*k8*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S)里面的k8不是參數(shù),是 k3*M+k4*GOS3+k5*GOS4,我之前簡寫了,沒有改過來, 正確的應(yīng)該是F'=(k1*k2* (k3*M+k4*GOS3+k5*GOS4) *S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S),能再幫忙擬合下嗎? 發(fā)自小木蟲Android客戶端 |
新蟲 (初入文壇)
|
不好意思,F(xiàn)'=(k1*k2*k8*S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S)里面的k8不是參數(shù),是 k3*M+k4*GOS3+k5*GOS4,我之前簡寫了,沒有改過來, 正確的應(yīng)該是F'=(k1*k2* (k3*M+k4*GOS3+k5*GOS4) *S*E0)/((k1+k2+k6)*(k3*M+k4*GOS3+k5*GOS4)+k1*k2*S),能再幫忙擬合下嗎? 發(fā)自小木蟲Android客戶端 |
新蟲 (初入文壇)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 一志愿中國海洋大學,生物學,301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 生物學調(diào)劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 考研化學學碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|