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

当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

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

  1. 用MATLAB做结构灵敏度分析怎么做,如何做灵敏度分析 - 数学 - 小木虫 - 学术 科研 互动社区...

    我有以下数据 年份        平均气温        最高气温        最低气温        相对湿度        日照时数        实际风速        日均ET 1953   ...

  2. matlab换热器优化设计,matlab用于求解一个简单的优化设计问题 - 数学 - 小木虫 - 学术 科研 互动社区...

    最近在学习工程优化设计与MATLAB实现,在模拟课本中的例子的时候,程序提示不正确. Error using ==> mrdivide Matrix dimensions must agree. ...

  3. matlab如何进行四维拟合,如何用matlab对四维数据进行非线性拟合？ - 数学 - 小木虫 - 学术 科研 互动社区...

    这个很简单啊,用lsqcurvefit或者lsqnonlin都可以.比如: function ydata = myfun(a,xdata) ydata=xdata(:,1).*xdata(:,2).* ...

  4. matlab imfill holes,到底怎么实现，bwfill和imfill - 程序语言 - 小木虫 - 学术 科研 互动社区...

    认真研读一下MATLAB的help文档吧,解释最权威: BWFILL Fill background regions in binary image. BWFILL is a grandfathere ...

  5. matlab中lower,【求助】matlab,这个错误究竟是什么？ - 数学 - 小木虫 - 学术 科研 互动社区...

    谢谢以上两位,程序改动了一下,以前的问题暂时没有出现,但是出现了下边的问题: 主程序中是这样写的,调用curvefunzscanAbs作为拟合方程,其中NonAbs是拟合参数,transIntUp,p ...

  6. matlab文件序号超出511,求教一段matlab的代码 - 数学 - 小木虫 - 学术 科研 互动社区...

    不知楼主说的没法求出优化值,是什么意思.是报错还是怎么. 我也是刚接触优化,觉得楼主的错误可能是在:主函数里 A=[-1;1];的括号中间应该用逗号,而不是分号.即改为A=[-1,1]; 我对你的程序 ...

  7. matlab 打开软件报错,matlab程序错误，提示如下【弄了半夜还是不行】 - 数学 - 小木虫 - 学术 科研 互动社区...

    function [x,y,Y]=levydistrintheorytest(t0,lambda,T,a,b); %  计算理论的二维levy分布的图象   画出来的3维曲面图 %   t0 :lev ...

  8. matlab级数求和阶乘,(阶乘)数列求和 - 数学 - 小木虫 - 学术 科研 互动社区

    翻看了一下坛子里的相关内容,觉得应该是发在这里最合适了.要是这里不合适,请大家告知我一下. 我在对一化学过程数学建模时,得到一数学关系,是一数列,其中各项表达如下: [latex]{x_1} = \f ...

  9. matlab读取超大txt文件,fortran如何读取超大文本文件 - 程序语言 - 小木虫 - 学术 科研 互动社区...

    有这样一个外部力的文件,需要每个时间步调用 Fbrk.txt 文件内容如下,第一列是时间 1.07800            13.73559           -16.39831        ...

最新文章

  1. YOLOv2/YOLO9000 《YOLO9000: Better, Faster, Stronger》论文笔记
  2. 第1节 yarn：15、关于yarn中常用的参数设置
  3. 括号生成Python解法
  4. 计算机专业网站的开题ppt,.计算机专业开题报告.ppt
  5. OpenCV_(Fit Line with points)用直线拟合一组点
  6. 为什么这么多的广告呢
  7. Delta3D（7）教程：在场景编辑器中加载角色
  8. USTC2017 writeup
  9. 自定义jsx解析工厂函数
  10. 史玉柱加持，Playtika上市能否圆昔日“巨人”翻身梦？
  11. 力扣刷题 DAY_62 回溯
  12. transition天坑
  13. 克鲁斯卡尔（kruskal）算法代码释义
  14. Ubuntu 22.04 LTS (Jammy Jellyfish) Daily Build安装镜像PC ARM Raspberry Pi
  15. 选择背光需要对比哪些因素呢？
  16. P1903 数颜色 (带修莫队)
  17. 哪个软件能代替斐讯路由_斐讯路由器控制软件 4.2.1 去广告版 - 用手机管理斐讯路由器...
  18. 【C语言】题目：输入某年某月某日，判断这一天是这一年的第几天？
  19. 信息技术应用 中职计算机教学,新形势下信息技术在中职计算机课程教学中的应用...
  20. 第四章 使用 matplotlib 绘制散点图

热门文章

  1. 笔记本有线网络共享为WiFi
  2. 医德：《大医精诚》和《希波克拉底誓言》
  3. PTA-链表学习-----奇数值结点链表
  4. 树莓派安装 Docker CE
  5. PYTHON地理出图配色及旁门左道
  6. 机械秒表的使用方法_秒表的使用方法！
  7. C++图片格式转换：BMP转JPEG
  8. 高速缓存对c程序的影响
  9. 自学c语言困难吗,你觉得这4块是C语言最难学的部分吗？
  10. safari浏览器找不到服务器