最小二乘法拟合的MATLAB和Excel实现

1、最小二乘法拟合的 MATLA盾口 Excel实现摘要:生活生产中我们会遇到各种各样的数据处理，然而这些数据并不像理想实验中得到的数据，有的是一元或多元函数的分布， 有的是一次或多次函数的分布，这就需要我们首先观察数据的散点图，进而选择合理的选择函数进行拟合，同时分析计算该拟合得到的误差，找MATLAEft Excel并对不同的拟合方MATLAE Excel 数据出最优的拟合方式。本文从数学上对最小二乘法原理进行了阐述，并通过 完成数据的拟合，在进行数据拟合中使用的一次函数拟合和多项式拟合， 式进行了比较，到了不同拟合方式下的拟合函数和拟合误差。同时对 拟合方式进行了对比。关键字：最小二乘法MA

2、TLAB Excel 数据拟合Abstract:we will encounter a variety of data processing in production life .However these data is not the data as we expect in ideal experiment;some distribution is a univariate or multivariate functions, some is one or more times function.So we should observe the scatter data chart,a

3、nd then choose the reasonable selection function fitting, make an error analysis and find out the best way of fitting. This paper expound the principle of least square mathematically,complete data fitting by MATLAB and Excel,and use a function fitting and polynomial fitting.we also compare the diffe

4、rent fitting methods,the fitting function and fitting error by the way of MATLAB and Excel.Keywords: Least squares MATLAB Excel Data fitting引言工程试验中我们常常遇到这样的问题，试验中我们会得到各种各样的数据， 不同的数据之间存在着这样那样的关系，如何把得到的试验数据用函数关系式来 得到不同组数据之间的关系，并且在经过数据处理后得到的函数能够客观准确的 描述数据与数据数据之间的关系。 如何选择数据的拟合方式，是线性拟合还是非 线性拟合？是一次还是多次拟合？

5、如何使相关系数 R2接近1?通常我们用最小二 乘法来确定拟合曲线和和该数据的经验公式。本文从数学角度给出最小二乘法的 推导过程，从现实生活中给出曲线拟合的实际应用，同时用 MATLAB日Excel两 种不同的方法对数据进行曲线拟合的实现。当由实验提供了大量数据时，不能要求拟合函数S(x)在数据点(x,y)处的偏差,即& =S(xD 一乂(i=1,2, -,m)严格为零，但为了使近似曲线尽量反映所给数据点的变化趋势，需对偏差有所要求.通常要求偏差平方和mmmZ 62 =，S (Xi) -yi2 =min Z S(Xi)- yi2(1-1 )i Ai =0i zQ这里S(x) =ao%(x

6、)+ai8i(x) + anJx) (n < m)(1-2)这就是一般的最小二乘逼近，用几何语言说，这就称为曲线拟合的最小二乘 法。1用最小二乘法求拟合曲线的问题，就是在形如(1-2 )式的S (x)中求一 函数y = S(x)使偏差的平方和最小。它转化为求多元函数mnI(a°,a1,an)=£ 功(x)£ aj气(为)f (x)2(1-3)i=0j =0的极小点(aQ,a*,a；)的问题。有求多元函数极值的必要条件，有.:I:ak(1-4)(1-5)(1-6)mn=2£ 切(为)£ aj气(x)-f(x)2q(x) =。，k=0,1 ，

7、- -ni =0j =0若记m("k) ='(xi)"xi)k(xi)i =0(f, k) ='(xi)f(x) lx广dk,k =0,1" ni =0则可改写为矩阵形式Ga = d其中 a =(&,涌,an)T,d =(d°,d1,dn)T,& ,平)财，中).(平，中JI000 '10 ' n(平中)(中中)(平中)1，01，1.1， n.仰平)(中中)(甲平)n，0 n，1. n， n2最小二乘法应用x1.22.84.35.46.87.9y2.111.528.141.972.391.4co :i111

