数值曲线拟合最小二乘法课程设计.docVIP

一、问题叙述 （一）、曲线拟合的最小二乘法定义：对给定数据 (i=0,1,…，m)，在取定的函数类中,求,使误差(i=0,1,…,m)的平方和最小，即 = 从几何意义上讲，就是寻求与给定点(i=0,1,…,m)的距离平方和为最小的曲线（图6-1）。函数称为拟合 函数或最小二乘解，求拟合函数的方法称为曲线拟合的最小二乘法。 在曲线拟合中，函数类可有不同的选取方法. （二）、多项式拟合 假设给定数据点(i=0,1,…,m)，为所有次数不超过的多项式构成的函数类，现求一,使得 (1) 当拟合函数为多项式时，称为多项式拟合，满足式（1）的称为最小二乘拟合多项式。特别地，当n=1时，称为线性拟合或直线拟合。 显然 为的多元函数，因此上述问题即为求的极值 问题。由多元函数求极值的必要条件，得 (2) 即  (3) （3）是关于的线性方程组，用矩阵表示为 (4) 可以证明，方程组（4）的系数矩阵是一个对称正定矩阵，故存在唯一解。从式（4）中解出(k=0,1,…，n)，从而可得多项式  (5) 为所求的拟合多项式。我们把称为最小二乘拟合多项式的平方误差，记作 由式(2)可得 (6) （三）、问题的提出 给出一组数据点列入表7-2中，试用线性最小二乘法求拟合曲线，并估计其误差，作出拟合曲线. 表7–2 一组数据 xi -2.5 -1.7 -1.1 -0.8 0 0.1 1.5 2.7 3.6 yi -192.9 -85.50 -36.15 -26.52 -9.10 -8.43 -13.12 6.50 68.04 问题的分析 多项式拟合的一般方法可归纳为以下几步： (1) 由已知数据画出函数粗略的图形——散点图，确定拟合多项式的次数n； (2) 列表计算和； (3) 写出正规方程组，求出； (4) 写出拟合多项式。 三、程序及实验结果 （1）在MATLAB工作窗口输入程序： x=[-2.5 -1.7 -1.1 -0.8 0 0.1 1.5 2.7 3.6]; y=[-192.9 -85.50 -36.15 -26.52 -9.10 -8.43 -13.12 6.50 68.04]; plot(x,y,r*), legend(实验数据(xi,yi)) xlabel(x), ylabel(y), title(数据点(xi,yi)的散点图) （2）运行后屏幕显示的散点图如下： （3）编写下列MATLAB程序计算在处的函数值，即输入程序 syms a1 a2 a3 a4 x=[-2.5 -1.7 -1.1 -0.8 0 0.1 1.5 2.7 3.6]; fi=a1.*x.^3+ a2.*x.^2+ a3.*x+ a4 运行后屏幕显示关于a1,a2, a3和a4的线性方程组 fi =[ -125/8*a1+25/4*a2-5/2*a3+a4, -4913/1000*a1+289/100*a2-17/10*a3+a4, -1331/1000*a1+121/100*a2-11/10*a3+a4, -64/125*a1+16/25*a2-4/5*a3+a4, a4, 1/1000*a1+1/100*a2+1/10*a3+a4, 27/8*a1+9/4*a2+3/2*a3+a4, 19683/1000*a1+729/100*a2+27/10*a3+a4, 5832/125*a1+324/25*a2+18/5*a3+a4] 编写构造误差平方和的MATLAB程序： y=[-192.9 -85.50 -36.15 -26.52 -9.10 -8.43 -13.12 6.50 68.04]; fi=[-125/8*a1+25/4*a2-5/2*a3+a4, -4913/1000*a1+289/100*a2-17/10*a3+a4, -1331/1000*a1+121/100*a2-11/10*a3+a4, -64/125*a1+16/25*a2-4/5*a3+a4, 1/1000*a1+1/100*a2+1/10*a3+a4, 27/8*a1+9/4*a2+3/2*a3+a4, 19683/1000*a1+729/100*a2+27/10*a3+a4, 5832/125*a1+324/25*a2+18/5*a3+a4]; fy=fi-y; fy2=fy.^2; J=sum(fy.^2) 运行后屏幕显示误差平方和如下