CN113899703A - 一种波长偏移下的线型展宽消除方法 - Google Patents

一种波长偏移下的线型展宽消除方法 Download PDF

Info

Publication number
CN113899703A
CN113899703A CN202111150439.9A CN202111150439A CN113899703A CN 113899703 A CN113899703 A CN 113899703A CN 202111150439 A CN202111150439 A CN 202111150439A CN 113899703 A CN113899703 A CN 113899703A
Authority
CN
China
Prior art keywords
signal
standard
harmonic
formula
calibration
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
CN202111150439.9A
Other languages
English (en)
Other versions
CN113899703B (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.)
Anhui Agricultural University AHAU
Original Assignee
Anhui Agricultural University AHAU
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 Anhui Agricultural University AHAU filed Critical Anhui Agricultural University AHAU
Priority to CN202111150439.9A priority Critical patent/CN113899703B/zh
Publication of CN113899703A publication Critical patent/CN113899703A/zh
Application granted granted Critical
Publication of CN113899703B publication Critical patent/CN113899703B/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/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/255—Details, e.g. use of specially adapted sources, lighting or optical systems
    • 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

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)
  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种波长偏移下的线型展宽消除方法，是基于一次谐波信号进行处理，包括：1去除背景噪声后再进行滤波处理，2标准一次谐波拟合，通过引入置信度并结合标准差对拟合的无吸收区域作为标准一次谐波自适应拟合是否结束的评判标准，得到标准一次谐波拟合信号，为线型展宽的消除提供了基础，3两次信号调整，先进行粗调，并对局部信号拟合，插值进行微调，消除展宽。本发明能在不增加硬件结构的基础上，通过一种波长偏移下的线型展宽消除方法，从而解决波长偏移下的线型展宽问题，并提高测量的准确性。

Description

