不确定度评估的基本方法

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

三、检测和校准实验室不确定度评估的基本方法
1、测量过程描述:
通过对测量过程的描述,找出不确定度的来源。

内容包括:测量内容;测量环境条件;测量标准;被测对象;测量方法;评定结果的使用。

不确定度来源:
● 对被测量的定义不完整; ● 实现被测量的测量方法不理想;
● 抽样的代表性不够,即被测样本不能代表所定义的被测量;
● 对测量过程受环境影响的认识不周全,或对环境的测量与控制不完善; ● 对模拟式仪器的读数存在人为偏移;
● 测量仪器的计量性能(如灵敏度、鉴别力、分辨力、死区及稳定性等)的局限性; ● 测量标准或标准物质的不确定度;
● 引用的数据或其他参量(常量)的不确定度; ● 测量方法和测量程序的近似性和假设性; ● 在相同条件下被测量在重复观测中的变化。

2、建立数学模型:
建立数学模型也称为测量模型化,根据被测量的定义和测量方案,确立被测量与有关量之间的函数关系。

● 被测量Y 和所有个影响量i X ),2,1(n i ,⋯=间的函数关系,一般可写为
),2,1(n
X X X f Y ,⋯=。

● 若被测量Y 的估计值为y ,输入量i X 的估计值为i x ,则有),x ,,x f(x y n ⋯=
21。

有时为简化
起见,常直接将该式作为数学模型,用输入量的估计值和输出量的估计值代替输入量和输出量。

● 建立数学模型时,应说明数学模型中各个量的含义。

● 当测量过程复杂,测量步骤和影响因素较多,不容易写成一个完整的数学模型时,可以分步评定。

● 数学模型应满足以下条件:
1) 数学模型应包含对测量不确定度有显著影响的全部输入量,做到不遗漏。

2) 不重复计算不确定度分量。

3) 选取合适的输入量,以避免处理较麻烦的相关性。

● 一般根据测量原理导出初步的数学模型,然后将遗漏的输入量补充,逐步完善。

