昆明理工大学数值分析上机报告7

《昆明理工大学数值分析上机报告7》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数值分析实验报告（数值积分）姓名：学号：2006202080专业：材料学学院：材料与冶金工程学院授课教师：昆明理工大学理06工科硕士《数值分析》上机实验报告专业：材料学姓名：学号：2006202080任课教师：作业完成实验室：实验成绩：理论描述(10)数值公式(10)程序流程图和程序结构(20)数据和结果(20)讨论(20)源程序(20)总分(100)曲线拟合的最小二乘法一、目的和意义在科学实验的统计方法研究中，往往要从一组实验数据中，寻找自变量x与因变量y之间的函数关系。由于观测数据往往不准确，因此不要求经过所有点，而只要求在给定点上

2、误差而只要求所在所有给定点上的误差按某种标准最小。若记，就是要求向量的范数最小。如果用最大范数，计算上困难较大，通常采用欧式范数作为误差度量的标准。的函数类型往往与实验的物理背景以及数据的实际分布有关，它一般含有某些待定参数。如果是所有待定参数的线性函数，那么相应的问题称为线性最小二乘问题，否则称为非线性最小二乘问题。最小二乘法还是实验数据参数估计的重要工具。这是因为这种方法比其他方法更容易理解，即使在其他方法失效的情况下，用最小二乘法还能提供解答，而且从统计学的观点分析，用该方法求得各项估计具有最优统计特征，因此这一方法也是系统识别的

3、重要基础。线性最小二乘问题可以借助多元微分学知识通过求解法方程组得到解答。用最小二乘法求拟合曲线时，首先要确定的形式。这不单纯是数学问题，还与所研究问题的运动规律以及所得观测数据有关；通常要从问题的运动规律以及给定数据描图，确定的形式，并通过实际计算选出较好的结果。为了使问题的提法更有一般性，通常把最小二乘法中的都考虑为加权平方和这里是上的加权函数，它表示不同点处的数据比重不同。二、计算方法分）051015202530354045505501.272.162.863.443.874.154.374.514.584.024.64在某冶炼过

4、程中，根据统计数据的含碳量与时间关系，试求含碳量与时间的拟合曲线。本题要求我们用对曲线进行拟合，这里故,,,,,,,由于,可以利用此式算出拟合曲线的，即所以求得，，，，误差为，而均方误差为下图可见实际测出值与拟合值的差别，下表可见拟合出的每一点的误差以及均方误差。ty拟合值误差误差平方0000051.271.20215-0.067850.004603623102.162.16620.00623.844E-05152.862.918550.058550.003428103203.443.48560.04560.00207936253.87

5、3.893750.023750.000564063304.154.16940.01940.00037636354.374.33895-0.031050.000964102404.514.4288-0.08120.00659344454.584.46535-0.114650.013144622504.024.4750.4550.207025554.644.48415-0.155850.024289222均方误差0.51293892三、结构程序设计在本题使用VisualStudioc#.NET编译程序。//在窗体的Load事件里调用Init

6、ialDeal和Deal函数来处理数据。privatevoidResultReport_Load(objectsender,System.EventArgse){this.InitialDeal();this.Deal();}//初始化各个变量。privatevoidInitialDeal(){inti=0;while(tString.Length>0){t[i++]=Convert.ToDouble(tString.Substring(0,tString.IndexOf(",")))*tUnit;tString=tString.Rem

7、ove(0,tString.IndexOf(",")+1);}num=i;i=0;while(yString.Length>0){y[i++]=Convert.ToDouble(yString.Substring(0,yString.IndexOf(",")))*yUnit;yString=yString.Remove(0,yString.IndexOf(",")+1);}i=0;while(i

8、0);0.00");i=0;while(i