测量不确定度的发展历史
2019-09-24 15:11 文章来源:检验医学网
由于测量不确定度及其评定不仅适用于计量领域,还适用于一切与测量有关的其它领域,因此1986年国际计量委员会要求国际计量局(BIPM)、国际电工委员会(IEC)、国际标准化组织(ISO)、国际法制计量组织(OIML)、国际理论和应用物理联合会(IUPAP)、国际理论和应用化学联合委员会(IUPAC)以及国际临床化学联合会(IFCCC)等七个国际组织成立了专门的工作组,起草了测量不确定度评定的指导性文件。经过工作组近七年的讨论,由ISO第四技术顾问组的第三工作组负责起草,并于1993年以七个国际组织名义联合发布了《测量不确定度表示指南》(Guide to the Expression of Uncertainty in Measurement,以下简称GUM)和第二版《国际通用计量学基本术语》International Vocabulary of Basic and General Terms in Metrology,以下简称VIM)。1995年又发布了GUM的修订版。这两个文件为在全世界统一采用测量结果的不确定度评定和表示奠定了基础。
除上述七个国际组织外,国际实验室认可合作组织(ILAC)也表示承认GUM。这就是说,在各国的实验室认可工作中,无论是校准实验室还是检测实验室,在进行测量结果的不确定度评定时均应以GUM为基础。这也表明GUM和VIM这两个文件的权威性。
GUM对所采用术语的定义和概念、测量不确定度的评定方法以及不确定度报告的表示方法都作了明确的统一规定。它使不同的国家和地区,以及不同的测量领域在表示测量结果及其不确定度时,具有相同的含义。
1998年我国发布了国家计量技术规范JJF1001-1998《通用计量术语及定义》,其中前六章的内容与VIM的第二版完全对应。
测量不确定度的概念以及不确定度的评定和表示方法的采用,是计量科学的一个新进展。
1、测量不确定度评定流程
测量不确定度评审的流程如图1所示。
图1 测量不确定度的评定流程
2、不确定度评定步骤
当被测量确定后,测量结果的不确定度仅仅和测量方法有关。因此,在进行不确定度评定之前必须首先确定被测量和测量方法。这里的测量方法包括测量原理、测量仪器、测量条件、测量程序以及数据处理程序等。确定这些条件后,测量不确定度评定步骤如下:
第1步:找出所有影响测量不确定度的来源
首先要找出所有对测量结果有影响的各种影响分量,即所有的测量不确定度来源。原则上,测量不确定度来源既不能遗漏,也不要重复计算,特别是对于比较大的不确定度分量。
测量过程中的随机效应和系统效应均会导致测量不确定度,数据处理中的修正也会导致不确定度。注意,这是从产生不确定度原因上所作的分类,与根据评定方法上所作的A、B分类之间不存在任何联系。
第2步:建立满足测量不确定度评定所需的数学模型
建立数学模型也称为测量模型化。目的在于要建立满足测量所要求准确度的数学模型,即被测量Y和所有各影响量Xi之间的具体函数关系,其一般形式可写成: (1)
影响量Xi也称为输入量,被测量Y也称为输出量。从原则上说,数学模型应该就是用以计算测量结果的计算公式。但多数情况下,公式都经过了一定程度的近似和简化,有些因素对测量结果影响很小,因此在计算时可以忽略,但对于测量不确定度来说可能是必须考虑的,因此数学模型和计算公式经常是有差别的。
数学模型并不是一成不变的。对于同样的被测量和同样的测量方法,当所要求的测量准确度不同时,需要考虑的不确定度分量数目可能不一样,此时的数学模型也可能会有所差别。
第3步:确定各输入量的估计值以及相应对应的标准不确定度u(xi)
测量结果是由各输入量的最佳估计值代入计算公式或数学模型后得到的,因此各输入量的不确定度显然会对测量结果的不确定度有影响。确定各输入量的最佳估计值的确定大体上分成两类:通过试验测量得到,或由诸如检定证书、校准证书、材料手册、文献资料以及实践经验等其它各种信息来源得到。那么对应的不确定度的评定方法也不同,即不确定度的A类评定和B类评定。A类评定是指通过一组观测列 进行统计分析,并以试验标准差来表征其标准不确定度的方法。与A类评定方法有所不同的其它方法均称为B类评定。它们是基于经验或其它信息的假定概率分布估算的,也用标准差表征。
第4步:确定对应于各输入量的标准不确定度分量ui(y)
若输入量估计值xi的标准不确定度u(xi), 则对应于输入量的标准不确定度分量为ui(y)为 (2)
(2)式中,ci称为灵敏度系数,它可以由数学模型对输入量xi求偏导数得到。当无法找到可靠的数学表达式时,灵敏度系数ci也可以由实验测量得到。在数值上它等于当输入量xi变化一个单位量时,被测量y的变化量。
第5步:列出不确定度分量汇总表
不确定度分量汇总表也称为不确定度概算。当然这一步并非是必须的。但经验表明,列出不确定度分量汇总表有利于对不确定度评定进行分析、检查、比较和交流。具有一目了然的效果。可以看出那些不确定度分量对测量结果起主要作用。如果最后得出的不确定度不满足要求,或者远远超过要求,可以结合成本考虑,对其主要影响分量进行合理的调整与改进,以能满足要求为目的。
第6步:将各标准不确定度分量u(xi)合成得到合成标准不确定度uc(y)
根据方差合成定理,当数学模型为线性模型,并且各输入量xi彼此独立无关时,合成标准不确定度为:
(3)式通常称为不确定度传播定律。不确定度传播定律实际上是将数学模型按泰勒级数展开后,对等式两边求方差得到的。对于线性数学模型,由于泰勒级数中二阶以上的偏导数项均为零,于是获得上述(3)式。当数学模型为非线性时,原则上(3)式不成立,而应考虑高阶项。但若非线性不很明显,通常有高阶项远小于一阶项,(3)式近似成立。若非线性很强时,则应考虑高阶项。
第6步:确定被测量Y可能值分布的包含因子k
根据被测量Y分布情况的不同,所要求的置信概率p,以及对测量不确定度评定具体要求的不同,分别采用不同的方法来确定包含因子k。得到各分量的标准不确定度后,应该先对被测量Y的分布进行估计。当被测量Y接近正态分布时,并且要求给出对应的置信概率为p的扩展不确定度Up时,需计算各分量的自由度和对应于被测量Y的有效自由度。并由有效自由度 和所要求的置信概率p查t分布表得到k值。如果确信自由度 不太小,例如15以上,可以直接取k=2。此时给出的扩展不确定度用U表示。
当无法判断被测量Y接近于何种分布时,一般直接取k=2,扩展不确定度用U表示。
第7步:确定扩展不确定度U或Up
扩展不确定度U=kuc。当包含因子k由被测量的分布以及规定的置信概率p得到时,扩展不确定度用Up=kpuc 表示。
第8步:给出测量不确定度报告。
报告中应给出尽可能多的信息,避免用户对所给不确定度产生错误的理解,以致错误地使用所给出的测量结果。报告中测量结果及其不确定度的表达方式应符合JJF1059-1999的规定,同时应注意测量结果及其不确定度的有效数字位数。