第十章 信号的时频分析.pptVIP

本章结束 青岛大学机电工程学院 第十章 信号的时域分析 一、短时傅里叶变换（STFT）的基本原理 1、基本原理 第一节 短时傅里叶变换 目的：为了克服FT没有时间分辨率的缺陷。通常采用两种方法来适应分析非平稳信号的需要。方法一：在FT中引入时间相关性而又保持线性不变。其思想是引入一个“局部频率”参数（在某时间内局部）。这样，“局部”FT便是通过一个窗口来观察信号，在这个窗口内信号接近平稳，见下图。通过移动时窗进行分段取样，将整个信号化为若干个局部“平稳”的信号；对这“平稳”信号进行FT，可得到一组原信号的“局部”频谱。方法二：将FT中所用的正弦基函数改为在时间上更集中而在频率上较分散的基函数。 2、定义与性质 信号x(t’)的短时傅里叶变换或短时谱定义为 t x(t’) t’ 以t为中心的“分析窗” ▲ 短时傅里叶变换（或称短时谱）——信号x(t’)在时间t的短时FT就是x(t’)乘上一个以t为中心的“分析窗” （——表示复共轭），所作的FT。x(t’)叫做基信号。 分析窗的选择：乘 等价于取出信号在分析点t’=t附近的一个切片，所以STFT是x(t’)在t附近的“局部谱”。对x(t’)的STFT很大程度上取决于分析窗 。作为分析窗的时窗函数 应该是紧支集（即为有限长度）的或很快趋近于零的函数。满足此条件的“最优”窗函数是任一高斯函数，其他一些函数也可作为窗函数，如矩形窗、三角窗、Hanning窗和Hamming窗等。 为了对 的位置和大小进行定量描述，定义其中心th和宽度半径Δth为 ▲ STFT(t, f)描述了x(t’)在时窗[t’+th-Δth, t’+th+Δth]内的局部信息。 如果x(t’)和 的FT为X(f’)和 ，则可得到STFT的另一种表示式 ej2πt(f ’-f) 可以认为，STFT是“加窗谱” 的IFT，窗谱Г(f)是时间窗γ(t)的FT. STFT的时域表达 式(10-1) STFT的频域表达 式(10-4) 根据用模拟滤波方式求取频率特征的原理，式(10-4)中的STFT可以采用让信号x(t’)通过一频率响应为 的滤波器得出，参见下图。 (a) 用带通滤波器的实现方式。在某个给定的频率f处，x(t’)先通过一中心频率为f的“分析带通滤波器”，然后再将滤波器的输出频移至零频率。由于最后的频移， 作为t的函数，对任意固定的f都是一低通信号。 (b) 采用低通滤波器的实现方法。 它与“带通方式”是等价的。而低通滤波器的冲激响应等于时间反向的分析窗 (a) (b) *共轭 依据FT的时延和频移特性， 对于短时傅里叶变换受基信号变换的影响，根据式(10-1)不难证明： 及 表明STFT保持了信号x(t’)的频移特性 表明STFT在一调制(相位因子)范围内可保存时间移位 二、短时傅里叶变换的时间-频率分辨率 一维信号x(t) STFT 映射成在时间-频率平面(t，f)内的二维函数 根据不确定性原理或Heisenberg不等式，如果用Δt和Δf分别表示时间分辨率和频率分辨率，则它们的乘积满足不等式 面临的问题：高的时间分辨率要求分析窗 尽量窄，而高的频率分辨率要求滤波器的带宽尽可能窄，即所对应的分析窗则尽量宽。二者相矛盾。 表明时间分辨率和频率分辨率不可兼顾 1、理想的时间分辨率 分析窗为无穷窄，可选择δ(t)函数作为窗函数，则 2、理想的频率分辨率 用不变窗γ(t)≡1，则 三、时频信号展开与STFT综合 STFT退化为信号x(t)，它保留了信号的所有时间变化，失去了频率分辨率 此时，STFT变为傅里叶变换，没提供任何时间分辨率. 短时傅里叶分析的另一种描述方式是：将信号表示成基信号经时频移位之后的线性迭加，即 此式可认为是x(t’)以gt,f(t’)=g(t’-t)ej2πft’为“基信号”的时频展开式。Tx(t, f)为时频展开式的系数函数，它描述了信号x(t’)在时-频点(t, f)的强度。 式(10-10) 系数函数Tx(t, f)也可选作STFT，即 时频展开式(10-10)对任何有限能量的信号x(t’)都是存在的。将式(10-11)代入式(10-10)，可得到： 条件是：选择STFT分析窗 满足 。这个条件并不是一个很苛刻的条件。另外，选择分析窗的自由度表明，展开式的系数函数Tx(t, f)并不是唯一被定义的。 式(10－11) 式(10-11)叫做STFT分析，而式(10-12)叫做STFT综合。综合式可以看作是分析式的逆变换。