張國(guó)城 張慶暖
(北京市計(jì)量檢測(cè)科學(xué)研究院,北京 100029)
?
對(duì)分析儀器兩點(diǎn)線性校準(zhǔn)引入不確定度的探討
張國(guó)城張慶暖
(北京市計(jì)量檢測(cè)科學(xué)研究院,北京 100029)
摘要:兩點(diǎn)線性校準(zhǔn)是分析儀器常用的校準(zhǔn)方法,該校準(zhǔn)過(guò)程所引入的不確定度是儀器測(cè)量結(jié)果不確定度的組成部分。但是對(duì)于兩點(diǎn)線性校準(zhǔn)引入的不確定度,卻未見(jiàn)有文獻(xiàn)進(jìn)行過(guò)評(píng)定或分析。本文在兩點(diǎn)校準(zhǔn)原理的基礎(chǔ)上,通過(guò)詳細(xì)的數(shù)學(xué)推導(dǎo),分析了各參數(shù)對(duì)測(cè)量結(jié)果不確定度的影響,發(fā)現(xiàn)響應(yīng)值y引入的不確定度對(duì)測(cè)量結(jié)果的不確定度沒(méi)有貢獻(xiàn),兩點(diǎn)線性校準(zhǔn)的過(guò)程是將標(biāo)準(zhǔn)溶液x的不確定度按一定規(guī)律傳遞給被測(cè)對(duì)象。
關(guān)鍵詞:計(jì)量;不確定度;兩點(diǎn)校準(zhǔn);線性;評(píng)定
分析儀器大部分是建立在間接測(cè)量的基礎(chǔ)之上,即通過(guò)檢測(cè)與被分析物的物理或化學(xué)屬性具有特定函數(shù)關(guān)系的物理量,來(lái)實(shí)現(xiàn)被測(cè)物相關(guān)物化參數(shù)的測(cè)量,例如利用利用電學(xué)或光學(xué)信號(hào)的強(qiáng)度來(lái)間接測(cè)量物質(zhì)的濃度或質(zhì)量的大小。這種特定的函數(shù)關(guān)系,在很多儀器中,都認(rèn)為是線性的。
那么在進(jìn)行未知樣品檢測(cè)之前,需要先用已知濃度或質(zhì)量(用x表示)的標(biāo)準(zhǔn)物質(zhì)建立起濃度或質(zhì)量與信號(hào)(用y表示)的線性關(guān)系(限于直線線性段),即:
y=a+bx
(1)
其中a和b分別為校準(zhǔn)曲線的截距和斜率。如果校準(zhǔn)曲線共有n個(gè)溶液,其中n≥3時(shí),稱為線性回歸,其線性回歸引入的不確定度評(píng)定已有大量文獻(xiàn)報(bào)道[1-3]。但是當(dāng)n=2時(shí),由于兩點(diǎn)確定一條直線,利用兩點(diǎn)線性校準(zhǔn)就可確定儀器的濃度-響應(yīng)值的線性關(guān)系。利用兩點(diǎn)實(shí)現(xiàn)儀器校準(zhǔn)的應(yīng)用很廣泛。對(duì)于兩點(diǎn)線性校準(zhǔn)引入的不確定度,由于不存在最小二乘法線性擬合,不能利用線性回歸的方法對(duì)其進(jìn)行不確定度的評(píng)定。而且此類儀器的不確定度評(píng)定案例中,鮮有考慮到線性校準(zhǔn)所引入的不確定度。那么兩點(diǎn)線性校準(zhǔn)引入的不確定度是否可忽略?本文將就此問(wèn)題進(jìn)行深入分析,為廣大計(jì)量工作者提供參考。
1不確定度的評(píng)定
假設(shè)用來(lái)校準(zhǔn)的兩點(diǎn)分別為(x1,y1)、(x2,y2),則直線斜率b和截距a分別為:
(2)
(3)
假設(shè)利用校準(zhǔn)后的儀器檢測(cè)某溶液得到的響應(yīng)值為yobs,則通過(guò)公式1可求得對(duì)應(yīng)的測(cè)量結(jié)果xpred:
(4)
將公式2和公式3代入公式4可得:
(5)
根據(jù)JJF 1059-1999[4]和JJF1135-2005[5],通過(guò)公式5,可得xpred的不確定度為:
(6)
其中:
(7)
(8)
(9)
(10)
公式(6)充分考慮了兩點(diǎn)線性校準(zhǔn)所引入的不確定度。由于y1、y2和yobs都是利用同一儀器測(cè)量的,則可認(rèn)為三者均為正相關(guān),即r(y1,y2)=1,r(yobs,y1)=1,r(yobs,y2)=1。
若x1和x2是利用相同方法定值的,則可認(rèn)為正相關(guān),r(x1,x2)=1。
則公式(8)和公式(9)可簡(jiǎn)化為:
(11)
(12)
從公式(5)可得到:
(13)
(14)
(15)
(16)
(17)
所以:
(18)
(19)
-(yobs-y1)u(y2)]2
(20)
-(yobs-y1)u(y2)]
(21)
若x1和x2是相互獨(dú)立的,即r(x1,x2)=0,則公式19應(yīng)為:
(22)
2討論
若u(y1)=u(y2) =u(yobs),則
(23)
(24)
(25)
由于同一臺(tái)儀器測(cè)量,u(y1)=u(y2) =u(yobs)。從公式(25)可看出,響應(yīng)值的不確定度對(duì)u(xpred)的貢獻(xiàn)為零。所以在兩點(diǎn)法線性校準(zhǔn)的儀器中,無(wú)需考慮響應(yīng)值不確定度的影響。此類儀器一般也可不給出響應(yīng)信號(hào)y的值,而直接給出測(cè)量結(jié)果x的值
(1)若x1和x2正相關(guān),且u(x1)=u(x2),則
(26)
u(xpred)=u(x1)
(27)
如果校準(zhǔn)溶液使用相同的方法定值,且使用相同的天平、移液管、容量瓶等進(jìn)行稀釋制備,則可認(rèn)兩種校準(zhǔn)溶液的不確定度相同。此時(shí)兩點(diǎn)線性校準(zhǔn)的結(jié)果是將標(biāo)液的不確定度直接傳遞給測(cè)量結(jié)果。
(2)若x1和x2正相關(guān),但u(x1)≠u(x2),則
(28)
這種情況可能對(duì)應(yīng)于使用不同濃度或廠家的GBW標(biāo)準(zhǔn)物質(zhì)配置,或使用不同量程的天平、移液管、容量瓶,導(dǎo)致配置的標(biāo)準(zhǔn)溶液定值不確定度有所區(qū)別。即若xobs離x1越近,則u(xpred)受u(x1)影響越大,u(x1)和u(x2)對(duì)u(xpred)的影響與xpred到x2和x1的距離成正比。
仍然以酸度計(jì)為例,如果儀器在25℃條件下,分別以6.86pH(混合磷酸鹽pH標(biāo)準(zhǔn)物質(zhì),GBW(E)130071,U=0.01)和9.18 pH(硼砂pH溶液標(biāo)準(zhǔn)物質(zhì),GBW(E)130077, U=0.01)進(jìn)行兩點(diǎn)校準(zhǔn),由于兩個(gè)校準(zhǔn)點(diǎn)的不確定度均為0.01pH,則校準(zhǔn)過(guò)程引入的U=0.01pH。如果選用的9.18 pH溶液標(biāo)準(zhǔn)物質(zhì)更改為GBW(E)130469, 其U=0.03,由于兩個(gè)校準(zhǔn)點(diǎn)的不確定度不同,根據(jù)公式(28),如果被測(cè)對(duì)象為7.00pH,則校準(zhǔn)過(guò)程引入的不確定度為0.011pH;如果被測(cè)對(duì)象為9.00pH,則校準(zhǔn)過(guò)程引入的不確定度為0.028pH。即在儀器量程范圍內(nèi),不同點(diǎn)評(píng)估得到的不確定度不同,離哪個(gè)校準(zhǔn)點(diǎn)越近,受該點(diǎn)不確定影響越大。
(3)若x1和x2是不相關(guān)的,則
(29)
這種情況可能對(duì)應(yīng)于兩個(gè)校準(zhǔn)點(diǎn)使用的標(biāo)準(zhǔn)物質(zhì)是利用兩種不同原理或組分制備的,使得x1和x2沒(méi)有任何關(guān)聯(lián)性,但這種情況基本不存在。
與多點(diǎn)線性回歸不確定度評(píng)定類似,也可以根據(jù)公式4,分別評(píng)定斜率b和截距a的不確定度,然后再來(lái)合成xpred的不確定度:
但是,由于是兩點(diǎn)校準(zhǔn),兩點(diǎn)確定一條直線,所以斜率b和截距a也是確定的,不存在不確定度,由統(tǒng)計(jì)學(xué)方法推導(dǎo)得到的不確定度公式不適用。因此,這種方法是不存在的。
3結(jié)論
分析結(jié)果表明,由于利用同一臺(tái)儀器進(jìn)行校準(zhǔn)和測(cè)量,所以響應(yīng)值y引入的不確定度互相抵消,對(duì)測(cè)量結(jié)果的不確定度沒(méi)有貢獻(xiàn),而起貢獻(xiàn)作用的是校準(zhǔn)溶液的不確定度。所以很多利用兩點(diǎn)法進(jìn)行校準(zhǔn)的儀器不給出響應(yīng)值的做法,不僅方便儀器使用者,而且具有科學(xué)合理性。若兩個(gè)校準(zhǔn)溶液的不確定度相同,如采用相同的定值方法,包括使用相同的稱量、稀釋等設(shè)備時(shí),兩點(diǎn)線性校準(zhǔn)的結(jié)果是將校準(zhǔn)溶液的不確定度直接傳遞給測(cè)量結(jié)果;若兩個(gè)校準(zhǔn)溶液的不確定度不同,則被測(cè)點(diǎn)離哪個(gè)校準(zhǔn)點(diǎn)越近,則該校準(zhǔn)點(diǎn)的不確定度貢獻(xiàn)越大,與距離成線性比例。該結(jié)論表明,利用兩點(diǎn)線性校準(zhǔn)法的儀器,一般不用考慮校準(zhǔn)過(guò)程所引入的不確定度,而只需要將校準(zhǔn)液的不確定度直接傳遞測(cè)量結(jié)果。該結(jié)論對(duì)于儀器使用者評(píng)定檢測(cè)結(jié)果不確定度具有重要的參考意義。
參考文獻(xiàn)
[1]中國(guó)合格評(píng)定國(guó)家認(rèn)可中心,等.材料理化檢驗(yàn)測(cè)量不確定度評(píng)估指南及實(shí)例(CNAS-GL 10:2006)[M]. 北京:中國(guó)計(jì)量出版社,2007:6-8.
[2]李慎安,等.化學(xué)實(shí)驗(yàn)室測(cè)量不確定度[M].北京:化學(xué)工業(yè)出版社,2006:76-78.
[3]中國(guó)合格評(píng)定國(guó)家認(rèn)可委員會(huì).化學(xué)分析中不確定度評(píng)估指南(CNAS-GL 2:2006)[M].北京:中國(guó)計(jì)量出版社,2002:133-134.
[4] JJF 1059-1999, 測(cè)量不確定度評(píng)定與表示技術(shù)規(guī)范[S].
[5] JJF1135-2005, 化學(xué)分析測(cè)量不確定度評(píng)定技術(shù)規(guī)范[S].
Uncertainty evaluation of analytical instrument caused by two-point linear calibration.
ZhangGuocheng,ZhangQingnuan
(BeijingInstituteofMetrology,Beijing, 100029,China)
Abstract:In this paper, we discuss the uncertainty of the measuring results introduced by the two-point linear calibration. It is found that the uncertainty of response intensity exerts no effect on the uncertainty of the measuring results. In comparison, the uncertainty caused by two-point linear calibration is mainly introduced by the uncertainty of the two standard solutions used in calibration according to specific regulation.
Key words:measurement; uncertainty; two-point calibration; linear; evaluation
收稿日期:2014-11-10
DOI:10.3936/j.issn.1001-232x.2015.02.014
作者簡(jiǎn)介:張國(guó)城,男,1980年出生,北京市計(jì)量檢測(cè)科學(xué)研究院工作,室主任,博士,高級(jí)工程師,主要從事化學(xué)分析儀器的計(jì)量檢測(cè)和研究工作,Email: zhanggc@bjjl.cn。