1stopt使用時(shí)擬合計(jì)算不穩(wěn)定
我在用1stopt對(duì)干燥方程進(jìn)行擬合時(shí),每次計(jì)算之后所得常數(shù)都會(huì)發(fā)生變化,數(shù)值不穩(wěn)定。所用方程為經(jīng)驗(yàn)方程可能出現(xiàn)過(guò)擬合情況。請(qǐng)大佬幫忙看一下,是我軟件版本低的情況嗎?馬上快畢業(yè)了,比較急求大佬指點(diǎn)萬(wà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
返回小木蟲(chóng)查看更多
京公網(wǎng)安備 11010802022153號(hào)
數(shù)據(jù)較少,可以嘗試設(shè)置參數(shù)的范圍進(jìn)行擬合。
版本較低也可能是一個(gè)原因
如用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
,
應(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
看目標(biāo)函數(shù)殘差平方和或均方差,越小越好。