上一篇讲到实际可以通过补零的方法增加序列长度,使得最后增大N来增加视在分辨率。但是它无法解决频率泄漏的问题。频率泄漏的根本原因在于窗的选择。
由于矩形窗突然截断,频谱旁瓣相对幅度过大,导致泄漏的分量很多。因此和FIR滤波器的思路一样,我们想到了其他窗。
接上次的例子,矩形窗:
ts = 0.01;
n = 0:24;
y = [sin(2*pi*20*n*ts),zeros(1,999)];
xk = abs(fft(y,1024));
stem(xk);
频谱如图:
我们换三角窗:yd = [y.*triang(25)‘,zeros(1,999)];注意先加权再补零吧(其实不是很确定的说)。频谱如下:
汉明窗:
虽然主瓣宽度加宽了,但咱可以继续加大N啊,所以不是问题,关键是现在频谱不泄漏了。
原文地址:http://blog.csdn.net/hunterlew/article/details/40123269