ar模型的功率譜估計burg算法的分析與仿真論文

《ar模型的功率譜估計burg算法的分析與仿真論文》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。

1、AR模型的功率譜估計BURG算法的分析與仿真錢平（信號與信息處理S101904010）一．引言現(xiàn)代譜估計法主要以隨機過程的參數(shù)模型為基礎(chǔ)，也可以稱其為參數(shù)模型方法或簡稱模型方法?，F(xiàn)代譜估計技術(shù)的研究和應(yīng)用主要起始于20世紀60年代，在分辨率的可靠性和濾波性能方面有較大進步。目前，現(xiàn)代譜估計研究側(cè)重于一維譜分析，其他如多維譜估計、多通道譜估計、高階譜估計等的研究正在興起，特別是雙譜和三譜估計的研究受到重視，人們希望這些新方法能在提取信息、估計相位和描述非線性等方面獲得更多的應(yīng)用。現(xiàn)代譜估計從方法上大致可分為參數(shù)模型譜估計和非參數(shù)模型譜估計兩種?；趨?shù)建摸的功率

2、譜估計是現(xiàn)代功率譜估計的重要內(nèi)容，其目的就是為了改善功率譜估計的頻率分辨率，它主要包括AR模型、MA模型、ARMA模型，其中基于AR模型的功率譜估計是現(xiàn)代功率譜估計中最常用的一種方法，這是因為AR模型參數(shù)的精確估計可以通過解一組線性方程求得，而對于MA和ARMA模型功率譜估計來說，其參數(shù)的精確估計需要解一組高階的非線性方程。在利用AR模型進行功率譜估計時，必須計算出AR模型的參數(shù)和激勵白噪聲序列的方差。這些參數(shù)的提取算法主要包括自相關(guān)法、Burg算法、協(xié)方差法、改進的協(xié)方差法，以及最大似然估計法。本章主要針對采用AR模型的兩種方法：Levinson-Durbi

3、n遞推算法、Burg遞推算法。實際中，數(shù)字信號的功率譜只能用所得的有限次記錄的有限長數(shù)據(jù)來予以估計，這就產(chǎn)生了功率譜估計這一研究領(lǐng)域。功率譜的估計大致可分為經(jīng)典功率譜估計和現(xiàn)代功率譜估計，針對經(jīng)典譜估計的分辨率低和方差性能不好等問題提出了現(xiàn)代譜估計，AR模型譜估計就是現(xiàn)代譜估計常用的方法之一。信號的頻譜分析是研究信號特性的重要手段之一，通常是求其功率譜來進行頻譜分析。功率譜反映了隨機信號各頻率成份功率能量的分布情況，可以揭示信號中隱含的周期性及靠得很近的譜峰等有用信息，在許多領(lǐng)域都發(fā)揮了重要作用。然而，實際應(yīng)用中的平穩(wěn)隨機信號通常是有限長的，只能根據(jù)有限長信號

4、估計原信號的真實功率譜，這就是功率譜估計。二．AR模型的構(gòu)建假定u(n)、x(n)都是實平穩(wěn)的隨機信號，u(n)為白噪聲，方差為，現(xiàn)在，我們希望建立AR模型的參數(shù)和x(n)的自相關(guān)函數(shù)的關(guān)系，也即AR模型的正則方程(normalequation)。由(1)由于u(n)是方差為的白噪聲，有（2）由Z變換的定義，，當(dāng)時，有h(0)=1。綜合（1）及（2）兩式，（3）在上面的推導(dǎo)中，應(yīng)用了自相關(guān)函數(shù)的偶對稱性。上式可寫成矩陣式：(4)（4）上述兩式即是AR模型的正則方程，又稱Yule-Walker方程。系數(shù)矩陣不但是對稱的，而且沿著和主對角線平行的任一條對角線上的元

5、素都相等，這樣的矩陣稱為Toeplitz矩陣。若x(n)是復(fù)過程，那么，系數(shù)矩陣是Hermitian對稱的Toeplitz矩陣。（4）式可簡單地表示為式中，為全零列向量，R是的自相關(guān)矩陣。可以看出，一個p階的AR模型共有p+1個參數(shù)，即，只要知道x(n)的前p+1個自相關(guān)函數(shù)，由(1)，(2)及(3)式的線性方程組即可求出這p+1個參數(shù)，即可求出x(n)的功率譜。三．AR模型階數(shù)的選擇AR模型的階次p一般事先是不知道的，需要事先選定一個稍大的值，在遞推的過程中確定。在使用Levinson遞推時，可以給出由低階到高階的每一組參數(shù)，且模型的最小預(yù)測誤差功率是遞減的

6、。直觀上講，當(dāng)達到所指定的希望值，或是不再發(fā)生變化時，其時的階次即是應(yīng)選的正確階次。因為是單調(diào)下降的，因此，的值降到多少才合適，往往不好選擇。為此，有幾個不同的準(zhǔn)則被提出，其中較常用的兩個是：最終預(yù)測誤差準(zhǔn)則：(1)(2)信息論準(zhǔn)則：式中N為數(shù)據(jù)的長度，當(dāng)階次k由1增加時，F(xiàn)PE(k)和AIC(k)都將在某一個k處取得極小值。將此時的k定為最合適的階次p。在實際運用時發(fā)現(xiàn)，當(dāng)數(shù)據(jù)較短時，它們給出的階次偏低，且二者給出的結(jié)果基本上是一致的。應(yīng)該指出，上面兩式僅為階次的選擇提供了一個依據(jù)，對所研究的某一個具體信號x(n)，究竟階次取多少為最好，還要在實踐中所得到的

7、結(jié)果作多次比較后，予以確定。四．Burg算法的理論分析Burg算法是較早提出的建立在數(shù)據(jù)基礎(chǔ)上的AR系數(shù)求解的有效算法[7]。其特點是：(1)令前后向預(yù)測誤差功率（5）為最小。(2)和的求和范圍從p至N-1，即，前后都不加窗，這時（6）在上式中，階次m由1至p時，（7）下式的遞推關(guān)系，即（8）（9）（10）式中。這樣，(5)式的僅是反射系數(shù)的函數(shù)。在階次m時，令相對為最小，即可估計出反射系數(shù)。將(6)、(7)及(8)式代入(5)式，令，可得使為最小的為式中。按此式估計出的滿足。按上式估計出后，在階次m時的AR模型系數(shù)仍然由Levinson算法遞推求出（11）（

8、12）式中。上面三式是假定在第(m-1