【答案】应助回帖

» 猜你喜欢

阳光

1楼2018-09-20 19:08:59

雨中奔跑

金虫 (正式写手)

应助: 6 (幼儿园)
金币: 1667.4
红花: 6
帖子: 520
在线: 291.7小时
虫号: 1899847
注册: 2012-07-18
专业: 无机非金属类光电信息与功

Title "Type your title here"
ComplexStr = i;
Variable x,y[realPart] ,y[imagPart] ;
Parameters a,b,c,d[0,1],e,f,g[0,1],h,k,l[0,1],m,n,o[0,1];
Function
y=a+(b^2)/(1+(i*2*pi*x*c)^(1-d))+(e^2)/(1+(i*2*pi*x*f)^(1-g))+(h^2)/(1+(i*2*pi*x*k)^(1-l))+(m^2)/(1+(i*2*pi*x*n)^(1-o));
17-a=(b^2)+(e^2)+(h^2)+(m^2);
Data;
1040000000 11.6592 3.2098
1120000000 11.5276 3.1712
1200000000 11.4288 3.1143
1280000000 11.329 3.0917
1360000000 11.231 3.0525
1440000000 11.1543 3.0323
1520000000 11.058 3.0114
1600000000 10.982 2.9778
1680000000 10.9125 2.9778
1760000000 10.8362 2.9415
1840000000 10.785 2.9225
1920000000 10.7126 2.9055
2000000000 10.6704 2.8747
2080000000 10.6182 2.8661
2160000000 10.5773 2.8389
2240000000 10.5489 2.8274
2320000000 10.5027 2.8186
2400000000 10.4917 2.7953
2480000000 10.4651 2.7943
2560000000 10.446 2.7695
2640000000 10.4378 2.7659
2720000000 10.414 2.7644
2800000000 10.4154 2.7531
2880000000 10.4098 2.7619
2960000000 10.4002 2.7591
3040000000 10.4015 2.7632
3120000000 10.3914 2.7789
3200000000 10.3948 2.7852
3280000000 10.3896 2.8121
3360000000 10.3793 2.8283
3440000000 10.3774 2.8464
3520000000 10.3595 2.8751
3600000000 10.3469 2.8921
3680000000 10.3316 2.9264
3760000000 10.3012 2.9473
3840000000 10.2838 2.9821
3920000000 10.2442 3.0138
4000000000 10.2155 3.0409
4080000000 10.1806 3.0829
4160000000 10.1354 3.1038
4240000000 10.1033 3.1417
4320000000 10.0461 3.1741
4400000000 10.0077 3.2004
4480000000 9.9586 3.241
4560000000 9.908 3.2592
4640000000 9.8687 3.2963
4720000000 9.8102 3.3238
4800000000 9.7696 3.3478
4880000000 9.7319 3.3801
4960000000 9.6863 3.401
5040000000 9.6633 3.4313
5120000000 9.6195 3.458
5200000000 9.5885 3.4708
5280000000 9.5667 3.5066
5360000000 9.5415 3.5154
5440000000 9.5243 3.5411
5520000000 9.4991 3.5715
5600000000 9.4985 3.5846
5680000000 9.4848 3.618
5760000000 9.472 3.6366
5840000000 9.4727 3.6674
5920000000 9.4536 3.6908
6000000000 9.4555 3.7171
6080000000 9.4456 3.7673
6160000000 9.4487 3.7746
6240000000 9.4325 3.8164
6320000000 9.4099 3.8545
6400000000 9.4149 3.9001
6480000000 9.3903 3.9319
6560000000 9.369 3.9482
6640000000 9.3489 4.0117
6720000000 9.3207 4.0436
6800000000 9.3032 4.09
6880000000 9.2669 4.1272
6960000000 9.2312 4.1595
7040000000 9.2116 4.2051
7120000000 9.1538 4.235
7200000000 9.1286 4.2686
7280000000 9.0805 4.3114
7360000000 9.0242 4.3348
7440000000 8.9918 4.3766
7520000000 8.9273 4.4089
7600000000 8.8799 4.4159
7680000000 8.8248 4.4629
7760000000 8.7661 4.4594
7840000000 8.7141 4.4651
7920000000 8.6207 4.4995
8000000000 8.5756 4.4939
8080000000 8.504 4.5157
8160000000 8.4285 4.51
8240000000 8.3682 4.5175
8320000000 8.2906 4.5381
8400000000 8.2364 4.512
8480000000 8.1507 4.52
8560000000 8.0714 4.5135
8640000000 8.0125 4.5121
8720000000 7.9271 4.5082
8800000000 7.8645 4.5065
8880000000 7.7978 4.4981
8960000000 7.7194 4.4708
9040000000 7.6638 4.4718
9120000000 7.5944 4.4655
9200000000 7.5436 4.435
9280000000 7.4844 4.4352
9360000000 7.444 4.4065
9440000000 7.4134 4.3954
9520000000 7.3516 4.3662
9600000000 7.3478 4.3547
9680000000 7.3196 4.3148
9760000000 7.2809 4.2941
9840000000 7.2807 4.2945
9920000000 7.2557 4.2549
10000000000 7.2343 4.226
10080000000 7.2279 4.2321
10160000000 7.2128 4.1977
10240000000 7.2051 4.1755
10320000000 7.1861 4.1788
10400000000 7.1784 4.1345
10480000000 7.1646 4.1463
10560000000 7.1408 4.1194
10640000000 7.1405 4.1036
10720000000 7.1093 4.129
10800000000 7.1082 4.1005
10880000000 7.0935 4.1307
10960000000 7.0672 4.133
11040000000 7.0739 4.1414
11120000000 7.032 4.1905
11200000000 7.0239 4.1807
11280000000 7.0109 4.2448
11360000000 6.9732 4.2669
11440000000 6.9939 4.2977
11520000000 6.9247 4.333
11600000000 6.9203 4.3798
11680000000 6.9111 4.4231
11760000000 6.8542 4.4322
11840000000 6.8583 4.5069
11920000000 6.772 4.5456
12000000000 6.7506 4.5541
12080000000 6.7079 4.6088
12160000000 6.6634 4.5479
12240000000 6.6526 4.593
12320000000 6.5891 4.5922
12400000000 6.5975 4.5558
12480000000 6.5271 4.6027
12560000000 6.5067 4.5452
12640000000 6.5046 4.5576
12720000000 6.4052 4.5334
12800000000 6.4364 4.477
12880000000 6.3728 4.5153
12960000000 6.3549 4.4214
13040000000 6.3586 4.446
13120000000 6.2792 4.4041
13200000000 6.3504 4.3582
13280000000 6.2785 4.363
13360000000 6.2927 4.2846
13440000000 6.3164 4.3063
13520000000 6.2558 4.2669
13600000000 6.3245 4.2256
13680000000 6.2869 4.2491
13760000000 6.3192 4.1794
13840000000 6.3548 4.2299
13920000000 6.3314 4.1556
14000000000 6.3987 4.1599
14080000000 6.3645 4.1712
14160000000 6.4232 4.099
14240000000 6.4457 4.1723
14320000000 6.4433 4.0781
14400000000 6.5009 4.1094
14480000000 6.4426 4.1376
14560000000 6.5533 4.063
14640000000 6.4958 4.1386
14720000000 6.4735 4.0682
14800000000 6.5502 4.1111
14880000000 6.4412 4.1262
14960000000 6.5073 4.075
15040000000 6.4498 4.1599
15120000000 6.3945 4.1076
15200000000 6.4377 4.1666
15280000000 6.3356 4.1649
15360000000 6.3656 4.1736
15440000000 6.3233 4.2555
15520000000 6.2634 4.2031
15600000000 6.3102 4.3022
15680000000 6.1832 4.279
15760000000 6.2274 4.3115
15840000000 6.1833 4.4105
15920000000 6.1591 4.3119
16000000000 6.1684 4.4256
16080000000 6.0823 4.3998
16160000000 6.1739 4.3946
16240000000 6.1045 4.4543
16320000000 6.1441 4.3789
16400000000 6.1645 4.4334
16480000000 6.1297 4.3888
16560000000 6.2595 4.3505
16640000000 6.1853 4.354
16720000000 6.2691 4.2791
16800000000 6.283 4.314
16880000000 6.2639 4.2019
16960000000 6.3661 4.1891
17040000000 6.3083 4.176
17120000000 6.4146 4.0681
17200000000 6.3832 4.1191
17280000000 6.3702 3.999
17360000000 6.4434 4.028
17440000000 6.346 3.9716
17520000000 6.4244 3.8973
17600000000 6.3648 3.9501
17680000000 6.3459 3.8387
17760000000 6.3668 3.9046
17840000000 6.2939 3.8368
17920000000 6.3263 3.81
18000000000 6.2377 3.8628

