| 3 | 1/1 | 返回列表 |
| 查看: 502 | 回復(fù): 2 | ||
yejianling新蟲 (正式寫手)
|
[求助]
求助高版本的1stOpt代跑個(gè)程序 已有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ù)組值不唯一。 計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:01:01:312 均方差(RMSE): 17.9229682919401 殘差平方和(SSE): 13491.7772805435 相關(guān)系數(shù)(R): 0.290562575469603 相關(guān)系數(shù)之平方(R^2): 0.0844266102635286 修正R平方(Adj. R^2): 0.0374741287385814 確定系數(shù)(DC): 0.0816922624187681 F統(tǒng)計(jì)(F-Statistic): 2.44181044608464 參數(shù) 最佳估算 -------------------- ------------- p2 20003.4712902898 p3 0.000151432517221469 |
新蟲 (正式寫手)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 290分調(diào)劑求助 +4 | 吉祥止止陳 2026-03-25 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 08工科,295,接受跨專業(yè)調(diào)劑 +6 | lmnlzy 2026-03-30 | 6/300 |
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 3/150 |
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +10 | 深V宿舍吧 2026-03-29 | 10/500 |
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
[考研] 292求調(diào)劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué),化工,生物,生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入! +3 | llss0711 2026-03-28 | 6/300 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|