高版本1stopt拟合求助,万分感谢
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@月只蓝 返回小木虫查看更多
三个问题:
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
,
公式里虚数符号都没有了,不就变成一个普通的实数拟合问题了:
均方差(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
c35.jpg
是data;
不是date;