第六章数据拟合方法.pptVIP

第1页，共27页，星期日，2025年，2月5日第六章数据拟合方法数据拟合的最小二乘法Bezier曲线第2页，共27页，星期日，2025年，2月5日例：考察某种纤维的强度与其拉伸倍数的关系,下表是实际测定的24个纤维样品的强度与相应的拉伸倍数的记录:纤维强度随拉伸倍数增加而增加。6.1 数据拟合的最小二乘法一、曲线拟合的数学描述与问题求解第3页，共27页，星期日，2025年，2月5日24个点大致分布在一条直线附近。故可认为强度y与拉伸倍数x的主要关系应为线性关系：必须找到一种度量标准来衡量什么曲线最接近所有数据点。第4页，共27页，星期日，2025年，2月5日xx1x2……xmf(x)y1y2……ym1、数据拟合问题研究内容：从一大堆看上去杂乱无章的数据中找出规律性来，即设法构造一条曲线（拟合曲线）反映所给数据点总的趋势，以消除其局部波动。这种要求曲线尽可能逼近给定数据的过程称“拟合”。给定一组值：求函数使得最小。第5页，共27页，星期日，2025年，2月5日据实验数据分布特点选取，可选幂函数类、指数函数类、三角函数类等。（1）若?(x)为一元函数，则函数曲线为平面图形，称曲线拟合。（2）?(x)为拟合函数，上式最小为拟合条件（即要求拟合曲线与各数据点在y方向的误差平方和最小）。（3）函数类的选取：说明：第6页，共27页，星期日，2025年，2月5日残差向量的各分量平方和记为：2、最小二乘法：以残差平方和最小问题的解来确定拟合函数的方法。令－－在回归分析中称为残差(i=1,2,…m)残差向量：第7页，共27页，星期日，2025年，2月5日由多元函数求极值的必要条件，有可得即第8页，共27页，星期日，2025年，2月5日上式为由n+1个方程组成的方程组，称正规方程组。由得即第9页，共27页，星期日，2025年，2月5日引入记号则由内积的概念可知显然内积满足交换律正规方程组便可化为第10页，共27页，星期日，2025年，2月5日将其表示成矩阵形式：其系数矩阵为对称阵。所以正规方程组的系数矩阵非奇异,即根据Crame法则,正规方程组有唯一解，称其为最小二乘解。第11页，共27页，星期日，2025年，2月5日作为一种简单的情况，常使用多项式函数Pn(x)作为(xi,yi)(i=1,2,…,m)的拟合函数。基函数之间的内积为：拟合函数φ(x)=Pn(x)的基函数为：第12页，共27页，星期日，2025年，2月5日即正规方程组为第13页，共27页，星期日，2025年，2月5日例.回到本节开始的实例，从散点图可以看出，纤维强度和拉伸倍数之间近似线性关系，故可选取线性函数为拟合函数建立正规方程组，其基函数为根据内积公式，可得正规方程组为第14页，共27页，星期日，2025年，2月5日解得残差平方和：拟合曲线与散点的关系如右图:即为所求的最小二乘解。故第15页，共27页，星期日，2025年，2月5日若mn+1，则此方程组称超定方程组（方程个数未知数个数）二、超定方程组的最小二乘解将拟合函数以向量表示：令（i=1,2,…,m）可得第16页，共27页，星期日，2025年，2月5日考虑正规方程组（k=0,1,…n）（1）未知数aj的系数为超定方程组中系数阵第k列与第j列对应积之和（即内积(φk,φj)）；（2）右端向量为系数阵第k列与m个函数对应积之和。可知：第17页，共27页，星期日，2025年，2月5日故正规方程组矩阵形式为：若有唯一解，称其为超定方程组的最小二乘解。注：最小二乘解并不能满足超定方程组中每个方程，但要求尽可能接近给定数据，即允许每个等式可以稍有偏差（即残差）。求一般超定方程组Ax=b的主要过程：（1）求出系数矩阵A的转置矩阵AT；（2）计算矩阵D=ATA和向量f=ATb；（3）求解正规方程组Dx=f。第18页，共27页，星期日，2025年，2月5日x12 3 4y4101826例1用多项式拟合函数：解：设得即第19页，共27页，星期日，2025年，2月5日