第八讲 扩展不确定度的计算
扩展不确定度

160 161
测量仪器名称
*酸度计 *可见分光光度计
校准参量
酸度 波长 透射比
校准规范
测量范围
扩展不确定度 说明 (k=2)
U=0.02pH U=0.6nm U=0.5% U=0.5nm U=0.5% Urel=0.16% Urel=0.3% Urel=1.6% 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效 有效
TCD:≥ 800mV•mL/mg
-10
Urel=4.2%
FID:≤5×10-10 g/s Urel=4.5% FPD:≤5×10 g/s(硫) ≤1×10 10 g/s(磷) NPD:≤5×10-12 g/s(氮) ≤1×10 11 g/s(磷) ECD:≤5×10-12 g/mL 176 177 178 179 *化学需氧量(COD) 浓度 测定仪 *离子计 *原子吸收分光光 度计 *原子荧光光度计 离子浓度 化学需氧量(COD)测定 (0.1~ 仪检定规程 JJG 975 1000)mg/L 离子计检定规程 JJG 757 原子吸收分光光度计 检定规程 JJG 694 原子荧光光度计检定 规程 JJG 939 (0~14)pX Cu:≤0.02μ g/mL Cd:≤4pg As:≤0.4ng Sb:≤0.4ng Urel=4.5%
有效
175
*气相色谱仪 检出限
气相色谱仪检定规程 JJG 700
Urel=4.2%
有效
Urel=4.7% U=2.8mg/L U=0.02pX U=0.005μg/mL U=0.4pg U=0.06ng U=0.06ng
有效 有效 有效 有效 有效 有效 有效
应用回归直线建立测长机扩展不确定度的估算公式

AE。 A 为 即期 望 E
( 的估计值 ; = E( 为变量 的中心化变 量 ;E ) ~ — ) Aa
( = ^ ) E( 为E ( 的期望 估 计 误差 。期 望 估 计误 ) E ( 一 ) )
由于未 定 系统误 差 的确 定 值是 未 知 的 , 因此 , 它 的评 估 对 实 际上是 对期望 估计 方法 可能存 在 的误差 进行评 估 。 以 所
【 随机变量的不确定度ucr itor dm v i l netn no r be a y fa aa
随机 变量 的不确 定 度是 表 征变 量 随机 部分 u 大小 ^ 的统计 特征 估计 值 。
维普资讯
T C N L GY S C I 技术篇 E H O O E TON
误差与不确 定度
术语 “ 不确定度"定义的剖析( ) 下
口王春艳 陆梅 高蔚 钱钟 泰
三、 随机变 量 的不确 定度 和测 量不 确定度
1 意 随机 变量 的“ 确 定度” . 任 不
/0 ( 10 mm)
5 0
相 同测 量 条 件 下 , 同被测 尺 寸 的扩展 不 确 定 度计 不
算值 , 列于表 2 。
0l 2 5 0
二、 回归 直 线 的 建立
被测件与测 长 机的分米尺热 B 膨胀系数差
2×1 0
三角
L t 4 × O 25 d A
5 0
在被测件 和测 长机充分恒温 , 温差不大于03 室 _ ℃, 温偏离2 ℃不超过土 ℃. 0 1 室温变化每小时不大于0 ℃等 _ 3 标准条件下 . 在测长机上对尺寸为3 0 m的被测钢件进 5m
行 直 接 测量 . 量 结果 的 不 确定 度 计 算步 骤列 于表 1 测 。
测量结果的扩展不确定度-解释

在理想状态下装置输出电能并测出其值 Wo,若被检表测量值为 Wx 则由计算机按
e机显 WX WO 100(%) Wo
计算出被检表误差并显示出来,作为实测结果. 在实际工作中,考虑到各种实际因素的影响,则有
e机显 W X WN e'CL WN
e结果 ex et ex (T , F , Q, R) [e0 e0 (T , F , Q, R)] ecl ' '
式中 e 机显—装置计算机计算并显示的被试表误差 e 结果—被试表误差的测量结果 WN—装置测出的电能值 Wx—被试表测出的电能值 et—时间控制误差 ex—被试表误差 eo—标准装置误差 ex(T,F,Q,R)—被试表由于受到温度、频率、功率稳定度及其它原因影响引入的误差 eo(T,F,Q,R)—标准装置由于受到温度、频率、功率稳定度及其它原因影响引入的误 差
(装置显示误差小数点后有 3 位) 5.4.2 结果化整引入不确定度分量 结果化整间距与被试表的级别有关,1.0 级被试表的化整间距为 0.1%,故:
U (ecl ) 2 0.05 3 0.02887 (%)
U ( ecl)1 0.005
1 结果化整带入的误差限为被试表误差限的20 5.5 被试表受影响量变化引入的不确定度分量 5.5.1 受温度变化影响引入的不确定度分量 根据对大量被试表数据和有关技术资料的分析,1 级表的温度系数小于 0.05%/℃, 实验室温度在(20±2)℃变化,且呈均匀分布,故:
根据
x i x
i 1
n
2
n 1
σ1.0=0.007(%) σ0.5L=0.008(%) 再由:n=10,得
u
(整理)不确定度的计算方法.

