CN111189817A - 激光诱导击穿光谱真假信号分辨数据分析方法 - Google Patents

激光诱导击穿光谱真假信号分辨数据分析方法 Download PDF

Info

Publication number
CN111189817A
CN111189817A CN202010022156.5A CN202010022156A CN111189817A CN 111189817 A CN111189817 A CN 111189817A CN 202010022156 A CN202010022156 A CN 202010022156A CN 111189817 A CN111189817 A CN 111189817A
Authority
CN
China
Prior art keywords
signal
sample
measured
laser
true
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010022156.5A
Other languages
English (en)
Other versions
CN111189817B (zh
Inventor
王祥丽
高智星
许小明
柏磊
程毅梅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Atomic of Energy
Original Assignee
China Institute of Atomic of Energy
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute of Atomic of Energy filed Critical China Institute of Atomic of Energy
Priority to CN202010022156.5A priority Critical patent/CN111189817B/zh
Publication of CN111189817A publication Critical patent/CN111189817A/zh
Application granted granted Critical
Publication of CN111189817B publication Critical patent/CN111189817B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/71—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light thermally excited
    • G01N21/718—Laser microanalysis, i.e. with formation of sample plasma
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Immunology (AREA)
  • Algebra (AREA)
  • Pathology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Plasma & Fusion (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明涉及一种激光诱导击穿光谱真假信号分辨数据分析方法，首先分别对本底信号和待测样品信号进行若干次重复测量，得到本底光谱数据和待测样品光谱数据，然后利用t检验法比较待测样品的测定结果的平均值与不存在标准试样的本底测量平均值之间是否存在显著性差异，如果存在显著性差异，则确定测得的信号值为样品信号，反之则为非有效信号。本发明通过将数理统计中的显著性检验与激光诱导击穿光谱技术相结合，显著提高测量数据的可信度，并且算法简便，无需复杂的程序，提高了数据处理效率。

Description

激光诱导击穿光谱真假信号分辨数据分析方法
技术领域
本发明属于光谱测量及分析技术，具体涉及一种激光诱导击穿光谱真假信号分辨数据分析方法。
背景技术
激光诱导击穿光谱(LIBS)是一种应用广泛的元素分析技术，通过利用激光与待测物表面少量样品的相互作用，使样品元素受到激发后发生原子化生成等离子体，在等离子体退激的过程中，待测元素将发射出特征光谱，连接光谱仪后对特征光谱信号进行分析，即可得出待测物中的元素信息。
LIBS因具有操作简便、快速、无需样品预处理、接近无损分析的特征，从上世纪70年代起受到关注，并在多个领域得到应用。进入21世纪，美国发射出的勇气号火星探测器因装备有LIBS检测装置，并成功探测到火星表面下覆盖冰层，说明了LIBS是一项极具潜力且定性探测能力较强的仪器分析技术。
但经过深入研究对比，LIBS在定量分析中具有明显弱势，主要在于检测下限无法满足需求，其中一项重要原因是测量时本底较高，无法分辨出真信号，因而难以进行下一步的数据处理。
根据国内外研究，LIBS的数据分析算法主要依靠测量数据信背比的计算来甄别真假信号，而当本底较强时，测量所得计数因无法满足较高的信背比，数据分析的准确性也随之受到影响。在此情况下，LIBS测量的重复性和稳定性一般用特征信号的相对标准偏差(RSD)来衡量，而当RSD数值较高时，LIBS测量的准确性会受到严重质疑。
在分析工作中，常常会遇到这样一些问题，如对标准试样或纯物质进行测定时，所得到的平均值与标准值的比较问题。不同分析人员、不同实验室和采用不同分析方法对同一试样进行分析时，两组分析结果的平均值之间的比较问题。革新、改造生产工艺后的产品分析指标与原指标的比较问题等。由于测量都有误差存在，数据之间存在的差异是由随机误差引起的，还是由系统误差引起的，若是前者，说明两者无显著性差异，在本质上并无不同，若是后者，说明比较的二者存在显著性差异，即认为它们之间有明显的区别。该类分析方法被称为显著性检验方法，t检验法是其中常用的显著性检验方法。
发明内容
本发明的目的是针对现有的激光诱导击穿光谱数据分析技术在测量本底较高时，难以分辨出真假信号，无法确定数据可信性的缺点，提供一种新的激光诱导击穿光谱真假信号分辨数据分析方法，从而提高测量数据的可信度。
本发明的技术方案如下：一种激光诱导击穿光谱真假信号分辨数据分析方法，首先分别对本底信号和待测样品信号进行若干次重复测量，得到本底光谱数据和待测样品光谱数据，然后利用t检验法比较待测样品的测定结果的平均值与不存在标准试样的本底测量平均值之间是否存在显著性差异，如果存在显著性差异，则确定测得的信号值为样品信号，反之则为非有效信号。
进一步，如上所述的激光诱导击穿光谱真假信号分辨数据分析方法，其中，所述的本底信号和待测样品信号的测量次数不小于20次。
进一步，如上所述的激光诱导击穿光谱真假信号分辨数据分析方法，其中，确定是否存在所述显著性差异的方法可以为，将测量得到的待测样品信号与本底信号进行显著性检验计算，得到t检验值，通过t检验显著性水平对照表查找比较计算值是否大于对应置信区间的t值。
更进一步，查表时所述的对应置信区间为90％。
另外，确定是否存在所述显著性差异的方法还可以为，在假设存在显著性差异的情况下，根据t检验显著性水平对照表中对应置信区间的t值，计算得出可判断为真信号的待测样品测定结果平均值的最低阈值，待测样品的测定结果平均值大于最低阈值时为样品信号。
进一步，如上所述的激光诱导击穿光谱真假信号分辨数据分析方法，其中，t检验法所利用的激光诱导击穿光谱信号数据，是横纵坐标分别为波长与波峰强度的若干个点的数据组，将对应波长下得到的激光诱导击穿光谱信号峰值进行显著性检验，且认为本底测量平均值为显著性检验中的已知样本平均值，待测样品的测量平均值为实验样本平均值。
本发明的有益效果如下：本发明通过将数理统计中的显著性检验与激光诱导击穿光谱技术相结合，提供了一种样品测量信号是否为真信号的判定方法，显著提高测量数据的可信度，并且算法简便，无需复杂的程序，提高了数据处理效率。
附图说明
图1为本发明激光诱导击穿光谱真假信号分辨数据分析方法流程图；
图2为具体实施例中激光诱导击穿光谱测量谱图中真假信号的比较示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明将激光诱导击穿光谱(LIBS)测量结果利用显著性检验方法中的t检验法进行真假信号的分辨，将t检验方法应用到LIBS的数据分析中，从而有效区分真假信号。为了检查分析数据是否存在较大的系统误差，可对待测试样和空白试样进行若干次分析，然后利用t检验法比较待测试样的测定结果的平均值与不存在标准试样的本底测量平均值之间是否存在显著性差异。
进行t检验时，按下式计算出t值
Figure BDA0002361199990000041
其中，
Figure BDA0002361199990000042
为测量所得的样品信号值，μ为空白试样的信号值，s为样品测量值的标准差，n为测量次数。
如图1所示，本发明的分析方法如下：
1)重复测量本底信号20次以上，得到本底光谱数据；
2)重复测量样品信号20次以上，得到相应的待测样品光谱数据；
3)利用t检验法比较待测样品的测定结果的平均值与不存在标准试样的本底测量平均值之间是否存在显著性差异。
可以通过两种方式进行是否存在显著性差异的判断：一是根据上述公式将测量得到的样品信号与本底信号进行显著性检验计算，得到t检验值，通过t检验显著性水平对照表(表1)查表比较计算值是否大于对应置信区间的t值。计算时，当测量次数为20次时，按照已知tα，f表中(见附表1)的数值，设定所测得信号为真信号的可能性(置信区间)为90％，即α＝0.1，此时的t值为1.325。因此，当计算得出的t值大于1.325时，可判定为真信号。二是在假设存在显著性差异的情况下，即已知t为1.325时，利用上述公式计算得出可判断为真信号的样品信号平均值的最低阈值，下文示例中使用该计算方式进行判定。
在该方法中所利用的LIBS信号数据，是横纵坐标分别为波长与波峰强度的若干个点的数据组。将对应波长下得到的LIBS信号峰值进行显著性检验，且认为本底测量平均值为显著性检验中的已知样本平均值，待测样品测量平均值为实验样本平均值。
表1t检验显著性水平对照表
Figure BDA0002361199990000051
实施例
在测量过程中，利用LIBS对空气本底测量200次，计算得出本底的信号强度平均值为3758。对样品实施测量100次，对样品信号值进行计算，得出样品信号标准偏差为1401。在该实施例中，样品测量次数n＝100>20，本底信号的平均值μ＝3758，样品信号的标准差s＝1401，取置信区间为90％，α＝0.1，查表1得t＝1.29。当测量数据平均值满足条件
Figure BDA0002361199990000061
时即表示有显著差异。
根据公式，当
Figure BDA0002361199990000062
时，可确定测得的信号值为样品信号，反之则非有效信号。在该实施例中，样品信号平均值为7591，因此可判断为获得了样品信号，而非本底，因此利用LIBS分析技术可对该浓度样品进行有效测量。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (6)

1.一种激光诱导击穿光谱真假信号分辨数据分析方法，其特征在于：首先分别对本底信号和待测样品信号进行若干次重复测量，得到本底光谱数据和待测样品光谱数据，然后利用t检验法比较待测样品的测定结果的平均值与不存在标准试样的本底测量平均值之间是否存在显著性差异，如果存在显著性差异，则确定测得的信号值为样品信号，反之则为非有效信号。
2.如权利要求1所述的激光诱导击穿光谱真假信号分辨数据分析方法，其特征在于：所述的本底信号和待测样品信号的测量次数不小于20次。
3.如权利要求2所述的激光诱导击穿光谱真假信号分辨数据分析方法，其特征在于：确定是否存在所述显著性差异的方法为，将测量得到的待测样品信号与本底信号进行显著性检验计算，得到t检验值，通过t检验显著性水平对照表查找比较计算值是否大于对应置信区间的t值。
4.如权利要求2所述的激光诱导击穿光谱真假信号分辨数据分析方法，其特征在于：确定是否存在所述显著性差异的方法为，在假设存在显著性差异的情况下，根据t检验显著性水平对照表中对应置信区间的t值，计算得出可判断为真信号的待测样品测定结果平均值的最低阈值，待测样品的测定结果平均值大于最低阈值时为样品信号。
5.如权利要求3或4所述的激光诱导击穿光谱真假信号分辨数据分析方法，其特征在于：查表时所述的对应置信区间为90％。
6.如权利要求1或2所述的激光诱导击穿光谱真假信号分辨数据分析方法，其特征在于：t检验法所利用的激光诱导击穿光谱信号数据，是横纵坐标分别为波长与波峰强度的若干个点的数据组，将对应波长下得到的激光诱导击穿光谱信号峰值进行显著性检验，且认为本底测量平均值为显著性检验中的已知样本平均值，待测样品的测量平均值为实验样本平均值。
CN202010022156.5A 2020-01-09 2020-01-09 激光诱导击穿光谱真假信号分辨数据分析方法 Active CN111189817B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010022156.5A CN111189817B (zh) 2020-01-09 2020-01-09 激光诱导击穿光谱真假信号分辨数据分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010022156.5A CN111189817B (zh) 2020-01-09 2020-01-09 激光诱导击穿光谱真假信号分辨数据分析方法

Publications (2)

Publication Number Publication Date
CN111189817A true CN111189817A (zh) 2020-05-22
CN111189817B CN111189817B (zh) 2021-09-03

Family

ID=70708773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010022156.5A Active CN111189817B (zh) 2020-01-09 2020-01-09 激光诱导击穿光谱真假信号分辨数据分析方法

Country Status (1)

Country Link
CN (1) CN111189817B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104807787A (zh) * 2015-05-20 2015-07-29 重庆邮电大学 一种基于激光诱导击穿光谱技术的茶叶分类鉴别方法
CN105203495A (zh) * 2015-09-11 2015-12-30 天津工业大学 一种基于希尔伯特-黄变换的光谱信号去噪方法
CN105222950A (zh) * 2015-09-18 2016-01-06 中国水利水电科学研究院 一种真假倍频压力脉动信号的区分方法
US20160077110A1 (en) * 2014-09-11 2016-03-17 Board Of Regents Of The University Of Texas System Detection of Misfolded Amyloid Beta Protein
CN106793917A (zh) * 2014-07-02 2017-05-31 新加坡国立大学 用于分析、表征和/或诊断如异常生长的样本或组织的类型或性质的拉曼光谱系统、装置和方法
CN106770189A (zh) * 2016-11-24 2017-05-31 浙江大学 一种基于激光诱导击穿光谱的烟叶铜元素快速检测方法
CN108226086A (zh) * 2018-01-22 2018-06-29 上海海洋大学 一种红葡萄酒酒精度和总糖含量定量分析方法
US20190043703A1 (en) * 2017-08-01 2019-02-07 Protein Metrics Inc. Interactive analysis of mass spectrometry data including peak selection and dynamic labeling
CN109738413A (zh) * 2019-01-08 2019-05-10 江南大学 基于稀疏非负最小二乘的混合物拉曼光谱定性分析方法
CN109900824A (zh) * 2019-03-16 2019-06-18 复旦大学 一种艾纳香中龙脑绝对含量的测定方法
CN110455737A (zh) * 2019-07-11 2019-11-15 中国林业科学研究院林产化学工业研究所 一种利用近红外光谱技术快速分析纸浆材中Klason木质素含量的方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106793917A (zh) * 2014-07-02 2017-05-31 新加坡国立大学 用于分析、表征和/或诊断如异常生长的样本或组织的类型或性质的拉曼光谱系统、装置和方法
US20160077110A1 (en) * 2014-09-11 2016-03-17 Board Of Regents Of The University Of Texas System Detection of Misfolded Amyloid Beta Protein
CN104807787A (zh) * 2015-05-20 2015-07-29 重庆邮电大学 一种基于激光诱导击穿光谱技术的茶叶分类鉴别方法
CN105203495A (zh) * 2015-09-11 2015-12-30 天津工业大学 一种基于希尔伯特-黄变换的光谱信号去噪方法
CN105222950A (zh) * 2015-09-18 2016-01-06 中国水利水电科学研究院 一种真假倍频压力脉动信号的区分方法
CN106770189A (zh) * 2016-11-24 2017-05-31 浙江大学 一种基于激光诱导击穿光谱的烟叶铜元素快速检测方法
US20190043703A1 (en) * 2017-08-01 2019-02-07 Protein Metrics Inc. Interactive analysis of mass spectrometry data including peak selection and dynamic labeling
CN108226086A (zh) * 2018-01-22 2018-06-29 上海海洋大学 一种红葡萄酒酒精度和总糖含量定量分析方法
CN109738413A (zh) * 2019-01-08 2019-05-10 江南大学 基于稀疏非负最小二乘的混合物拉曼光谱定性分析方法
CN109900824A (zh) * 2019-03-16 2019-06-18 复旦大学 一种艾纳香中龙脑绝对含量的测定方法
CN110455737A (zh) * 2019-07-11 2019-11-15 中国林业科学研究院林产化学工业研究所 一种利用近红外光谱技术快速分析纸浆材中Klason木质素含量的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孟宪纾主编: "第五节 假设检验", 《中成药分析》 *

