合成標準不確定度評定中相關系數的實例分析
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的引言中就已明確:“測量不確定度能從對測量結果有影響的不確定度分量導出,且與這些分量怎樣分組無關,也與這些分量如何進一步分解為下一級分量無關”。
從上述各例分析可以看出,處理合成標準不確定度評定中的相關問題時,不應忽視協方差項的影響,對協方差項中的相關系數的性質及計算等問題的正確處理是十分重要的。
|