当前位置: 首页 > 计算模拟 >求1stOPT高版本无参数限制的大牛 帮跑一组数据

求1stOPT高版本无参数限制的大牛 帮跑一组数据

作者 kennydota2
来源: 小木虫 200 4 举报帖子
+关注

求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 返回小木虫查看更多

今日热帖
  • 精华评论
  • 独孤神宇

    请人代跑最起码将代码写好吧。方便应助者就是方便自己。

  • kennydota2

    引用回帖:
    2楼: Originally posted by 独孤神宇 at 2020-06-10 20:09:32
    请人代跑最起码将代码写好吧。方便应助者就是方便自己。

    好的,是我疏忽了,谢谢提醒!!!代码如下:

    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

  • 独孤神宇

    引用回帖:
    3楼: Originally posted by kennydota2 at 2020-06-11 09:31:51
    好的,是我疏忽了,谢谢提醒!!!代码如下:

    ComplexStr = i;
    Variable x,y,y;
    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,y;
    0.5 ...

    均方差(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

  • kennydota2

    引用回帖:
    4楼: Originally posted by 独孤神宇 at 2020-06-11 11:22:13
    均方差(RMSE): 0.4116268853424
    残差平方和(SSR): 37.2760724020708
    相关系数(R): 0.934350795935571
    相关系数之平方(R^2): 0.873011409865435
    修正R平方(Adj. R^2): 0.872124889383887
    确定系数(DC): 0.8430 ...

    大佬能把那个图和计算的数据给我吗?感激不尽了!!

猜你喜欢