标准测量不确定度的B类评定李慎安 测量不确定度表述讲座 第六讲 标准测量不确定度的B类评定 □国家质量技术监督局 李慎安 61 什么叫不确定度的B类评定? 测量不确定度的评定方法主要分成两大类一类是用统计方法进行评定,称之为A类评定(参阅12),而其他的非统计方法,统称之为B类评定,又称之为非统计方法的评定,由此评定出来的不确定度一般称为B类不确定度或称为B类标准不确定度要注意的是INC—1(1980)(参阅12)中以及《JJF1027》中都曾规定A类不确定度分量用符号si,而B类不确定度分量用符号uj表示,这一方式在《导则》以及《JJF1059》中已作了更改,s只是实验标准偏差的符号,当它作为不确定度时,则不论是A类还是B类方法所得到,一律用u作为符号而只以数字序号作为下标相区别,一般则写作ui从量值上说s=u,但含义不同 B类不确定度只是其评定方法与A类不同,如此而已在合成过程以及对测量不确定度的贡献中完全一样,它们都以标准偏差给出,也都可以评定其自由度 62 用于评定B类不确定度的信息一般有哪一些? 由于B类评定方法不是按统计方法进行的,一般不需要对被测量在统计控制状态下(或是重复性条件下或复现性条件下)进行重复观测,而是按现有信息加以评定。
所用信息一般有: (1)以前的观测数据 例如,对某一型号的测量仪器的重复性(参阅《JJF1001—1998》第7.27条)按A类评定方法,重复了20次观测,得出了其单次示值的分散性,即重复性标准偏差sr由于这个sr的自由度υ=20-1=19,一般来说,也是充分可靠的了所以,这个数据可以用作该测量仪器进行一次、或重复几次测量结果的不确定度评定信息 但是,同一型号的某测量仪器的重复性如果彼此并不一定接近,例如,有1/3或1/4左右的差,那么,如果我们对例如只有三台这样的测量仪器,分别各进行了20次观测试验,并分别得出它们的重复性分别为sr1=34,sr2=24,sr3=29这三个标准偏差可认为充分可靠,而这样的差别则是反映了仪器本身重复性的不同这一现象,说明不能用一台仪器的sr代替同类型的其他仪器的sr使用如果,如上例,把这三个sr取平均值sr=29,那么其标准偏差按贝塞尔公式 s(sr)= = 008 =0.29自由度为2,因此,拿2.9作为其他同类型测量仪器的单次测量的分散性标准差也是不行的这个例子说明:对以前的观测数据应加以分析,看其是否可用于当前的测量结果。
(2)对有关技术资料和测量仪器特性的了解和经验 例如50 mm的量块,其中心长度最大允许示值误差,对于零级来说是±025μm,一级是±050μm、二级是±100μm仅是这样的信息是不够的,还应了解在二级量块中是不会出现示值误差在±050μm内的量块;同样,在一级量块中,也决不会出现示值误差在±025μm内的量块原因是在成批生产出同一标称尺寸的量块后,按最大允许示值误差把±025μm者作为零级挑出,把从±026μm至±050μm者挑出作为一级,如此类推因而虽然有最大允许示值误差的信息,按上述特性,只是接近两点分布,有了这一信息,就可以评定其标准不确定度了 (3)生产部门提供的技术说明文件 例如对某些装备了玻璃短标尺的测长光学仪器,一般,生产厂是给出了该尺的线膨胀系数的,我们就可以依据它对由于温度测量导致的不确定度进行评定 (4)校准证书、检定证书或其他文件提供的数据、准确度的等别或级别,包括目前还暂在使用的极限误差等 校准是指在规定条件下,为确定测量仪器或测量系统所指示的量值,或实物量具或参考物质所代表的量值,与对应的由标准所复现的量值之间关系的一组操作简而言之,即为确定测量仪器示值误差的一组操作,或给出示值的一组操作。
例如对标准硬度块的定度,对标准砝码、量块、标准电池的赋值等当然,在校准证书中必定要给出其校准结果的不确定度,甚至还要给出其自由度给出该测量仪器所属的等别(order)或(和)级别(class),这些都是据以评定标准不确定度的极为常见的信息 至于极限误差,目前还在JJG2009《射频与微波功率计量器具检定系统》、JJG2010《射频与微波衰减计量器具检定系统》、JJG2016《粘度计量器具检定系统》、JJG2027《磁感应强度计量器具检定系统》等为数尚不少的技术规范中使用,这是我国50年代以来的作法,由于对极限误差的评定缺乏统一的规范,导致其含义不确切当前不应再使用这一概念来表述测量结果,但对早已存在的规范,未修订以前,还应可以作为评定标准不确定度的依据 (5)手册或某些资料给出的参考数据及其不确定度 例如手册上给出的基本物理常量、阿伏加德罗常数L=(6022 136 7±0000 003 6)×1023 mol-1,并声明±号后之值为标准偏差,其自由度为17 又如国际上1992年公布的相对原子质量例如对于碳原子Ar(C)=12011(1),括弧中的1表示相对原子质量值的标准偏差为0001。
可惜的是没有指明不确定度的自由度但这可用于对不确度分量的评定 (6)技术规范中对某些测量方法所规定的重复性限r或复现性限R 在一些检测方法的国家标准或其他技术规范中给出的r或R一方面用于查明实验过程是否处于所规定的状态,是否出现过大的误差,另一方面,说明了这一方法的不确定度 B类标准不确定度u(qi)的二次方,可简称为B类方差 63 当已知扩展不确定度U及包含因子k时,如何评定其标准不确定度? 例如当前对圆锥量规、锥度仪器、热电偶、维氏硬度计、色温度仪器、发光强度计、核素活度计等类,在其检定系统中均明确规定校准结果要按k=3给出扩展不确定度U而例如超声功率计、橡胶硬度计、黑白密度计等,采用k=2给出扩展不确定度U在这些情况下,只要用U除以k即可得出标准不确定度例如:某交流数字电流表,流量上限为10A,相对扩展不确定度为5×10-4,k=3,则其相对标准不确定度urel=5×10-4/3=17×10-4 64 如果检定、校准证书上给出总体标准偏差σ的倍数时,如何评定标准偏差? 当前还有不少测量仪器仍沿用50年代以来的习惯,用总体标准偏差σ及其倍数表示校准结果的可靠程度。
例如:真空测量仪器、燃烧热测量仪器、电导流量仪器、水流量测量仪器、某些压力测量仪器等本来σ只是一种理想概念,在计量学中的总体指重复性条件下,对同一被测量进行了无限多次的重复观测,通过这无限多个观测结果按下式: σ= ∑(qi-μ)2/N计算出来的,式中μ为对被测量Q进行无限多次的算术平均值,称为总体均值,N是重复观测次数,为无限大由于实际测量中的重复次数n是极为有限的,按贝塞尔公式算出的实验标准偏差s只是σ的一个估计,应该说,把s等同于σ是不合理的,s本身还有不确定度,而σ的不确定度为零但是,现在在一些测量仪器既已用了σ的情况下,我们在评定标准不确定度时,可以把σ直接作为标准不确定度u处理,即从量值上说σ=u又如,射频与微波功率计、脉冲参数计量仪器、真空测量仪器等,它们用3σ给出时,采用其三分之一即获得标准不确定度对于那些燃烧热、气体流量、石油螺纹测量仪器等,它们用2σ给出时,采用其一半即获得相应标准不确定度 对于例如液体闪烁放射性活度测量仪器、质量测量仪器等,给出置信概率为9973%的情况下,由于原假设为理想的正态分布,而且所得到的合成标准不确定度又十分可靠,在这一情况下给出了上述置信概率,尽管从实际上说不十分合理,虽原为3σ的含义在进行不确定度评定时,可按U99处理,即除以26,但也未必不可以除3,前者偏保守即稍可靠一些。
作这种评定时,自由度可估计大一些,而用除以3来评定u时,自由度就不能估计得太大 65 已知Up的情况下,对标准不确定度的评定如何? 从现在起,将会越来越多地在校准证书上给出校准结果的扩展不确定度U95,或U99,这样就十分明确地交代了置信概率p的大小如果没有特殊说明,一般总按正态分布考虑,并据以评定其标准不确定度在正态分布情况下,置信概率p与包含因子kp之间的对应关系为: p=090 kp=164; p=095 kp=196; p=099 kp=258例如肖氏硬度计用U95,而(27315~90315)K温度计用U99显然,如给出的是U95,则u=U95/196,如给出的是U99,则u=U99/258 66 如果校准证书上既给出了Up,又给出了其自由度υ,应如何评定其标准不确定度u? 在《JJF1059》中,推荐在给出扩展不确定度U99或U95时,同时给出其自由度当这些信息已在校准证书中给出后,使用者不仅可以更有把握地给出标准不确定度u,而且这个u之值会比65所评定出的要小一些,对不确定度评定是有利的 当有了自由度υ以后,可以按以下t值在表格中查出kp=tp(υ)这个值。
例如:已知的自由度υ=10,采用的置信概率p=099,表中给出kp=317,这个值比上述65中的258大了不少,约五分之一这样,通过Up/kp所得的标准不确定度就会比用k=258除得出的值要小约五分之一油水不少! 67 当已知量X之值x分散区间的半宽为a,且x落在x-a至x+a区间的概率p为100%,即全部落在此范围内,如何评定标准不确定度u(x)? 在x±a内包含了X的全部可能值a处于x的两侧,而且,x处于已知区间x-a至x+a的中央拿a作为扩展不确定度U(k=3)或是U99未必不可以,只是不太确切,也不很可靠p=100%可认为接近99%,而k=3则不一定确切在这一情况下要评定其标准不确定度u(x)时,与X可能值的分布类型关系很大,因而必须对其分布作出一个近似的估计,需要有一定的经验分布的情况与包含因子k的关系是:越接近(向正态分布趋近)正态分布,k值越大从两点分布到正态分布,k值由1增加到3 接近正态分布: k=3 三角分布: k= 6 ≈245 梯形分布: k=2 矩形分布: k= 3 ≈173 反正弦分布: k= 2 ≈141 两点分布: k=1以上梯形分布指其上底与下底之比β=071的较为多见的标准状态。
当包含因子k增大,评定出的标准不确定度u(x)=a/k就减小 分散区间的半宽a可以理解为可能出现的误差极限值之模即绝对值最大的误差,一般用误差限这个概念但表述为误差限时,带有正负号(±),例如测量仪器的最大允许误差 一般用a-,表示小于x的那个下界,而用a+表示大于x的那个上界,或称为x的下限、上限值 68 如何对分散区间x-a至x+a范围内,x可能值的分布类别进行估计? 如果测量结果x是几个观测值的平均值(例如三个以上)(经常用算术平均值作为测量结果),则其分布必定是正态的如果被测量Y是由两三个以上的分量按线性合成时(例如相加),而这些分量彼此的大小又比较接近,则y的分布是接近正态的,如果Y是由两个分量线性合成,而这两个分量接近,且是三角分布,则y的分布也会是正态的自由度的大小与分布无关,不能作为考虑分布的依据 在缺乏任何其他信息的情况下,一般估计为矩形分布(均匀分布。