数值分析曲线拟合的最小二乘法实验报告精选 .doc

数值分析曲线拟合的最小二乘法实验报告 一、目的和意义 在科学实验的统计方法研究中，往往要从一组实验数据?xi,yi??i?0,1,2,?,m?中，寻找自变量x与因变量y之间的函数关系y?F?x?。由于观测数据往往不准确，因此不要求y?F?x?经过所有点?xi,yi?，而只要求在给定点xi上误差而只要求所在所有给定点xi上的误差?i?F?yi ?i?0,1,2,?,m?按某种标准最小。若记????0,?1,?2,?,?m?，就是要求向量?的范数 如果用最大范数，计算上困难较大，通常采用欧式范数?最小。 2T 作为误差度量的标准。F?x?的函数类型往往与实验的物理背景以及数据的实际分布有关，它一般含有某些待定参数。如果F?x?是所有待定参数的线性函数，那么相应的问题称为线性最小二乘问题，否则称为非线性最小二乘问题。最小二乘法还是实验数据参数估计的重要工具。这是因为这种方法比其他方法更容易理解，即使在其他方法失效的情况下，用最小二乘法还能提供解答，而且从统计学的观点分析，用该方法求得各项估计具有最优统计特征，因此这一方法也是系统识别的重要基础。线性最小二乘问题可以借助多元微分学知识通过求解法方程组得到解答。 用最小二乘法求拟合曲线时，首先要确定S?x?的形式。这不单纯是数学问题，还与所研究问题的运动规律以及所得观测数据?xi,yi?有关；通常要从问题的运动规律以及给定数据描图，确定S?x?的形式，并通过实际计算选出较好的结果。为了使问题的提法更有一般性，通常把最小二乘法中的? 22 都考虑为加权平方和 22 ? ????xi???S?xi??f?xi??? i?0 m 2 这里??xi??0是?a,b?上的加权函数，它表示不同点?xi,f?xi?处的数据比重不同。 ? 二、计算方法 在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量y与时间t的拟合曲线。 本题要求我们用 ??t??a1t?a2t2?a3t3对曲线进行拟合，这里 m?11,?1?t,?2?t2,?3?t3,故 ??1,?1???t i?0 11 2i ?12650,??1,?2????2,?1???ti3?544500, i?011 11 ??2,?2???t i?0 11 4i ??1,?3????3,?1? ??ti4 i?0 11 5i ??3,?2????2,?3???t i?011 ?1193362500,??3,?3???ti6?58593218750 i?0 11 11 ??1,y???tiyi?,??2,y???ti2yi?, i?0 i?0 ??3,y???ti3yi? i?0 11 由于 ???,??a k j j?0 n j ?dj?k?0,1,?,n?, 可以利用此式算出拟合曲线的ai，即 ???1,?1???1,?2???1,?3???a1????1,y1?????????,??,??,?a??,y?????????2122?22? 23??2? ??,??,??,???a???,y???31??32??33???3???33?? 所以求得 5 ，a1?2.?6?61a02???10?7 ， a3??10?9 ， ??t???10?5t??10?7t2??10?9t3，误差为?i?yi???ti??i?0,1,?11? max?i? ? 下图可见实际测出值与拟合值的差别，下表可见拟合出的每一点的误差以及均方误差。 t 0 5 10 15 20 25 30 35 40 45 50 55 y 0 拟合值 0 误差 0 - - - - - 误差平方 均方误差 三、结构程序设计 在本题使用Visual Studio c# .NET编译程序。 //在窗体的Load事件里调用InitialDeal和Deal函数来处理数据。 private void ResultReport_Load { ; ; } //初始化各个变量。 private