3、不确定度的A 类评定:
(1)基本方法——贝塞尔公式(实验标准差)方法
在重复性条件下对被测量X 做n 次独立重复测量,得到的测量结果为
i x )
,2,1(n i ,⋯=。

则X 的
最佳估计值可以用n 次独立测量结果的算术平均值来表示:
n x
x n
i i
∑==
1。

根据定义,用标准差表示的不确定度为标准不确定度。

于是单次测量结果的标准不确定度可用贝塞尔公式表示:
若在实际工作中,采用n 次测量结果的算数平均值作为测量结果的最佳估计值,则平均值的标准不确定度为:2()
1
1
=()(1)n
A Y i i u S Y Y n ==--∑()。

◆ )(i x u 和)(x u 的自由度都为1-n 。

◆ 显然,采用m 次测量结果的算数平均值作为测量结果的最佳估计值,比单次测量结果更可靠,因此,算术平均值的标准不确定度(实验标准差)比单次测量结果的标准不确定度(实验标准差)小。

◆ 在使用贝塞尔公式时,要求n 应比较大。

JJF1033—《计量标准考核规范》中规定,在进行计量标准的重复性测量时,要求测量次数n ≥10。

◆ 如果通过n 次重复测量得到的单次测量结果的标准不确定度(实验标准差),可以保持相当长时间不变,若出现测量结果是m (m 可能比较小)次重复测量的算术平均值,则该平均值的标准不确定度(实验标准差)为:2()
11
=()(1)n
Y A i i S u Y Y m n m
==--∑()。

(2)合并样本标准差)(i P x s 方法
若在实际工作中,在重复性条件下,对被测量X 做n 次独立测量,并有k 组这样的测量结果。

由于各组之间的测量条件可能会稍有不同,因此不能直接用贝塞尔公式对总共k n ⨯次测量计算标准不确定度(实验标准差),而必须使用合并样本标准差)(i P x s ,公式可表示为:
2()()
11
1()(1)k k k n
j Y P Y ji j i u s Y Y k n ====--∑∑ 式中ji x 是第j 组的第i 次测量结果,j x 是第j 组的n 个测量结果的算术平均值。

◆ 合并样本标准差也称为组合实验标准差。

◆ 若已分别算出k 组测量结果的实验标准差)(i j x s ,而且每组包含的测量次数相同,合并样本
标准差)(i P x s 可表示为:k
x s x s k
j i j i P ∑==
1
2
)
()(。

◆ 合并样本标准差)(i P x s 应该采用方差的平均值,即合并样本方差)(2
i P x s 等于各组样本方差
)(2
i j x s 的平均值。

◆ 若各组所包含的测量次数不完全相同,合并样本标准差)(i P x s 表示为:
∑∑==--=
k
j j
k
j i j j
i P n
x s n
X s 1
1
2
)
1()
()1()(。

式中j n 为第j 组的测量次数。

◆ 以上计算得到的合并样本标准差仍是单次测量结果的实验标准差。

◆ 若实际工作中最后给出的测量结果是由h 次测量结果的算术平均值,则该平均值的实验标准差为:h
x s x s i p )()(=。

(3)极差法
在重复性条件下,对被测量X 做n 次独立测量,n 个测量结果中最大值与最小值之差R 称为极差,在可以估计被测量X 接近正态分布的前提下,单次测量结果i x 的标准不确定度(实验标准差)可表示为:
C
R
x s x u i i =
=)()( 式中级差系数C 如下表,其值与测量次数有关:
◆ 一般在测量次数较少时采用该法。

n
2 3 4 5 6 7 8 9 10 15 20 C
1.13
1.69
2.06
2.33
2.53
2.70
2.85
2.97
3.08
3.47
3.73
(4)最小二乘法
当被测量X 的估计值是由实验数据通过最小二乘法拟合的直线或曲线得到时,则任意预期的估计最,或拟合曲线参数的标准不确定度均可以利用已知的统计程序计算得到。

一般来说,两个物理量X 和Y 之间的关系问题,且估计值y x 和之间有线性关系bx a y +=。


y x 和独立测得n 组数据,其结果为),(,),,(),,(2211n n y x y x y x ⋯,且2>n 。

同时假定x 的测量不确定度远小于y 的测量不确定度(即x 的测量不确定度)(x u 可以忽略不计),则可利用最小二乘法得到参数b a ,(拟合直线方程的截距和斜率)以及它们的标准不确定度)()(b u a u 和。

由于测得的i y 存在误差,因而通常i i bx a y +≠,于是bx a y +=的误差方程可以写为:
)(111bx a y v +-= )(222bx a y v +-=
……
)(n n n bx a y v +-=
将上列各等式两边平方后相加,可得残差i v 的平方和为:∑∑==+-=n
i i i n
i i bx a y v 1
21
2)]([
为使残差i v 的平方和∑=n
i i v 1
2达到最小值,必须使上式对b a 和的偏导数同时为零。

于是由
01
2
1
2
=∂∂=∂∂∑∑==b
v a
v n
i i
n
i i

可得
[]
222)(2)(1
1
2
=-+=---=∂+-∂∑∑==y n x nb na bx a y a
bx a y n
i i i n
i i i 和
0222])[(2)]([1
1
2
1
1
2
=-+=---=∂+-∂∑∑∑∑====n
i i i n i i n i i i i n
i i i y x x b x na x bx a y b
bx a y
得到联立方程:{
1
1
2=-+=-+∑∑==n
i i i n i i y x x b x na y n x nb na
对a 求解得:x b y a -=;
对b 求解得:∑∑∑∑=∙∙
===∙∙--=
=-+-n
i i
n
i i
i n i n
i i i i x
x n x
y
x
n y x b y x x b x x nb y n 1
2
11
1
20,于是
假设∑∑∑=∙=∙--=+-=--=n
i i i n i i i i n i i i xy y x n y x y x y x y x y y x x S 1
1
1
)())((
∑∑∑==∙=-=+-=-=n
i i n i i i n i i xx x n x x x x x x x x S 1
22
1
2
1
2
)()2()(。

最后得到xx
xy S S b =
将b a ,的值代回误差方程,可求得残差i v 和残差的平方和∑=n
i i v 1
2。

于是y 的实验标准差)(y s 为:
2
-)(1
2
-=

=n v v y s n
i i
)(。

通过计算b a 和的方差,可以得到它们的标准不确定度为:
2
1
()()n
i
i xx
x
u a s a s nS =∙
==∑
而参数b a 和是由同一组测量结果计算得到的,因此两者之间理应存在一定的相关性,由于
x b a y +=,对等式两边求方差后得到:),(2)()()(2
x b a x b V a V x b a V n
s σ++=+=
)()(),(2)()()(222b xs a s b a r b s x a s ∙∙++=
n
x S s x b a r S s x nS x s
n
i i
xx
xx xx
n
i i
∑∑=∙
∙∙∙=∙++=1
2
2
22
12
2
),(2)(
于是b a 和之间的相关系数),(b a r 为:∑∑∑∑∑==∙==∙=-=
--=
--=
n
i i
n
i i
n
i i xx n
i i
xx
xx xx n
i i
x n x n x n x x n x S n
x
S x S x nS x n b a r 12
1
2
2
1
21
2
2
12
2)(2)(1),(
在Y 轴上拟合值0y 的标准不确定度
当对x 进行测量,测得值为0x ,并通过参数b a 和得到拟合值0y 时,可以计算出0y 的标准不确
定度)(0y u 。

测得值0x 与拟合值0y 之间满足关系:00bx a y +=。

其方差为:)()(),(2)()()(02
00b s a s b a x b V x a V y V ++==
由于xx xx
xx
n
i i
n
i i
S x x s S s
nS x
s
x n x n x b s a s b a r x 021
212
0022)()(),(2-=-=∙
=∙
∙=∙
∑∑
于是:xx
xx xx xx xx xx xx
n
i i
S x
x S x nS x n S s S x x s S x s nS x
s
y u 02
02022021
2
2
02)(2)(-++=-+=
∙=∙
∑ 将上式简化后得到:
在X 轴上拟合值0x 的标准不确定度
当对y 重复测量p 次,得到y 的平均值y ,并通过参数b a 和得到拟合值0x 时,同样可以求出0
x 的标准不确定度020()
()11x xx
x x S u b
p n S -=++。

4、不确定度的B 类评定
获得B 类评定标准不确定度的信息来源:
● 以前的观测数据;
● 对有关技术资料和测量仪器特性了解和经验; ● 生产部门提供的技术说明文件;
● 校准证书、检定证书或其他文件提供的数据、准确度的等级或级别、误差限等; ● 手册或某些资料给出的参考数据及不确定度;
● 规定实验方法的国家标准或类似文件中给出的重复性限或复现性限。

(1)信息来自校准证书或检定证书
自校准证书或检定证书给出的误差为扩展不确定度,根据扩展不确定度和标准不确定度之间的关系,可求出标准不确定度:()
()U x u x k
=
(2)信息来自测量仪器的误差 标准不确定度为:()3
A
u x =
,式中A 为仪器的误差。

(3)信息来自测量仪器的分辨力
标准不确定度为:()0.29x u δ=∙,式中δ为仪器的分辨力。

(4)信息来自数据修约
标准不确定度为:()0.29u x δ=,式中δ为数字修约。

(5)信息来自方法中的重复性限 标准不确定度为:() 2.83
r
u x =
,式中r 为重复性限。

(6)信息来自方法中的复现性限 标准不确定度为:() 2.83
R
u x =,式中R 为复现性限。

5、合成标准不确定度
(1)灵敏系数i c 和不确定度分量
根据各输入量的标准不确定度)(i x u ,以及由数学模型或实际测量得到的灵敏系数i c ,就可以得到对应于各输入量的标准不确定度分量)(y u i 。

)()(i i i x u c y u =。

灵敏系数i c 可由数学模型对输入量i x 求偏导数得到:i
i x y c ∂∂=。

当无法得到灵敏系数的可靠数学表达式时,灵敏系数也可以有实验测量得到。

在数值上它等于输入量i x 变化一个单位时,被测量y 的变化量,即后者与前者的比值。

(2)输出量等于各输入量加和的数学模型的合成标准不确定度
输出量合成标准不确定度)(y u c 可表示为各输入量标准不确定度分量)(i y u 的合成方差的正平方根:
2
1
()()n
c i
i u y u
y ==
∑。

(3)输出量等于各输入量相乘的数学模型的合成标准不确定度
输出量合成相对标准不确定度)(y u crel 可表示为各输入量相对标准不确定度)(i rel x u 的合成方差的正平方根:∑==
n
i i rel crel x u y u 1
2
)()(。

由于()y
y u y u c crel =)(,i i i rel x x u x u )()(=,则∑=∙=n
i i
i c x x u y y u 1
2
))((
)(。

(4)输出量与各输入量成幂函数的数学模型的合成标准不确定度 若n p
n p
p
x x bx y ⋯=2121,式中b 为比例常数,如22==b z x y p 中,。

则∑==
n
i i rel
crel x u
p y u 1
22)()(,导出∑=∙
=n
i i c x x u p y y u 1
22))(()(。

可以看出,若指数1=p ,第(4)种情况即为第(3)种情况。

(5)输出量与各输入量既有加成关系,又有相乘的关系时的数学模型的合成标准不确定度 出现该种情况,先计算相乘关系的不确定度分量(即用相对标准不确定度计算),再计算加成关系的标准不确定度分量。

(6)合成标准不确定度中相关性的处理
当各输入量之间存在不可忽略的相关性时,合成标准不确定度为:
∑∑∑∑∑-=+=∙∙=∙==∙+=∂∂∂∂=111
12
112
),(2)(),()(n i n
i j j i j i n j i i j i n
i n
j j i c x x u c c x u c x x u x f x f y u
式中)(j i x x u 为输入量j i x x 和之间的协方差。

由于相关系数定义为:)
()(),(),(j i j i j i x u x u x x u x x r =
,也可以用相关系数来表达成为:
),()()(2)()(111
1
2
2
j i n i n
i j j i j i n j i i c x x r x u x u c c x u c y u ∙-=+=∙∙∙=∑∑∑+=。

用不确定度分量表示为:
),()()(2)()(111
1
2j i n i n
i j j i n
j i
c x x r y u y u y u
y u ∙-=+=∙=∑∑∑+=
若考虑仅有两个输入量的情况21x x y +=:
● 若21x x 和之间不相关,即相关系数02,1=r ,此时合成标准不确定度等于两个不确定度分量之方和根,即2
221u u u c +=。

● 若21x x 和之间完全正相关,即相关系数12,1=r ,此时合成不确定度等于两个不确定度分量之和,即21u u u c +=。

● 若21x x 和之间完全负相关,即相关系数12,1-=r ,此时合成不确定度等于两个不确定度分量之差的绝对值,即21u u u c -=。

● 对于一般情况,21x x 和之间部分相关,即—1<2,1r <1,此时合成不确定度表示为:
2,1212
22
12r u u u u u c ++=。

若考虑仅有三个输入量的情况321x x x y ++=:
● 若321x x x 、、之间不相关,此时合成标准不确定度表示为:
222123c u u u u =++。

● 若21x x 和之间存在相关性,此时合成不确定度表示为:
2,1212
32
22
12r u u u u u u c +++=。

● 若三个输入量321,x x x 和之间均存在相关性,此时合成不确定度表示为:
3,1313,2322,1212
32
22
1222r u u r u u r u u u u u u c +++++=。

从原则上说,必须要知道相关系数后,才能求出合成标准不确定度。

相关系数:
● 通过实验,同时测量n 组输入量21x x 和之值,由公式得到输入量之间的相关系数和协方差:
1
)()(),(,)
()()
,(),(221
1121212121---=
=
∑=n x x x x
x x s x s x s x x s x x r i n
i i
● 输入量21x x 和的n 组测量结果的平均值21x x 和之间的相关系数和协方差为:
)
1()()(),(),,()()()
,(),(221
112121212121---===∑=n n x x x x x x s x x r x s x s x x s x x r i n
i i ● 由于相关系数的实验测量比较麻烦,因此在进行测量不确定度评定中除非确有必要,一般应尽量避免处理相关性。

相关行动处理有以下几种方法:
1) 采用合适的测量方法和测量程序,因可能避免输入量之间的相关性。

2) 如果可以选择测量不确定度评定中所采用的输入量,应尽量选用不相关的输入量。

3) 如果已知两个输入量之间存在相关性,若相关性较弱,则可以忽略其相关性。

