| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 917 | 回復: 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 相關系數(shù)(R): 0.981929994660137 相關系數(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,我之前簡寫了,沒有改過來, 正確的應該是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,我之前簡寫了,沒有改過來, 正確的應該是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ā)表 | |
|---|---|---|---|---|
|
[考研] 求調劑院校信息 +4 | CX 330 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 考研化學學碩調劑,一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|
[考研] 297求調劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 317求調劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 306求0703調劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 265求調劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 303求調劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研] 求調劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 311求調劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 261求B區(qū)調劑,科研經歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 086500 325 求調劑 +3 | 領帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 0856調劑,是學校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 333求調劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 中科院材料273求調劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 0856專碩279求調劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|