24小时热门版块排行榜    

查看: 1069  |  回复: 3

clownxxl

新虫 (小有名气)

[求助] 1stopt程序运行 已有1人参与

想拟合复数函数中的各个参数值,求大神帮忙平跑一下程序,v是自变量,epsilon_re是函数实部,epsilon_im是函数虚部,我编程使用的是手册中介绍的分解式的拟合方法,下面是写好的程序代码:
Variable v,epsilon_re,epsilon_im;
SharedModel;
Function epsilon_re = epsmax-w_pe1^2/((600*pi/v)^2+r_e1^2)+derta*w_pe2^2*(w_oe^2-(600*pi/v)^2)/((w_oe^2-(600*pi/v)^2)^2+(r_e2*(600*pi/v))^2);
epsilon_im=w_pe1^2*r_e1/((600*pi/v)^3+r_e1^2*(600*pi/v))+derta*w_pe2^2*r_e2*(600*pi/v)/((w_oe^2-(600*pi/v)^2)^2+(r_e2*(600*pi/v))^2);   
Data;
4        -79.1178        8.23761
4.02222        -80.1449        8.5075
4.04444        -82.0825        8.40739
4.06667        -83.6393        9.23381
4.08889        -84.8806        9.10384
4.11111        -87.0815        9.41879
4.13333        -88.3183        10.1001
4.15556        -90.0435        9.79016
4.17778        -92.3154        10.5519
4.2        -93.4154        11.0448
4.22222        -95.5231        10.7695
4.24444        -97.7177        11.8124
4.26667        -98.8034        12.0807
4.28889        -101.257        11.7945
4.31111        -103.588        13.0228
4.33333        -104.767        13.3184
4.35556        -107.417        13.0851
4.37778        -109.919        14.4726
4.4        -111.115        14.899
4.42222        -113.831        14.5706
4.44444        -116.701        15.9846
4.46667        -117.974        16.7
4.48889        -120.65        16.1836
4.51111        -124.187        17.4619
4.53333        -125.783        18.8641
4.55556        -128.08        18.4646
4.57778        -132.006        19.2723
4.6        -134.24        21.2338
4.62222        -136.16        21.1082
4.64444        -140.402        21.191
4.66667        -143.997        23.699
4.68889        -145.505        24.7699
4.71111        -149.066        24.1406
4.73333        -154.281        25.9221
4.75556        -156.672        28.8856
4.77778        -158.795        28.7495
4.8        -164.255        28.6672
4.82222        -169.551        32.2478
4.84444        -171.259        34.8206
4.86667        -174.79        34.081
4.88889        -181.992        35.2695
4.91111        -186.99        40.4723
4.93333        -188.619        42.5843
4.95556        -193.887        41.7276
4.97778        -202.998        44.6473
5        -207.777        51.88
5.02222        -209.404        53.9681
5.04444        -216.512        53.1831
5.06667        -228.333        58.217
5.08889        -233.314        69.0166
5.11111        -234.418        72.2202
5.13333        -243.105        72.124
5.15556        -258.021        80.4418
5.17778        -262.249        96.9575
5.2        -261.398        102.141
5.22222        -270.623        103.63
5.24444        -287.108        116.436
5.26667        -288.464        138.895
5.28889        -284.299        146.413
5.31111        -291.433        150.212
5.33333        -305.267        165.207
5.35556        -305.505        189.387
5.37778        -299.086        200.509
5.4        -302.175        206.562
5.42222        -312.792        220.077
5.44444        -317.37        244.714
5.46667        -311.056        264.572
5.48889        -308.418        276.665
5.51111        -313.404        291.239
5.53333        -319.548        315.848
5.55556        -316.074        346.506
5.57778        -306.184        371.954
5.6        -299.343        394.687
5.62222        -295.56        423.091
5.64444        -286.818        461.573
5.66667        -264.483        503.8
5.68889        -231.941        540.58
5.71111        -195.46        573.387
5.73333        -153.34        607.112
5.75556        -97.4825        641.048
5.77778        -22.5489        663.196
5.8        61.7332        660.105
5.82222        139.56        632.392
5.84444        203.672        589.706
5.86667        255.714        538.738
5.88889        296.274        479.126
5.91111        319.506        413.018
5.93333        323.411        350.362
5.95556        315        298.804
5.97778        302.618        258.075
6        289.911        223.756
6.02222        275.875        192.81
6.04444        259.082        165.891
6.06667        240.766        144.579
6.08889        223.331        128.527
6.11111        208.105        115.838
6.13333        194.775        104.617
6.15556        182.14        94.1435
6.17778        169.46        84.9086
6.2        157.042        77.4912
6.22222        145.576        71.7663
6.24444        135.405        67.0798
6.26667        126.324        62.7677
6.28889        117.814        58.5607
6.31111        109.484        54.6394
6.33333        101.329        51.2995
6.35556        93.5821        48.6175
6.37778        86.4301        46.4231
6.4        79.8729        44.4587
6.42222        73.7474        42.5481
6.44444        67.8536        40.6797
6.46667        62.0888        38.9587
6.48889        56.4787        37.4833
6.51111        51.1055        36.2663
6.53333        46.0258        35.2446
6.55556        41.2362        34.3312
6.57778        36.6829        33.4632
6.6        32.2975        32.6255
6.62222        28.0332        31.8415
6.64444        23.8795        31.1429
6.66667        19.8508        30.5457
6.68889        15.9662        30.0433
6.71111        12.2364        29.6138
6.73333        8.65784        29.23
6.75556        5.21472        28.8681
6.77778        1.88333        28.5142
6.8        -1.36104        28.1656
6.82222        -4.53808        27.8294
6.84444        -7.6589        27.5183
6.86667        -10.7248        27.2452
6.88889        -13.7282        27.0184
6.91111        -16.6567        26.8369
6.93333        -19.4999        26.6877
6.95556        -22.258        26.5498
6.97778        -24.9464        26.4032
7        -27.5914        26.2408
7.02222        -30.2192        26.0731
7.04444        -32.8434        25.925
7.06667        -35.4577        25.8246
7.08889        -38.0367        25.7895
7.11111        -40.5465        25.8141
7.13333        -42.9635        25.8678
7.15556        -45.2913        25.9077
7.17778        -47.5629        25.9024
7.2        -49.8255        25.8494
7.22222        -52.1168        25.7778
7.24444        -54.447        25.7349
7.26667        -56.7918        25.765
7.28889        -59.1003        25.8861
7.31111        -61.3185        26.0756
7.33333        -63.4197        26.2757
7.35556        -65.424        26.4225
7.37778        -67.3905        26.4802
7.4        -69.3881        26.4596
7.42222        -71.4634        26.4128
7.44444        -73.6206        26.411
7.46667        -75.8167        26.5147
7.48889        -77.9756        26.7447
7.51111        -80.0208        27.0652
7.53333        -81.9159        27.3941
7.55556        -83.6883        27.6428
7.57778        -85.4171        27.7608
7.6        -87.1945        27.7585
7.62222        -89.087        27.6996
7.64444        -91.11        27.6751
7.66667        -93.221        27.7702
7.68889        -95.3302        28.0295
7.71111        -97.3334        28.4314
7.73333        -99.1595        28.8884
7.75556        -100.81        29.2841
7.77778        -102.36        29.5287
7.8        -103.924        29.5984
7.82222        -105.603        29.5405
7.84444        -107.455        29.4511
7.86667        -109.473        29.4426
7.88889        -111.585        29.6077
7.91111        -113.674        29.9833
7.93333        -115.615        30.5264
7.95556        -117.333        31.1213
7.97778        -118.84        31.6275
8        -120.233        31.9422
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
均方差(RMSE):26.1269048149089
残差平方和(SSE):247106.686185046
相关系数(R): 0.98500778531148
相关系数之平方(R^2): 0.970240337124227
确定系数(DC): 0.970084420040968
F统计(F-Statistic): 948.644795216808