4) 如果已知两个输入量之间存在相关性,若其本身在合成标准不确定度中不起主要作用,在可以忽略其相关性。

5) 如果已知两个输入量之间存在相关性,若相关性较强,则假定其相关系数为1。

6) 如果已知两个输入量之间存在相关性,若相关系数为负值,则可以忽略其相关性,只要最后得到的扩展不确定度满足要求。

7) 仅在以上方法全部都不适用的情况下,才考虑由实验测量并计算相关系数。

6、扩展不确定度
扩展不确定度U 等于合成标准不确定度c u 乘以包含因子k 。

因此必须先确定被测量y 可能得分布,以确定包含因子。

应建立不确定度分量一览表。

(1)被测量y 的分布接近于正态分布的判定(中心极限定理)及扩展不确定度计算 中心极限定理:如果一个随机变量是大量相互独立的随机变量之和,则不论这些随机变量具有何种类型的分布,该随机变量的分布近似于正态分布。

随着独立随机变量个数的增加,它们的和就越接近正态分布。

对被测量y 判断是否为正态分布的依据:
● 被测量y 用扩展不确定度p U 给出,而对其分布又没有特殊指明时,估计值y 的分布。

● 被测量y 的合成标准不确定度)(y u c 中相互独立的分量)(y u i 中,存在两个界限值接近的三角分布,或4个界限值接近的均匀分布时。

● 被测量y 的合成标准不确定度)(y u c 的相互独立的分量中,量值较大的分量(起决定作用的分量)接近正态分布时。

● 如果∑==n
i i i x c y 1,即被测量y 是各输入量i x 的线性函数,且各输入量均为正态分布并相互独
立,则被测量y 服从正态分布。

也就是说,正态分布的线性叠加仍是正态分布。

● 即使输入量i x 不是正态分布,根据中心极限定理,只要被测量y 的方差)(2y σ比各输入量i
x 的分量的方差)(22i x c σ大得多,或各分量的方差)(22i x c σ相互接近,则被测量y 近似的满足正态分布。

● 若取多次测量的被测量y 算数平均值y 作为最佳估计值(或结果),此时不论被测量y 为何种分布,随着测量次数的增大,y 的分布趋于正态分布。

● 对于正态分布,包含因子k 与置信概率p 的关系如下表: P% 50 68.27 90 95 95.45 99 99.73 k p
0.67
1
1.645
1.960
2
2.576
3
● 被测量接近正态分布时,原则上应计算各分量的自由度和合成标准不确定度的有效自由度,根据置信概率p 由t 分布表得到包含因子
p
k ,此时扩展不确定度为:c p c p p u v t u k U ∙∙==)(。

● 若有些测量程序确保了被测量y 自由度不会太小(15以上),也可以不计算自由度,仍可以估计其置信概率大体上分别为)3%(99)2%(95==k k 和,此时扩展不确定度为:
c c c c u u k U u u k U 32%99%95====∙∙或
● 若有些测量程序使被测量y 自由度较小,选用32或=k 时对应的置信概率可能与95%或99%相差甚远。

ISO/TS14253中提供了一种补偿办法:在计算扩展不确定度时,增加一个安全因子h 作为乘数,即c c c c hu u h k U hu u h k U 32%99%95====∙∙∙∙或。

2=k 时对应的安全因子如下:
测量次数
2
3
4
5
6
7
8
9
10
n h (2=k )
7.0
2.3
1.7
1.4
1.3
1.3
1.2
1.2
1
(2)被测量y 的分布为某种非正态分布的判定及扩展不确定度计算 ● 被测量y 接近于矩形(均匀)分布的判定 1) 数据修约导致的不确定度。

