高等数字信号处理课件作者吴正国尹为民第5章节短时富里叶变换.ppt

第五章 短时富里叶变换 第一节 时频分析的基本概念 第二节 短时富里叶变换 第三节 离散短时富里叶变换及其计 算 第四节 基于离散短时富里叶变换的信号重构 第五节 戈勃(Gabor)展开 第六节 循环平稳信号处理简介 第一节 时频分析的基本概念一、时频分析 ⒈连续时间函数的富氏变换 第二节 短时富里叶变換 第三节 离散短时富里叶变换及其计算 一、离散短时富里叶变换的定义 第四节 基于DSTFT的信号重构 一、滤波器组求和法 第六节 循环平稳信号处理简介 一、 循环平稳的基本概念 二、谱相关密度函数 1、谱相关密度函数 令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)为高斯窗函数。 试用MATLAB软件，取不 同时间长度的窗函数， 分别求x(n)的离散短时 富里叶变换，并进行信号 重构。试讨论窗函数 长度对时频分辨率、重构 精度的影响。 对题5的信号x(n)，信号长度N=1024，取高斯窗 函数，取不同的过取样率q=1、2、4，分别求 x(n)的戈勃展开系数。 1、一阶循环平稳信号 循环平稳信号：统计特性呈现周期性或多周期 (各周期不能通约)性平稳变化的非平稳随机信号。 时变均值的计算-------同步平均法 循环均值与循环频率 循环均值是将信号频谱左移循环频率后，再取时间平均。 只要信号频谱存在循环频率谱线，就不为零。 循环均值表征了对应的循环频率分量的幅值与相位。 可以检测 是否为零( )作为信号是否具有一阶循环平稳的判据。 存在多个不可通约的循环频率(多循环平稳) 时变均值： 循环均值： 循环均值可将频率不等于循环频率α的分量 (含零均值白噪声)置零 尚辅网