24小时热门版块排行榜    

查看: 1315  |  回复: 4

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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

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

发自小木虫Android客户端
数值计算
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[基金申请] 申请科研项目纯靠实力能成功吗? 200+3 David1982 2026-10-13 6/300 2026-10-14 11:14 by 北京莱茵编辑
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见