 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1934 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
976246540新蟲 (初入文壇)
|
[交流]
1stopt使用時(shí)擬合計(jì)算不穩(wěn)定 已有3人參與
|
||
|
我在用1stopt對(duì)干燥方程進(jìn)行擬合時(shí)，每次計(jì)算之后所得常數(shù)都會(huì)發(fā)生變化，數(shù)值不穩(wěn)定。所用方程為經(jīng)驗(yàn)方程可能出現(xiàn)過擬合情況。請(qǐng)大佬幫忙看一下，是我軟件版本低的情況嗎？馬上快畢業(yè)了，比較急求大佬指點(diǎn)萬分感謝。 Variable t,m; Parameters a,k,b,g,c,h; Function m=a*exp(-k*t)+b*exp(-g*t)+c*exp(-h*t); Data; 1 0.761468058 2 0.572266754 3 0.429168622 4 0.306584911 5 0.221415875 6 0.153590152 7 0.104075754 8 0.066698653 9 0.038105337 9.5 0.028177319 第一次計(jì)算擬合 "Type your title here" 迭代數(shù): 144 計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:00:05:264 優(yōu)化算法: 準(zhǔn)牛頓法(BFGS) + 通用全局優(yōu)化法 計(jì)算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 均方差(RMSE): 0.00167811805641684 殘差平方和(SSE): 2.81608021127223E-5 相關(guān)系數(shù)(R): 0.999974701263212 相關(guān)系數(shù)之平方(R^2): 0.999949403166451 決定系數(shù)(DC): 0.99994940316645 卡方系數(shù)(Chi-Square): 4.48676514663589E-5 F統(tǒng)計(jì)(F-Statistic): 158104.6611854 參數(shù) 最佳估算 ---------- ------------- a -3.91637356524864 k 0.409125168268213 b 4.93177468581318 g 0.374949664548045 c -0.0267488857371107 h -0.0180087835359877 ====== 結(jié)果輸出 ===== No 實(shí)測(cè)值m 計(jì)算值m 1 0.761468058 0.7611145 2 0.572266754 0.5741946 3 0.429168622 0.4253856 4 0.306584911 0.3095410 5 0.221415875 0.2208491 6 0.153590152 0.1538050 7 0.104075754 0.1036216 8 0.066698653 0.0663422 9 0.038105337 0.0388043 10 0.028177319 0.0278935 第二次計(jì)算擬合 "Type your title here" 迭代數(shù): 196 計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:00:05:617 優(yōu)化算法: 準(zhǔn)牛頓法(BFGS) + 通用全局優(yōu)化法 計(jì)算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 均方差(RMSE): 0.00167811998989461 殘差平方和(SSE): 2.81608670048389E-5 相關(guān)系數(shù)(R): 0.999974701204915 相關(guān)系數(shù)之平方(R^2): 0.999949403049858 決定系數(shù)(DC): 0.999949403049858 卡方系數(shù)(Chi-Square): 4.48675487738193E-5 F統(tǒng)計(jì)(F-Statistic): 158104.296839769 參數(shù) 最佳估算 ---------- ------------- a -0.0267533848791293 k -0.0179991625982107 b -3.88507793264088 g 0.409264735572213 c 4.90048287680716 h 0.374841529171177 ====== 結(jié)果輸出 ===== No 實(shí)測(cè)值m 計(jì)算值m 1 0.761468058 0.7611144 2 0.572266754 0.5741946 3 0.429168622 0.4253856 4 0.306584911 0.3095410 5 0.221415875 0.2208491 6 0.153590152 0.1538051 7 0.104075754 0.1036217 8 0.066698653 0.0663422 9 0.038105337 0.0388043 10 0.028177319 0.0278935 第三次計(jì)算擬合 "Type your title here" 迭代數(shù): 201 計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:00:05:610 優(yōu)化算法: 準(zhǔn)牛頓法(BFGS) + 通用全局優(yōu)化法 計(jì)算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 均方差(RMSE): 0.00167810849402698 殘差平方和(SSE): 2.81604811772551E-5 相關(guān)系數(shù)(R): 0.999974701551534 相關(guān)系數(shù)之平方(R^2): 0.99994940374308 決定系數(shù)(DC): 0.999949403743079 卡方系數(shù)(Chi-Square): 4.48681005484139E-5 F統(tǒng)計(jì)(F-Statistic): 158106.463142428 參數(shù) 最佳估算 ---------- ------------- a -0.0267168568795642 k -0.0180845713175149 b -4.08185733520843 g 0.408430070251824 c 5.09722825774971 h 0.37550527538769 ====== 結(jié)果輸出 ===== No 實(shí)測(cè)值m 計(jì)算值m 1 0.761468058 0.7611144 2 0.572266754 0.5741945 3 0.429168622 0.4253857 4 0.306584911 0.3095411 5 0.221415875 0.2208491 6 0.153590152 0.1538050 7 0.104075754 0.1036216 8 0.066698653 0.0663422 9 0.038105337 0.0388044 10 0.028177319 0.0278935 |
回復(fù)此樓» 猜你喜歡
幫我的英語口語老師找學(xué)生
已經(jīng)有4人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有109人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有1人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
獨(dú)孤神宇: 金幣+2, 鼓勵(lì)交流 2020-08-05 09:31:02
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
獨(dú)孤神宇: 金幣+2, 鼓勵(lì)交流 2020-08-05 09:31:02
|
應(yīng)該一致的，只是順序不同而已（公式中三個(gè)部分形式是一樣的），另外指數(shù)數(shù)值計(jì)算，有可能出現(xiàn)計(jì)算精度差異。 Root of Mean Square Error (RMSE): 0.00167798458106697 Sum of Squared Residual: 2.81563225429849E-5 Correlation Coef. (R): 0.999974705287561 R-Square: 0.999949411214945 Adjusted R-Square: 0.999934957276358 Determination Coef. (DC): 0.999949411214944 Parameter Best Estimate ---------- ------------- a 84.9030954163188 k 0.390313878805315 b -83.8883211096539 g 0.392092581899161 c -0.0261161738829758 h -0.0196386740358328 |
4樓2020-08-05 09:27:01
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時(shí)
- 蟲號(hào): 3522474
- 注冊(cè): 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動(dòng)力學(xué)
- 管轄: 計(jì)算模擬
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
數(shù)據(jù)較少，可以嘗試設(shè)置參數(shù)的范圍進(jìn)行擬合。 版本較低也可能是一個(gè)原因 發(fā)自小木蟲Android客戶端 |
2樓2020-08-04 21:25:35
hzlhm
至尊木蟲 (著名寫手)
- 應(yīng)助: 387 (碩士)
- 金幣: 17808.5
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時(shí)
- 蟲號(hào): 1517335
- 注冊(cè): 2011-11-30
- 性別: GG
- 專業(yè): 常微分方程與動(dòng)力系統(tǒng)
★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
獨(dú)孤神宇: 金幣+5, 鼓勵(lì)交流 2020-08-05 08:29:42
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
獨(dú)孤神宇: 金幣+5, 鼓勵(lì)交流 2020-08-05 08:29:42
|
如用matlab做的話，結(jié)果變化不大，這個(gè)跟初值有關(guān)系。 第1次計(jì)算擬合 a b c g h k -5.48855631156569 4.4979065232128 1.92482409022818 0.469010073268637 0.44404365198402 0.511719978797606 第2次計(jì)算擬合 -5.48855631156569 4.4979065232128 1.92482409022818 0.469010073268637 0.44404365198402 0.511719978797606 第3次計(jì)算擬合 -5.48855631156569 4.4979065232128 1.92482409022818 0.469010073268637 0.44404365198402 0.511719978797606 123.jpg |
3樓2020-08-05 00:00:17
976246540
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 27
- 帖子: 7
- 在線: 1.8小時(shí)
- 蟲號(hào): 14375310
- 注冊(cè): 2019-02-26
- 專業(yè): 食品加工技術(shù)
5樓2020-08-05 14:26:02
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[文學(xué)芳草園] 伙伴們，祝我生日快樂吧 +4 | myrtle 2026-03-10 | 5/250 |
|
|---|---|---|---|---|
|
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +7 | 11啦啦啦 2026-03-11 | 7/350 |
|
|
|
[考研] 0857 資源與環(huán)境 285分 +6 | 未名考生 2026-03-09 | 6/300 |
|
|
|
[考博] 售SCI一區(qū)T0P文章，我:8.O.5.5.1.O.5.4,科目齊全,可+急 +3 | 7x4iz20zqa 2026-03-11 | 8/400 |
|
|
|
[考研] 材料專碩288分求調(diào)劑 一志愿211 +3 | 在家想你 2026-03-11 | 3/150 |
|
|
|
[考研] 085600 材料與化工 295 求調(diào)劑 +9 | dream…… 2026-03-10 | 9/450 |
|
|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 13/650 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +7 | 小奶蒼 2026-03-07 | 7/350 |
|
|
|
[考研] 一志愿南大化學(xué)339分求調(diào)劑，四六級(jí)已過，有比賽，有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
|
[考研] 0856材料與化工290求調(diào)劑 +7 | Nebulala 2026-03-08 | 8/400 |
|
|
|
[考研] 085701環(huán)境工程317分求調(diào)劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級(jí)有競(jìng)賽求調(diào)劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
|
[考研] 材料科學(xué)（0805）338 求調(diào)劑 +7 | xiaokang3286 2026-03-07 | 7/350 |
|
|
|
[考研] 化學(xué)290求調(diào)劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
|
[考研] 第一志愿上海大學(xué)，專業(yè)化學(xué)工程與技術(shù)，總分288，求調(diào)劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
|
[考研]
 材料工程調(diào)劑
+6
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|
|
[考研] 材料與化工354調(diào)劑 +4 | Lucy-xiao 2026-03-06 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|