5测量结果的不确定度估算

1、4 大学物理实验91.5测量结果的不确定度估算不确定度的概念一般来说，真值是无法测得的，因此误差也就无法得到。我们只能通过一定的方法对测 量误差进行估计，这就需要引入不确定度的概念。不确定度是指由于测量误差的存在而对被测量值不能肯定的程度，是对被测量的真值所 处的量值范围的评定。我们在表示完整的测量结果时，除给出被测量Xo的量值(一般用被测量的算术平均值来表示)，还要同时标出测量的总不确定度&写成x=x0 A ( P = P)(1-11)式中P为置信概率，式(1-11)的含义是：区间(Xo -A,Xo十)内包含被测量X的真值的 可能T是P。为了直观地评定测量结果，也常采用相对不确定度的概念。用

2、Ur表示相对不确定度，则有Ur =X100%(1-12)Xo根据估计方法的不同，总不确定度可分为两类分量，一类是可以通过多次重复测量用统 计学方法估算出的 A类分量Aa,另一类是用非统计方法估算出的B类分量Ab。将两类分量按方和根的方法合成，就得到测量结果的总不确定度： 二 ,Aa2 +42(1-13)1.5.2 A类不确定度分量的估算A类不确定度分量是指可以用统计学方法估算的分量，一般指随机误差。具体估算的方 法如下：根据误差理论，当重复测量次数足够多时，可求得置信概率为0.95的A类不确定度分量4 =1.96sx(1-14)式中sx是算术平均值的标准偏差。但当重复测量次数较少时，随机误差不

3、再符合正态分布。这样，需对式(1-14)做一个修正。即a =t&(1-15)式中t是由测量次数决定的修正系数，它的取值与测量次数和置信概率有关。置信概率为0.95时，t与不同测量次数n之间的关系如表1-1所示。表1-1t与不同测量次数n的对应关系n23456789101520oot12.74.303.182.782.572.452.362.312.262.142.091.96根据重复测量的次数，从表 1-1中查出相应的t值，就可得到修正后的置信概率为0.95的A类不确定度分量4第1章测量误差与数据处理的基础知识51.5.3 B类不确定度分量的估算.仪器误差测量仪器和量具本身总是存在一定误差，我

4、们习惯上称之为仪器误差，用符号仪表示。它是指仪器在规定条件下使用时，所允许的误差限值。仪器误差是一个统称，对于具体的各类仪器和量具具有不同的表示方式。例如，游标卡尺和螺旋测微计的仪器误差用示值误差表示。国家标准规定，量程为 0300mm以下的游标 卡尺，其示值误差在数值上等于该尺的最小分度值。螺旋测微计分零级、一级和二级3种精度级别，通常实验室使用的为一级螺旋测微计，其示值误差随量程而异。量程为025mm的一级螺旋测微计，示值误差为M.004mm。对物理天平而言，仪器误差用指示值变动性误差来表示天平称衡结果的可靠程度，这是 由于天平调节、操作、温差、气流以及振动等原因，致使重复称衡时各次平衡位

5、置产生差异。 按规定，合格天平的示值变动性误差不应大于该天平的最小分度值。对于实验室常用的 WL-1型物理天平，其误差限值一般取为0.05g。电表及电阻箱的仪器误差用准确度等级 K表示。在规定条件下使用电表测量，其示值的 误差限为电表量程与准确度等级百分数（即K%）的乘积。可见，电表的仪器误差大小由电表准确度等级和电表量程二者决定。为可在仪器出厂说明书或仪器标牌上查到，对于精度较低的仪器，A仪可取为其最小分度值 的一半，如米尺的示值误差为0.5mm。在工业和商业用途上，仪器误差的置信概率一般为0.9S. B类不确定度分量不能用统计学方法估算的分量为 B类不确定度分量，一般指系统误差。若不特别说

6、明，则 品4= -1-16CC称为置信因子，置信概率取 0.95时，C=1.05o1.5.4 直接测量量的不确定度估算及测量结果表示1 .单次直接测量在实际测量过程中，有的被测量是随时间变化着的，我们无法对其进行重复测量，只能 进行单次测量。还有些被测量，对它们的测量精度要求不高，只要进行单次测量就可以了。在单次测量中，用单次测量值x测作为被测量的最佳估计值。 测量值的不确定度与所用测量仪器的精度、测量者的估读能力及测量条件等很多因素有关，因此它的合理估计是比较复 杂的。在一般情况下，对随机误差很小的测量，可以只估计不确定度的B类分量，用仪器误差仪作为x测的总不确定度，测量结果表示为*=冲九4

7、仪（1-17）有的测量随机误差可能比较大，此时可以估计一个误差限来作为单次测量的不确定度A例如，用0.1s分度的秒表计时，由于人的感官灵敏度的限制与技术上的不熟练，常常造成“启6 大学物理实验9动”和“停止”秒表所用的时间超过0.1s,这必然使测量误差限超过秒表的仪器误差限。这时可依据实际情况来估计误差限，如可取 = 0.2s。又如，用钢卷尺测量较长的距离，不可能 保证尺子拉直拉平，则可依实际情况取 = 5mm或更大。总之要根据测量的不同情况以及观 测者实验技巧的高低来对单次测量的总不确定度做出估计。2.多次直接测量对多次直接测量的数据 Xi,X2, -Xn进行处理的一般步骤是:n 计算被测量