一种波长偏移下的线型展宽消除方法
技术领域
本发明涉及环境光学领域，具体的说是针对波长偏移所导致线型展宽影响，提出一种消除波长偏移所导致线型展宽的方法。
背景技术
在利用激光光谱技术进行气体检测时，温度变化，电流扫描信号直流电平漂移，驱动电流与温度的漂移，热敏电阻、激光器芯片的非完全耦合，光路、电路老化等会造成激光器波长偏移，导致线型展宽，影响气体浓度测量的准确性。因此有效的消除线型展宽问题，提高测量的长期稳定性与准确性有重要研究意义。
现有的研究主要从改造硬件系统例如控温、锁相等减少漂移问题，很少讨论如何消除波长偏移导致的线型展宽的影响。B.Tuzson等人(Atmospheric MeasurementTechniques,2013,6:927)指出抑制波长漂移引起的谱线展宽是保证QCL系统长期稳定的基础。陈等人(光谱学与光谱分析,2018)设计了一种高精度温度控制系统，以消除背景信号漂移的影响，增强测量系统的稳定性。现有的技术虽然在一定程度上解决波长偏移的问题，但在实际的测量过程中，引起线型展宽复杂多样，仅通过单一的改变硬件结构的角度，并不能从根本上解决。
发明内容
鉴于现有技术方法存在的不足，本发明提出一种波长偏移下的线型展宽消除方法，以期能在不增加硬件结构的基础上，解决波长偏移下的线型展宽问题，并提高测量的准确性。
本发明解决上述问题的技术方案如下：
本发明一种波长偏移下的线型展宽消除方法，是应用于由激光器、控制器、函数发生器、校准通路、探测通路、采集处理器依次连接所构成的线型展宽消除气体检测系统中，其特点是，所述线型展宽消除方法是按如下步骤进行：
步骤1.所述采集处理器获取当前周期下所述校准通路的信号序列{C(m)|m＝1,2,…,N}与所述探测通路的信号序列{Y(m)|m＝1,2,…,N}，其中，C(m)为所述校准通路在当前周期下第m位置上的校准电压信号；Y(m)为所述探测通路在当前周期下第m位置上的探测电压信号；
步骤2.分别对信号序列{C(m)|m＝1,2,…,N}和{Y(m)|m＝1,2,…,N}去除背景噪声后再进行滤波处理，得到滤波后的信号序列{C′(m)|m＝1,2,…,N}和{Y′(m)|m＝1,2,…,N}；其中，C′(m)为所述校准通路在当前周期下第m位置上滤波后的校准电压信号；Y′(m)为所述探测通路在当前周期下第m位置上滤波后的探测电压信号；
步骤3.选用洛伦兹展宽主导线型，分别对滤波后的信号序列{C′i(m)|m＝1,2,…,N}和{Y′i(m)|m＝1,2,…,N}进行标准谐波拟合后相应得到标准一次谐波信号序列{XC1f(m)|m＝1,2,…,N}与{XY1f(m)|m＝1,2,…,N}；其中，XC1f(m)为所述校准通路在当前周期下第m位置上标准一次谐波拟合后的校准信号；XY1f(m)为所述探测通路在当前周期下第m位置上标准一次谐波拟合后的探测信号；
步骤4.对标准一次谐波信号序列{XY1f(m)|m＝1,2,…,N}求取极大值位置的信号XY1f(m0)与极小值位置的信号XY1f(m2)；
对标准一次谐波信号序列{XC1f(m)|m＝1,2,…,N}求取极大值位置的信号XC1f(m1)与极小值位置的信号XC1f(m3)；
利用式(1)得到极大值点的偏移量Δmax与极小值点的偏移量Δmin：
Figure BDA0003286852650000021
式(1)中，m0和m2分别表示标准一次谐波信号序列{XY1f(m)|m＝1,2,…,N}的极大值和极小值的位置；m1和m3分别表示标准一次谐波信号序列{XC1f(m)|m＝1,2,…,N}的极大值和极小值的位置；
利用式(2)所得到的偏移量Δ，对标准一次谐波信号序列{XY1f(m)|m＝1,2,…,N}进行粗调，从而得到粗调后的标准一次谐波信号序列{XY1f(m′)|m′＝1,2,…,N′}，其中，XY1f(m′)表示对第m位置上探测信号XY1f(m)进行粗调后的信号，且m′＝m±Δ：
Figure BDA0003286852650000022
步骤5.以极大值的位置m0为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(3)对粗调后的信号XY1f(m′)进行二次低阶拟合，得到局部峰值信号
Figure BDA0003286852650000023
Figure BDA0003286852650000024
式(3)中，b0、b1、b2表示三个拟合系数，k表示步进值，且m0-3Δ≤k≤m0+3Δ；
以极大值的位置m1为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(4)对校准信号XC1f(m)进行二次低阶拟合，得到局部校准峰值信号
Figure BDA0003286852650000031
Figure BDA0003286852650000032
式(4)中，a0、a1、a2表示三个拟合系数，k表示步进值，且m1-3Δ≤k≤m1+3Δ；
利用插值法对局部峰值信号
Figure BDA0003286852650000033
求取极大值位置的信号
Figure BDA0003286852650000034
利用插值法对局部校准峰值信号
Figure BDA0003286852650000035
求取极大值位置的信号
Figure BDA0003286852650000036
利用式(5)得到极大值点微调偏移量ΔL：
Figure BDA0003286852650000037
式(5)中，
Figure BDA0003286852650000038
表示信号
Figure BDA0003286852650000039
的极大值位置，
Figure BDA00032868526500000310
表示信号
Figure BDA00032868526500000311
的极大值位置；
步骤6.以极小值的位置m2为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(6)对粗调后的信号XY1f(m′)进行二次低阶拟合，得到局部谷值信号
Figure BDA00032868526500000312
Figure BDA00032868526500000313
式(6)中，c0、c1、c2表示三个拟合系数，k表示步进值，且m2-3Δ≤k≤m2+3Δ；
以极小值的位置m3为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(7)对校准信号XC1f(m)进行二次低阶拟合，得到局部校准谷值信号
Figure BDA00032868526500000314
Figure BDA00032868526500000315
式(7)中，a0、a1、a2表示三个拟合系数，k表示步进值，且m3-3Δ≤k≤m3+3Δ；
利用插值法对局部谷值信号
Figure BDA00032868526500000316
求取极大值位置的信号
Figure BDA00032868526500000317
利用插值法对局部校准谷值信号
Figure BDA00032868526500000318
求取极大值位置的信号
Figure BDA00032868526500000319
利用式(8)得到极小值点微调偏移量ΔM：
Figure BDA00032868526500000320
式(8)中，
Figure BDA0003286852650000041
表示信号
Figure BDA0003286852650000042
的极大值位置，
Figure BDA0003286852650000043
表示信号
Figure BDA0003286852650000044
的极大值位置；
根据所述的极大值点微调偏移量ΔL与极小值点微调偏移量ΔM计算均值求得微调偏移量ΔA，用于对粗调后的信号XY1f(m′)进行微调，得到展宽消除后的信号XY1f(m″)；
步骤7.对展宽消除后的信号XY1f(m″)与信号XC1f(m1)进行反演，并得到所述探测通路中的气体浓度。
本发明所述的一种波长偏移下的线型展宽消除方法的特点也在于，所述步骤3中的标准谐波拟合是判断滤波后的信号序列中的每个电压信号到零线的距离是否处于所设置的置信度阈值范围内，若是，在保持相应的电压信号，否则，相应的电压信号设为零，从而得到相应的一次谐波信号的临时序列；
再利用最小二乘拟合法对临时序列进行处理，并在信号序列及其对应的临时序列之间的标准差之和最小时，得到最终的标准一次谐波信号序列。
与现有的技术相比较，本发明的有益效果体现在：
1.本发明中在不增加硬件结构的基础上，通过对一次谐波信号标准拟合，先进行粗调，并对局部信号拟合，插值进行微调，通过两次信号调整，从而消除了展宽，提高了测量准确性，解决了信号的波长偏移引起的线型展宽，影响测量准确性的问题。
2.本发明中通过引入置信度结合标准差对拟合的无吸收区域作为标准一次谐波自适应拟合是否结束的评判标准，得到标准一次谐波拟合信号，为线型展宽的消除提供了基础。
附图说明
图1为本发明实施案例中的线型展宽消除气体检测系统的示意图；
图2本发明的一种波长偏移下的线型展宽消除方法的流程图；
图3为本发明的展宽消除前后图；
图4为本发明实施案例中的浓度修正前后图。
具体实施方式
本实施例中，参见图1，一种线型展宽消除气体检测系统，包括：激光器3、控制器2、函数发生器1、校准通路、探测通路、采集处理器。其中，校准通路由准直器4、标准参考池5、光电探测器6组成；探测通路由准直器7、多次反射池8、光电探测器9组成；
激光器1作为检测光源进行测量，通过控制器2进行温度控制和电流控制，从而改变激光器1的输出波长；函数发生器1向控制器2的电流驱动器提供扫描信号和调制信号，共同作用于激光器3,调谐激光稳定输出；激光分路后分别经过校准通路和探测通路，其中校准通路的激光光束经准直器5准直后，通过标准参考池5后，由光电探测器6接收；探测通路的激光光束经准直器7后，进入多次反射池8后，汇聚到光电探测器9；两种信号都由采集处理器完成信号采集与处理。
参见图2，本实施例中，一种基于线型展宽消除气体检测系统波长偏移下的线型展宽消除方法是按如下步骤进行：
步骤1.采集处理器获取当前周期下校准通路的信号序列{C(m)|m＝1,2,…,N}与探测通路的信号序列{Y(m)|m＝1,2,…,N}，其中，C(m)为校准通路在当前周期下第m位置上的校准电压信号；Y(m)为探测通路在当前周期下第m位置上的探测电压信号；
步骤2.分别对信号序列{C(m)|m＝1,2,…,N}和{Y(m)|m＝1,2,…,N}去除背景噪声后再进行滤波处理，得到滤波后的信号序列{C′(m)|m＝1,2,…,N}和{Y′(m)|m＝1,2,…,N}；其中，C′(m)为校准通路在当前周期下第m位置上滤波后的校准电压信号；Y′(m)为探测通路在当前周期下第m位置上滤波后的探测电压信号；
步骤3.选用洛伦兹展宽主导线型，分别对滤波后的信号序列{C′i(m)|m＝1,2,…,N}和{Y′i(m)|m＝1,2,…,N}进行标准谐波拟合后相应得到标准一次谐波信号序列{XC1f(m)|m＝1,2,…,N}与{XY1f(m)|m＝1,2,…,N}；其中，XC1f(m)为校准通路在当前周期下第m位置上标准一次谐波拟合后的校准信号；XY1f(m)为探测通路在当前周期下第m位置上标准一次谐波拟合后的探测信号；
具体实施中，标准谐波拟合是判断滤波后的信号序列中的每个电压信号到零线的距离是否处于所设置的置信度阈值范围内，若是，在保持相应的电压信号，否则，相应的电压信号设为零，从而得到相应的一次谐波信号的临时序列；
再利用最小二乘拟合法对临时序列进行处理，并在信号序列及其对应的临时序列之间的标准差之和最小时，得到最终的标准一次谐波信号序列。
步骤4.对标准一次谐波信号序列{XY1f(m)|m＝1,2,…,N}求取极大值位置的信号XY1f(m0)与极小值位置的信号XY1f(m2)；
对标准一次谐波信号序列{XC1f(m)|m＝1,2,…,N}求取极大值位置的信号XC1f(m1)与极小值位置的信号XC1f(m3)；
利用式(1)得到极大值点的偏移量Δmax与极小值点的偏移量Δmin：
Figure BDA0003286852650000061
式(1)中，m0和m2分别表示标准一次谐波信号序列{XY1f(m)|m＝1,2,…,N}的极大值和极小值的位置；m1和m3分别表示标准一次谐波信号序列{XC1f(m)|m＝1,2,…,N}的极大值和极小值的位置；
利用式(2)所得到的偏移量Δ，对标准一次谐波信号序列{XY1f(m)|m＝1,2,…,N}进行粗调，从而得到粗调后的标准一次谐波信号序列{XY1f(m′)|m′＝1,2,…,N′}，其中，XY1f(m′)表示对第m位置上探测信号XY1f(m)进行粗调后的信号，且m′＝m±Δ：
Figure BDA0003286852650000062
步骤5.以极大值的位置m0为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(3)对粗调后的信号XY1f(m′)进行二次低阶拟合，得到局部峰值信号
Figure BDA0003286852650000063
Figure BDA0003286852650000064
式(3)中，b0、b1、b2表示三个拟合系数，k表示步进值，且m0-3Δ≤k≤m0+3Δ；
以极大值的位置m1为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(4)对校准信号XC1f(m)进行二次低阶拟合，得到局部校准峰值信号
Figure BDA0003286852650000065
Figure BDA0003286852650000066
式(4)中，a0、a1、a2表示三个拟合系数，k表示步进值，且m1-3Δ≤k≤m1+3Δ；
利用插值法对局部峰值信号
Figure BDA0003286852650000067
求取极大值位置的信号
Figure BDA0003286852650000068
利用插值法对局部校准峰值信号
Figure BDA0003286852650000069
求取极大值位置的信号
Figure BDA00032868526500000610
利用式(5)得到极大值点微调偏移量ΔL：
Figure BDA00032868526500000611
式(5)中，
Figure BDA0003286852650000071
表示信号
Figure BDA0003286852650000072
的极大值位置，
Figure BDA0003286852650000073
表示信号
Figure BDA0003286852650000074
的极大值位置；
步骤6.以极小值的位置m2为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(6)对粗调后的信号XY1f(m′)进行二次低阶拟合，得到局部谷值信号
Figure BDA0003286852650000075
Figure BDA0003286852650000076
式(6)中，c0、c1、c2表示三个拟合系数，k表示步进值，且m2-3Δ≤k≤m2+3Δ；
以极小值的位置m3为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(7)对校准信号XC1f(m)进行二次低阶拟合，得到局部校准谷值信号
Figure BDA0003286852650000077
Figure BDA0003286852650000078
式(7)中，a0、a1、a2表示三个拟合系数，k表示步进值，且m3-3Δ≤k≤m3+3Δ；
利用插值法对局部谷值信号
Figure BDA0003286852650000079
求取极大值位置的信号
Figure BDA00032868526500000710
利用插值法对局部校准谷值信号
Figure BDA00032868526500000711
求取极大值位置的信号
Figure BDA00032868526500000712
利用式(8)得到极小值点微调偏移量ΔM：
Figure BDA00032868526500000713
式(8)中，
Figure BDA00032868526500000714
表示信号
Figure BDA00032868526500000715
的极大值位置，
Figure BDA00032868526500000716
表示信号
Figure BDA00032868526500000717
的极大值位置；
根据的极大值点微调偏移量ΔL与极小值点微调偏移量ΔM计算均值求得微调偏移量ΔA，用于对粗调后的信号XY1f(m′)进行微调，得到展宽消除后的信号XY1f(m″)；
步骤7.对展宽消除后的信号XY1f(m″)与信号XC1f(m1)进行反演，并得到探测通路中的气体浓度。
为了验证本发明方法的效果，本例利用所建立的线型展宽消除气体检测系统进行实验，本发明的设置扫描频率为100Hz，调制频率为50kHz，采样速度为200kHz的采集卡，探测通路路采用20m的多次反射池，在反射池内密封浓度约为40ppm的甲烷进行测量。利用本发明方法修正前后结果图如图3所示，修正后波长偏移下的线型展宽得到了抑制。进一步验证本发明方法对浓度反演计算中的准确性影响，对固定浓度30ppm的甲烷进行测量。测量结果如图4示，展宽修正后标准差为0.047，未考虑展宽影响的标准差为0.203，实验结果证明了本发明提供的方法能在波长发生展宽时有效提高气体浓度测量的准确性。

