功率譜估計
自然界出現的許多現象都可以在統計平均意義上很好的表現出來。例如，氣象學中的氣溫與氣壓的變動等，均可以以統計的方式表示為隨機過程。在電阻器和電子設備中生成的熱噪音電壓，也是被抽象為隨機過程模型的物理訊號的例子。由於這些訊號為隨機訊號，我們必須採用一種統計觀點來處理隨機訊號的平均特徵。特別的是隨機訊號的自相關函數很適合用於代表時域中的隨機訊號，並且自相關函數的傅立葉轉換可生成功率譜密度，也可提供時域到頻域的轉換。
對於有限時間長度的訊號，數據序列的有限記綠長度是對功率譜估計(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方法包含了三個步驟。首先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。