CN105300920B - 一种基于太赫兹反射光谱提取固体薄片复折射率的方法 - Google Patents
一种基于太赫兹反射光谱提取固体薄片复折射率的方法 Download PDFInfo
- Publication number
- CN105300920B CN105300920B CN201510367771.9A CN201510367771A CN105300920B CN 105300920 B CN105300920 B CN 105300920B CN 201510367771 A CN201510367771 A CN 201510367771A CN 105300920 B CN105300920 B CN 105300920B
- Authority
- CN
- China
- Prior art keywords
- sample
- signal
- terahertz
- domain
- time
- 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
本申请公开了一种基于太赫兹反射光谱提取固体薄片复折射率的方法，制备同一种材料的薄厚不同的两个固体薄片样品，其中薄样品为被测样品，厚样品作为时域信号基线校正时的辅助样品，测量二者的太赫兹时域反射信号；用被测样品和辅助样品的太赫兹时域反射信号之差作为基线，将被测样品的第一、二个峰分别作为基线校正后的参考信号和样品信号；对基线校正后的样品信号和参考信号分别进行复合数字滤波，以去除系统回波和随机噪声；根据滤波后的样品信号和参考信号，用自参考方法提取被测样品的复折射率。
Description
技术领域
本申请涉及光子学技术与电子学技术领域，尤其涉及一种基于太赫兹反射光谱提取固体薄片复折射率的方法。
背景技术
太赫兹波通常指的是频率在0.1THz～10THz(1THz＝1012Hz)，波长在3mm～30μm范围内的电磁辐射。太赫兹波在电磁波谱上位于毫米波和红外线之间，属于远红外波段，太赫兹频段在电磁波谱中的位置如图1所示，是光子学技术与电子学技术、宏观与微观的过渡区域。太赫兹辐射具有非常重要的科学价值，目前已成为世界性的研究热点，是2l世纪科学研究最前沿的领域之一。太赫兹频段所具有的独特性质决定了它具有广阔的应用前景。
利用一些光谱技术可以对药物进行检测，目前已采用拉曼散射光谱、红外光谱对药物进行了一些研究。相比于这些技术，太赫兹时域光谱技术具有独特的优势，例如，它对化合物晶型有很高的灵敏度，能够反映整个分子的低频振动和分子间的弱相互作用，能量低不会造成药物发生化学变化等。因此，利用太赫兹光谱在化学药品的生产、储存和流通环节进行质量的检查等，成为太赫兹辐射最具前景的应用之一。
化学药品通常以固体薄片、颗粒和液态形式存在。通过对化学药品薄片的太赫兹时域信号进行快速傅立叶变换，可以同时获得它的太赫兹信号的振幅和相位信息，从而提取到物质在太赫兹频段的复折射率，即获得药品的“指纹谱”，从而达到分析和鉴别药品的目的。
现有技术中，通常采用全反射镜参考方法提取固体薄片的复折射率。图2(a)所示为太赫兹脉冲垂直入射到固体薄片样品上时光线的传播情况，图中表示被测固体薄片放置在干燥空气中。若忽略空气对太赫兹辐射的吸收，则空气的折射率为实数，且近似为1；设薄片的复折射率为其中ns为样品的折射率，ks为消光系数。
现有的全反射镜参考方法的处理流程如图3所示，包括如下步骤：
步骤301：基于反射式太赫兹时域光谱系统，测量全反射镜的时域信号并将其作为参考信号对其做傅里叶变换得到参考信号的频谱
步骤302：将全反射镜取下换成样品，并使得二者的前表面位置完全重合，测量样品的时域信号同样再对其做傅里叶变换，得到样品的频谱
步骤303：根据公式(6)计算得到样品的折射率ns和消光系数ks，从而获得复折射率
假设全反射镜无损耗反射，则全反射镜的反射信号等于入射信号以图2(a)中作为样品信号，将作为参考信号。则参考信号和样品信号的频谱可以分别表示为
其中，可由菲涅尔(Fresnel)公式得到，通常采用P分量探测，则
用样品信号频谱除以参考信号频谱可得到传递函数的表达式
联立公式(3)和(4)，得到样品复折射率和传递函数的关系
是复数，令则由公式(5)可以推导出样品的折射率和消光系数的表达式
吸收系数αs＝2ωks/c，其中c是真空中的光速。传递函数的振幅ρ(ω)和相位N(ω)都是可测量的量，因此，由公式(6)便可以获得折射率ns和消光系数ks，从而获得复折射率
这种全反射镜参考方法存在如下问题：首先，传递函数的相位φ(ω)不可避免地会存在一个机械相移误差。如图2(b)所示，样品表面和全反射镜的表面很难实现完全复位，二者的位置会存在△L的光程差。因此，参考信号与样品信号会存在一个固有的相移使得提取的固体薄片的复折射率不够精确；其次，该方法要求测量系统的机构精密度足够高，因此，增加了系统硬件制作的难度；最后，由于需要分别测量参考与样品的时域信号，因此，测量效率低。
发明内容
本申请提供了一种基于太赫兹反射光谱提取固体薄片复折射率的方法，可以提高测量得到的复折射率的精确度，并提高测量效率。
本申请实施例提供的一种基于太赫兹反射光谱提取固体薄片复折射率的方法，包括：
A、制备同一种材料的薄厚不同的两个固体薄片样品，其中薄样品为被测样品，厚样品作为时域信号基线校正时的辅助样品，测量二者的太赫兹时域反射信号；
B、用被测样品和辅助样品的太赫兹时域反射信号之差作为基线，将被测样品的第一、二个峰分别作为基线校正后的参考信号和样品信号；
C、对基线校正后的样品信号和参考信号分别进行复合数字滤波，以去除系统回波和随机噪声；
D、根据滤波后的样品信号和参考信号，用自参考方法提取被测样品的复折射率，其中包括：采用遗传算法对公式
进行优化求解，以传递函数为目标函数，进行N次种群迭代，最终得到优化后的复折射率的精确解；其中，为传递函数的相位；为传递函数的振幅；菲涅尔系数ras、tas、rsa、tsa可由菲涅耳公式得到：
表示复折射率，c是真空中的光速，d表示被测样品厚度，下标as表示由空气到样品的界面，下标sa表示由样品到空气的界面。
较佳地，步骤B包括：
B1、将被测样品和辅助样品的太赫兹时域反射信号的第一个峰在时间轴上对齐；
B2、用被测样品的太赫兹时域反射信号减去辅助样品的太赫兹时域反射信号，得到样品信号；
B3、将被测样品的太赫兹时域反射信号中的第二个峰置为前后数据点的平均值，将辅助样品的太赫兹的时域反射信号中的第一个峰置为前后数据点的平均值，再用被测样品的太赫兹时域反射信号减去辅助样品的太赫兹的时域反射信号，得到参考信号。
较佳地，所述复合数字滤波包括解卷积和小波变换。
较佳地，所述解卷积包括：
通过实验得到第一个回波的振幅A1和第一个回波与主峰的时间延迟(t1-t0)；
依据公式滤波，获得主峰的频谱；其中，A(ω)和A0(ω)分别为原始时域信号和主峰的频谱。
较佳地，所述小波变换包括：将基线校正后的样品信号和参考信号多尺度分解，使噪声按频率区域划分；采用Matlab中的wdencmp函数，对高频系数进行阈值量化去噪处理，最后根据小波分解的各层系数和经过量化处理后的高频系数，进行一维信号的重构，得到去除噪声后的样品信号和参考信号。
较佳地，所述N大于或等于100。
较佳地，由公式计算得到的n0(ω)、k0(ω)作为遗传算法的初代种群。
较佳地，所述被测样品和辅助样品的厚度差的绝对值等于或大于0.8mm。
从以上技术方案可以看出，与目前的全反射镜参考方法相比，本发明的方法具有以下优点：
(1)仅需一次性制备同种材料的厚样品，以后再测相同材料时可重复使用作为基线校正的辅助样品，从而提高了测量效率；
(2)采用数字滤波技术，可以通过编程灵活地调节参数，不受硬件的限制，并且可以采用最先进的数字滤波算法；
(3)无需用全反射镜作为参考，避免了样品表面与参考物不共面所带来的机械相移误差，提高了复折射率的测量精度；
(4)将被测固体薄片的前表面反射的光作为参考信号，将后表面反射的光作为样品信号，光经过样品内部后将携带更多的样品信息，便于提取样品的特征吸收峰。
附图说明
图1为太赫兹频段在电磁波谱中的位置示意图；
图2(a)为太赫兹脉冲垂直入射到固体薄片上的光线传播示意图，图2(b)为相应的反射相移示意图；
图3为现有的全反射镜参考方法的处理流程示意图；
图4为本申请提供的基于太赫兹反射光谱提取固体薄片复折射率的方法的流程示意图；
图5所示为全反射镜与厚度为2.588mm、直径为13mm的葡萄糖薄片的太赫兹时域信号；
图6所示为采用现有的全反射镜参考方法提取的葡萄糖薄片的折射率(a)和消光系数(b)；
图7所示为厚度分别为2.588mm和3.750mm的葡萄糖薄片的太赫兹反射时域信号；
图8所示为厚度为2.588mm的葡萄糖薄片去除基线前后的时域信号：(a)样品信号；(b)参考信号；
图9所示为厚度为2.588mm的葡萄糖薄片第一个峰的频谱复合滤波前后的比较：(a)测量的原始频谱；(b)小波去噪后的频谱；(c)解卷积后的频谱；(d)复合滤波后的频谱；
图10所示为厚度为2.588mm的葡萄糖薄片的折射率和吸收系数随频率变化曲线。
具体实施方式
为使本申请技术方案的技术原理、特点以及技术效果更加清楚，以下结合具体实施例对本申请技术方案进行详细阐述。
图4为本申请提供的基于太赫兹反射光谱提取固体薄片复折射率的方法的流程示意图，该流程包括：
步骤401：制备同一种材料的薄厚不同的两个固体薄片样品，测量二者的太赫兹时域反射信号，其中薄样品为被测样品，厚样品作为时域信号基线校正时的辅助样品，薄厚样品的厚度差的绝对值应等于或大于0.8mm；
步骤402：用被测样品和辅助样品的太赫兹时域反射信号之差作为基线，将被测样品的第一、二个峰分别作为基线校正后的参考信号和样品信号；
步骤403：对基线校正后的样品信号和参考信号分别进行复合数字滤波，以去除系统回波和随机噪声；
步骤404：根据滤波后的样品信号和参考信号，用自参考方法提取被测样品的复折射率，其中包括：采用遗传算法进行优化求解，以传递函数为目标函数，进行N次种群迭代，最终得到优化后的复折射率的精确解。其中，N大于等于100。
其中，步骤402的作用是消除基线噪声。具体可以包括如下子步骤：
子步骤402-1：将被测样品和辅助样品的太赫兹时域反射信号的第一个峰在时间轴上对齐；
子步骤402-2：用被测样品的太赫兹时域反射信号减去辅助样品的太赫兹时域反射信号，得到样品信号；
子步骤402-3：将被测样品的太赫兹时域反射信号中的第二个峰置为前后数据点的平均值，将辅助样品的太赫兹的时域反射信号中的第一个峰置为前后数据点的平均值，再用被测样品的太赫兹时域反射信号减去辅助样品的太赫兹的时域反射信号，得到参考信号。
步骤403中的复合数字滤波包括解卷积和小波变换。
太赫兹波在砷化镓(GaAs)基底和碲化锌(ZnTe)探测晶体中的多次反射会产生系统回波，给太赫兹光谱带来系统误差。本申请方案中利用解卷积算法去除系统回波。
假定考虑到测量是在干燥空气条件下进行，认为THz脉冲在真空中传播的过程中无耗散，因此系统回波脉冲等于主脉冲的相对强度乘以时间延迟：
其中，A(t)和A0(t)分别为原始时域信号和主脉冲；Ai和ti分别表示第i个回波相对于主脉冲的幅值和延迟时间，当i＝0表示主脉冲；为卷积运算符。理论上说，系统回波的数目是无限的，其振幅依次衰减，计算时考虑的回波数目越多计算结果越精确，但是计算过程越复杂。本申请实施例中，只考虑较强的第一个系统回波，即取i＝1。将时域信号表达式(8)进行傅里叶变换，并取i＝1，可以得到
其中，A(ω)和A0(ω)分别为原始时域信号和主峰的频域表达式；(t1-t0)为系统的第一个回波与主峰的时间延迟。由公式(9)可以看出，A(ω)相当于在A0(ω)基础上加正弦振荡。因此，只要通过实验得到第一个回波的振幅A1和其与主峰的时间延迟(t1-t0)，就可以去除系统回波噪声。
本申请方案利用小波变换去除随机噪声，小波变换是把函数在特定的空间内按照小波基展开和逼近，可以同时进行时域和频域的分析，并具有很强的表征信号局部特征的能力。因为小波变换去噪对于具有尖峰和突变特性的太赫兹时域信号保真度高。根据小波变换理论将噪声信号进行多尺度分解，使噪声按频率区域划分，再分别用滤波器对信号进行滤波，用得到的近似系数和精确系数进行小波重构，最后得到去除随机噪声的真实时域信号。对于信号f(t)，其小波变换的形式为
其中，ψ(a,b)(x)为母小波ψ(x)经过伸缩和平移得到的小波函数，a表示伸缩因子，b表示平移因子。小波变换ωf(a,b)体现的是原信号在x＝b处、在x-b≤a范围内对信号进行分析，具有时频局域化和多分辨的特性。
一个含噪声的一维时域信号可以认为是真实信号叠加了噪声信号，以如下方式表示：
S(t)＝f(t)+N(t),t＝0,1,...,n-1, (11)
其中，S(t)为测量信号，f(t)为真实信号，N(t)为噪声信号。在实际应用中，噪声多表现为高频信号。因此，在用小波分析时，将一维时域信号S(t)(基线校正后的时域样品信号和时域参考信号)进行分解后，采用Matlab中的wdencmp函数，对高频系数进行阈值量化去噪处理，最后根据小波分解的各层系数和经过量化处理后的高频系数，进行一维信号的重构，得到去除噪声后的时域样品信号和时域参考信号。
以下给出步骤404所述的一种提取固体薄片太赫兹频段复折射率的自参考方法的具体实现过程，该方法以图2中的作为时域参考信号，而以作为时域样品信号，将二者分别进行傅里叶变换，得到各自的频谱复折射率的求解公式推导过程如下
其中，P分量菲涅尔系数ras、tas、rsa、tsa可由菲涅耳公式得到
下标as表示由空气到样品的界面，下标sa表示由样品到空气的界面。
传播因子ps[ω,Ls(ω)]的表达式为，
d表示被测样品厚度。
传递函数为
将代入上式，并化简得到
若令则有
由上述等式可推导出
abs()表示对括号内求绝对值，ln[]表示对括号内求自然对数。
对于弱吸收介质，即ns>>ks，可以忽略样品复折射率的虚部ks，则 从而得到样品复折射率的近似解
利用公式(20)，采用遗传算法，通过MATLAB编程计算，便可以提取出样品的复折射率。
较佳地，将由公式(21)计算得到的n0(ω)、k0(ω)作为遗传算法的初代种群，这样便于提高计算效率，但是，最终的优化结果与赋给的初值无关。
该自参考方法不仅避免了反射太赫兹时域光谱系统的固有相移问题，而且由于光程增加而提取了更强的物质属性信息。另外，作为基线校正的厚样品一次性制备完成后，或者当测量基线完好，并不需要校正时，仅需测量一次样品的时域信号，便可以获得复折射率的值。因此，实验效率提高了。这对于太赫兹光谱技术的实际应用具有重要意义。本发明的方法不局限于固体薄片，也可以推广到液态样品复折射率的提取。若将颗粒样品压制成薄片，同样可以应用本发明方法提取颗粒样品的复折射率。
以下以厚度为2.588mm、直径为13mm的无水D-葡萄糖粉末的固体薄片的实验为例，说明本申请方法与传统方法的区别。
I.现有的全反射镜参考方法
现有的全反射镜参考方法要求参考全反射镜表面与样品前表面的位置完全重合。图5为全反射镜与厚度为2.588mm、直径为13mm的葡萄糖薄片的太赫兹时域信号，小图为第一个峰放大后的曲线。如图5中的小图可见，从波峰位置看到，金镜和各样品的表面并没有完全重合。利用公式(1)-(6)提取的葡萄糖薄片的折射率和消光系数，如图6所示。
II.本发明方法
1)光谱相减法基线校正
在制备了被测样品(厚度为2.588mm、直径为13mm的葡萄糖粉末的固体薄片)后，还制备了厚度为3.750mm、直径为13mm的同一种葡萄糖粉末的固体薄片，用于时域信号基线校正时的辅助样品。
图7展示的分别是纯葡萄糖不同厚度薄片的两条时域信号，粗线、细线分别是厚度为3.750mm和2.588mm的薄片的时域信号，图中的虚线为零线。可以明显看到，随着测量时间的增加，谱线向下漂移，表现出基线存在误差。本发明对原始数据进行如下基线校正：
(1)将厚度不同的纯葡萄糖薄片的时域信号的第一个峰位通过调节时间手动重合；然后，根据数据点统一将其时间截到60ps；最后，用薄样品的时域信号减去厚样品的时域信号，结果如图8(a)中细线所示，粗线为原始时域信号；
(2)将薄样品时域信号中的第二个峰置为前后数据点的平均值，将厚样品时域信号中的第一个峰置为前后数据点的平均值；然后，用薄样品的时域信号减去厚样品的时域信号，即得到参考信号。如图8(b)中的细线所示，粗线为测量的原始时域信号。
综上所述，采取不同厚度的样品时域谱来消除基线噪声方法的具体使用步骤可以归纳为：首先，将不同厚度样品的两条时域谱的第一个峰手动对齐；然后，用薄样品的时域信号减去厚样品的时域信号，得到样品信号；最后，将薄样品时域信号中的第二个峰置为前后数据点的平均值，将厚样品时域信号中的第一个峰置为前后数据点的平均值，用薄样品的时域信号减去厚样品的时域信号，得到参考信号。
由图8可见，虽然经过基线校正后，系统回波信号幅度有所减弱，但是依然存在。因此，还要考虑消除系统回波。
2)复合滤波法去噪
图9是厚度为2.588mm厚的葡萄糖薄片的复合滤波前后的比较。按照上述解卷积和小波变换的原理进行滤波去噪，得到如图9所示的第一个峰的傅里叶变换结果。
3)自参考法提取复折射率
将参考信号和样品信号进行傅里叶变换获得各自的频谱图，两式相除获得传递函数，按照自参考的原理计算获得厚度为2.588mm厚的葡萄糖薄片的折射率和吸收系数，结果如图10所示。可以看出，葡萄糖的折射率在1.3THz后呈现反常色散，且吸收系数随着频率的增加而增大，并存在特征吸收峰。
III.本发明方法与传统全反射镜参考方法结果比较
图6和图10分别为用传统的全反射镜参考方法和本发明方法得到的葡萄糖薄片的折射率和消光系数。由图6(a)可见，折射率变化比较平缓，图6(b)中的消光系数也没有表现出明显的特征吸收，二者没有满足的介电模型。然而，从图10(a)中可见，折射率有多处突变，且突变位置与图10(b)中的消光系数的峰位相关，峰位在图中标出，二者满足洛伦兹介电模型。这些吸收峰是无水D-葡萄糖的特征吸收峰，从特征吸收峰我们可以鉴别是否为无水D-葡萄糖这种物质。
以上所述仅为本申请的较佳实施例而已，并不用以限制本申请的保护范围，凡在本申请技术方案的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请保护的范围之内。
Claims (7)
1.一种基于太赫兹反射光谱提取固体薄片复折射率的方法，其特征在于，包括：
A、制备同一种材料的薄厚不同的两个固体薄片样品，其中薄样品为被测样品，厚样品作为辅助样品，测量二者的太赫兹时域反射信号；
B1、将被测样品和辅助样品的太赫兹时域反射信号的第一个峰在时间轴上对齐；
B2、用被测样品的太赫兹时域反射信号减去辅助样品的太赫兹时域反射信号，得到样品信号；
B3、将被测样品的太赫兹时域反射信号中的第二个峰置为前后数据点的平均值，将辅助样品的太赫兹的时域反射信号中的第一个峰置为前后数据点的平均值，再用被测样品的太赫兹时域反射信号减去辅助样品的太赫兹的时域反射信号，得到参考信号；
C、对所述样品信号和参考信号分别进行复合数字滤波，以去除系统回波和随机噪声；
D、根据滤波后的样品信号和参考信号，用自参考方法提取被测样品的复折射率，其中包括：采用遗传算法对公式
进行优化求解，以传递函数为目标函数，进行N次种群迭代，最终得到优化后的复折射率的精确解；其中，为传递函数的相位；为传递函数的振幅；菲涅尔系数ras、tas、rsa、tsa可由菲涅耳公式得到：
表示复折射率，c是真空中的光速，d表示被测样品厚度，下标as表示由空气到样品的界面，下标sa表示由样品到空气的界面，j表示虚数单位，ω表示频率，ns(ω)表示样品的折射率，ks(ω)表示样品的消光系数。
2.根据权利要求1所述的方法，其特征在于，所述复合数字滤波包括解卷积和小波变换。
3.根据权利要求2所述的方法，其特征在于，所述解卷积包括：
通过实验得到第一个回波的振幅A1和第一个回波与主峰的时间延迟(t1-t0)；
依据公式滤波，获得主峰的频谱；其中，A(ω)和A0(ω)分别为原始时域信号和主峰的频谱。
4.根据权利要求2所述的方法，其特征在于，所述小波变换包括：将基线校正后的样品信号和参考信号多尺度分解，使噪声按频率区域划分；采用Matlab中的wdencmp函数，对高频系数进行阈值量化去噪处理，最后根据小波分解的各层系数和经过量化处理后的高频系数，进行一维信号的重构，得到去除噪声后的样品信号和参考信号。
5.根据权利要求1所述的方法，其特征在于，所述N大于或等于100。
6.根据权利要求1所述的方法，其特征在于，由公式计算得到的n0(ω)、k0(ω)作为遗传算法的初代种群。
7.根据权利要求1所述的方法，其特征在于，所述被测样品和辅助样品的厚度差的绝对值等于或大于0.8mm。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510367771.9A CN105300920B (zh) | 2015-06-29 | 2015-06-29 | 一种基于太赫兹反射光谱提取固体薄片复折射率的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510367771.9A CN105300920B (zh) | 2015-06-29 | 2015-06-29 | 一种基于太赫兹反射光谱提取固体薄片复折射率的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105300920A CN105300920A (zh) | 2016-02-03 |
| CN105300920B true CN105300920B (zh) | 2018-05-08 |
Family
ID=55198430
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510367771.9A Expired - Fee Related CN105300920B (zh) | 2015-06-29 | 2015-06-29 | 一种基于太赫兹反射光谱提取固体薄片复折射率的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105300920B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105823756B (zh) * | 2016-03-24 | 2018-11-13 | 西安电子科技大学 | 一种金属太赫兹至远红外复折射率的联合反演方法 |
| CN106644073A (zh) * | 2016-11-07 | 2017-05-10 | 北京师范大学 | 一种消除太赫兹光谱的水汽噪声的方法 |
| CN106841095A (zh) * | 2017-01-04 | 2017-06-13 | 北京环境特性研究所 | 一种使用太赫兹脉冲测量材料参数和材料厚度的方法 |
| CN108489917B (zh) * | 2018-02-28 | 2020-09-01 | 浙江大学 | 一种去除太赫兹吸收光谱中低频误差和高频噪声的方法 |
| CN109932338B (zh) * | 2019-04-03 | 2021-06-04 | 北京环境特性研究所 | 基于太赫兹频段测量样品复折射率的方法和装置 |
| CN111998783B (zh) * | 2020-07-13 | 2021-12-31 | 北京工业大学 | 基于进化优化算法的反射式太赫兹时域光谱厚度测量方法 |
| CN112162295B (zh) * | 2020-09-23 | 2022-07-15 | 青岛青源峰达太赫兹科技有限公司 | 一种基于时频分析的太赫兹厚度检测优化方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101419157A (zh) * | 2008-09-05 | 2009-04-29 | 中国计量学院 | 一种利用太赫兹时域谱精确测定食用油光学参数的方法 |
| CN103196868A (zh) * | 2013-04-10 | 2013-07-10 | 上海理工大学 | 光刻胶折射率的测定方法 |
| CN103645154A (zh) * | 2013-11-26 | 2014-03-19 | 东莞理工学院 | 一种利用太赫兹光谱信号提取材料光学常数的方法 |
| CN104076009A (zh) * | 2014-07-02 | 2014-10-01 | 中国人民解放军电子工程学院 | 一种生物颗粒远红外波段复折射率的快速测量方法 |
| CN104237143A (zh) * | 2013-06-08 | 2014-12-24 | 青岛科技大学 | 一种基于太赫兹光谱的固体农药鉴别方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6683681B2 (en) * | 2002-04-10 | 2004-01-27 | Baker Hughes Incorporated | Method and apparatus for a downhole refractometer and attenuated reflectance spectrometer |
| JP2013174548A (ja) * | 2012-02-27 | 2013-09-05 | Otsuka Denshi Co Ltd | 測定装置および測定方法 |
-
2015
- 2015-06-29 CN CN201510367771.9A patent/CN105300920B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101419157A (zh) * | 2008-09-05 | 2009-04-29 | 中国计量学院 | 一种利用太赫兹时域谱精确测定食用油光学参数的方法 |
| CN103196868A (zh) * | 2013-04-10 | 2013-07-10 | 上海理工大学 | 光刻胶折射率的测定方法 |
| CN104237143A (zh) * | 2013-06-08 | 2014-12-24 | 青岛科技大学 | 一种基于太赫兹光谱的固体农药鉴别方法 |
| CN103645154A (zh) * | 2013-11-26 | 2014-03-19 | 东莞理工学院 | 一种利用太赫兹光谱信号提取材料光学常数的方法 |
| CN104076009A (zh) * | 2014-07-02 | 2014-10-01 | 中国人民解放军电子工程学院 | 一种生物颗粒远红外波段复折射率的快速测量方法 |
Non-Patent Citations (2)
| Title |
|---|
| A reliable method for extraction of material parameters in terahertz time-domain spectroscopy;Lionel Duvillaret et al.;《IEEE JOURNAL OF SELECTED TOPICS IN QUANTUM ELECTRONICS》;19960930;第2卷(第3期);第739-746页 * |
| 太赫兹时域光谱技术;孟田华 等;《山西大同大学学报(自然科学版)》;20071231;第23卷(第3期);第14-17页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105300920A (zh) | 2016-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105300920B (zh) | 一种基于太赫兹反射光谱提取固体薄片复折射率的方法 | |
| US5939721A (en) | Systems and methods for processing and analyzing terahertz waveforms | |
| Mittleman et al. | Recent advances in terahertz imaging | |
| CN109100325B (zh) | 一种基于光谱吸收率二次谐波特征提取的气体浓度测量方法 | |
| JP6022456B2 (ja) | 時間領域信号の減衰定数を決定する方法とシステム | |
| Naftaly et al. | A method for removing etalon oscillations from THz time-domain spectra | |
| Bernier et al. | Precise determination of the refractive index of samples showing low transmission bands by THz time-domain spectroscopy | |
| CN112067574A (zh) | 不等光程法太赫兹光谱测量方法及其系统 | |
| CN103645154A (zh) | 一种利用太赫兹光谱信号提取材料光学常数的方法 | |
| Zhai et al. | Terahertz time-of-flight tomography beyond the axial resolution limit: autoregressive spectral estimation based on the modified covariance method | |
| CN107831157A (zh) | 基于除谱的拉曼光谱荧光背景扣除方法 | |
| Fujiyama et al. | Vibrational intensities XXV: some systematic errors in infrared absorption spectrophotometry of liquid samples | |
| Hirsch et al. | Techniques for cancellation of interfering multiple reflections in terahertz time-domain measurements | |
| Arbab et al. | Application of wavelet transforms in terahertz spectroscopy of rough surface targets | |
| Vermare et al. | Density fluctuation measurements using X-mode fast sweep reflectometry on Tore Supra | |
| Yang et al. | A background reduction method based on empirical mode decomposition for tunable diode laser absorption spectroscopy system | |
| CN106908422B (zh) | 一种非共线光参量放大荧光光谱仪的数据采集方法 | |
| Jördens et al. | Fibre-coupled terahertz transceiver head | |
| Angeluts et al. | Reconstruction of the spectrum of the relative permittivity of the plane-parallel plate from the angular dependences of its transmission coefficients | |
| Gu et al. | Extraction of optical constants in the terahertz band using material dispersion models | |
| Zhang et al. | The elimination of the FP effect in terahertz signals from a parameter-optimized decomposition | |
| CN114152589A (zh) | 一种基于双光频梳光源的太赫兹高光谱成像系统 | |
| Guan et al. | Terahertz time-domain spectral signal noise reduction analysis based on DG-SSA-VMD algorithm | |
| Wang et al. | A novel method to process the THz time-domain spectrum and software to make data visualization | |
| Zhang et al. | Quasi-Optical Measurement and Complex Refractive Index Extraction of Flat Plate Materials Using Single Time-Domain Transmission Model in Y-Band |
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 |
Granted publication date: 20180508 Termination date: 20210629 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |