《相关系数的计算》由会员分享,可在线阅读,更多相关《相关系数的计算(6页珍藏版)》请在金锄头文库上搜索。
1、相关系数的计算相关系数的计算相关系数一般可以通过实验测量的方法或理论经验的分析得到。即一是用同时观测两个量的方法确定相关系数估计值;二是当两个量或以上均因与同一个量有关而相关时,依据相关系数定义公式,计算相关系数的估计值。1根据对x和y两个量同时测量的n组测量数据,相关系数的估计值按公式(4)计算:(4)1() () ( , )(1) ( ) ( )nii ixXyY r x yns x s y 式中,s(x),s(y)-为X和Y的实验标准偏差。公式(4)还可以表示为:(5) nknkjjkiiknkjjkiikjijixxxxxxxxxuxuxxuyxr11221,示例示例1:用同一钢卷尺测
2、量某矩形的面积,对矩形的长( )和l宽(d ) 各测量10次,其测量列如表1 所示。表1 矩形长和宽的测量数据mml40.140.240.040.140.140.040.140.140.240.1 =40.10lmmd20.020.220.020.120.120.020.020.120.120.1=20.07d矩形面积的数学模型:,因为对长和宽采用了同一测量lds仪器,则它们的估计值会出现相关,根据表1有 和d 算术平均值l的标准不确定度为:mmlllsii021. 01101011012 mmdddsii021. 01101011012 =0.03 ddlli ii 1012mm04. 02
3、101 iill2mm041. 02101 iidd2mm所以相关系数 74. 0 041. 0040. 003. 021012101101 ii iii iiddllddll r面积 S= =804.81mm2ld则考虑相关系数 r 得: %15. 0221 22 dds llsrdds lls ssuc当不考虑相关系数r时, %12. 021 22 dds lls ssuc从以上两式的结果可以看出考虑相关系数与不考虑相关系数存在明显的区别,不考虑相关系数时,明显使评定的不确定度偏小。2.当两个量均因与同一个量有关而相关时,计算相关系数的估计值。假如在得到两个输入量的估计值xi和xj时,是使
4、用了同一个测量标准、测量仪器或参考数据或采用了相同的具有相当大不确定度的测量方法,则xi和xj两个量均因与同一个量有关而相关。示例示例 2 2:2014 年度一级注册计量师考试测量数据处理及计量专业实务科目中的单项选择题第 26 题为“用 1k的标准电阻Rs校准标称值均为 1k的两个电阻器,校准值,.已知SRR11SRR22标准电阻Rs的标准不确定度为,若 1 . 0sRu, 11, 12,假设、R1 互不相关,则R1与R2的相关 4 21101uu12系数为( ) 。A.1.0 B.0.75 C.0.5 D.0.25 ”解 1)每个电阻Ri校准时与标准电阻Rs比较得到比值i,校准值为:Ri
5、=iRs2) 根据不确定度传播定律,每个Ri的标准不确定度:u(Ri) =22)()(SiiSRuuR式中的u(i)对每一个校准值近似相等,且i1,由比较仪的不确定度为u(i)= , 则: 4 21101uuu(Ri) =)()(222 SiSRuuR3) 任意两个电阻校准值的相关系数:;SiiRRSjjRRRi、Rj之间协方差的估计值:u(Ri,Rj)=)()()(2222 SSjiS SjSiRuRuRuRRRR由于i j= 1,协方差u(Ri,Rj) = u2(RS)Ri、Rj之间相关系数:=ijr SiSSjiji jiRuuRRu RuRuRRuRRr2222)( )()(),(),
6、(= 12 1 SSRRuu由题意知,; 代入上式,得410SSRRu 4 21101uu=r(Ri,Rj)=0.5ijr本题正确选项为: C.0.5分析可知,; ; =0.5410SSRRu 6 2110100uuijr如果 ; 0.990 6 211010uuijr如果 ; 1.000 6 21101uuijr因此当 和 , 0u1r SiRuRu一般来说,在与校准值比较时,如本示例,已校项的估计值间是相关的,其相关的程度取决于校准过程(比对过程)引入的不确定度与参考标准的不确定度之比。仅当与参考标准的不确定度相比,校准过程(比对过程)的不确定度可以忽略不计时,相关系数等于+1,并且每个校
7、准项的不确定度与其参考标准的不确定度相同。示例示例 3 3:在示例 2 的条件中,若将R1和R2串联成Rref =R1 +R2 =2k的电阻,试确定串联后的Rref =2k电阻的合成标准不确定度uc(Rref)为多少?解法一1)每个电阻Ri校准时与标准电阻Rs比较得到比值i,校准值为:Ri =iRs2) 根据不确定度传播定律,每个Ri的标准不确定度:u(Ri) =22)()(SiiSRuuR式中的u(i)对每一个校准值近似相等,且i1,由比较仪的不确定度为u(i)= , 则: 4 21101uuu(Ri) =)()(222 SiSRuuR=22431 . 0101101=0.12即:u(R1)
8、=u(R2)= 0.123)根据示例 2 的计算结果,R1 、R2两个电阻校准值的相关系数:r(R1,R2)=0.54)串联电阻Rref的合成标准不确定度:根据Rref的测量模型: Rref= R1 +R2Rref的合成标准不确定度为: uc(Rref)=),()()(2)()(212122 12RRrRuRuRuRu=5 . 0) 1 . 02() 1 . 02(2) 1 . 02() 1 . 02(22=0.2521 . 06解法二由于输出量 Rref= R1 +R2输入量Ri的不确定度由两个分量构成,其一来自标准电阻u(Rs)=0.1,作为uc(Rref)的分量,其灵敏系数ci均为+1;
9、其二来自校准(比对)过程,已知为Rsu()=1103110-4=0.1,灵敏i系数也ci均为+1。可以认为Rref一共有 4 个标准不确定度分量,其中 2 个为标准器Rs引入的分量均为 0.1,它们之间为强相关,可设定r=+1,这 2个 0.1取代数和合成为:0.1+0.1=0.2;另 2 个为校准过程(比对过程)引入的分量Rsu()=1103110-4=0.1,主要是i随机效应引起,可以设定彼此独立,r=0,而按方和根合成为:0.1;来源于Rs与来源于校准过程彼此不相关,因而2uc(Rref)= =0.2522) 1 . 02(2 . 021 . 06解法二的特点是把相关和不相关的不确定度分量分别合成后再合成。这样的分组合成方法,在 JJF1059.1-2012 的引言中就已明确:“测量不确定度能从对测量结果有影响的不确定度分量导出,且与这些分量怎样分组无关,也与这些分量如何进一步分解为下一级分量无关” 。