最小二乘法原理及算例.ppt

勒让德多项式(Legendre)[-1,1],?(x)=1递推关系:P0(x)=1,P1(x)=x,第31页,共33页，2024年2月25日，星期天Tn(x)=cos(narccosx)切比雪夫多项式(Chebyshev)递推关系:T0(x)=1,T1(x)=x,T2(x)=2x2-1,T3(x)=4x3-3x,………第32页,共33页，2024年2月25日，星期天感谢大家观看第33页,共33页，2024年2月25日，星期天****关于最小二乘法原理及算例实例讲解某种合成纤维的强度与其拉伸倍数有直接关系，下表是实际测定的24个纤维样品的强度与相应拉伸倍数的记录。提示：将拉伸倍数作为x,强度作为y,在座标纸上标出各点，可以发现什么?第2页,共33页，2024年2月25日，星期天数据表格第3页,共33页，2024年2月25日，星期天第4页,共33页，2024年2月25日，星期天从上图中可以看出强度与拉伸倍数大致成线形关系，可用一条直线来表示两者之间的关系。解：设y*=a+bxi，令δ=yi-y*i=yi-a-bxi，根据最小二乘原理，即使误差的平方和达到最小，也就是令

nQ=∑δi2i=1

为最小，即求使

(a,b)=

有最小值的a和b的值。第5页,共33页，2024年2月25日，星期天计算出它的正规方程得解得：a=0.15,b=0.859

直线方程为：y*=0.15+0.859x第6页,共33页，2024年2月25日，星期天一问题的提出插值法是使用插值多项式来逼近未知或复杂函数的，它要求插值函数与被插函数在插值节点上函数值相同，而在其他点上没有要求。在非插值节点上有时函数值会相差很大。若要求在被插函数的定义区间上，所选近似函数都能与被插函数有较好的近似，就是最佳逼近问题。最佳逼近是在函数空间M中选P(x)满足但由于绝对值函数不宜进行分析运算,常将上式化为来讨论，于是最佳逼近问题变为最佳平方逼近问题,而离散的最佳平方逼进问题就是常说的曲线拟合它们都可用最小二乘法求解。 主页第7页,共33页，2024年2月25日，星期天曲线拟合的最小二乘法最小二乘原理当由实验提供了大量数据时,不能要求拟合函数在数据点处的偏差,即(i=1,2,…,m)严格为零,但为了使近似曲线尽量反映所给数据点的变化趋势,需对偏差有所要求.通常要求偏差平方和最小,此即称为最小二乘原理第8页,共33页，2024年2月25日，星期天?最小二乘法的求法第9页,共33页，2024年2月25日，星期天第10页,共33页，2024年2月25日，星期天?最小二乘法的几种特例第11页,共33页，2024年2月25日，星期天第12页,共33页，2024年2月25日，星期天例题第13页,共33页，2024年2月25日，星期天第14页,共33页，2024年2月25日，星期天二线性最小问题的存在与唯一在科学实验中，很多情况数据间存在线性或可转化为线性的关系。线性最小二乘是最基本也是最重要的一种。1线性最小二乘问题与线性最小二乘求解 设 Ax=b 其中A?Rm?n，b?Rm，x?Rn 当m?n时，上方程超定方程组 令r=b-Ax，一般，超定方程无通常意义下解，既无x使t=0。对这类方程求解意义是求x，使 ??r??22=??b-Ax??22 为最小， 称x为Ax=b的最小二乘解。 主页第15页,共33页，2024年2月25日，星期天2最小二乘解的存在性与唯一性 定理：x*为Ax=b的最小二乘解充要条件 ATAX*=ATb证明：充分性：若存在X*，使ATAX*=ATb则对任意向量令