数学实验报告tzx..docVIP

数学实验报告 实验序号： 日期： 2012 年 5 月 2 日 班级 数学系09B班 姓名 童周祥 学号 20091611241 实验名称 使用matlab编写程序计算最小二乘法的拟合 问题背景描述： 实验目的： １、学习使用matlab编写数值计算程序。 ２、了解最小二乘法的拟合原理和编程方法。 ３、根据最小二乘法的拟合原理编写matlab程序，并运行的出结果。 实验原理： 设A是阶矩阵，称线性方程组： 　　　　　　　　　　　　　　　　　　(1) 为超定方程组，这里。如果Ａ的秩r(A)=n，称Ａ为列满秩矩阵。记残向量 ｒ=b-Ax，考虑确定一个向量ｘ，使 　　　　　　　　　　　　　　　　 达到最小的问题称为线性最小二乘问题，这样的ｘ称为方程组（１）的最小二乘解。线性最小二乘问题与欧式空间的正交性等有密切联系。 设超定方程组（１）的系数矩阵A是列满秩矩阵，由定理（当mn时，超定方程组（１）的最小二乘解总是存在。最小二乘解惟一的充分必要条件是）知，方程组（１）的最小二乘解总是存在而且惟一的。 设x是最小二乘解，。由定理（当mn时，超定方程组（１）的最小二乘解总是存在。最小二乘解惟一的充分必要条件是），残向量 ， 因此，　　（２） 也就是说。X是方程组（１）的最小二乘解等价于是方程组（２）的解。 方程组（２）称为最小二乘问题的正则方程组或法方程组，由于是列满秩的，因此，方程组（２）是正定方程组。 用正则方程组求最小二乘解的计算步骤为如下： ①形成； ②计算的分解，记； ③求解下三角方程组； ④求解上三角方程组。 主要内容： 依照实验原理编写最小二乘法的程序。 实验过程记录： 首先，在电脑上安装matlab，然后，启动matlab，新建一个M文件。 实验程序如下： 程序代码 function c=ercheng(x,y,m) for i=1:m+1 a(i,:)=x.^(i-1); end =a*a; b=a*y; c=inv()*b; x0=-100:100; y0=0; for i=1:m+1 y0=y0+c(i).*x0.^(i-1); end plot(x,y,+,x0,y0,k:) 四、实验数据及结果： １、实验送的初值为： 　　　　　　　x=[-100 -75 -50 -25 0 25 50 75 100]; 　　　　　　　y=[4.06 6.78 9.49 16.27 40.67 97.62 146.43 151.85 162.7]; 　　　　　　　m=4; 　　实验结果为： ＝49.2061　　　　　　＝1.4692　　　　　　＝0.0104 ＝-0.0001　　　　　　＝-0.0000 　　拟合波形为： ２、实验送的初值为： 　　　　　　　x=[-100 -75 -50 -25 0 25 50 75 100]; 　　　　　　　y=[ 3.855 7.555 25.554 36.245 45.547 65.558 98.587 129.548 256.3525]; 　　　　　　　m=4; 　　实验结果为： ＝51.2424　　　　　　＝0.4511　　　　　　＝-0.0007 ＝0.0001　　　　　　 ＝0.0000 　　拟合波形为： ３、实验送的初值为： 　　　　　 　x=[-100 -75 -50 -25 0 25 50 75 100]; 　　　 　　　y=[25.6 65.4 75.58 50.63 17.58 47.57 12.55 35.25 12.35]; 　　　 　　　m=4; 　　实验结果为： ＝32.8234　　　　　　＝-0.5284　　　　　　＝-0.0087 ＝0.0000　　　　　　 ＝-0.0000 　　拟合波形为： 五、实验分析： １、在本次实验中所输的初值ｘ与ｙ对应的不太好时，就会出现像初值３那样的，拟合点和 　　拟合曲线之间有很大的差距，但基本上拟合点还是均匀分布在曲线周围的，初值１和初 　　值２的拟合曲线和拟合点均对应的很好。基本在曲线上，或是曲线附近很小的范围内。 ２、这种算法只适用于拟合函数是多项式的情形。 　　用正则方程组求最小二乘解得方法，有两个注意的问题： 　　　　　　①形成正则方程组时，计算出来的可能不是正定的，也可能是奇异的，以致 　　　　　　　后面不能用分解来求解正则方程组。比如说： 　　　　　　　　　　　　　　　　　　 　　　　　　　采用六位有限数位计算，得，这时是奇异的。 　　　　②矩阵Ａ的状态数与的状态数密切相关，对列满秩的阶矩阵Ａ， 　　　　　有因此，用正则方程组求最