观察最小二乘多项式不稳定现象.doc

观察最小二乘多项式的数不稳定现象实验1实验任务1.1 在-1,1区间上取n=20个等距节点，计算出以相应节点上的ex的值作为数据样本，以1，x，x2，xl为基函数做出l=3,5,7,9次的最小二乘拟合多项式。1.2 画出lncondA-l曲线，其中A是确定最小二乘多项式系数的矩阵。1.3 计算出不同阶最小二乘多项式给出的最小偏差l。1.4 将基函数改为1，P1x，P2x，Plx，其中Pix是勒让德多项式，结果如何？2实验原理与理论基础2.1一般线性最小二乘拟合的法方程组为：0,0 0,1 0,n1,0 1,1 1,n n,0 n,1 n,na0a1an=0,y1,yn,y由于以1，x，x2，xl为基函数，所以i=xi ， i=0,1l。把i=xi ， i=0,1l代入一般线性最小二乘拟合的法方程组中可得多项式拟合的法方程组。2.2最小二乘拟合多项式的存在唯一性：定理1：设节点x0，x1，xn互异，则多项式拟合的法方程组的解存在唯一性。定理2：设a0，a1，an是多项式拟合的法方程组的解，则Pnx=k=0nakxk是最小二乘拟合多项式。2.3 矩阵A的条件数是cond(A)=AA-1。根据范式的不同矩阵的条件数也有3中，这里选取A为A的2-范式。2.4勒让德多项式L0=1，L1=x ，Ln+1=2n+1n+1xLnx-nn+1Ln-1x，n13实验内容及实验结果1.基函数为1，x，x2，xl3.1在-1,1区间上取n=20个等距节点，计算出以相应节点上的ex的值作为数据样本。xi=-1+i*h，h=219，i=0，1，19yi=exi，i=0，1，193.2计算法方程组对应得系数矩阵，及增广矩阵。系数矩阵的的第i行，第j列的元素为ai,j=xi-1,xj-1=k=019xki-1*xkj-13.3求系数矩阵的条件数。cond(A)=AA-1。3.4化简增广矩阵为简化行阶梯型矩阵，得出法方程组的解。其中第一行是三阶最小二乘多项式的系数，第二行是五阶最小二乘多项式的系数等等。所以三阶最小二乘多项式为：y=0.9955+0.9976x+0.5404x2+0.1770x3五阶最小二乘多项式：y=1+x+0.4992x2+0.1665x3+0.0438x4+0.0087x5七阶最小二乘多项式：y=1+x+0.5+0.1667+0.0416x4+0.0083x5+0.0014x6+0.00020457x7九阶最小二乘多项式同样代入系数可得。3.5 cond(A)的取值如下图所示，分别是3,5,7,9,11,13,15阶行列式对应得A的条件数：画出lncondA-l的曲线：分析lncondA-l的曲线可看出，随着l增大，cond(A)迅速增大。这意味着当l越大时，正规方程组的病态越严重。3.6计算出不同阶最小二乘多项式给出的最小偏差l。l=i=019yxi-yi2，l=3,5,7,93=0.000313692373446139，5=2.24566107559548e-087=4.01160377773610e-13，9=2.25819731366317e-18可以看出随着多项式阶数的增大，其最小偏差在迅速减小。这似乎是和3.5中的结论“l越大时，正规方程组的病态越严重”相悖。3.7为了进一步研究这个问题，我们取l=3,5,7,9,11,13下图是cond(A)的取值，分别是3,5,7,9,11,13,15阶行列式对应得A的条件数。 下图是lncondA-l曲线，可以发现曲线走势基本与l=3,5,7,9时确定的曲线走势一致。计算11=7.5119e-24，13=9.90515e-23，15=3.4205e-21对比13，5，7，9的值我们可以得出结论，随着拟合多项式次数的增大A的条件数迅速增大，确定最小二乘多项式的方程组的病态程度也随之增加。这在偏差中反映出来的便是，刚开始时次数增大偏差减小，当次数达到一定程度后，次数越大，偏差反而越大。2基函数为勒让德多项式3.8把基函数改为勒让德多项式组成的基函数：L0=1，L1=x ，Ln+1=2n+1n+1xLnx-nn+1Ln-1x，n1所以定义函数function out = L(n,x) if n=0 out=1; elseif n=1 out=x; else out=(2*n+1)/(n+1)*x*L(n-1,x)-n/(n+1)*L(n-2,x); endend3.9 用L(n,x)取代原程序中的 xn 求出最小二乘拟合函数的系数，其中第一行为三阶最小二乘拟合函数的系数，第二行为五阶最小二乘拟合函数的系数，等等每一行从左到右分别为a0，a1，al。把系数代入下式，即可得到各阶最小二乘拟合函数。y=k=0lakLk(x)3.10 cond(A)的取值如下图所示，分别是3,5,7,9,11,13,15阶行列式对应得A的条件数：作出画出lncondA-l的曲线：可以看出以勒让德多项式为基函数得到的法方程组的系数矩阵的条件数增长，相对于多项式拟合的法方程组的系数矩阵的条件数增长要缓慢得多。所以可以说随着阶数增长，以勒让德多项式为基函数得到的法方程组的病态，相对于多项式拟合的法方程组的病态，不算严重。3.11计算出不同阶最小二函数给出的最小偏差l。l=i=019yxi-yi2，l=3,5,7,93=0.49328797611472