Welch功率谱估计中窗函数选择与算法分析.docVIP

Welch功率谱估计中窗函数选择与算法分析 　　摘 要： 以Matlab为平台对Welch算法中窗函数的选择与使用进行了分析。从频谱分辨率、噪声水平等方面对窗函数特性进行分析，其中矩形窗与凯撒窗频谱分辨率较高，但信号频率附近噪声水平较高。分析认为矩形窗与凯撒窗适用于高信号信噪比的高精度频谱估计。汉宁窗与切比雪夫窗对频谱泄漏抑制效果较好，信号频率附近噪声水平未受影响，但其分辨率相对较低。分析认为汉宁窗与切比雪夫窗适用于信噪比较低的信号频率的大致估计。Welch功率谱估计算法针对短信号使用的研究认为，直接对短信号进行Welch功率谱估计不仅频谱分辨率低，信号也容易被噪声干扰导致探测失败。短信号功率谱的估计无法直接通过Welch功率谱估计法来实现。 　　关键词： Welch功率谱估计； 窗函数； 频谱分辨率； 频谱泄漏； 信噪比 　　中图分类号：TP391.9 文献标志码：A 文章编号：1006-8228（2018）02-01-03 　　Abstract： The selection and use of window function in Welch algorithm are analyzed with Matlab as the platform. The characteristics of window function are analyzed from spectrum resolution and noise level. The spectral resolution of Rectangular window and Kaiser window is higher， but the noise level near the signal frequency is higher. It is considered that Rectangular window and Kaiser window are suitable for high precision spectrum estimation of signals with high signal to noise ratio. Hanning window and Chebyshev window have better effect on spectrum leakage suppression， and the noise level near the signal frequency is not affected， but the frequency resolution is relatively low， suitable for general frequency estimation of signals with low signal to noise ratio. The study of using Welch algorithm in frequency estimation of short signals shows that it is not only the low spectrum resolution， but also the signal is easily disturbed by the noise， which leads to the failure of the detection. 　　Key words： Welch power spectrum estimation； window function； spectral resolution； spectrum leakage； signal to noise ratio 　　0 引言 　　功率?V估计技术是通过信号的相关性，对接收信号功率随频率的变化关系进行估计的一种频谱估计方法，其基本功能是实现宽带噪声中窄带信号的检测。功率谱估计可以实现对信号的分析与识别。周期图法利用相关函数与傅里叶变换的卷积关系求功率谱，由于周期图法不满足一致性估计条件，分辨率与方差难以同时取得较高的参数。Welch算法是在周期图法基础上进行修正的一种比较简单，便于计算的方法。通过对数据进行分段和加窗，可以在保证分辨率不受影响的前提下尽量减小谱估计的方差。在算法应用中，窗函数的选择将直接影响谱估计的结果。 　　1 周期图及Welch算法 　　1.1 周期图法 　　周期图法是将观测得到的有限个样本数据x（n）={x（0），x（1），…，x（N-1）}视为能量有限信号，利用FFT直接计算出其傅里叶变换，然后取其复制的平方并除以数据长度N，作为功率谱估计： 　　⑴ 　　其中X（k）为有限长序列x（n）（n=0，1，……，N-1）的傅里叶变换。 　　通过周