蒙特卡洛法评定永停滴定法测定磺胺甲氧嘧啶钠注射液含量的不确定度.docxVIP

蒙特卡洛法评定永停滴定法测定磺胺间甲氧嘧啶钠注射液含量的不确定度[摘要]介绍用蒙特卡洛法（MCM）评定永停滴定法测定磺胺间甲氧嘧啶钠注射液含量的不确定度，与引用文献测量不确定度表示指南（GUM）评定结果比较。[关键词]蒙特卡洛法；不确定度；磺胺间甲氧嘧啶钠注射液；含量测定不确定度的评定通常都采用测量不确定度表示指南（GUM）[1]，GUM法在输入量概率分布不对称、测量模型为非线性等情况不适用；蒙特卡洛法（MCM）是一种新的评定不确定度的方法，输入量概率分布不对称、测量模型为非线性等亦适用，适用范围比GUM法更广。MCM不仅可补充GUM法的不足，并且能作为验证GUM法的一种途径[1-2]。本文介绍MCM评定磺胺间甲氧嘧啶钠注射液含量的不确定度方法，为避免重复及比较研究，完全引用杨星等的GUM法数据资料[5]。MCM评定不确定度，需要复杂的计算，采用美国MathWorks公司的Matlab软件（2012b）[3-4]。1、MCM评定不确定度原理简介[2]MCM是实现概率分布传播的一种数值方法，它通过对输入量Xi的概率密度函数（PDF）离散抽样，由测量模型传播输入量的分布，计算获得输出量Y的PDF的离散抽样值，进而由输出量的离散分布数值直接获得输出量的最佳估计值、标准不确定度和包含区间。2、Matlab常用函数说明MCM很关键一点是对所有输入量的PDF进行随机抽样，用数学方法产生随机数，Matlab是一个很强大的数学软件，完全能胜任产生随机数的功能。MCM主要用到正态分布、均匀分布、三角分布等随机函数，Matlab中产生相应随机数的函数及对应语法简介如下：正态分布normrnd(mu,sigma,1,m)产生平均数为mu，标准差为sigma，1行m列的正态分布随机数；均匀分布unifrnd(a,b,1,m)产生1行m列在[a,b]区间的均匀分布随机数；三角分布a+(b-a)/2*(rand(1,m)+rand(1,m))产生1行m列在[a,b]区间的三角分布随机数。2、测量模型及经验参数MCM评定不确定度建立数学模型，与GUM法相同。为了简化过程，NaNO2滴定液的浓度标定MCM评定不确定度省略，仅引用其GUM法评定的不确定度。经化简后测定磺胺间甲氧嘧啶钠注射液含量的数学模型为：式中:V滴为滴定消耗的NaNO2的体积，mL；C为NaNO2滴定液的浓度，mol/L；M为磺胺间甲氧嘧啶钠的摩尔质量，g/moL；V供为移取磺胺间甲氧嘧啶钠的体积，mL；rep为方法重复性因子。据杨星等的文献[5]引用相应参数列于表1。表1 输入量的大小，不确定度和设定的分布输入量描述单位大小标准不确定度或半宽度分布v_NaNO2滴定消耗NaNO2的体积mL15.180.03三角c_NaNO2滴定用NaNO2的浓度moL/L0.10390.00016正态v_SS移取磺胺间甲氧嘧啶钠的体积mL5.000.01三角M_SS磺胺间甲氧嘧啶钠的摩尔质量g/moL302.284810.00589均匀rep重复性因子11.00000.00288正态3、输入量概率密度函数Matlab随机抽样实现3.1、方法的重复性，按正态分布，以因子总量为1，标准不确定度必须用相对不确定度表示，抽样的概率分布为：rep=normrnd(1.0000,0.00288,1,m)3.2、NaNO2的浓度，按正态分布，抽样的概率分布为：c_NaNO2=normrnd(0.10390,0.00016,1,m)3.3、磺胺间甲氧嘧啶钠的摩尔质量，按均匀分布，抽样的概率分布为：M_SS=unifrnd(302.28481-0.00589,302.28481+0.00589,1,m)3.4、滴定消耗NaNO2的体积，按三角分布，抽样的概率分布为：a=15.18-0.03;b=15.18+0.03;v_NaNO2=a+(b-a)/2*(rand(1,m)+rand(1,m))3.5、移取磺胺间甲氧嘧啶钠的体积，按三角分布，抽样的概率分布为：c=5.00-0.01;d=5.00+0.01;v_SS=c+(d-c)/2*(rand(1,m)+rand(1,m))4、MCM评定不确定度Matlab实现据周桃庚的MCM实施的Matlab程序[3]，更换输入量的抽样分布参数及数学模型，略作改动及注释后如下：clear;format long;p=0.95;%概率m=1000000;%设定试验数rep=normrnd(1.0000,0.00288,1,m);%重复性因子c_NaNO2=normrnd(0.10390,0.00016,1,m);%滴定用NaNO2的浓度M_SS=unifrnd(302.28481-0.00589,302.28481+0.00589,1,m);%磺胺间甲氧嘧啶钠的摩尔质量a=15.18-0