碩士生現(xiàn)代信號處理-經(jīng)典譜分析 已讀

碩士生現(xiàn)代信號處理-經(jīng)典譜分析已讀第六章經(jīng)典譜分析第六章經(jīng)典譜分析6.2周期圖法6.3平滑的周期圖法6.1預備知識6.4小結(jié)6.1預備知識1.確定信號的能量譜密度設x(n)為能量有限確定離散信號，即 那么對x(n)可做付氏變換：

由Parseval定理，有：6.1預備知識

表示信號分布在頻率區(qū)間 上的能量，稱為x(n)的能量譜密度(EnergyDensitySpectrum)。2.隨機信號的功率譜i)隨機過程的功率譜密度(PowerDensitySpectrum)

定義6.1預備知識那么由Parseval定理，有：6.1預備知識信號的平均功率：定義功率譜密度：對廣義平穩(wěn)隨機過程，有6.1預備知識可以證明，對于零均值平穩(wěn)隨機過程，其自相關(guān)函數(shù) 與功率譜密度為一對付氏變換，即6.1預備知識ii)離散時間隨機過程的功率譜密度 設x(n)為零均值廣義平穩(wěn)離散時間過程，其自相關(guān)函數(shù)為： 那么x(n)的功率譜密度為的SFT,即x(n)的平均功率為6.1預備知識3.譜估計(分析)的定義利用給定的有限個樣本數(shù)據(jù)來估計一個平穩(wěn)隨機過程的功率譜密度(PowerDensitySpectrum)的過程稱為譜分析或譜估計。6.1預備知識4.隨機信號自相關(guān)函數(shù)的計算方法1i)估計式： 設觀察到N個零均值樣本序列x(n)的值：

假設當n<p及n>p+N-1時，x(n)=0。令6.1預備知識那么：6.1預備知識ii)估計性能當時，上式為一致估計，即當N一定，時，較大。自相關(guān)函數(shù)估計的一致性證明

估計偏差:自相關(guān)函數(shù)估計的一致性證明

估計方差:6.1預備知識5.隨機信號自相關(guān)函數(shù)的計算方法2i)估計式：6.1預備知識ii)估計性能當時，,是一致估計。6.2周期圖法1.定義設x(n)為零均值平穩(wěn)隨機序列。定義N點序列 及其付氏變換 采用下式計算x(n)的功率譜密度6.2周期圖法自相關(guān)的估計方法2：6.2周期圖法用得到的功率譜估計稱為周期圖(Periodogram)，記作，即6.2周期圖法周期圖求出的功率譜與基于自相關(guān)估計方法二計算的功率譜等價。假設觀察到x(n)的N個值x(p),x(p+1),…x(p+N+1)，可通過SFT直接求得，然后求得，而不必先估計自相關(guān)函數(shù)，計算方便。6.2周期圖法2.周期圖譜估計的性能i)周期圖的均值 設A(n)為長度為N的矩形序列，即：那么6.2周期圖法

其中是一個三角函數(shù)，又稱Bartlett窗函數(shù)：于是有：周期圖的均值是真實頻譜與三角窗頻譜的卷積。6.2周期圖法6.2周期圖法周期圖是功率譜的有偏估計，即由于

即周期圖在時是的無偏估計。6.2周期圖法ii)周期圖的方差 設序列是一個零均值、方差為的白色高斯過程的樣本，那么6.2周期圖法

，即周期圖非一致估計。對于非白高斯樣本，有6.2周期圖法iii)周期圖的協(xié)方差假設，那么6.2周期圖法以的整數(shù)倍為頻率間隔的周期圖的值是不相關(guān)的。當N增大時，協(xié)方差為零的功率譜樣本之間的間距減小，因此周期圖的起伏增快。6.2周期圖法周期圖的起伏現(xiàn)象實例N=16N=32N=646.2周期圖法3.周期圖譜估計的特點頻譜泄漏現(xiàn)象：當真實頻譜為譜線時，周期圖估計出的頻譜不再是譜線。原因：是與窗函數(shù)頻譜的卷積。本質(zhì)原因是觀察數(shù)據(jù)為有限長。觀察數(shù)據(jù)越短，頻譜泄漏越嚴重。6.2周期圖法頻率分辨率頻率分辨率：區(qū)分兩個相鄰頻率分量的能力。頻率分辨率與觀察長度成反比。周期圖的頻率分辨率為1/N。三角窗函數(shù)的頻譜：6.2周期圖法例如：信號的真實頻譜周期圖估計頻譜6.3平滑的周期圖法1.巴特利特平均周期圖法i)引出：設序列是不相關(guān)的隨機變量，每一個具有期望值，方差。 定義，那么

屢次測量平均可以減小測量方差。6.3平滑的周期圖法ii)方法描述：數(shù)據(jù)分段；計算各段的周期圖；平均各段的周期圖。將序列分成L段，每段M(N=LM)個樣本，那么第i段樣本序列為：的周期圖為6.3平滑的周期圖法 假設很小，那么可假定各段周期圖是不相關(guān)的。將各段周期圖的平均作為對功率譜的估計，即：6.3平滑的周期圖法iii)估計性能偏差：

Bartlett估計是有偏估計，但。6.3平滑的周期圖法估計方差

iv)結(jié)論：Bartlett估計從減少方差的角度對周期圖法進行了改進，得到了頻譜的一致估計。6.3平滑的周期圖法2.窗化法i)思路: 用適宜的窗函數(shù)譜對周期圖平滑，即：ii)窗函數(shù)的選擇 6.3平滑的周期圖法iii)估計性能

其中為Bartlett窗。6.3平滑的周期圖法 假設，于是：iv)結(jié)論：窗化法使頻譜估計結(jié)果起伏變小，但分辨率變差。6.3平滑的周期圖法3.Welch法思想:綜合分段平均法和窗化法修正1：選擇窗函數(shù)直接對信號窗化U為歸一化因子：修正2：分段時各段之間有重疊。6.3平滑的周期圖法4.譜估計例如信號真實頻譜6.3平滑的周期圖法周期圖估計的頻譜(N=64)6.3平滑的周期圖法周期圖平均估計的頻譜(M=64，L=50)6.3平滑的周期圖法周期圖估計的頻譜(N=128)6.3平滑的周期圖法周期圖平均估計的頻譜(M=128，L=50)6.3平滑的周期圖法(續(xù)15)平滑的周期圖法總結(jié)估計的期望值:真值與窗函數(shù)頻譜卷積。分辨率變差：窗函數(shù)頻譜的主瓣寬度增