频谱泄漏

版权声明:欢迎转载欢迎学习 https://blog.csdn.net/dreamandxiaochouyu/article/details/45393319

使用DFT进行信号频谱分析导致频谱泄露问题的探究

 

一,引论

    这篇文章里,我将要叙述一个在利用DFT去估计信号频谱时很容易引发的一个问题--频谱泄露。这个名词算是比较专业的词。之所以前面我用了两篇文章去讲时域窗和DTFT/DFT就是为这个问题做铺垫。我决定不在开始就讲他的概念是因为这个概念直接讲出来会很困难,也难以被接受,因为它确实需要模型去论证。初学者如果没有信号理论,DFT,窗的基本常识,直接看这个可能会有点困难,所以建议先打好基础。

 

二,有限点数DFT的实现过程

    通常的DFT用来分析一小段长度有限的离散信号,只是大多数信号肯定不是一小段这么简单,所以欲要进行DFT分析,先要将长信号截断,如何截?很简单嘛,在时域乘以一个窗不就完成了截断嘛。

    我们知道,时域乘积对应频域卷积,因此必有:


    基于频域卷积这个事实,我们来分析一下最简单的情况,在时域找一个单频信号x(n)=sin(2pi*w0*n),窗就选矩形窗来做个实验。

    首先声明,这个图解过程并非直接用的DFT,而是用DTFT来展示,这样做的好处是连续的谱线看起来更直观,事实上当我们对最后的结果Y(w)进行抽样就得到了DFT的离散值。这样稍微的变动并不会影响分析,反而会使分析变得简单。



     这个图很具体的展示出了经过频域运算后原信号x(n)的频谱X(w)变成新的截断信号的频谱Y(w)的过程,从图中可以很明显的感觉到Y(w)产生的变化:

    (1),原来的X(w)在w0处单一的频率谱线经过窗之后形成的Y(w)变成了被展宽的带状区域,在w0处明显变胖了,毫无疑问这个就是矩形窗固有的主瓣。
    (2),因为矩形窗旁瓣的存在,使得运算后的Y(w)引入了很多拖尾(浪花)。
   

    我们似乎能感觉到:窗的引入会使得原信号频谱产生一种可能会影响我们分析DFT的反面功效。试想一下,如果我们选取了一个窗,他的主瓣比较胖,而我们待分析的原始信号有两个相隔很近的频率点,例如:x(n)=sin(0.3pi*n)+sin(0.302pi*n),请问最终的DFT图谱上我们还能看到两个峰值点吗?答案肯定是NO,因为那个太胖的主瓣完全包容了0.3和0.302这两个频率点的微小距离,还能如何把他们区分开来?

    再试想另外一种情况,如果我们选了一个窗,它的旁瓣衰减不够猛,旁瓣的高度比较高,正好我们分析的信号其中有一个频率点的功率比较小,例如:x(n)=sin(0.3pi*n)+0.001cos(0.45pi*n),那么请问0.45pi那个频率点的信号会在最终的DFT图谱上展示出来吗?答案必然是NO,因为如果那个旁瓣比较高,那么他完全就把你这个小功率的0.5pi这个点给完全淹没掉。

 

    以上这段解释是我发挥想象来讲述的,事实上这个需要理论论证或者仿真,总之他揭示了窗的反面功效,从侧面告诉我么一个道理--截断信号,进行DFT分析之前,针对不同类型的信号,时域窗类型的选择至关重要。如果我们待截断的信号包含一些相隔很近的频率点,那么必须选取的窗的主瓣不能太胖。如果我们待截断的信号包含一些功率很小的频率点,那么我们选取的窗的旁瓣衰减必须要足够大。

 

    下面我附上信号x(n)=0.01cos(0.45pi*n)+sin(0.3pi*n)+sin(0.302pi*n)通过矩形窗截取,和通过汉明窗截取之后再进行DFT计算之后的异同。


                                                                
    可以看到,矩形窗主瓣过宽,导致有微小偏差0.002pi的两个正弦信号显示不出来,并且由于旁瓣衰减不够大,所以那个小功率的余弦信号也未能展现。

 

 

 


    选取汉明窗之后,通过增加窗口点数,使得它主瓣变窄窄,所以具有0.002pi微小频率偏差的两个正弦信号放大之后显形,同时该窗旁瓣衰减的给力,使得小功率点的正弦信号也浮出水面。

 

 

三,频谱泄露概念的引出

    之所以现在才说什么是频谱泄露也是基于上面的论证过程,其实看到这里想必大家已经知道了究竟什么是频谱泄露了,那就是----采用窗的截取方式去截断信号之后进行DFT分析频谱,可能会因为窗的选取不当或者窗的长度N不够,导致最终DFT频谱图中丢失了原信号应该有的频率成分。

 

 

四,结论

    欲采用窗的截取方式进行DFT分析离散信号频谱,窗的选取是有原则的。一是对于原有信号存在微小频率偏差的情况,窗的主瓣宽度必须小于信号频率偏差的大小,当然,这可以通过改变窗的长度N去实现。二是,对于原信号中存在小功率频率点,必须考虑窗的旁瓣衰减要给力才行,当然这就要考验窗的本身性能了。本博客在前面已经给出了常见的窗的频域参数,这就是在分析信号的过程中很有必要查阅的一部分表格之一。

展开阅读全文

没有更多推荐了,返回首页