不确定度测定汇总

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

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

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

4、测量不确定度 简称不确定度，是指 根据用到的信息，表征赋予被测量值分散性的非负 参数。”测量不确定度一般由若干分量组成。其中一些分量可根据一系列测量值的统计分布，按测量不确定度的A类评定（随机效应引起的）进行评定，并用标准偏差表征；而另一些分量 则可根据基于经验或其它信息所获得的概率密度函数，按测量不确定度的B类评定（系统效应引起的）进行评定，也用标准偏差表征。标准不确定度（全称为标准测量不确定度）可采用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(Xi)相对标准不确定度Urelrel 表示 相对”的英文字母的缩写。测量不确定度的A类评定UA 或 UA (Xi)测量不确定度的B类评定UB 或 UB (Xi)合成标准不确定度Uc或 Uc (y)相对合成标准不确定度Llcrel 或 Ucrel (y)扩展不确定度U或UpUp包含概率为p的扩展不确定度相对扩展不确定度Urel或 Uprel包含因

7、子k或 kpkp包含概率为p的包含因子包含概率p如，p :=95% , p = 99%。有效自由度Veffeff -表示 有效”的英文字母的缩写。注： 表中 A、B、c、rel、eff为正体；x、y、k、i、p、n、u、U 为斜体。 表中大写U表示扩展不确定度；小写 u表示标准不确定度，如：标准不确定度A类评定：UA标准不确定度B类评定：UB合成标准不确定度，uc或Uc (y) 扩展或相对扩展不确定度，U或Up、Sei或Uprei2.11测量模型是指测量中涉及的所有已知量间的数学关系。测量模型简称模型。测量模型的通用形式是方程：f (Y, Xi,，Xn) = 0,其中测量模型中的 输出量Y是被

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

9、下对同一被测量(一个被测件)独立重复观测n次，得到n个观测值Xi (i=1,2，,n),被测量X的最佳估计值是n个独立测得值的算术平均值X，按公式(1-1)计算:1 nXXi(1-1)n i 4单个测得值Xk的实验方差s2(xk),按公式(1-2)计算：2 1 n 2s (Xk)(Xi -X)(1-2)n 1 y单个测得值Xk的实验标准偏差s(Xk)，按公式(1-3)计算:nn 二s(Xk) J/ (Xi-x)(1-3)Y n-1 y公式(1-3)是贝塞尔公式，自由度 v为n-1。实验标准偏差s(xq表征了测得值x的分散性, 测得重复性用s(Xk)表征。被测量估计值X的A类标准不确定度为：uX

10、)二 s(X)二 s(Xk) . , n( 1-4)A类标准不确定度u丄(X)的自由度为实验标准偏差 s(Xk)的自由度，即v=n-1。实验标准偏差s(X)表征了被测量估计值 X的分散性。3.2在规范化的常规检定、校准或检测中评定合并样本标准偏差若对每个被测件的被测量Xi在相同条件下进行n次独立测量，测得值为Xi1,X2,Xin,其平均值为Xi ;若有m个被测件，则有m组这样的测得值，可按公式(1-5)计算单个测得值的合成样本标准偏差Sp(Xk):sXk)=m(nm n玮汪jX)2式中： i组数，i=1,2，,m;j 每组测量的次数,(1-5)j=1,2，n。公式(1-5)给出的sp(Xk),

11、其自由度为m(n-1)。若对每个被测件已分别按n次重复测量算出了其实验标准偏差环则m组测得值的合并样本标准偏差Sp(xk)可按公式(1-6)计算:sm(Xk)=(1-6)当实验标准偏差s的自由度均为v0时，公式(1-6)给出的Sp(xk)的自由度为mv0。若对m个被测量Xi分别重复测量的次数不完全相同，设各为ni，而Xi的实验标准偏差s(x)的自由度为Vi，通过m个Si与w可得Sp(xk)按公式(1-7)计算：(1-7)SXk)二m公式(1-7)给出Sp(Xk)的自由度为7 =、Vi oi J由上述方法对某个被测件进行n次测量时，所得测量结果最佳估计值的A类标准不确定度为：U X)二 s(X)

