最小二乘法实验报告.docxVIP

最小二乘法实验报告【实验目的】：观察最小二乘多项式的数值不稳定现象【实验内容】：1 在[-1,1]区间上取n=20个等距节点，计算出以相应节点上的值做为数据样本，以为基函数作出次的最小二乘多项式，画出～之间的曲线，其中A是确定最小二乘多项式的系数矩阵。计算出不同阶最小二乘多项式给出的最小误差 2 在[-1,1]区间上取n=20个等距节点，计算出以相应节点上的值做为数据样本，以为基函数作出次的最小二乘多项式，其中，是勒让德多项式。画出～之间的曲线，其中A是确定最小二乘多项式的系数矩阵。计算出不同阶最小二乘多项式给出的最小误差，把结果与1比较【实验步骤及结果】：在[-1,1]区间上取n=20个等距节点，步长h=2/19，计算出以相应节点上的值做为数据样本，数据如表格1。表格 1 数据样本值ixyixy1-1011012-0012013-0013014-0014015-0015016-0016017-0017018-0018029-00190210-002012以为基函数拟合 在matlab中编写函数lsmex（x, y, l）,生成最小二乘法的系数矩阵A、右端向量d，求出系数，得不同阶数下的最小二乘多项式，其中，。计算系数的结果如下：① l=3， 0.995548986716942 0.997579001689689 0.540351495874469 0.176998749075551② l=5， 1.000038552490093 1.000020023273578 0.499246808308251 0.166497111425869 0.043753810227387 0.008681888224899 ③ l=7， 0.999999832859033 0.999999916830632 0.500005765944358 0.166667853917488 0.041636431725214 0.008328850726684 0.001438513264655 0.000204569079642④ l=9， 1.000000000409680 1.000000000193609 0.499999977559871 0.166666662393254 0.041666859391586 0.008333359277458 0.001388319252092 0.000198349327547 0.000025478122502 0.000002822438546⑤ l=11， 1.000000000008706 1.000000000002729 0.500000000050966 0.166666666671517 0.041666666004638 0.008333333185874 0.001388891754987 0.000198411522433 0.000024795348762 0.000002756249160 0.000000281543798 0.000000025378540⑥ l=13， 0.999999999770700 0.999999998584826 0.500000000002793 0.166666666744277 0.041666666993633 0.008333334699273 0.001388894421398 0.000198394060135 0.000024799884896 0.000002682209015 0.000000272755242 0.000000014901161 0.000000003426662 0⑦ l=15， 1.000000006592917 1.000000012805685 0.500000000015796 0.166666664183140 0