功率譜估計

自然界出現的許多現象都可以在統計平均意義上很好的表現出來。例如，氣象學中的氣溫與氣壓的變動等，均可以以統計的方式表示為隨機過程。在電阻器和電子裝置中生成的熱噪音電壓，也是被抽象為隨機過程模型的物理訊號的例子。由於這些訊號為隨機訊號，我們必須採用一種統計觀點來處理隨機訊號的平均特徵。特別的是隨機訊號的自相關函式很適合用於代表時域中的隨機訊號，並且自相關函式的傅立葉轉換可生成功率譜密度，也可提供時域到頻域的轉換。基於有限長訊號觀察的功率譜估計.對於有限時間長度的訊號，資料序列的有限記綠長度是對功率譜估計(power spectrum estimation)的主要限制。當處理統計平穩訊號時，資料記錄越長，可從資料提取的訊號估計就越好。另一方面若訊號統計是非平穩的，我們不能選擇任意長度記綠對譜進行估計，在這種情況下我們選擇的資料記綠長度是由訊號統計上的時變速度決定的。最後我們要選擇盡可能短且能解析不同訊號分量譜特徵的資料記綠，所表現的這些訊號分量具有相近空間譜。對基於有限長度資料記綠的古典功率譜估計方法，所面臨的問題之一是我們試圖要估計出的譜會有失真。這一問題無論在確定性訊號的譜計算方面還是在隨機訊號的功率譜估計方面都會遇到，既然很容易觀察到有限長度資料記錄對確定性訊號的效應，我們就先考察確定性訊號的情況然後再討論隨機訊號及其功率譜估計。能量密度譜計算.考慮有限長度資料序列確定性訊號的譜計算，可功率譜密度。隨機訊號的自相關和功率譜估計:週期圖.有限能量訊號可進行傅立葉轉換，並在譜域用它們的能量密度譜來表現。另一方面代表為平穩隨機過程的重要型別訊號不具有有限能量，因此不能進行傅立葉轉換。這類訊號具有有限平均功率，因此表現為功率譜密度。如果formula_1是一個平穩隨機過程，它的自相關函式是:formula_2其中formula_3代表統計平均。於是由維納-辛欽定理(Wiener–Khinchin theorem)，平穩隨機過程的功率譜是自相關函式的傅立葉轉換，即:formula_4實際上我們處理隨機過程的單個實現並從中估計該過程的功率密度譜。由於不知到真實的自相關係數formula_5，導致我們不能按上式計算傅立葉轉換來得到formula_6。從隨機過程的單個實現，可以計算時間平均自相關函式:formula_7其中formula_8是觀察期間。如果平穩隨機過程的一階和二階矩(均值和自相關函式)是各態歷經的，那麼formula_9這一關係證實了時間平均自相關函式formula_10可用做對統計自相關函式formula_5的估計。更進一步formula_10的傅立葉轉換提供了對功率密度譜的估計formula_13,即:formula_14formula_15formula_16實際功率密度譜是formula_13在極限formula_18時的期望值:formula_19formula_20我們將從隨機過程的單個實現樣本考慮功率密度譜估計。假設formula_21以formula_22取樣，其中B是隨機過程功率密度譜包含的最高頻率。因此透過對formula_21取樣，我們得到有限長序列formula_24。從這些樣本我們計算時間平均自相關序列:formula_25並且對於m的負值，我們有formula_26。於是我們計算傅立葉轉換:formula_27上上式中的歸一化因子formula_28導致了均值估計:formula_29其中formula_30是formula_31的真實(統計的)自相關序列。因此formula_32是自相關函式formula_30的無偏差估計。估計formula_32的方差近似為:formula_35這是由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_47formula_48其偏移為formula_49，因為其均值是:formula_50然而該估計具有較小的方差，近似為:formula_51注意到formula_30是漸進無偏的，即formula_53。並且當formula_39時其方差收斂於零。因此估計formula_30也是formula_56的一致估計。在處理功率譜估計問題時，我們將使用由式formula_47和formula_48給出的估計formula_30。相應的功率譜密度是:formula_60我們把得到的formula_30代入上式，估計formula_62可以表示為:formula_63其中formula_64是樣本序列formula_31的傅立葉轉換。這種常見形式的功率密度譜估計稱為週期圖。它最初是由Schuster於1898年引入用來檢測和測量存在於資料中的"隱藏週期"。從式formula_60可推出週期圖估計formula_62的均值是formula_68formula_69對上兩式的解釋是，估記譜的均值是窗自相關函式的傅立葉轉換formula_70，其中窗函式是(三角形的)巴特利特(Bartlett)窗。因此估計譜的均值是:formula_71其中formula_72是巴特利特窗的譜特徵。上式說明瞭估計譜的均值是真實功率譜密度formula_73與巴特利特窗傅立葉轉換formula_72的旋積。結果估計譜的均值是真實譜的平滑版，受損於有限資料點導致的相同的譜洩漏。注意到估計譜是漸進無偏的，即:formula_75然而一般來說當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個資料段:formula_95對於每一段，可計算週期圖:formula_96最後對K段的週期圖進行平均得到Bartlett功率譜估計:formula_97該估計的統計特性很容易得到。首先均值是formula_98。由式formula_68和formula_69及formula_71可得到單個週期圖的期望值是:formula_102其中formula_103是巴特利特窗formula_104的頻率特性。從formula_105的公式注意到，現在真實譜與巴特利特窗的頻率特性formula_72有關。資料長度從formula_43點減小到formula_108，導致視窗的譜寬度增加formula_92因子。結果頻率解析度得到減小因子formula_92。解析度降低的結果使得方差減小。Bartlett估計的方差是:formula_111如果我們利用式formula_78和上式，得到:formula_113因此，Bartlett功率譜估計的方差減小因子為formula_92。