相关系数的计算相关系数一般可以通过实验测量的方法或理论经验的分析得到 即一是用同时观测两个量的方法确定相关系数估计值；二是当两个量 或以上均因与同一个量有关而相关时，依据相关系数定义公式，计算 相关系数的估计值1.根据对X和y两个量同时测量的n组测量数据,相关系数的估计值按公式(4)计算：h(x - X) (y - Y)iir ( x, y) = 4^(n - 1)s( x) s( y)4)式中, s(x), s(y)---为 X 和 Y 的实验标准偏差公式(4)还可以表示为：()丫()uk, xAr,y =十十㈠=-u\x )• uh 丿,j起k=1 vX —X ik ik=1 八X —yik ik =1 ' 丿k=1X —Xjk j丿X —Xjk j丿5)示例1:用同一钢卷尺测量某矩形的面积，对矩形的长(1)和宽d ) 各测量10次,其测量列如表1 所示表1 矩形长和宽的测量数据1 (mm )40.140.240.040.140.140.040.140.140.240.11=40.10d (mm 丿20.020.220.020.120.120.020.020.120.120.1d =20.07矩形面积的数学模型：s = d X1，因为对长和宽采用了同一测量仪器，则它们的估计值会出现相关，根据表1有l和d算术平均值的标准不确定度为：,、艺（ _ l） = 0.021（mm） 10 x（10 — 1丿 i丫 i=1冷）= 丄 2 C — d） =0.021（mm）10x（10-1丿］i兰（一 /'d 一 d ）=0.03 Cmm2）iii=1芳（-1L 0.04 Cmm 2） ii 二1d — d ） = 0.041 Cm 2）ii=1所以相关系数艺（—l）（— d ）r i 0.03 0 74„■£ d — l）昱 C - d） ®040X 0.041i ir i=1 i=1面积 S=/ • d =804.81mm2则考虑相关系数r得：uc（s）』曲2 I-站2 + 2止塾2 = 015%s [L l J L d J l d J当不考虑相关系数r时，uS =J-曲2 +L站2F = 012%s [L 门 L d」J ■从以上两式的结果可以看出考虑相关系数与不考虑相关系数存在明显的区别，不考虑相关系数时，明显使评定的不确定度偏小。
2.当两个量均因与同一个量有关而相关时，计算相关系数的估计 值假如在得到两个输入量的估计值 x 和 x 时，是使用了同一个测 ij 量标准、测量仪器或参考数据或采用了相同的具有相当大不确定度的 测量方法，则 x 和 x 两个量均因与同一个量有关而相关ij示例 2：2014 年度一级注册计量师考试《测量数据处理及计量专业实务》科目中的单项选择题第26题为“用lk0的标准电阻R校准s标称值均为lk0的两个电阻器，校准值R R , R二a R .已知标准1 1 S 2 2 S电阻 R 的标准不确定度为u(R L 0.10 ， 若ssa〜1, a〜1, uCx )= uCx )= 1 x 10一4，假设a、a、R互不相关，则R与1 2 1 2 1 2 1 1R 的相关系数为( )2A.1.0 B.0.75 C.0.5 D.0.25 ”[解]1)每个电阻R校准时与标准电阻R比较得到比值a ,校准值i s i为：R=a Ri i s2) 根据不确定度传播定律，每个 R 的标准不确定度：iu(R) = A；i RSu(ai )]2 + [aiu( RS )]2式中的u(a )对每一个校准值近似相等，且a心1,由比较仪的ii不确定度为 u(a )= u(a)= u(a )= 1 x 10-4,则:i 1 2u(Ri) =iR 2u 2 (a ) + u 2(R )S i S3) 任意两个电阻校准值的相关系数：Ri =a iRS ； Rj j jRSRi、 R 之间协方差的估计值： ijdR ORu(R，R)= i Lu2(RS) = aa u2(R ) = a2u2(R )' i / OR OR S 1 j S SSS由于a ~a二a 心1,协方差 u(R , R)二 U2(R)i j i j SR 、 R 之间相关系数：ij= u(R ,R ) u2(R )r..= r (R , R ) = i j = —ij i j u(R )u(R ) R 2u 2 Q 丿+ u 2 \R )i j S i Sr u (x)]2< 1 +u(R ) ■ R1- S • S」>由题意知，u(R》R二10一4 ； u(d )= uCx )= 1 x 10一4代入上式，得S ' S 1 2r =r(R ， R )=0.5ij i j本题正确选项为： C.0.5分析可知，u(R ) R = 10-4; u(a )= u(a )= 100x 10-6; r =0.5S ' S 1 2 ij如果 u(x )= u(x )= 10x 10-6; r 心0.9901 2 ij如果 u(x )= u(x )= 1 x 10-6; r 心1.0001 2 ij因此当 u(x )t 0, r T1 和 u(R )t u(R )iS一般来说，在与校准值比较时，如本示例，已校项的估计值间是 相关的，其相关的程度取决于校准过程(比对过程)引入的不确定度 与参考标准的不确定度之比。
仅当与参考标准的不确定度相比，校准 过程(比对过程)的不确定度可以忽略不计时，相关系数等于+1，并 且每个校准项的不确定度与其参考标准的不确定度相同示例3：在示例2的条件中，若将R和R串联成R =R +R =2k01 2 ref 1 2的电阻，试确定串联后的R =2k0电阻的合成标准不确定度u (R )ref c ref 为多少？[解法一]1 )每个电阻R校准时与标准电阻R比较得到比值x ,校准i s i 值为：R=x Ri i s2) 根据不确定度传播定律，每个 R 的标准不确定度：iu(R) = vi RSu(x. )]2 + [x.u( RS )]2式中的u(x)对每一个校准值近似相等，且x心1,由比较仪的ii不确定度为 u© )二 uCx )= uCx )= 1 x 10-4,则:iu(Ri) =iR 2u2(x ) + u2(R )S i S=\；( X 103 X 1 X 10 - 4 ) + 0.12 0=.込 X0.10即：u(R)=u(R)=迈 X0.101 2 琴3) 根据示例2的计算结果，R、R两个电阻校准值的相关系数:12r(R ， R)=0.5124) 串联电阻 R 的合成标准不确定度：ref根据 R 的测量模型： R = R+Rref ref 1 2R 的合成标准不确定度为：refU (R )= 'u 2(R ) + u 2(R ) + 2u (R )u (R ) - r(R , R )c ref 1 2 1 2 1 2=話(迈 X 0.1)2 + (迈 X 0.1)2 + 2 X (迈 X 0.1) X (迈 X 0.1) X 0.5 0「6x0.12 05.250［解法二］ 由于输出量 R = R+Rref 1 2输入量 R 的不确定度由两个分量构成，其一来自标准电阻iu(R)=0.10，作为U(R )的分量，其灵敏系数C均为+1；其二来自校 s c ref i准(比对)过程，已知为Ru(a )=1X103X1X10-40=0.10,灵敏系数si也C均为+1。
i可以认为R 一共有4个标准不确定度分量，其中2个为标准器refR引入的分量均为0.10,它们之间为强相关，可设定r=+1，这2个 s0.10取代数和合成为： 0.10+0.10=0.20；另2个为校准过程(比对过程)引入的分量Ru(a )=1X103X1X10-40=0.10，主要是随机效si应引起，可以设定彼此独立，r=0，而按方和根合成为：迈X0.10；来源于R与来源于校准过程彼此不相关，因而su (R )二\：'0.22 + (*：2 x 0.1)2 0二斗 6 x 0.12 0心0.250c ref "解法二的特点是把相关和不相关的不确定度分量分别合成后再合成这样的分组合成方法，在JJF1059.1-2012的引言中就已明确: “测量不确定度能从对测量结果有影响的不确定度分量导出，且与这 些分量怎样分组无关，也与这些分量如何进一步分解为下一级分量无 关”。