63连续Gabor变换_定义与窗函数.docVIP

6.3 连续Gabor变换：定义与窗函数 6.3.1 Gabor变换：定义与说明 一、(函数与广义函数 在式（6.2.1）中的函数，数学上称为第一个广义函数①，是狄喇克在量子物理的研究中提出的一个怪函数： 它的一种定义形式为 ．　　　　　　（6.3.1） 将式（6.3.1）代入式（5.2.1）： 　　　　　　　 （6.3.2） （对时间进行一定精度定位）， 再重写式（5.1.1a）于下 　　　　　　 （6.3.3） （对频率准确定位）， 对比分析以上两个式子：一个在一定精度下定位时间（即式（6.3.2）），其方法是在时间轴上用低通函数开出一个窗口——时间窗口，窗口持续时间用参量来控制，其定位功能则通过该窗口对信号在时间上进行扫描；另一个准确定位频率（即式（6.3.3））．对频率的准确定位需要宽时窗口：，而这又不能对时间精确定位，对时间的准确定位需要窄时窗口：． 我们可否折衷考虑取为一恰当的数值 放宽对时间准确定位要求，只要达到一定的精度即可，也就是说，将（6.3.5）中的固定——我们认为可以接受的某一确定值，在保证时间窗口函数满足低通性（low-pass property） （6.3.4） 的条件下，将分析式（6.3.2）改写成 　　　　　　　　　（6.3.5） （对时间具有一定定位能力） 对比式（6.3.3）与式（6.3.5），这两个变换式在形式又何其地相似！将两个式子融和在一起，就能达到同时对时间、对频率的一定定位能力呢？由此而来我们引入如下的Gabor变换． 二、连续Gabor变换 Gabor变换定义 对频率定位，必需保留频率变量，对时间定位必需保留扫描时间变量，因此，信号的一种时-频表示应当是将式（6.3.3）与（6.3.5）结合在一起： 　　　　　　　（6.3.6） 这就是连续Gabor变换定义，在许多地方称为窗口付里叶变换（windowed Fourier Transform），或短时傅里叶变换（Short Time Fourier Transform(STFT）——如果我们取窗口函数为矩形情形． 式（6.3.6）的物理意义是：Gabor变换在扫描时间附近局部地测量了频率为的正弦分量的幅度．也就是说，表示在窗口内观察到的那部分信号的频率结构情况．实际上就是用窗口函数先对信号在时间轴上进行扫描并定位，然后进行富利叶变换．对信号的开窗过程我们可用图6.3.1来说明．图6.3.1(a)表示了用小时宽窗口函数在信号时间和处开窗，得到该两时刻附近的信号和；图6.3.4(b)表示了用大时宽窗口函数在相同两时刻开窗情形． 显然，信号的Gabor变换精确地按窗口宽度分解了信号的频谱，提取出它的局部频谱信息，当在整个信号时间轴上扫描时，就给出了的完整频谱．因此，并没有损失信号在频域上的任何信息． 时-频窗口函数 式（6.3.6）中g通常称为时间窗口函数（time window function），简称时窗函数．除满足式（6.3.4）低通性外，为了同时具有一定的时间、频率定位与分辨率，时窗函数g应当具有好的时-频局域化特征，其时-频胞元面积越小越好，基于这一特性，将g称为时-频窗函数（time-frequency window function）更为贴切． (a)小时宽窗口t (b)大时宽窗口 图6.3.1　Gabor变换对信号的开窗 最初，Gabor选取g为高斯函数： （6.3.7） (　　　　　　　　　　( 　 　 　（6.3.8） 我们知道这是最佳的时-频窗口函数，它具有最小的胞元面积: ． 总结以上分析得出：时-频窗口函数g应当具有一定时-频局域化特征的低通函数．但并没有限制它的取值是实数还是复数，因此，更为一般化的连续Gabor变换定义成如下的形式 ．　　　　　（6.3.9a） GT的基函数——积分核、时-频分析胞元 将Gabor变换式（6.3.12a）写成内积形式则为 　　　　（6.3.9b） 其基函数，也即式（6.3.9a）中的积分核函数，当，时，有 ． （6.3.10） 改变与就可对整个时-频相平面中的任意信号进行扫描，得到信号的一种时-频表示，因而人们常常称为GT分析元（analysis cell），用它就可对时-频平面上分布的任何信号进行分析．在这里我们仿照的称谓——扫描时间，也叫为扫描频率（scanning frequency），如图6.3.2所示． 图6.3.2 GT时-频分析胞元： 一旦选定时-频窗函数g，其时-频胞元形状及其面积就被确定了，用它来分析信号，对时间、频率的定位准确性和分辨率也随之被确定了．当时-频指标（，）在内变化时，分析元的Heisenberg盒将覆盖整个