不确定度的表示方法

1、不确定度数据表示方法一.不确定度概述：在科学实验、产品生产、商业贸易及日常生活的各个领域，我们都要 进行测量工作。测量的目的是确定被测量的值，测量不确定度表示测量结果的不确 定或不肯定的程度，也就是不可信度。定义：不确定度是与测量结果相关联的，用于合理表征被测量值分散性大小 的参数。分类及表示：标准不确定度：以标准差表示的不确定度，以M表示。扩展不确定度：以标准不确定度的倍数表示的不确定度，以U表示。（扩展不确定 度表明了具有较大置信概率的区间的半宽）合成标准不确定度：各标准不确定度分量的合成，以M c表示（测量结果标准差的 估计值）1.1.合成标准不确定度被测量y由N个其他量xi的函数确定时

2、，假设其函数关系为y=f（x1，x2，xN）/ co=jx * a）+2sx普畀“吃内）心）酿匕）V i=多曰尸i+1气0XjdXi上式称为不确定度传播率。为灵敏系数，r （xi, xj）上式称为不确定度传播率。为灵敏系数，r （xi, xj）为相关系数。当被测量的函数形式为：y=A1x1+A2x2+ANxN，且各输入量之间 不相关时，合成标准不确定度为：“（y）= 片=四片。）若用灵敏系数表示:当被测量的函数形式为:-（X）/ 的不合成标准不确定度为:合成标准不确定度为:所有输入量都相关，且相关系数为1时，合成标准不确定度为:uc(y)：合成标准不确定度ui(x)：各输入量的标准不确定度：

3、ui(x)的自由度eff越大表明评定的合成标准不确定度uc(y)越可靠。自由度的含义：自由度是方差之不确定度的度量，由于测量不确定度用标准 偏差(方差的正平方根)表示，自由度也就是“测量不确定度的不确定度”。自由度 大表示测量不确定度的不确定度小，即测量结果之不确定度的可信度高，反之亦 然。用第一张ppt的例子来说明，当自由度很大时，表示“被测量的值落在831。 9 839.1 区间的置信水平约为95 %”的可信度高，对于自由度v= 12, 3.6 的不可信度大约是21 %。1.2扩展不确定度分为两种U和Up。U：就是合成标准不确定度的倍数，U=kuc,即由合成标准不确定度直 接乘以包含因子k

4、 ( k的典型值为23)Up：对于给定的置信概率P,扩展不确定度记为Up=kpuc,此时包含 因子kp的选择如下如果组成uc的不确定度分量较多，且各分量对不确定度的影响不大时，据中 心极限定理，合成不确定度uc的分布接近正态分布。若有效自由度充分大，按正态分布计算若有效自由度较小，按t分布计算(按有效自由度查表)如果uc的概率分布为非正态分布时，应根据相应的分布确定kp。.不确定度的评定测量不确定度的评定方法分为两类，即A类和B类，两者之间无主次之分， 享有同等地位。A类不确定度评定A类不确定度是采用观察列进行统计分析的方法来评定标准不确定度的，用 标准误差来表示。测量列算术平均值的标准误差为

5、当测量次数较少时，其估算值会偏大，这是，从理论上可得A类不确定度的 估算值为UA=t (nT)式中，t (n-1)是一个大于1的修正值(被称为t分布临界值)。测量次数 n不同，修正量t (n-1)不同。下表给出了不同测量次数n对应的修正量。n2345610208t (n-1)1.841.321.201.141.111.051.031.00B类不确定度的评定实则基于对一个事件发生的信任程度。很多不确定度分量实际上还必须用别 的非统计方法来评定。B类不确定度评定的信息来源主要有六项：:以前的测量数据。 :对有关数据资料和测量仪表特性的了解和经验。:生产部门提供的技术说明文件。:校准证书、检定证书或

6、其他文件提供的数据、准确度的等别或级别，包括目前还 在使用的极限误差等。 :手册或某些资料给出的参考数据及其不确定度。:规定实验方法的国家标准或类似技术文件中给出的重复性限r或重复性限R。B类不确定度的评定方法采用不同于A类的其他方法估算。首先，根据仪器、仪表说明书，国家标准，材料特性等来确定测量误差限A ，例如，已知仪表精度等级和量程可计算出误差限。其次，确定测量误差的分布，常见的有正态分布和均匀分布。最后，将测量误差限（对应的置信度1）换算成相似的标准误差 （对应一倍的标准误差置信度）。对于均匀分布的误差，其B类不确定度估算为叫金对于服从正态分布的误差，其B类不确定度估算为.测量不确定度评

