线性回归标准曲线法不确定度(检验检疫)

仪器 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 中线性回归 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 曲线法分析结果不确定度评估 一、前言 对测试方法制定不确定度评估程序是ISO/IEC 17025对实验室的要求[1]，也是检验工作的需要。由ISO等7个国际组织联合发布的《测量不确定度表达指南》[2]采用当前国际通行的观点和方法，使涉及测量的技术领域和部门可以用统一的准则对测量结果及其质量进行评定、表示和比较，满足了不同学科之间交往的需要[3]。采用《测量不确定度表达指南》对测试结果不确定度进行评估，也是检验工作同国际标准接轨的需要。 线性回归标准曲线法是仪器分析中最常用的方法，这类仪器包括原子吸收分光光度计、发射光谱仪、分光光度计、气相（液相）色谱仪等。这类分析测定结果的不确定度都有相似的来源，可概括为仪器精密度、标准物质不确定度及溶液制备过程中带来的不确定度等。因此，可用相似的方法对它们进行评估。本文以ICP-AES法测定钢铁中磷为例，推导了仪器分析中线性回归标准曲线法测定不确定度的计算方法，并提供了计算过程所需的各参数的采集和计算方法，评估了标准不确定度、自由度和扩展不确定度的数值。 二、测定过程和数学模型 仪器分析中线性回归标准曲线测定方法，利用被测物质相应的信号强度与其浓度成正比关系，通过测定已知浓度的溶液（即标准溶液）的信号强度，回归出浓度-信号强度标准曲线，从标准曲线上得到被测定溶液信号强度相应的浓度。计算过程的数学模型如下： 用yi和yt分别表示标准溶液和被测溶液的信号线强度，以xi和xt分别表示第i个标准溶液和被测样品溶液的浓度，i=1~n，n表示标准溶液个数，则： (1) 其中， (2) (3) （1）式也可表示成： (4) 把式（2）、（3）代入式（4）得： （5） 式（5）表明了被测量xt与输入量x1,x2...xn和y1,y2...yn、yt的函数关系，可简写成： 由上式可知，样品溶液浓度测定结果不确定度可分成标准溶液浓度不确定度分量及其信号强度不确定度分量和被测定溶液信号强度不确定度分量，其中标准溶液浓度不确定度分量可由标准样品标称含量不确定度和配制过程引入的不确定度合成得到，而信号强度不确定度分量是由仪器测量的误差引起的，可从仪器的精密度数据得到。如果标准溶液是从各自不同的母标准物质配制，并且在配制过程使不同的器皿，则各标准溶液之间浓度估计值之间无关，因此xt各输入量估计值之间无关，由文献[2]5.1款式10，得样品溶液浓度的合成标准不确定度如下： （6） 由式（5）得传播系数的表达式如下： (7) (8) (9) 我们以测定钢铁中磷的浓度为例，选择5个标准钢铁样品，按SN/T0750-1999《进出口碳钢、低合金中铝、砷、铬、钴、铜、磷、钼、镍、硅、锡、钛、钒含量的测定-电感耦合等离子体原子发射光谱（ICP-AES）法》的方法制备标准溶液和被测样品溶液。在ICP-AES上选用P213.6nm做分析谱线，分别测定标准溶液和被测溶液的分析谱线强度，用最小二乘线性回归法做出浓度-分析谱线强度标准曲线，从标准曲线上得到与被溶液分析谱线强度相应的浓度，再把被测溶液浓度转换成样品的百分比含量。测试过程的各种量如表1所示。 表1 钢铁中磷的测定有关量值 样品编号 标称含量（%） 称样量（g） 溶液浓度（ug/ml） 谱线强度（cps） BH0507-9A 0.011 0.4755 0.523 57.6 30#第39号 0.022 0.4864 1.07 139.6 GB01202 0.029 0.5167 1.39 190 GB01203 0.038 0.5053 1.92 266.7 BH0506-3A 0.063 0.4524 2.58 351.9 被测样品   0.4805 2.30 316.7           三、各量值不确定度的计算 如前所述，仪器分析中标准曲线法测定结果不确定度来源包括仪器精密度、标准物质不确定度及制备各种溶液过程中引入的不确定度，这些不确定度来源可从仪器说明书及其它有关文件中得到，通过这些数据及合适的合成方法，可得到测定结果的不确定度。 我们以上面测定钢铁中磷的例子，各种量值不确定度计算如下： 1.标准溶液浓度不确定度u(xi) 标准溶液溶度由下式得出： （10） 其中，xi为标准溶液浓度（单位为ug/ml），ci为标准样品质量比浓度（单位为1），mi为标样质量(单位为ug)，vi为溶液体积(单位为ml)。 由此得出标准溶液浓度不确定度来自三个分量：第一，钢铁标样标值的不确定度uc(xi)；第二，称量过程引入的不确定度um(xi)；第三，定容过程引入的不确定度uv(xi)。由文献[3]公式10得： (11) 式中传播系数： （12） （13） （14） （1）标值的标准不确定度uc(xi) 标准钢铁采用国家标准，其相对扩展不确定度为1%，按正态分布95%置信概率计算，标准不确定度： uc(xi)=0.01/1.96ci =0.0051ci 国标钢样由国家标准研究所提供，其标不准确定度可信度很高，自由度： vc= （2）称量过程引入的不确定度um(xi)  称量不确定度来自两个方面，第一，天平校正产生的不确定度，按检定证书上给出的在95%置信概率时为0.1mg，换算成标准偏差为u1=0.1/1.96=0.052mg自由度为；第二称量变动性，根据历史记录，在50g以内，变动性标准偏差为u2=0.07mg，自由度为。 样品称量过程中，样品质量m=mG-mT得出，其中，mG是样品及容器总质量，mT是容器质量，由此得出样品质量不确定度的四个分量：第一，系统误差导致的毛重不确定度uG1=0.052mg；第二，系统误差导致的皮重不确定度ut1=0.052mg；第三，由随机误差导致的毛重不确定度uG2=0.007mg；第四，由随机误差导致的皮重不确定度uT2=0.007mg。由于毛重和皮重由同一天平称出，因此，系统误差导致的毛重和皮重值不确定度是强相关的，即相关系数为1，而随机误差导致的不确定度是无关的，由此得： 由文献[2]式（G2），自由度为： （3）定容过程引入的不确定度uv(xi)[4] 使用100ml容量瓶配制标准和样品溶液，其不确定度包括三部分：第一，容量瓶体积的不确定度，按制造商给定为0.10ml，按均匀分布换算成标准偏差为： ，自由度为：；第二，充满液体至容量瓶刻度的变动性，可通过重复称量进行统计，例如重复10次统计出标准偏差为0.012ml，自由度：10-1=9；第三，容量瓶和溶液温度与校正时温度不同引起的体积不确定度，假设差3C，水体积膨胀系数为2.1X10-4/C，则95%置信概率时体积变化的区间为100X3X2.1X10-4=0.063ml，转换成标准偏差：0.063/1.96=0.032ml，自由为：。以上三项合成得出： 自由度为： 把5点标准溶液体积、称样量和标称含量代入式（12）、（13）、（14），得各量值的不确定度传播系数。对标准溶液浓度不确定度有贡献的量、其不确定度及传播系数如表2所示。把表2中有关量值代入式（11），得第i个标准溶液不确定度u(xi)，其自由度v(xi)=，列于表3中。 表2 标准溶液浓度不确定度u(xi)计算的有关量值 量 量值 标准不确定度 自由度 传播系数 不确定度分量(ug/ml) c1 0.00011 0.00000056 4755ug/ml 0.00269 c2 0.00022 0.00000112 4864ug/ml 0.005466 c3 0.000269 0.00000137 5167ug/ml 0.00709 c4 0.00038 0.00000194 5053ug/ml 0.00979 c5 0.00063 0.00000321 4524ug/ml 0.0145 m1 475500ug 10ug 0.0000011ml-1 0.000011 m2 486400ug 10ug 0.0000022ml-1 0.000022 m3 516700ug 10ug 0.0000029ml-1 0.000030 m4 505300ug 10ug 0.0000038ml-1 0.000038 m5 452400ug 10ug 0.0000063ml-1 0.000063 v1 100ml 0.067ml -0.00523ug/ml -0.00035 v2 100ml 0.067ml -0.0107ug/ml -0.00072 v3 100ml 0.067ml -0.0150ug/ml -0.00093 v4 100ml 0.067ml -0.0192ug/ml -0.00129 v5 100ml 0.067ml -0.0285ug/ml -0.00191             表3标准溶液浓度不确定度u(xi)及自由度v(xi) 量xi 量值 标准不确定度u(xi) 自由度v(xi) x1 0.523ug/ml 0.00269ug/ml x2 1.07ug/ml 0.00551ug/ml x3 1.39ug/ml 0.00716ug/ml x4 1.92ug/ml 0.00989ug/ml x5 2.58ug/ml 0.0147ug/ml         2.分析谱线强度的标准不确定度 分析谱线强度标准不确定度来源于ICP-AES仪器本身的长时间精密度。以我们采用的美国PE公司OPTIMA 3300 RL为例，其仪器说明书上列明的短时间精密度为10次连续重复测量结果的相对标准偏差不超过1%，长时间精密度是每10分钟重复测量一次相同溶液，4小时内得到数据的相对极差小于10%，即相对标准偏差小于： 我们认为取长时间稳定性做仪器不确定度计算依据较为可信，其标准不确定度： u(yi)=0.0202yi 把5点标准溶液的分析谱线强度代入上式，得其标准不确定度，如表4所示。 自由度： v(yi)=4X60/10-1=23。 相似地，可得样品溶液分析谱线的强度不确定度： u(yt)=0.0202yt, 自由度：v(yt)=23 如表4所示。 表4溶液分析谱线强度不确定度u(yi)和自由度v(yi) 量 量值 标准不确定度u(yi) 自由度v(yi) y1 57.6 cps 1.16cps 23 y2 139.6 cps 2.82cps 23 y3 190.0 cps 3.84cps 23 y4 266.7 cps 5.39cps 23 y5 351.9 cps 7.10cps 23 yt 316.7 cps 6.40cps 23         四、样品溶液浓度合成标准不确定度的计算 计算完上述有关量值，把它们分别代入式（7）、（8）、（9），得出各不确定度分量的传播系数。以上述例，结果如表5所示。 表5 被测样品溶液浓度不确定度u(xt)计算的有关量值 量 量值 标准 不确定度 自由度 传播系数 不确定度分量(ug/ml) x1 0.523ug/ml 0.00269 ug/ml -0.0445 -0.0012 x2 1.07 ug/ml 0.00551 ug/ml 0.0904 0.000498 x3 1.39 ug/ml 0.00716 ug/ml 0.165 0.00118 x4 1.92 ug/ml 0.00989 ug/ml 0.301 0.00297 x5 2.85 ug/ml 0.0147 ug/ml 0.488 0.00717 y1 57.6 cps 1.16cps 23 0.000356ug/ml/cps 0.00414 y2 139.6 cps 2.82cps 23 -0.000630ug/ml/cps -0.00178 y3 190.0 cps 3.84cps 23 -0.00120ug/ml/cps -0.00461 y4 266.7 cps 5.39cps 23 -0.00215ug/ml/cps -0.0116 y5 351.9 cps 7.11cps 23 -0.00334ug/ml/cps -0.0237 yt 316.7 cps 6.33cps 23 0.00696ug/ml/cps 0.0445