高等数字信号处理吴正国尹为民第5章短时富里叶变换课件教学.pptVIP

第五章 短时富里叶变换 第一节 时频分析的基本概念 第二节 短时富里叶变换 第三节 离散短时富里叶变换及其计 算 第四节 基于离散短时富里叶变换的信号重构 第五节 戈勃(Gabor)展开 第六节 循环平稳信号处理简介 第一节 时频分析的基本概念一、时频分析 ⒈连续时间函数的富氏变换 第二节 短时富里叶变換 第三节 离散短时富里叶变换及其计算 一、离散短时富里叶变换的定义 第四节 基于DSTFT的信号重构 一、滤波器组求和法 第六节 循环平稳信号处理简介 一、 循环平稳的基本概念 二、谱相关密度函数 1、谱相关密度函数 3、循环自相关函数 时变自相关函数： 循环自相关函数： 循环谱密度函数： 为矩形窗(时宽为T)下， 的富里叶变换，它仍 为t的函数。 对上述卷积应用卷积定理，其富里叶变换即 应为两频谱的乘积，考虑到不平稳性，代入 上述时变谱后，还应对其取时间平均，结果得： 为 的频谱在频域分别移动 后的两频谱间的互相关函数(谱相关密度函数) 谱相关系数 在频率 处的频谱自相关函数 谱相关系数 α为信号的循环频率时为1，其余小于1。 2、谱支撑域 令带限信号： 或 或 3、乘积的谱相关密度函数 例： X(t)的功率谱如图(b) 三、 循环统计量的估计 1、循环统计量的离散表示 二、离散短时富里叶变换的性质 性质1： 性质2：若x(n)的长度为N，w(n)的长度为M， 则的时域长度为N+M-1； 性质3： 的带宽总小于等于w(n)的带宽； 性质4：若x(n)、w(n)皆为因果序列，则在时 域上也是因果序列； 令m=l+n，l=k+sN，k=0,1,….,N-1，即将n时刻 窗内的信号分段，每段长度为N，则得： 使用N点FFT进行计算 三、 离散短时富里叶变换的计算 二、精确重构条件 时，仅有r=0项不为零 时，需对任何r不等于零的整数有： 频域等价条件： 第三节 戈勃展开 (Gabor’s Expansion ) 一、连续信号的戈勃展开 ⒈定义： 式中 称之戈勃基， 为戈勃展开系数，T和Ω分别 时间、频率取样周期。 TΩ= 2π为临界取样，TΩ 2π为过取样。 ⒉展开系数的计算 ①难点：戈勃基不一定是正交基。 ②完全重构条件 (Perfect Reconstruction Condition )：用于计算展开系数的辅助 函数γ( t ) 应满足如下约束条件 ①难点：戈勃基不一定是正交基。 ②完全重构条件 ： 用于计算展开系数的辅助 函数γ( t ) 应满足 如下约束条件 ③计算公式： 展开系数是STFT的二维取样值。 3.举例 ①母函数为高斯函数 特点：母函数在时、频域上皆集中分布，但其辅助函数在时域上却很分散。这样，戈勃展开系数的时间分辨率就很差。 ②母函数为矩形函数： 特点：母函数和辅助函数在时域上皆能集中分布，但在频域上却很分散。 二、离散信号的戈勃展开 1、离散信号戈勃展开的定义：长度为 有限 长度序列可以看作以 为周期的周期序列x( k )的主 值序列，x ( k )满足 则戈勃展开为 为给定窗序列 h ( k ) 的周期 延拓 戈勃展开式一般不唯一。并有双正交条件： 的长度等于 x ( k ) 的周期，无冗余。 但M、N多解， 2、临界取样： 3、举例 三、过取样条件下离散信号的戈勃展开 MATLAB实现：[tfr,dgr,gam]=tfrgabor (sig,N,q,h,trace) 信号长度为Nx， Gabor系数的时域长度为N， 频域长度为M，q=NM/Nx，h为戈勃基函数，长度 为N+1，默认为高斯窗，gam为对偶窗，tfr为戈勃 系数模平方，dgr为戈勃系数，trace非零，则显示 算法进程。 举例： 作业： 试证明(3.1.18)式。 试证明(3.2.3)式和(3.2.46)式。 试证明(3.3.21)式与(3.3.22)式一对变换对。 试证明(3.3.80)式和(3.3.81)式。 5. 令离散时间信号 ， w(n)为高斯窗函数。