太赫兹波检测颗粒样品频谱散射线的去除算法制造技术

本发明专利技术涉及一种太赫兹波检测颗粒样品频谱散射线的去除算法，针对太赫兹波检测颗粒样品时由于散射效应的存在，导致频率谱线随着频率增大而缓慢变高，在后期应用中提取特征谱线中特征峰的位置和峰高时会引入较大误差，从而无法准确地描述被测样品在太赫兹波段内的光谱特性这一问题，从太赫兹系统采集得到的时域或者频域信号，建立对应的数学模型，推导频谱散射线与频率之间的对应关系，最终得到去除散射线之后的特征谱线。可以提高被测样品的特征峰和峰高的提取精度。该方法简单易行，适用范围广，准确度高。

【技术实现步骤摘要】

本专利技术涉及一种太赫兹波检测技术，特别涉及一种太赫兹波检测颗粒样品频谱散射线的去除算法。
技术介绍
太赫兹波一般是指频率在0.1到10THz(1THz＝1012Hz)范围内的电磁波，在电磁波谱中位于微波与红外波之间。由于在这一波段电磁波具有光子能量低、透射性高、携带信息量丰富等特性，太赫兹波在很多领域都显示了其巨大的应用价值，如在天文观测、物质识别、安检以及通信等。目前，主要的应用研究有太赫兹时域光谱技术、太赫兹成像技术、安全检查、太赫兹雷达、天文学、通信技术。目前，太赫兹系统可以分为采集得到的信息是时域信息和频域信息两类。采集得到的信息是时域信息时，将得到的太赫兹时域波形经过傅里叶变换后可以得到频域信息。根据两种系统得到的频域信息进行处理，获得材料在太赫兹波段的复介电常数，即色散及吸收等信息。但是，当太赫兹波照射到颗粒状样品上时会发生散射现象，导致频域谱线会随着频率的增大而缓慢上升，在后期应用中提取特征谱线中特征峰的位置和峰高时会引入较大误差，从而无法准确地描述被测样品在太赫兹波段内的光谱特性。但是目前，国内外使用太赫兹波对颗粒状样品进行物质成分测定分析时，只考虑了物质本身的吸收特性，而忽略了散射效应的影响，物质成分测定分析的精确度低。
技术实现思路
本专利技术是针对太赫兹波检测颗粒样品时由于散射效应的存在，导致频率谱线随着频率增大而缓慢变高，在后期应用中提取特征谱线中特征峰的位置和峰高时会引入较大误差，从而无法准确地描述被测样品在太赫兹波段内的光谱特性的问题，提出了一种太赫兹波检测颗粒样品频谱散射线的去除算法，从太赫兹系统采集得到的时域或者频域信号，建立对应的数学模型，推导频谱散射线与频率之间的对应关系，最终得到去除散射线之后的特征谱线。本专利技术的技术方案为：一种太赫兹波检测颗粒样品频谱散射线的去除算法，太赫兹系统测试信号为时域信号时，去除频谱散射线具体步骤包括为：1)检测系统在未放置样品的情况下，扫描得到太赫兹时域参考信号RTime(t)，对其进行傅里叶变换后，得到参考频域信息RTime(ω)，根据平均频率计算公式求得参考频谱信号平均频率μr_Time；2)在相同的检测系统中，放置被测样品后，扫描得到太赫兹时域测试信号SS_Time(t)，对其进行傅里叶变换后，得到测试频域信息SS_Time(ω)，根据平均频率计算公式求得测试频谱信号平均频率μS_Time；3)计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的特征谱线αTime(ω)；4)计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的频谱散射线αS_Time(ω)；5)用步骤3)所得实际测得的被测样品特征谱线减去频谱散射线αS_Time(ω)，得到最终的样品特征谱。所述太赫兹波检测颗粒样品频谱散射线的去除算法，太赫兹系统测试信号为频域信号时，去除频谱散射线具体步骤包括为：A：检测系统在未放置样品的情况下，直接获得参考频谱信号RFreq(ω)，根据平均频率计算公式，求得参考频谱信号平均频率μr_Freq；B：在相同的检测系统中，放置被测样品后，获得被测频谱信息SS_Freq(ω)；同样，根据平均频率计算公式，求得测试频谱信号平均频率μs_Freq；C：计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的特征谱线αFreq(ω)；D：计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的频谱散射线αS_Freq(ω)；E：用步骤C所得实际测得的被测样品特征谱线减去频谱散射线αS_Freq(ω)，得到最终的样品特征谱。所述太赫兹系统检测颗粒样品时的频谱散射线αS(ω)公式如下：αS(ω)=π(μr2-μS2)4dμr2μS2ω2,]]>其中μr为根据平均频率计算公式求得参考频谱信号平均频率，μS为根据平均频率计算公式求得测试频谱信号平均频率，d为被测样品的厚度。本专利技术的有益效果在于：本专利技术太赫兹波检测颗粒样品频谱散射线的去除算法，可以提高被测样品的特征峰和峰高的提取精度。该方法简单易行，适用范围广，准确度高。附图说明图1为本专利技术太赫兹时域光谱系统对样品进行太赫兹波检测时，未放置样品时，得到的参考时域图；图2为本专利技术太赫兹时域光谱系统对样品进行太赫兹波检测时，放置样品时，得到的参考时域图；图3为本专利技术太赫兹系统检测颗粒样品时的频谱散射线图；图4为本专利技术得到采集信号为时域信息的太赫兹系统进行太赫兹波检测时去除频谱散射线后的被测样品特征谱线图；图5为本专利技术使用红外光谱系统对样品进行太赫兹波检测时，未放置样品时，直接获得的参考频谱图；图6为本专利技术红外光谱系统对样品进行太赫兹波检测时，放置样品时，直接获得的参考频谱图；图7为本专利技术得到采集信号为频域信息的太赫兹系统检测颗粒样品时的频谱散射线图；图8为本专利技术使用红外光谱系统对颗粒样品进行太赫兹波检测时去除频谱散射线后的被测样品特征谱线图。具体实施方式实施例1:本专利技术所用的采集信号为时域信息的太赫兹系统，以太赫兹时域光谱系统(TDS)为例，其它采集信号为时域信息的太赫兹系统处理方法与之类似。所用的被测样品成分为C5爆炸物，厚度为1mm。TDS系统对样品进行太赫兹波检测时，未放置样品时，得到参考时域信息RTDS(t)，如图1。经过傅里叶变换后，可以得到参考频谱信号RTDS(ω)。根据下面平均频率计算公式，求得参考频谱信号平均频率μr_TDS：μ=∫0∞ωR(ω)dω∫0∞R(ω)dω.]]>在同一检测系统中，放入被测样品后，得到测试时域信息SS_TDS(t)，如图2。经过傅里叶变换后，可以得到测试频谱信号SS_TDS(ω)。根据相同的平均频率计算公式，求得测试频谱信号平均频率μS_TDS。采集信号为时域信息的太赫兹系统检测颗粒样品时的特征谱线为：αTDS(ω)=-1dln[SS_TDS(ω)RTDS(ω)]]]>其中，d＝1mm，为被测样品的厚度。根据太赫兹信号数学模型可以计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的频谱散射线，如图3：αS_TDS(ω)=π(μr_TDS2-μS_TDS2)4dμr_TDS2μr_TDS2ω2]]>最后，得到采集信号为时域信息的太赫兹系统对颗粒样品进行太赫兹波检测时去除频谱散射线后的被测样品特征谱线，如图4：α0(ω)＝αTDS(ω)-αS_TDS(ω)。实施例2:本专利技术所用的采集信号为频域信息的太赫兹系统，以傅里叶变换红外光谱系统(FTIR)为例，其它采集信号为频域信息的太赫兹系统处理方法与之类似。所用的被测样品成分为谷氨酸，厚度为1mm。使用FTIR对样品进行太赫兹波检测时，在不放置样品的情况下，可以直接获得参考频谱信号RFTIR(ω),如图5,根据平均频率计算公式，求得参考频谱信号平均频率μr_FTIR。在同一检测系统中，放入被测样品后，可以获得被测频谱信息SS_FTIR(ω),如图6。同样，根据平均频率计算公式，求得测试频谱信号平均频率μ本文档来自技高网...

