陜西省計量測試研究所  許倩鈺

    在不確定度評定中，當全部輸入量估計值彼此獨立或不相關時，各分量的合成采用方和根計算，但當輸入量估計值出現(xiàn)相關時，由JJF1059-1999《測量不確定度評定與表示》第(22)、(23)、(24)式可知，不確定度合成運算中出現(xiàn)了不為零的協(xié)方差或相關系數(shù)，這使得計算變得更加復雜和難以進行(非理想狀態(tài))。在實際工作中，有些相關問題求解并不都是那么難以操作，在此以砝碼檢測為例來闡明。
    

　　一、砝碼不確定度相關實例

    依據(jù)JJG99-1990《砝碼》試行檢定規(guī)程，砝碼量傳方法有兩種，即直接比較法和組合比較法。直接比較法是用一個被檢砝碼與一個標稱值相同的標準砝碼比較，在這個量傳過程中，通常把各輸入量之間看成是彼此獨立的，與此有關的實例已見到很多。組合比較法是用一個砝碼與其標稱質量相當?shù)囊蝗喉来a比較，依次向上或向下累積，組合傳遞。在量傳過程中，這是個非常重要的不能用其它方法替代的方法。在組合比較法中部分輸入量之間就存在著無法回避的相關。如在大砝碼校準中，首先用25個20kg的二等砝碼組合成500kg，做為標準傳給兩個500kg的F₂等的砝碼，然后再由這兩個500kg的砝碼組合傳給1000kg的M₁₁級砝碼。本文就以M₁₁級1000kg砝碼校準為例，分析輸入量存在相關時不確定度的合成。
    1.不確定度來源
    在砝碼的校準中，不確定度的來源主要是:
    (1)標準砝碼帶來的不確定度
    (2)測量過程帶來的不確定度
    (3)環(huán)境影響帶來的不確定度
    在砝碼檢定時，不應有容易察覺的振動和氣流，并應遠離振源和磁源，根據(jù)JJG99-1990規(guī)程，在常溫、相對濕度不大于80%，并滿足上述檢定條件時，環(huán)境影響帶來的不確定度很小且對合成不確定度的貢獻可忽略不計，故在下面的分析中對此不做考慮。
    2.數(shù)學模型
    依據(jù)JJG99-1990規(guī)程，被測量砝碼的折算質量m^*_A由下面方程給出:
    

    式中:m^*_A——被檢砝碼的折算質量，g；m^*_B——標準砝碼的折算質量，g；V_A——被檢砝碼的體積，cm³；V_B——標準砝碼的體積，cm³；Δρ_K——室內的空氣密度ρ_K相對于約定的標準空氣密度，ρ_1.2=1.2mg/cm³的偏離量，Δρ_K=ρ_K-ρ_1.2=ρ_K-1.2mg/cm³；I_A——被檢砝碼在天平上的讀數(shù)，格；I_B——標準砝碼在天平上的讀數(shù)，格；d——天平的實際分度值，mg/格。
    3. 求靈敏系數(shù)
    為了計算方便，對公式(1)進行如下簡化
    

    則(1)式變?yōu)?m^*_A=m^*_B+m_ρ+L    (2)
    其關系式為:m^*_A=f(m^*_B，m_ρ，L)    (3)
    在(1)式中，由于m^*_Bi 是由同一標準傳遞的，故這幾個輸入量的估計值是相關的，但在(2)式中m^*_B，m_ρ，L三個量的估計值則是相互獨立的。    
    合成方差
    u²(m^*_A)=c₁²u²(m^*_B)+c₂²u²(m_ρ)+c₃²u²(L)    (4)
    

    則(4)式可寫為:
    u²(m^*_A)=u²(m^*_B)+u²(m_ρ)+u²(L)    (5)
    4.不確定度分量的評定
    (1)標準砝碼引起的不確定度分量u(m^*_B)
    在1000kgM₁₁級砝碼的檢定中，標準砝碼為兩個500kg的F₂等砝碼，即，其標準不確定度分別是u_B1、u_B2 ，因它們是由同一組二等標準砝碼傳遞的，這使m^*_B1 、m^*_B2可能同時偏大或偏小，即兩個砝碼的估計值明顯強相關，而且是正相關。這里的難題就是求相關系數(shù)，第一種方法是統(tǒng)計法，一般情況下(非理想狀態(tài))，用統(tǒng)計法求相關系數(shù)很困難；第二種方法是物理(實驗)判斷法，根據(jù)經(jīng)驗，當輸入量m^*_B1、m^*_B2為正強相關時，可近似取r(m^*_B1 ，m^*_B2 )=1，此時不確定度各分量的合成就是線性相加而不是方和根(JJF1059-1999中式(18)、(19))。
    即u(m^*_B)=u(m^*_B1 )+u(m^*_B2 )
    根據(jù)JJG99-1990規(guī)程，500kg質量允差為±7.5×10³mg，質量允差的1/3為擴展不確定度，則500kg F₂等砝碼的擴展不確定度是2.5×10³mg(k=3)，其標準不確定
    u(m^*_B1)=u(m^*_B2)=(2.5×10³/3)mg=(2.5/3)g    
    則標準砝碼引起的不確定度分量u(m^*_B)
    u(m^*_B)=u(m^*_B1)+u(m^*_B2)=(2.5/3)g+(2.5/3)g=(5/3)g
    (2)空氣浮力引起的標準不確定度分量u(m_ρ)
    根據(jù)JJF99-1990，在M₁₁級砝碼的檢定中，空氣浮力對砝碼檢定結果帶來的影響可以忽略不計，空氣浮力影響帶來的不確定度對合成不確定度的貢獻很小，在這里也忽略不計。
    (3)校準過程中的不確定度u(L)
    校準過程中的不確定度包括人員和天平帶來的不確定度，人員因素主要是天平的開、關以及讀數(shù)，這已在測量結果中反映，故人員因素帶來的不確定度在以下分析中不再單獨考慮。機械天平性能主要包括天平分度值、示值變動性、不等臂性(因砝碼檢定用的是替代法，不等臂性影響已消除)等，將這些因素綜合到一起，其對測量結果帶來的不確定度不超過5g，即u(L)=5g(有關這方面的文章較多，故此處不再做詳細說明)。
    5.不確定度合成
    

    6.擴展不確定度
    取包含因子k=2，則U=2u(m^*_A)=2×5.77g=11.5g，修約為U=12g
    

　　二、結束語

    在不確定度評定中，相關是個比較復雜的問題，但也不是無法操作，通過適當?shù)姆纸饨M合，可以將復雜問題簡單化。