CN113567366B - 一种基于光强时间积分的光谱测量衰荡时间的推算方法 - Google Patents
一种基于光强时间积分的光谱测量衰荡时间的推算方法 Download PDFInfo
- Publication number
- CN113567366B CN113567366B CN202110766267.1A CN202110766267A CN113567366B CN 113567366 B CN113567366 B CN 113567366B CN 202110766267 A CN202110766267 A CN 202110766267A CN 113567366 B CN113567366 B CN 113567366B
- Authority
- CN
- China
- Prior art keywords
- ring
- time
- signal
- integral
- test data
- 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.)
- Active
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/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了一种基于光强时间积分的光谱测量衰荡时间的推算方法,包括以下步骤:获取腔衰荡吸收光谱测试中的包含完整衰荡过程的测试数据,生成衰荡曲线;从测试数据中选择多个用于推算光谱测量衰荡时间的分析时间点;从衰荡曲线分别确定多个分析时间点对应的测试数据,利用光强‑时间积分方法分别计算多个分析时间点的事件积分值;利用衰荡时间计算公式处理多个事件积分值,以确定完整衰荡过程的衰荡时间值;本发明选择衰荡过程中的多个衰荡数据段进行积分处理,利用积分值简单地推算衰荡事件的衰荡时间,使用过程简单、计算机运行时间相对较短,且其推算精度不低于传统的最小均方差拟合方法。
Description
技术领域
本发明涉及腔衰荡光谱测试技术领域,具体涉及一种基于光强时间积分的光谱测量衰荡时间的推算方法。
背景技术
在腔衰荡吸收光谱测量技术中,根据衰荡过程产生的衰荡过程数据来推算衰荡时间,进而推测腔内气体的吸收吸收。通常,理想的衰荡过程为单指数衰荡,而在实际的测试中由于光电探测器的特点,衰荡信号会被探测器的零输入时的偏置电压和系统噪声所干扰,衰荡信号偏离单指数衰荡的特点。
在现有技术中,通常利用最小均方拟合来推测衰荡时间,即在腔衰荡光谱测试(CRDS)中,其中的关键数据处理环节是利用实验得到的在衰荡过程中的输出强度随时间的演化推算衰荡时间,进而求得腔内介质的吸收参数,一方面其推算方法或过程直接决定着测量的精度,另一方面,基于测试过程会产生大量的数据和对计算机数据处理的快速要求,多次迭代和矩阵求解的拟合方法中需要较为复杂的数学处理,这样通常需要较长的计算机处理过程和计算时间。
发明内容
本发明的目的在于提供一种基于光强时间积分的光谱测量衰荡时间的推算方法,以解决现有技术中多次迭代和矩阵求解的拟合方法中需要较为复杂的数学处理,这样通常需要较长的计算机处理过程和计算时间的技术问题。
为解决上述技术问题,本发明具体提供下述技术方案:
一种基于光强时间积分的光谱测量的衰荡时间推算方法,包括以下步骤:
步骤100、从腔衰荡吸收光谱测试实验中获取包含完整衰荡过程的测试数据,生成衰荡曲线;
步骤200、从所述衰荡曲线中选择多个用于推算光谱测量衰荡时间的分析时间点;
步骤300、从所述衰荡曲线分别确定多个所述分析时间点对应的测试数据,利用光强-时间积分方法分别计算多个所述分析时间点的测试数据积分值;
步骤400、处理多个所述测试数据积分值得到衰荡时间计算公式,以确定所述完整衰荡过程的衰荡时间值。
作为本发明的优选方案,在步骤100中,所述测试数据在完整衰荡过程中具体为单指数衰荡信号、光电探测器的噪声信号以及偏置信号的叠加数据。
作为本发明的优选方案,在步骤100中,所述衰荡曲线为包含偏置电压和系统噪声信号的单指数衰荡曲线;
所述衰荡曲线的函数关系式为:
其中,I0exp(-t/τRD)单指数衰荡信号;IP为偏置信号;ΔI为噪声信号,且噪音信号的时间积分等于零。
作为本发明的优选方案,在步骤200中,选择分析时间点的筛选规则为:
从所述完整衰荡过程的衰荡前期选择一组前期分析时间点,所述前期分析时间点对应的注入光被完全关断,且所述前期分析时间点的测试数据保持较高的信噪比;
从所述完整衰荡过程的衰荡后期选择一组后期分析时间点,所述后期分析时间点的测试数据包含光电探测器的噪声信号和腔衰荡吸收光谱测试系统的偏置电压信号。
作为本发明的优选方案,所述后期分析时间点对应的所述单指数衰荡信号的强度接近于零,所述后期分析时间点的测量数据只包含偏置信号IP和噪声信号ΔI;
所述前期分析时间点的测量数据包含单指数衰荡信号、偏置信号IP和噪声信号ΔI,任意两个相邻的前期分析时间点之间存在大于零的积分时间段。
作为本发明的优选方案,在步骤300中,利用光强-时间积分方法分别计算多个所述分析时间点的测试数据积分值的实现方法为:
将两个所述前期分析时间点为一组,确定一组前期分析时间点之间的时间差值作为积分时间段,计算每组所述前期分析时间点的测试数据积分值;
计算两个所述后期分析时间点之间的测试数据积分值。
作为本发明的优选方案,所述分析时间点分别为tA、tB、tA+NΔt、tB+NΔt、tC和tD;其中tA、tB、tA+NΔt和tB+NΔt为所述前期分析时间点,tC和tD为后期分析时间点
生成的四个所述测试数据积分值分别为:
其中,处于衰荡后期的所述后期分析时间点的测量数据只包含偏置信号和噪声信号,且噪音信号的时间积分等于零;
PA,N为从tA至tA+NΔt时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tA至tA+NΔt时刻的光强度之和;
PB,N为从tB至tB+NΔt时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tB至tB+NΔt时刻的光强度之和;
PA,B为从tA至tB时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tA至tB时刻的光强度之和;
PC,D为从tC至tD时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tC至tD时刻的光强度之和;
I0为对应测量值的第一个数据;N为正整数;Δt为两个数据点的时间间隔。
作为本发明的优选方案,在步骤400中,所述衰荡时间计算公式中具体为测试数据积分值的等式变换;
基于上述四个所述测试数据积分值的等式,消去偏置信号IP,生成以下等式:
转换等式衰荡时间计算公式,以获得完整衰荡过程的衰荡时间值为:
本发明与现有技术相比较具有如下有益效果:
本发明基于去除系统噪声和偏置电压的基础上,选择衰荡过程中的多个衰荡数据段进行积分处理,利用积分值简单地推算衰荡事件的衰荡时间,使用过程简单、计算机运行时间相对较短,且其推算精度不低于传统的最小均方差拟合方法。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。
图1为本发明实施例提供的衰荡时间推算方法的流程示意图;
图2为本发明实施例提供的衰荡曲线的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供了一种基于光强时间积分的光谱测量的衰荡时间推算方法,本实施方式选择衰荡过程中的多个衰荡数据段进行积分处理,利用积分值简单地推算衰荡事件的衰荡时间,以去除系统噪声和偏置电压,因此推算精度不低于传统的最小均方差拟合方法。
具体包括以下步骤:
步骤100、获取腔衰荡吸收光谱测试实验中的包含完整衰荡过程的测试数据,生成衰荡曲线,所述测试数据在完整衰荡过程中具体为单指数衰荡信号、光电探测器的噪声信号以及偏置信号的叠加数据,即衰荡曲线为包含偏置电压和系统噪声信号的单指数衰荡曲线。
所述衰荡曲线的函数关系式为:
其中,I0exp(-t/τRD)单指数衰荡信号;IP为偏置信号;ΔI为噪声信号,需要特别说明的是,噪音信号的时间积分等于零。
步骤200、从所述测试数据中选择多个用于推算光谱测量衰荡时间的分析时间点。
在步骤200中,选择分析时间点的筛选规则为:
从所述完整衰荡过程的衰荡前期选择一组前期分析时间点,所述前期分析时间点对应的注入光被完全关断,且所述前期分析时间点的测试数据保持较高的信噪比。
从所述完整衰荡过程的衰荡后期选择一组后期分析时间点,所述后期分析时间点的测试数据包含光电探测器的噪声信号和腔衰荡吸收光谱测试系统的偏置电压信号。
所述后期分析时间点对应的所述单指数衰荡信号的强度接近于零,所述后期分析时间点的测量数据只包含偏置信号IP和噪声信号ΔI;
所述前期分析时间点的测量数据包含单指数衰荡信号、偏置信号IP和噪声信号ΔI,任意两个相邻的前期分析时间点之间存在大于零的积分时间段。
步骤300、从所述衰荡曲线分别确定多个所述分析时间点对应的测试数据,利用光强-时间积分方法分别计算多个所述分析时间点的测试数据积分值。
根据步骤100和步骤200可知,其一,后期分析时间点的测试数据中的单指数衰荡信号的强度接近于零,且噪声信号ΔI的积分数据为零,因此后期分析时间点的测试数据积分值仅仅为偏置信号IP的积分值。
其二,当选择的前期分析时间点大于等于三个时,在计算前期分析时间点的测试数据积分值,处于中间的前期分析时间点与前一个前期分析时间点之间的时间差大于零,且小于前一个前期分析时间点与后一个前期分析时间点之间的时间差。
其中,处于衰荡后期的所述后期分析时间点的测量数据只包含偏置信号IP和噪声信号ΔI,且噪音信号的时间积分等于零。
因此利用光强-时间积分方法分别计算多个所述分析时间点的测试数据积分值的实现方法为:
将两个所述前期分析时间点为一组,确定一组前期分析时间点之间的时间差值作为积分时间段,计算每组所述前期分析时间点的测试数据积分值;
计算两个所述后期分析时间点之间的测试数据积分值。
为了具体说明上述衰荡时间推算过程,本实施方式基于腔衰荡吸收光谱测试的实际衰荡曲线,如图2所示,举例说明计算衰减时间的步骤如下:
1、从腔衰荡吸收光谱测试的衰荡事件中得到测试数据,测试数据包含了一个完整的衰荡过程,在衰荡过程的尾部,测试数据只剩下光电探测器的噪声和偏置信号的叠加数据;
2、按照上述分析时间点的筛选规则,选择四个分析时间点,分别为tA、tB、tC和tD;tA和tB选择包含尽量多的有效数据的衰荡过程,因此保证注入光被完全关断以及保持较高的信噪比,tC和tD尽量选择在注入光关断的后期,使这个时间域的测量信号尽量只含有系统的噪声信号成分和探测器的偏置电压信号;
3、在tA和tB的时间段附近选择两个分析时间点tA+NΔt,以及tB+NΔt,其中NΔt的时间长度小于tA和tB的时间跨度,且利用光强-时间积分方法计算tA和tB的测试数据积分值,以及tC和tD的测试数据积分值;
此时所述分析时间点分别为tA、tB、tA+NΔt、tB+NΔt、tC和tD;其中tA、tB、tA+NΔt和tB+NΔt为所述前期分析时间点,tC和tD为后期分析时间点;
生成的四个所述测试数据积分值分别为:
处于衰荡后期的所述后期分析时间点的测量数据只包含偏置信号和噪声信号,且噪音信号的时间积分等于零;
其中,PA,N为从tA至tA+NΔt时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tA至tA+NΔt时刻的光强度之和;
PB,N为从tB至tB+NΔt时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tB至tB+NΔt时刻的光强度之和;
PA,B为从tA至tB时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tA至tB时刻的光强度之和;
PC,D为从tC至tD时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tC至tD时刻的光强度之和;
I0为对应测量值的第一个数据;N为正整数;Δt为两个数据点的时间间隔。
步骤400、处理多个所述测试数据积分值得到衰荡时间计算公式,以确定所述完整衰荡过程的衰荡时间值。
在步骤400中,所述衰荡时间计算公式中具体为测试数据积分值的等式变换;
基于上述四个所述测试数据积分值的等式,消去偏置信号IP,生成以下等式:
转换等式衰荡时间计算公式,以获得完整衰荡过程的衰荡时间值为:
因此,在本实施方式基于得到的实验数据的特点,利用数据的时间平均直接滤掉系统带来的噪声数据的简单推算公式和数据处理方法,即利用衰荡过程多个衰荡数据段的求和积分方法,去除系统噪声和偏置电压,简单地推算衰荡事件的衰荡时间,计算机运行时间相对较短,由于在计算过程中充分考虑到系统噪声和偏置电压,因此衰荡曲线更符合腔衰荡光谱测试的实际数据,并通过积分方式消除系统噪声和偏置电压,因此其推算精度至少不低于传统的最小均方误差拟合方法得到的衰荡时间。
以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。
Claims (5)
1.一种基于光强时间积分的光谱测量的衰荡时间推算方法,其特征在于,包括以下步骤:
步骤100、从腔衰荡吸收光谱测试实验中获取包含完整衰荡过程的测试数据,生成衰荡曲线;
所述衰荡曲线为包含偏置电压和系统噪声信号的单指数衰荡曲线;
所述衰荡曲线的函数关系式为:
其中,I0exp(-t/τRD)单指数衰荡信号;IP为偏置信号;ΔI为噪声信号,且噪音信号的时间积分等于零;
步骤200、从所述衰荡曲线中选择多个用于推算光谱测量衰荡时间的分析时间点;
在步骤200中,选择分析时间点的筛选规则为:
从所述完整衰荡过程的衰荡前期选择一组前期分析时间点,所述前期分析时间点对应的注入光被完全关断,且所述前期分析时间点的测试数据保持高信噪比;
从所述完整衰荡过程的衰荡后期选择一组后期分析时间点,所述后期分析时间点的测试数据包含光电探测器的噪声信号和腔衰荡吸收光谱测试系统的偏置电压信号;
步骤300、从所述衰荡曲线分别确定多个所述分析时间点对应的测试数据,利用光强-时间积分方法分别计算多个所述分析时间点的;
计算多个所述分析时间点的测试数据积分值的实现方法为:
将两个所述前期分析时间点为一组,确定一组前期分析时间点之间的时间差值作为积分时间段,计算每组所述前期分析时间点的测试数据积分值;
计算两个所述后期分析时间点之间的测试数据积分值;
步骤400、处理多个所述测试数据积分值得到衰荡时间计算公式,以确定所述完整衰荡过程的衰荡时间值,所述衰荡时间计算公式具体为测试数据积分值的等式变换。
2.根据权利要求1所述的一种基于光强时间积分的光谱测量衰荡时间的推算方法,其特征在于,
在步骤100中,所述测试数据在完整衰荡过程中具体为单指数衰荡信号、光电探测器的噪声信号以及偏置信号的叠加数据。
3.根据权利要求1所述的一种基于光强时间积分的光谱测量衰荡时间的推算方法,其特征在于,
所述后期分析时间点对应的所述单指数衰荡信号的强度接近于零,所述后期分析时间点的测量数据只包含偏置信号IP和噪声信号ΔI;
所述前期分析时间点的测量数据包含单指数衰荡信号、偏置信号IP和噪声信号ΔI,任意两个相邻的前期分析时间点之间存在大于零的积分时间段。
4.根据权利要求1所述的一种基于光强时间积分的光谱测量衰荡时间的推算方法,其特征在于,
所述分析时间点分别为tA、tB、tA+NΔt、tB+NΔt、tC和tD;其中tA、tB、tA+NΔt和tB+NΔt为所述前期分析时间点,tC和tD为后期分析时间点;
生成的四个所述测试数据积分值分别为:
其中,处于衰荡后期的所述后期分析时间点的测量数据只包含偏置信号和噪声信号,且噪音信号的时间积分等于零;
PA,N为从tA至tA+NΔt时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tA至tA+NΔt时刻的光强度之和;
PB,N为从tB至tB+NΔt时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tB至tB+NΔt时刻的光强度之和;
PA,B为从tA至tB时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tA至tB时刻的光强度之和;
PC,D为从tC至tD时刻实验测量,包括偏置信号和噪声信号的强度积分值,即在离散状态下,从tC至tD时刻的光强度之和;
I0为对应测量值的第一个数据;N为正整数;Δt为两个数据点的时间间隔。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110766267.1A CN113567366B (zh) | 2021-07-07 | 2021-07-07 | 一种基于光强时间积分的光谱测量衰荡时间的推算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110766267.1A CN113567366B (zh) | 2021-07-07 | 2021-07-07 | 一种基于光强时间积分的光谱测量衰荡时间的推算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113567366A CN113567366A (zh) | 2021-10-29 |
CN113567366B true CN113567366B (zh) | 2022-09-30 |
Family
ID=78163889
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110766267.1A Active CN113567366B (zh) | 2021-07-07 | 2021-07-07 | 一种基于光强时间积分的光谱测量衰荡时间的推算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113567366B (zh) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050012931A1 (en) * | 2003-07-17 | 2005-01-20 | Sze Tan | Method and apparatus for rapidly and accurately determining a time constant from cavity ring-down data |
CN101261182B (zh) * | 2008-03-26 | 2010-06-23 | 中国科学院光电技术研究所 | 基于频率选择性光反馈光腔衰荡技术的高反射率测量方法 |
CN101644673A (zh) * | 2009-09-09 | 2010-02-10 | 中国科学院光电技术研究所 | 基于量子级联激光器的红外光腔衰荡光谱痕量气体检测方法 |
CN102768198A (zh) * | 2012-07-05 | 2012-11-07 | 中国科学技术大学 | 一种利用锁频激光的光腔衰荡光谱技术测量气体组分含量的系统和方法 |
CN104458586B (zh) * | 2014-12-05 | 2017-02-22 | 中国科学院光电技术研究所 | 一种在光反馈光腔衰荡技术中筛选衰荡信号的方法 |
CN110672554B (zh) * | 2019-11-11 | 2020-09-25 | 北京航空航天大学 | 一种随机振动驱动衰荡腔免标定气体浓度测量系统 |
CN111122496B (zh) * | 2019-12-20 | 2021-03-23 | 清华大学 | 一种免标定的气体浓度测量装置及方法 |
-
2021
- 2021-07-07 CN CN202110766267.1A patent/CN113567366B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113567366A (zh) | 2021-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112881959B (zh) | 一种用于磁共振成像的梯度涡流补偿方法及系统 | |
CN103475369B (zh) | 基于信号源误差一次性校准识别的高精度adc测试方法 | |
CN111487213A (zh) | 多光谱融合的化学需氧量测试方法及装置 | |
Tseitlin et al. | Uncertainty analysis for absorption and first‐derivative electron paramagnetic resonance spectra | |
Ables et al. | A 1024− channel digital correlator | |
CN110768661B (zh) | 一种基于神经网络的锁相放大器 | |
CN113567366B (zh) | 一种基于光强时间积分的光谱测量衰荡时间的推算方法 | |
CN111337746A (zh) | 一种射频信号检测定位方法 | |
Max | Testing high speed high accuracy analog to digital converters embedded in systems on a chip | |
US20050230632A1 (en) | Spectra acquisition system with threshold adaptation integrator | |
Mitzkus et al. | Surface brightness fluctuation spectrum: a new probe of evolved stars in unresolved stellar populations | |
Geraci et al. | Shared baseline restoration at minimum noise for high resolution spectroscopy | |
Werle | Time domain characterization of micrometeorological data based on a two sample variance | |
Bittanti et al. | High-accuracy fit of the poles of spectroscopy amplifiers designed for mixed analog-digital filtering | |
CN110865057A (zh) | 一种应用于荧光寿命成像的非均匀时间数字转换器 | |
Capofreddi et al. | The use of linear models for the efficient and accurate testing of A/D converters | |
Carbone et al. | Statistical efficiency of the ADC sinewave histogram test | |
Van den Bogaert et al. | Evaluation and correction of signal model errors in a matched filter for the quantification of chromatographic data | |
Ajith et al. | Null-stream veto for two co-located detectors: implementation issues | |
Woehl et al. | User-independent nonlinear modeling using adjusted spline-interpolated knots (UNMASK) and indirect hard modeling for deriving compositions from spectra with background signals | |
Diversi | A three-step identification procedure for ARARX models with additive measurement noise | |
US20060176979A1 (en) | Wave detection device, method, program, and recording medium | |
US7010439B2 (en) | Digital circuit for measuring the power of a signal | |
Gardner | 6. Uncertainty Estimates in Radiometry | |
RU2104495C1 (ru) | Способ измерения физической величины |
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 |