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

1、精品文档观察最小二乘多项式的数不稳定现象实验1 实验任务1.1 在-1,1区间上取 n=20 个等距节点，计算出以相应节点上的?的值作为数2，?据样本，以 1，x，x，?为基函数做出 l = 3,5,7,9次的最小二乘拟合多项式。1.2 画出 ln(?) ) - ?曲线，其中 A 是确定最小二乘多项式系数的矩阵。1.3 计算出不同阶最小二乘多项式给出的最小偏差(?) 。1.4 将基函数改为 1，?1( ?)，?2( ?)， ? ，?(?)，其中 ?(?)是勒让德多项式，结果如何？2 实验原理与理论基础2.1 一般线性最小二乘拟合的法方程组为：(? ,? )( ?,?)?(?,?)?( ?, y

2、)000010 ?10(?1 ,?0 )( ?1,?1)?(?1, ?)( ?1,y) ? ?=?(?,?)(?, ?)?( ?,?)?(? , y)? 0?1? ?2?由于以 1， x，x，? ， ?为基函数，所以 ?= ?， ?= 0,1 ? l 。?把 ?= ?， ?= 0,1 ? l 代入一般线性最小二乘拟合的法方程组中可得多项式拟合的法方程组。2.2 最小二乘拟合多项式的存在唯一性：定理 1：设节点 ?，?， ? ， ?互异，则多项式拟合的法方程组的解存在唯一性。01?定理 2：设?，?， ? ，?是多项式拟合的法方程组的解， 则?01?n(x) =?=0?是最小二乘拟合多项式。-1

3、。根据范式的不同矩阵的条件数也2.3 矩阵 A 的条件数是 cond(A)= ? ? ?有 3 中，这里选取 ? 为 A的 2- 范式。2.4 勒让德多项式?2?+1( )?( )，0 = 1，?1= ?， ?+1 =?+1 ?-1n1?+1 ? -?。1欢迎下载精品文档3 实验内容及实验结果?1. 基函数为 1， x，?，? ，?3.1 在-1,1区间上取 n=20 个等距节点，计算出以相应节点上的ex 的值作为数据样本。?= -1 + ? ?，? =2，i = 0，1， ? ，1919?， i = 0， 1， ? ， 19?= ?3.2 计算法方程组对应得系数矩阵，及增广矩阵。系数矩阵的的

4、第i 行，第 j 列的元素为19?-1 ?-1) = ?-1?-1?,?= (? , ?=03.3 求系数矩阵的条件数。cond(A)= ?-1 。3.4 化简增广矩阵为简化行阶梯型矩阵，得出法方程组的解。其中第一行是三阶最小二乘多项式的系数， 第二行是五阶最小二乘多项式的系数等等。所以23三阶最小二乘多项式为： y = 0.9955 + 0.9976x + 0.5404? + 0.1770?2345五阶最小二乘多项式： y = 1 + x + 0.4992?+ 0.1665? + 0.0438? + 0.0087?45七 阶 最 小 二 乘 多 项 式 ： y = 1 + x + 0.5 +

5、 0.1667 + 0.0416? + 0.0083? +670.0014? + 0.00020457?九阶最小二乘多项式同样代入系数可得。3.5 cond(A) 的取值如下图所示，分别是 3,5,7,9,11,13,15 阶行列式对应得 A的条件数：。2欢迎下载精品文档画出 ln(?) -?的曲线：分析 ln(?) -?的曲线可看出，随着l 增大， cond(A) 迅速增大。这意味着当 l 越大时，正规方程组的病态越严重。3.6 计算出不同阶最小二乘多项式给出的最小偏差(?) 。19(?)= y ( ?) -?2 ，l = 3,5,7,9?=0(3) = 0.0003136923734461

6、39 ，(5) = 2.24566107559548? - 08 (7) = 4.01160377773610? - 13， (9) = 2.25819731366317? - 18 可以看出随着多项式阶数的增大，其最小偏差在迅速减小。这似乎是和 3.5 中的结论“ l 越大时，正规方程组的病态越严重”相悖。3.7 为了进一步研究这个问题，我们取l = 3,5,7,9,11,13下图是 cond(A) 的取值，分别是 3,5,7,9,11,13,15阶行列式对应得 A 的条件数。下图是 ln(?) - ?曲线，可以发现曲线走势基本与l = 3,5,7,9时确定的曲线走势一致。3欢迎下载精品文档

7、计算 (11) = 7.5119?- 24， (13) = 9.90515?- 23，(15) = 3.4205?- 21对比 (13 )， (5)，( 7) ，(9) 的值我们可以得出结论，随着拟合多项式次数的增大 A 的条件数迅速增大，确定最小二乘多项式的方程组的病态程度也随之增加。这在偏差中反映出来的便是，刚开始时次数增大偏差减小，当次数达到一定程度后，次数越大，偏差反而越大。2 基函数为勒让德多项式3.8 把基函数改为勒让德多项式组成的基函数：? = 1，? = ?，?=2?+ 1(?) -?(?) ，n 101?+1?+ 1 ?+ 1?-1所以定义函数functionout = L(

8、n,x)ifn=0out=1;elseifn=1out=x;elseout=(2*n+1)/(n+1)*x*L(n-1,x)-n/(n+1)*L(n-2,x);endend3.9 用L(n,x) 取代原程序中的? 求出最小二乘拟合函数的系数，其中第一行为三阶最小二乘拟合函数的系数，第二行为五阶最小二乘拟合函数的系数，等等每一行从左到右分别为 ?，?， ? ， ?。01?。4欢迎下载精品文档把系数代入下式，即可得到各阶最小二乘拟合函数。?y = ?(?)?=03.10 cond(A) 的取值如下图所示，分别是 3,5,7,9,11,13,15 阶行列式对应得 A的条件数：作出画出 ln(?) - ?的曲线：可以看出以勒让德多项式为基函数得到的法方程组的系数矩阵的条件数增长，相对于多项式拟合的法方程组的系数矩阵的条件数增长要缓慢得多。所以可以说随着阶数增长，以勒让德多项式为基函数得到的法方程组的病态，相对于多项式拟合的法方程组的病态，不算严重。3.11 计算出不同阶最小二函数给出的最小偏差(?)。19(?)= y ( ?) -?2 ，l = 3,5,7,9?=0(3) = 0.493287976114723 ， (5) = 0.493844031116877 (7) = 0.493805766121517 ， (9) = 0.493805