24小时热门版块排行榜

返回列表

duandacheng1

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 2
帖子: 2
在线: 24分钟
虫号: 3189795
注册: 2014-05-08

[交流] 1stopt拟合求助，好难啊已有1人参与

求大神帮忙拟合一下啊，本人用matlab非线性拟合出来的效果不理想，据说1stopt可以拟合带复数的公式，我尝试用1stopt拟合，结果定义的复数符号j也被拟合出来了

待拟合的方程式电力线传输函数幅频特性的一个4阶方程（4个只变gp和dp的表达式之和），所以我取了绝对值换算成对数使其单位变成dB，测试的数据放在附件的excel里，第一列是变量x，第二列是对应因变量y.
$\"1stopt拟合求助，好难啊\"$ 就是这个函数取得对数，vp=2e8,
跪求大神帮帮忙，万分感激啊！（如果拟合效果不好可增加阶数，但是待求拟合系数也会增加）

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"

迭代数: 0
计算用时(时:分:秒:微秒): 00:00:42:814
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 3.82641829064014
残差平方和(SSE): 14656.1184118803
相关系数(R): 0.89236767319423
相关系数之平方(R^2): 0.796320064162084
决定系数(DC): 0.79632005624803
卡方系数(Chi-Square): -453.261902781698
F统计(F-Statistic): 3905.75409809133

参数最佳估算
---------- -------------
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

回复此楼