(整理)不确定度的计算⽅法.精品⽂档测量结果的正确表达被测量X 的测量结果应表达为:)(单位U X X ±= 其中X 是测量值的平均值,U 是不确定度。
例如:⽤最⼩刻度为cm 的直尺测量⼀长度最终结果为:L =(0.750±0.005)cm ;测量⾦属丝杨⽒模量的最终结果为:E = (1.15±0.07)×1011Pa 。
1. 不确定度的计算⽅法直接测量不确定度的计算⽅法22仪?+=S U其中: 1)(2--=∑n X XS i为标准差;仪?是仪器误差,⼀般按仪器最⼩分度的⼀半计算,但是游标卡尺和⾓游标按最⼩分度计算。
也可按仪器级别计算或查表。
间接测量不确定度的合成⽅法间接测量)??=,,,(z y x f N 的平均值公式为:)??=,,,(z y x f N ;不确定度合成公式为: +++=222222)()()(Z Y X N U ZN U Y N U X N U 。
也可根据表1中的公式计算间接测量的不确定度。
表1 常⽤函数不确定度合成公式函数表达式合成公式2γβαZY X N =222222)()()(ZUY U X U N U Z Y X N γβα++= 注:1. 在函数关系是乘除法时,先计算相对不确定度(NU N)⽐较⽅便.例如表中第⼆⾏的公式.2. 不确定度合成公式可以联合使⽤.例如:若φθτ3sin =,令θsin =u ,φ3=w 则wu=τ.精品⽂档根据表中第⼆⾏公式,有:22)()(wUu U U w u +=ττ; 根据表中第⼀⾏公式,有: φφU U U w 332 2==; 根据表中第三⾏公式,有: θθU U u ?=cos . 所以, 2222)( )sin cos ()33()sin cos (φθθτφθθτφθφθτU U U U U +??=+??=。
不确定度计算

不确定度计算2、不确定度各分量的评定根据测量步骤可知,测量氨氮质量的不确定度来源有几个方面,一是由标准曲线配制所产生的不确定度,二是测试过程所产生的不确定度。
按《化学分析中不确定度的评估指南》,对于只涉及积或商的模型,例如:c N=m/v,合成标准不确定度为:式中,u(c)为质量m和体积v的合成标准测量不确定度,mg/L;u(m)为质量m的标准测量不确定度,ug;u(v)为体积v的标准测量不确定度,mL。
2.1 取样体积引入的相对不确定度u rel(v)所取水样用50mL单标线吸管移取。
查JJG 196-2006《常用玻璃量器检定规程》,A级50mL 单标线吸管的容量允差为0.05mL,根据JJF 1059-1999《测量不确定度评定与表示》的规定,标定体积为三角分布,则容量允差引入的不确定度为:u(△V)=0.050/√6 。
根据制造商提供的信息,吸量管校准温度为20℃,设实验室内温度控制在±5℃范围内波动,与校准时的温差为5℃,由膨胀系数(以水的膨胀系数计算)为2.1×10-4/℃得到50mL水样的标准不确定度为(假定为均匀分布):2.2重复性测定引入的相对不确定度u rel(rep)采用A类方法评定,与重复性有关的合成标准不确定度均包含其中。
对某水样进行7次重复性测定,所得结果如下:1.33、1.35、1.34、1.34、1.35、1.38、1.35mg/L,平均值1.35 mg/L。
重复测量数据的标准不确定度为:2.3 铵(以氮计)的绝对量m引入的不确定度u rel(m)配制过程中引入的不确定度u rel(1)a.) 标准贮备液的不确定度u rel(1-1):包括纯度、称量、体积及摩尔质量计算4个部分,其中,摩尔质量计算不确定度可省略不计(与其它因素相比,其对标准浓度计算相差1-2个数量级)。
纯度p:按供应商提供的参考数据,分析纯氯化铵[NH4Cl]纯度为≥99.5%,将该不确定度视为矩形分布,则标准不确定度为u(p) =0.5/√3=28.9×10-4;称量m:经检定合格的天平最大允许误差±0.1mg,将该不确定度视为矩形分布,标准偏差为0.058mg,称量3.819g时的相对标准偏差为u(m) =0.152×10-4;体积v:影响体积的主要不确定度有校准及温度。
扩展不确定度包含因子的简单计算

·46·质量技术监督研究Quality and Technical Supervision Research2010年第1期(总第7期)NO.1.2010General No.7一、引言在测量不确定度评定过程,需要选择包含因子k或kp,在JJF 1059-1999《测量不确定度评定与表示》(以下简称“JJF 1059”)中7.1条,对计算扩展不确定度有两种方法:一种直接选择k=2或3,乘合成标准不确定度uc(y)得扩展不确定度U;另一种在给定概率p下,查“t分布临界值”的kp,然后乘合成标准不确定度uc(y)得扩展不确定度U。
但没有明确在什么情况下选择那种包含因子评定扩展不确定度,只提到“但veff充分大时,可以近似认为k95=2、k99=3”。
因此本文主要根据ISO、IEC等7个国际组织联合制定并发布《测量不确定度表达指南》〔(GUM)ISO:1993(E)〕有关内容阐述如何简单选取包含因子计算扩展不确定度。
二、对《测量不确定度表达指南》中推荐的简单方法的理解根据《测量不确定度表达指南》附录G选取包含因子简单选择方法有:在Y和uc(y)表征的概率分布近似为状态分布,且uc(y)的有效自由度veff较大时,可取k=2,其所形成的区间±2uc(y)具有置信概率约为95%;若取k=3,其所形成的区间±3uc(y)具有置信概率约为99%。
具体成立条件为:1、被测量Y的估计值y,是由适当的输入量Xi 的估计值xi 得到的,Xi 由良好概率分布描述,如正态分布或矩形分布;2、可由A类或B类评定输入量的估计值的标准不确定度u(xi ),它们对测量结果合成不确定度uc(y)有一定的贡献;3、由不确定度传播律要求的线性近似是合适的;4、uc(y)的不确定度是很小的,因其有效自由度veff具有足够的大小,如veff>10。
实际评定测量不确定度对上述4个条件中大部分满足前3条。
不确定度的计算方法

测量结果的正确表达
被测量X的测量结果应表达为:'- !
其中是测量值的平均值,是不确定度。
例如:用最小刻度为cm的直尺测量一长度最终结果为:L=(0.750 ±.005cm;测量金属丝杨氏模量的最终结果为:E=(1.15 0.07 >1011Pa。
1.不确定度的计算方法
直接测量不确定度的计算方法
; ------ _
其中:为标准差;
是仪器误差,一般按仪器最小分度的一半计算,但是游标卡尺和角游标按最小分度计算。
也可按仪器级别计算或查表。
间接测量不确定度的合成方法
间接测量’的平均值公式为:•;
不确定度合成公式为:。
也可根据表1中的公式计算间接测量的不确定度。
表1常用函数不确定度合成公式
1.在函数关系是乘除法时,先计算相对不确定度(>比较方便.例如表中第二行的公式•
2.不确定度合成公式可以联合使用.
sin & u
例如:若,令,则根据表中第二行公式,有:;
U = iVu1 = W7
根据表中第一行公式,有:;根据表中第三行公式,有:•
U t 所以,
普八(与
sin 3 0
皿i
1」
1sin0 0
U
、
U^= cos X U。
不确定度的扩展与合成

uC4 N ui4
(1)
i1 i
N
且 eff i i 1
eff
E(uC )4 N Pi E(ui )4
(2)
i 1
i
相对不确定度
E(uC ) uC y
E(ui ) ui xi
B类不确定度的有效自由度
uB
5%
相对不确定度
有效自由度vB
200
10% 50
20% 12
25% 8
Page . 5
【例题1】产品性能的检定
某工厂产品上需要使用1000的精密电阻,
设计要求其允许误差极限在土0.
2%
以内。为此,
用一台4
1 2
数字多用表直接测量随机抽取的10个电阻,请报告测量结果。该仪表的技术指标如下:
最大允许误差=(0.01%*读数+1个LSB)
Page . 4
B类不确定度uB的有效自由度可按下式来近似计算
1 ������ 2
������ ������B ������B
−2
有效自由度的计算结果应取整数;
式中,
[uB]代表不确定度uB的不确定度,������
������B ������B
则是不确定度uB的相对不确定度;
常用的相对不确定度取值及其有效自由度如表所示:
N
uC ( y) ciu(xi )2 i 1
传播系数ci
f xi
Page . 10
【例题2】间接测量中不确定度的评定举例
已知电阻 R1=50.2±0.5(Ω), R2 =149.8±0.5(Ω), 求它们串联的电阻R及其合成不确定度。
解:串联电阻的阻值为 合成不确定度
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第八讲扩展不确定度的计算
减小字体增大字体作者:李慎安来源: 发布时间:2007-05-08 10:33:45
计量培训:测量不确定度表述讲座
国家质量技术监督局李慎安
8.1 什么叫扩展不确定度?
按《JJF1001》扩展不确定度定义为:确定测量结果区间的量,合理赋予被测量之值分布的大部分可望含于此区间。
也称展伸不确定度或范围不确定度。
符号为大写斜体U,U P。
当除以被测量之值后,称为相对扩展不确定度,符号为U rel,U prel。
符号中的p为置信概率,一般取95%,99%,这时其符号成为U95,U99,U95rel或U99rel。
定义中所指大部分,最常用的是95%和99%。
扩展不确定度过去曾称总不确定度(overall uncertainty),这一名称已为《导则》所禁止使用,因其从含义上易与合成不确定度混淆。
扩展不确定度是比合成标准不确定度大的一个参数,它等于合成标准不确定度乘以包含因子k后的值,对于合成标准不确定度而言,它是成倍地被扩大了的一个值。
8.2 扩展不确定度分成几种?
扩展不确定度根据所乘的包含因子k的不同,分成两大类。
当包含因子k之值取2或3时,扩展不确定度U只是合成标准不确定度u C的k倍。
在给出U时,必须指明k的取值。
实际上,这时的U所包含的信息与u C一样,并未因乘以k后,其信息有所增多。
此外,还有一种包含因子k p,它是为了使扩展不确定度所给出的区间内能有概率为p的合理赋予被测量之值含于其中所必须有的因子。
所得到的扩展不确定度为U p。
一般,只在被测量Y可能值y的分布类型可估计为正态时才给出U P。
这时的k p之值,按u c(y)的有效自由度υeff,通过本讲座6.6中的表得出,即t p值,k p=t p(υ)。
随υ的增大,k有所降低,随p的增大,k p有所增加。
与上述类似,相对扩展不确定度亦有两种。
8.3 什么情况下使用U,什么情况下使用U p来说明测量结果的不确定度?
(1)根据有关测量仪器校准的技术规范。
例如,以下技术规范规定取k=3,JJF2002,2003,2004,2018,2019,2025,2026,2030,2032~2041,2045,2446等,不一一例举。
而以下技术规范规定取k=2,JJF2049,2050,2072,2089等。
也有一些技术规范规定用U95,如JJF2006,2061,等。
规定采用U99的如JJF2020,2056,146等。
(2)可以估计被测量Y估计值y之分布接近正态时,可给出U p,否则只能给出U。
8.4 什么情况下可用包含因子k95=2及k99=3?
如果y的分布是比较理想的正态分布,那么,当合成标准不确定度u C(y)的有效自由度充分大时,即可做出这样较简单的处理,例如,在p=95%时,自由度为12,这时,按本讲座6.6,k p=2.18,如取k p=2,其值小了不到十分之一,应该说就无足轻重了。
当p=99%时,υeff无穷大的k p=2.58≈2.6,整化为k99=3,已较保守;而当υeff=20时,k99之值为2.85,它比2.6大约大十分之一,因此,这时如不用2.85而用2.6,所得U99也只小十分之一左右,应可忽略。
因此,在《JJF1059》中所要求的有效自由度应充分大,拿十分之一作为可忽略的标准,则对于p=95%时,υeff应大于12,对于p=99%,应大于20。
8.5 什么情况下,虽未计算合成标准不确定度u c(y)的有效自由度,取包含因子k=2给出的扩展不确定度U可以估计是置信区间在p=95%的半宽,可否在检定证书中给出其值为U95?
虽未算出υeff,但其值估计不太小,例如,大于12,而且,可以估计Y的估计值的分布接近正态,这时,一般可以认为U=2u c(y)的置信概率p大约为95%。
但是不能在证书上给出其值为U95之值。
对于取k=3,认为U=3u c(y)近似地有p=99%亦类似。
可作如此估计,但不能在证书上给出。
8.6 给出校准测量能力时,包含因子k应取多少?
校准测量能力(calibration mea-surement capability)定义为:通常提供给用户的最高校准测量水平,规定用包含因子k=2的扩展不确定度表示。
校准测量能力有时又称最佳测量能力(best measurement capa-bility)。
在国际比对中,有时也特别指明采用k=2的扩展不确定度。
8.7 已确知被测量Y可能值y的分布并非正态分布而是某种其他分布,例如:三角分布、梯形分布、均匀分布、两点分布等较为典型而且规则的分布时,可否给出U=ku c(y)或U p?
当我们根据经验,一般可以按本讲座6.8对被测量Y可能值做出其分布的评定。
例如说,我们完全可以认定是均匀分布,那么,取包含因子k=2或k=3是不合理的。
例如,某个被测量Y可能值y的分布,主要决定于某一个均匀分布的输入量,其他输入量的影响可忽略不计时,我们按这个主要的输入量的最大
允许误差,按均匀分布取了其标准偏差,即a/k=≈a/1.7。
然后又取k=2得出一个扩展不确定度2×(a/1.7)≈1.2a(关于a的含义见本讲座6.7),这比原来的分散区间半宽a(p=100%)还大了约五分之一。
如果取k=3,得到的扩展不确定度U=3×≈1.8a,(即较a大了约五分之四。
十分不合理,导致误解。
如果是两点分布,情况就更糟。
当然,三角分布比较接近正态分布,情况会好些。
因此,在可以
确定Y可能值的分布接近某种其他分布时,不应取k=2或k=3。
对于均匀分布来说,对于U95,k p=1.65;对于U99,k p=1.7。
本讲座6.7中给出的k值,近似地为相应分布的U99之值。
8.8 标准测量仪器(或校准装置)的扩展不确定度中,是否应包含被校准测量仪器的示值重复性?
测量仪器的重复性定义为在相同测量条件下,重复测量同一个被测量,测量仪器提供相近示值的能力。
这些条件即重复性条件。
而测量仪器的重复性是用示值分散性定量表述的,即重复性标准偏差。
当我们用校准设备对其进行校准时,一般来说,校准设备的重复性标准偏差大大地小于被校准仪器的重复性标准偏差。
校准过程中出现的重复性标准偏差主要是被校准仪器的。
因此,这一部分不应包含在标准测量仪器的扩展不确定度之中。
8.9 单侧检验中,扩展不确定度U95应如何计算?
在统计检验中,当量是一维的情况时,以小于(或大于)某一给定值的所有值的集合,作为拒绝域的检验称为单侧检验。
计量学中,通过测量以确定被测量的真值以给定置信概率p(p一般取95%,偶也有99%或90%),不大于(或不小于)某值的检验亦称为单侧检验。
例如,通过测量得出样品中砷的质量浓度以95%的概率不大于0.1mg/L,这就是一种单侧检验。
当一个被校准的砝码质量m经校准后表达为,例如:
m=100.006g±0.004g或100.002g≤m≤100.010g,很明显,m的可能值有上、下两个界,即100.002g与100.010g,这就是所谓双侧检验。
得出真值以某个置信概率处于某两极限值之中,而单侧检验则只有一个上界或一个下界。
上例中的0.1mg/L即为一个上界而不存在下界。
单侧检验与双侧检验在合成标准不确定度评定中的方法相同,只是给出扩展不确定度(p=95%)所乘的包含因子k95,与双侧检验不同,即不是查本讲座6.6中的t值表而是下表:
当计算出被测量Y的扩展不确定度U95及其最佳估计值y之后,其单侧的上界(或称单侧置信上限): T2>y-U95
其下界(或称单侧置信下限):
T1<y+U95
根据需要,上述T1与T2计算中的U95可用U99代替。
ISO2602《检测数据的统计处理:平均值的估算、置信区间》中给出的单侧置信区间及上、下限的计算方法,由于只考虑了测量在重复性条件下的分散性,没有考虑其他不确定度分量,不应再作为评定的依据。
(未完)。