計量培訓:測量不確定度表述講座
國家質(zhì)量技術(shù)監督局 李慎安
3.1 實(shí)驗方差s2(qk)是方差σ2的無(wú)偏估計的含義為何?標準偏差s是否也是總體標準差σ的無(wú)偏估計?
在用貝塞爾公式
計算任一次測量結果qk的實(shí)驗標準偏差s(qk)時(shí),未開(kāi)方前以及用本講座2.12問(wèn)題中給出的式子計算時(shí),未開(kāi)方前,均稱(chēng)為實(shí)驗方差s2(qk)。σ稱(chēng)為總體標準偏差,σ2則稱(chēng)為總體方差或簡(jiǎn)稱(chēng)方差,在計量學(xué)中,特別是測量不確定度評定中,總體是指被測量Y在重復性條件下或復現性條件下無(wú)限多次的測量結果。根據這無(wú)限多次測量結果計算出的標準偏差就是σ。由于實(shí)驗中,重復的次數n總是有限的,計算出的實(shí)驗方差s2只是σ2的一個(gè)估計值。n越大,這個(gè)估計值越可靠。所謂無(wú)偏估計,可以簡(jiǎn)單地理解為:s2比σ2大的概率與s2比σ2小的概率相等,即均為50%。而且當次數n越大時(shí),差值(s2-σ2)的總和越趨近為零,當n為無(wú)窮大時(shí),s2-σ2就等于零。當s2是σ2的無(wú)偏估計時(shí),s就不是σ的無(wú)偏估計而是有偏的了,s是σ的偏小估計,即s-σ是負值的概率大于s-σ是正值的概率。在測量不確定度評定中,可以不去考慮這種偏小,因為隨n的增大它們會(huì )趨于相等。
3.2 為什么在按貝塞爾公式計算的實(shí)驗標準偏差時(shí),次數n應充分大(開(kāi)方后為什么只取正值)復現性條件下的重復測量結果可否采用貝塞爾公式計算一次測量結果的實(shí)驗標準偏差?
次數n越大,計算出來(lái)的實(shí)驗標準偏差s(qk)越可靠。一般文獻均提出應充分大,當然是越大越好,盡可能多地重復測量。不過(guò)一般來(lái)說(shuō),次數n≥30就認為充分了。因為n等于40或50雖比n=30好一點(diǎn),但好不了多少。當我們研究測量?jì)x器的特性,特別是其示值分布的情況時(shí),則是另一種目的,次數n往往要超過(guò)100甚至200。
數字中的平方根,總是帶有正負號的,例如=±10。但是,在標準差的計算中,開(kāi)方后的值只取正值,原因在于標準差表示的是分散性,而分散性所給出的是一個(gè)區間或理解為一個(gè)范圍。作為物理量的區間的大小,用負值是沒(méi)有意義的。這也就是測量不確定度只有正值而不存在負值的原因。
復現性條件下的重復測量結果之間,也存在分散性,這種條件下的任一次測量結果的實(shí)驗標準偏差,也無(wú)例外地可以用3.1或2.12中給出的式子進(jìn)行計算。復現性條件似乎不是等精度測量所要求的條件,但是,所謂等精度是個(gè)定性的概念,重復性條件下的測量結果有大有小,它們的測量誤差也各不相同,但應該說(shuō)是有限程度的等精度。復現性條件下出現了某些條件的變化,導致測量結果分散性某種程度的擴大,但仍可以用實(shí)驗標準偏差來(lái)定量表述,也可稱(chēng)之為等精度測量。參閱2.9。
3.3 測量不確定度的定義如何理解?
測量不確定度定義的英文為:Parameter,associated with the result of a measurement,that characterizes the dispersion of the values that could reasonably be attributed to the measurand。一般譯為:與測量結果相聯(lián)系的參數,用來(lái)表征合理地賦予被測量之值的分散性。上述譯文把“associated with”譯為“相聯(lián)系”不太貼切,英文的含義是“與…一起”,而“相聯(lián)系”一詞在漢語(yǔ)中,特別是在科技文獻中,往往令人要問(wèn),如何聯(lián)系,函數形式如何?其實(shí),在這里測量不確定度與測量結果之間的“聯(lián)系”,只不過(guò)是“在一起”,除此以外無(wú)其他含義。
當我們在重復性條件下,對一穩定的被測量X獨立進(jìn)行了n次重復測量,在這一測量列中,通過(guò)n個(gè)結果按貝塞爾公式計算出的,第i次結果xi的實(shí)驗標準差s(xi),與xi之間有怎樣的聯(lián)系?這里的xi雖指第i次測量結果,而其實(shí)際含義則為:任一次的測量結果。因此,s(xi)=u(xi)表明這個(gè)不確定度(分散性)是這個(gè)測量列中任意一次的結果的不確定度。當然,如果在相同的重復性條件下再測一次,得到的結果xi的標準不確定度同樣也是s(xi)。我們能看出這里的一次測量結果的標準不確定度u(xi)與xi之間有怎樣的聯(lián)系呢?沒(méi)有。
怎么叫合理?怎么是非合理?在《導則》中未予交代。有人說(shuō),這里“合理”一詞妙極,合理就是合理。沒(méi)有,也不必要有任何解釋。只要賦予被測量之值的分散性不能用不確定度來(lái)表征,則賦予被測量之值就不合理。如果是這樣,我們如何理解不確定度的概念呢?
定義中所謂的合理,是指處于統計控制狀態(tài)下的測量。當測量是處于統計控制狀態(tài)下時(shí),其結果的分散性才能用不確定度這一參數表征,否則不行。
國際上對實(shí)驗標準差的定義是表征結果分散性的,還有,1994年12月公布的國際標準ISO 5725—《測量方法與測量結果的準確度》其中對重復性標準差以及復現性標準差(sr與sR),都是明確規定重復性條件下和復現性條件下,對同一被測量獨立測量若干次的測量列,按貝塞爾公式所得到的分散性用標準差定量地給出的值,ISO分別用了標準化的符號sr與sR以示區別。
所謂統計控制狀態(tài)的含義,在統計學(xué)中就是指隨機控制過(guò)程狀態(tài)。在計量學(xué)中,一般來(lái)說(shuō),可以具體化為:重復性條件和復現性條件可以保證下的狀態(tài)。
當我們把測量過(guò)程中所用的標準測量?jì)x器,按證書(shū)所給的修正量或修正曲線(xiàn),對其示值(某些情況下就是測量結果)進(jìn)行修正后,由于修正值的不確定度導致的誤差,其期望是可以,而且往往只能,作為零來(lái)估計的,這就是一種統計控制狀態(tài),因它處于隨機過(guò)程之中。
不確定度是否就是測量結果的可能誤差?答復是肯定的。不確定度的含義雖為賦予被測量之值的分散性,但是,分散性的形成:一是隨機效應;二是系統效應。系統效應導致的誤差分量其期望(指對那些已知系統誤差進(jìn)行過(guò)修正后的)與隨機效應導致的誤差分量一樣,都是為零。因此,只要沒(méi)有遺漏重大的不確定度分量,最后給出的擴展不確定度,無(wú)論是U還是Up,都是一種可能誤差(possible error)的量度。事實(shí)上,在計量學(xué)中,過(guò)去給測量不確定度曾經(jīng)有過(guò)一個(gè)定義:由測量結果所給出的被測量估計值中,可能誤差的量度。這個(gè)定義雖已為1995年的《導則》放棄,但是,其概念與當前所采用的定義并不矛盾,可能誤差在大多數情況下,表達為一種誤差限,或最大允許誤差等。因此,我們在按檢定證書(shū)或某些儀器的技術(shù)規范中的這一指標,來(lái)估算其所導致的不確定度分量時(shí),就有理由把它們作為U或Up來(lái)對待。例如:證書(shū)上給出了最大允許誤差不超出±18μA,就可認為U99=18μA。而其標準不確定度在正態(tài)分布的前提下可估算為U99/3=18μA/3=6μA。
測量不確定度是否仍可理解為被測量真值所處范圍的量度?答案也是肯定的,JJF1001-1991中,曾對測量不確定度按當時(shí)國際上的意見(jiàn)定義為:表征被測量的真值所處量值范圍的評定。這一定義也為國際計量學(xué)界所放棄,原因是這兩個(gè)定義中均涉及到“真值”、“誤差”這樣的理論上的概念而不具有“可操作性”。雖然如此,其所表達的概念并未被國際計量學(xué)界所否定。德國于1996年3月所公布的標準DIN1319—3《單一被測量測量結果不確定度的估算》中,對測量不確定度的定義卻是采用了:和測量結果一起,用于說(shuō)明被測量真值所處范圍的一個(gè)參數。
不確定度與測量結果有多大的聯(lián)系?
例如:1個(gè)三等砝碼,交給某個(gè)實(shí)驗室,按檢定規程的要求進(jìn)行了測量,得到其質(zhì)量為m1。然后,把這個(gè)砝碼交給另一個(gè)實(shí)驗室,同樣按檢定規程進(jìn)行測量,得到其質(zhì)量為m2,這兩個(gè)實(shí)驗室各自使用自己的二等標準砝碼與天平,m1≠m2是常見(jiàn)的。但是,這兩個(gè)測量結果的不確定度是十分接近的,都不超過(guò)檢定規程的三等砝碼的要求。因此,只要測量程序、條件相同,不同的測量結果可以有相同的不確定度。反之,如果測量程序、條件并不相同,雖然測量結果相同,也未必有相同的不確定度。從這個(gè)意義上來(lái)看,測量不確定度獨立于測量結果。
應該認為:測量不確定度主要決定于測量程序與條件,而測量結果應是這一測量程序與條件下的測量結果而非其他。其聯(lián)系僅此而已。
不確定度指測量結果的可疑程度,即對測量結果正確性的可疑程度。其值大則表示不可靠,其值較小,則表示較為可靠,其準確度較高。
測量不確定度無(wú)例外地只用正值表述。例如:擴展不確定度U95=0.45 mA。如與測量結果用數學(xué)符號聯(lián)系起來(lái),則另加正負號(±)。例如:電流I=(70.000±0.054)A。
[page_break]
3.4 不確定度分為哪些類(lèi)?
不確定度的定義與概念已如3.3所述,當不確定度除以真值(或測量結果)時(shí),稱(chēng)之為相對不確定度。這是個(gè)無(wú)量綱量,通常用百分數或10的負數冪表示(例如10-6,10-9)等,而其符號則加下標rel,例如:U95rel。
不帶形容詞的不確定度指一般概念,當需要明確某一測量結果的某種不確定度時(shí),要適當采用一個(gè)形容詞,常用的形容詞有:標準、擴展(展伸或范圍)。在這兩個(gè)形容詞前,還可再加“相對”。例如:相對擴展不確定度Urel。
當不確定度是用標準偏差表述時(shí),稱(chēng)為標準不確定度。
當不確定度是采用統計方法(例如用2.12或3.1給出的計算式)得到時(shí),稱(chēng)A類(lèi)評定,而得出的不確定度稱(chēng)為A類(lèi)標準不確定度。當用不同于統計方法的其他方法得到時(shí),稱(chēng)為B類(lèi)評定,所得出的不確定度稱(chēng)為B類(lèi)標準不確定度。由各個(gè)不確定度的方差和協(xié)方差之和算出的標準不確定度,稱(chēng)為合成標準不確定度,它是測量結果標準偏差的估計。
擴展不確定度則是合成標準不確定度乘一個(gè)包含因子(見(jiàn)3.5)之后擴大了兩倍或兩倍以上的不確定度,它給出的區間能包含被測量可能值的大部分(例如:95%,99%等)。但也可以不必給出其百分數。
不確定度的分類(lèi)可用下圖表明:
上圖中所給出的各類(lèi)不確定度均可對應地給出它們的相對不確定度,即再除以測量結果。
決不可以用隨機不確定度和系統不確定度這樣的概念和術(shù)語(yǔ)。因不確定度只是一個(gè)分散性區間,這個(gè)區間沒(méi)有什么隨機與系統性的問(wèn)題,更不存在系統的分散性或系統的分散區間。如需要表明某不確定度是怎樣的原因所導致的,可以用隨機效應導致的不確定度或系統效應導致的不確定度,以避免概念上的混淆。必須注意,隨機效應導致的不確定度既可以是A類(lèi),也可以是B類(lèi),而系統效應導致的不確定度也不一定都是B類(lèi)。
3.5 包含因子k的概念如何?用在不確定度評定中的什么場(chǎng)合?
包含因子又稱(chēng)覆蓋因子,它等于擴展不確定度與合成標準不確定度uc之比。由于擴展不確定度有U與Up兩類(lèi)(下標p為置信概率,Up給出置信概率為p的置信區間半寬)。包含因子亦有k與kp兩類(lèi),k與kp在名稱(chēng)上沒(méi)有區別,使用時(shí),如需區別,必須給出符號。k一般為2,有時(shí)為3。而kp則是在給定概率p時(shí)所要求的因子,對于被測量Y之估計值接近正態(tài)分布的情況下,kp就是t分布給出的t值。在擴展不確定度Up的計算中,有了p與自由度ν即可查表得到kp,從而得到Up。如只要求給出擴展不確定度U,則無(wú)需kp而只用k即可。
3.6 自由度的含義如何?不確定度評定中起何作用?什么情況下可以不必進(jìn)行自由度的估算?
自由度一詞在不同學(xué)科有不同定義,例如在運動(dòng)學(xué)中,指物體活動(dòng)空間是幾維。質(zhì)點(diǎn)只能沿固定軌道進(jìn)退為1維,在一個(gè)面上自由運動(dòng)則是2維,在空間自由運動(dòng)則是3維,分別稱(chēng)為自由度等于1,2,3。剛體加上旋轉的自由,最大等于6。自由度從字面上看是指松動(dòng)程度。計量學(xué)中,如一個(gè)被測量只測了一次,有一個(gè)結果,不存在選擇余地,自由度為零,但有了兩個(gè)測量結果,就多了一個(gè)選擇。不確定度的評定中,自由度用于表明所得到的標準偏差的可靠程度,它被定義為方差計算中和的項數減對和的限制數。按貝塞爾公式(見(jiàn)3.1)計算方差s2(qi)時(shí),∑符號后的項數等于n,因為重復了n次,被測量Q有n個(gè)結果,與其平均值之差就有n個(gè),成為n項之和。但有一個(gè)限制,就是由于這n個(gè)殘差之和必為零,即∑(qi-
)=0,這算一個(gè)限制條件。自由度ν=n-1。
一般,我們可以認為自由度等于測量次數n減被測量的個(gè)數m,即ν=n-m。
自由度越大,這個(gè)標準差越可靠,自由度ν與這個(gè)標準差s的相對不確定度的平方成反比。
自由度只用于包含因子kp的獲得,如果在擴展不確定度的評定中只要求U而不是Up,則不必進(jìn)行自由度的評定及有效自由度的計算。合成標準不確定度的自由度稱(chēng)為有效自由度νeff。
3.7 置信概率的含義如何?與置信區間有何關(guān)系?
按測量不確定度的定義,合理賦予被測量之值的分散區間是包括全部被測量的測量結果的,即測量結果100%存在于這一區間。這一分散區間的半寬一般用a表示。但是如只要求某個(gè)區間只包含其95%的賦予被測量之值,這個(gè)區間就稱(chēng)為概率p=95%的置信區間,其半寬就是擴展不確定度U95,如要求99%的概率,則為U99。相應的概率稱(chēng)為置信概率,有:
U95<U99<a
至于大多少,與賦予被測量之值的分布情況有關(guān)。
3.8 測量誤差的定義、分類(lèi)以及使用時(shí)應注意哪些問(wèn)題?
測量誤差(簡(jiǎn)稱(chēng)誤差)的定義從20世紀70年代以來(lái)沒(méi)有改變,定義為:測量結果減被測量的真值。但是,長(cháng)期以來(lái),國內外計量學(xué)界常錯誤地使用誤差一詞。從定義看,誤差與測量結果有關(guān)而與測量方法無(wú)關(guān)。不同的測量結果有不同誤差,相同測量結果有相同誤差,而不論測量結果是來(lái)自何種測量方法。合理賦予被測量之值,各有其誤差而并不存在一個(gè)共同的誤差。一個(gè)測量結果的誤差,如不是正值就是負值,取決于這個(gè)結果是大于還是小于真值。因此,誤差決不會(huì )帶有正負號(±)。
測量結果的誤差往往由若干分量組成,這些分量按其特性分成隨機誤差與系統誤差兩大類(lèi),而測量結果的誤差無(wú)例外地是全部分量的代數和。即,對誤差的合成只有代數和這一種方式。
隨機誤差的定義在1993年以來(lái)作了原則性的改變,它被定義為:測量結果減重復性條件對同一量進(jìn)行無(wú)限多次測量結果的平均值。測量結果是真值、系統誤差與隨機誤差這三者的代數和,而無(wú)限多次結果的平均值則只是真值與系統誤差的代數和。它們的差則是這一測量結果的隨機誤差分量。
注意:不再有偶然誤差這一術(shù)語(yǔ),也不再有另外的定義。
系統誤差也有了原則性的改變,它被定義為:重復性條件下對被測量的無(wú)限多次測量結果的算術(shù)平均值減被測量的真值。由于只能有限次數的重復,而真值只能用給定真值代替,因此,所得到的系統誤差只是個(gè)估計值,并具有一定的不確定度。這個(gè)不確定度也就是修正值的不確定度,與其他來(lái)源的不確定度分量一樣進(jìn)入合成標準不確定度,僅此而已。而不是把系統誤差分成為已知系統誤差和未知系統誤差,也不能說(shuō)未知系統誤差按隨機誤差處理。因為這里所謂的未知系統誤差并非誤差分量而是不確定度,而且,所謂按隨機誤差處理,概念是不清的。
至于誤差限、最大允許誤差、可能誤差、引用誤差等術(shù)語(yǔ),它們前面帶有正負號(±),是一種可能誤差的分散區間,而并非一個(gè)測量結果的誤差。
過(guò)去所謂的誤差傳播定律,所傳播者并非誤差而是不確定度?,F在已改稱(chēng)為不確定度傳播定律。
應該注意,誤差一詞只能按其定義使用。今后不應用它來(lái)定量表明測量結果的可靠程度。
3.9 測量誤差與測量不確定度之間存在哪些主要不同之處?
下表給出它們之間的主要不同之處:
3.10 什么叫變量之間的相關(guān)?如何定量表達?不確定度評定中如何給出其估計值,什么情況下發(fā)生相關(guān)?有沒(méi)有B類(lèi)估算?
一切被測量的估計值,由于諸多不穩定因素的影響,它們不是固定的,因而稱(chēng)之為變量。也就是說(shuō)量的測量結果都是變量。
當某些被測量的估計值有相同的不確定度來(lái)源,特別是受相同的系統效應的影響,例如用了同一個(gè)標準器,則這樣的變量間存在相關(guān)。均可能偏大或均可能偏小,稱(chēng)為正相關(guān);一個(gè)偏大而另一個(gè)偏小,稱(chēng)為負相關(guān)。這種相關(guān)性導致的方差稱(chēng)為協(xié)方差,進(jìn)入合成方差的計算,從而擴大了所得的合成方差(合成方差是合成標準差的二次方)。
協(xié)方差的評定既有A類(lèi)評定也有B類(lèi)評定。往往也可通過(guò)測量的操作程序來(lái)避免相關(guān)的產(chǎn)生,即使其協(xié)方差小到可忽略不計,例如通過(guò)改變所用的標準器等。