标准曲线不确定度的分析.doc

标准曲线不确定度的分析曹 涛（深圳市通量检测科技有限公司，广东 深圳 518102）摘 要：本文阐述了标准曲线不确定度分析的通用方法，讲解了标准曲线不确定度的分析步骤和方法。并且通过检测过程中经典的几个项目的标准曲线的不确定度实例分析，来论证影响不确定度分析的因素；通过实验数据证明：参与标准曲线的点的个数和样品测定的次数越多，标准不确定度越小；并且样品测定的结果越靠近标准曲线的重点，标准不确定度越小； 标准曲线的线性越好，斜率越高，标准不确定度越好；关键词：标准曲线，不确定度中图分类号：R155.5 文献标识码：A 文章编号：作者简介：曹涛（1987-06-04），男，汉族，陕西省汉中市，中山大学本科毕业，研究食品中营养成分、添加剂、农兽药残留的检测分析；E-mail： Analysis of Standard Curves UncertaintyCao Tao（Shenzhen Total-Test Technology Co., Ltd, Shenzhen 518102）Abstract: This paper expounds a general method of the analysis of standard curves uncertainty and explains the steps of the analysis of standard curves 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 tests times of the standard curve and sample can make the standards uncertainty small. And The results of sample are more close to the center of the standard curve, the standards 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 标准曲线的不确定度忽略标准溶液的不确定度的引入用线性最小二乘法拟合曲线程序的前提是假定横坐标（标准溶液的浓度）量的不确定度远小于纵坐标的量的不确定度，因此通常的C0不确定度计算程序仅仅与响应值不确定度有关，而与校准溶液不确定度无关，也不与从同一溶液中逐次稀释产生必然的相关性1。1.2 标准曲线方差和残余方差的概念和计算假定一条标准曲线 （1.2.1）响应值y1，y2yn的偏差（）主要由两方面原因引入：一是自变量x取值不同导致的响应值的偏离均值（），另一方面是测量误差导致的响应值偏离标准曲线（）；如下图图1：标准曲线响应值偏差分析yi第xi点对应的响应值第xi点带入标准曲线（1.2.1）得到的理论响应值所有浓度响应值的平均值n参与标准曲线的浓度点和响应值的个数；（注：并不仅仅是标准曲线的点的个数，比如标准曲线有五个点，其中有两个点平行做了3次，则，同样计算的时候也需要把平行做的点计算在内，即x和y的平均值是参与本次标准曲线所有点的几何重心）在量化的过程中，用n个取值的偏离平方和来描述，这里分别记为：总偏差平方和S（反映了的总的分散程度）、回归平方和T（反映了回归值的分散程度，仅与标准曲线斜率相关）和残余平方和Q（反映了观测值偏离回归直线的程度，真实描述标准曲线不确定度）；由于交叉项结果为0（可参考最小二乘法计算线性回归方程公式）；因此 (1.2.3)其中， 因此S的自由度vS=vT+vQ；总偏差平方和的自由度由于受平均值的约束为n -1，而回归偏差平方和的自由度只与斜率b有关，自由度为1，因此残余偏差平方和的自由度vQ=n -2通过残余偏差平方和除以自由度得到对应的残余偏差 所以标准曲线的残余标准偏差 （1.2.4）1.3 标准曲线斜率和截距的不确定度分析由于得到的标准曲线中斜率b和截距a参与到未知样品的浓度的计算，因此必须对斜率b和截距a进行不确定度的分析；斜率的标准不确定度： （1.3.1）扩展不确定度U（b）= （1.3.2）式中tP是指选定置信空间P水平内，根据残余偏差自由度vQ查t分布表所得的t值；截距的标准不确定度： （1.3.3）扩展不确定度U（a）= （1.3.4） 式中tP是指选定置信空间P水平内，根据残余偏差自由度vQ查t分布表所得的t值1.4 代入未知样品响应值计算浓度过程不确定度分析标准曲线相关不确定度都已通过公式计算得到，而未知样品根据响应值再带入标准曲线得到浓度值的过程也会引入相应的不确定度，计算如下:被测物含量标准偏差： （1.4.1） 扩展不确定度U（x0）= （1.4.2）p被测物测定平行次数；被测物平行测定的平均值；xi标准曲线中各个浓度值；tP选定置信空间P水平内，根据残余偏差自由度vQ查t分布表所得的t值2 以原子吸收法测铅的标准曲线不确定度分析进行实例分析2表1 原子吸收法测铅标准曲线原始记录浓度值X01020304050吸光度yi0.01180.06290.11130.16230.20830.25310.01200.06300.11440.16280.20540.2523由上述原始记录可知下列数据：2.1：参与标准曲线的点n=122.2：线性曲线的方程：y=0.00482x+0.01458 ； 相关系数：r=0.99962.3：标准曲线的重心（25，0.1350），即，2.4：观测值与理论值一一对应关系如下： 表2 实测值与理论值对应关系表浓度值X观测值yi理论值观测值yi理论值00.01180.01460.01200.0146100.06290.06280.06300.0628200.11130.11100.11440.1110300.16230.15920.16280.1592400.20830.20740.20540.2074500.25310.25560.25230.25562.5：残余偏差 2.6：斜率标准不确定度 残余偏差自由度vQ=12-2=10，取置信空间为95%查t分布表的tP=2.228，斜率扩展不确定度U（b）=2.2280.00045=0.0010032.7：截距标准不确定度同上，截距扩展不确定度U（a）=2.2280.00136=0.003032.8：样品测得一次，吸光度分别为：0.1223，代入标准曲线过程得到的浓度为22.35，其标准不确定度如下：2.9：样品平行测得两次，吸光度分别为：0.1223和0.1315，代入标准曲线过程得到的浓度分别为22.35和24.26，其标准不确定度如下：2.10：样品平行测得三次，吸光度分别为：0.1223、0.1315、0.1254，代入标准曲线过程得到的浓度分别为22.35、24.26、22.99，其标准不确定度如下：2.11：样品平行测得四次，吸光度为：0.1223、0.1315、0.1254、0.1229，带入标准曲线得到的浓度分别为22.35、24.26、22.99、22.47，其标准不确定度如下2.12：样品平行测得五次，吸光度为：0.1223、0.1315、0.1254、0.1229、0.1304，带入标准曲线得到的浓度分别为22.35、24.26、22.99、22.47、24.03，其标准不确定度如下2.13：样品平行测得六次，吸光度分别为：0.1223、0.1315、0.1254、0.1229、0.1304、0.1309，带入标准曲线得到的浓度分别为22.35、24.26、22.99、22.47、24.03和24.13，其标准不确定度如下：2.14：样品平行测得十次，吸光度分别为：0.1223、0.1315、0.1254、0.1229、0.1304、0.1309、0.1244、0.1312、0.1321、0.1277，带入标准曲线得到的浓度分别为22.35、24.26、22.99、22.47、24.03和24.13，其标准不确定度如下：总结样品测定次数和不确定度的关系曲线如下：因此：样品测定一定要测定2-3次平行，超过四次对不确定度的减少意义并不大；3 紫外分光光度法测亚硝酸盐标准曲线的不确定度分析3表3 紫外分光光度法测亚硝酸盐标准曲线Q1原始记录浓度值X01.02.03.04.05.07.510.012.5吸光度yi00.01910.04810.07820.11320.14010.20940.28970.3596表4 紫外分光光度法测亚硝酸盐标准曲线Q2原始记录浓度值X01.02.03.04.05.07.510.012.5吸光度yi00.01410.04830.07490.10130.16100.22990.29410.3463由上述原始记录可知下列数据：3.1：参与标准曲线的点n=93.2：线性曲线Q1的方程：y=0.02932x-0.00691 ； 相关系数：r=0.9995 线性曲线Q2的方程：y=0.02938x-0.00570 ； 相关系数：r=0.99523.3：标准曲线Q1的重心（5.0，0.1397），即， 标准曲线Q2的重心（5.0，0.1412），即，3.4：观测值与理论值一一对应关系如下：表5 实测值和理论值关系表浓度值X标准曲线Q1标准曲线Q2观测值yi理论值观测值yi理论值00.0-0.00690.0-0.00571.00.01910.02240.01410.02412.00.04810.05170.04830.05383.00.07820.08110.07490.08364.00.11320.11040.10130.11335.00.14010.13970.16100.14317.50.20940.21300.22990.217410.00.28970.28630.29410.291812.50.35960.35960.34630.36623.5：标准曲线Q1残余偏差 标准曲线Q2残余偏差 3.6：标准曲线Q1斜率标准不确定度 残余偏差自由度vQ=9-2=7，取置信空间为95%查t分布表的tP=2.365， 标准曲线Q1斜率扩展不确定度U（b）=2.3650.000336=0.000795 标准曲线Q2斜率标准不确定度 同上，标准曲线Q2斜率扩展不确定度U（b）=2.3650.00112=0.002653.7：标准曲线Q1截距标准不确定度同上，标准曲线Q1截距扩展不确定度U（b）=2.3650.00215=0.00508标准曲线Q2截距标准不确定度同上，标准曲线Q2截距扩展不确定度U（b）=2.3650.00715=0.01693.8：测得样品平行两次，吸光度分别为：0.1554和0.1469，代入标准曲线Q1过程得到的浓度分别为5.54和5.34，其标准不确定度如下： 代入标准曲线Q2过程的到浓度分别为5.42和5.22，其标准不确定度如下：总结以上数据得出如下结论：线性较好的标准曲线不确定度较低，而且线性的好坏对不确定度的影响非常大！4 液相色谱法测三聚氰胺的标准曲线的不确定度分析4表6 液相色谱法测三聚氰胺的标准曲线原始记录浓度值X5.010.025.050.0100.0峰面积yi23959656316314579472632863552759424014557132214689712641009/由上述原始记录可知下列数据：4.1：参与标准曲线的点n=94.2：线性曲线的方程：y=55184x+10123 ； 相关系数：r=0.99964.3：标准曲线的重心（31.11，1726957），即，4.4：观测值与理论值一一对应关系如下：表7 实测值和理论值关系表浓度值x观测值yi理论值观测值yi理论值523959628604324014528604310.056316356196357132256196325.0145794713897231468971138972350.02732863276932327410092769323100.055275945528523/4.5：标准曲线残余偏差 4.6：标准曲线斜率标准不确定度 残余偏差自由度vQ=9-2=7，取置信空间为95%查t分布表的tP=2.365， 标准曲线斜率扩展不确定度U（b）=2.365565.24=1336.794.7：标准曲线截距标准不确定度同上，标准曲线截距扩展不确定度U（b）=2.36524202=572374.8：测得样品平行两次，峰面积分别为：1798791和1810244，代入标准曲线得到浓度分别为32.41和32.62，其标准不确定度如下： 4.9：测得另一组样