最小二乘法拟合方法的应用

1、CENTRAL SOUTH UNIVERSITY数值分析实验报告最小二乘法拟合的方法的应用一、问题背景在实际的工程应用中，人们往往只能测得一些分散的数据点，为了从这些分散的数据点中找到内在的规律性，就需要利用这些分散的数据点，运用最小二乘法、多项式或其他的已知函数等方法来生成一个新的多项式或是函数来逼近这些已知点。由于MATLAB语言强大的计算功能和绘图功能，使得用户可以很方便的进行曲线拟合并绘制曲线拟合图。二、数学模型我们常说的最小二乘拟合通常指最小二乘多项式拟合。比多项式更一般的拟合函数形式为 其中，为m个函数(多项式拟合中取为幂函数)。假设有n组观测数据，将它们代入上面设定的拟合函数形式

2、中得到 ， 上面的方程组不一定有解，故写成约等号。这里的拟合就是确定参数，的一组值，记为，使得由其中计算得到的数值与观测数据 尽可能接近，这组，可通过解最小化问题得到，即求这种使与的误差平方和在最小二乘意义下最小所确定的函数 称为最小二乘拟合函数。如果定义的拟合模型是关于参数的线性函数，则称为线性模型；如果拟合模型关于参数是非线性函数，则称为非线性模型。在多数情况下，可以通过函数变换的方式将非线性模型转化为线性模型。三、算法及流程对于给定的一组数据，假设拟合函数其中基函数为线性无关的函数系。最小二乘拟合问题就是求系数，使下面的表达式取极小值。根据多元函数极值问题，E取极值的必要条件是设以及，那

3、么-如此，最小二乘拟合问题就转化为超定线性方程组的最小二乘问题。四、计算结果及分析求数据x-3-2-10123y4230-1-2-5的最小二乘拟合。用MATLAB编写主函数如下：% 最小二乘法拟合的主函数 x=-3 -2 -1 0 1 2 3' y=4 2 3 0 -1 -2 -5' % c为要求的变量% 形成线性方程的系数矩阵a= zx_nh_f(x(1) zx_nh_f(x(2)zx_nh_f(x(3) zx_nh_f(x(4) zx_nh_f(x(5) zx_nh_f(x(6) zx_nh_f(x(7); b=y; % 法方程A=a'*a; B=a'*b;

4、 c=AB % 在MATLAB中实际上可以直接使用c=ab，其实以及进行了最小二乘% c就是拟合基函数系数% 作出原离散数据与你和函数图像 x_n=-3:0.02:3; y_n=c(1)*1+c(2)*x_n+c(3)*x_n.2; plot(x,y,'*',x_n,y_n) grid 将文件以文件名zx_nh.m保存。再编写拟合的基函数程序：% 最小二乘法的拟合基函数 function f=zx_nh_f(x) f(1)=1; f(2)=x; f(3)=x2; 将文件以文件名zx_nh_f.m保存。在MATLAB命令窗口输入：>>zx_nh回车得到：c = 0.666666666666666 -1.392857142857143 -0.130952380952381则由最小二乘法拟合得到的函数为：运行程序还得到原离散数据与拟合函数图像，如下图所示，数据拟合与插值的最大区别在