2) 数字式测量仪器对示值量化(分辨力)导致的不确定度。

3) 测量仪器由于滞后、摩擦效应导致的不确定度。

4) 按级使用的数字式仪表、测量仪器最大允许误差导致的不确定度。

5) 用上、下界给出的线膨胀系数。

6) 测量仪器度盘或齿轮回差引起的不确定度。

7) 平衡指示器调零不准导致的不确定度。

此时扩展不确定度为:c c c c u u k U u u k U 71.165.1%99%95====∙∙或 ● 被测量y 接近于三角分布的判定
1) 相同修约间隔给出的两独立量之和或差,由修约导致的不确定度。

2) 因分辨力引起的两次测量结果之和或差的不确定度。

3) 用替代法检定标准电子元件或测量衰减时,调零不准导致的不确定度。

4) 两相同均匀分布的合成。

此时扩展不确定度为:95%99%1.90 2.20c c c c U k u u U k u u ∙∙====或 (3)无法判定被测量y 分布的扩展不确定度计算
由于无法判断被测量y 的分布,也就无法根据规定的置信概率求出包含因子k ,此时假设k 值取2(认为置信概率95%)或3(认为置信概率99%)。

JJF1059—2012中规定:k 值一般在2~3之间,在大多数情况下2=k ,当取其他值时,应说明其来源。