Also Published As

Publication number Publication date
CN111189817B (zh) 2021-09-03

Similar Documents

Publication Publication Date Title
CN104251846B (zh) 一种结合判别分析的激光诱导击穿光谱定量分析方法
CN108398416A (zh) 一种基于激光拉曼光谱的混合物成分测定方法
CN1831516A (zh) 用可见光和近红外光谱技术无损鉴别卷烟品种及真假的方法
CN112345759B (zh) 一种用于荧光强度峰检测的方法
US7768639B1 (en) Methods for detecting and correcting inaccurate results in inductively coupled plasma-atomic emission spectrometry
CN110108658B (zh) 一种污染气体红外光谱识别方法及系统
CN104297206A (zh) 一种液态牛奶品牌快速鉴别方法
CN110231328B (zh) 一种基于半峰高距离法的拉曼光谱定量分析方法
JP2018040779A (ja) 未知試料判定方法、未知試料判定装置及び未知試料判定プログラム
JPH1038807A (ja) プラスチックの判別方法およびプラスチックの判別装置
CN111189817B (zh) 激光诱导击穿光谱真假信号分辨数据分析方法
Gornushkin et al. Insights into linear and rank correlation for material identification in laser-induced breakdown spectroscopy and other spectral techniques
CN110987899B (zh) 一种测定白酒中总酸总酯含量及鉴别白酒质量等级的方法
JP6565801B2 (ja) 質量分析データ処理装置、質量分析装置、質量分析データ処理方法、及び質量分析データ処理用プログラム
CN105717093A (zh) 一种基于大数据库辨识的水泥特性分析方法
CN105717094A (zh) 一种基于大数据库辨识的金属元素含量分析方法
CN110609030B (zh) 一种基于检出概率模型的拉曼快检性能综合评价方法
CN102564965B (zh) 基于光谱检测技术的检测方法
CN114739925A (zh) 基于水拉曼峰强度和波长校准的激光光谱白酒检测方法
CN111044504B (zh) 一种考虑激光诱导击穿光谱不确定性的煤质分析方法
CN113791062A (zh) 一种基于拉曼光谱判断固定物质类别的方法
CN111307753B (zh) 一种基于太赫兹技术的纤维材料检测系统及方法
CN111239237B (zh) 一种高通量质谱检测试剂盒及其质检方法
AU2019100481A4 (en) Quick and nondestructive cubilose authenticity discrimination method
CN112113954A (zh) 一种用于光谱仪的线阵cmos数据处理方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant