|
|
[交流]
1stopt擬合求助,好難啊 已有1人參與
求大神幫忙擬合一下啊,本人用matlab非線(xiàn)性擬合出來(lái)的效果不理想,據(jù)說(shuō)1stopt可以擬合帶復(fù)數(shù)的公式,我嘗試用1stopt擬合,結(jié)果定義的復(fù)數(shù)符號(hào)j也被擬合出來(lái)了
待擬合的方程式電力線(xiàn)傳輸函數(shù)幅頻特性的一個(gè)4階方程(4個(gè)只變gp和dp的表達(dá)式之和),所以我取了絕對(duì)值換算成對(duì)數(shù)使其單位變成dB,測(cè)試的數(shù)據(jù)放在附件的excel里,第一列是變量x,第二列是對(duì)應(yīng)因變量y.
就是這個(gè)函數(shù)取得對(duì)數(shù),vp=2e8,
跪求大神幫幫忙,萬(wàn)分感激。。ㄈ绻麛M合效果不好可增加階數(shù),但是待求擬合系數(shù)也會(huì)增加)
Title "H_1";
Parameters gp01,a0,a1,K,dp01,gp02,dp02,gp03,dp03,gp04,dp04;
ComplexStr = j;
Variable x , y;
Function y=20*log10(abs(gp01*exp(-(a0+a1*x^K+2*pi*j*x/2e8)*dp01)+gp02*exp(-(a0+a1*x^K+2*pi*j*x/2e8)*dp02)+gp03*exp(-(a0+a1*x^K+2*pi*j*x/2e8)*dp03)+gp04*exp(-(a0+a1*x^K+2*pi*j*x/2e8)*dp04))) ;
data;
500000 -1.582418505
519500 -1.620108522
539000 -1.657877333
558500 -1.695696353
.....
"H_1"
迭代數(shù): 0
計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:00:42:814
優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法
計(jì)算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn)
均方差(RMSE): 3.82641829064014
殘差平方和(SSE): 14656.1184118803
相關(guān)系數(shù)(R): 0.89236767319423
相關(guān)系數(shù)之平方(R^2): 0.796320064162084
決定系數(shù)(DC): 0.79632005624803
卡方系數(shù)(Chi-Square): -453.261902781698
F統(tǒng)計(jì)(F-Statistic): 3905.75409809133
參數(shù) 最佳估算
---------- -------------
gp01 0.0905058486354724
a0 -3.21061454592761
a1 0.0303224962664343
K 0.00164251742836774
dp01 0.24425464628118
gp02 0.162608273800758
dp02 0.244245680941401
gp03 0.000818566218494661
dp03 0.244738688504121
gp04 0.000143396963924432
dp04 2.87024769969935
j 18.9404153768973 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋:小木蟲(chóng)僅提供交流平臺(tái),不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶(hù)自主發(fā)布,如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題,其責(zé)任在于用戶(hù)本人,如對(duì)版權(quán)有異議,請(qǐng)聯(lián)系郵箱:xiaomuchong@tal.com
- 附件 1 : code.xlsx
2016-01-16 11:44:07, 32.59 K
» 猜你喜歡
|