阳光

2楼2018-09-20 19:12:02

雨中奔跑

金虫 (正式写手)

应助: 6 (幼儿园)
金币: 1667.4
红花: 6
帖子: 520
在线: 291.7小时
虫号: 1899847
注册: 2012-07-18
专业: 无机非金属类光电信息与功

感觉没人回复，可以有偿的。

阳光

3楼2018-09-25 15:07:28

独孤神宇

版主 (知名作家)

应助: 488 (硕士)
贵宾: 0.008
金币: 29532.8
散金: 802
红花: 120
沙发: 1
帖子: 5337
在线: 850.2小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★
雨中奔跑: 金币+5, ★有帮助, 谢谢你 2018-09-25 20:27:45

// Title "Type your title here"
ComplexStr = i;
Variable x,y[realPart] ,y[imagPart] ;
Parameter a,b,c,e,f,h,k,m,n;
Parameter d=[0,1],g=[0,1],l=[0,1],o=[0,1];
Function y=a+(b^2)/(1+(i*2*pi*x*c)^(1-d))+(e^2)/(1+(i*2*pi*x*f)^(1-g))+(h^2)/(1+(i*2*pi*x*k)^(1-l))+(m^2)/(1+(i*2*pi*x*n)^(1-o));
17-a=(b^2)+(e^2)+(h^2)+(m^2);
Data;