8、111y=S (x)例1已知数据试用最小二乘法拟合这组数据的曲线。解：根据所给的数据，在坐标纸上标出，见右 图,从途中看到各点在一条直线附近，故选择线性函数作拟合曲线，即令S(x) =ao +ax，这里 m=5 n=1, 9°(x)=1,%x)=x,故5(*0,0)=' i = 6,i =05(、0, 1)= ("o)=寸，为=284i =055( ""=' 点=165.58,( ", f) 5 =247.3, i =0i=05(%, f)=寸 Bfj =1595.51i =0由法方程得线性方程组 6a0 28.4a1 =247

9、.3' 28.4 165.58% =1595.51解得a°=-23.3498 , a1=13.6408。丁是所求拟合曲线为S1(x) - -23.3498 13.6408 x同理，可以分别求出二次和三次拟合的方程S2(x )和&(x )分别为:25 (x) =1.4424* x2 0.4740* x -0.83286 (x) =0.0899* x3 =2.6690* x2 -4.2774* x 3.81433数据拟合结果的实现和误差的分析在例1中我们运用数学公式进行了最小二乘法的拟合运算2,这种方法在计算过程中太过丁复杂，下面介绍两种方法对上例进行数据拟合，首先用MA

10、TLAB进行编程处理。拟合的程序代码见附录。该实验程序得到的拟合图形如下：10Q80604020100数据散点图疗一朋91 伊 1*x'+2 669*/4277除+3一814100。原始数据1二次曲线£/ - rysz挪-1y=1 442*A 4740*x- 83288060402062S806040200La。原始数据三次曲线图3-1数据的散点图和不同类型的拟合曲线在MATLAB的得到的拟合函数及均方偏差结果如下:yi = 13 .64 * x - 23 .352 一 一 一、1 = 2.734957e+ 002y2_ 一 2_ _ ._ =1 .442 * x 0.474

11、0 * x - 0.8328、2 = 1.645645e+ 0013y3 = 0.0899 * x 2.669 *、32 = 1.363550e+ 0012- _x - 4.2774 * x 3.814从得到的拟合函数和均方偏差可知，用一次拟合得到的一次函数的均方偏差 比较的大，也就是用多项式函数拟合得到的均方偏差更小， 随着多项式的次数越 多，均方偏差趋近丁零。数据散点图曲线拟合（一次）曲线拟合（二次）曲线拟合（三次）图3-2用Excel表格进行数据拟合结果图用Excel表格拟合得到的拟合函数及相关系数结果如下yi = 13 .641 * x - 23 .350R12 = 0.9549 ,_

12、.2y2 = 1 .4424 * x 0.4740 * x - 0 .8328R22 = 0.9973* x 3 .8143y3 = 0.0899 *x3 = 2.6690 *x 2 - 4.2774R； = 0.9978同理，用Excel得到的拟合函数和 MATLAB得到的函数结果相同，在Excel 中进行数据拟合时用相关系数 F2表示误差，（0<R2<1）当相关系数越趋丁 1,说明 得到的试验结果越准确，误差越小。可以看出运用MATLAB Excel进行对数据进行曲线拟合比人工计算大大减 少了计算量，运用软件处理数据已经是当前时代解决数学问题的重要方法之一。参考文献1 李庆阳，

13、王能超，易大义.数值分析（第五版）.北京：清华大学出版社，20082 姜健飞，胡良剑，唐俭.数值分析及其 MATLAB式验.北京：科学出版社，2004学院:学号:完成人签名:研究方向：导师姓名：完成时间：2014年12月5日附录MATLAB程代码如下：x=1.2 2.8 4.3 5.4 6.8 7.9;y=2.1 11.5 28.1 41.9 72.3 91.4;subplot(2,2,1)stem(x,y);xlabel('x')ylabel('y')title('数据散点图')p=polyfit(x,y,1); %拟合出的一次函数的系数yc=

14、y-polyval(p,x); % 计算误差s2=sum(yc.A2); %误差的平方和disp(sprintf('一 次误差的平方和=%d',s2);x1=linspace(min(x),max(x); %绘图用到的点的横坐标y1=polyval(p,x1); % 拟合曲线的纵坐标subplot(2,2,2);plot(x,y,'o',x1,y1); %绘图，原始数据+拟合曲线legend(原始数据','一次曲线',2); %图示s=char(vpa(poly2sym(p,'x'),4); % 一次函数式转换为字符申，v

15、pa转换小数，保 留4位有效数字title('y=' s);p=polyfit(x,y,2); %拟合出的二次函数的系数yc2=y-polyval(p,x); % 计算误差s2=sum(yc2.A2); %误差的平方和disp(sprintf('二次误差的平方和=%d',s2);x2=linspace(min(x),max(x); %绘图用到的点的横坐标y2=polyval(p,x2); %拟合曲线的纵坐标subplot(2,2,3);plot(x,y,'o',x2,y2); %绘图，原始数据+拟合曲线legend(原始数据,'二次曲线',2); %图示s=char(vpa(poly2sym(p,'x'),4); %二次函数式转换为字符申，vpa转换小数，保 留4位有效数字title('y=' s);p=polyfit(x,y,3); %拟合出的三次函数的系数yc3=y-polyval(p,x); % 计算误差s2=sum(yc3.A2); %误差的平方和disp(sprintf('三次误差的平方和=%d',s2);x3=lin