7、自由度
(1)不确定度的A 类评定的自由度
● 用贝塞尔公式计算实验标准差时,若测量次数为n ,则自由度为1-n v =。

● 当同时测量t 个被测量(输入量)时,自由度为t n v -=。

● 若t 个被测量之间另有m 个约束条件,自由度为m t n v +=-。

● 对于合并样本标准差,其自由度为各组的自由度之和。

● 当用极差法估算实验标准差时,其自由度与测量次数n 的关系如下表:
(2)不确定度的B 类评定的自由度 不确定度的B 类评定的自由度:2
)()]([21⎭

⎫⎩⎨⎧=
x u x u u v ,式中
)(x u 为被测量x 的标准不确定度。

)]([x u u 为标准不确定度)(x u 的标准不确定度。

)
()]
([x u x u u 为标准不确定度)(x u 的相对标准不确定度。

相对标准不确定度
)
()]
([x u x u u 与自由度v 的关系如下: )()]
([x u x u u 0.10 0.20 0.25 0.30 0.40 0.50
v
50 12.5 8 5.5 3.1 2
● 自由度越大,说明不确定度越可靠。

主要依靠不确定度信息来源确定自由度大小。

● 当相对不确定度的评定有严格的数字关系(如设备的误差、数字修约引起的不确定度)时,自由度较大。

● 当计算相对不确定度的数据来自校准证书、检定证书、手册等时,自由度一般在20—50之间。

● 当相对不确定度的计算带有一定的主观因素时,自由度较小。

n 2 3 4 5 6 7 8 9 10 15 20 v
0.9
1.8
2.7
3.6
4.5
5.3
6.0
6.8
7.5
10.5
13.1
(3)合成标准不确定度的有效自由度
合成标准不确定度的自由度称为有效自由度,以eff v 表示。

当合成标准不确定度有两个或两个以上不确定度分量合成,并且其分布接近于正态分布时,合成标准不确定度的自由度计算公式为:
∑==n i i i c eff v y u y u v 1
4
4
)()
(。

当用相对标准不确定度来评定时的公式为:∑∑∑======n
i i i i c n i i i rel i el r c n i i rel i el r c eff v x x u r y u v x u r y u v y u y u v 14
41
444
144])/([)/y]([)()()()( 四、测量结果的处理 1、被测量估计值的计算表达
(1)∑∑==⋯===n
i ki i i n
i i
x x x f n n y
y y 1
211
),,(1,
根据对所有输入量k x 的第i 次观察结果,计算出第i 次的测量结果i y ,然后再对y 取平均值。

(2))
,,,(1211
k n
i ki
k x x x f y x n x ⋯==∑= 先将每个输入量k x 的n 次独立观测值取平均值,然后再由各输入量的平均值得到被测量。

2测量结果中离群值的剔除
(1)σ3准则
对被测量做n 次独立测量,计算实验标准差。

在正态分布情况下,当某一残差x x v i i -=的绝对值超过三倍实验标准差,即)(3i i x s v >,则认为该测量结果属于离群值而应予剔除。

将离群值剔除后,重新计算实验标准差,对剩余数据进行判别、处理,直到测量结果中不包含离群值为止。

(2)格拉布斯准则
对被测量做n 次独立测量,计算实验标准差。

在正态分布情况下,当某一残差x x v i i -=超过实验标准差与临界系数)(n g 的乘积,即g(n))(∙>i i x s v ,则认为该测量结果属于离群值而应予剔除。

将离群值剔除后,重新计算实验标准差,对剩余数据进行判别、处理,直到测量结果中不包含离群值为止。

临界系数)(n g 表 n g(n) n g(n) n g(n) 3 1.155 10 2.290 17 2.260 4 1.481 11 2.355 18 2.651 5 1.715 12 2.412 19 2.681 6 1.877 13 2.462 20 2.709 7 2.020 14 2.507 30 2.908 8 2.126 15 2.549 40 3.306 9 2.215
16
2.585
50
3.128
3、数据修约
(1)有效数字
当一个近似数所引入的误差的绝对值小于该近似数末位数的0.5,从该近似数左边第一个非零数字算起,直到最后末位数为止是有效数字。

例如:
1) 对3.14159265截取到百分位,为3.14,引入的误差绝对值为:
2
01
.00015926.03.14159265-3.14<
=,所以近似数3.14为原数值的3位有效数字。

2) 对3.14159265截取到千分位,为3.141,引入的误差绝对值为:
201
.005926000.03.14159265-3.141>
=,所以近似数3.141末位数不是原数值的有效数字。

