标准曲线不确定度的分析
曹 涛
(深圳市通量检测科技有限公司,广东 深圳 518102)
摘 要:本文阐述了标准曲线不确定度分析的通用方法,讲解了标准曲线不确定度的分析步骤和方法。并且通过检测过程中经典的几个项目的标准曲线的不确定度实例分析,来论证影响不确定度分析的因素;通过实验数据证明:参与标准曲线的点的个数和样品测定的次数越多,标准不确定度越小;并且样品测定的结果越靠近标准曲线的重点,标准不确定度越小; 标准曲线的线性越好,斜率越高,标准不确定度越好; 关键词:标准曲线,不确定度
中图分类号:R155.5 文献标识码:A 文章编号:
作者简介:曹涛(1987-06-04),男,汉族,陕西省汉中市,中山大学本科毕业,研究食品中营养成分、添加剂、农兽药残留的检测分析;E-mail:vitcy@yahoo.cn
Analysis of Standard Curve’s Uncertainty
Cao Tao
(Shenzhen Total-Test Technology Co., Ltd, Shenzhen 518102)
Abstract: This paper expounds a general method of the analysis of standard curve’s uncertainty and explains the steps of the analysis of standard curve’s uncertainty. Prove the effect of the factors of uncertainty analysis through the testing process in classical several items of standard curves of uncertainty analysis. Experiments demonstrate: Improving the test’s times of the standard curve and sample can make the standard’s uncertainty small. And The results of sample are more close to the center of the standard curve, the standard’s uncertainty is smaller. The standard curve is linear in the better, and the slope is higher, the standard uncertainty can be smaller. Key word: the standard curve. uncertainty
前言:仪器分析中线性回归标准曲线测定方法,利用被测物质相应的信号强度与其浓度成正比关系,通过测定已知浓度的溶液(即标准溶液)的信号强度,通过最小二乘法将响应值和浓度的对应的线性关系拟合成一条直线,再根据未知样品的响应值推算出对应的浓度;然而测得的所有点未必全部都落在标准曲线上(除非线性相关系数r=1),因此得到的标准曲线本身具备相应的不确定性,而通过标准曲线去计算得到的浓度值就不可避免的具备不确定性,而且这个不确定性往往是整个实验不确定度的最大来源;因此对标准曲线计算不确定度非常有必要; 1 标准曲线不确定度分析的概念和计算
1.1 标准曲线的不确定度忽略标准溶液的不确定度的引入
1
用线性最小二乘法拟合曲线程序的前提是假定横坐标(标准溶液的浓度)量的不确定度远小于纵坐标的量的不确定度,因此通常的C0不确定度计算程序仅仅与响应值不确定度有关,而与校准溶液不确定度无关,也不与从同一溶液中逐次稀释产生必然的相关性。
[1]
1.2 标准曲线方差和残余方差的概念和计算
假定一条标准曲线
y?a?bx (1.2.1)
响应值y1,y2…yn的偏差(yi-y)主要由两方面原因引入:一是自变量x取值不同导致的响应值的
?);如下图 -y?-y)偏离均值(y,另一方面是测量误差导致的响应值偏离标准曲线(yi
图1:标准曲线响应值偏差分析
yi——第xi点对应的响应值
?——第xi点带入标准曲线(1.2.1)得到的理论响应值 yy——所有浓度响应值的平均值
n——参与标准曲线的浓度点和响应值的个数;(注:并不仅仅是标准曲线的点的个数,比如标准曲线有五个点,其中有两个点平行做了3次,则n?3?2?1?1?1?9,同样计算x的时候也需要把
平行做的点计算在内,即x和y的平均值是参与本次标准曲线所有点的几何重心)
在量化的过程中,用n个取值的偏离平方和来描述,这里分别记为:总偏差平方和S(反映了yi的总的分散程度)、回归平方和T(反映了回归值y?i的分散程度,仅与标准曲线斜率相关)和残余平方和Q(反映了观测值yi偏离回归直线的程度,真实描述标准曲线不确定度);
2
?i)?i-y)???i)???i-y)?2??i)?i-y)S??(yi-y)??(yi-y?(y(yi-y(y(yi-y?(y222i?1i?1i?1i?1i?1nn??2nnn由于交叉项
(y?i?1nni?i)?i-y)结果为0(可参考最小二乘法计算线性回归方程公式); -y?(yn因此
2?i)???i-y)S??(yi-y(y?Q?T (1.2.3)
2i?1i?1其中Q2?i), T?? ?i-y)??(yi-y(yn2ni?1i?1因此S的自由度vS=vT+vQ;总偏差平方和的自由度由于受平均值的约束为n -1,而回归偏差平方和的自由度只与斜率b有关,自由度为1,因此残余偏差平方和的自由度vQ=n -2
通过残余偏差平方和除以自由度得到对应的残余偏差
s2?Q n-2所以标准曲线的残余标准偏差s?(y?i?1ni2?i)-yn-2 (1.2.4)
1.3 标准曲线斜率和截距的不确定度分析
由于得到的标准曲线中斜率b和截距a参与到未知样品的浓度的计算,因此必须对斜率b和截距a进行不确定度的分析;
斜率的标准不确定度:s(b)?s(x?i?1ni2-x) (1.3.1)
扩展不确定度U(b)=tp (1.3.2) ?s(b)式中tP是指选定置信空间P水平内,根据残余偏差自由度vQ查t分布表所得的t值;
n截距的标准不确定度:
s(a)?s??xi?1ni?12i2n??(xi-x)1?s??nx(x?i?1n2 (1.3.3)
2-x)i扩展不确定度U(a)=tp (1.3.4) ?s(a)式中tP是指选定置信空间P水平内,根据残余偏差自由度vQ查t分布表所得的t值
3
1.4 代入未知样品响应值计算浓度过程不确定度分析
标准曲线相关不确定度都已通过公式计算得到,而未知样品根据响应值再带入标准曲线得到浓度值的过程也会引入相应的不确定度,计算如下: 被测物含量标准偏差:
(xs?x)2 (1.4.1) 扩s11s(x0)????2nbpn?(xi?x)i?1展不确定度U(x0)=tp (1.4.2)p?s(x0)——被测物测定平行次数;
xs——被测物平行测定的平均值;
xi——标准曲线中各个浓度值;
tP——选定置信空间P水平内,根据残余偏差自由度vQ查t分布表所得的t值 2 以原子吸收法测铅的标准曲线不确定度分析进行实例分析
表1 原子吸收法测铅标准曲线原始记录
[2]
浓度值X 0 0.0118
10 0.0629 0.0630
20 0.1113 0.1144
30 0.1623 0.1628
40 0.2083 0.2054
50 0.2531 0.2523
吸光度yi
0.0120
由上述原始记录可知下列数据: 2.1:参与标准曲线的点n=12
2.2:线性曲线的方程:y=0.00482x+0.01458 ; 相关系数:r=0.9996 2.3:标准曲线的重心(25,0.1350),即x?25,y?0.1350
2.4:观测值与理论值一一对应关系如下:
表2 实测值与理论值对应关系表
浓度值X
0 10 20 30 40 50
观测值yi 0.0118 0.0629 0.1113 0.1623 0.2083 0.2531
? 理论值y0.0146 0.0628 0.1110 0.1592 0.2074 0.2556
观测值yi 0.0120 0.0630 0.1144 0.1628 0.2054 0.2523
? 理论值y0.0146 0.0628 0.1110 0.1592 0.2074 0.2556
4