参数                  最佳估算
--------------------        -------------
epsmax        40.3456192867918
w_pe1        -3732.89238502787
r_e1        7.83476811178322
derta        118.992611095546
w_pe2        -202.103223597636
w_oe        -328.288770160549
r_e2        23.0734958438302
1stopt程序运行
c154.jpg


1stopt程序运行-1
c155.jpg

» 本帖已获得的红花(最新10朵)

2楼2015-05-06 21:50:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clownxxl

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by dingd at 2015-05-06 21:50:33
均方差(RMSE):26.1269048149089
残差平方和(SSE):247106.686185046
相关系数(R): 0.98500778531148
相关系数之平方(R^2): 0.970240337124227
确定系数(DC): 0.970084420040968
F统计(F-Statistic): 948.644795 ...

谢谢了,可不可以再帮我拟合一组参数,下面是程序的代码以及数据:
Variable v,mu_re,mu_im;
SharedModel;
Function mu_re=mumax+w_pm^2*(w_om^2-(600*pi/v)^2)/((w_om^2-(600*pi/v)^2)^2+(r_m*(600*pi/v))^2);
mu_im=w_pm^2*r_m*(600*pi/v)/((w_om^2-(600*pi/v)^2)^2+(r_m*(600*pi/v))^2);   
Data;
3        4.57288        0.0759301
3.02222        4.49991        0.0548232
3.04444        4.44365        0.0393381
3.06667        4.37966        0.0360459
3.08889        4.31531        0.0194135
3.11111        4.27023        0.00900236
3.13333        4.21228        0.00883828
3.15556        4.16094        -0.00867656
3.17778        4.12197        -0.00979982
3.2        4.06726        -0.0149305
3.22222        4.03075        -0.0280318
3.24444        3.98877        -0.0221584
3.26667        3.94411        -0.0367971
3.28889        3.91603        -0.0339028
3.31111        3.86734        -0.0397817
3.33333        3.84388        -0.0458781
3.35556        3.79876        -0.0421289
3.37778        3.77256        -0.0547031
3.4        3.73445        -0.045073
3.42222        3.7042        -0.0596571
3.44444        3.67287        -0.048591
3.46667        3.63869        -0.0630235
3.48889        3.61188        -0.0506645
3.51111        3.57783        -0.0645219
3.53333        3.55331        -0.0532245
3.55556        3.5193        -0.0671069
3.57778        3.49522        -0.055298
3.6        3.46309        -0.0704342
3.62222        3.43657        -0.0570564
3.64444        3.40912        -0.07273
3.66667        3.37805        -0.0582038
3.68889        3.35933        -0.0722893
3.71111        3.32241        -0.062504
3.73333        3.30986        -0.070685
3.75556        3.26876        -0.0684952
3.77778        3.25955        -0.0654166
3.8        3.22215        -0.0753562
3.82222        3.20608        -0.0617597
3.84444        3.18221        -0.0784924
3.86667        3.15253        -0.0648104
3.88889        3.14324        -0.0751943
3.91111        3.10405        -0.0747662
3.93333        3.09655        -0.0667566
3.95556        3.06784        -0.082823
3.97778        3.04321        -0.0669068
4        3.03535        -0.0789418
4.02222        2.99671        -0.0791421
4.04444        2.98954        -0.0679622
4.06667        2.9672        -0.0868192
4.08889        2.93557        -0.0735485
4.11111        2.93251        -0.0760179
4.13333        2.89746        -0.087218
4.15556        2.87853        -0.0680756
4.17778        2.87357        -0.0813029
4.2        2.83398        -0.0846037
4.22222        2.82307        -0.0669745
4.24444        2.81219        -0.0853345
4.26667        2.77129        -0.0816589
4.28889        2.7647        -0.0646179
4.31111        2.75209        -0.0841374
4.33333        2.71184        -0.0786363
4.35556        2.70618        -0.062124
4.37778        2.69385        -0.0819788
4.4        2.65273        -0.07776
4.42222        2.64493        -0.0593346
4.44444        2.63586        -0.0775863
4.46667        2.59375        -0.077395
4.48889        2.58198        -0.0552388
4.51111        2.58026        -0.071177
4.53333        2.53589        -0.0817569
4.55556        2.511        -0.0535364
4.57778        2.51916        -0.0539617
4.6        2.48909        -0.0794501
4.62222        2.44715        -0.06233
4.64444        2.4468        -0.0439814
4.66667        2.43716        -0.0676352
4.68889        2.38924        -0.0689909
4.71111        2.37088        -0.0394767
4.73333        2.37967        -0.0444589
4.75556        2.34619        -0.0702937
4.77778        2.30279        -0.0513386
4.8        2.30309        -0.0282768
4.82222        2.30332        -0.0527525
4.84444        2.2541        -0.0672764
4.86667        2.2217        -0.0375104
4.88889        2.23263        -0.0255576
4.91111        2.2202        -0.060819
4.93333        2.16412        -0.0636594
4.95556        2.14157        -0.0291256
4.97778        2.1598        -0.0276287
5        2.1369        -0.0721174
5.02222        2.07458        -0.0681137
5.04444        2.05704        -0.0291716
5.06667        2.08442        -0.0335868
5.08889        2.05743        -0.092769
5.11111        1.98542        -0.0906444
5.13333        1.96605        -0.0500232
5.15556        1.99976        -0.0603113
5.17778        1.96473        -0.139202
5.2        1.87577        -0.141331
5.22222        1.84582        -0.101056
5.24444        1.87235        -0.115388
5.26667        1.82127        -0.201939
5.28889        1.71499        -0.205001
5.31111        1.66335        -0.166003
5.33333        1.65896        -0.167384
5.35556        1.60316        -0.22491
5.37778        1.49711        -0.22708
5.4        1.42376        -0.187157
5.42222        1.38963        -0.15937
5.44444        1.34919        -0.17286
5.46667        1.26506        -0.179238
5.48889        1.18082        -0.148305
5.51111        1.12207        -0.106022
5.53333        1.07997        -0.0797836
5.55556        1.02455        -0.0718895
5.57778        0.949235        -0.0528813
5.6        0.877152        -0.0145843
5.62222        0.819476        0.0288422
5.64444        0.77021        0.0632198
5.66667        0.713903        0.0883375
5.68889        0.646683        0.118401
5.71111        0.579067        0.159282
5.73333        0.518345        0.205236
5.75556        0.463561        0.248368
5.77778        0.407634        0.285506
5.8        0.344983        0.321443
5.82222        0.277448        0.36189
5.84444        0.209405        0.407499
5.86667        0.143219        0.455555
5.88889        0.0780535        0.502376
5.91111        0.00972264        0.546307
5.93333        -0.0657924        0.590157
5.95556        -0.14865        0.638487
5.97778        -0.235954        0.694001
6        -0.324088        0.75611
6.02222        -0.411879        0.821064
6.04444        -0.503039        0.885868
6.06667        -0.602685        0.952794
6.08889        -0.712        1.02757
6.11111        -0.827519        1.1147
6.13333        -0.943819        1.21422
6.15556        -1.05837        1.3214
6.17778        -1.17511        1.4318
6.2        -1.30073        1.5473
6.22222        -1.43792        1.675
6.24444        -1.58336        1.82173
6.26667        -1.73037        1.98952
6.28889        -1.87444        2.17458
6.31111        -2.0179        2.37197
6.33333        -2.16781        2.58267
6.35556        -2.32864        2.81498
6.37778        -2.49783        3.07952
6.4        -2.66686        3.38344
6.42222        -2.82608        3.72746
6.44444        -2.97058        4.10797
6.46667        -3.10178        4.52395
6.48889        -3.22185        4.98236
6.51111        -3.326        5.49635
6.53333        -3.3992        6.0786
6.55556        -3.41833        6.7348
6.57778        -3.35794        7.46041
6.6        -3.19614        8.24276
6.62222        -2.91578        9.06665
6.64444        -2.49946        9.91709
6.66667        -1.92439        10.7744
6.68889        -1.16439        11.6044
6.71111        -0.201276        12.3519
6.73333        0.958514        12.9441
6.75556        2.2711        13.3068
6.77778        3.65689        13.3896
6.8        5.02037        13.184
6.82222        6.27685        12.7234
6.84444        7.37026        12.0662
6.86667        8.27506        11.2762
6.88889        8.98781        10.4095
6.91111        9.51693        9.5116
6.93333        9.87665        8.61962
6.95556        10.0862        7.76529
6.97778        10.1712        6.97413
7        10.1626        6.26177
7.02222        10.0917        5.63158
7.04444        9.98403        5.07625
7.06667        9.85602        4.58269
7.08889        9.71454        4.13771
7.11111        9.55982        3.73263
7.13333        9.39039        3.36516
7.15556        9.20786        3.03751
7.17778        9.01873        2.7519
7.2        8.83223        2.5065
7.22222        8.65638        2.29458
7.24444        8.4949        2.10652
7.26667        8.34632        1.93298
7.28889        8.20512        1.76803
7.31111        8.06468        1.61086
7.33333        7.92093        1.46534
7.35556        7.7748        1.33679
7.37778        7.6315        1.22808
7.4        7.49741        1.13754
7.42222        7.37685        1.05972
7.44444        7.27027        0.987574
7.46667        7.17419        0.915115
7.48889        7.08247        0.839608
7.51111        6.98891        0.762736
7.53333        6.89023        0.68975
7.55556        6.78782        0.626501
7.57778        6.68673        0.576137
7.6        6.5928        0.537643
7.62222        6.50996        0.506629
7.64444        6.43885        0.477274
7.66667        6.37673        0.444411
7.68889        6.31837        0.405348
7.71111        6.25791        0.361075
7.73333        6.19142        0.316083
7.75556        6.11892        0.276277
7.77778        6.04423        0.245974
7.8        5.97287        0.225964
7.82222        5.90953        0.213578
7.84444        5.85638        0.204083
7.86667        5.81268        0.192387
7.88889        5.77499        0.174599
7.91111        5.73814        0.149361
7.93333        5.69699        0.118644
7.95556        5.64862        0.0872583
7.97778        5.59378        0.0607795
8        5.53649        0.0429587
3楼2015-05-07 08:47:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

Root of Mean Square Error (RMSE):0.387372592569377
Sum of Square Error:67.826001514212
Correlation Coef. (R): 0.993131339824686
R-Square: 0.986309858141975
Determination Coef. (DC): 0.985250159765273
F-Statistic: 3138.76454601051

Parameter                  Best Estimate
--------------------        -------------
mumax        3.16728622593316
w_pm        -255.229864047357
w_om        278.719272152532
r_m        17.2556124790563
4楼2015-05-07 09:52:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 clownxxl 的主题更新
信息提示
请填处理意见