%%

13个参数，跑了1个小时都没有结果，数据点太多了，又或者是电脑性能弱了

。神经网络考虑一下？

可以考虑每隔5个或者10个重新取一组数据试试

赞一下(1人)

» 本帖已获得的红花（最新10朵）

雨中奔跑

数值计算

4楼2018-09-25 18:21:35

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15356.8
散金: 101
红花: 234
帖子: 3409
在线: 1222.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
独孤神宇: 金币+1, 鼓励交流 2018-09-25 19:25:35
雨中奔跑: 金币+15, ★★★★★最佳答案, 谢谢你 2018-09-25 20:27:24
月只蓝: 金币+20, 感谢热心帮助！ 2018-09-26 07:52:55

参考下，确实比较费时间：

均方差(RMSE): 0.224777666596414
残差平方和(SSR): 21.5236497446252
相关系数(R): 0.969583595816403
相关系数之平方(R^2): 0.940092349276266
修正R平方(Adj. R^2): 0.939581573343302
确定系数(DC): 0.938516870471897
F统计(F-Statistic): 697.58700082705

参数最佳估算
-------------------- -------------
a -310.616619796388
b 1.59261674662206
c 3.61859362315012E-10
d 1.86183669970379E-20
e 2.60341995044174
f -1.63603943690019E-11
g 3.17836609940756E-13
h 17.7457110247239
k -3.36380940892916E-14
l 9.62892686773623E-13
m 11786555.614163
n -6156.94549402471
o 3.82364086608296E-15

赞一下(1人)

» 本帖已获得的红花（最新10朵）

雨中奔跑

5楼2018-09-25 18:32:53

雨中奔跑

金虫 (正式写手)

应助: 6 (幼儿园)
金币: 1667.4
红花: 6
帖子: 520
在线: 291.7小时
虫号: 1899847
注册: 2012-07-18
专业: 无机非金属类光电信息与功

送红花一朵

引用回帖:

5楼: Originally posted by dingd at 2018-09-25 18:32:53
参考下，确实比较费时间：

均方差(RMSE): 0.224777666596414
残差平方和(SSR): 21.5236497446252
相关系数(R): 0.969583595816403
相关系数之平方(R^2): 0.940092349276266
修正R平方(Adj. R^2): 0.93958157 ...

你好！可以把拟合的数据分享一下吗？我可以加你qq吗？我的qq1225912760

阳光

6楼2018-09-25 20:29:08

雨中奔跑

金虫 (正式写手)

应助: 6 (幼儿园)
金币: 1667.4
红花: 6
帖子: 520
在线: 291.7小时
虫号: 1899847
注册: 2012-07-18
专业: 无机非金属类光电信息与功

引用回帖:

4楼: Originally posted by 独孤神宇 at 2018-09-25 18:21:35
// Title "Type your title here"
ComplexStr = i;
Variable x,y ,y ;
Parameter a,b,c,e,f,h,k,m,n;
Parameter d=,g=,l=,o=;
Function y=a+(b^2)/(1+(i*2*pi*x*c)^(1-d))+(e^2)/(1+(i*2*p ...

神经网络怎么搞没用过，希望指点一下谢谢了

阳光

7楼2018-09-25 20:29:51

雨中奔跑

金虫 (正式写手)

应助: 6 (幼儿园)
金币: 1667.4
红花: 6
帖子: 520
在线: 291.7小时
虫号: 1899847
注册: 2012-07-18
专业: 无机非金属类光电信息与功

送红花一朵

引用回帖:

神经网络需要什么的计算机配置。

阳光

8楼2018-09-25 20:31:47

独孤神宇

版主 (知名作家)

应助: 488 (硕士)
贵宾: 0.008
金币: 29532.8
散金: 802
红花: 120
沙发: 1
帖子: 5337
在线: 850.2小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

引用回帖:

7楼: Originally posted by 雨中奔跑 at 2018-09-25 20:29:51
神经网络怎么搞没用过，希望指点一下谢谢了...

5楼已经给出结果了，你将参数代入公式不就可以计算出预测值了么。。

神经网络我基本没用过。而且也不是几句话就可以解释清楚的，matlab里面就有神经网络工具箱。

数值计算

9楼2018-09-25 20:51:10

雨中奔跑

金虫 (正式写手)

应助: 6 (幼儿园)
金币: 1667.4
红花: 6
帖子: 520
在线: 291.7小时
虫号: 1899847
注册: 2012-07-18
专业: 无机非金属类光电信息与功

引用回帖:

9楼: Originally posted by 独孤神宇 at 2018-09-25 20:51:10
5楼已经给出结果了，你将参数代入公式不就可以计算出预测值了么。。

神经网络我基本没用过。而且也不是几句话就可以解释清楚的，matlab里面就有神经网络工具箱。...

好的，谢谢你。

阳光

10楼2018-09-25 21:14:11