应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟 » beefly » 求1stOPT高版本无参数限制的大牛 帮跑一组数据
51/1返回列表
查看: 1022  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kennydota2

新虫 (著名写手)

[求助] 求1stOPT高版本无参数限制的大牛 帮跑一组数据 已有1人参与

求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
回复此楼

» 猜你喜欢
1楼 2020-06-10 19:31:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
请人代跑最起码将代码写好吧。方便应助者就是方便自己。

发自小木虫Android客户端
赞一下(2人) 回复此楼
数值计算
2楼2020-06-10 20:09:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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楼2020-06-11 09:31:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
kennydota2: 金币+100, ★有帮助 2020-06-20 13:11:15
引用回帖:
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
赞一下 回复此楼
数值计算
4楼2020-06-11 11:22:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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 ...

大佬能把那个图和计算的数据给我吗？感激不尽了！！
赞一下 回复此楼
5楼2020-06-11 13:36:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kennydota2 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