功率谱估计

自然界出现的许多现象都可以在统计平均意义上很好的表现出来。例如，气象学中的气温与气压的变动等，均可以以统计的方式表示为随机过程。在电阻器和电子设备中生成的热噪音电压，也是被抽象为随机过程模型的物理讯号的例子。由于这些讯号为随机讯号，我们必须采用一种统计观点来处理随机讯号的平均特征。特别的是随机讯号的自相关函数很适合用于代表时域中的随机讯号，并且自相关函数的傅立叶转换可生成功率谱密度，也可提供时域到频域的转换。
基于有限长讯号观察的功率谱估计.
对于有限时间长度的讯号，数据序列的有限记绿长度是对功率谱估计(power spectrum estimation)的主要限制。当处理统计平稳讯号时，数据记录越长，可从数据提取的讯号估计就越好。另一方面若讯号统计是非平稳的，我们不能选择任意长度记绿对谱进行估计，在这种情况下我们选择的数据记绿长度是由讯号统计上的时变速度决定的。最后我们要选择尽可能短且能解析不同讯号分量谱特征的数据记绿，所表现的这些讯号分量具有相近空间谱。对基于有限长度数据记绿的古典功率谱估计方法，所面临的问题之一是我们试图要估计出的谱会有失真。这一问题无论在确定性讯号的谱计算方面还是在随机讯号的功率谱估计方面都会遇到，既然很容易观察到有限长度数据记录对确定性讯号的效应，我们就先考察确定性讯号的情况然后再讨论随机讯号及其功率谱估计。
能量密度谱计算.
考虑有限长度数据序列确定性讯号的谱计算，可参见:功率谱密度。
随机讯号的自相关和功率谱估计:周期图.
有限能量讯号可进行傅立叶转换，并在谱域用它们的能量密度谱来表现。另一方面代表为平稳随机过程的重要类型讯号不具有有限能量，因此不能进行傅立叶转换。这类讯号具有有限平均功率，因此表现为功率谱密度。如果formula_1是一个平稳随机过程，它的自相关函数是:
其中formula_3代表统计平均。于是由维纳-辛钦定理(Wiener–Khinchin theorem)，平稳随机过程的功率谱是自相关函数的傅立叶转换，即:
实际上我们处理随机过程的单个实现并从中估计该过程的功率密度谱。由于不知到真实的自相关系数formula_5，导致我们不能按上式计算傅立叶转换来得到formula_6。从随机过程的单个实现，可以计算时间平均自相关函数:
其中formula_8是观察期间。如果平稳随机过程的一阶和二阶矩(均值和自相关函数)是各态历经的，那么
这一关系证实了时间平均自相关函数formula_10可用做对统计自相关函数formula_5的估计。更进一步formula_10的傅立叶转换提供了对功率密度谱的估计formula_13,即:
实际功率密度谱是formula_13在极限formula_18时的期望值:
我们将从随机过程的单个实现样本考虑功率密度谱估计。假设formula_21以formula_22取样，其中B是随机过程功率密度谱包含的最高频率。因此通过对formula_21取样，我们得到有限长序列formula_24。从这些样本我们计算时间平均自相关序列:
并且对于m的负值，我们有formula_26。于是我们计算傅立叶转换:
上上式中的归一化因子formula_28导致了均值估计:
其中formula_30是formula_31的真实(统计的)自相关序列。因此formula_32是自相关函数formula_30的无偏差估计。估计formula_32的方差近似为:
这是由Jenkins和Watts于1968年给出的结果，显然formula_36有formula_37。
因为formula_38,并且当formula_39时估计的方差收敛于零，所以估计formula_32是相容的。对于较大值的滞后参数formula_41，特别当formula_41逼近于formula_43时，由formula_25给出的估计formula_32具有较大方差。这是由于很少的数据点数进入大的滞后情况下的估计。作为式formula_25的备用方法，我们使用如下的估计:
其偏移为formula_49，因为其均值是:
然而该估计具有较小的方差，近似为:
注意到formula_30是渐进无偏的，即formula_53。并且当formula_39时其方差收敛于零。因此估计formula_30也是formula_56的一致估计。在处理功率谱估计问题时，我们将使用由式formula_47和formula_48给出的估计formula_30。相应的功率谱密度是:
我们把得到的formula_30代入上式，估计formula_62可以表示为:
其中formula_64是样本序列formula_31的傅立叶转换。这种常见形式的功率密度谱估计称为周期图。它最初是由Schuster于1898年引入用来检测和测量存在于数据中的"隐藏周期"。从式formula_60可推出周期图估计formula_62的均值是
对上两式的解释是，估记谱的均值是窗自相关函数的傅立叶转换formula_70，其中窗函数是(三角形的)巴特利特(Bartlett)窗。因此估计谱的均值是:
其中formula_72是巴特利特窗的谱特征。上式说明了估计谱的均值是真实功率谱密度formula_73与巴特利特窗傅立叶转换formula_72的旋积。结果估计谱的均值是真实谱的平滑版，受损于有限数据点导致的相同的谱泄漏。注意到估计谱是渐进无偏的，即:
然而一般来说当formula_39时估计formula_62的方差不会衰减到零。例如当数据序列是一个随机过程时，方差是formula_78,当formula_39时，极限为formula_80。 
因此我们认为周期图不是真实谱密度的一致估计(即不收敛于真正的功率谱密度)。概括来说估计的自相关formula_30是真实自相关函数formula_56的一致估计。然而它的傅立叶转换formula_62即周期图不是真实功率谱密度的一致估计。我们注意到formula_62是formula_85的渐进无偏差估计，但是对于一个有限长序列，从式formula_71得到的formula_62均值包含了偏移，说明真实功率谱密度产生了失真。于是估计谱受损于巴特利特窗的平滑效应和具体的泄漏，平滑和泄漏最终限制了我们准确分析紧密谱的能力。
功率谱估计的非参数化方法.
非参数方法是由Bartlett(1948)、Blackman和Tukey(1958)及Welch(1967)开发的方法，这些方法没有假定数据是如何生成的，故被称作非参数方法。既然这些估计全部基于数据的有限记录，这些方法的频率分辨率最好等于长度为formula_43的矩形窗的宽度，也就是在formula_89近似为formula_90。这将会更加精确地指定具体方法的频率分辨率。为了减小谱估计的方差，非参数方法会降低频率分辨率。
Bartlett方法:平均周期图.
减小周期图方差的Bartlett方法包含了三个步骤。首先formula_43点序列被划分为formula_92个不重叠段，每段的长度为formula_93。这样就生成了formula_92个数据段:
对于每一段，可计算周期图:
最后对K段的周期图进行平均得到Bartlett功率谱估计:
该估计的统计特性很容易得到。首先均值是formula_98。由式formula_68和formula_69及formula_71可得到单个周期图的期望值是:
其中formula_103是巴特利特窗formula_104的频率特性。
从formula_105的公式注意到，现在真实谱与巴特利特窗的频率特性formula_72有关。数据长度从formula_43点减小到formula_108，导致窗口的谱宽度增加formula_92因子。结果频率分辨率得到减小因子formula_92。分辨率降低的结果使得方差减小。Bartlett估计的方差是:
如果我们利用式formula_78和上式，得到:
因此，Bartlett功率谱估计的方差减小因子为formula_92。