原创力文档- 1、本文档共15页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
11.2拟合
测量数据的数据量较大或者测量值与真实值有误差，这时一般用曲线拟合的方法来解决问题.已知离散点上的数据集，即已知在点集上的函数值，构造一个解析函数f(x)（其图形为一曲线），使f(x)在原离散点xi上的值尽可能接近给定的yi值，这一构造函数f(x)的过程称为曲线拟合.最常用的曲线拟合方法是最小二乘法，该方法是寻找函数f(x)使得最小.在MATLAB中，利用polyfit和lsqcurvefit函数进行拟合.
11.2.1多项式拟合polyfit函数调用格式如下：p=polyfit(x,y,n)其中x和y是节点数据向量，n是多项式的次数；p是返回的多项式p(x)=p1xn+p2xn?1+...+pnx+pn+1系数向量，系数按降幂排列，p的长度为n+1.
例11.8求如下给定数据的拟合曲线：x=[0.5,1.0,1.5,2.0,2.5,3.0]，y=[1.75,2.45,3.81,4.80,7.00,8.60].解MATLAB命令如下：x=[0.5,1.0,1.5,2.0,2.5,3.0];y=[1.75,2.45,3.81,4.80,7.00,8.60];p=polyfit(x,y,2)%选取次数的时候，次数长度不得超过length(x)-1x1=0.5:0.05:3.0;y1=polyval(p,x1);%求以p为系数的多项式的值plot(x,y,*r,x1,y1,-b)计算结果为p=0.56140.82871.1560此结果表示拟合函数为：用此函数拟合数据的效果如图11-7所示.图11-7多项式拟合
例11.9某乡镇企业2015-2021年的生产利润如表11-3.?试预测2022年的利润.?解作已知数据的散点图，x0=[2014201520162017201820192020];y0=[70122144152174196202];plot(x0,y0,*)发现该乡镇企业的年生产利润几乎直线上升（见图11-8）.因此，我们可以y=a1x+a0用作为拟合函数来预测该乡镇企业未来的利润.年份2015201620172018201920202021利润（万元）70122144152174196202表11-3生产利润
?图11-8数据的散点图
11.2.2非线性拟合lsqcurvefit函数是MATLAB提供的用于求解非线性最小二乘拟合问题的专用函数，用作各种类型曲线的拟合，调用格式如下：x=lsqcurvefit(fun,x0,xdata,ydata)其中，fun是非线性函数模型；x0是求解初始值；xdata和ydata是原始输入和输出数据向量；返回值x为拟合输出的未知参数.
例11.10用最小二乘法求形如的函数拟合下列数据.x0.10.20.150-0.20.3y0.950.840.861.061.500.72表11-4拟合数据解先将参数，合写为，编写如下程序：fun=@(c,x)c(1)*exp(c(2)*x)x=[0.1,0.2,0.15,0,-0.2,0.3];y=[0.95,0.84,0.86,1.06,1.50,0.72];c=lsqcurvefit(fun,[0,0],x,y)norm(feval(fun,c,x)-y)^2
其中，[0,0]是初始值，最后一句是计算残差的平方和，也就是拟合函数在给定点的值和原始数据的差的平方和，运行结果为：Optimizationterminated:relativefunctionvaluechangingbylessthanOPTIONS.TolFun.c=1.0997-1.4923ans=0.0031可得a=1.0997，b=-1.4923，拟合函数为.ans=0.0031说明残差很小.
例11.11我国从1954年到2005年全国总人口的数据如表11-5，请预测2019年和2030年的全国人口数.表11-5各年份全国总人数（单位：千万）
解因
文档评论（0）