8、的算术平均值 X =Z X / n ,把X作为被测量的最佳估计值。i 3求出各测量值的残差 Vi =Xi -Xn V2用贝塞尔公式求出测量列的标准偏差S=1。n1审查测量数据，如发现有异常数据，应予以舍弃。舍弃异常数据后，再重复步骤 ，直至完全剔除异常数据。求出算术平均值的标准偏差SSn,并查表1-1求出总不确定度的A类分量Aa =1Sx 。求出总不确定度 A = JKa+联表不出最后测量结果 x =X A,U r =二父100%X注：利用计算器的统计计算功能，将多次测量结果输入后，可直接求得X及S。【例1-1】有一组以cm为单位的长度测量数据为2.20, 2.25, 2.30, 2.15,

9、2.10, 2.15,2.25, 2.10, 2.20, 2.20, 2.10, 2.15, 2.25, 2.20, 2.20, 2.15, 2.25, 2.20, 2.20, 3.50,据 “3/ 准则判断其中是否有异常数据需剔除。解：因为测量次数 n10,所以可以用“ 3,准则剔除异常数据。 20 - . x X = = 2.26cm i 土 20根据贝塞尔公式，用标准偏差S代替标准误差aS = 0.3cm3S =0.9cm根据“3,准则，因为 &0 =1.24cm 3S,所以3.50应舍去。舍去3.50后，再重新计算 测量数据，得XX = = 2.19cmii19S = 0.06cm3S

10、 =0.18cm在这19个数据中，没有一个数据的残差大于此3S,所以，根据准则“ 3，这19个数据中没有异常数据。【例1-2】用量程为025mm的一级螺旋测微计（仪=0.004mm ）对一铁板的厚度进行第1章测量误差与数据处理的基础知识了 8 次重复测量，以 mm 为单位，测量数据为：3.784, 3.779, 3.786, 3.781, 3.778, 3.782,3.780, 3.778,同时读得螺旋测微计的零位为+ 0.004mm,请给出完整的测量结果。解：可求得L =3.781mmS =0.0029mm没有一个数据的残差大于3S,所以，根据准则“ 3，这8个数据中没有异常数据，不用剔除。

11、但考虑到零位修正L =3.781 -0.004 =3.777(mm)一 SSL =0.0011mmn查表1-1可知，n=8时，t=2.36,计算得：A类不确定度分量Aa =tSL =0.0025mmB类不确定度分量仪 0 004b = mm = 0.0038mmC 1.05总不确定度l =八2 A + A2b =0.005mm测量结果为L =二、l = 3.777 , 0.005 mmUr 100% =0.13%1.5.5间接测量量的不确定度估算及测量结果表示在很多实验中，我们进行的测量都是间接测量。因为间接测量量是各直接测量量的函数，所以直接测量量的误差必定会给间接测量量带来误差，这被称为误

12、差的传递。这样一来，直 接测量结果的不确定度就必然会影响到间接测量结果，这种影响的大小可以由相应的数学公 式估算出来。.间接测量的不确定度合成公式设间接测量量y是各自相互独立的直接测量量Xi,&,Xm的函数，其函数形式为y = f (Xi,X2, ，Xm )(1-18)设各直接测量量Xi,X2, Xm的测量结果分别为Xi AXi,X2iX2, - Xm 土加,则间接测量量 y的最佳估计值为y = f (xi,X2, Xm )(1-19)由于不确定度都是微小的量，相当于数学中的“增量”，因此间接测量量的不确定度的计 算公式与数学中的全微分公式基本相同。不同之处是：要用不确定度 ,等替代微分dXi

13、等;要考虑到不确定度合成的统计性质。具体做法如下。大学物理实验9首先对函数式（1-18）求全微分:dy =dxi + dx2+，+ dxm（1-20）然后用不确定度 如取，跳，&m替代dy, dxi, dx2,dxm并将等式右端进行方和根合成， 得到间接测量量的不确定度方和根合成公式：(1-21)(1-22)(1-23)用不确定度，y,:出，必，4xm替代dy, dx1, 是间接测量量的相对不确定度的方和根合成公式：改2,一，ckm后，再进行方和根合成，得到的(1-24)对于积商形式的函数，为计算方便，可先对函数式（1-18）取对数，得 TOC o 1-5 h z In y =ln f x,x

14、2,xm再对上式求全微分:dy _ f dx1 f d&f dxmy f 乂 ff注意：用式（1-21）式（1-24）估算间接测量量的不确定度时，应使式中各直接测量 量的不确定度具有相同的置信概率。作粗略的不确定度估算时，也可采用间接测量量的不确定度算术合成公式:a句=27 &1+ :网+1+ 再%21y _行及1廿+且y&1f依2 f仅mf (1-25)(1-26)用算术合成公式估算出的间接测量量的不确定度偏大。AR【例1-3】 求函数式（1） N=A+BC和（2） N =里 的不确定度传递公式，式中A、CB、C为变量。解：（1）对函数式N=A + B-C求全微分得dN =dA dB -dC

15、用不确定度代替微分，再方和根合成便得到不确定度传递公式：二手人2b , ；2CAB（2）因为函数式 N=是积商的形式，所以先对其取对数，可得CIn N =ln A In B -In C再求全微分：第1章测量误差与数据处理的基础知识dN dA dB dC一=十_一N A B C用不确定度代替微分，再方和根合成便得到相对不确定度传递公式:2.间接测量的数据处理步骤按照直接测量量的数据处理程序求出各直接测量量的结果:X1 =X1 二 Xi,X2 =X2 二 X2, Xm = Xm 二Xm将各直接测量量的最佳估计值代人函数关系式中，求得间接测量量的最佳估计值:y = f Xi,X2, Xm求出间接测量不确定度的方和根合成公式： TOC o 1-5 h z y = /伍加2 ；令m丁一1方f; 同J鬲求出间接测量值的不确定度 表示出最后测量结果y=y 土 Ay, Ur=yxl00%y【例1-4】已测得金属环的外径 Di =2.880 0.004cm ,内径D2 =3.600 0