插值拟合逼近.ppt

第1页，共37页，星期日，2025年，2月5日§2.1曲线拟合的最小二乘法常见的数据拟合问题可描述为，给定离散数据表x x1 x2 ?? xm f(x)y1 y2 ?? ym 求拟合函数使得达到最小。这里，函数称为拟合函数，式（2）称为拟合条件。通常假设（1）中是已给定的函数。（1）（2）数据拟合问题第2页，共37页，星期日，2025年，2月5日拟合函数是一元函数时，所对应的函数图形是平面曲线。数据拟合的几何背景是寻找一条近似通过给定离散点的曲线，故称为曲线拟合问题。为了确定数据拟合问题，首先要选取适当的函数类例如，选用幂函数类,则这就是多项式拟合函数。在实际问题中常用的拟合函数类有指数函数类、三角函数类等等，可以根据实验数据的分布特点来选取。第3页，共37页，星期日，2025年，2月5日曲线拟合的最小二乘法设函数已选定，根据拟合条件（2）确定拟合函数（1）中系数的方法称为最小二乘法。拟合函数与数据表中函数在各结点上的差值组成的向量称为残差，记为残差向量r的分量平方和为现在确定使残差平方和最小，可令第4页，共37页，星期日，2025年，2月5日即由于是未知数，将上式整理为方程组（3）称为正规方程组，由它的解可以确定拟合函数（3）第5页，共37页，星期日，2025年，2月5日线性拟合与多项式拟合将拟合函数取线性函数或多项式函数是一种简单的数据拟合方法。拟合函数称为对数据的线性拟合。对于线性拟合问题，需要求函数的最小值。对两个变量求导数，得令其等于零，得正规方程组第6页，共37页，星期日，2025年，2月5日方程组系数矩阵方程组右端项超定方程组:GX=F?正规方程组:GTGX=GTF超定方程组第7页，共37页，星期日，2025年，2月5日例1.已知实验数据如下,求线性拟合函数。x 1 2 3 45f(x)4 4.5 6 8 9第8页，共37页，星期日，2025年，2月5日||r||2=0.7583残差向量:?(1)－4=－0.40?(2)－4.5=0.45?(3)－6=0.30?(4)－8=－0.35?(5)－9=0?(x)=2.25+1.35x第9页，共37页，星期日，2025年，2月5日拟合函数称为对数据的多项式拟合。对于多项式拟合问题，需要求函数采用同样的方法，可得到正规方程组第10页，共37页，星期日，2025年，2月5日方程组系数矩阵方程组右端项超定方程组:GX=F?正规方程组:GTGX=GTF第11页，共37页，星期日，2025年，2月5日例2.求数据的二次拟合函数P(x)=a0+a1x+a2x2x12345f(x)44.5689第12页，共37页，星期日，2025年，2月5日二次拟合误差:||r||2=0.6437比较线性拟合误差:||r||2=0.7583第13页，共37页，星期日，2025年，2月5日超定方程组的最小二乘解给定数据表xx1x2··········xmf(x)y1y2··········ym拟合函数?(x)=a0?0(x)+a1?1(x)+······+an?n(x)为了进一步理解数据拟合的最小二乘方法，考虑如下线性方程组问题。第14页，共37页，星期日，2025年，2月5日以矩阵符号记为：GX=F··············mn+1超定方程组?GTGX=GTF第15页，共37页，星期日，2025年，2月5日系数矩阵按列分块GX=F?GTGX=GTF·····正规方程组的解称为超定方程组的最小二乘解第16页，共37页，星期日，2025年，2月5日取求拟