《Hilbert-Huang变换分析THz脉冲信号的时频特性》.pdf

第 26卷第 7期 强 激 光 与 粒 子 束 Vo1．26，NO．7 2014年 7月 HIGH POW ER LASER AND PARTICLE BEAM S Ju1．，2014 Hilbert—Huang变换分析 THz脉冲信号的时频特性 谢 静 ，谭佐军 ， 陈 阳 ， 石舒宁。， 陈 璐。 (1．华中农业大学 理学院 ，武汉 430070； 2．华 中农业大学 工学 院，武汉 430070) 摘 要： 提出一种基于 Hilbert—Huang变换 的THz时域光谱分析方法 ，将 THz时域脉冲信号分解成有 限数 目的单分量信号之和，利用 Hilbert变换求得瞬时频率来获得幅值 的时频分布——Hilbert—Huang变换谱 ， 实现 了通过水蒸气 的THz脉冲信号的时频分析 ，揭示了THz波与水蒸气相互作用 的频谱时域分布特性，并与 基于小波变换 的时频 图进行 了对 比分析 。结果表 明，该方法可以同时提高 THz脉冲时频分布 的时间分辨率和 频率分辨率 ，具有局部化分析和 自适应选择 的特点，还能直观地表现 出各频率成分之间的相对时间延迟 。 关键词 ： Hilbert—Huang变换 ； THz脉冲 ； 时频分析 ； Hilbert边 际谱 中图分类号 ： 0433．4 文献标志码 ： A doi：10．11884／HPLPB201426．073103 太赫兹 (THz)时域光谱信号口 是一种非平稳信号 ，在时间上的延迟 、信号强度的衰减、频谱的展宽等都反 映了被测物体的吸收、色散等特性 。传统的THz时域光谱分析技术是基于 Fourier变换_3]，将整个时域信号 变换到频域进行频谱分析，时频分离的做法掩盖了THz脉冲时间分辨的光谱特性 ，不能显示 出某个时间区域 的频率特性或某个频率特性发生的时间区域 。Gabor变换_6。在一定程度上解决了不能局部化分析的问题，可 以对 THz脉冲进行时间一频率联合分析，然而由于时频窗 口形状 固定 ，无法根据 THz脉冲频率的变化来改变 时频分辨率 ，且不能同时获得高的时间分辨率和频率分辨率 。小波变换l8具有可调的时频窗，可以非均匀地 划分 THz脉冲信号的时间轴和频率轴 ，在一定程序上克服了时间和频率分辨率不能兼顾的问题。但是小波变 换无法从 THz脉冲的局部特征出发选择小波基 ，且不同的小波基具有不同的性质 ，对 THz脉冲的分析能力也 不同，对同一 THz脉冲采用不同的小波基获得的结果没有可 比性 ，这造成了选择最优小波基的困难 。本文提 出一种基于 Hilbert—Huang变换 (HHT)的THz时域光谱分析方法。该方法将 THz时域脉冲信号分解成有 限数 目的单分量信号之和，再利用 Hilbert变换和求得 的瞬时频率来获得 HHT谱时频分布 ，可以同时提高 THz脉冲时频特性分析的时间分辨率和频率分辨率 。 1 Hilbert-Huang变换分析方法 HHT[1叩是一种针对非线性非平稳信 号的两步骤处理方法。第一步利用实验模态分解 (EMD)方法将 THz脉冲信号分解成有限数 目的单分量信号之和，即 ” -厂()一z(￡)+r() (1) i一 1 式中： ()为第 i个单分量信号，称为固有模态函数 (IMF)，即局部仅含有一个频率成分或者以某个频率为中 心的窄带信号； (￡)为残留分量 ，具有单调性或为常数。每个 IMF通过筛分来获取 ，从特征尺度出发 ，一步步 的把信号 中所包含的最精细的模态分离出来 ，即 z (f)=h 一h姒一1)一 m* (2) 式中：m 为第k次筛分的上和下包络线的均值函数 ；h 为第k次筛分的结果 ；h。为经过 i一1次分解后的残 余信号，