欢迎来到天天文库
浏览记录
ID:47074401
大小:627.69 KB
页数:12页
时间:2019-07-17
《误差理论与大数据处理 实验报告材料》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、《误差理论与数据处理》实验指导书姓名学号机械工程学院2016年05月11实验一误差的基本性质与处理一、实验内容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程
2、序:l=[24.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,
3、bh<0,故以上计算正确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.单次测量的标准差',n
4、um2str(bz)]);p=sort(l);%用格罗布斯准则判断粗大误差,先将测量值按大小顺序重新排列g0=2.03;%查表g(8,0.05)的值11g1=(x1-p(1))/bz;g8=(p(8)-x1)/bz;%将g1与g8与g0值比较,g1和g8都小于g0,故判断暂不存在粗大误差ifg15、36;%查表t(7,0.05)值jx=t*sc;%算术平均值的极限误差disp(['8.算术平均值的极限误差为:',num2str(jx)]);%l1=x1+jx;%写出最后测量结果%l2=x1-jx;%写出最后测量结果disp(['9.测量结果为:(',num2str(x1),'±',num2str(jx),')']);11实验二测量不确定度二、实验内容1.由分度值为0 .01mm的测微仪重复6次测量直径D和高度h,测得数据如下:/mm8.0758.0858.0958.0858.0808.060/mm86、.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);%高度平均值disp(['2.高度平均值为:',num2str(h)]);V=pi*D*D*h/4;7、%体积测量结果估计值disp(['3.体积测量结果估计值为:',num2str(V)]);s1=std(A);%直径标准差disp(['4.直径标准差为:',num2str(s1)]);u1=pi*D*h*s1/2;%直径测量重复性引起的不确定度分量11disp(['5.直径测量重复性引起的不确定度分量为:',num2str(u1)]);v1=5;%自由度s2=std(B);%高度标准差disp(['6.高度标准差为:',num2str(s2)]);u2=pi*D*D*s2/4;%高度测量重复性引起的不确8、定度分量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时对
5、36;%查表t(7,0.05)值jx=t*sc;%算术平均值的极限误差disp(['8.算术平均值的极限误差为:',num2str(jx)]);%l1=x1+jx;%写出最后测量结果%l2=x1-jx;%写出最后测量结果disp(['9.测量结果为:(',num2str(x1),'±',num2str(jx),')']);11实验二测量不确定度二、实验内容1.由分度值为0 .01mm的测微仪重复6次测量直径D和高度h,测得数据如下:/mm8.0758.0858.0958.0858.0808.060/mm8
6、.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);%高度平均值disp(['2.高度平均值为:',num2str(h)]);V=pi*D*D*h/4;
7、%体积测量结果估计值disp(['3.体积测量结果估计值为:',num2str(V)]);s1=std(A);%直径标准差disp(['4.直径标准差为:',num2str(s1)]);u1=pi*D*h*s1/2;%直径测量重复性引起的不确定度分量11disp(['5.直径测量重复性引起的不确定度分量为:',num2str(u1)]);v1=5;%自由度s2=std(B);%高度标准差disp(['6.高度标准差为:',num2str(s2)]);u2=pi*D*D*s2/4;%高度测量重复性引起的不确
8、定度分量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时对
此文档下载收益归作者所有