以前谈到序列的实际长度可以通过零填充方法加入,使得最终增加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);

频谱如图:

余弦信号DFT频谱分析(继续)-LMLPHP

我们换三角窗:yd = [y.*triang(25)',zeros(1,999)];注意先加权再补零吧(事实上不是非常确定的说)。

频谱例如以下:

watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvaHVudGVybGV3/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center" alt="">

汉明窗:

watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvaHVudGVybGV3/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center" alt="">

尽管主瓣宽度加宽了,但咱能够继续加大N啊,所以不是问题。关键是如今频谱不泄露。

版权声明:本文博客原创文章,博客,未经同意,不得转载。

05-08 14:59