当前位置: 首页 > 计算模拟 >高版本1stopt拟合求助,万分感谢

高版本1stopt拟合求助,万分感谢

作者 虾米一只111
来源: 小木虫 250 5 举报帖子
+关注

Title "Z1";
ComplexStr = i;
Parameters Z0,m,t,c;
Variable f,Z1;
Function Z1=Z0*(1-m*(1-1/(1+(i*6.28*t)^c)));
date;
//f,Z1;
1        443756
1.063        429749
1.131        419324
1.202        410917
1.278        403661
1.359        397516
1.445        392408
1.536        387741
1.633        383226
1.737        376250
1.847        361775
1.963        358622
2.088        355508
2.22        352868
2.36        350175
2.509        347738
2.668        345227
2.837        342783
3.016        340612
3.207        338454
3.41        336397
3.626        334375
3.855        332334
4.099        330332
4.358        328422
4.634        326564
4.927        324825
5.239        323205
5.57        321568
5.923        319993
6.297        318219
6.696        316639
7.119        314986
7.569        313276
8.048        311859
8.557        310149
9.099        308495
9.674        306832
10.286        305369
10.937        303939
11.629        302441
12.364        300800
13.147        299535
13.978        298043
14.862        296489
15.802        295235
16.802        293748
17.865        292449
18.995        291079
20.197        289796
21.474        288547
22.833        287272
24.277        286248
25.813        284624
27.445        283488
29.182        281835
31.028        280981
32.99        279672
35.077        279557
37.296        277560
39.655        276725
42.164        275719
44.831        272656
47.667        266930
50.682        278738
53.888        271849
57.297        269797
60.921        266744
64.775        267350
68.872        266589
73.229        264085
77.861        265795
82.787        263371
88.024        262529
93.592        261369
99.512        260271
105.81        259093
112.5        258666
119.62        256814
127.18        257089
135.23        255139
143.78        254975
152.88        252915
162.55        252496
172.83        251582
183.76        250884
195.39        249707
207.75        248618
220.89        247906
234.86        246220
249.72        246206
265.52        244551
282.31        243170
300.17        242259
319.16        241559
339.35        240035
360.81        239375
383.64        238206
407.9        237374
433.71        236602
461.14        235460
490.31        234354
521.33        233147
554.31        231975
589.37        230975
626.65        229600
666.29        228638
708.44        227104
753.25        226191
800.9        224234
851.56        222483
905.43        221009
962.71        220065
1023.6        218528
1088.4        216901
1157.2        215441
1230.4        213857
1308.2        211998
1391        210249
1479        208709
1572.5        206718
1672        205051
1777.8        203223
1890.2        201217
2009.8        199413
2136.9        197447
2272.1        195517
2415.8        193347
2568.7        191256
2731.2        189369
2903.9        187234
3087.6        185282
3282.9        183018
3490.6        180902
3711.4        178725
3946.2        176663
4195.8        174523
4461.2        172469
4743.4        170445
5043.5        168390
5362.5        166381
5701.7        164353
6062.4        162416
6445.9        160419
6853.6        158492
7287.2        156555
7748.1        154680
8238.3        152775
8759.4        150965
9313.5        149092
9902.6        147333
10529        145516
11195        143812
11903        142135
12656        140317
13457        138607
14308        136974
15213        135245
16176        133629
17199        131896
18287        130301
19443        128789
20673        127023
21981        125490
23372        123808
24850        122189
26422        120667
28093        119018
29870        117313
31760        115980
33769        114113
35905        112620
38176        110979
40591        109205
43159        107277
45889        105884
48792        104128
51878        102391
55160        100793
58649        98973.9
62359        97412.7
66304        95785.3
70498        93933.4
74958        92271.8
79699        90550.8
84741        88700.3
90101        87115.7
95801        85353
101860        83534.2
108300        81915.5
115160        80110.5
122440        78560.4
130190        76937.8
138420        74759.4
147180        73019.2
156490        71283.6
166390        69591.8
176910        67770.5
188100        65973.3
200000        64205.6@月只蓝 返回小木虫查看更多

今日热帖
  • 精华评论
  • dingd

    三个问题:
    1:Z1是实部值还是虚部值?
    2:参数Z0,m,t,c都是实数?或是复数?
    3:公式“ Z1=Z0*(1-m*(1-1/(1+(i*6.28*t)^c)));”根本没出现自变量f,正确吗?

  • 虾米一只111

    引用回帖:
    2楼: Originally posted by dingd at 2017-04-10 09:18:32
    三个问题:
    1:Z1是实部值还是虚部值?
    2:参数Z0,m,t,c都是实数?或是复数?
    3:公式“ Z1=Z0*(1-m*(1-1/(1+(i*6.28*t)^c)));”根本没出现自变量f,正确吗?

    Z1是复数,测出来的Z1是模值,参数都是实数,大神在用下面的试一下吧,为什么我用1stopt1.5就算不了下面的呢
    Title "Type your title here";
    Parameters Z0,m[0,1],t,c[0,1];
    Variable f,Z1;
    Function Z1=Z0*((1+2*(1-m)*((6.28*f*t)^c)*cos(3.14*c/2)+((1-m)^2)*(6.28*f*t)^(2*c))/(1+(2*(6.28*f*t)^c)*cos(3.14*c/2)+(6.28*f*t)^(2*c)))^(1/2);
    date;
    //f,Z1;
    1        2672280
    1.063        2637580
    1.131        2603820
    1.202        2572840
    1.278        2541090
    1.359        2507950
    1.445        2478280
    1.536        2452960
    1.633        2429570
    1.737        2406790
    1.847        2384150
    1.963        2364720
    2.088        2344850
    2.22        2326000
    2.36        2306890
    2.509        2289380
    2.668        2268680
    2.837        2252920
    3.016        2236680
    3.207        2222060
    3.41        2205420
    3.626        2189910
    3.855        2173360
    4.099        2159710
    4.358        2144220
    4.634        2131120
    4.927        2117340
    5.239        2104890
    5.57        2091420
    5.923        2078460
    6.297        2063730
    6.696        2051670
    7.119        2039140
    7.569        2025270
    8.048        2007670
    8.557        1996740
    9.099        1987110
    9.674        1974880
    10.286        1961470
    10.937        1948590
    11.629        1934940
    12.364        1924080
    13.147        1909720
    13.978        1892840
    14.862        1879130
    15.802        1868170
    16.802        1855690
    17.865        1845950
    18.995        1830630
    20.197        1812600
    21.474        1807850
    22.833        1803980
    24.277        1793930
    25.813        1789510
    27.445        1755860
    29.182        1757580
    31.028        1705720
    32.99        1724630
    35.077        1684100
    37.296        1680680
    39.655        1715190
    42.164        1641390
    44.831        1739690
    47.667        1646190
    50.682        1692550
    53.888        1793590
    57.297        1617400
    60.921        1506230
    64.775        1502970
    68.872        1558190
    73.229        1463570
    77.861        1561710
    82.787        1496610
    88.024        1458560
    93.592        1437150
    99.512        1423670
    105.81        1396740
    112.5        1373040
    119.62        1364220
    127.18        1339260
    135.23        1335800
    143.78        1276450
    152.88        1290290
    162.55        1240340
    172.83        1213290
    183.76        1179790
    195.39        1163060
    207.75        1155230
    220.89        1100840
    234.86        1096520
    249.72        1063860
    265.52        1055530
    282.31        1012090
    300.17        970230
    319.16        959770
    339.35        935819
    360.81        913507
    383.64        884057
    407.9        860688
    433.71        828063
    461.14        813422
    490.31        794543
    521.33        769208
    554.31        740689
    589.37        721380
    626.65        700036
    666.29        685124
    708.44        661148
    753.25        647662
    800.9        627651
    851.56        607534
    905.43        591570
    962.71        579160
    1023.6        562435
    1088.4        547797
    1157.2        532802
    1230.4        516786
    1308.2        508478
    1391        493107
    1479        483909
    1572.5        474241
    1672        462584
    1777.8        452253
    1890.2        443462
    2009.8        435123
    2136.9        426321
    2272.1        418628
    2415.8        410716
    2568.7        403914
    2731.2        396138
    2903.9        390006
    3087.6        383816
    3282.9        377615
    3490.6        371902
    3711.4        365772
    3946.2        360513
    4195.8        355408
    4461.2        350125
    4743.4        345059
    5043.5        340532
    5362.5        335705
    5701.7        331151
    6062.4        326533
    6445.9        322325
    6853.6        317968
    7287.2        313705
    7748.1        309473
    8238.3        305276
    8759.4        301016
    9313.5        297117
    9902.6        292686
    10529        288990
    11195        284157
    11903        280747
    12656        276693
    13457        272381
    14308        268361
    15213        264550
    16176        260134
    17199        256043
    18287        251799
    19443        247506
    20673        243872
    21981        238974
    23372        234782
    24850        230811
    26422        226763
    28093        222279
    29870        218160
    31760        213460
    33769        210108
    35905        204966
    38176        198866
    40591        196629
    43159        192043
    45889        188583
    48792        183884
    51878        179425
    55160        175130
    58649        171387
    62359        166847
    66304        162712
    70498        158643
    74958        154292
    79699        150412
    84741        145853
    90101        141803
    95801        137847
    101860        133918
    108300        130054
    115160        126031
    122440        122288
    130190        118386
    138420        114729
    147180        110928
    156490        107557
    166390        103713
    176910        100136
    188100        96895.8
    200000        93390.1

  • dingd

    引用回帖:
    3楼: Originally posted by 虾米一只111 at 2017-04-14 12:00:32
    Z1是复数,测出来的Z1是模值,参数都是实数,大神在用下面的试一下吧,为什么我用1stopt1.5就算不了下面的呢
    Title "Type your title here";
    Parameters Z0,m,t,c;
    Variable f,Z1;
    Function Z1=Z0*( ...

    公式里虚数符号都没有了,不就变成一个普通的实数拟合问题了:

    均方差(RMSE): 62837.5235514104
    残差平方和(SSE): 789710873214.811
    相关系数(R): 0.99697395680469
    相关系数之平方(R^2): 0.993957070546801
    修正R平方(Adj. R^2): 0.993895721009205
    确定系数(DC): 0.993956797613149
    卡方系数(Chi-Square): 344523.815157506
    F统计(F-Statistic): 10736.7575976312

    参数        最佳估算
    ----------        -------------
    z0        2642790.06747986
    m        0.972182633784368
    t        0.00223886517828685
    c        0.492671333535272
    高版本1stopt拟合求助,万分感谢
    c35.jpg

  • lipenggg

    是data;
    不是date;

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