欢迎来到天天文库
浏览记录
ID:59494245
大小:617.50 KB
页数:12页
时间:2020-11-03
《误差理论与数据处理-实验报告.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=[2
2、4.674,24.675,24.673,24.676,24.671,24.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,故以上计算
3、正确ifbh<0disp('3.经校核算术平均值及计算正确');elsedisp('算术平均值及误差计算有误');endxt=sum(v(1:4))-sum(v(5:8));%判断系统误差(算得差值较小,故不存在系统误差)ifxt<0.1disp(['4.用残余误差法校核,差值为:',num2str(x1),'较小,故不存在系统误差']);elsedisp('存在系统误差');endbz=sqrt((sum(v.^2)/7));%单次测量的标准差disp(['5.单次测量的标准差',num2str(bz)]);p
4、=sort(l);%用格罗布斯准则判断粗大误差,先将测量值按大小顺序重新排列g0=2.03;%查表g(8,0.05)的值g1=(x1-p(1))/bz;g8=(p(8)-x1)/bz;%将g1与g8与g0值比较,g1和g8都小于g0,故判断暂不存在粗大误差ifg15、*sc;%算术平均值的极限误差disp(['8.算术平均值的极限误差为:',num2str(jx)]);%l1=x1+jx;%写出最后测量结果%l2=x1-jx;%写出最后测量结果disp(['9.测量结果为:(',num2str(x1),'±',num2str(jx),')']);实验二测量不确定度二、实验内容1.由分度值为0 .01mm的测微仪重复6次测量直径D和高度h,测得数据如下:/mm8.0758.0858.0958.0858.0808.060/mm8.1058.1158.1158.1108.1158.6、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);%高度平均值disp(['2.高度平均值为:',num2str(h)]);V=pi*D*D*h/4;%体积测量结果估计值disp(['3.体积测量结果估计值为:7、',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;%高度测量重复性引起的不确定度分量disp(['7.高度测量重复性引起的不确定度分量为:',num8、2str(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=(
5、*sc;%算术平均值的极限误差disp(['8.算术平均值的极限误差为:',num2str(jx)]);%l1=x1+jx;%写出最后测量结果%l2=x1-jx;%写出最后测量结果disp(['9.测量结果为:(',num2str(x1),'±',num2str(jx),')']);实验二测量不确定度二、实验内容1.由分度值为0 .01mm的测微仪重复6次测量直径D和高度h,测得数据如下:/mm8.0758.0858.0958.0858.0808.060/mm8.1058.1158.1158.1108.1158.
6、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);%高度平均值disp(['2.高度平均值为:',num2str(h)]);V=pi*D*D*h/4;%体积测量结果估计值disp(['3.体积测量结果估计值为:
7、',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;%高度测量重复性引起的不确定度分量disp(['7.高度测量重复性引起的不确定度分量为:',num
8、2str(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=(
此文档下载收益归作者所有