部分内容参考
[1] 基于Welch算法的经典功率谱估计的MATLAB分析(http://www.doc88.com/p-3734310145747.html)
功率谱估计是利用有限长的数据估计信号的功率谱,分为经典谱估计和现代谱估计。下面我们从经典谱估计中的周期图法入手,进一步讨论其改进Welch算法。
1、周期图法
周期图法有傅立叶变换得到的,因此也称为直接法。若我们有离散的 N N N点随机序列 x N ( n ) , n = 0 , 1 , … , N − 1 x_N(n),n=0,1,\ldots,N-1 xN(n),n=0,1,…,N−1,则可以得到其 N N N点离散傅立叶变换(DFT) X N ( f ) X_N(f) XN(f),进一步可以得到对 x ( n ) x(n) x(n)真实功率谱密度 P X ( f ) P_X(f) PX(f)的估计:
(1) P ^ X ( f ) = 1 N ∣ X N ( f ) ∣ 2 = 1 N ∣ ∑ n = 0 N − 1 x ( n ) e − j 2 π f n ∣ 2 . \tag{1} \hat P_X(f)=\frac{1}{N}|X_N(f)|^2=\frac{1}{N}{\Large |}\sum_{n=0}^{N-1}x(n)e^{-j2\pi fn}{\Large |}^2. P^X(f)=