新MATLAB辅助现代工程数字信号处理 第二版 教学课件 李益华 第6 10章_ 第8章.pptVIP

第8章 非平稳信号分析与处理 　　　　　8.1 短时傅里叶变换8.1.1 时域窗法　　对于时变信号，需了解它的局部傅里叶变换，最简单的方法就是对信号乘上一个滑动窗来研究它的傅里叶变换。　　设离散信号为x(n)，在n0时刻局部的傅里叶变换为 式中，w(n)为窗函数，其长度为Nw，通常可选为中心对称。定义n时刻的短时段信号为　　　　　　　　　　xn(m)=x(m)w(n－m)　　　　　　(8.3)　　因此，时域窗法就是用一个中心对称的滑动窗函数截取观测信号，对不同时刻的短时段信号进行傅里叶变换，最后得到各段信号构成的时变谱阵。短时傅里叶变换对信号分析的示意图如图8.1所示。 图 8.1 短时傅里叶变换示意图 　　通常对于短时傅里叶变换X(n，ω)进行时域抽取，值得注意的是抽取的间隔L必须满足L≤Nw/2，这样才能保证信号不丢失，即可以从时域信号恢复原信号。抽取后的短时傅里叶变换为 8.1.2 频域窗法　　如果用W(ω)表示窗函数w(n)的傅里叶变换，则信号的STFT还可表示为 图 8.2 短时傅里叶变换频域窗法原理图 　　由图可以看出，对每个频率值，STFT是序列x(n)e－jω0n与窗函数w(n)的卷积，从频率解释就是先将信号频谱X(ω)下变频为X(ω+ω0)，然后通过滤波器滤波。　　因此，频域窗法就是将信号通过具有不同中心频率、带宽为Δω的窄带滤波器W(ω)进行滤波。为使其具有较高的频率分辨率， W(ω)带外衰减越大越好。 8.1.3 窗函数的选取　　1. 不同窗的选择　　矩形窗具有最窄的主瓣宽度，但同时也具有最大的旁瓣峰值和最慢的旁瓣峰值衰减速度，矩形窗存在的边界效应会造成频谱泄漏，大大影响分析性能; 汉宁窗和汉明窗具有较小的旁瓣峰值和较大的旁瓣峰值衰减速度，但主瓣稍宽于矩形窗，在计算机短时傅里叶变换时，将大大减少频谱泄漏，具有较优的分析性能。　　理想的最佳分析窗函数应具有最窄的主瓣宽度、最小的旁瓣峰值和最大的衰减速度，因此，应在不同条件下通过比较选择最合适的分析窗函数。 　　2. 不同的窗长度选择　　在选定了分析窗函数之后，合理地选择分析窗长度也有利于分析精度的提高。若窗口过长，则信号不能近似看做平稳信号，傅里叶变换将失去作用，也会加大运算量; 若窗口过小，又会丢失信息。因此，慎重选取窗函数的长度也是很有必要的。　　MATLAB信号处理工具箱提供了函数specgram，用于计算短时傅里叶变换的频谱图，其调用格式为 　　B=specgram(a) 　　B=specgram(a， nfft， Fs) 　　[B， f]= specgram(a， nfft， Fs) 　　[B， f， t]= specgram(a， nfft， Fs) 　　　B=specgram(a， nfft， Fs， windows) 　　 B=specgram(a， nfft， Fs， windows， noverlap) 　　 specgram(a) 　　 B=specgram(a， f， Fs， windows， noverlap)其中，B=specgram(a)为返回信号向量a的谱图，该调用缺省时，参数配置如下: 　　(1) FFT计算点数nfft＝min(256，length(a)); 　　(2) 采样频率Fs＝2; 　　(3) 窗函数windows＝hanning; 　　(4) 分段重叠的点数noverlap＝length(windows)/2; 　　当a为实数时，specgram只在正频域内计算离散傅里叶变换。当a为偶数时，B的行数为(nfft/2+1); 当a为奇数时，B的行数为nfft/2，B的列数为fix((n－noverlap)/(length (windows)－noverlap)。当a为复数时，specgram函数要计算正、负频域的离散傅里叶变换，这时B为nfft行的复矩阵，从B中第一点开始，时间随列数线性增加，而频率从0开始随行数线性下降。　　B＝specgram(a，nfft)采用指定的nfft点FFT，点数为2的幂次时，运算速度较快。　　[B， f]＝specgram(a， nfft， Fs)除了得到谱图B外，还返回离散傅里叶变换向量f，Fs决定了f的频率分辨率，并不影响B的结果。 　　[B， f， t]= specgram(a， nfft， Fs)附加地得到记录时域窗函数与a相交时刻的向量t。　　B=specgram(a， nfft， Fs， windows)中windows指定所采用的窗函数和窗的宽度。窗函数的宽度应该小于或等于nfft，否则补零。　　B=specgram(a， nfft， Fs， wi