7、定实例用电压表测量稳压电源的输入电压测量方法及测量的数学模型用已经校准的电压表测量一台稳压电源的输出电压U。电压表的分辨力为 0.01V。电压表校准的不确定度和表的分辨力引起的不确定度可以忽略不计。因 此，多次直接测量，数据的平均值即为输出电压的最佳估计值。故测量的数学模型 可以表示为：U=U 测（1.1）测量数据进行了 10次测量，测量数据及相关计算列于表1.1表1.1输出电压测量数据及相关计算序号i输出电压Ui/V残差u i/V残差平方u 2i/V21200.63+ 0.070.00492200.28-0.280.07843200.08-0.480.23044201.36+ 0.800.6

8、4005200.47-0.090.00816200.82+ 0.260.06767199.73-0.830.68898201.09+ 0.530.28099200.41-0.150.022510200.73+ 0.170.0289测量值平均：U = 200.56VN +q=183, N-q=-183=2.0506检查平均值和残差的计算是否有误，可将正残差与负残差分别相加，若两个和的绝对值不相等，且两者之差大于末位的1/2,则可判定计算有误。本例中|S+| =183，再复核计算，表明计算正确。也可直接求残差的代数和看是否为零，或小于末位的半个单位来进行判断。10次测量值的平均即为输出电压U的最佳

9、估计值。3.1.3根据贝塞尔公式计算测量列的实验标准差单次测量值的实验标准差S(Ui)表征测量列中测量数据的分散性。假定测量值服从正态分布，就可以 估计，大约有68.3%的测量值处在(200.560.48) V区间内，95%的测量值处在 (200.56 + 2X0.48) V 区间内，99.7%的测量值处在(200.563X0.48) V 区间 内。残差绝对值大于3X0.48V的测量值不应该出现(小概率事件)。如果出现， 可判定为粗大误差。10次测量的每一个测量值的实验标准差均为0.48V。这10个测量值仅是测量值总体的一个样本。由此计算的标准差仅是这个样本 的标准差，而不是总体标准差。总体标

