| 3 | 1/1 | 返回列表 |
| 查看: 501 | 回復: 2 | ||
yejianling新蟲 (正式寫手)
|
[求助]
求助高版本的1stOpt代跑個程序 已有1人參與
|
|
代碼如下: Parameter p1,p2,p3,p4; Variable t,x1,x2,x3,x4; ODEFunction x2'=(x3/270)*(x1-x2)-p1*exp((-p2+p3*96485*(x4-p4))/(8.314*363))*x2^2; data; t,x1,x2,x3,x4; 0.1 165 69 100.00 361.00 0.3 115 71 110.00 356.00 0.5 103 62 112.50 358.00 0.7 150 69 118.75 380.00 0.9 128 95 112.50 380.00 1.1 140 63 100.00 383.00 1.3 110 36 100.00 383.00 1.5 103 33 118.75 370.00 1.7 96 31 137.50 367.00 1.9 70 25 100.00 356.00 2.1 83 31 125.00 373.00 2.3 90 81 137.50 363.00 2.5 100 88 125.00 353.00 2.7 112 60 137.50 346.00 2.9 119 35 137.50 353.00 3.1 66 30 112.50 353.00 3.3 70 35 125.00 373.00 3.5 160 63 112.50 370.00 3.7 147 74 112.50 370.00 3.9 140 86 112.50 373.00 4.1 108 76 107.50 373.00 4.3 109 79 75.00 363.00 4.5 76 50 125.00 353.00 4.7 79 47 131.00 363.00 4.9 82 43 112.00 362.00 5.1 136 61 135.00 373.00 5.3 145 69 125.00 363.00 5.5 143 52 110.00 366.00 5.7 159 63 125.00 370.00 5.9 178 76 100.00 370.00 6.1 112 56 100.00 360.00 6.3 121 70 81.25 364.00 6.5 152 47 112.50 363.00 6.7 157 60 100.00 373.00 6.9 162 82 93.75 383.00 7.1 170 63 118.75 363.00 7.3 62 49 112.50 291.00 7.5 83 78 112.50 356.00 7.7 227 95 137.50 359.00 7.9 226 75 125.00 360.00 8.1 223 47 125.00 355.00 8.3 116 44 143.75 366.00 8.5 130 70 143.75 370.00 |

鐵桿木蟲 (職業(yè)作家)
|
樓主的公式明顯有問題: p1*exp((-p2+p3*96485*(x4-p4))/(8.314*363)) 完全等同于: exp((-p2+p3*96485*(x4))/(8.314*363)) p1和p4多余,去掉后效果是一樣的,否則過擬合,參數(shù)組值不唯一。 計算用時(時:分:秒:微秒): 00:01:01:312 均方差(RMSE): 17.9229682919401 殘差平方和(SSE): 13491.7772805435 相關系數(shù)(R): 0.290562575469603 相關系數(shù)之平方(R^2): 0.0844266102635286 修正R平方(Adj. R^2): 0.0374741287385814 確定系數(shù)(DC): 0.0816922624187681 F統(tǒng)計(F-Statistic): 2.44181044608464 參數(shù) 最佳估算 -------------------- ------------- p2 20003.4712902898 p3 0.000151432517221469 |
新蟲 (正式寫手)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學與工程調劑 +11 | 深V宿舍吧 2026-03-30 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 求調劑 +3 | 研s. 2026-03-25 | 3/150 |
|
|
[考研] 335求調劑 +3 | 321* 2026-03-31 | 3/150 |
|
|
[考研] 一志愿085600中科院寧波所276分求調劑 +6 | 材料學257求調劑 2026-03-28 | 7/350 |
|
|
[考研] 求調劑 +3 | 圖鑒212 2026-03-30 | 3/150 |
|
|
[考研] 085701環(huán)境工程求調劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
[考研] 317求調劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
[考研] 367求調劑 +5 | 芋泥啵啵… 2026-03-28 | 5/250 |
|
|
[考研] 085601材料工程找調劑 +17 | oatmealR 2026-03-29 | 18/900 |
|
|
[考研] 求化學調劑 +11 | wulanna 2026-03-28 | 11/550 |
|
|
[考研] 348求調劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
[考研] 086000生物與醫(yī)藥調劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 316求調劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 086502化學工程342求調劑 +6 | 阿姨復古不過 2026-03-27 | 6/300 |
|
|
[考研] 調劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 考研調劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 中國科學院深圳先進技術研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 309求調劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 290分調劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-25 | 3/150 |
|