1stOpt非線性擬合自定義函數(shù)(gamma函數(shù))
顯示函數(shù)表達(dá)式錯(cuò)誤,多次修改未成功,函數(shù)如附件所示
Title "Type your title here";
Parameters a,b;
Variable t,y;
Function y= t^(a-1) * exp(-t/b)*b^(-a) /(int(t^(a-1)*exp(-t),0,inf) ;
Data;
0
12 -6.08
79 -7.32
120 -8.27
144 -9.38
172 -7.26
206 -9.03
230 -8.19
255 -7.00
351 -5.10
395 -5.70
425 -6.15
445 -6.66
486 -7.28
516 -7.77
545 -7.31
573 -7.02
606 -8.52
632 -8.07
665 -4.24
699 -3.24
723 -4.84
752 -5.97
789 -6.06
815 -7.00
851 -7.22
881 -7.65
912 -6.96
940 -6.86
966 -7.80
988 -6.78

3.png
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號(hào)
Parameters a,b;
Variable t,y;
Function y= t^(a-1)*exp(-t/b)/b^a/Gamma(a);
Data;
0
12 -6.08
......
***********************
均方差(RMSE): 6.51708182055846
殘差平方和(SSR): 1316.64301913146
相關(guān)系數(shù)(R): 0.00732870945702214
相關(guān)系數(shù)之平方(R^2): 5.37099823054458E-5
修正R平方(Adj. R^2): -0.0713710250189585
確定系數(shù)(DC): -12.2545266162522
卡方系數(shù)(Chi-Square): -171.789745491722
F統(tǒng)計(jì)(F-Statistic): 26.8504257160124
參數(shù) 最佳估算
---------- -------------
a -0.01320043171917
b 1.2315680815176E308
====== 結(jié)果輸出 ======
No. 實(shí)測值 y 計(jì)算值 y
1 0 0
2 -6.08 -12.3230025991806
3 -7.32 -1.82585740146398
4 -8.27 -1.19540786357906
5 -9.38 -0.993778589151446
6 -7.26 -0.830051534835844
7 -9.03 -0.691404461863923
8 -8.19 -0.618357711796463
9 -7 -0.556975246943513
10 -5.1 -0.402937024657334
11 -5.7 -0.357495143350713
12 -6.15 -0.331939278614345
13 -6.66 -0.316828278804579
14 -7.28 -0.289762657873544
15 -7.77 -0.272700287427244
16 -7.31 -0.258003339305755
17 -7.02 -0.245233609702048
18 -8.52 -0.231707974303466
19 -8.07 -0.22205251255645
20 -4.24 -0.210891625758833
21 -3.24 -0.20050164587487
22 -4.84 -0.19375962591301
23 -5.97 -0.186190829461552
24 -6.06 -0.177346970499099
25 -7 -0.171615811595814
26 -7.22 -0.164262167534759
27 -7.65 -0.158596128921691
28 -6.96 -0.153135329168843
29 -6.86 -0.148514555472775
30 -7.8 -0.144465229327475,
謝謝大神!還有一個(gè)小問題,我該怎么定義參數(shù)a和b的范圍呢?他們都是大于0的,
模型和數(shù)據(jù)對(duì)不上,可能模型不適合。
Parameters a,b;
ParameterDomain = [0, ];
Variable t,y;
Function y= t^(a-1)*exp(-t/b)/b^a/Gamma(a);
Data;
0 0
.......
好的,謝謝大神!那么在擬合過程中,擬合出現(xiàn)了參數(shù)值、相關(guān)性和均方根誤差,但是擬合的y值都為0,是不是就是說明模型不適合,沒有擬合成功呢?
好的,謝謝大神!