频谱泄露的分析及其处理方法.pdfVIP

频谱泄露的分析及其处理方法 在现代信号处理中，由于信号的频域分析比时域分析具有更加清 晰的物理概念和深刻含义，因而在信息技术领域中，FFT 运算 频谱 分析是一种常用的分析手段。对信号进行频谱分析首先需要通过信号 的傅里叶变换计算出信号对应的频谱函数，但是由于实际应用中接触 到的大量非周期连续信号x(t) 的频谱函数X(jω)是连续函数，利用计算 机对其进行频谱分析时往往需要对信号进行离散化处理以近似分析 相应的频谱。在离散化处理过程中由于被处理信号的有限记录长度 时域、频域的离散性往往造成在频谱分析中会出现一些特殊的效应， 例如混叠现象、泄漏现象以及栅栏现象，频谱泄漏就是这样出现的。 一．频谱泄漏的分析 所谓频谱泄露，就是信号频谱中各谱线之间相互影响，使 量结 果偏离实际值，同时在谱线两侧其他频率点上出现一些幅值较小的假 谱，导致频谱泄漏的原因是采样频率和信号频率的不同步，造成周期 采样信号的相位在始端 终端不连续。 设X(t)为实际信号，T0 为信号周期，f0=1/T0 为信号频率，Ts 为采 样周期，fs=1/Ts 为采样频率，L 是截取的周期数，N 是采样点数，L、N 均为正整数，X(t)经过长度为 LT0 的时间窗后得到离散序列 X(n) ，必 须满足采样频率和信号频率同步，即同步采样的要求: LT0/Ts=Nfs/f0。 当信号 X(t) 的频率 f0 是 fs/N 的整数倍时，这说明在处理长度 NT 内有信号的 K 个整周期。这时由 X(t)构成的以 NT 为周期的周期性信 号是连续的。当信号 X(t) 的频率 f0 不是 fs/N 的整数倍时，则在 NT 的 处理长度内，就不是恰好为信号周期的整数倍，有 X(t) 以 NT 为周期 进行周期延拓所得到的周期性信号就出现了不连续点，造成了频谱分 量从其正常频谱扩展开来，就这样形成了频谱泄漏现象。 在对信号做 FFT 分析时，如果采样频率固定不变，由于被采样信 号自身频率的微小变化以及干扰因素的影响，就会使数据窗记录的不 是整数个周期。从时域来说，这种情况在信号的周期延拓时就会导致 其边界点不连续，使信号附加了高频分量；从频域来说，由于 FFT 算法 只是对有限长度的信号进行变换，有限长度信号在时域相当于无限长 信号和矩形窗的乘积，也就是将这个无限长信号截短，对应频域的傅 里叶变换是实际信号傅里叶变换与矩形窗傅里叶变换的卷积。 当信号被截矩后的频谱不同于它以前的频谱。例如，对于频率为 fs 的正弦序列，它的频谱应该只是在 fs 处有离散谱。但是，在对它的 频谱做了截短后，结果使信号的频谱不只是在 fs 处有离散谱，而是在 以fs 为中心的频带范围内都有谱线出现，它们可以理解为是从 fs 频率 上泄漏出去的，这种现象就是频谱泄漏。泄漏现象对功率谱估计及正 弦分量的检 均带来有害的影响，因为弱信号的主瓣很容易被强信号 泄漏到邻近的副瓣所淹没及畸变的，从而造成谱的模糊与失真。 通 过 LABVIEW 信号处理实验室可以看到当边界点不连续时出现的频 谱泄漏的情况如下图 1 所示： 图 1 信号边界点不连续时 接下来举例说明以上的情况。假设连续信号 X(t) 的周期为 T，现 在对它进行采样，采样时间为 t ，采样 N 个点，那么 T=N*t,因为 f(t)的 频率 f0=2*pi/T,同时又有 T=N*t、fs=2*pi/t,则有 f0=2*pi/N*t=fs/N。接 着我们假设对一个周期采样 N=32 个点，则有 f0=fs/N ；当对一个周期 采样 N1=64 个点，那么 N1=2*N ，有 f0=fs/N=fs/N1/2,即 f0=2*fs/N1 ；同 理当 N2=128 ，f0=4*fs/N2… 也就是说如果采样的不是整数倍的信号周期，那么这 32 个点、64 个点、128 个点就不是在一个整周期内采到的，那么上面的等式也 就不成立了，因此也就发生了频谱泄漏。如果原始信号的频谱成份与 FFT 中的谱线完全一致，这种情况下采样数据的长度为信号周期的整 数倍，频谱中只有主瓣。没有出现旁瓣的原因是旁瓣正处在窗函数主 瓣两侧采样频率间隔处的零分量点。如果时间序列的