大家好!我想对两组数据(数据和图见下)同时进行拟合,得到一组参数。不知哪位高手能编写出程序,最后生成:原始数据,拟合数据,参数值,每个参数的误差以及总方差。
当x较小时,y变化较快,需要对x求一下对数。
待求参量:
#A, #B, #C, #D, #E, #F, #G (y1含有前6个参量,y2含有7个参量)
对参量的约束条件:
1<#A<30;
20<#B<200,#B>#A;
10^(-10)<#C<10^(-7);
0<=#D<=1;
10^(-12)<#E<10^(-5);
0<=#F<=1;
10^(-12)<#G<10^(-5)
拟合函数:
y1=#A+(#B-#A)*(1+(2*3.1416*x*#C)^(1-#D)*sin(#D*3.1416/2))/(1+2*(2*3.1416*x*#C)^(1-#D)*sin(#D*3.1416/2)+(2*3.1416*x*#C)^(2-2*#D))+(2*3.1416*x)^(#F-1)/(8.854*10^-12)*#E*tan(3.1416/2*#F)
y2=(#B-#A)*(2*3.1416*x*#C)^(1-#D)*cos(#D*3.1416/2)/(1+2*(2*3.1416*x*#C)^(1-#D)*sin(#D*3.1416/2)+(2*3.1416*x*#C)^(2-2*#D))+#G/(2*3.1416*x*8.854*10^-12)+(2*3.1416*x)^(#F-1)*#E/(8.854*10^-12)
观测值:
x y1 y2
40 130.9108 88.32217
42.799 111.5734 79.52563
45.794 135.18311 67.34402
48.999 123.1284 84.34343
52.428 130.7267 69.48528
56.097 112.3516 78.91016
60.022 121.4998 65.06214
64.222 98.48351 73.76794
68.717 102.7466 67.00543
73.525 94.70784 64.25732
78.671 110.2479 67.14107
84.176 96.17056 63.70906
90.067 89.13078 64.87036
96.369 83.95633 56.52533
103.113 83.2327 56.5466
110.329 80.07961 53.27859
118.05 80.67174 55.32287
126.311 77.64321 55.53664
135.14999 75.29877 51.31252
144.60699 71.00301 47.85907
154.72701 69.6681 48.15019
165.55499 68.74583 47.59437
177.14 67.43265 45.89816
189.536 66.16602 42.72585
202.8 61.29647 42.18708
216.991 60.87629 41.72636
232.17599 60.04741 41.9826
248.424 55.9209 39.09031
265.80801 55.70895 38.088
284.409 54.11443 37.93357
304.31201 52.39143 35.2491
325.608 51.25287 34.43814
348.39301 50.44999 33.30919
372.77399 48.69473 31.95148
398.85999 47.92783 31.00966
426.772 46.7323 29.805
456.63699 45.95324 29.28115
488.59201 44.60682 28.21328
522.78302 43.6366 27.10598
559.367 42.65884 25.81369
598.51099 41.59007 25.16921
640.39398 40.77101 24.17984
685.20898 39.69152 23.49695
733.159 38.98373 22.45583
784.46503 38.16704 21.86946
839.36102 37.49696 20.99112
898.099 36.84111 20.10395
960.94702 36.07808 19.54207
1028.19299 35.38936 18.71043
1100.14502 34.77517 17.99798
1177.13196 34.15973 17.39159
1259.50696 33.59982 16.71645
1347.646 33.05671 16.13231
1441.953 32.58645 15.63731
1542.85999 32.01151 14.94919
1650.828 31.51963 14.4221
1766.35205 30.96424 13.92618
1889.95898 30.60521 13.38352
2022.21704 30.12229 12.95935
2163.72998 29.70016 12.4626
2315.146 29.32128 12.0107
2477.15796 28.91598 11.56298
2650.50708 28.53142 11.16361
2835.98706 28.17922 10.75587
3034.44702 27.83117 10.37639
3246.79492 27.50975 10.01584
3474.00293 27.17833 9.65472
3717.11108 26.86233 9.3266
3977.23096 26.54549 8.96851
4255.5542 26.27776 8.67039
4553.354 25.99566 8.35968
4871.99414 25.72867 8.06178
5212.93213 25.45255 7.77694
5577.72803 25.2096 7.5081
5968.05322 24.95931 7.25553
6385.69189 24.72685 6.99624
6832.55713 24.50187 6.75443
7310.69385 24.28611 6.51664
7822.28906 24.07368 6.28817
8369.68555 23.86149 6.08031
8955.38965 23.66551 5.87051
9582.08008 23.46948 5.67001
10252.62598 23.2845 5.47798
10970.0957 23.10101 5.28829
11737.77344 22.93066 5.10582
12559.17285 22.7633 4.93827
13438.05273 22.59646 4.76626
14378.43555 22.44618 4.60521
15384.62598 22.2944 4.45443
16461.22852 22.14483 4.30725
17613.17188 22.00053 4.16272
18845.72656 21.86432 4.0239
20164.5332 21.73107 3.88866
21575.63086 21.6034 3.75837
23085.47266 21.47595 3.63318
24700.97461 21.35446 3.51537
26429.52734 21.23675 3.39919
28279.04297 21.12457 3.28676
30257.98438 21.01587 3.1775
32375.41211 20.90948 3.07392
34641.01563 20.80532 2.97374
37065.16406 20.70618 2.88081
39658.95313 20.60925 2.7862
42434.25 20.51352 2.69559
45403.76172 20.42471 2.60942
48581.07813 20.33565 2.52619
51980.73828 20.24922 2.44531
55618.30859 20.16652 2.36399
59510.42969 20.08609 2.28895
63674.91406 20.00891 2.21597
68130.82813 19.93334 2.14577
72898.57031 19.85959 2.07717
77999.94531 19.78893 2.01078
83458.3125 19.72021 1.94762
89298.65625 19.65268 1.88451
95547.69531 19.58786 1.82502
102234.0391 19.52355 1.76602
109388.2891 19.4645 1.71045
117043.1875 19.40484 1.65523
125233.7656 19.34647 1.60227
133997.5156 19.29044 1.54989
143374.5469 19.23703 1.50043
153407.7656 19.18514 1.45196
164143.1094 19.13492 1.39906
175629.7031 19.08501 1.35566
187920.125 19.03795 1.31159
201070.6094 18.99183 1.2743
215141.3594 18.94748 1.23305
230196.75 18.90616 1.19191
246305.7188 18.86405 1.15433
263542 18.82389 1.11703
281984.4375 18.78548 1.07926
301717.4375 18.74832 1.04396
322831.375 18.71346 1.01068
345422.8438 18.67629 0.97772
369595.2188 18.64261 0.94535
395459.1563 18.61007 0.91493
423133.0625 18.57844 0.88441
452743.5313 18.55013 0.85544
484426.125 18.52 0.82496
518325.8438 18.49192 0.79865
554597.875 18.46391 0.77155
593408.125 18.43856 0.74664
634934.3125 18.41295 0.72168
679366.4375 18.38894 0.69756
726907.9375 18.364 0.67463
777776.3125 18.34106 0.65122
832204.4375 18.31869 0.63035
890441.375 18.29826 0.60946
952753.6875 18.27835 0.58924
1.02E+06 18.25747 0.56853
1.09E+06 18.23779 0.54971
1.17E+06 18.21692 0.52943
1.25E+06 18.19978 0.5139
1.34E+06 18.18398 0.49688
1.43E+06 18.16855 0.4808
1.53E+06 18.15184 0.46522
1.64E+06 18.1368 0.4505
1.75E+06 18.12179 0.43518
1.87E+06 18.1055 0.42126
2.00E+06 18.09197 0.40866
2.15E+06 18.07855 0.39695
2.30E+06 18.06528 0.38488
2.46E+06 18.05082 0.3717
2.63E+06 18.03846 0.3607
2.81E+06 18.02646 0.34861
3.01E+06 18.01395 0.33598
3.22E+06 18.00267 0.32527
3.44E+06 17.99192 0.31493
3.69E+06 17.97955 0.30545
3.94E+06 17.96903 0.29605
4.22E+06 17.95798 0.28646
4.51E+06 17.94747 0.27755
4.83E+06 17.93674 0.26857
5.17E+06 17.92351 0.26092
5.53E+06 17.91276 0.25269
5.92E+06 17.90242 0.24409
6.33E+06 17.89152 0.23572
6.77E+06 17.88139 0.22712
7.25E+06 17.87085 0.21797
7.76E+06 17.86081 0.20954
8.30E+06 17.85117 0.20016
8.88E+06 17.84113 0.19098
9.50E+06 17.8315 0.18221
1.02E+07 17.81958 0.16899
1.09E+07 17.81006 0.16002
1.16E+07 17.80062 0.15037
1.25E+07 17.79126 0.14172
1.33E+07 17.78176 0.13266
1.43E+07 17.77221 0.124
1.53E+07 17.76218 0.11459
1.63E+07 17.75181 0.10532
1.75E+07 17.74104 0.09607
1.87E+07 17.73037 0.08655
2.00E+07 17.71898 0.07692
2.14E+07 17.71209 0.07477
2.29E+07 17.69966 0.06699
2.45E+07 17.68639 0.05993
2.62E+07 17.67342 0.05409
2.80E+07 17.65588 0.04905
3.00E+07 17.63785 0.04541
两曲线同时拟合.png