误差理论与数据处理实验报告.doc

《误差理论与数据处理实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《误差理论与数据处理》实验指导书姓名学号机械工程学院2016年05月实验一误差的基本性质与处理一、实验内容1．对某一轴径等精度测量8次，得到下表数据，求测量结果。序号（10-4）1234567824.67424.67524.67324.67624.67124.67824.67224.674-0.00010.0009-0.00110.0019-0.00310.0039-0.0021-0.00010.00020.00770.01270.03520.09770.15020.04520.0002Matlab程序：l=[24.674,24.675,24.673,24.676,24.671,24.

2、678,24.672,24.674];%已知测量值x1=mean(l);%用mean函数求算数平均值disp(['1.算术平均值为：',num2str(x1)]);v=l-x1;%求解残余误差disp(['2.残余误差为：',num2str(v)]);a=sum(v);%求残差和ah=abs(a);%用abs函数求解残差和绝对值bh=ah-(8/2)*0.001;%校核算术平均值及其残余误差,残差和绝对值小于n/2*A,bh<0，故以上计算正确ifbh<0disp('3.经校核算术平均值及计算正确');elsedisp('算术平均值及误差计算有误');endxt=sum(v(1:4))

3、-sum(v(5:8));%判断系统误差（算得差值较小，故不存在系统误差）ifxt<0.1disp(['4.用残余误差法校核，差值为：',num2str(x1),'较小，故不存在系统误差']);elsedisp('存在系统误差');endbz=sqrt((sum(v.^2)/7));%单次测量的标准差disp(['5.单次测量的标准差',num2str(bz)]);p=sort(l);%用格罗布斯准则判断粗大误差，先将测量值按大小顺序重新排列g0=2.03;%查表g(8,0.05)的值g1=(x1-p(1))/bz;g8=(p(8)-x1)/bz;%将g1与g8与g0值比较，g1和g8

4、都小于g0，故判断暂不存在粗大误差ifg1

5、确定度二、实验内容1．由分度值为0 .01mm的测微仪重复6次测量直径D和高度h，测得数据如下：/mm8.0758.0858.0958.0858.0808.060/mm8.1058.1158.1158.1108.1158.110请按测量不确定度的一般计算步骤，用自己熟悉的语言编程完成不确定度分析。MATLAB程序及分析如下：A=[8.0758.0858.0958.0858.0808.060];B=[8.1058.1158.1158.1108.1158.110];D=mean(A);%直径平均值disp(['1.直径平均值为：',num2str(D)]);h=mean(B);%高度平均值

6、disp(['2.高度平均值为：',num2str(h)]);V=pi*D*D*h/4;%体积测量结果估计值disp(['3.体积测量结果估计值为：',num2str(V)]);s1=std(A);%直径标准差disp(['4.直径标准差为：',num2str(s1)]);u1=pi*D*h*s1/2;%直径测量重复性引起的不确定度分量disp(['5.直径测量重复性引起的不确定度分量为：',num2str(u1)]);v1=5;%自由度s2=std(B);%高度标准差disp(['6.高度标准差为：',num2str(s2)]);u2=pi*D*D*s2/4;%高度测量重复性引起的不

7、确定度分量disp(['7.高度测量重复性引起的不确定度分量为：',num2str(u2)]);v2=5;%自由度ue=0.01/(3^0.5);%均匀分布得到的测微仪示值标准不确定度u3=(((pi*D*h/2)^2+(pi*D*D/4)^2)^0.5)*ue;%示值引起的体积测量不确定度disp(['8.示值引起的体积测量不确定度为：',num2str(u3)]);v3=1/(2*0.35^2);%取相对标准差为0.35时对应自由度uc=(