10、准差可表示为：这无法实际测得，只是理论上存在，又叫理论标准差。而样本标准差仅是理 论标准差的有偏估计值。样本方差S2 (Ui)才是总体理论方差。2 (Ui)的无偏估 计值，即最佳估计值。所以在统计分析中，多用方差作为数据分散性的度量。标准 差是方差的正平方根，在实际工作中使用更为方便。3.1.4计算平均值的实验标准差s(q)=s(uj/薪口 015(v)可以这样理解S (U测)的含义：再进行若干组测量，每组n个测量值均可求 得一个平均值和相应的测量列的标准差。若测量条件不变，各组测量列的标准差相 互很接近。各组测量值的平均值不会完全一样，它们也构成一个数列。由这个数列 也可用贝塞尔公式计算它的

11、标准差。由统计理论可证明，平均值数列的标准差等于 单个测量列标准差的1/n倍，即为(1.5)式。这表明，平均值数列的分散性比单 一测量列的分散性小，即各平均值比测量列中的测量值相互之间更靠近。v = n a叫做自由度，它是求和的项数n,减去对和的限制数a。上例中，n = 10,限制条件仅有残差之和为零，即Zu i = 0,故a=1,因而v=10 1 = 9。 自由度越大，计算的样本标准差就越接近总体标准差，所得结果的可信度就越高。3.1.5测量结果报告平均值的实验标准差就是它的标准不确定度。本例中，这是唯一的不确定度 来源。认为近似服从正态分布，测量结果可以表示为：输出电压U = 200.56

12、V, u = S (U 测)=0.15V,(1.6)或者 U = 200.56V2Xu= 200.56V0.30V, p = 95%(1.7)0.30V即为扩展不确定度。说明：测量值的分散是输出电压的随机变化和测量仪表读数随机起伏的综合 反映。3.2用发光强度标准灯校准光照度计的示值校准方法如图(2.1)所示，将光强标准灯和照度计的光度探测器安置在测光导轨上， 并调整好它们的状态，设定两者之间的距离为1。则标准灯在接收面上产生的标准 照度值为ES = I / 12(2.1)式中，I为标准灯的发光强度，I = 268.8 cd, 1 = 1.600 m,故有：ES = 268.8 cd/ (1.

13、600 m) 2 = 105.0 lx(2.2)由照度计测得的照度示值为Et图2.1发光强度标准灯校准光照度计示意图数学模型如上所述，校准的数学模型可以表示为= & _ 耳=& (2.3)式中，为照度计的示值误差。输入量的标准不确定度评定(1)由4E分别对Et, I和l求偏导数，即得相应的灵敏系数1弭C /(回=213 di f由此可得：u =耳)的引他(?)=(，)Et为照度计10次测量所得示值Eti的平均值，即Et = Z Eti/10。故 用A类方法评定其标准不确定度。10次测量值及相关计算列于表示2.1表2.1测量数据一览表序号测量值Eti / lx残差u i / lxi?.2 th2

14、198.93 0.580.3364299.83+ 0.320.1024399.36 0.150.02254100.10+ 0.590.3481599.55+ 0.040.0016699.20-0.310.0961799.50-0.010.0001899.48-0.030.0009999.64+ 0.130.01691099.54+ 0.030.0009平均Et99.51= 0.9259根据贝塞尔公式，计算得实验标准差s (Eti)为则平均值的实验标准差s (Et)为“是叫味氏二。0 bt因此，A类方法评定的标准不确定度为(2.8)u1 = s (Et)=0.10 lx(2.8)其自由度为v=n

15、-1=9u1反映了各种随机因素，如标准灯供电的随机起伏引起发光的随机起伏，灯 丝本身的分子热运动产生的发光起伏，空气绕动使接收面上的照度发生随机起伏， 可能出现的微小震动使接收面和灯丝面在平衡位置附近摆也使接收面上的的照度产 生随机起伏，以及若干没有认识到的随机因素影响的综合作用，使得接收面上的照 度不是恒定的，因而照度测量值不重复。此外，照度计本身在各种随机因素影响 下，其响应度也不是恒定不变的，也是测量值不重复性的原因。上面用统计方法求 得的标准差或标准不确定度正是测量值不重复性的一个评定，它表征测量值的分散性。（3）I的不确定度有两部分，均用B类方法评定其标准不确定度a）检定证书上给出的

16、光强值的扩展不确定度为1.0%,并注明包含因子k = 3 （近似服从正态分布）。由此可得光强值的相对标准不确定度为与=等1 = 0.33%(止 26&K cdx0.0033 J 0,89 cd所求得的u （I1）本身的相对标准不确定度估计为25%,它大致相应于正态分 布9次观测值的平均值的标准差的相对标准差。在确定标准灯光强值时，观测次数 一般不少于9次，在没有更准确的资料时，就取故相应的自由度为I1的不确定度对的不确定度贡献为以(4)_ 0.89 cd以(4)_ 0.89 cd2(1.600姆-0.35 lxb）供给标准灯的电流在数值上与检定时供给的电流完全一样。由于电测系 统的不确定度，在

17、这两种情况下，供给灯的实际电流，一般说来不会完全相同。估 计两者的最大差异可达0.03%，即电测系统的扩展不确定度为0.015%，置信水准为 1。两者的差别从0到0.03%都以等概率出现，认为服从均匀分布，则电流值的相 对标准不确定度为Wr(i)= 0.0087%di7di7一 7故，由于灯电流的不确定度引起光强的相对标准不确定度为uf （/2） = 6x0.0087% 0.052%则w（Z2） = 268.8 cdx0.00052n 0.14 cd因而对E的不确定度的贡献为0.14 cd0140.14 cd014 I 0.055 lx2.56病由于电流的置信限是准确知道的，可以认为求得的un

18、2是准确的，故V（/?） oo如果标准灯工作和检定时使用同一电测系统，则这一项不确定度分量就不予 考虑。（4）测量距离l的标准不确定度也用B类方法评定根据对导轨长度标尺的校准和接收面，灯丝平面的调整资料，可以判定l值可能变动的范围不会超过1mm，而且出现在区间中心的概率比出现在两端的概率要大得多，可以看作服从三角分布。因而测量l的标准不确定度为对照度计示值误差标准不确定度的贡献为：4 =陵? x 0.41 x 10-$ 附=0,054 lx（1.600同样，可以认为u（1）是准确的，故（5）合成标准不确定度综上所述，将各不确定度分量和相关信息列于表2.2表2.2标准不确定度分量一览表不确定度来源标准不确定度值灵敏系数 d