JJF1059_1999_测量不确定度
U NCERTAINT Y OF MEASUREM ENT 测量不确定度
JJF1059) 19995测量不确定度评定与表示6讨论之六
测量仪器不确定度
李慎安
(国家质量监督检验检疫总局, 北京 100088)
[摘 要]讨论了测量仪器不确定度这一术语的不同定义及其一般评定的方法。[关键词]测量仪器; 测量不确定度[中图分类号]TB9
[文献标识码]B
[文章编号]1002-1183(2006) 04-0039-03
1 测量仪器不确定度的概念
从JJF1001) 19985通用计量术语及定义6中, 以及当前国际上所普遍认同的不确定度定义中, 明确地表明了/不确定度0是与测量结果相联系的一个参数。它和被测量的测量结果(最佳估计值) 一起, 说明被测量之值可能出现的一个区间或说明所获得的最佳估计值的可能出现的误差限。总之, 它只与测量结果相联系而不用于测量仪器。因而可以认为:不存在测量仪器不确定度这一概念。
2002年末, 国内一些专家在北京的一个研讨会上提出了:测量仪器不确定度是测量仪器示值不确定度的简称。而/示值0则可以是测量结果, 即令一些实物量具, 例如砝码、量块等, 其示值就是其标称值, 也存在不确定度。这样的解释, 可以认为存在测量仪器不确定度这一概念。
2004年, 国际上包括ISO 在内的7个国际组织对5International Vocabulary of Basic and General T erms in M etrology, 19932nd ed 6进行修订的征求意见稿中, 增加了过去所不曾列入的/测量仪器不确定度0这一术语, 并定义为:测量结果不确定度的分量中, 来源于测量仪器的分量。但正式定稿的第三版尚未见到。如果按这一概念, 其含义可理解为测量仪器示值不确定度(用标准偏差给出) 乘以它的灵敏系数(决定于所给出的量值在数学模型中的形式) 。而当灵敏系数为正1或负1时, 测量仪器不确定度就是测量仪器示值不确定度。
当前, 在第三版尚未颁布之前, 如有必要, 为了
避免混淆, 建议分别采用:测量仪器导致的不确定度和测量仪器示值不确定度这样的表述。当它们都是标准偏差时, 也可在/不确定度0前加/标准0二字。本文以下只讨论测量仪器示值标准不确定度, 简称仪器不确定度。
2 仪器不确定度一般具有的两种分量
211 系统效应导致的分量
设某个被测量Q , 在规定的条件下, 采用某测量仪器, 得到的某个示值q i 或若干个示值的平均值q , 它们都会存在一个修正值$q 。如果这一仪器经过校准, 在其校准证书上, 给出了校准值或修正值, 必然, 还应同时给出它们的扩展不确定度U , 包含因子k 或是扩展不确定度U p 。通过这样的信息, 可以很方便地计算出修正值$q 的标准不确定度u ($q ) , 而不论$q 是否等于零(参见JJF1059) 1999中512~513节) 。
如果这一仪器未经校准而只是进行了检定, 在检定证书上通常没有校准值或修正值而只给出了是否合格以及其准确度级别的结论, 那么, 这时的$q 是未知的, 当然, 也不需要知道, 也就是说, 对仪器的示值不加修正, 也不必去修正。我们说, 这一情况下, $q =0, 这个/零0当然并不可靠, 因为仪器存在最大允许误差M PE, $q 以100%的概率处于M PE 所给出的区间。按JJF1059中516节, 通过对其分布的估计, 可以得出标准不确定度u ($q ) , 通常, 按矩形分布处理。MPE 的绝对值MPEV 成为分散区间的半
[收稿日期]2006-05-10
[作者简介]李慎安(1926-) , 男, 湖北武汉人, 教授级高工, 毕业于同济大学, JJF1059) 1999起草执笔人, 为发展我国科学技术事业做出了突出贡献, 获得国务院发给的政府特殊津贴。 4#
测量不确定度 UNCERTAINTY OF M EASU REM ENT
宽a , 乘以转换系数b =016(如按JJF1059, 则是除以k =
3, 这里, b 是根据ISO 14253-2:1999的
u ($q ) =ab =016MPEV
例如:标称值m =200g 的4级(F 2级) 砝码, 按JJG99) 1990, 其最大允许误差MPE =? 310mg, 分散区间半宽a =310mg, 估计为矩形分布, 得其修正值$m =0的标准不确定度:
u ($m ) =016@310mg=118mg 。
212 随机效应导致的分量
测量仪器在重复性条件下, 对同一被测量Q 进行测量所得的某一个示值q i 或是若干次结果的平均值q 往往均存在一个重复性标准偏差s r (q i ) 或s r (q ) , 它们之间有
s r (q ) =s r (q i ) /
式中:n 为q 的重复观测次数。
例如:化学发光法氮氧化物分析仪, 在检定规程JJG 801) 2004中, 给出了仪器的任一次检测结果的重复性标准偏差为所测得的体积分数U i 的1%, 即
-2
s r (U U i ) =10i
应超过的R 之值, 视作为复重性限r 来对待。r 的定义是, 在重复性条件下, 任意两示值之差, 以95%的概率, 不会超出的值。而r 与s r (q i ) 之间有:s r (q i ) =r/218。这个式子在JJF1059中的5110节提到过。
例如:金属洛氏硬度计的示值重复性, 按JJG112) 2003, 就是对标准硬度块上, 均匀分布的5点所得示值中的最大者与最小者之差, 该规范对A 标尺(20~75) H RA 硬度范围内, 示值重复性规定不超过018HR, 则其示值的重复性标准不确定度可简单地计算:s r =018H R/218=013H R
(2) 重复性要求给出的值为若干次的重复条件下的示值与它们的算术平均值之差中的最大差值。
测量仪器示值q i 与平均值q 之差, 称之为残差v i =q i -q 。通过其中的最大残差v max 不得超出之值, 即所谓最大残差法, 用于评定s r (q i ) 也是一种简化的方法, 其缺点与极差法类似。
一个较为简单可行的办法是将规程所给出的这一指标, 往往也称为重复性, 作为(1/2) 重复性限r , 即015r 来对待。
例如:可见分光光度计的透射比S 的重复性要求, 在JJG178) 1996中, 对于棱镜型要求为015%, 实际上可以认为是015r 。
s r (S i ) =r /218=015r /(218@015) =
015%/114=0136%
(3) 同时给出多种与重复性有关的指标。例如, 拉力、压力和万能试验机, 在检定规程JJG139) 1999中, 同时给出了:示值重复性相对误差、示值进回程相对误差、零点相对误差以及示值相对分辨力。显然, 它们对重复性标准差s r 都直接有关。但是, 作为仪器不确定度中的随机效应分量, 可只考虑其中最大的那一项。例如:只有重复性标准差和分辨力带来的标准差两项, 这两项中, 哪一个的值较大就用哪一项而不必把这两项合成为随机效应导致的分量。上述JJG139中给出的4项, 实际上也可只采用其中的最大的那一项。例如:2级试验机, 所规定的示值重复性相对误差为210%, 示值进回程相对误差310%, 零点相对误差012%, 示值相对分辨力110%, 其中的前3项, 都相当于给定操作中的重复性限, 相对标准不确定度分别为:017%, 110%, 011%。而分辨力所导致的相对标准不确定度为0129@110%=
114
简化建议的计算方法) 即:
n
在现行的一些技术规范中, 对测量仪器的重复性, 存在一些不同形式的技术要求, 以下分别讨论在这样一些情况下, 如何计算s r (q i ) 。
(1) 提出在若干次的检测示值中, 极差R 不得超过的值。
极差R 的定义是, 重复性条件下若干不同示值中, 其最大者q max 减最小者q min 之差。如何通过R 计算s r (q i ) , 在JJF1059的414节中已有完整的叙述。但是, 采用极差法对s r (q i ) 进行评定, 要求所得出的这若干个q i 接近正态分布。事实上, 我们不能通过几个q i 来判断, 即令有20~30个q i 也不可能。一般最少应有50个。如要求较准确的判断, 应有200个以上。对我们来说很难。其次, 极差法所得到的s r (q i ) 可靠性较差, 即其自由度小于贝塞尔法所得出的s (q i ) 的自由度, 如果并不接近正态分布而采用了这一方法, 所得s r (q i ) 必然偏大。如较大地偏离正态, 则会偏大得很多。此外, 如这一方法中所用的重复次数较小, 例如, 少于4次, 就会太不可靠。而要得到一个较可靠的s r (q i ) , 与之相乘的安全因子h 就会很大(参见ISO14253-2:1999) 。
一个较为简单可行的办法是将规程中所给出的不
U NCERTAINT Y OF MEASUREM ENT 测量不确定度
013%。因此, 作为仪器不确定度中的随机效应导致的分量, 只采用110%即可, 其余3项均应不计。
(4) 一些检测标准中给出的重复性限, 只能作为方法确认的重复性标准差的依据而不是仪器示值的重复性标准差的依据。
例如:H G/T 2771) 96规定对同一样品中氯化镍所作的任两次定量分析结果, 其质量分数w 之差不得大于012%(以NiCl 2#6H 2O 计) 。这个重复性限r 除以218所得出的012%/218=0107%是该方法所确认的任一次分析结果的重复性标准偏差, 而不是仪器的。
来评定它所导致的标准不确定度(系统效应) 导致的分量。但如果MPE 是用相对误差给出时, 按211节, 乘以转换系数016后, 并非相对标准不确定度。
例如, 显微硬度计, 按JJG260) 1991, 在硬度在200~300HV 内范围内, 示值允许误差规定为? 510%。这里给出的值的含义是:仪器的示值与标准块给出的示值之差(即示值与标准值之差) 除以标准块硬度值, 不得超出的值。设标准硬度块之值为1H V, 最大允许误差M PE=? 510%HV 。由此而导致的系统效应分量(参见211节) :
u ($HV) =016@5%HV=3%HV
这个式子的确切含义是, 仪器示值的修正值$HV=0的标准不确定度为3%HV 。这里用标准块的值H V 还是用仪器的示值与3%相乘, 关系不大。因为它们是接近的。
但是, 我们不能认为$HV 的相对标准不确定度:
u rel ($H V) =3%
按相对标准不确定度的定义:
u rel (q ) =u (q ) /q
那么,
u rel ($H V) =u ($H V) /$HV
3 仪器不确定度中两类分量的合成
311 对于实物量具, 只有一类分量
按JJF1001) 1998中512节, 实物量具的示值就是它所标出的值, 其中, 如砝码、量块、标准电阻、量器、参考物质等。在其检定证书中, 给出了其级别或最大允许误差及是否合格的结论, 据此我们可以按JJF1059第516节评定其标准不确定度或相对标准不确定度, 作为仪器不确定度而不存在另外的分量。312 仪器不确定度的两类分量彼此独立, 且灵敏系数均为1
一般仪器对某量Q 的测量结果q i 或q , 在给出其数学模型时为:
Q =q i +$q
或是
Q =q +$q
Q 的两个输入量估计值是彼此独立的。在这一线性关系式中, q i 与$q 两个输入量的灵敏系数均为1。因此, 在评定不确定度时, 就只需要将这两个分量u (q i ) 和u ($q ) 平方相加后开方即可。这就是仪器示值不确定度。u (q i ) 或u (q ) 是随机效应导致的重复性标准偏差s r ; 而u ($q ) 是修正值$q =0的标准不确定度, 是系统效应导致的分量。本文在第2节中, 对s r 的评定, 都是据有关JJG 给出的不同技术要求计算的, 所谓B 类评定(非统计方法的评定) 。实际上, 如果在JJG 或有关仪器的技术规范和标准中, 缺乏这类信息时, 当然也可采用A 类评定(统计方法的评定) 来获得s r , 具体方法从略。
313 尽管不确定度的分量往往来源于最大允许相对误差, 但往往并非相对不确定度
在上述211节中已提到通过最大允许误差M PE
4由于$H V=0而上式不成立。
上述评定所给出的只是:u ($H V) /1HV=3%
这一点经常被忽略而导致错误, 值得注意。314 由于分量均为与示值之比或者均为与满量程(FS) 之比, 在合成时, 可直接采用方和根
例如:显微硬度计JJG260) 1991给出上述硬度范围内, 示值重复性为6%, 但其含义为仪器在标准硬度块上5个均匀分布的点所得的5个示值的极差R 除以平均示值, 本文前面已讨论过, 这种情况下的值相当于重复性限r /HV, 由此而评定出的任一次测量结果HV 的相对重复性标准偏差(也就是相对标准不确定度)
u rel (HV) =s r (HV) /HV=r /218H V=6%/218=2%
这里的u rel (H V) 和313节提到的u ($HV) /HV 之值3%可以采用方和根得出仪器相对不确定度, 因为它们都是对HV 的比值。即
3+2%=316%。
[编辑:薛 敏]
#