版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

976246540

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 27
帖子: 7
在線: 1.8小時
蟲號: 14375310
注冊: 2019-02-26
專業(yè): 食品加工技術(shù)

[交流] 1stopt使用時擬合計(jì)算不穩(wěn)定已有3人參與

我在用1stopt對干燥方程進(jìn)行擬合時，每次計(jì)算之后所得常數(shù)都會發(fā)生變化，數(shù)值不穩(wěn)定。所用方程為經(jīng)驗(yàn)方程可能出現(xiàn)過擬合情況。請大佬幫忙看一下，是我軟件版本低的情況嗎？馬上快畢業(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ì)算用時(時:分:秒:微秒): 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í)測值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ì)算用時(時:分:秒:微秒): 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í)測值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ì)算用時(時:分:秒:微秒): 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í)測值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ù)此樓

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
幫我的英語口語老師找學(xué)生已經(jīng)有4人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有162人回復(fù)
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有7人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有2人回復(fù)

1樓 2020-08-04 16:51:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機(jī)械動力學(xué)
管轄: 計(jì)算模擬

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

數(shù)據(jù)較少，可以嘗試設(shè)置參數(shù)的范圍進(jìn)行擬合。

版本較低也可能是一個原因

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

數(shù)值計(jì)算

2樓2020-08-04 21:25:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hzlhm

至尊木蟲 (著名寫手)

應(yīng)助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時
蟲號: 1517335
注冊: 2011-11-30
性別: GG
專業(yè): 常微分方程與動力系統(tǒng)

★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
獨(dú)孤神宇: 金幣+5, 鼓勵交流 2020-08-05 08:29:42

如用matlab做的話，結(jié)果變化不大，這個跟初值有關(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

贊一下

回復(fù)此樓

QQ：2120156492

3樓2020-08-05 00:00:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
獨(dú)孤神宇: 金幣+2, 鼓勵交流 2020-08-05 09:31:02

應(yīng)該一致的，只是順序不同而已（公式中三個部分形式是一樣的），另外指數(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

贊一下

回復(fù)此樓

4樓2020-08-05 09:27:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

976246540

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 27
帖子: 7
在線: 1.8小時
蟲號: 14375310
注冊: 2019-02-26
專業(yè): 食品加工技術(shù)

引用回帖:

4樓: Originally posted by dingd at 2020-08-05 09:27:01
應(yīng)該一致的，只是順序不同而已（公式中三個部分形式是一樣的），另外指數(shù)數(shù)值計(jì)算，有可能出現(xiàn)計(jì)算精度差異。

Root of Mean Square Error (RMSE): 0.00167798458106697
Sum of Squared Residual: 2.81563225429 ...

嗯，好像只是順序變了。但是你的a值和b值比我擬合的要大

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

5樓2020-08-05 14:26:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

976246540

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 27
帖子: 7
在線: 1.8小時
蟲號: 14375310
注冊: 2019-02-26
專業(yè): 食品加工技術(shù)

引用回帖:

3樓: Originally posted by hzlhm at 2020-08-05 00:00:17
如用matlab做的話，結(jié)果變化不大，這個跟初值有關(guān)系。
第1次計(jì)算擬合
a b c g h k
-5.48855631156569 4.4979065232128 1.92482409022818 0.469010073268637 0.4 ...

初值是什么不太懂。求賜教，謝謝

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

6樓2020-08-05 14:33:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

看目標(biāo)函數(shù)殘差平方和或均方差，越小越好。

[ 發(fā)自手機(jī)版 http://m.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

7樓2020-08-05 15:11:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hzlhm

至尊木蟲 (著名寫手)

應(yīng)助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時
蟲號: 1517335
注冊: 2011-11-30
性別: GG
專業(yè): 常微分方程與動力系統(tǒng)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

6樓: Originally posted by 976246540 at 2020-08-05 14:33:48
初值是什么不太懂。求賜教，謝謝
...

初值，就是擬合系數(shù)計(jì)算最初的值(開始的第一個值)。

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

QQ：2120156492

8樓2020-08-05 15:14:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 976246540 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 考研一志愿長安大學(xué)材料與化工309分請求調(diào)劑 +6	dtdxzxx 2026-03-06	8/400	2026-03-10 08:31 by Linda Hu
[考研] 085600材料與化工 326 求調(diào)劑 +4	熱愛生活ing 2026-03-09	4/200	2026-03-10 08:25 by Linda Hu
[考研] 315食品工程求調(diào)劑 +3	Oreov0 2026-03-06	5/250	2026-03-10 00:39 by wdwdoo
[考研] 296求調(diào)劑 +4	Xinyu Wu311 2026-03-09	4/200	2026-03-09 11:51 by 勇敢太監(jiān)王公公
[考研] 293一志愿華東理工 0817化學(xué)工程與技術(shù) 調(diào)劑 +5	fjj0912 2026-03-07	5/250	2026-03-09 09:13 by 30660438
[考研] 298求調(diào)劑 +9	fjj0912 2026-03-03	11/550	2026-03-08 22:19 by qingfeng258
[考研] 303求調(diào)劑 +8	forgman95 2026-03-05	10/500	2026-03-08 12:41 by 藍(lán)莓都給你
[考研] 材料化工求調(diào)劑 +7	博斯特525 2026-03-06	8/400	2026-03-08 03:51 by lfhuang
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3	月山斜 2026-03-06	3/150	2026-03-07 18:56 by houyaoxu
[考研] 267化工調(diào)劑求助 +7	聰少OZ 2026-03-04	7/350	2026-03-07 17:33 by 無際的草原
[考研] 求調(diào)劑推薦 +4	微辣不吃 2026-03-06	4/200	2026-03-07 00:28 by leaiy
[考博] 2026年博士名額撿漏 +4	科研ya 2026-03-04	7/350	2026-03-06 16:05 by 科研ya
[考研] 材料085601一志愿哈工大317 +4	壓迫感行 2026-03-04	4/200	2026-03-05 20:08 by 黑衣饅頭人
[考研] 347求調(diào)劑 +6	啊歐歐歐 2026-03-03	8/400	2026-03-05 11:40 by 0202liuyan
[考研] 材料328求調(diào)劑 +10	一個蘿卜02 2026-03-03	10/500	2026-03-05 09:31 by oxidpl
[考研] 化學(xué) 0703求調(diào)劑總分293 一志愿211 +9	土土小蟲 2026-03-03	11/550	2026-03-05 09:26 by vfdbghn
[考研] 085600 材料與化工 298 +14	小西笑嘻嘻 2026-03-03	14/700	2026-03-04 16:23 by sslc1985
[考研] 281求調(diào)劑 +3	我是小小蔥蔥 2026-03-03	5/250	2026-03-04 14:23 by kakakapanpan
[考研] 347分電子信息專碩求調(diào)劑 +3	咕嘟咕嘟1 2026-03-03	3/150	2026-03-04 12:32 by 熱情沙漠
[考研] 化工專碩調(diào)劑 +4	利好利好. 2026-03-03	7/350	2026-03-03 21:30 by L135790

24小時熱門版塊排行榜

[交流] 1stopt使用時擬合計(jì)算不穩(wěn)定 已有3人參與

» 猜你喜歡

[交流] 1stopt使用時擬合計(jì)算不穩(wěn)定已有3人參與