测量不确定度的评定与表示2015.5.28课件.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,测量不确定度评定与表示 JJF 1059.1-2012,雷雷雷阿雷,2015.5.28,一、,不确定度的基本概念,二、,评定测量不确定度的一般流程,分析不确定度来源和建立测量模型,标准不确定度的A类评定（预评估重复性）,标准不确定度的B类评定,合成标准不确定度（灵敏系数）,扩展不确定度,测量不确定度的表示,三、,声明检测或校准结果与规范符合性,一、不确定度的基本概念,为何与报告的结果不一致呢？,测量不确定度是对任何测量的结果存有怀疑。,世间本无“绝对”，每一次测量过程，操作人员，使用设备，测量方法，所处环境等因素都会对测量结果有影响。,每一次测量过程，我们都得到了一个测得值，但上述等等因素对我们得到的这个结果有没有影响，有多大影响，都值得我们去怀疑，也就是不确定度。,在日常说话中，这可以表述为“出入”。例如一根绳子可能2米长，有1cm的“出入”,“出入”也就是所谓的不确定。,一、不确定度的基本概念,一、不确定度的基本概念,不确定度值：40V0.2%=0.08V测得值：39.99V可能分布的范围：39.99V0.08V即：39.91V40.07V,一、不确定度的基本概念,本无“绝对”，即便是这里的可能分布的范围（不确定度大小）也是有概率的。,k,=2，表示有95%的信心（包含因子，包含概率）。,某种意义上也就是说，在该数字多用表40V这个测量点，后续进行重复测量，有95%的概率都是落在39.91V40.07V这个范围内。,一、不确定度的基本概念,测量不确定度定义,测量不确定度（uncertainty of measurement）：,根据所用到的信息，表征赋予被测量值分散性的非负参数。,注：1.,此参数可以是诸如标准测量不确定度的标准偏差（或其特定倍数），或说明了包含概率的区间半宽度。,2.测量不确定度由若干分量组成。其中的一些分量可根据一系列测量结果的统计分布，按测量不确定度的A类评定进行评定，并可用标准偏差表征。另一些分量可用基于经验或其他信息获得的概率密度函数，按测量不确定度的B类评定进行评定,也用标准偏差表征。,一、不确定度的基本概念,经验？信口开河？,测量不确定度由若干个分量组成,注2,“测量不确定度由若干分量组成。其中的一些分量可根据一系列测量结果的,统计分布,，按测量不确定度的A类评定进行评定，并可用,标准偏差,表征。另一些分量可用基于经验或其他信息获得的,概率密度函数,，按测量不确定度的B类评定进行评定,也用,标准偏差,表征。”,一、不确定度的基本概念,标准不确定度（standard uncertainty）：,以,标准偏差,表示的测量不确定度。,实验标准偏差（experimental standard deviation）：,对同一被测量进行n次测量，表征测量结果分散性的量。用符号,s,表示。（,）,一、不确定度的基本概念,测量不确定度的A类评定（type A evaluation of measurement uncertainty）：,对在规定测量条件下测得的量值用,统计分析,的方法进行的测量不确定度分量的评定。,测量不确定度的B类评定（type B evaluation of measurement uncertainty）：,用不同于测量不确定度A类评定的方法对测量不确定度分量进行的评定。（,基于经验或其他信息获得的,概率密度函数,）,一、不确定度的基本概念,合成标准不确定度（combined standard uncertainty）：,由在一个测量模型中各输入量的标准测量不确定度获得的输出量的标准测量不确定度。,（简单的说就是每一个分量求出来之后合成一个总的不确定度，但并非简单的加法）,相对标准不确定度（relative standard measurement uncertainty）：,标准不确定度除以测得值的绝对值。,一、不确定度的基本概念,一、不确定度的基本概念,包含因子（coverage factor）：（,k,）,为获得扩展不确定度，对合成标准不确定度所乘的大于1的数。,包含区间（coverage interval）：,基于可获得的信息确定的包含被测量一组值的区间，被测量以一定概率落在该区间内。,包含概率（coverage factor）：（,p,）,在规定的包含区间内包含被测量的一组值的概率。,预备知识,分布数据散布的“形状”,一组数值的散布会取不同的形式，或称为服从不同的概率分布。,正态分布,在一组读数中，较多的读数值靠近平均值，少数读数值离平均值较远。这就是正态分布或高斯分布的特征。,正态分布是最常见的一种分布，简单的理解就是中间群体是占大多数的，极高和极低都是少部分。,正态分布,P(%),68.27,90,95,95.45,99,99.73,K,1,1.65,1.96,2,2.58,3,正态分布,右图的,即标准偏差，也可以理解为（合成）不确定度，可以看出，在2（,k,）倍标准差时，概率由68.26%提升至95.45%，但此时是中间值（测得值）加减2倍标准差（扩展不确定度）。,一、不确定度的基本概念,扩展不确定度（expanded uncertainty）：,合成标准不确定度与一个大于1的数字因子的乘积。,一般来说，包含因子越大，所求得的扩展不确定度值会变大，包含概率（信心）增加。,包含概率增加是我们希望看到的，意味着发生在包含区间之外的概率极低，也就是说不在我们可控范围（信心范围）内的情况概率极低。,但扩展不确定值变大是我们不希望看到的，意味着包含区间变大。即便有极高的包含概率，如果扩展不确定度过大，那也是没有意义的。,实际中，一般去寻求概率和扩展不确定度值的平衡，比如取包含因子,k,=2，也就是包含概率大概为,p,=95%，这样求得的扩展不确定度为合成标准不确定度的2倍，概率由约68%提升至约95%，两者都达到一个可以接受的范围。,二、评定测量不确定度的一般流程,分析不确定度来源和建立测量模型,标准不确定度,计算合成不确定度,确定扩展不确定度,报告测量结果,分析不确定度来源,在实际测量中，有许多可能导致测量不确定度的来源。例如：,a)被测量的定义不完整,b)复现被测量的测量,方法,不理想,c),取样,的代表性不够，即被测样本不能代表定义的被测量,d)对测量过程受,环境影响,的认识不恰如其分或对环境的测 量与控制不完善,e),对模拟式仪器的,读数,存在人为偏差,分析不确定度来源,在实际测量中，有许多可能导致测量不确定度的来源。例如：,f)测量仪器的,计量性能,（如,灵敏度,，鉴别力阈，,分辨力,，死区及,稳定性,等）的局限性,g),测量标准或标准物质的不确定度（准确度）,h)引用的数据或其他,参量的不确定度,i)测量方法和测量程序的近似和假设,j)在相同条件下被测量在,重复观测,中的变化,测量模型的建立,测量中，当被测量（即输出量）Y由N个其他量X,1,X,2,，X,N,(即输入量），通过函数f来确定时，则下式称为测量模型：,Y,=,f,（,X,1,X,2,，X,N,）,例：一个随温度,t,变化的电阻器两端的电压为V，在温度为,t,0,（20）时的电阻为R,0,，电阻器的温度系数为a，则电阻器的损耗功率P（被测量）取决于V，R,O,，a和t，即测量模型为：,P=,f,（V，R,0,，a，t）=,测量模型的建立,在简单的直接测量中测量模型可能简单到下式：,Y=X,1,-X,2,（如：,=X,1,-X,2,）,甚至简单到下式：,Y=X,例如：用压力表测量压力，被测量（压力）的估计值y就是仪器（压力表）的示值x。测量模型为y=x。,标准不确定度的A类评定,贝塞尔公式法,在重复性条件或复现性条件下对同一被测量独立重复观测n次，得到n个测得值 （i=1,2，n),被测量X的最佳估计值是n个独立测得值的算术平均值 ，计算公式为：,标准不确定度的A类评定,贝塞尔公式法,单个测得值 的实验标准偏差 ，按下式计算：（贝塞尔公式）,此式是单次测量的实验标准偏差(,)，也就是标准不确定度,u,(,x,)。自由度（反应了相应实验标准偏差的可靠程度）,v,=,n,-1。,标准不确定度的A类评定,贝塞尔公式法,m次测量平均值 的A类标准不确定度 ，按下式计算：,如果测量是以m次测量的平均值为报告结果时，则,标准不确定度的A类评定,例：,在我们平时校准耐电压测试仪的电压参数时，按检定规程要求，取2次测得值的平均值为报告结果，则此时m=2.,即我们在评定耐电压测试仪电压参数由重复性引入的标准不确定时，可先进行n=10次重复测量，求出单次测量的实验标准差,s,(x,k,),因报告时是取m=2次测量值的平均值，则由重复性引入的标准不确定度为,【例】,某实验室事先对某一电流量进行,n,10,次重复测量，测量值列于下表。按下表的计算步骤得到单次测量的估计标准偏差,s,(,x,),0.074mA,。,在同一系统中在以后做单次（,m,1,）测量，测量值,x,46.3mA,，求这次测量的标准不确定度,u,(,x,)。,在同一系统中在以后做,3,（,m,3,）次测量，,mA,，求这次测量的标,准不确定度,。,表 对某一电流量进行,n,10次重复测量的测量值,次数,i,1,2,3,4,5,测量值,mA,46.4,46.5,46.4,46.3,46.5,次数,i,6,7,8,9,10,测量值,mA,46.3,46.3,46.4,46.4,46.4,平均值,46.39mA,单次测量的,标准偏差,s,(,x,),0.074mA,【解】,对于单次测量，则其标准不确定度等于1倍单次测量的标准偏差：,x,46.3mA,u,(,x,)=,s,(,x,)=0.074mA,【解】,对于,m,3测量，测量结果为：,的标准不确定度为：,标准不确定度的B类评定,B类评定的方法是根据有关的信息或经验，判断被测量的可能值区间,x,-,a,x,+,a,假设被测量的概率分布，根据概率分布和要求的概率,p,确定,k,，则B类标准不确定度 可由下式得到：,a,-被测量可能值区间的半宽度,k,-包含因子,预备知识,分布数据散布的“形状”,一组数值的散布会取不同的形式，或称为服从不同的概率分布。,（,1,）,正态分布,在一组读数中，较多的读数值靠近平均值，少数读数值离平均值较远。这就是正态分布或高斯分布的特征。,（,2,）,均匀分布（矩形分布）,当测量值非常平均地散布在最大值和最小值之间的范围内时，就产生了矩形分布或称为均匀分布。,（,3,）,其他分布,还有其他分布形状，但较少见，例如三角分布、反余弦分布（,U,型分布）等。表,2.1,给出了几种概率分布及其包含因子。,2,a,(=,a,),x,矩形（均匀）分布,矩形(均匀)分布,标准不确定度(标准偏差)：,特征：估计值以,p,100的概率均匀散布在,a,区间内，落在该区间外的概率为零。,B类评定方法求得的标准不确定度也是标准偏差,P(%),68.27,90,95,95.45,99,99.73,K,1,1.65,1.96,2,2.58,3,正态分布,三角分布,标准不确定度(标准偏差),：,特征,：,估计值以,p,100的概率落在,a,区间内，靠近,x,的数值比接近边界的值多，落在该区间外的概率为零。,x,-,a,x,2,a,(=,a,),三角分布,x,+,a,标准不确定度的B类评定,概率分布按以下不同情况假设：,a）被测量受许多随机影响量的影响，当它们各自的效应同等量级时，不论各影响量的概率分布是什么形式，被测量的随机变化近似正态分布。,b）,如果有证书或报告给出的不确定度是具有包含概率为0.95、0.99的扩展不确定度,U,，此时除非另有说明，可按正态分布来评定。,标准不确定度的B类评定,概率分布按以下不同情况假设：,c）当利用有关信息或经验估计出被测量可能值区间的上限和下限，其值在区间外的可能几乎为零时，若被测量落在该区间内的任意值处的可能性相同，则可假设为均匀分布；若被测量落在该区间内中心的可能性最大，则假设为三角分布；若落在该区间中心的可能性最小，而落在该区间上限和下限的可能性最大，则可假设为反正弦分布。,标准不确定度的B类评定,概率分布按以下不同情况假设：,d）已知被测量的分布是两个不同大小的均匀分布合成时，则可假设为梯形分布。,e）,对被测量的可能值落在区间内的情况缺乏了解时，一般假设为均匀(矩形)分布。,f）实际工作中，可依据同行专家的研究结果或经验来假设概率分布。,【,例】,由校准证书数据确定标准不确定度,用于校准的某台设备的校准证书中说明，在它的校准范围内的测量不确定度为,U,(,x,)0.5V，,k,=2。在其校准范围内，由该设备引起的标准不确定度,u,(,x,)为：,（V）,【例】,由准确度误差限数据确定标准不确定度,查JJG596-1999 电子式电能表,检定规程，,0.2级,A型,三相（平衡负载）标准电能表，基本误差限为,0.2。则区间半宽度为,a,=0.2，服从矩形分布，包含因子 。,由此引起的标准不确定度为：,【例】,数字显示的分辩力引起的标准不确 定度,如果仪器的分辩力为,，则示值,x,将等概率处于,x,/2，,x,/2的区间任何位置。即被测量的可能值服从均匀分布，区间半宽度,a,/2，包含因子 。,因此分辩力引起的标准不确定度为：,1、不管是A类评定还是B类评定，都是指评定不确,定度的方法，即求得各个分量的标准不确定度,时所用的方法，并非指不确定度的分类。,2、求每一个分量的不确定度时，注意其单位的统,一。不管是每个分量的标准不确定度，还是合,成标准不确定度和扩展不确定度，都是有单位,的，且必须有一致性。,需要注意的一些问题,预评估重复性,每一次测量都会评定一次测量不确定度？,在日常开展同一类被测件的常规检定、校准或检测工作中，如果测量系统稳定，测量重复性无明显变化，则可用该测量系统以与被测件相同的测量程序、操作者、操作条件和地点，预先对典型的被测件的典型被测量值进行n次测量（一般n不小于10），由贝塞尔公式计算出单个测得值的实验标准偏差 ，即测量重复性。,预评估重复性,预评估重复性,在对某个被测件实际测量时可以只测量 次（），并以 次独立测量的算术平均值作为被测量的估计值，则该被测量估计值由于重复性导致的A类标准不确定度按下式计算：,用这种方法评定的标准不确定度的自由度仍,v=n-,1,。,应注意，当怀疑测量重复性有变化时，应及时重新测量和计算实验标准偏差 。,预评估重复性,目前，我们采用的都是预评估重复性，但是，有一个问题需要思考。常见的校准项目，比如数字多用表，温湿度计，数显卡尺等，校准的数量都是比较大的，但我们在评定测量不确定度时，基本都是选取一个进行预评估重复性。,这样，取样的代表性够么？,合并样本标准偏差,两种状况,一、同一被测品采用核查标准和控制图的方法使测量过程处于统计控制状态。如数字多用表连续测量30天，每天测10个数据（或每天测不同个数据，如5月1号得到10个数据，5月2号得到8个数据，5月3号得到9个数据）。,二、同一类型的多个被测品进行测量。如实验室有30个同一型号数字多用表，每一个数字多用表测10个数据（或每个样品测不同个数据，如第1个得到10个数据，第2个得到6个数据，第3个得到7个数据）。,合成标准不确定度,影响测量结果不确定度的因素很多，为了计算总不,确定度，需要将各不确定度分量进行合成。在计算合成,标准不确定度之前，,需要确定各输入量的标准不确定度,是否彼此相关，如相关，则需要考虑其带来的影响。,对,于大多数情况，我们可以认为输入量的标准不确定度是,彼此互不相关的（,本PPT仅考虑输入量不相关的况,）。,合成标准不确定度,各输入量不相关时合成标准不确定度的计算,当输入量不相关时，相关系数为零。被测量的估计值（测得值）,y,的合成标准不确定度,u,c,(,y,)按下式计算：,公式中的,u,(,x,i,)即各个分量的标准不确定度。,灵敏系数,在求出各个输入量的不确定度分量,u,i,(,x,)之后，还需要计算灵敏系数,c,i,：,每一个分量的不确定度对合成标准不确定度的贡献（传播）并不是均等的，可能有的分量对合成标准不确定度影响大，有的影响小，这样，当在合成时，就需要一个加权系数，也就是灵敏系数。,灵敏系数,某公司，甲乙丙丁四人分别持公司股份的1/2、1/4、1/8、1/8，某种意义上讲，甲的权力最大，乙丙丁分别次之；也可以说，甲的决议，对公司的运作影响最大，乙丙丁分别次之。,也就是说，甲乙丙丁的行为都会对公司产生影响，即便是同样一个行为，影响的程度都不一致，而多大的程度，就可以理解为加权系数，即灵敏系数。,甲,乙,丙,丁,灵敏系数,灵敏系数,测量模型：,R,=R-U/I,例：,使用多功能校准器FLUKE 5720A和数字多用表FLUKE 8508A来校准电阻值1m的电阻，如下图，即用FLUKE 5720A输出电流I通过被校电阻R，同时用FLUKE 8508A在被校电阻R两端测得电压U，则实测电阻值=U/I，已知电阻误差公式为标称值-实测值，试建立电阻测量不确定度评估测量模型，并给出各不确定度分量的灵敏系数。,电阻R,FLUKE 5720A输出,电流I,FLUKE 8508A 测得电压U,灵敏系数：,扩展不确定度,扩展不确定度是测量可能值包含区间的半宽度。扩展不确定度分为,U,和,U,p,两种。在给出测量结果时，一般情况下报告扩展不确定度,U,。,1、扩展不确定度,U,2、扩展不确定度,U,p,扩展不确定度,1、扩展不确定度,U,扩展不确定度,U,由合成标准不确定度,u,c,乘包含因子,k,得到，按下式计算：,U,=,ku,c,通常结果可用下式表示：,Y,=y,U,y,是被测量Y的估计值，被测量Y的可能值以较高的包含概率落在,y,-,U,y,+,U,区间内，即,y,-,U,Y,y,+,U,。被测值的值落在包含区间内的包含概率取决于所取的包含因子,k,的值，,k,的值一般取2或3。,扩展不确定度,1、扩展不确定度,U,当y和,u,c,(,y,)所表征的概率分布近似为正态分布时，且,u,c,(,y,)的有效自由度较大情况下，若,k,=2，则由,U,=2,u,c,所确定的区间具有的包含概率约为95%。若,k,=3，则由,U,=3u,c,所确定的区间具有的包含概率约为99%。,在通常的测量中，一般取,k,=2。当给出扩展不确定度,U,时，一般应注明所取的,k,值。,扩展不确定度,2、扩展不确定度,U,p,当要求扩展不确定度所确定的区间具有接近于规定的包含概率,p,时，扩展不确定度用符号,U,p,表示，当,p,为0.95或0.99时，分别表示为,U,95,和,U,99,。,U,p,由下式获得：,U,p,=,k,p,u,c,k,p,是包含概率为,p,时的包含因子，由下式获得：,k,p,=,t,p,(,v,eff,),根据合成标准不确定度的有效自由度和需要的包含概率，查t分布表得到。,t分布表（节选）,测量不确定度的表示,测量不确定度的表示,1、合成标准不确定度,u,c,(,y,)的报告可用以下三种形式之一。,例如，标准砝码的质量为,m,s,，被测量的估计值为100.02147g，合成标准不确定度,u,c,(,m,s,)=0.35mg,则报告为：,测量不确定度的表示,a),m,s,=100.02147g;合成标准不确定度,u,c,(,m,s,)=0.35mg。,b),m,s,=100.02147(35)g;括号内的数是合成标准不确定度的值，其末位与前面结果内末位数对齐。,c),m,s,=100.02147(0.00035)g;括号内是合成标准不确定度的值，与前面结果有相同计量单位。,测量不确定度的表示,2、当用扩展不确定度报告测量结果的不确定度时，,U,=,ku,c,(,y,)的报告可用以下四种形式之一。,例如，标准砝码的质量为,m,s,，被测量的估计值为100.02147g，合成标准不确定度,u,c,(,y,)=0.35mg,取包含因子,k,=2，,U,=20.35mg=0.70mg,则报告为：,测量不确定度的表示,a),m,s,=100.02147g，,U,=0.70mg；,k,=2。,b),m,s,=(100.021470.00070)g；,k,=2。,c),m,s,=100.02147(70)g；括号内为,k,=2的,U,值，其末位与前面结果内末位数对齐。,d),m,s,=100.02147(0.00070)g,；括号内为,k,=2的,U,值，与前面结果有相同计量单位。,测量不确定度的表示,3、当用扩展不确定度报告测量结果的不确定度时，,U,p,=,k,p,u,c,(,y,)的报告可用以下四种形式之一。,例如，标准砝码的质量为,m,s,，被测量的估计值为100.02147g，合成标准不确定度,u,c,(,y,)=0.35mg，,v,eff,=9，按,p,=95%，查表得,k,p,=t,95,(9)=2.26，,U,95,=2.260.35mg=0.79mg，则报告为：,测量不确定度的表示,a),m,s,=100.02147g，,U,95,=0.79mg，,v,eff,=9。,b),m,s,=(100.021470.00079)g，,v,eff,=9，括号内第二项为,U,95,之值。,c),m,s,=100.02147(79)g，,v,eff,=9，括号内为,U,95,之值，其末位与前面结果内末位数对齐,。,d),m,s,=100.02147(0.00079)g,，,v,eff,=9，括号内为,U,95,之值，与前面结果有相同计量单位。,三、,声明检测或校准结果与规范符合性,如何作出判定，需要考虑测量不确定度么？,三、,声明检测或校准结果与规范符合性,CNAS条文,5.3 在多数情况下，规范要求在证书或报告中做出符合性声明，但没有指 明进行符合性评价时需考虑不确定度的影响。在这种情况下，用户可以在不考虑不确定度的情况下，根据测得值是否在规定限值范围内做出符合性判断。（CNAS-GL27）,这就是通常所说的“风险共担”，通过使用约定的测量方法进行测试后的产品还有可能不符合规范要求，最终用户可能会承担此风险。在这种情况下，一般假设约定的测量方法的不确定度是可以接受的，而且重要的是其不确定度在必要时是可以评估的。国家法规可以否决“风险共担”的原则，并使不确定度引起的风险由其中一方承担。,三、,声明检测或校准结果与规范符合性,CNAS条文,5.4 在用户与实验室之间的协议或实施准则或规范中可能已声明其采用方法的准确性是足够的，且在判断符合性时，可以忽略不确定度。上文中对“风险共担”的考虑也适用于此情况。（CNAS-GL27）,三、,声明检测或校准结果与规范符合性,三、,声明检测或校准结果与规范符合性,END,谢谢各位！,