功率谱估计MATLAB实现.docVIP

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