Claims (2)

1.一种波长偏移下的线型展宽消除方法，是应用于由激光器、控制器、函数发生器、校准通路、探测通路、采集处理器依次连接所构成的线型展宽消除气体检测系统中，其特征是，所述线型展宽消除方法是按如下步骤进行：
步骤1.所述采集处理器获取当前周期下所述校准通路的信号序列{C(m)|m＝1，2，…，N}与所述探测通路的信号序列{Y(m)|m＝1，2，…，N}，其中，C(m)为所述校准通路在当前周期下第m位置上的校准电压信号；Y(m)为所述探测通路在当前周期下第m位置上的探测电压信号；
步骤2.分别对信号序列{C(m)|m＝1，2，…，N}和{Y(m)|m＝1，2，…，N}去除背景噪声后再进行滤波处理，得到滤波后的信号序列{C′(m)|m＝1，2，…，N}和{Y′(m)|m＝1，2，…，N}；其中，C′(m)为所述校准通路在当前周期下第m位置上滤波后的校准电压信号；Y′(m)为所述探测通路在当前周期下第m位置上滤波后的探测电压信号；
步骤3.选用洛伦兹展宽主导线型，分别对滤波后的信号序列{C′i(m)|m＝1，2，…，N}和{Yi′(m)|m＝1，2，…，N}进行标准谐波拟合后相应得到标准一次谐波信号序列{XC1f(m)|m＝1，2，…，N}与{XY1f(m)|m＝1，2，…，N}；其中，XC1f(m)为所述校准通路在当前周期下第m位置上标准一次谐波拟合后的校准信号；XY1f(m)为所述探测通路在当前周期下第m位置上标准一次谐波拟合后的探测信号；
步骤4.对标准一次谐波信号序列{XY1f(m)|m＝1，2，…，N}求取极大值位置的信号XY1f(m0)与极小值位置的信号XY1f(m2)；
对标准一次谐波信号序列{XC1f(m)|m＝1，2，…，N}求取极大值位置的信号XC1f(m1)与极小值位置的信号XC1f(m3)；
利用式(1)得到极大值点的偏移量Δmax与极小值点的偏移量Δmin：
Figure FDA0003286852640000011
式(1)中，m0和m2分别表示标准一次谐波信号序列{XY1f(m)|m＝1，2，…，N}的极大值和极小值的位置；m1和m3分别表示标准一次谐波信号序列{XC1f(m)|m＝1，2，…，N}的极大值和极小值的位置；
利用式(2)所得到的偏移量Δ，对标准一次谐波信号序列{XY1f(m)|m＝1，2，…，N}进行粗调，从而得到粗调后的标准一次谐波信号序列{XY1f(m′)|m′＝1，2，…，N′}，其中，XY1f(m′)表示对第m位置上探测信号XY1f(m)进行粗调后的信号，且m′＝m±Δ：
Figure FDA0003286852640000021
步骤5.以极大值的位置m0为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(3)对粗调后的信号XY1f(m′)进行二次低阶拟合，得到局部峰值信号
Figure FDA0003286852640000022
Figure FDA0003286852640000023
式(3)中，b0、b1、b2表示三个拟合系数，k表示步进值，且m0-3Δ≤k≤m0+3Δ；
以极大值的位置m1为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(4)对校准信号XC1f(m)进行二次低阶拟合，得到局部校准峰值信号
Figure FDA0003286852640000024
Figure FDA0003286852640000025
式(4)中，a0、a1、a2表示三个拟合系数，k表示步进值，且m1-3Δ≤k≤m1+3Δ；
利用插值法对局部峰值信号
Figure FDA0003286852640000026
求取极大值位置的信号
Figure FDA0003286852640000027
利用插值法对局部校准峰值信号
Figure FDA0003286852640000028
求取极大值位置的信号
Figure FDA0003286852640000029
利用式(5)得到极大值点微调偏移量ΔL：
Figure FDA00032868526400000210
式(5)中，
Figure FDA00032868526400000211
表示信号
Figure FDA00032868526400000212
的极大值位置，
Figure FDA00032868526400000213
表示信号
Figure FDA00032868526400000214
的极大值位置；
步骤6.以极小值的位置m2为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(6)对粗调后的信号XY1f(m′)进行二次低阶拟合，得到局部谷值信号
Figure FDA00032868526400000215
Figure FDA00032868526400000216
式(6)中，c0、c1、c2表示三个拟合系数，k表示步进值，且m2-3Δ≤k≤m2+3Δ；
以极小值的位置m3为中心，取±3Δ范围内粗调后的信号作为待拟合区，利用式(7)对校准信号XC1f(m)进行二次低阶拟合，得到局部校准谷值信号
Figure FDA0003286852640000031
Figure FDA0003286852640000032
式(7)中，a0、a1、a2表示三个拟合系数，k表示步进值，且m3-3Δ≤k≤m3+3Δ；
利用插值法对局部谷值信号
Figure FDA0003286852640000033
求取极大值位置的信号
Figure FDA0003286852640000034
利用插值法对局部校准谷值信号
Figure FDA0003286852640000035
求取极大值位置的信号
Figure FDA0003286852640000036
利用式(8)得到极小值点微调偏移量ΔM：
Figure FDA0003286852640000037
式(8)中，
Figure FDA0003286852640000038
表示信号
Figure FDA0003286852640000039
的极大值位置，
Figure FDA00032868526400000310
表示信号
Figure FDA00032868526400000311
的极大值位置；
根据所述的极大值点微调偏移量ΔL与极小值点微调偏移量ΔM计算均值求得微调偏移量ΔA，用于对粗调后的信号XY1f(m′)进行微调，得到展宽消除后的信号XY1f(m″)；
步骤7.对展宽消除后的信号XY1f(m″)与信号XC1f(m1)进行反演，并得到所述探测通路中的气体浓度。
2.根据权利要求1所述的一种波长偏移下的线型展宽消除方法，其特征是，所述步骤3中的标准谐波拟合是判断滤波后的信号序列中的每个电压信号到零线的距离是否处于所设置的置信度阈值范围内，若是，在保持相应的电压信号，否则，相应的电压信号设为零，从而得到相应的一次谐波信号的临时序列；
再利用最小二乘拟合法对临时序列进行处理，并在信号序列及其对应的临时序列之间的标准差之和最小时，得到最终的标准一次谐波信号序列。
CN202111150439.9A 2021-09-29 2021-09-29 一种波长偏移下的线型展宽消除方法 Active CN113899703B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111150439.9A CN113899703B (zh) 2021-09-29 2021-09-29 一种波长偏移下的线型展宽消除方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111150439.9A CN113899703B (zh) 2021-09-29 2021-09-29 一种波长偏移下的线型展宽消除方法

