功率谱估计性能分析及其MATLAB实现
一、经典功率谱估计分类简介
1.间接法
根据维纳-辛钦定理,1958年Blackman和Turkey给出了这一方法的具体实现,即先由N个观察值,估计出自相关函数,求自相关函数傅里叶变换,以此变换结果作为对功率谱的估计。
2.直接法
直接法功率谱估计是间接法功率谱估计的一个特例,又称为周期图法,它是把随机信号的N 个观察值直接进行傅里叶变换,得到,然后取其幅值的平方,再除以N,作为对功率谱的估计。
3.改进的周期图法
将N点的观察值分成L个数据段,每段的数据为M,然后计算L个数据段的周期图的平均,作为功率谱的估计,以此来改善用N点观察数据直接计算的周期图的方差特性。根据分段方法的不同,又可以分为Welch法和Bartlett法。
Welch法
所分的数据段可以互相重叠,选用的数据窗可以是任意窗。
Bartlett法
所分的数据段互不重叠,选用的数据窗是矩形窗。
二、经典功率谱估计的性能比较
1.仿真结果
为了比较经典功率谱估计的性能,本文采用的信号是高斯白噪声加两个正弦信号,采样率Fs=1000Hz,两个正弦信号的频率分别为f1=200Hz,f2=210Hz。所用数据长度N=400.
仿真结果如下: