法律状态公告日
法律状态信息
法律状态
2019-10-25
专利权质押合同登记的生效 IPC(主分类):G01N27/48 登记号:Y2019430000024 登记生效日:20190926 出质人:长沙华时捷环保科技发展股份有限公司 质权人:长沙银行股份有限公司科技支行 发明名称:线性扫描极谱曲线重叠峰分离方法 授权公告日:20140709 申请日:20120616
专利权质押合同登记的生效、变更及注销
2017-05-03
专利权的转移 IPC(主分类):G01N27/48 登记生效日:20170410 变更前: 变更后: 申请日:20120616
专利申请权、专利权的转移
2014-07-09
授权
授权
2012-12-19
实质审查的生效 IPC(主分类):G01N27/48 申请日:20120616
实质审查的生效
2012-10-24
公开
公开
技术领域] 本发明涉及一种线性扫描极谱曲线重叠峰分离方法,特别是离子之间呈高浓度比时的重叠峰分离。
背景技术]
线性扫描极谱法是金属离子浓度在线检测中应用最广泛的一种测定方法,它分别通过测定极谱波的半波电位和波高实现金属离子浓度的定性和定量分析,当使用该方法同时测定半波电位相近且呈高浓度比的几种金属离子时,由于测试体系不可能完全消除高含量基体的影响,高含量的重金属离子的极谱峰极易掩蔽低含量的重金属离子极谱峰;同时,极谱特性相近的两种重金属组分之间也存在极谱峰重叠的问题,严重影响测定方法的测量精度和分辨力。而对产生重叠峰的金属离子进行分析时,通常采用实验室离线分析的方法,滞后时间长,无法为生产过程提供实时信息,影响工业生产的优化控制运行。
重叠峰分离方法的研究和应用普遍存在于可以用标准Gauss函数等单一函数描述的色谱或光谱分离中,然而,线性扫描极谱法生成的极谱曲线难以用单一函数进行描述,因此,这些方法无法应用于此类重叠峰的分离,且目前尚无有效方法对线性扫描极谱波重叠峰实现分离。这一难题已成为制约线性扫描极谱法在离子浓度在线检测领域应用的主要障碍之一。
发明内容]
针对极谱特性相近的金属离子,其线性扫描极谱曲线容易产生重叠峰,而影响其极谱分析检测的问题,本发明的目的是提供一种新型类Gauss函数,以描述金属离子非对称性的线性扫描极谱曲线,然后,通过基于小波变换的非线性加权最小二乘法对重叠信号曲线进行拟合,最后,将各金属离子的极谱图进行分离重构,进而通过计算可以得到各离子的极谱电流,从而实现线性扫描极谱波重叠信号的分离。
所述的类Guass函数构造方法是将标准Gauss函数、反正切函数和线性函数叠加形成,用来精确近似描述线性扫描极谱曲线,利用该函数的参数与峰高、半峰宽以及半波电位等极谱曲线的特征参数之间的对应关系,将计算得到的该函数的参数值用于求解线性扫描极谱曲线的各特征参数值,实现线性扫描极谱曲线的精确定量近似描述。
单个线性扫描极谱检测波谱曲线如图1所示,其具有如下特点:
a. 峰形非对称性;
b. 前峰可以近似为标准Gauss曲线;
c. 除去标准Gauss曲线的部分有一个迅速上升之后趋于恒定值的过程,可以近似为一个S形曲线,因此可以选用反正切函数进行描述。
将两个基本函数叠加之后可以近似描述单个线性扫描极谱检测波谱曲线,如图2 所示。
单个线性扫描极谱检测波谱曲线可以表示为:
式中,k1、k2、…、k6为该极谱波的重构参数;x表示扫描电压,i残余表示对单个空白溶液施加扫描电压时生成的电流。
则多个极谱曲线可以表示为:
式中,n表示波谱信号有n个离子产生极谱波构成,ki1、ki2、…、ki6为第i个极谱波的重构参数; I残余表示对n种空白溶液施加扫描电压时生成的电流的和。
假设残余电流I残余为:
式中,k、b分别表示残余电流I残余的重构参数。
由公式(2)、(3)可得零阶导数波谱信号的表达式:
对公式(4)求导,可得到一阶导数波谱信号δ(x):
对公式(5)求导,可得到二阶导数波谱信号γ(x):
对原始极谱曲线求一阶和二阶导数,可以提高检测曲线的分辨率。为提高拟合精度和拟合曲线的物理意义,需要对极谱波形进一步分析,确定重叠在一起的极谱波峰的原始位置。小波变换的多尺度分析可以将不同频率特性的波形分开,而极谱重叠波在波峰或者波谷位置重叠时会导致该点的频率发生改变。因此,可以利用小波变换的多尺度分解准确找出重叠波原始波形的波峰和波谷位置。
为保证小波变换的多尺度分解可以将重叠峰完全分离,须选择最佳的离散细节D(j) ,根据Nyquist临界采样频率可知:
可得,最佳离散细节D(j)的频率要大于或等于信号频率f0,从而,有2jfc≈f0,这里将极谱波信号拓展为周期信号,可认为其频率满足:
为表示重叠峰的重叠程度,引入色谱分析领域分离度R的概念,其表达式如式(9)所示:
式中,pA和pB分别表示重叠峰中前后两个峰的位置,wA和wB分别表示两个峰的半峰宽度,对于两个重叠峰来说:
由公式(7)~(10)可以得到:
对于同一信号而言,最佳分解次数只与分离度有关,而与小波基的选取无关。
由此可以得到各阶波谱曲线组成重叠信号的单元素极谱电位
对于公式(5)存在:
对于公式(6)存在:
公式(13)~(14)中
由于杂波和残余电流的影响,在极谱波峰前后的波谱曲线存在较多的干扰因素,会影响曲线拟合的质量,本发明的重点在于极谱波形,为提高曲线拟合的物理意义,采用非线性加权最小二乘法(Weighted Least Squares, WLS)拟合。则拟合函数的偏差平方和可以表示为:
考虑到本发明拟合曲线专注于信号重叠峰的分离,极谱波与高斯曲线类似,因此,设计权函数表达式为:
式(15)中,
通过非线性WLS以及公式(12)~(14)约束条件,对分别对零阶、一阶以及二阶极谱曲线进行参数辨识,即可得到公式(4)~(6)中的参数。最后,对利用拟合得到的类Gauss函数的参数重构单个离子的极谱曲线,实现线性扫描极谱波重叠峰的分离。
本发明通过小波变换能够精确的确定各极谱峰的位置,对曲线拟合起到了很好的导向作用,半波电位的相对误差基本维持在1%以下,基于小波变换的新型类Gauss拟合重叠峰分离方法,起到了很好的分离效果。
附图说明]
图1 线性扫描极谱法电流-电位曲线;
图2 构造类Gauss函数曲线;
图3 原始线性扫描极谱波重叠峰;
图4 一阶导数线性扫描极谱波重叠峰;
图5 二阶导数线性扫描极谱波重叠峰;
图6 Co2+和Zn2+重叠峰WLS拟合的权值分布曲线;
图7 基于小波变换的类Gauss拟合原始极谱曲线;
图8 基于小波变换的类Gauss拟合一阶极谱曲线;
图9 基于小波变换的类Gauss拟合二阶极谱曲线。
具体实施方式]
以极谱电位相差约50mV,浓度比1:105的Co2+和Zn2+为检测实例,通过线性扫描极谱法获取钴锌的标准极谱曲线如图3,对原始数据进行一阶、二阶求导后得到一阶极谱曲线如图4和二阶极谱曲线如图5,针对产生的重叠峰具体分离方法实施如下:
1、确定各个组分极谱波波峰和波谷的位置
依据公式(11)可以得到图3~图5的三种极谱曲线的最佳分解层数为8,对三种极谱波数据通过小波变换进行压缩,然后分别按最佳分解次数原则进行小波分解,确定钴离子极谱波的位置。
可以从D(8)层重构的波形得到钴波的峰值位电位为
表2 小波变换判定各阶极谱波峰、谷位置
2、确定拟合函数
由公式(4)~(6)得到钴锌重叠极谱零阶导数谱峰信号的表达式:
一阶导数谱峰信号δ(x)可以表示为:
同时,待定参数满足如下条件:
公式(11)中,
同理,二阶导数谱峰信号γ(x)可以表示为:
待定参数满足如下条件:
3、确定非线性WLS拟合权函数和拟合参数
钴锌极谱检测为两个离子极谱波重叠,由此
式(15)中,
由图6中,可以看出在Co2+波和Zn2+波位置权值较高,拟合曲线将尽可能逼近实际测量值,而在两个极谱波之外权值较小,其在WLS拟合寻优过程中所占的比重较小,由非线性WLS拟合的零阶导数极谱波拟合曲线如图7所示。
为进一步提高运算速度和拟合质量,本发明引入含约束条件的非线性WLS用于导数极谱波重叠峰的曲线拟合。一阶导数极谱波拟合函数参数的约束条件可以表示为:
式中
二阶导数极谱波拟合函数参数的约束条件可以表示为公式(23),式中
由于线性扫描极谱法的波形不对称特性,因此,本发明提出类高斯函数以实现该类重叠峰的分离,由于引入了不对称函数,因此,相应的几个波峰和波谷的位置参考点在拟合前,并不能准确描绘出各个待定参数的具体关系,只能为拟合提供初始值,以提高拟合速度。
表3 拟合数据的相似度和均方差统计表
由图7至图9的拟合结果以及表3的评价指标可以看出,零阶导数的拟合效果最好,二阶次之,一阶最差。那么,影响拟合结果出现偏差的因素有很多种,可能性最大的有:其它痕量杂质离子的干扰,噪声的干扰,残余电流的影响等。为进一步验证拟合结果的好坏,进行了单组分对比分析实验,分别测量了与钴、锌混合溶液相同浓度的
表5-5 拟合结果与实验对比
由表3和表4可以看出一个问题,半波电位的相对误差基本维持在1%以下,说明小波变换能够精确的确定各极谱峰的位置,对曲线拟合起到了很好的导向作用。另外,虽然从图中的拟合效果看,随着阶次的升高拟合效果逐渐变差,然而,实验对比的结果却显示随着阶次的升高极谱峰电流的相对误差逐渐减小,一方面说明导数极谱波可以提高检测的灵敏度,另一方面说明,重叠峰检测信号可能由于杂质等因素影响了极谱波形,然而,曲线拟合的过程反而起到了修复波形的作用。综上所述,基于小波变换的Gauss拟合重叠峰分离方法,起到了很好的分离效果。
机译: 用于记录极谱仪电压曲线的装置
机译: 用komp enso图记录极谱电流-电压-曲线的装置
机译: 分离光谱中重叠峰的装置和使用相同方法的X射线荧光分析装置