STFT短时傅里叶变换

短时傅里叶变换唐娜主要内容1.短时傅里叶变换简介2.测不准定理3.短时傅里叶变换缺陷背景:傅里叶变换的缺陷1）不适用于非平稳信号的处理2）没有局域性3）时域与频域的分割提出与基本思想鉴于傅里叶变换的缺陷提出了窗函数的概念，提出一个灵活可变的时间—频率窗，使得在这个窗内能够体现频率的信息，这种信号分析方法称为时间—频率分析。而窗固定的时间—频率分析方法即为短时傅里叶变换。短时傅里叶变换（STFT，short-timeFouriertransform）。其主要思想是将信号加窗，将加窗后的信号再进行傅里叶变换，加窗后使得变换为时间r附近的很小时间上的局部谱，窗函数可以根据r的位置变化在整个时间轴上平移，利用窗函数可以得到任意位置附近的时间段频谱实现时间局域化STFT定义1946年，Gabor就提出了STFT，给定一信号，其STFT定义为：公式涵义：在时域用窗函数去截信号，对截下来的局部信号作傅立叶变换，即在t时刻得该段信号得傅立叶变换，不断地移动t，也即不断地移动窗函数的中心位置，即可得到不同时刻的傅立叶变换，这些傅立叶变换的集合，即是STFT可以看成是用基函数来代替傅里叶变换中的基函数jjxetgxdetgxtSTFT)(),()()(),(),(tSTFTxjtetgg)()(,窗函数（1.1）0200400600800100012001400160018002000-0.4-0.3-0.2-0.100.10.20.3（1.1）式内积的结果即可实现对进行时－频定位的功能。对两边做傅里叶变换，有式中和是等效的频率变量jtetgg)()(,deetgGjjt)()(,tdetgetjtj)()()(tjeG)()(该式指出，对在时域加窗，引导出在频域对加窗。)(x)(tg)(X)(),()(),(,21,ttGXgtxdeGXetSTFTtjtjx)()(),(21所以)(G信号谱窗谱deGXtj)(21)()(窗函数的中心和半径：定义非平凡函数称为一个窗函数如果也是属于的，一个窗函数的中心定义为：半径定义为：RL2xxRL2dxxxt222*12122*2})()({1dxxtx这样我们可以认为函数集中定义在以为中心，为半径长为的区间上，取此区间为的有效区间是合适的。对函数其中心为半径为，对的Fourier变换设其中心为半径为矩形，称为函数的时-频窗。该窗的面积为w*t*,*tt2w*t*ww]*,*[tt×**,w221.短时傅里叶变换简介2.测不准定理3.短时傅里叶变换缺陷以Gabor函数为例，令Gabor数函为窗函数，已知Gabor函数的表达式如下：以Gabor函数为窗函数的STFT称为Gabor变换，其定义为则Gabor函数的中心和半径为：ataeatg4221dttgtfesati),(dxxgxgt222||1*根据窗函数半径公式，可知道Gabor窗函数的半径为：因为窗的面积为，对于Gabor窗函数说因为为的傅里叶变换，则对于Gabor函数就要求出g的傅里叶变换，再代入上式得出则经计算得adxxgtxgbabagba2122*2}{122baggbagbagaabag2141则有：可以证明，不论采用何种函数作为窗函数，其时间窗和频率窗宽度的乘积的最小值都是2，这就是测不准原理，此定理告诉我们，不可能在时间和频率两个空间同时以任意精度逼近被测信号，因此就必须在信号的分析上对时间或者频率的精度做取舍。当利用STFT时，若我们希望能得到好的时－频分辨率，或好的时－频定位，应选取时宽、带宽都比较窄的窗函数，遗憾的是，由于受不定原理的限制，我们无法做到使同时为最小。221422aabagag)(g当我们对信号作时－频分析时，一般，对快变的信号，我们希望它有好的时间分辨率以观察其快变部分（如尖脉冲等），即观察的时间宽度要小，受时宽－带宽积的影响，这样，对该信号频域的分辨率必定要下降。由于快变信号对应的是高频信号，因此对这一类信号，我们希望有好的时间分辨率，但同时就要降低高频的分辨率。反之，对慢变信号，由于它对应的是低频信号，所以我们希望在低频处有好的频率分辨率，但不可避免的要降低时域的分辨率。主要内容1.短时傅里叶变换简介2.测不准定理3.短时傅里叶变换缺陷前面推导了测不准定理，知道STFT不具备自动调节能力窗函数选定形状不会发生改变时频窗在时间轴频率轴方向上的宽度确定时频分辨率确定不随时间、频率的变化而变化从上面的分析我知道，如果要改变分辨率，则需要重新选择窗函数。因为受到不确定准则的限制，时频窗的面积不小于2，故不能兼顾频率与时间分辨率的需求，这也就从另一个侧面说明了短时傅里叶变换窗函数的时间与频率分辨率不能同时达到最优，我们对时间分辨率和频率分辨率只能取一个折中，一个提高了，另一个就必然要降低，反之亦然。短时傅里叶变换用来分析分段平稳信号或者近似平稳信号犹可，但是对于非平稳信号，对频率分辨率和时间分辨率的要求是要按照一定的规律变化的，但短时傅里叶变换的函数一旦选定时频分辨率是确定不随时间、频率的变化而变化。加窗实例[]xn[]wn[]yn频谱频谱加窗后频谱产生失真正弦序列接收下行转换接收下行转换接收下行转换A/DA/DA/DW1W2Wp采用自适应算法调整系数解调控制部分▽▽▽天线阵列波束形成网络id期望信号散射媒介基站移动台)(ty)(ts