第四节 直接测量的数据处理
在物理实验教学中,直接测量数据处理的一般程序是:
(1)以测量列x1,x2,…,xn为样本;
(2)求出样本均值作为测量结果的最佳值;
(3)根据样本进行A类评定求出uA;
(4)根据测量仪器的性能进行B类评定,求出uB;
(5)再求出标准不确定度;
(6)最后求出扩展不确定度U=2u,写出结果表达式。
一、测量结果的最佳估计值
对被测量进行直接测量时,通常是对被测量进行n次测量,得到一个测量列
x1,x2,…,xn
它的样本均值
(2-5)
就是测量结果的最佳值。
二、不确定度A类评定
n个测量值中,任一测量值xi的大小都具有随机性,即它分布在某一个区间内的任何一点都是可能的,由贝塞尔公式给出实验标准差
(2-6)
表明了任一测量值xi的分散性。样本均值是通过随机取出的样本计算出来的,所以也有分散性,比如对同一被测量进行若干组重复测量,每一组都得到n个测量值,则每一组都能计算出一个平均值这些平均值也是各不相同的,所以也有分散性。只不过的分散性比S(xi)更小,平均值的分散性的大小用平均值的标准差
(2-7)
来表征。
由于实验中测量次数较少(一般5~8次),的分散性还要大一些,还需要乘以t分布因子tn(p),其大小与测量次数n和置信概率p有关(参见阅读材料)。于是用统计分析方法评定出的关于的A类标准不确定度为
uA=tn(p)·S (2-8)
uA的大小反映了测量结果由重复性引起的不确定度分量。
表2-2给出置信概率为0.683时tn(p)随测量次数变化的关系。
表2-2 tn(0.683)与n对应关系表
三、不确定度B类评定
下面讨论由系统效应引起的不确定度分量的评定方法。在多数情况下,由系统效应引起的不确定度分量需采用不同于A类评定的其他方法来评定,称为B类评定。
在进行B类评定时,首要问题是要知道测量仪器的“最大允许误差”。所谓最大允许误差,是指对给定的测量仪器,有关规范、规程允许的误差极限值,本书以Δ仪来表示。仪器的型号不同,其最大允许误差也不同。有些仪器可以通过查询国家计量检定规程而得到,如卡尺、千分尺、天平等。有些仪器可以在其铭牌和使用说明书中查到,如直流电桥、直流电位差计等。还有些仪器,在铭牌上给出了准确度等级,它可以换算成Δ仪。总之,在进行B类评定时,要通过查阅相关资料以获得测量仪器的性能参数。
仪器的最大允许误差,是仪器示值相对测量结果之差的最大允许值,也可以理解为这个差值的分布不会超过(-Δ仪,+Δ仪)这一区间,或者说在不考虑随机误差时被测量值以100%的置信概率包含在这一区间之内。因为uB的置信概率小于Δ仪的置信概率,所以在量值上uB<Δ仪,本书约定
(2-9)
四、标准不确定度和扩展不确定度的计算
uA是由A类评定得到的标准不确定度分量,uB是由B类评定得到的标准不确定度分量,将uA和uB按方和根合成,得到的标准不确定度
(2-10)
将式(2-8)、式(2-9)代入上式得
(2-11)
这样,的扩展不确定度为
U=2u (2-12)
五、测量结果表达式
对直接测量来说,测量结果为,扩展不确定度为U,则测量结果表达式为
或
(2-13)
下面对式(2-13)做几点说明
(1)式(2-13)说明被测物理量的量值以较大概率落于区间之内。当k=2时,这一概率可望达到95%以上。或者反过来说,被测量的量值落在上述区间之外的可能性较小,其概率不会超过5%。
(2)测量结果表达式应附有被测物理量的单位。
(3)关于和S的计算,应借助带有统计运算功能的袖珍计算器,它可以使繁重的计算变得轻松和快捷。袖珍计算器的操作方法见附录1。
(4)关于有效数字取位的原则:扩展不确定度U只取1位有效数字,第2位按四舍五入处理。当U的首位数为1或2时,可取2位有效数字,相对扩展不确定度应取2位有效数字。平均值的有效数字首先按运算规则确定,结果表达时的末位要与U所在位对齐(有效数字末位是不确定度所在位),最末一位按“四舍六入五凑偶”原则舍入。
六、直接测量的数据处理举例
【例2-1】 用Δ仪=0.02mm的游标卡尺测某物长度,测量数据(单位:mm)是
29.18,29.24,29.28,29.26,29.22,29.24
求:(1)样本均值;(2)单次测量值的实验标准差S(xi);(3)平均值的标准差S;(4)A类评定的不确定度分量uA;(5)B类评定的不确定度分量uB;(6)扩展不确定度U;(7)写出测量结果表达式。
解 将测量值逐一输入计算器,经统计运算后,屏上分别显示出
(1):29.23666667
(2)S(xi):0.0344448134
(3)S():0.014063392
(4)A类分量
1.11×0.014063392=0.0156mm
(5)B类分量
=0.0115mm
(6)扩展不确定度
(7)下面对U和进行有效数字修约,U的首位数是3,可取1位有效数字,所以U=0.04,末位与U所在位对齐,即的末位应在百分位上,因为U就是在百分位上。这样结果表达式应写成
x=(29.24±0.04)mm, k=2
也可以写成x=29.24(1±0.14%)mm, k=2
【例2-2】 改正下面给出的结果表达式错误。
(1)x=17000±1000m;
(2)x=10.8000±0.25cm;
(3)x=567.84±8.47(s);
(4)x=3.29×103±86.3(cm3)
解 (1)x=(1.70±0.10)×104m;
(2)x=(10.80±0.25)cm;
(3)x=(568±8)s;
(4)x=(3.29±0.09)×103cm3
【例2-3】 某一数字多用表,最大允许误差为Δ仪=0.005%×读数+3×最小步进值。
用此仪表测高值电阻共测量10次,数据如下(单位kΩ)
999.31,999.41,999.59,999.26,999.54,999.23,999.14,999.06,999.92,999.62
试写出测量结果表达式。
解 测量结果
实验标准差
平均值的标准差
不确定度分量A类评定
由系统效应(仪表准确性)引入的不确定度分量B类评定
测量结果的标准不确定度
扩展不确定度为
U=2u=0.196=0.20kΩ
因为U的首位数是1,所以取两位有效数字。结果表达式为
R=(999.41±0.20)kΩ, k=2
相对扩展不确定度 也可以写成
R=999.41(1±0.020%)kΩ, k=2