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

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

24小時(shí)熱門版塊排行榜

返回列表

976246540

新蟲 (初入文壇)

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

[交流] 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論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有254人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會(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ù)

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小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

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

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

版本較低也可能是一個(gè)原因

發(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小時(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

如用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

贊一下

回復(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小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

★ ★ ★
小木蟲: 金幣+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

贊一下

回復(fù)此樓

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

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

976246540

新蟲 (初入文壇)

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

引用回帖:

4樓: Originally posted by dingd at 2020-08-05 09:27:01
應(yīng)該一致的，只是順序不同而已（公式中三個(gè)部分形式是一樣的），另外指數(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小時(shí)
蟲號(hào): 14375310
注冊(cè): 2019-02-26
專業(yè): 食品加工技術(shù)

引用回帖:

3樓: Originally posted by hzlhm at 2020-08-05 00:00:17
如用matlab做的話，結(jié)果變化不大，這個(gè)跟初值有關(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小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

看目標(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小時(shí)
蟲號(hào): 1517335
注冊(cè): 2011-11-30
性別: GG
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

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

引用回帖:

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

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

發(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 的主題更新

返回列表

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

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

[考研] 22408 275分求調(diào)劑 +4	宇智波比 2026-03-03	5/250	2026-03-10 14:50 by 在風(fēng)落中
[考研] 278求調(diào)劑 +8	Gale1314 2026-03-06	8/400	2026-03-10 14:18 by Linda Hu
[考研] 收調(diào)劑 +4	調(diào)劑的考研學(xué)生 2026-03-10	4/200	2026-03-10 14:09 by 賈3723
[考研] 求調(diào)劑一志愿蘇州大學(xué)，0856化工323分 \| 本科應(yīng)化 \| 有專利/競(jìng)賽/科研助手經(jīng)歷 \| +7	橙子cyx 2026-03-06	9/450	2026-03-10 12:03 by kiokin
[考研] 材料調(diào)劑，307分 +10	張泳銘1 2026-03-09	10/500	2026-03-10 11:47 by peike
[考研] 304求調(diào)劑（085602一志愿985） +6	化工人999 2026-03-09	6/300	2026-03-10 11:40 by peike
[考研] 化工284求調(diào)劑 +12	小黑想變白 2026-03-07	12/600	2026-03-10 09:17 by Wahxp
[碩博家園] 2026級(jí)碩士研究生招生/調(diào)劑 +4	lbj6746988 2026-03-03	7/350	2026-03-10 08:15 by lbj6746988
[考研] 334求調(diào)劑 +5	Trying] 2026-03-07	5/250	2026-03-09 14:47 by ruiyingmiao
[考研] 286求調(diào)劑 +12	Faune 2026-03-06	14/700	2026-03-08 14:54 by hmn_wj
[考研] 085701環(huán)境工程專業(yè)，初試305，均過國家A區(qū)線 +7	卡卡來了@ 2026-03-07	8/400	2026-03-08 06:40 by 劉兵
[考研] 346分材料求調(diào)劑 +5	snow_反季節(jié)版 2026-03-07	5/250	2026-03-07 22:40 by Leeding1356
[考研] 276求調(diào)劑本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4	路lyh123 2026-03-06	6/300	2026-03-07 11:37 by shdliugang
[考研] 化學(xué)專業(yè)調(diào)劑 +5	好好好1233 2026-03-04	6/300	2026-03-06 15:38 by @颯颯颯颯
[考研] 287求調(diào)劑 +3	看看我. 2026-03-05	6/300	2026-03-06 10:40 by Iveryant
[考研] 271求調(diào)劑 +7	月色c 2026-03-05	8/400	2026-03-05 19:51 by wangjihu
[考研] 301求調(diào)劑 +4	李LJR 2026-03-04	4/200	2026-03-05 13:31 by Iveryant
[考研] 274環(huán)境工程求調(diào)劑 +6	扶柳盈江 2026-03-05	6/300	2026-03-05 13:16 by 夢(mèng)天888
[考研] 一志愿314求調(diào)劑 +7	202111120625 2026-03-03	7/350	2026-03-04 15:56 by zhukairuo
[考研] 298求調(diào)劑一志愿中海洋 +3	lour. 2026-03-03	3/150	2026-03-03 20:41 by gxg2025

24小時(shí)熱門版塊排行榜

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

» 猜你喜歡

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