CN105181624B - 一种基于散射类比的太赫兹光谱定量分析方法 - Google Patents
一种基于散射类比的太赫兹光谱定量分析方法 Download PDFInfo
- Publication number
- CN105181624B CN105181624B CN201510557368.2A CN201510557368A CN105181624B CN 105181624 B CN105181624 B CN 105181624B CN 201510557368 A CN201510557368 A CN 201510557368A CN 105181624 B CN105181624 B CN 105181624B
- Authority
- CN
- China
- Prior art keywords
- terahertz
- sample
- mixture
- scattering
- composition
- 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.)
- Expired - Fee Related
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种基于散射类比的太赫兹光谱定量分析方法,属于太赫兹光谱应用技术领域。本发明对三个不同浓度的纯聚乙烯(PE)压片样品进行太赫兹时域光谱检测,获取其太赫兹吸收谱后进行平均计算,用平均计算后的PE太赫兹吸收谱来表征散射效应;分别获取混合物样品及混合物中的各成分单质的太赫兹吸收谱;实验测得的样品的太赫兹吸收谱实际由样品本身对太赫兹波的吸收以及散射效应造成的衰减两部分组成,据此建立实验测得的混合物样品的太赫兹吸收谱模型;依据类比散射的原理,建立混合物样品的太赫兹吸收谱与混合物中各成分的太赫兹吸收谱及PE太赫兹吸收谱之间的关系;利用最小二乘算法,计算得到混合物中各成分的浓度。
Description
技术领域
本发明涉及一种基于散射类比的太赫兹光谱定量分析方法,属于太赫兹光谱应用技术领域。
背景技术
近年来随着超快激光技术的发展,为太赫兹波的产生和探测提供了一系列有效的方式。太赫兹(THz,1THz=1012Hz)波作为电磁波谱中的一个波段,长久以来由于缺乏有效的产生源和探测装置,一直没有得到深度的开发利用,以至于在很长一段时间内都被称为“太赫兹间隙”。由于太赫兹波的独特特性,使其在多个领域展示出了巨大的潜力,如物质识别、医药、安检以及天文观测等。其中,基于太赫兹吸收谱的定量分析是太赫兹技术最为重要的应用方向之一。
目前,应用于太赫兹定量分析领域的方法主要有最小二乘、偏最小二乘等。早在2007年,北京科技大学的张朝晖教授等人就通过分析半胱氨酸、天冬酰胺和苏氨酸等三种L型氨基酸及其混合物的太赫兹吸收谱,利用最小二乘法对混合物中各成分的浓度进行了计算,取得了良好的效果。桂林电子科技大学的陈涛等人研究了基于偏最小二乘的特征谱区筛选法在太赫兹光谱定量分析中的应用。实验表明,采用移动窗口偏最小二乘法的交互验证均方根误差可达0.9803,证明该方法能够有效的对多元混合物的成分含量进行快速检测。
但是现有技术在太赫兹定量分析领域中均忽略了散射效应的影响。尤其是当样品的颗粒尺寸与入射波的波长相接近的时候,会发生严重的散射现象。此时测得的吸收谱由两部分组成,即样品本身对入射光的吸收以及由散射效应导致的衰减。以往的分析方法仅考虑了样品本身的吸收情况,而忽略了散射效应的作用,使得定量分析的精确度不高。
发明内容
本发明的目的是提供一种基于散射类比的太赫兹光谱定量分析方法,以解决目前太赫兹光谱定量分析中仅考虑了样品本身的吸收情况,而忽略了散射效应的作用,从而导致定量分析的结果准确度不高的技术问题。
本发明为解决上述技术问题而提供一种基于散射类比的太赫兹光谱定量分析方法,该方法包括以下步骤:
1)对三个不同浓度的纯聚乙烯(PE)压片样品进行太赫兹时域光谱检测,获取其太赫兹吸收谱并进行平均计算,并用平均计算后的PE太赫兹吸收谱来表征散射效应;
2)分别对混合物样品及混合物中的各成分单质进行太赫兹时域光谱检测,获取其太赫兹吸收谱;
3)实验测得的样品的太赫兹吸收谱包含两部分,即样品本身对太赫兹波的吸收以及由于散射效应造成的衰减,据此建立实验测得的混合物样品的太赫兹吸收谱模型;
4)依据类比散射的原理,对所述步骤3)中的混合物样品的太赫兹吸收谱模型进行变换,得到变换后的混合物样品的太赫兹吸收谱模型,从而建立混合物样品的太赫兹吸收谱与混合物中各成分的太赫兹吸收谱及PE太赫兹吸收谱之间的关系;
5)利用最小二乘算法,计算得到混合物中各成分的浓度。
所述步骤3)中的混合物样品的太赫兹吸收谱模型为:
(1)
其中ω i 是频率,α mix_e 是通过实验得到的混合物样品的太赫兹吸收系数,α mix_s 是混合物样品的散射系数,下标1,2,…,n表示混合物中的各种成分,下标e表示通过实验所得,下标s表示样品的散射系数,即α 1_e 是通过实验得到的第一种成分的太赫兹吸收系数,α 2_e 是通过实验得到的第二种成分的太赫兹吸收系数,α 1_s 是第一种成分的散射系数,α 2_s 是第二种成分的散射系数,以此类推。c 1 是第一种样品在混合物中的浓度,c 2 是第二种样品在混合物中的浓度,以此类推。
所述步骤4)中的类比散射的原理为:
(2)
其中α PE 是平均后的PE的吸收系数,k 1 ,k 2 ,…,k n ,k mix 是混合物中各成分及混合物与散射相关的无量纲的系数。
所述步骤4)中的变换后的混合物样品的太赫兹吸收谱模型为:
(3)
(4)
其中K是一个与散射相关的无量纲的系数。
所述步骤5)中的最小二乘算法为:
令,,,
则公式(3)可表示为:
(5)
其中MIX为混合物样品的吸收系数矩阵,A为混合物中各成分单质及PE的吸收系数矩阵,C为包含散射系数K的混合物中各成分的浓度矩阵。则有:
(6)
本发明的有益效果是:本发明在对混合物样品的太赫兹波段散射效应进行分析的基础上,采用类比的方法,提出了混合物的太赫兹吸收谱模型。在利用最小二乘法对混合物进行定量分析时,考虑了散射效应的影响,从而提高了定量分析的准确度。取得了良好的应用效果。同时通过定量分析实验进行了验证,本发明能够实现对混合物浓度的精确定量分析。
附图说明
图1. 本发明算法流程图
图2. 混合物样品(谷氨酰胺,组氨酸)的扫描电镜图
图3. 纯PE压片的太赫兹吸收谱(虚线)及其平均值(实线)
图4. 谷氨酰胺、组氨酸以及苏氨酸的摩尔吸收系数谱
图5. 谷氨酰胺、组氨酸混合样品(虚线)和组氨酸、苏氨酸混合样品(实线)的太赫兹吸收谱。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明。本发明提供一种基于散射类比的太赫兹光谱定量分析方法,在对太赫兹光谱定量分析时,考虑散射相应造成的影响。通过一系列实验证明了在对混合物进行定量分析时将散射效应考虑进来,能够提高定量分析的准确度。
当样品的颗粒尺寸与入射波的波长相接近的时候,会发生严重的散射现象。此时测得的样品吸收谱由两部分组成,即样品本身对入射光的吸收以及由散射效应导致的衰减:
其中α e 是样品的太赫兹吸收系数,α in 是样品本身的吸收,α s 是散射效应造成的衰减,也称为散射系数。
根据朗博比尔定律的描述,物质对光的吸收与其浓度成正比例关系,由多种成分组成的混合物样品的吸收系数就可以用下面的公式表示:
(1)
其中ω i 是频率,α mix_e 是通过实验得到的混合物样品的太赫兹吸收系数,α mix_s 是混合物样品的散射系数,下标1,2,…,n表示混合物中的各种成分,下标e表示通过实验所得,下标in表示样品本身对太赫兹波的吸收,下标s表示样品的散射系数,即α 1_e 是通过实验得到的第一种成分的太赫兹吸收系数,α 2_e 是通过实验得到的第二种成分的太赫兹吸收系数,α 1_in 是第一种成分本身对太赫兹波的吸收系数,α 2_in 是第二种成分本身对太赫兹波的吸收系数,α 1_s 是第一种成分的散射系数,α 2_s 是第二种成分的散射系数,以此类推。c 1 是第一种样品在混合物中的浓度,c 2 是第二种样品在混合物中的浓度,以此类推。公式(1)也可称为混合物的太赫兹吸收谱模型。
所有的样品都由聚乙烯(PE)进行稀释后压成压片然后检测。由于PE对太赫兹波没有吸收,因此可以认为实验测得的PE的吸收谱完全是由散射效应造成的。由于样品颗粒与PE颗粒的尺寸相接近,而散射效应主要与颗粒的尺寸有关,因此我们测试了三个由纯PE组成的压片,并将得到的太赫兹吸收谱进行了平均计算。然后采用类比的方法用平均计算后的PE的吸收谱来表征样品在太赫兹波段的散射效应,即将样品的散射效应表示为带有系数的PE吸收谱的一次函数:
(2)
其中αPE是平均后的PE的吸收系数,k1,k2,…,kn,kmix是混合物中各成分及混合物与散射相关的无量纲的系数。因此,公式(1)就可以写成下列形式:
(3)
(4)
其中K是一个与散射相关的无量纲的系数。利用公式(3),并结合最小二乘算法,就可以对成分已知的混合物进行定量分析,并取得准确结果。
令,,,
则公式(3)可表示为:
(5)
其中MIX为混合物样品的吸收系数矩阵,A为混合物中各成分单质及PE的吸收系数矩阵,C为包含散射系数K的混合物中各成分的浓度矩阵。则有:
(6)
为了验证本方法的优越性,设计了一系列定量分析的实验。样品选取了3种氨基酸,分别是谷氨酰胺,组氨酸和苏氨酸,并将其进行了两两混合,然后用指标qe来表征定量分析的误差。qe的定义如下所示:
其中qe表征定量分析的误差,n是混合物样品中不同成分种类的个数,c i_real 和c i_cal 分别是各组分真实的以及计算得到的浓度,均由百分比表示。
样品的组成以及定量分析的结果如下表所示。
表1. 样品的组成以及定量分析的误差
上述实验数据表明,利用该算法能够实现对混合物浓度的精确定量分析,误差大致在3%以下,取得了优异的效果。
Claims (1)
1.一种基于散射类比的太赫兹光谱定量分析方法,其特征在于,该分析方法包括以下步骤:
1)对三个不同浓度的纯聚乙烯PE压片样品进行太赫兹时域光谱检测,获取其太赫兹吸收谱并进行平均计算,并用平均计算后的PE太赫兹吸收谱来表征散射效应;
2)分别对混合物样品及混合物中的各成分单质进行太赫兹时域光谱检测,获取其太赫兹吸收谱;
3)实验测得的样品的太赫兹吸收谱包含两部分,即样品本身对太赫兹波的吸收以及由于散射效应造成的衰减,据此建立实验测得的混合物样品的太赫兹吸收谱模型;
4)依据散射类比的原理,对所述步骤3)中的混合物样品的太赫兹吸收谱模型进行变换,得到变换后的混合物样品的太赫兹吸收谱模型,从而建立混合物样品的太赫兹吸收谱与混合物中各成分的太赫兹吸收谱及PE太赫兹吸收谱之间的关系;
5)利用最小二乘算法,计算得到混合物中各成分的浓度;
所述步骤3)中的混合物样品的太赫兹吸收谱模型为:
αmix_e(ωi)=c1·α1_e(ωi)+c2·α2_e(ωi)+…+cn·αn_e(ωi)-(c1·α1_s(ωi)+c2·α2_s(ωi)+…+cn·αn_s(ωi)-αmix_s(ωi)) (1)
其中ωi是频率,αmix_e是通过实验得到的混合物样品的太赫兹吸收系数,αmix_s是混合物样品的散射系数,下标1,2,...,n表示混合物中的各种成分,下标e表示通过实验所得,下标s表示样品的散射系数,即α1_e是通过实验得到的第一种成分的太赫兹吸收系数,α2_e是通过实验得到的第二种成分的太赫兹吸收系数,α1_s是第一种成分的散射系数,α2_s是第二种成分的散射系数,以此类推;c1是第一种样品在混合物中的浓度,c2是第二种样品在混合物中的浓度,以此类推;
所述步骤4)中的散射类比的原理为:
α1_s(ωi)=k1·αPE(ωi)
α2_s(ωi)=k2·αPE(ωi)
……
αn_s(ωi)=kn·αPE(ωi)
αmix_s(ωi)=kmix·αPE(ωi) (2)
其中αPE是平均后的PE的吸收系数,k1,k2,...,kn,kmix是混合物中各成分及混合物与散射相关的无量纲的系数;
混合物样品的太赫兹吸收谱模型为:
αmix_e(ωi)=c1·α1_e(ωi)+c2·α2_e(ωi)+…+cn·αn_e(ωi)
-(c1·k1·αPE(ωi)+c2·k2·αPE(ωi)+…+cn·kn·αPE(ωi)-kmix·αPE(ωi))
=c1·α1_e(ωi)+c2·α2_e(ωi)+…+cn·αn_e(ωi)-K·αPE(ωi) (3)
K=c1·k1+c2·k2+…+cn·kn-kmix (4)
其中K是一个与散射相关的无量纲的系数;
所述步骤5)中的最小二乘算法为:
令
则公式(3)可表示为:
MIX=A·C (5)
其中MIX为混合物样品的吸收系数矩阵,A为混合物中各成分单质及PE的吸收系数矩阵,C为包含散射系数K的混合物中各成分的浓度矩阵;则有:
C=(A′·A)-1·A′·MIX (6)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510557368.2A CN105181624B (zh) | 2015-09-06 | 2015-09-06 | 一种基于散射类比的太赫兹光谱定量分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510557368.2A CN105181624B (zh) | 2015-09-06 | 2015-09-06 | 一种基于散射类比的太赫兹光谱定量分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105181624A CN105181624A (zh) | 2015-12-23 |
CN105181624B true CN105181624B (zh) | 2017-11-07 |
Family
ID=54903857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510557368.2A Expired - Fee Related CN105181624B (zh) | 2015-09-06 | 2015-09-06 | 一种基于散射类比的太赫兹光谱定量分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105181624B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572076B (zh) * | 2016-01-13 | 2018-11-30 | 南京邮电大学 | 基于散射效应的太赫兹波谱测量装置及其测量方法 |
CN106295614A (zh) * | 2016-08-23 | 2017-01-04 | 上海理工大学 | 太赫兹波检测颗粒样品频谱散射线的去除算法 |
CN114965843A (zh) * | 2022-05-11 | 2022-08-30 | 浙江树人学院 | 一种多通道检测的离子色谱仪及其使用方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104730025A (zh) * | 2015-04-01 | 2015-06-24 | 河南工业大学 | 一种基于太赫兹光谱的混合物定量分析方法 |
CN104849233A (zh) * | 2015-05-13 | 2015-08-19 | 河南工业大学 | 一种检测谷物新陈度的方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009236565A (ja) * | 2008-03-26 | 2009-10-15 | Jasco Corp | 複数成分の自動連続定量分析方法およびその装置 |
ES2875957T3 (es) * | 2012-12-20 | 2021-11-11 | Amgen Inc | Agonistas del receptor APJ y usos de los mismos |
-
2015
- 2015-09-06 CN CN201510557368.2A patent/CN105181624B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104730025A (zh) * | 2015-04-01 | 2015-06-24 | 河南工业大学 | 一种基于太赫兹光谱的混合物定量分析方法 |
CN104849233A (zh) * | 2015-05-13 | 2015-08-19 | 河南工业大学 | 一种检测谷物新陈度的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105181624A (zh) | 2015-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105181624B (zh) | 一种基于散射类比的太赫兹光谱定量分析方法 | |
Kahraman et al. | The effect of rock classes on the relation between uniaxial compressive strength and point load index | |
WO2018168939A8 (ja) | 蛍光x線分析方法、蛍光x線分析プログラムおよび蛍光x線分析装置 | |
WO2011027613A1 (ja) | 蛍光x線分析方法 | |
CN110208241B (zh) | 基于受激拉曼散射的大气单颗粒物快速三维化学成像方法 | |
NZ703875A (en) | Spectroscopic analysis | |
EA201270260A1 (ru) | Способ и устройство анализа материалов с помощью сфокусированного пучка электронов с использованием характеристического рентгеновского излучения и обратноотраженных электронов | |
CN105651801B (zh) | 一种矿浆矿物在线分析方法 | |
CN104914089A (zh) | 用表面增强拉曼光谱对痕量混合物进行半定量分析的方法 | |
CN102213682A (zh) | 一种干涉不敏感太赫兹波透射测量方法 | |
King et al. | A micro-reflectance IR spectroscopy method for analyzing volatile species in basaltic, andesitic, phonolitic, and rhyolitic glasses | |
CN105092436A (zh) | 一种沉积物粒度光谱分析方法及装置 | |
Sitko et al. | Fast and simple method for determination of fatty acid methyl esters (FAME) in biodiesel blends using X-ray spectrometry | |
CN102410998A (zh) | 激光光散射-拉曼光谱联用装置 | |
CN105004641A (zh) | 一种血液分析仪及分析方法 | |
CN104076054B (zh) | 一种基于x射线荧光快速检测烟草中重金属元素铅的方法 | |
JP2019158560A5 (zh) | ||
CN102854177B (zh) | 用拉曼光谱法对待测样品中高铁血红蛋白含量MetHb%进行测定的方法 | |
Xu et al. | Quantitative mineral analysis by FTIR spectroscopy | |
Bertucci et al. | A new model for x‐ray fluorescence autoabsorption analysis of pigment layers | |
CN103226100A (zh) | 一种量子产率的测试方法 | |
Powell et al. | Evaluation of uncertainties in X‐ray photoelectron spectroscopy intensities associated with different methods and procedures for background subtraction. I. Spectra for monochromatic Al X‐ray | |
CN106092989B (zh) | 一种通过计数量子点团聚比例定量肝素的方法 | |
CN106769996A (zh) | 一种小麦粉中山梨酸钾检测方法 | |
CN105699357B (zh) | 一种基于二维拉曼光谱的颗粒粒度测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171107 Termination date: 20180906 |