[发明专利]一种基于小波变换改进激光探针定量分析的方法无效

专利信息
申请号: 201410183871.1 申请日: 2014-05-04
公开(公告)号: CN103983617A 公开(公告)日: 2014-08-13
发明(设计)人: 李祥友;邹孝恒;曾晓雁;郭连波;郝中骐;沈萌;陆永枫 申请(专利权)人: 华中科技大学
主分类号: G01N21/63 分类号: G01N21/63
代理公司: 华中科技大学专利中心 42201 代理人: 曹葆青
地址: 430074 湖北*** 国省代码: 湖北;42
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 基于 变换 改进 激光 探针 定量分析 方法
【权利要求书】:

1.一种激光探针定量分析,该方法利用激光探针技术得到该待测样品的LIBS光谱,然后利用对应的定标样品得到的优化后的小波函数、分解层数和比例因子γ对光谱进行校正,最后通过单变量线性回归模型预测待测样品中分析元素的含量;所述定标样品是指元素含量已知,且与待测样品具有相同基体元素的样品;

其中,所述优化后的小波函数、分解层数和比例因子γ是按照下述过程得到:

第1步利用激光探针技术对一组定标样品进行检测,得到定标样品的LIBS光谱图;

第2步然后利用小波变换对定标样品的光谱进行背景校正,将背景校正后的特征谱线强度和分析元素浓度分别作为自变量和因变量建立单变量线性回归模型,根据定标均方根误差依次优化小波函数、分解层数和比例因子γ,得到定标样品优化后的小波函数、分解层数和比例因子γ。

2.根据权利要求1所述的激光探针定量分析,其特征在于,第2步具体包括下述过程:

第2.1步对光谱信号进行小波分解:

选用Daubechies小波函数系列,每种小波函数对应的分解层数为4~13,分别对光谱信号进行小波分解,假设小波函数为w,分解层数为l,进行小波分析得到各分解层的高频细节系数D1,D2,……,Dl,以及最高层的低频近似系数C;

第2.2步.第一次重构光谱信号:

将最高层的低频近似系数C置零,进行小波分析,利用修改后的低频分量与其它高频分量重构光谱信号;

第2.3步.分析元素特征谱线的获取:

以选用的标准原子光谱数据库为参考,结合实际光谱,选取无重叠峰干扰和无自吸收现象的分析元素的特征谱线,提取各个定标样品的特征谱线强度I1,I2,……,In,其中,n是定标样品个数;

第2.4步.建立分析元素的单变量线性回归模型:

根据分析元素的浓度矩阵c=[c1,c2,……,cn]T和谱线强度矩阵I=[I1,I2,……,In]T建立单变量线性回归模型,回归方程为:

ci=bIi+a+εi

其中,T为矩阵的转置,b为回归方程的斜率,a为回归方程的截距,ε=[ε12,……,εn]T是残差;

第2.5步.根据第2.1步中选取的小波函数w和分解层数l,γ取定值:γini=1,计算在w和l这两种参数组合下的定标均方根误差RMSEC:

RMSEC=Σi=1nϵi2n]]>

第2.6步.确定使RMSEC取到最小值的小波函数W和分解层数L的组合,即满足:

RMSEC(W,L,γini)=RMSECmin(w,l,γini)

第2.7步.第二次进行小波分解和重构,确定最佳的γ值:

利用第2.6步求得的小波函数W和分解层数L再次对光谱信号进行小波分解。对分解后的最高层低频近似系数C作用一个比例因子γ,即C′=(1-γ)C,利用修改后的低频分量与其它高频分量重构光谱信号,求解γ的最佳值γopt,使RMSEC再次达到最小值,即满足:

RMSEC(W,L,γopt)=RMSECmin(W,L,γ)。

3.根据权利要求1或2所述的激光探针定量分析,其特征在于,所述样品的材料为钢铁。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于华中科技大学,未经华中科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201410183871.1/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top