1、误差理论与数据处理实验报告误差理论与数据处理实 验 指 导 书 姓名 学号 机械工程学院2016年05月实验一 误差的基本性质与处理一、实验内容1对某一轴径等精度测量8次,得到下表数据,求测量结果。序号(10-4)1234567824.67424.67524.67324.67624.67124.67824.67224.674-0.0001 0.0009 -0.0011 0.0019 -0.0031 0.0039 -0.0021 -0.00010.0002 0.0077 0.0127 0.0352 0.0977 0.1502 0.0452 0.0002Matlab程序:l=24.674,24.6
2、75,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,bh0,故以上计算正确if bh0 disp(3.经校核算术平均值及计算正确);else disp(算术平均值及误差计算有误);endxt
3、=sum(v(1:4)-sum(v(5:8);%判断系统误差(算得差值较小,故不存在系统误差)if xt0.1 disp(4.用残余误差法校核,差值为:,num2str(x1),较小,故不存在系统误差);else disp(存在系统误差);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都小于g0,故判断
4、暂不存在粗大误差if g1g0&g8=max(Falpha) table2,7=高度显著;elseif (Fb=min(Falpha) table2,7=显著;else table2,7=不显著;endtable3、程序x=12 13 14 15 16 18 20 22 24 26;y=52.0 55.0 58.0 61.0 65.0 70.0 75.0 80.0 85.0 91.0;plot(x,y,*k)title(散点图);X=ones(size(x), x;b= regress(y,X,0.05);disp(y随x变化的经验公式为:y=,num2str(b(1),+,num2str(b(2),x)