2025年测量不确定度与数据分析实验研究报告.docxVIP

S海工程技术大学

ShanghaiUniversityofEngineeringScience

《误差理论与数据处理》

实验指导书

姓名

学号

机械工程学院05月

试验一误差的基本性质与处理

一、试验内容

1.对某一轴径等精度测量8次，得到下表数据，求测量成果。

序号

l,/mm

v;/mm

v2/mm2(10?)

1

24.674

2

24.675

3

24.673

-0.0001

0.00020.0077

4

24.676

0.0009-0.0011

0.0019-0.0031

0.01270.0352

5

24.671

0.0039-0.0021

0.09770.1502

6

24.678

0.04520.0002

7

24.672

-0.0001

8

24.674

Matlab程序：

l=[24.674,24.675,24.673,24.676,24.671,24.678,24.672,24.674];%已知测量值

xl=mean(1);%用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,故以上计算

对的ifbh0

disp(3.经校核算术平均值及计算对的);else

disp(算术平均值及误差计算有误);end

xt=sum(v(1:4))-sum(v(5:8));%判断系统误差(算得差值较小，故不存在系统误差)ifxt0.1

disp([4.用残存误差法校核，差值为：,num2str(x1),较小，故不存在系统误差J);else

disp(存在系统误差);end

bz=sqrt(sum(v.^2)/7);%单次测量的原则差disp([5.单次测量的原则差，num2str(bz)]);

命令行窗口

disp([5.单次测量的标准差’,num2str(bz)]):

p=sort(1):%用格罗布斯准则判断粗大误差，先将测量值按大小顺序重新排列

g0=2.03:%查表g(8,0.05)的值g1=(x1-p(1))/bz;

g8=(p(8)-x1)/bz;%将g1与g8与g0值比较，g1和g8都小于g0,故判断暂不存在粗大误差ifg1gOg8g0

disp(6.用格罗布斯准则判断，不存在粗大误差):

end

sc=bz/(sqrt(8));%算数平均值的标准差

disp(7.算术平均值的标准差为：,num2str(sc))):t=2.36:%查表t(7,0.05)值

jx=t*sc:%算术平均值的极限误差

disp([8.算术平均值的极限误差为：,num2str(jx)]):%11=x1+jx:%写出最后测量结果

%12=x1-jx:%写出最后测量结果

disp([9.测量结果为：,num2str(x1),±,num2str(jx),)]):1.算术平均值为：24.6741

2.残余误差为：-0.0001250.00875-0.0011250.001875-0.0031250.003875-0.002125-0.000125

3.经校核算术平均值及计算正确

4.用残余误差法校核，差值为：24.6741较小，故不存在系统误差5.单次测量的标准差0.0022321

6.用格罗布斯准则判断，不存在粗大误差7.算术平均值的标准差为：08.算术平均值的极限误差为：0.00186249.测量结果为：(24.6741±0.0018624)

试验二测量不确定度

二、试验内容

1.由分度值为0.01mm的测微仪反复6次测量直径D和高度h,测得数据如下：

D,/mm

8.075

8.085

8.095

8.085

8.080

8.060

h,/mm

8.105

8.115

8.115

8.110

8.115

8.110

请按测量不确定度的一般计算环节，用自己熟悉的语言编程完毕不确定度分析。MATLAB程序及分析如下：

A=[8.0758.0858.0958.0858.0808.060];

B=[8.10