多種功率譜估計(jì)的比較.doc

《多種功率譜估計(jì)的比較.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、.多種功率譜估計(jì)的比較1．實(shí)驗(yàn)?zāi)康模篴.了解功率譜估計(jì)在信號(hào)分析中的作用；b.掌握隨機(jī)信號(hào)分析的基礎(chǔ)理論，掌握參數(shù)模型描述形式下的隨機(jī)信號(hào)的功率譜的計(jì)算方法；c.掌握在計(jì)算機(jī)上產(chǎn)生隨機(jī)信號(hào)的方法；d.了解不同的功率譜估計(jì)方法的優(yōu)缺點(diǎn)。2．實(shí)驗(yàn)準(zhǔn)備：有三個(gè)信號(hào)源，分別代表三種隨機(jī)信號(hào)（序列）。信號(hào)源1：其中，z(n)是一個(gè)一階AR過程，滿足方程：e(n)是一高斯分布的實(shí)白噪聲序列，方差信號(hào)源2和信號(hào)源3：都是4階的AR過程，它們分別是一個(gè)寬帶和一個(gè)窄帶過程，滿足方程：e(n)是一高斯分布的實(shí)白噪聲序列，方差，參數(shù)如下：參數(shù)信號(hào)源a(1)a(2)a(3)a(4)

2、信號(hào)源2－1.3001.200－0.6000.2501.000信號(hào)源3－2.7503.799－2.6500.9281.0003．實(shí)驗(yàn)內(nèi)容：a.描繪出這三個(gè)實(shí)驗(yàn)信號(hào)的真實(shí)功率譜波形。b.在計(jì)算機(jī)上分別產(chǎn)生這個(gè)三個(gè)信號(hào)，令所得到的數(shù)據(jù)長(zhǎng)度256=N。注意：產(chǎn)生信號(hào)的時(shí)候注意避開起始瞬態(tài)點(diǎn)。例如，可以產(chǎn)生長(zhǎng)度為512的信號(hào)序列，然后取后面256個(gè)點(diǎn)作為實(shí)驗(yàn)數(shù)據(jù)。c.分別用如下的譜估計(jì)方法，對(duì)三個(gè)信號(hào)序列進(jìn)行譜估計(jì)。1、經(jīng)典譜估計(jì)z周期圖法z自相關(guān)法z平均周期圖法（Bartlett法）zWelch法（可選每段64點(diǎn)，重疊32點(diǎn)，用Hamming窗）2、現(xiàn)代譜估計(jì)zY

3、ule-Walker方程（自相關(guān)法）z最小二乘法注：階次p可在3－20之間，由自己給定。4.實(shí)驗(yàn)結(jié)果及分析頁腳.1分析信號(hào)源11>周期圖法周期圖法又稱直接法，是直接建立在功率譜的定義式上的。圖1一次周期圖法圖250次周期圖法平均（紅色：平均，藍(lán)色：?jiǎn)未危┲芷趫D法在理論上漸近無偏的非一致估計(jì)，數(shù)據(jù)量越大，周期圖沿頻率軸的起伏變化愈激烈所求得的功率譜振蕩劇烈，信號(hào)方差較大，不利于對(duì)功率信號(hào)的分析。但是從上圖上可以看出，對(duì)于確定信號(hào)成分較大信號(hào)源里，周期圖法能夠較好的表現(xiàn)出信號(hào)的頻譜成分。本信號(hào)為確定信號(hào)和白噪聲相加，且性噪比較大，因此信號(hào)有較好的平均性。2>自相

4、關(guān)法原理：由維納-辛欽公式，經(jīng)自相關(guān)函數(shù)間接獲得頁腳.圖3一次自相關(guān)法圖450次自相關(guān)法平均（紅色：平均，藍(lán)色：?jiǎn)未危┳韵嚓P(guān)法是由維納-辛欽公式出發(fā)的，本質(zhì)上是對(duì)周期圖法的插值，因此而這本質(zhì)上來說是一致的。由于采用差值的方法，會(huì)使譜線相對(duì)平滑一些。同樣，對(duì)于本信號(hào)源（確定信號(hào)+白噪聲），自相關(guān)法較好的跟隨了譜線峰值，且也有較好的平均性。3>平均周期圖法（Bartlett法）原理：K個(gè)獨(dú)立同分布的隨機(jī)變量的均值之方差，等于單個(gè)變量方差的1/K。方法：長(zhǎng)數(shù)據(jù)N分成K段，每段M=N/K，每段用周期圖法求譜：頁腳.K段平均后的譜估計(jì)為：圖5一次BT法圖650次BT法

5、平均（紅色：平均，藍(lán)色：?jiǎn)未危┢骄芷趫D法是通過改進(jìn)的經(jīng)典譜估計(jì)法，是以犧牲分辨率和偏差來改善方差性能的。但是在本信號(hào)源（確定信號(hào)+白噪聲）分析中，很好的分辨出了相鄰的兩根譜線值。4>Welch法原理：1把N個(gè)數(shù)據(jù)分成K段，每段可以互相獨(dú)立(如平均周期圖法),也可以互相交疊,例如交疊一半,即K=N/L或K=(N-L/2)/(L/2)2.再把每段數(shù)據(jù)乘上窗函數(shù)w(n)(如加窗平滑法)后作DFT。頁腳.圖7一次welch法圖850次welch法平均（紅色：平均，藍(lán)色：?jiǎn)未危￤elch法結(jié)合了平均周期法和加窗平滑法的優(yōu)點(diǎn)，保證了方差性能，但是由于每段數(shù)據(jù)點(diǎn)數(shù)較少，

6、所以在跟隨峰值的方面較差。對(duì)于本信號(hào)，使用welch法有較好的平均性。5>Yule-Walker方程（自相關(guān)法）目標(biāo)：找到已知參數(shù)和未知參數(shù)的關(guān)系，以便求解未知參數(shù)。未知參數(shù)：頁腳.已知參數(shù)：圖9一次Y-W法圖1050次Y-W法平均（紅色：平均，藍(lán)色：?jiǎn)未危℡ule-Walker方程（自相關(guān)法）是現(xiàn)代譜估計(jì)方法。從圖中可以看出，相對(duì)于welch頁腳.法有較好的分辨率。對(duì)于本信號(hào)（確定信號(hào)+白噪聲），使用該方法的平均性較差。6>最小二乘法是零均值方差為的平穩(wěn)隨機(jī)過程，把上式寫成則參數(shù)Φ的線性最小二乘估計(jì)：圖11一次LS法頁腳.圖1250次LS法平均（紅色：平均

7、，藍(lán)色：?jiǎn)未危┳钚《朔ㄊ蔷烙?jì)，是在Y-W方法上改進(jìn)，采用的數(shù)據(jù)量大，故有較好的分辨率。從圖上可以看出，分辨率和方差性能都較好。2分析信號(hào)源21>周期圖法圖1350次周期圖法（紅色：平均，藍(lán)色：?jiǎn)未危┲芷趫D法的方差性能比較差，因此譜線的起伏比較大。2>自相關(guān)法圖1450次自相關(guān)法平均（紅色：平均，藍(lán)色：?jiǎn)未危┳韵嚓P(guān)法在本質(zhì)上和周期圖法是相同的，方差性能比較差，因此譜線的起伏比較大。3>平均周期圖法頁腳.圖1550次周期法平均平均周期圖法從理論上改善了方差性能，是以犧牲分辨率來提高方差性能的，該方法的估計(jì)與自相關(guān)法和周期圖法相比較，曲線的方差特性明顯變好，但

8、是分辨率較低。4>Welch法圖1550次welch