合成标准不确定度评定中相关系数的实例分析
JJF1059.1-2012《测量不确定度评定与表示》与JJF 1059-1999版相比,主要修订内容之一就是“合成标准不确定度评定中增加了各输入量间相关时协方差和相关系数的估计方法,以便规范处理相关的问题”。
当输入量Xi明显相关时,其合成标准不确定度uc(y)必须考虑协方差项,且需预先求得相关系数后,才能计算合成标准不确定度。本文针对合成标准不确定度评定中相关系数的计算进行实例分析,对相关系数的计算有一定的适用性。
一、不确定度传播律和相关系数性质
当被测量Y由N个其它量X1,X2,…,XN通过测量函数f确定时,被测量的估计值y为:
被测量的估计值y的合成标准不确定度uc(y)按下式计算:
(1)
上式称为不确定度传播律,式中: r(xi, xj) u(xi) u(xj)= u(xi, xj) 是输入量xi与xj的协方差,r(xi, xj)为输入量xi与xj之间的相关系数。
相关系数是两个变量之间相互依赖性的度量,它等于两个变量间的协方差除以各自方差之积的正平方根,用符号(x, y)表示
(2)
定义的相关系数是在无限多次测量条件下的理想概念。有限次测量时相关系数的估计值用r(x,y)表示,
(3)
相关系数性质:
1)相关系数r(x ,y)= r(y ,x)∈[-1,1];
2)当r >0时,表示两变量正相关,当r<0时,表示两变量为负相关;当|r|=1时,表示两变量为完全线性相关,当r=1时,称为完全正相关(正强相关),而当r =-1时,称为完全负相关(负强相关);当r =0时,表示两变量间无线性相关关系。
二、相关系数的计算
相关系数一般可以通过实验测量的方法或理论经验的分析得到。即一是用同时观测两个量的方法确定相关系数估计值;二是当两个量或以上均因与同一个量有关而相关时,依据相关系数定义公式,计算相关系数的估计值。
1.根据对x和y两个量同时测量的n组测量数据,相关系数的估计值按公式(4)计算:
(4)
式中,s(x),s(y)---为X和Y的实验标准偏差。
公式(4)还可以表示为:
(5)
示例1:用同一钢卷尺测量某矩形的面积,对矩形的长( )和宽(d ) 各测量10次,其测量列如表1 所示。
表1 矩形长和宽的测量数据
40.1 40.2 40.0 40.1 40.1 40.0 40.1 40.1 40.2 40.1 =40.10
20.0 20.2 20.0 20.1 20.1 20.0 20.0 20.1 20.1 20.1 =20.07
矩形面积的数学模型: ,因为对长和宽采用了同一测量仪器,则它们的估计值会出现相关,根据表1有 和d 算术平均值的标准不确定度为:
=0.03
所以相关系数
面积S= • =804.81mm2
则考虑相关系数 r 得:
当不考虑相关系数r时,
从以上两式的结果可以看出考虑相关系数与不考虑相关系数存在明显的区别,不考虑相关系数时,明显使评定的不确定度偏小。
2.当两个量均因与同一个量有关而相关时,计算相关系数的估计值。假如在得到两个输入量的估计值xi和xj时,是使用了同一个测量标准、测量仪器或参考数据或采用了相同的具有相当大不确定度的测量方法,则xi和xj两个量均因与同一个量有关而相关。
示例2:2014年度一级注册计量师考试《测量数据处理及计量专业实务》科目中的单项选择题第26题为“用1k的标准电阻Rs校准标称值均为1k的两个电阻器,校准值 , .已知标准电阻Rs的标准不确定度为 ,若 ,假设 、 、R1 互不相关,则R1与R2的相关系数为( )。
A.1.0 B.0.75 C.0.5 D.0.25 ”
[解] 1)每个电阻Ri校准时与标准电阻Rs比较得到比值i,校准值为:
Ri =iRs
2) 根据不确定度传播定律,每个Ri的标准不确定度:
u(Ri) =
式中的u(i)对每一个校准值近似相等,且i≈1,由比较仪的不确定度为u(i)= , 则:
u(Ri) =
3) 任意两个电阻校准值的相关系数:
;
Ri、Rj之间协方差的估计值:
u(Ri,Rj)=
由于i ≈j= ≈1,协方差u(Ri,Rj) = u2(RS)
Ri、Rj之间相关系数:
=
=
由题意知, ; 代入上式,得
=r(Ri,Rj)=0.5
本题正确选项为: C.0.5
分析可知, ; ; =0.5
如果 ; ≈0.990
如果 ; ≈1.000
因此当 和
一般来说,在与校准值比较时,如本示例,已校项的估计值间是相关的,其相关的程度取决于校准过程(比对过程)引入的不确定度与参考标准的不确定度之比。仅当与参考标准的不确定度相比,校准过程(比对过程)的不确定度可以忽略不计时,相关系数等于+1,并且每个校准项的不确定度与其参考标准的不确定度相同。
示例3:在示例2的条件中,若将R1和R2串联成Rref =R1 +R2 =2k的电阻,试确定串联后的Rref =2k电阻的合成标准不确定度uc(Rref)为多少?
[解法一]1)每个电阻Ri校准时与标准电阻Rs比较得到比值i,校准值为:
Ri =iRs
2) 根据不确定度传播定律,每个Ri的标准不确定度:
u(Ri) =
式中的u(i)对每一个校准值近似相等,且i≈1,由比较仪的不确定度为u(i)= , 则:
u(Ri) =
=
= ×0.1
即:u(R1)=u(R2)= ×0.1
3)根据示例2的计算结果,R1 、R2两个电阻校准值的相关系数:r(R1,R2)=0.5
4)串联电阻Rref的合成标准不确定度:
根据Rref的测量模型: Rref= R1 +R2
Rref的合成标准不确定度为:
uc(Rref)=
=
= ≈0.25
[解法二]由于输出量 Rref= R1 +R2
输入量Ri的不确定度由两个分量构成,其一来自标准电阻u(Rs)=0.1,作为uc(Rref)的分量,其灵敏系数ci均为+1;其二来自校准(比对)过程,已知为Rsu( )=1×103×1×10-4=0.1,灵敏系数也ci均为+1。
可以认为Rref一共有4个标准不确定度分量,其中2个为标准器Rs引入的分量均为0.1,它们之间为强相关,可设定r=+1,这2个0.1取代数和合成为:0.1+0.1=0.2;另2个为校准过程(比对过程)引入的分量Rsu( )=1×103×1×10-4=0.1,主要是随机效应引起,可以设定彼此独立,r=0,而按方和根合成为: ×0.1;来源于Rs与来源于校准过程彼此不相关,因而
uc(Rref)= = ≈0.25
解法二的特点是把相关和不相关的不确定度分量分别合成后再合成。这样的分组合成方法,在JJF1059.1-2012的引言中就已明确:“测量不确定度能从对测量结果有影响的不确定度分量导出,且与这些分量怎样分组无关,也与这些分量如何进一步分解为下一级分量无关”。
从上述各例分析可以看出,处理合成标准不确定度评定中的相关问题时,不应忽视协方差项的影响,对协方差项中的相关系数的性质及计算等问题的正确处理是十分重要的。
|