matlab 比较两个函数,Matlab同时拟合两个函数 - 数学 - 小木虫 - 学术 科研 互动社区...

大家好!我想对两组数据(数据和图见下)同时进行拟合,得到一组参数。不知哪位高手能编写出程序,最后生成:原始数据,拟合数据,参数值,每个参数的误差以及总方差。

当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

w120h2494731_1457423759_330.png

两曲线同时拟合.png

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值