不确定度测定汇总

1、测量不确定度评定与表示测量的目的是确定被测量值或获取测量结果.有测量必然存在测量误差,在 经典的误差理论中,由于被测量自身定义和测量手段的不完善,使得真值不可知, 造成严格意义上的测量误差不可求.而测量不确定度的大小反映着测量水平的高 低,评定测量不确定度就是评价测量结果的质量.图11识别测量不确定度的来源测量不确定度来源的识别应从分析测量过程入手,即对测量方法、测量系统 和测量程序作详细研究,为此必要时应尽可能画出测量系统原理或测量方法的方 框图和测量流程图.检测和校准结果不确定度可能来自：(1)对被测量的定义不完善；(2)实现被测量的定义的方法不理想；(3)取样的代表性不够,即被测量的样本

2、不能代表所定义的被测量；(4)对测量过程受环境影响的熟悉不全,或对环境条件的测量与限制不完善；对模拟仪器的读数存在人为偏移；(6)测量仪器的计量性能(如最大允许误差、灵敏度、鉴别力、分辨力、死区及稳 定性等)的局限性,即导致仪器的不确定度；(7)赋予计量标准的值或标准物质的值不准确；(8)引用于数据计算的常量和其它参量不准确；(9)测量方法和测量程序的近似性和假定性；(10)在外表上看来完全相同的条件下,被测量重复观测值的变化.分析时,除了定义的不确定度外,可从 测量仪器、测量环境、测量人员、测量方法等方面全面考虑,特别要注意对测量结果影响较大的不确定度来源,应尽量做 到不遗漏、不重复.2定义

3、2.1 测量误差简称误差,是指“测得的量值减去参考量值.2.2 系统测量误差简称系统误差,是指“在重复测量中保持恒定不变或按可预见的方式变化的测量误差的分量.系统测量误差的参考量值是真值,或是测量不确定度可忽略不计的测量标准 的测量值,或是约定量值.系统测量误差及其来源可以是的或未知的.对于的系统 测量误差可以采用修正来补偿.系统测量误差等于测量误差减随机测量误差.2.3 随机测量误差简称随机误差,是指“在重复测量中按不可预见的方式变化的 测量误差的分量.随机测量误差的参考量值是对同一个被测量由无穷屡次重复测量得到的平均 值.随机测量误差等于测量误差减系统测量误差.图2测量误差示意图2.4 测

4、量不确定度简称不确定度,是指“根据用到的信息,表征赋予被测量值分 散性的非负参数.测量不确定度一般由假设干分量组成.其中一些分量可根据一系列测量值的统计分布,按测量不确定度的 A类评定随机效应引起的进行评定,并用标准偏差 表征；而另一些分量那么可根据基于经验或其它信息所获得的概率密度函数,按测 量不确定度的B类评定系统效应引起的进行评定,也用标准偏差表征.2.5 标准不确定度 是“以标准偏差表示的测量不确定度.标准不确定度全称为标准测量不确定度可采用A类标准不确定度、B类标准不确定度及合成标准不确定度、相对合成标准不确定度等表示.测量不确定度的A类评定,简称A类评定,是指“对在规定测量条件下测

5、得的量值用统计分析的方法进行的测量不确定度分量的评定.测量不确定度的B类评定,简称B类评定,是指“用不同于测量不确定度 A类评定的方法进行的测量不确定度分量的评定.2.6 合成标准不确定度 全称合成标准测量不确定度,是指“由在一个测量模型中各输入量的标准测量不确定度获得的输出量的标准测量不确定度.2.7 相对标准不确定度全称相对标准测量不确定度,是指“标准不确定度除以测得值的绝对值.2.8 自由度是指“在方差的计算中,和的项数减去对和的限制数.2.9 扩展不确定度 全称扩展测量不确定度,是指“合成标准不确定度与一个大于1的数字因子的乘积.2.10 包含区间是指“基于可获信息确定的包含被测量一组

6、值的区间,被测量值以一定概率落在该区间内.包含概率 是指“在规定的包含区间内包含被测量的一组值的概率.包含因子是指“为获得扩展不确定度,对合成标准不确定度所乘的大于1的数.包含因子有时也称扩展因子,用符号 k表示.表1表示测量不确定度常用的名称及符号名称符 号说明标准不确定度u 或 u(x i)相对标准不确定度Urelrel一表示“相对的英文字母的缩写.测量不确定度的A类评定Ua 或 UA(Xi)测量不确定度的联评定Ub 或 Ub (Xi)合成标准不确定度Uc或 Uc (y)相对合成标准不确定度Ucrel 5 Ucrel( y)扩展正确定度urnLP包含概率为p的扩展不确定度相对扩展正确定度U

7、el 或 LPrel包含因子k或kpkp包含概率为p的包含因子包含概率P如,p = 95% , p = 99%.启效自由度Veffeff 表示“有效的英文字母的缩写.注： 表中A、B、c、rel、eff为正体；x、y、k、i、p、n、u、以斜体. 表中大写 竦示扩展不确定度；小写u表示标准不确定度,如：标准不确定度A类评定：Ua标准不确定度B类评定：ub合成标准不确定度,uc或uc (y)扩展或相对扩展不确定度,L® Lb Uel或LPrel2.11测量模型是指测量中涉及的所有量间的数学关系.测量模型简称模型测量模型的通用形式是方程：f (Y, X,Xn) = 0,其中测量模型中的

8、输 出量Y是被测量,其量值由测量模型中 输入量X,Xn的有关信息推导得到.在测量模型中,输入量与输出量间的函数关系又称测量函数.建立测量模型,即被测量与各输入量之间的函数关系.假设Y的测量结果为y,输入量X的估计值为Xi,那么y=f (Xi,X2,Xn).在建立模型时要注意有一些潜在的不确定度来源不能明显地呈现在上述函数关系中,它们对测量结果本身有影响,但由于缺乏必要的信息无法写出它们与被 测量的函数关系,因此在具体测量时无法定量地计算出它们对测量结果影响的大 小,在计算公式中只能将其忽略而作为不确定度处理.图3测量不确定度评定内容3标准不确定度的A类评定(分量)3.1贝塞尔公式法在重复性条件

9、下或复现性条件下对同一被测量(一个被测件)独立重复观测n次,得到n个观测值Xi (i=1,2,.,n),被测量X勺最正确估计值是n个独立测得值的算术平均值X,按公式(1-1)计算：_ 1n,、X 1 X (1-1)n i 1单个测得值Xk的实验方差s2(Xk),按公式(1-2)计算：21 n _ 2s (Xk) (Xi X) (1-2)n 1 i 1单个测得值Xk的实验标准偏差s(Xk),按公式(1-3)计算：n s(Xk) J； (X X) (1-3),n 1 i 1公式(1-3)是贝塞尔公式,自由度v为n-1.实验标准偏差s(Xk)表征了测得值x的分散性,测得重复性用s(Xk)表征被测量估

10、计值x的A类标准不确定度为：u (x) s(x) s(xk)/Vn(1-4)A类标准不确定度u (x)的自由度为实验标准偏差s(xk)的自由度,即丫=氏1.实 验标准偏差s(x)表征了被测量估计值x的分散性.3.2在标准化的常规检定、校准或检测中评定合并样本标准偏差假设对每个被测件的被测量 X在相同条件下进行n次独立测量,测得值为xixi2,疝,其平均值为无；假设有册被测件,那么有ma这样的测得值,可按公式(1-5)计算单个测得值的合成样本标准偏差sp( xk):1 imn2"s(xk) , , 八(xij x)2Vm(n 1) i 1 j 1(1-5 )式中： i 一组数,i=1,

11、2,mj 每组测量的次数,j=1,2,n.公式(1-5 )给出的sp(x0 ,其自由度为mn-1).假设对每个被测彳已分别按 n次重复测量算出了其实验标准偏差si,那么m组测得值的合并样本标准偏差 sp(xk)可按公式(1-6)计算：1 ms (xk) J- s2(1-6),m i 1当实验标准偏差si的自由度均为Vo时,公式(1-6)给出的sp(xk)的自由度为mv.假设又m个被测量X分别重复测量的次数不完全相同,设各为 n,而X的实验标准偏差s(x)的自由度为Vi,通过m个si与Vi可得sp(xk)按公式(1-7)计算：1 m二,一s (xk) p- Vis(1-7)仆1公式(1-7)给出

12、sp(xk)的自由度为VVi.由上述方法对某个被测件进行 n'次测量时,所得测量结果最正确估计值的A类标准不确定度为：u (x) s(x) s (%)/布(1-8)用这种方法可以增大评定的标准不确定度的自由度,也就提升了可信程度.3.3 预评估重复性在日常开展同一类被测件的常规检定、校准或检测工作中,如果测量系统稳 定,测得重复性无明显变化,那么可用该测量系统以与测量被测件相同的测量程序、 操作者、操作条件和地点,预先对典型的被测件的典型被测量值进行n次测量(一般n不小于10),由贝塞尔公式计算出单个测得值的实际标准偏差s(Xk),即测量重复性.在对某个被测件实际测量时可以只测量n&#

13、39;次(1&n' <n),并以n'次独立测量的算术平均值作为被测量的估计值,那么该被测量估计值由于重复性导致的A类标准不确定度按公式(1-9)计算：u(x) s(X) s(Xi)/Vn-(1-9)用这种方法评定的标准不确定度的自由度仍为v=n-10注意：当疑心被测量重复性有变化时,应及时重新测量和计算实验标准偏差s(xk) o3.4 A类评定流程4标准不确定度的E类评定(分量)4.1 B类评定的一般表示B 类评定的方法是根据有关的信息或经验,判断被测量的可能值区间X-a, X+a,假设被测量值的概率分布,根据概率分布和要求的包含概率p估计因子k,那么B类标准不确

14、定度 3可由(2-1 )式得到：u 旦(2-1 )k式中：a为被测量可能值区间的半宽度当k为扩展不确定度的倍乘因子时称包含因子,其他情况下根据概率论获得的k称置信因子4.2 峻评定(来源)通常基于诸如以下信息：(1)权威机构发布的量值；(2)有证标准物质的量值；(3)校准证书；(4)仪器的漂移；(5)经检定的测量仪器的准确度等级；(6)根据人员经验推断的极限值等.4.3 确定B类评定的区间半宽度 a(1)生产厂提供的测量仪器的最大允许误差为±4,或由手册查出所用的参考数据误差限为±4,或当测量仪器或实物量具给出准确度等级等,并经计量部门检定合格,那么评定仪器的不确定度时,可

15、能值区间的半宽度为：a =(2)校准证书提供的校准值,给出了其扩展不确定度为U,那么区间的半宽度为：a=U(3)由有关资料查得某参数的最小可能值为a-和最大值为a,最正确估计值为该区间的中点,那么区间半宽度可以用下式估计：a =( a+ a-)/2(4)必要时,可根据经验推断某量值不会超出的范围,或用实验方法来估计可能的区间.4.4 k确实定方法(1)扩展不确定度是合成标准不确定度的假设干倍时,该倍数就是包含因子k值.2得到k值(2)假设被测量值服从正态分布时,根据要求的概率查表表2正态分布情况下概率 p与k值间的关系p0.500.680.900.950.95450.990.9973k0.67

16、511.6451.96022.5763(3)假设为非正态分布时,根据要求的概率查表3得到k值表3常用非正态分布时的k值及B类标准不确定度Ub(x)分布类别P 衿kUbX三角100a/ v16梯形B =0.71 1002a/2矩形(均匀)100a/V3反正弦100a/ v12两点1001a注：表3中B为梯形的上底与下底之比, 对于梯形分布来说,k <6/(1 2),特别当B等于1时,梯形分布变为矩形分布；当 B等于0时,变为三角分布.4.5 B类评定概率分布的假设(1)被测量受许多随机影响量的影响,当它们各自的影响都很小时,不管各影响量的概率分布是什么形式,被测量的随机变化服从正态分布.如

17、证书或报告给出的不确定度是具有包含概率为0.90、0.95、0.99的扩展不确定度(即给出 以°、LU、U99),此时,除非另有说明,可按正态分布来评定B类标准不确定度.(2)当利用有关信息或经验,估计出被测量可能值区间的上限和下限,其值在区间外的可能几乎为零时,假设被测量值落在该区间内的任意值处的可能性相同,那么可假设为均匀分布(或称矩形分布、等概率分布).如数据修约、测量仪器最大允许误差或分辨力、参考数据的误差限、度盘或齿轮的回差、平衡指示器调零不准、测量仪器的滞后或摩擦效应导致的不确定度及对被测量的可能值落在区间内的情况缺乏了解等,一般假设为均匀分布.(3)当利用有关信息或经验

18、,假设被测量值落在该区间中央的可能性最大,那么假设 为三角分布.如两相同均匀分布的合成、两个独立量之和值或差值服从三角分布.(4)当利用有关信息或经验,假设落在该区间中央的可能性最小,而落在该区间上限和下限的可能性最大,那么可假设为反正弦分布(即U形分布).如度盘偏心引起的测角不确定度、正弦振动引起的位移不确定度、无线电测量中失配引起的不确定度、随时间正弦或余弦变化的温度不确定度等.(5)按级使用量块时,中央长度偏差的概率分布可假设为两点分布.(6)安装或调整测量仪器的水平或垂直状态导致的不确定度常假设为投影分布.(7)实际工作中,可依据同行共识确定概率分布.4.6 分辨力导致的B类不确定度分

19、量假设数字显示器的分辨力为Sx,由分辨力导致的标准不确定度分量 u(x)采用B类评定,那么区间半宽度为 a=Sx/2,假设可能值在区间内为均匀分布,查表 3得k =73,因此由分辨力导致的标准不确定度分量u(x)为：u(x) - 一 0.29 xk 2、3(2-2)4.7 B类标准不确定度分量的自由度22v 1 U(xi)- u(xi)(2-3)2 u(x.2 u(x.根据经验,按所依据的信息来源的可信程度来判断u(xi)的相对标准不确定度 u(x)/ u(xi)0按上式计算出的自由度列于表4.表 4 u(x)/ u(xi)与 Vi 的关系 U(X)/u( xi)Vi u(xi)/u(xi)V

20、i0oo0.3060.10500.4030.20120.5020.2584.8 B类评定流程5合成标准不确定度评定5.1合成标准不确定度表示被测量Y的估计值y=f (xi,x 2,x n)的标准不确定度是由相应输入量 xi,X2,XN的标准不确定度合理合成求得的,其表示式的符号为Uc (y) o合成标准不确定度Uc (y)表征合理赋予被测量之值 Y的分散性,是一个估计标准偏差.求各个输入分量标准不确定度对输出量y的标准不确定度的奉献在求出各个输入量的不确定度分量Ui(x)之后,还需要计算传播系数(灵敏系数)最后计算由此引起的被测输出量y的标准不确定度分量：5 (y) cu(x) u(xi)(3

21、-1)xi式中传播系数或灵敏系数Ci上的含义是,输入量的估计值 xi的单位变化引x起的输出量的估计值y的变化量,即起到了不确定度的传播作用.合成标准不确定度的Uc (y)的计算公式：NN 1 NUc(y) J ()2u2(xi) 2r(xi,xj) u(xi)u(xj)(3-2),i 1 xii 1 j i 1 xi xj在实际工作中,假设各输入量之间均不相关,或有局部输入量相关,但其相关系数较小(弱相关)而近似为 r(xi,xj)=0,于是便可以化简为:N fUc(y) J ()2u2(x)(3-3)i 1 x当上1,那么可进一步化简为：xNUc(y)u2(Xi)(3-4).i 1此即计算合

22、成不确定度一般采用的方和根法,即将各个标准不确定分量平方后求其和再开跟.5.2常用的表达形式5.2.1 当简单直接测量,测量模型为y=x时,应该分析和评定测量时导致测量不确定度的各分量Ui,假设相互间不相关,那么合成标准不确定度按公式(3-5)计算：iN U：( 3-5 )1 15.2.2 当测量模型为 Y= A1X1+ A2X+ AnXn且各输入量间 互不相关 时,合成标准 不确定度可以用公式(3-6)计算：NUc(y) J A-2(x)(3-6), i 15.2.3当测量模型为Y AX1P1X；2XPN且各输入量间 互不相关 时,合成标准不确定度可使用公式(3-7)计算：JNI NPu(x

23、i)/xi2 .PUcrel(xi)2(3-7)i 1: i 1当测量模型为Y AX1X2 Xn且各输入量间互不相关时,公式(3-7)变换为公式(3-8 ):NNUc(y)/y J 爪为)/为 2 J Ucr；(xi)(3-8)-i 1 i 1注：只有在测量函数是各输入量的乘积时,可由输入量的相对合成标准不确定度Ucrel (Xi) U(X)/ Xi计算输出量的相对标准不确定度.5.2.4各输入量间正强相关,相关系数为1时,合成标准不确定度应按公式(3-9)计算：,、N f ,、Uc(y)u(x)i 1 XiCU(X) i 1(3-9)假设灵敏系数为1,那么公式(3-9)变换为(3-10):N

24、Uc(y)u(Xi)(3-10)i 15.3 关于相关性的说明5.3.1 对大局部检测工作(除涉及航天、航空、兴奋剂检测等特殊领域中要求较高的场合外),只要无明显证据证实某个分量有强相关时,均可按不相关处理,如发现分量存在强相关,如采用相同仪器测量的量之间,那么尽可能改用不同仪器分 量测量这些量使其不相关.5.3.2 如证实某些分量之间存在强相关, 那么首先判断相关性是正相关还是负相关, 并分别取相关系数为+1或-1 ,然后将这些相关分量算术相加后得到一个“净分量,再将它与其他独立无关分量用方和根求得Uc (y)05.3.3 如发现各分量中有一个占支配地位时(该分量大于其次那个分量三倍以上),

25、合成不确定度就决定于该分量.5.4 有效自由度有效自由度是指合成标准不确定度Uc(y)的自由度,用符号 Veff表示.Veff反映了 Uc(y)的可靠程度,Veff越大,Uc(y)越可靠.以下情况需要计算有效自由度 Veff ： (1)当评定某包含概率下的扩展不确定度UP时,为求得包含因子kp需要计算Uc(y)的有效自由度Veff ；(2)当客户需要了解不确定度的可靠程度而提出要求时.当各分量间相互独立且输出量接近正态分布或t分布(测量模型为线性函数)时,合成标准不确定度的有效自由度通常可按公式(3-11 )计算:u4( )一NVeffnUC4(3-11)且VeffMUi (y)i 1i 1

26、Vi当测量模型为Y AX1P1X22XPN时,有效自由度可用相对标准不确定度的形式计算,见公式(3-12):Veffuc(y)/y4N Pu(x)/x4(3-12)实际计算中,得到的有效自由度Veff不一定是一个整数,可采用将Veff数字舍位到最接近的一个较低的整数.如计算得到Veff =12.65 ,那么取Veff =12.5.5 合成标准不确定度计算流程6扩展不确定度评定6.1扩展不确定度：是被测量可能值包含区间的半宽度.扩展不确定度分为U和UP 两种.一般情况下,在给出测量结果时报告扩展不确定度Uo(1)扩展不确定度 U由合成标准不确定度 Uc乘包含因子k得到：U=kuc (4-1)当y

27、和uc( y)所表征的概率分布近似为正态分布(不确定度分量较多且其大小也比拟接近,可估计为正态分布)时,且uc(y)的有效自由度较大情况下,假设 k=2,那么由U=2uc所确定的区间具有的包含概率约为95%假设k=3,那么由U=3uc所确定的区间具有的包含概率约为 99%在通常的测量中,一般取 k=20当取其他值时,应说明其来源.当给出扩展不确定度U时,一般应注明所取的 k值；假设未注明k值,k=2o(2)当要求扩展不确定度所确定的区间具有接近于规定的包含概率p时,扩展不确定度用符号UP表示,当p为0.95,0.99时,分别表示为 5和U90U = kpUc (4-2)kp是包含概率为p时的包

28、含因子.kp= t p( Veff)(4-3 )根据合成标准不确定度 Uc(y)的有效自由度Veff和需要的包含概率,查?t分 布在不同概率p与自由度v时的tp(v)值(t值)表?得到tp(Veff)值,该值即包含 概率为p时的包含因子kp值.如果合成不确定度中 A类分量占比重较大,如把23而且作A类评估时重复 u测量次数n较少,那么包含因子k必须查t分布表获得.扩展不确定度UP = k pUc(y)提供了一个具有包含概率为p的区间y±U.在给出UP时,应同时给出有效自由度Veff 0(3)如果可以确定 Y可能值的分布不是正态分布,而是接近于其他某种分布,那么 不应按kp= t p(Veff )计算 UP.例如Y可能近似为矩形分布,取 p=0.95时kp=1.65弋3 0.95;取p=0.99时 kp=1.71 Q氢 0.99 ;取 p=1 时 kp=1.73 弋任 o图4正态分布概率分布图6.2扩展不确定度的有效位数估计值y的数值和它的合成标准不确定度Uc(y)或扩展不确定度U的数值均不应给出过多的有效位数.通常最终报告的Uc(y)和U最多为两位有效数字.对各标准不