Publications (2)

Publication Number Publication Date
CN113899703A true CN113899703A (zh) 2022-01-07
CN113899703B CN113899703B (zh) 2022-07-05

Family

ID=79189298

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111150439.9A Active CN113899703B (zh) 2021-09-29 2021-09-29 一种波长偏移下的线型展宽消除方法

Country Status (1)

Country Link
CN (1) CN113899703B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114970639A (zh) * 2022-06-10 2022-08-30 安徽农业大学 一种开放空间外界环境的气体闪烁噪声消除方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7969576B1 (en) * 2007-03-23 2011-06-28 The Regents Of The University Of California Optical sensing based on wavelength modulation spectroscopy
US20130222801A1 (en) * 2010-08-16 2013-08-29 The University Of Chicago Real-Time Mapping of Electronic Structure with Single-Shot Two-Dimensional Electronic Spectroscopy
US20130311110A1 (en) * 2012-05-18 2013-11-21 Konstantin AIZIKOV Methods and Apparatus for Obtaining Enhanced Mass Spectrometric Data
CN104126106A (zh) * 2011-12-20 2014-10-29 小利兰·斯坦福大学托管委员会 一种用于气体感测的免校准扫描波长调制光谱的方法
CN108061722A (zh) * 2017-12-05 2018-05-22 中国科学院合肥物质科学研究院 一种一氧化碳浓度的检测装置及检测方法
US20180275049A1 (en) * 2015-09-30 2018-09-27 Consiglio Nazionale Delle Ricerche - Cnr Method for measuring the concentration of trace gases by scar spectroscopy
CN109270027A (zh) * 2018-09-19 2019-01-25 清华大学 一种基于正弦调制时域拟合的气体吸收率在线测量方法
CN111829980A (zh) * 2020-07-23 2020-10-27 安徽农业大学 一种基于谐波技术的线型非线性修正的检测系统及方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7969576B1 (en) * 2007-03-23 2011-06-28 The Regents Of The University Of California Optical sensing based on wavelength modulation spectroscopy
US20130222801A1 (en) * 2010-08-16 2013-08-29 The University Of Chicago Real-Time Mapping of Electronic Structure with Single-Shot Two-Dimensional Electronic Spectroscopy
CN104126106A (zh) * 2011-12-20 2014-10-29 小利兰·斯坦福大学托管委员会 一种用于气体感测的免校准扫描波长调制光谱的方法
US20130311110A1 (en) * 2012-05-18 2013-11-21 Konstantin AIZIKOV Methods and Apparatus for Obtaining Enhanced Mass Spectrometric Data
US20180275049A1 (en) * 2015-09-30 2018-09-27 Consiglio Nazionale Delle Ricerche - Cnr Method for measuring the concentration of trace gases by scar spectroscopy
CN108061722A (zh) * 2017-12-05 2018-05-22 中国科学院合肥物质科学研究院 一种一氧化碳浓度的检测装置及检测方法
CN109270027A (zh) * 2018-09-19 2019-01-25 清华大学 一种基于正弦调制时域拟合的气体吸收率在线测量方法
CN111829980A (zh) * 2020-07-23 2020-10-27 安徽农业大学 一种基于谐波技术的线型非线性修正的检测系统及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIANG, QIJUN 等: "《Temperature Dependence of Free Carriers and Lattice Vibrations in Novel Nonlinear Optical Crystals in the THz Frequency Regime》", 《2017 CONFERENCE ON LASERS AND ELECTRO-OPTICS EUROPE & EUROPEAN QUANTUM ELECTRONICS CONFERENCE (CLEO/EUROPE-EQEC)》 *
刘立拓 等: "《激光诱导击穿光谱数据特征自动提取方法研究》", 《光谱学与光谱分析》 *
尹松林: "《被动调Q激光诱导击穿光谱探测研究》", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅰ辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114970639A (zh) * 2022-06-10 2022-08-30 安徽农业大学 一种开放空间外界环境的气体闪烁噪声消除方法
CN114970639B (zh) * 2022-06-10 2024-04-16 安徽农业大学 一种开放空间外界环境的气体闪烁噪声消除方法

Also Published As

Publication number Publication date
CN113899703B (zh) 2022-07-05

Similar Documents

Publication Publication Date Title
CN113899703B (zh) 一种波长偏移下的线型展宽消除方法
CN107462551B (zh) 一种用于气体检测的激光波长锁定的方法及装置
CN111398216A (zh) 气体浓度检测的多频调制方法、气体浓度检测方法及系统
JPH10510355A (ja) 情報取得方法
CN110160976B (zh) 基于二次谐波干涉包络补偿的气体浓度检测方法和系统
CN111537470A (zh) 一种基于数字调制的tdlas气体浓度检测方法
CN103868884A (zh) 一种基于调制系数的气体吸收率在线测量方法
CN111462707B (zh) 一种基于建模补偿的液晶稳功率系统温漂抑制方法
US20130057429A1 (en) Positioning apparatus and signal processing method thereof
TW201336309A (zh) 符號率偵測裝置及符號率偵測方法
KR102286360B1 (ko) 에칭 처리 장치, 에칭 처리 방법 및 검출기
CN118031823A (zh) 晶圆的反射谱处理方法、装置、设备及存储介质
Dowsett et al. Quantification of secondary-ion-mass spectroscopy depth profiles using maximum entropy deconvolution with a sample independent response function
CN112729544A (zh) 一种扫描吸收光谱的激光波长锁定系统与方法
CN116519629A (zh) 一种高压环境下的tdlas测量信号谱线参数获取方法
CN114970639B (zh) 一种开放空间外界环境的气体闪烁噪声消除方法
CN111350494A (zh) 一种硬件聚焦阵列侧向自适应恒功率控制方法
CN115682964A (zh) 一种晶圆片厚度大量程、高精度快速解算方法
CN115389457A (zh) 基于变区域锁相放大的开放式激光气体检测方法及系统
CN109470647A (zh) 一种水蒸气太赫兹吸收谱的测量方法
TWI780618B (zh) 電漿處理裝置及電漿處理方法
CN113049120B (zh) 红外波段化学激光线宽的测量装置及其测量方法
KR101883298B1 (ko) 항공 수심 라이다의 수심 보정 방법
CN104776795B (zh) 基于可变衰减器的光学非线性一阶误差消除方法
JP2000292128A (ja) 半導体薄膜の膜厚測定装置および膜厚測定方法ならびに半導体ウェハおよびその製造方法

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