必须
将其进位成3.142,此时所引人的误差绝对值为:201
.0004073500.03.14159265-3.142<
=,所以近似
数3.142为原数值的4位有效数字。

JJF1059—1999规定,合成标准不确定度、扩展不确定度以及输入量的估计值的标准不确定度通常为两位。

在实际计算过程中,为了避免过大的数据修约误差,可以多保留数值的位数。

(2)修约间隔
修约间隔是确定保留位数的一种方式,也称为修约区间。

修约间隔一经确定,修约数只能是修
约间隔的整数倍。

修约间隔一般以n k 10⨯的形式表示,称为以”“k 间隔修约,并由n 确定修约到哪
一位。

数据会引入不确定度,其大小与修约间隔有关,若修约间隔为x σ,则修约后可能引入的最大误差为2/x σ,由于数据修约引起的不确定度满足矩形分布,固由修约引入的标准不确定度为:
x x
u σσ29.03
2==
(3)修约规则
1) 对于“1”间隔修约,若舍去的数值小于所保留数值末位的0.5,则保留数值的末位数字不变。

2) 对于“1”间隔修约,若舍去的数值大于所保留数值末位的0.5,则保留数值的末位数字加1。

3) 对于“1”间隔修约,若舍去的数值等于所保留数值末位的0.5,则保留数值的末位数字按奇偶规则进行修约,即当末位为偶数时,末位数字不变;当末尾数字为奇数时,末位数字加1.
4) 对于非“1”间隔修约,例如“2”或“5”间隔修约,可先将拟修约数除以2或5,然后按“1”间隔修约,然后再将修约数乘以2或5.
5) 负数的修约按其绝对值进行,修约后再加上负号。

6) 数据修约应一步到位,不得连续修约。

连续修约会导致修约不确定的增大。

7) 在计算合成标准不确定度的有效自由度并取整时,习惯上采取只舍不进的规则(较小的自由度对应较大的包含因子)。

五、测量不确定度的报告和表示
1、测量不确定度报告
比较重要的测量,不确定度报告一般包括以下内容: ● 数学模型和对应于各输入量的灵敏系数i c 。

● 修正值和常熟的来源及不确定度。

● 输入量i x 的实验观测数据及其估计值i x ,标准不确定度)(i x u 的评定方法及其量值和自由度,并列表表示。

● 对所有有相关性的输入量给出协方差或相关系数及其获得方法。

● 测量结果的数据处理程序,该程序应易于重复。

2、合成标准不确定度)(y u c 的报告形式
在报告以下测量结果时,使用合成标准不确定度)(y u c ,同时给出有效自由度eff v : ● 基础计量学研究。

● 基本物理常数测量。

● 复现国际单位制单位的国际对比。

3、扩展不确定度的报告形式
(1)当用U 报告扩展不确定度时
当扩展不确定度用)(y ku U c =表示时,应给出k 值,测量结果表示为: m=100.111g ;U=11mg ;k=2;或m=(100.111±0.011)g ;k=2。

(2)当用p U 报告扩展不确定度时
● 当扩展不确定度用)(y ku U c p =表示时,应给出p 值和有效自由度eff
v 。

测量结果表为:
m=100.111g ;U 95=11mg ;eff v =8。

m=(100.111±0.011)g ;eff v =8;括号内第二项为U 95之值。

m=100.111(11)g;
v=8;括号内第二项为U95之值,其末位前面结果末尾数对齐。

eff
m=100.111(0.011)g;
v=8;括号内第二项为U95之值,与前面结果有相同的计量单位。

eff
●当被测量接近于非正态分布时,应给出p值和k值,同时指出其分布类型。

●。

相关文档
最新文档