【技术保护点】
一种太赫兹波检测颗粒样品频谱散射线的去除算法，其特征在于，太赫兹系统测试信号为时域信号时，去除频谱散射线具体步骤包括为：1)检测系统在未放置样品的情况下，扫描得到太赫兹时域参考信号RTime(t)，对其进行傅里叶变换后，得到参考频域信息RTime(ω)，根据平均频率计算公式求得参考频谱信号平均频率μr_Time；2)在相同的检测系统中，放置被测样品后，扫描得到太赫兹时域测试信号SS_Time(t)，对其进行傅里叶变换后，得到测试频域信息SS_Time(ω)，根据平均频率计算公式求得测试频谱信号平均频率μS_Time；3)计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的特征谱线αTime(ω)；4)计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的频谱散射线αS_Time(ω)；5)用步骤3)所得实际测得的被测样品特征谱线减去频谱散射线αS_Time(ω)，得到最终的样品特征谱。

【技术特征摘要】
1.一种太赫兹波检测颗粒样品频谱散射线的去除算法，其特征在于，太赫兹系统测试信号为时域信号时，去除频谱散射线具体步骤包括为：1)检测系统在未放置样品的情况下，扫描得到太赫兹时域参考信号RTime(t)，对其进行傅里叶变换后，得到参考频域信息RTime(ω)，根据平均频率计算公式求得参考频谱信号平均频率μr_Time；2)在相同的检测系统中，放置被测样品后，扫描得到太赫兹时域测试信号SS_Time(t)，对其进行傅里叶变换后，得到测试频域信息SS_Time(ω)，根据平均频率计算公式求得测试频谱信号平均频率μS_Time；3)计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的特征谱线αTime(ω)；4)计算得采集信号为时域信息的太赫兹系统检测颗粒样品时的频谱散射线αS_Time(ω)；5)用步骤3)所得实际测得的被测样品特征谱线减去频谱散射线αS_Time(ω)，得到最终的样品特征谱。2.根据权利要求1所述太赫兹波检测颗粒样品频谱散射线的去除算法，其特征在于，太赫兹系统测试信号为频域信号时，去除频谱散射线具体步骤包括为：A：...

【专利技术属性】
技术研发人员：王凯亮，彭滟，朱亦鸣，张秀平，徐博伟，张腾飞，王俊炜，庄松林，
申请(专利权)人：上海理工大学，
类型：发明
国别省市：上海;31