12、二 sXk). . n(1-8)用这种方法可以增大评定的标准不确定度的自由度，也就提高了可信程度。3.3预评估重复性在日常开展同一类被测件的常规检定、校准或检测工作中，如果测量系统稳定，测得重复性无明显变化，则可用该测量系统以与测量被测件相同的测量程序、操作者、操作条件和地点，预先对典型的被测件的典型被测量值进行n次测量(一般n不小于10),由贝塞尔公式计算出单个测得值的实际标准偏差S(Xk),即测量重复性。在对某个被测件实际测量时可以只测量门次(K nl当1 （只测一次）A类标;隹不确定度u（x） = 5（）自由度为4V 丿4标准不确定度的B类评定（分量）4.1 B类评定的一般表示B类评定的

13、方法是根据有关的信息或经验，判断被测量的可能值区间X -a, X +a,假设被测量值的概率分布， 根据概率分布和要求的包含概率p估计因子k,则B类标准不确定度UB可由（2-1）式得到：a u 扩 k（2-1）式中：a为被测量可能值区间的半宽度。当k为扩展不确定度的倍乘因子时称包含因子，其他情况下根据概率论获得的k称置信因子。4.2 B类评定（来源）通常基于诸如以下信息：（1 ）权威机构发布的量值；（2）有证标准物质的量值；（3）校准证书；（4）仪器的漂移；（5）经检定的测量仪器的准确度等级；（6）根据人员经验推断的极限值等。4.3确定B类评定的区间半宽度 a（1）生产厂提供的测量仪器的最大允许

14、误差为址，或由手册查出所用的参考数据误差限为，或当测量仪器或实物量具给出准确度等级等，并经计量部门检定合格，则评定仪器的不确定度时，可能值区间的半宽度为：a =(2) 校准证书提供的校准值，给出了其扩展不确定度为U，则区间的半宽度为：a =U(3) 由有关资料查得某参数的最小可能值为a-和最大值为a+,最佳估计值为该区间的中点,则区间半宽度可以用下式估计：a =(a+- aj/2(4) 必要时，可根据经验推断某量值不会超出的范围，或用实验方法来估计可能的区间。4.4 k的确定方法(1 )已知扩展不确定度是合成标准不确定度的若干倍时，该倍数就是包含因子k值。(2)假设被测量值服从正态分布时，根据

15、要求的概率查表2得到k值。表2正态分布情况下概率 p与k值间的关系p0.500.680.900.950.95450.990.9973k0.67511.6451.96022.5763(3)假设为非正态分布时，根据要求的概率查表3得到k值。表3常用非正态分布时的k值及B类标准不确定度Ub(x)分布类别p (%)kUb(X)三角100a/v6梯形(3 =0.71)1002a/2矩形(均匀)100a/反正弦100a/力两点1001a注：表3中B为梯形的上底与下底之比，对于梯形分布来说，k= 6/(1 1 2)，特别当B等于1时，梯形分布变为矩形分布；当 B等于0时，变为三角分布。4.5 B类评定概率分

16、布的假设(1) 被测量受许多随机影响量的影响，当它们各自的影响都很小时，不论各影响量的概率 分布是什么形式，被测量的随机变化服从正态分布。如证书或报告给出的不确定度是具有包含概率为0.90、0.95、0.99的扩展不确定度(即给出U90、U95、U99),此时，除非另有说明，可按正态分布来评定B类标准不确定度。(2 )当利用有关信息或经验，估计出被测量可能值区间的上限和下限，其值在区间外的可能几乎为零时，若被测量值落在该区间内的任意值处的可能性相同，则可假设为均匀分布(或称矩形分布、等概率分布)。如数据修约、测量仪器最大允许误差或分辨力、参考数据的误差限、度盘或齿轮的回差、平衡指示器调零不准、

17、测量仪器的滞后或摩擦效应导致的不确定 度及对被测量的可能值落在区间内的情况缺乏了解等，一般假设为均匀分布。(3 )当利用有关信息或经验，若被测量值落在该区间中心的可能性最大，则假设为三角分布。如两相同均匀分布的合成、两个独立量之和值或差值服从三角分布。(4 )当利用有关信息或经验，若落在该区间中心的可能性最小，而落在该区间上限和下限的可能性最大，则可假设为反正弦分布(即U形分布)。如度盘偏心引起的测角不确定度、正弦振动引起的位移不确定度、无线电测量中失配引起的不确定度、随时间正弦或余弦变化的温度不确定度等。(5 )按级使用量块时，中心长度偏差的概率分布可假设为两点分布。(6) 安装或调整测量仪

18、器的水平或垂直状态导致的不确定度常假设为投影分布。(7) 实际工作中，可依据同行共识确定概率分布。4.6分辨力导致的B类不确定度分量若数字显示器的分辨力为由分辨力导致的标准不确定度分量u(x)采用B类评定，则3得k =-.3 ，因此由分辨力导致的标准不确定度分量u(x)为：u(x)二k0.29 x2*3(2-2)区间半宽度为a= & 12，假设可能值在区间内为均匀分布，查表4.7 B类标准不确定度分量的自由度(2-3)1 u2(Xi)1 ：u(x)Vi22 二 U(xjJ 2 I u(xj根据经验，按所依据的信息来源的可信程度来判断u(Xi)的相对标准不确定度 u(Xi)/u(Xi)。按上式计

19、算出的自由度列于表4。表 4 u(Xi)/ u(Xi)与 Vi 的关系 u(Xi)/ u(Xi)Vi u(Xi)/ u(Xi)Vi0oo0.3060.10500.4030.20120.5020.2584.8 B类评定流程5合成标准不确定度评定5.1合成标准不确定度表示被测量Y的估计值y=f ( Xi,X2,XN)的标准不确定度是由相应输入量 X1,X2,XN的标准不 确定度合理合成求得的， 其表示式的符号为Uc(y)。合成标准不确定度Uc(y)表征合理赋予被 测量之值Y的分散性，是一个估计标准偏差。求各个输入分量标准不确定度对输出量y的标准不确定度的贡献在求出各个输入量的不确定度分量Ui(X)

20、之后，还需要计算传播系数(灵敏系数)Ci,最后计算由此引起的被测输出量y的标准不确定度分量：cfU(y) = c u(xj =-u(Xi)(3-1)cx式中传播系数或灵敏系数ci-的含义是，输入量的估计值Xi的单位变化引起的输出cXi量的估计值y的变化量，即起到了不确定度的传播作用。合成标准不确定度的 (y)的计算公式：Cf 2 2 亠NT冴话Uc(y).(一 ) u (Xi) 2-r(Xi,Xj) u(Xi)u(Xj)(3-2)V iCXii# j=fc十 CXi CXj在实际工作中，若各输入量之间均不相关，或有部分输入量相关，但其相关系数较小(弱相关)而近似为r(Xi,Xj)=O,于是便可

21、以化简为:uc(y)石(如x);i 4:-X(3-3)当f =1，则可进一步化简为:;:Xi%(y)= ： u2(Xi)(3-4)此即计算合成不确定度一般采用的方和根法，即将各个标准不确定分量平方后求其和再 开跟。5.2常用的表达形式5.2.1当简单直接测量，测量模型为y=x时，应该分析和评定测量时导致测量不确定度的各分量Ui，若相互间不相关，则合成标准不确定度按公式(3-5)计算：NUc(y)u2(3-5)5.2.2当测量模型为 Y= A1X1+ A2X2+ AnXn且各输入量间 互不相关时，合成标准不确定度 可以用公式(3-6)计算：:N(3-6)瓦 A2u2(x)i =15.2.3当测量

22、模型为Y二AXXJXN且各输入量间 互不相关 时，合成标准不确定度可 使用公式(3-7)计算：fN%(y)/y = 送 Pu(Xi)/X 2 =严 Pucrel(N)F(3-7)当测量模型为 Y =AX1X2 XN且各输入量间 互不相关 时，公式(3-7)变换为公式(3-8 )：rNrNuc(y)/ yMxJ/Xi2ucrel2(Xi)(3-8)i i =1i =1注：只有在测量函数是各输入量的乘积时，可由输入量的相对合成标准不确定度u“el(Xi )二u(Xj )/ Xj计算输出量的相对标准不确定度。5.2.4各输入量间正强相关，相关系数为1时，合成标准不确定度应按公式(3-9)计算:NNU

23、c(y) = Zk4-xiu(x) = E Gu(Xi)7(3-9)若灵敏系数为 1则公式（3-9）变换为（3-10）:N(3-10)%(y)八 u(Xi)i 二5.3关于相关性的说明5.3.1对大部分检测工作（除涉及航天、航空、兴奋剂检测等特殊领域中要求较高的场合外），只要无明显证据证明某个分量有强相关时，均可按不相关处理，如发现分量存在强相关， 如采用相同仪器测量的量之间，则尽可能改用不同仪器分量测量这些量使其不相关。5.3.2如证实某些分量之间存在强相关，则首先判断相关性是正相关还是负相关，并分别取相关系数为+1或-1，然后将这些相关分量算术相加后得到一个“净”分量，再将它与其他 独立无

24、关分量用方和根求得uc （y）。5.3.3如发现各分量中有一个占支配地位时（该分量大于其次那个分量三倍以上），合成不 确定度就决定于该分量。5.4有效自由度有效自由度是指合成标准不确定度Uc（y）的自由度，用符号 veff表示。veff反映了 Ujy）的可靠程度，Veff越大，Uc（y）越可靠。以下情况需要计算有效自由度Veff：（1 ）当评定某包含概率下的扩展不确定度Up时，为求得包含因子kp需要计算Uc（y）的有效自由度veff；（2）当客户需要了解不确定度的可靠程度而提出要求时。当各分量间相互独立且输出量接近正态分布或t分布（测量模型为线性函数）时，合成标准不确定度的有效自由度通常可按公

25、式（3-11）计算:Ui4(y)(3-11 )且Nve f fvii T当测量模型为 Y = AX1PXX,时，有效自由度可用相对标准不确定度的形式计算，见公式（3-12）：(3-12)Uc(y)/y4N Ui(N)/Xi F, vi实际计算中，得到的有效自由度Veff不一定是一个整数，可采用将Veff数字舍位到最接近的一个较低的整数。如计算得到Veff =12.65，则取Veff =12。5.5合成标准不确定度计算流程6扩展不确定度评定6.1扩展不确定度：是被测量可能值包含区间的半宽度。扩展不确定度分为U和Up两种。一般情况下，在给出测量结果时报告扩展不确定度U。（1） 扩展不确定度 U由合

26、成标准不确定度 Uc乘包含因子k得到：U = kUc（ 4-1）当y和Uc（y）所表征的概率分布近似为正态分布（不确定度分量较多且其大小也比较接近，可估计为正态分布）时，且uy）的有效自由度较大情况下，若k=2，则由U=2uc所确定的区间具有的包含概率约为95%。若k=3，则由U=3uc所确定的区间具有的包含概率约为99%。在通常的测量中，一般取k=2。当取其他值时，应说明其来源。当给出扩展不确定度U时，一般应注明所取的k值；若未注明k值，k=2。（2） 当要求扩展不确定度所确定的区间具有接近于规定的包含概率p时，扩展不确定度用符号Up表示，当p为0.95, 0.99时，分别表示为 U95和U

27、99。Up = kpuc（4-2）根据合成标准不确定度 Uc(y)的有效自由度Veff 和需要的包含概率，查t分布在不同概 率p与自由度v时的tp(v)值(t值)表得到tp(Veff)值，该值即包含概率为 p时的包含因子 kp值。如果合成不确定度中A类分量占比重较大，如Uc(y) : 3而且作A类评估时重复测量次UA数n较少，则包含因子 k必须查t分布表获得。扩展不确定度 Up = kpUc(y)提供了一个具有包含概率为p的区间yJp。在给出Up时，应同时给出有效自由度Veff。(3) 如果可以确定 Y可能值的分布不是正态分布，而是接近于其他某种分布，则不应按kp= tp(Veff)计算 Up。例如Y可能近似为矩形分布，取p=0.95时kp=1.653 0.95 ;取p=0.99时kp=1.713 0.99 ；取 p=1 时 kp=1.73&. 3。图4正态分布概率分布图6.2扩展不确定度的有效位数估计值y的数值和它的合成标准不确定度Uc(y)或扩展不确定度U的数值均不应给出过多的有效位数。通常最终报告的Uc(y)和U最多为两位有效数字。对各标准不确定度分量u(Xi),为了在连续计算中避免修约误差导致不确定度，可以适当保留多余的位数。在报告最终结果时，一般采用GB/T 8170-2008数值修约规则与极限数值的表示和判定修约到需要的有效数字。如U=28.05kHz经修约写成28k