最小二乘抛物曲线拟合

课程：数值方法班级：2009级数学与应用数学（师范）1班姓名：廖璐学号：200902114020最小二乘多项式拟合摘要本次实验通过利用最小二乘法的原理以及多项式讨论最小二乘抛物曲线拟合可以得到最小二乘多项式拟合，给出了实例的具体做法。关键词曲线拟合最小二乘法多项式抛物线正文一.最小二乘法原理设已知某过程的一组观测数据要求在某特定函数类寻找一个函数φ（x）作为的近似函数，使其在xk上的误差ek=φ（xk）-f(xk),k=1,2,...,N按某种形式为最小这就是拟合问题。要求误差的平方和最小的拟合为曲线拟合的最小二乘法。也可以表示为均方根误差：的值最小。从几何意义上讲：就是数据点到曲线的垂直距离平方和的最小值。多项式拟合假设给定数据点，ψ为所有次数不超过m(m<N)的多项式构成的函数类，现求一。当拟合曲线为多项式时，称为多项式拟合，满足I的fn(x)称为最小二乘拟合多项式。特别的当n=1时，称为线性拟合或者直线拟合。本文具体探讨当n=2的情形，也就是抛物曲线拟合三.最小二乘抛物曲线拟合实验目的：1.掌握拟合曲线的最小二乘法的原理和基本内容2.会利用最小二乘法求出抛物线实验要求：1.了解最小二乘拟合曲线y=Ax2+Bx+C的原理和方法；2.根据具体点求解最小二乘抛物线系数3.绘制出曲线拟合图实验原理：由最小二乘抛物线y=f(x)=Ax2+Bx+C，通过使的值最小（即是令偏导数ɑE/ɑA,ɑE/ɑB,ɑE/ɑC为零）得到的函数就是拟合函数（也称最小二乘解）实验基本步骤：1.设有N个点，横坐标是确定的。首先建立最小二乘抛物线的系数表示为y=f(x)=Ax2+Bx+C。通过求表达式E（A,B,C）=的最小值可得到A，B和C。令偏导数ɑE/ɑA,ɑE/ɑB,ɑE/ɑC为零，可得利用加法分配律，可将上述中的A,B,C移到求和的外面，以得到正规方程为了在matlab中实现可以表示成以下矩阵形式：左边=右边=因此，=１\*GB3①首先建立一个3×N的零矩阵B（N=length(x)以存放；=２\*GB3②接着令F等于一个Ｎ×３的零矩阵用来存放；=３\*GB3③然后给F赋值F(:,k)=x'.^(k-1)；计算F和F'；=4\*GB3④最后实验程序：functionC=lspoly(X,Y,M)n=length(X);B=zeros(1:M+1);F=zeros(n,M+1);fork=1:M+1F(:,k)=Ｘ'.^(k-1);endA=F'*F;B=F'*Y';C=A\B;C=flipud(C);输入数据：X=[-3024]X=-3024>>Y=[3113]Y=3113>>M=2M=2输出结果：>>C=lspoly(X,Y,M)C=0.1785-0.19250.8505因此结束语：通过实验可以更深刻的理解曲线拟合最小二乘法的基本原理，也更深刻的理解了数值方法应用的广泛性。本学期学习《数值方法》这门课，虽然在听课的时候很多东西都没有接触过感觉学习比较困难许多地方都不懂，但是我觉得自己还是有一些收获（比如在matlab的操作上有了进步；了解了数