误差理论实验报告2.docx

《误差理论与数据处理》实验报告实验名称：线性函数的最小二乘法处理实验目的线性函数的最小二乘法是解决有关组合测量最佳估计问题的典型的数据处理方法。本实验要求学生编写最小二乘数据处理程序并对组合测量数据进行处理，求出最佳估计值并进行精度分析。实验原理1.最小二乘法原理指出，最可信赖值应在是残差误差平方和的条件下求得。?2.最小二乘法可以将误差方程转化为有确定解的代数方程组（其方程组的数目正好等于未知数的个数），从而可求解出这些未知参数。这个有确定解的代数方程组称为最小二乘法的正规方程。?3.线性参数的最小二乘法处理程序为：首先根据具体问题列出误差方程式；再按最小二乘原理，利用求极值的方法将误差方程转化为正规方程；然后求解正规方程，得到代求的估计量；最后给出精度估计。?4.正规方程又转化为残差方程，残差方程可用矩阵方法求出方程的解。因此可用Matlab求解最小二乘法参数。?5.求出最小二乘法的参数后，还要对参数进行精度估计。?相应的标准差为ttxtxxddd222111，其中ttddd..2211称为不定乘数。实验内容和结果程序及流程在MATLAB环境下建立一个命令M-文件，编写解答以下组合测量问题数据处理的程序：现要检定刻线A,B,C,D间的距离x1,x2,x3，采用组合测量方法，直接测量刻线间的各种组合量，得到数据如下测量数据：l1=1.051mm; l2=0.985; l3=1.020mm; l4=2.016mm; l5=1.981mm; l6=3.032mm编程求x1,x2和x3的最小二乘估计值；对直接测量数据进行精度估计对x1,x2和x3的最小二乘估计值进行精读估计。程序： A=[1 0 0;0 1 0;0 0 1;1 1 0;0 1 1;1 1 1] A*A C=A*A inv(C) l=[1.015; 0.985; 1.020; 2.016; 1.981;3.032]; X=inv(C)*A*l V=l-A*X V*V STD1=sqrt(V*V/3) inv(C) STDX1=sqrt(0.5)*STD1实验结果（数据或图表）结果分析心得体会通过本次实验，我掌握等精度测量线性参数最小二乘法的处理，并能够应用Matlab用矩阵的方法求出拟合方程的参数，及能够对各个参数进行精度估计。同时能根据等精度线性参数理解不等精度线性参数及非线性参数情况下的最小二乘法处理。对以后的学习有了很大的帮助《误差理论与数据处理》实验报告实验名称：一元/多元回归数据分析实验目的回归分析是对实验数据进行处理的重要方法。通过本实验使学生掌握一元线性回归方程的求解和方差分析、显著性检验方法；掌握一元非线性回归方程的求解和显著性检验的方法；掌握多元线性回归方程的求解和方差分析、显著性检验方法；掌握回归数据处理的程序设计方法。实验原理回归分析是研究随机现象中变量之间相关关系的一种统计方法。一元线性回归一元线性回归就是研究两个具有线性相关关系的随机变量之间的依存关系。即求取它的经验公式。1.一元线性回归的数学模型： (i=1,2,3……n)其中表示随机因素对影响总和，一般假设他们是一组相互独立，并服从同一正态分布N（0、）的随机变量。Xi是可以严格控制的变量：yi是服从正态分布N（b0+b1xi,）的随机变量。b0,b1是待估参数。一元线性回归方程：x利用最小二乘法可求得方差分析：误差来源平方和自由度方差F显著性回归U1S^2=F=0.010.050.1或其他残差Qn-2总计Sn-12.多元线性回归多元线性回归的数学模型假设因变量y与另外m个自变量的内在联系是现行的通过实验得到n组观测数据：（xi1,xi2……,xim;yi） (i=1,2,……,n)那么这批数据有如下的结构形式：其中(i=1,2,3,……,n)是一组相互独立，并服从同一正态分布N（0,）的随机变量。（i=1,2,3,……,n）是可以严格控制的变量；bi（i=0,1,2,……,m）是待估参数。多元线性回归方程:其中X=；ｂ＝利用最小二乘法可求得：ｂ＝（　方差分析误差来源平方和自由度方差F显著性回归UＭS^2=F=0.010.050.1或其他残差Qn-ｍ－１总计Sn-1实验内容和结果程序及流程用ＭＡＴＬＡＢ编写程序解答下面各题材料的抗剪强度与材料承受的正应力有关。某种材料实验数据如下表：正应力ｘ（Ｐａ）２６．８２５．４２８．９２３．６２７．７２３．９２４．７２８．１２６．９２７．４２２．６２５．６抗剪强度ｙ（Ｐａ）２６．５２７．３２４．２２７．１２３．６２５．９２６．３２２．５２１．７２１．４２５．８２４．９编写程序，要求求出ｙ对ｘ一元钱线性回归方程和ｘ对ｙ的一元线性回归方程并在同一张图上绘出测量数据的散点以及所拟合的两条直线图。x=[26.8;25.4;28.9;23.6;