求1stOPT高版本无参数限制的大牛 帮跑一组数据
求1stOPT高版本无参数限制的大牛帮忙拟合一组数据。大牛跑完可以贴一下计算后的实部和虚部的数据,本人感激不尽!
复数拟合
拟合公式:y=y1-i*y2=a1+a2/((1+(i*x/a3)^a4)^a5)-i*((a6+i*a7)+(2*pi*x)*(a8+i*a9))/(8.85e-12*(2*pi*x).^a10)
十个参数(a1,a2,a3>0 0<a4<1, 0<a5<1, a6,a7,a8,a9无范围限制, 0<a10<1)
Data: x, 实部y1,实部y2
0.5 23.439 0.592
0.5875 23.495 0.593
0.675 23.445 0.603
0.7625 23.455 0.607
0.85 23.424 0.602
0.9375 23.376 0.642
1.025 23.386 0.647
1.1125 23.359 0.669
1.2 23.364 0.694
1.2875 23.374 0.703
1.375 23.348 0.744
1.4625 23.37 0.766
1.55 23.352 0.79
1.6375 23.364 0.839
1.725 23.386 0.872
1.8125 23.361 0.911
1.9 23.394 0.96
1.9875 23.374 0.989
2.075 23.382 1.059
2.1625 23.407 1.103
2.25 23.4 1.153
2.3375 23.423 1.215
2.425 23.443 1.276
2.5125 23.438 1.341
2.6 23.468 1.406
2.6875 23.478 1.476
2.775 23.495 1.55
2.8625 23.506 1.635
2.95 23.496 1.702
3.0375 23.507 1.815
3.125 23.516 1.901
3.2125 23.511 2.012
3.3 23.521 2.105
3.3875 23.499 2.229
3.475 23.518 2.359
3.5625 23.521 2.495
3.65 23.53 2.631
3.7375 23.509 2.76
3.825 23.471 2.928
3.9125 23.457 3.07
4 23.41 3.261
4.0875 23.408 3.413
4.175 23.36 3.592
4.2625 23.333 3.724
4.35 23.243 3.896
4.4375 23.227 4.011
4.525 23.1 4.115
4.6125 23.006 4.207
4.7 22.806 4.219
4.7875 22.684 4.217
4.875 22.511 4.165
4.9625 22.374 4.158
5.05 22.171 3.992
5.1375 22.022 3.915
5.225 21.947 3.749
5.3125 21.751 3.612
5.4 21.705 3.445
5.4875 21.544 3.289
5.575 21.532 3.143
5.6625 21.385 2.996
5.75 21.389 2.879
5.8375 21.305 2.731
5.925 21.322 2.642
6.0125 21.247 2.514
6.1 21.242 2.408
6.1875 21.171 2.279
6.275 21.188 2.234
6.3625 21.138 2.092
6.45 21.139 2.078
6.5375 21.087 1.928
6.625 21.061 1.954
6.7125 21.107 1.827
6.8 20.945 1.773
6.8875 21.075 1.72
6.975 20.918 1.644
7.0625 21.016 1.621
7.15 20.808 1.489
7.2375 20.914 1.522
7.325 20.851 1.457
7.4125 20.856 1.547
7.5 20.789 1.507
7.5875 20.778 1.612
7.675 20.814 1.6
7.7625 20.717 1.662
7.85 20.759 1.694
7.9375 20.684 1.691
8.025 20.734 1.799
8.1125 20.631 1.734
8.2 20.713 1.895
8.2875 20.547 1.79
8.375 20.704 1.992
8.4625 20.535 1.869
8.55 20.655 2.035
8.6375 20.514 1.963
8.725 20.622 2.052
8.8125 20.563 2.045
8.9 20.471 2.024
8.9875 20.549 2.164
9.075 20.418 2.077
9.1625 20.601 2.246
9.25 20.337 2.045
9.3375 20.573 2.236
9.425 20.312 2.026
9.5125 20.473 2.205
9.6 20.315 2.093
9.6875 20.374 2.155
9.775 20.289 2.071
9.8625 20.235 2.068
9.95 20.289 2.094
10.0375 20.212 2.012 @beefly
返回小木虫查看更多
京公网安备 11010802022153号
请人代跑最起码将代码写好吧。方便应助者就是方便自己。
好的,是我疏忽了,谢谢提醒!!!代码如下:
ComplexStr = i;
Variable x,y[realPart],y[imagPart];
Function y=a1+a2/((1+(i*x/a3)^a4)^a5)-i*((a6+i*a7)+(2*pi*x)*(a8+i*a9))/((8.85e-12)*(2*pi*x)^a10);
Data; //x,y[realPart],y[imagPart];
0.5 23.439 0.592
0.5875 23.495 0.593
0.675 23.445 0.603
0.7625 23.455 0.607
0.85 23.424 0.602
0.9375 23.376 0.642
1.025 23.386 0.647
1.1125 23.359 0.669
1.2 23.364 0.694
1.2875 23.374 0.703
1.375 23.348 0.744
1.4625 23.37 0.766
1.55 23.352 0.79
1.6375 23.364 0.839
1.725 23.386 0.872
1.8125 23.361 0.911
1.9 23.394 0.96
1.9875 23.374 0.989
2.075 23.382 1.059
2.1625 23.407 1.103
2.25 23.4 1.153
2.3375 23.423 1.215
2.425 23.443 1.276
2.5125 23.438 1.341
2.6 23.468 1.406
2.6875 23.478 1.476
2.775 23.495 1.55
2.8625 23.506 1.635
2.95 23.496 1.702
3.0375 23.507 1.815
3.125 23.516 1.901
3.2125 23.511 2.012
3.3 23.521 2.105
3.3875 23.499 2.229
3.475 23.518 2.359
3.5625 23.521 2.495
3.65 23.53 2.631
3.7375 23.509 2.76
3.825 23.471 2.928
3.9125 23.457 3.07
4 23.41 3.261
4.0875 23.408 3.413
4.175 23.36 3.592
4.2625 23.333 3.724
4.35 23.243 3.896
4.4375 23.227 4.011
4.525 23.1 4.115
4.6125 23.006 4.207
4.7 22.806 4.219
4.7875 22.684 4.217
4.875 22.511 4.165
4.9625 22.374 4.158
5.05 22.171 3.992
5.1375 22.022 3.915
5.225 21.947 3.749
5.3125 21.751 3.612
5.4 21.705 3.445
5.4875 21.544 3.289
5.575 21.532 3.143
5.6625 21.385 2.996
5.75 21.389 2.879
5.8375 21.305 2.731
5.925 21.322 2.642
6.0125 21.247 2.514
6.1 21.242 2.408
6.1875 21.171 2.279
6.275 21.188 2.234
6.3625 21.138 2.092
6.45 21.139 2.078
6.5375 21.087 1.928
6.625 21.061 1.954
6.7125 21.107 1.827
6.8 20.945 1.773
6.8875 21.075 1.72
6.975 20.918 1.644
7.0625 21.016 1.621
7.15 20.808 1.489
7.2375 20.914 1.522
7.325 20.851 1.457
7.4125 20.856 1.547
7.5 20.789 1.507
7.5875 20.778 1.612
7.675 20.814 1.6
7.7625 20.717 1.662
7.85 20.759 1.694
7.9375 20.684 1.691
8.025 20.734 1.799
8.1125 20.631 1.734
8.2 20.713 1.895
8.2875 20.547 1.79
8.375 20.704 1.992
8.4625 20.535 1.869
8.55 20.655 2.035
8.6375 20.514 1.963
8.725 20.622 2.052
8.8125 20.563 2.045
8.9 20.471 2.024
8.9875 20.549 2.164
9.075 20.418 2.077
9.1625 20.601 2.246
9.25 20.337 2.045
9.3375 20.573 2.236
9.425 20.312 2.026
9.5125 20.473 2.205
9.6 20.315 2.093
9.6875 20.374 2.155
9.775 20.289 2.071
9.8625 20.235 2.068
9.95 20.289 2.094
10.0375 20.212 2.012
,
均方差(RMSE): 0.4116268853424
残差平方和(SSR): 37.2760724020708
相关系数(R): 0.934350795935571
相关系数之平方(R^2): 0.873011409865435
修正R平方(Adj. R^2): 0.872124889383887
确定系数(DC): 0.843048505755613
F统计(F-Statistic): 277.416655834306
参数 最佳估算
-------------------- -------------
a1 19.9594394545544
a2 3.34963854221475
a3 -4.08919598037546
a4 1.58896355660606
a5 0.589175439529943
a6 123077987206.485
a7 -44972049724.8562
a8 2504510700.0565
a9 -3107519528.08157
a10 127.716672294193
大佬能把那个图和计算的数据给我吗?感激不尽了!!