短时傅里叶变换stft.pptVIP

短时傅里叶变换 ? FT 在信号处理中的局限性 : ? 用傅立叶变换提取信号的频谱需要利用信号的全部 时域信息。 ? 傅立叶变换没有反映出随着时间的变化信号频率成 分的变化情况。 1. 短时傅里叶变换简介 ? 提出与基本思想 ? 鉴于傅里叶变换的缺陷提出了窗函数的概念，提出一个灵活 可变的时间 - 频率窗，使得在这个窗内能够体现频率的信息， 这种信号分析方法称为时间 - 频率分析。而窗固定的时间 - 频 率分析方法即为短时傅里叶变换。 ? 短时傅里叶变换（ STFT ， short-time Fourier transform ）。 其主要思想是将信号加窗，将加窗后的信号再进行傅里叶变 换，加窗后使得变换为时间 t 附近的很小时间上的局部谱， 窗函数可以根据 t 的位置变化在整个时间轴上平移，利用窗 函数可以得到任意位置附近的时间段频谱实现时间局域化。 ? STFT 定义 : 1946 年， Gabor 就提出了 STFT ，给定一信号，其 STFT 定义为 : ? 短时谱的特点： 1) 时变性：既是角频率 ω 的函数又是时间 t 的函数。 2) 周期性：是关于 ω 的周期函数，周期为 2 π 。 ? ? ?? ? ? ? ? ?? ?? ? ? ? ? ? ? j j x e t m x d e t m x t S ) ( ), ( ) ( ) ( ) , ( 窗函数 （ 1.1 ） ? 公式涵义： ? 在时域用窗函数去截信号，对截下来的局部信号 作傅里叶变换，即在 t 时刻得该段信号的傅里叶变 换，不断地移动 t ，也即不断地移动窗函数的中心 位置，即可得到不同时刻的傅里叶变换，这些傅 里叶变换的集合，即是 。 ? STFT 可以看成是用基函数 来 代替傅里叶变换中的基函数 。 ) , ( t S x ? ?? ? ? ? j t e t m m ) ( ) ( , ? ? 0 200 400 600 800 1000 1200 1400 1600 1800 2000 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 ) ( ? x ) ( t m ? ? ? （ 1.1 ）式内积的结果即可实现对进行时－频 定位的功能。对 两边 做傅里叶变换，有 式中 和 是等效的频率变量 ?? ? ? ? j t e t m m ) ( ) ( , ? ? ? ? ? ? ? ? ? ? ?? ?? ? d e e t m M j j t ) ( ) ( , ? ? ? ? ? ? ? ? ? t d e t m e t j t j ) ( ) ( ) ( ? ? ? ? t j e M ) ( ) ( ? ? ? ? ? ? ? ? ? 该式指出，对 在时域加窗 ，引导在 频域对 加窗 。 ) ( ? x ) ( t m ? ? ) ( ? X ? ? ?? ? ) ( ), ( ) ( ), ( , 2 1 , ? ? ? ? ? ? t t M X m t x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? d e M X e t S t j t j x ) ( ) ( ) , ( 2 1 所以 ) ( ? ? ? M 信号谱 窗谱 ? ? ? ? ? ? ? d e M X t j ) ( 2 1 ) ( ) ( ? ? ? ? ? ? ? Parseval 定理 窗函数的中心和半径： 定义非平凡函数 称为一个窗函数，如果 也是属于 的，这个窗函数的中心定义为 ： 半径定义为： ? ? R L w 2 ? ? ? x xw ? ? R L 2 ? ? ? ? ? ? ? dx x w x w t 2 2 2