第三章_短时傅里叶变换与Gabor展开 文件.DOC

第三章 短时傅里叶变换与Gabor展开 问题： 怎样用FT对时间定位—平移加窗后再进行FT 连续STFT高度冗余，需进行采样以降低冗余度 选择高斯窗函数可实现最高时频分辨率—Gabor展开 怎样求窗函数的对偶函数？对偶函数时频定位性与采样率关系？ 拟正交Gabor展开法 周期Gabor展开：要求窗函数与数据长度一致 离散Gabor展开：降低了窗函数与数据长度的一致限制 3.1 短时傅里叶变换 一．短时傅里叶变换的引入 1.STFT信号处理步骤 ①将 ②对乘积结果作FT，获得t=0附近的频率特性 ③移动窗口，重复①②步骤，以获得其他位置的频率特性 2.定义： = 基本函数 3.物理意义： ①信号与 ② 4.特征 ①的时宽与带宽恒定，但两者分辨率需折中考虑 ②STFT谱与分析窗选择密切相关 ③二维时频面的标度由平移量和调制量决定 ④STFT的平方构成二维STFT谱 例3.1 ① 当时，上式，即时功率 STFT无新信息 ② 只要归一，则STFT谱的平均即为信号能量 ③平均即时频率 例3.2高斯函数 高斯分析窗 ①谱中心在（0,0） ② 比WVD分辨率要差 二．反STFT 1.表示方式 ①一维反变换： 令 ②从对偶函数h(t)进行反变换 若S(t)用另一窗函数h(t)处理后形成 = = 2.启示 ① ②STFT高度冗余，选其离散样本即可生成反变换 新问题 分析与综合函数不再相同 ③STFT可看作从时域到时频联合域的映射，但任一二维时频函数B(t,w)不一定存在反变换 三．离散信号的STFT 1.定义 令采样间隔 = 2.性质 ①离散STFT频域内具有周期性 ②任一二维离散函数不一定存在反变换 3.2 Gabor展开 一．Gabor展开 1.物理意义：存在占据有最小时频面积的基本信号（信息原子），任一信号可看成是由时频面内不同位置的基本信号叠加而成 2.基本函数 高斯函数（占据最小时频面积） ②任何函数均可作为基本函数（条件） 二．Gabor展开与STFT 问题：如何计算Gabor系数 1.1980 Bastias找到了高斯函数的对偶函数的解析解 2.系数求解—Gabor变换 3.两种Gabor展开形式—视具体应用选择 4.Gabor展开的缺陷 Gabor初始函数集无法形成正交基，即 不能直接求解相应对偶函数 三．对偶函数 利用Poisson求和公式 不能直接求出，很难得到解析解，时频分辨率也不好数值解法 3.3 周期离散Gabor展开 一．引入 1.Gabor展开：位移量，调频量离散了，信号分析与综合窗时间量还未离散化 离散化问题 时间离散 2.采样引起的问题：时域，频域均周期化 3.离散Gabor展开：（时频域进入步长） 采样 = 注：均为周期为L的序列（加“～”标志） 4.参数说明 过采样率 ，临界采样 ＞1，过采样 Gabor展开冗余 5.周期离散Gabor展开的DFT实现 同理： N—数字频率间隔数 6.周期离散Gabor展开的滤波器组解释 抽取（内插因子）小于频道数N，即＜N 二．离散周期Gabor展开对偶函数求解 1.对偶满足条件，将 要求上式双重级数求和项为函数 而其中 0，其它 利用离散Poisson求和公式，简化为 令 （p—频率调制量,q—时间平移量） 矩阵简化 离散Wexler—Raz等式 2.对偶的求解情况 ①临界采样：H为L×L方阵，对偶有唯一解，且与h(k)呈双正交关系 缺点：时频聚集性不好（Fig3–8, Fig3–9） ②过采样：,求解方程为欠定方程，对偶解不唯一。 3. 对时频分析影响 ①问题：Gabor基本函数不能构成正交基。 良好的时频聚集性 ②解决办法： a.选择好h(k) b.在采样情况下，选择最优 ③求解手段： a.Zak变换（局限：窗函数长度=数据长度） b．拟正交展开法（Qian） 3.4 拟正交Gab