資源描述:
《維納濾波器課件.ppt》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、現(xiàn)代數(shù)字信號處理第一章維納濾波主要內(nèi)容1.1維納濾波問題描述1.2維納FIR濾波器1.3維納非因果IIR濾波器1.4維納因果IIR濾波器1.5深入了解的預(yù)備知識1.6維納IIR濾波器計(jì)算例子1.7維納預(yù)測器真實(shí)信號觀察/測量數(shù)據(jù)加性噪聲/干擾線性估計(jì)問題最小均方誤差(MMSE)估計(jì)(minimummean-squareerror)估計(jì)誤差1.1維納濾波問題描述維納濾波->對真實(shí)信號的最小均方誤差估計(jì)問題.平滑濾波預(yù)測這里我們主要考慮濾波問題,即……線性估計(jì)根據(jù)其取值范圍不同通常有下面幾種情況:問題在于估計(jì)濾波器的參數(shù)/單位沖激響應(yīng)序列正交方程:標(biāo)準(zhǔn)方程(Wi
2、ener-Hopfequations):維納-霍夫(Wiener-Hopf)/標(biāo)準(zhǔn)方程任何時(shí)刻的估計(jì)誤差都與用于估計(jì)的所有數(shù)據(jù)(即濾波器的輸入)正交注意:Adata-dependantlinearleastsquareerrorestimationWiener-Hopfequation-solutionsOrthogonalequation-decorrelationWienerFilters下標(biāo)i的取值范圍決定了FIR,非因果IIR,因果IIRFIR(FiniteImpulseResponse)WienerFilter1.2求解Wiener-HopfEqu
3、ations--FIR濾波器維納-霍夫方程展開為矩陣形式Solution:練習(xí)完成習(xí)題2.9(1)和2.9(2)后如何利用上面給出的公式計(jì)算FIR維納濾波器的參數(shù)?data?P=E[s(n)x(n)],R=E[x(n)xT(n)]hop=R-1PFIR維納濾波器結(jié)構(gòu)WienerFilters標(biāo)準(zhǔn)方程:1.3求解Wiener-HopfEquations--非因果IIR濾波器Solution:1.4求解Wiener-HopfEquations--因果IIR濾波器標(biāo)準(zhǔn)方程:WhiteningfilterOptimumcausalfilterforwhitein
4、put將IIR濾波器分解為兩部分第一部分為白化濾波器(將輸入信號變?yōu)榘自肼暎┑诙糠譃橐园自肼暈榧?lì)的最優(yōu)因果濾波器。計(jì)算步驟因式分解(譜分解定理)(2)分解為因果部分和非因果部分(3)計(jì)算系統(tǒng)函數(shù)(4)計(jì)算沖激響應(yīng)(逆Z變換)(5)計(jì)算最小均方誤差維納濾波器的均方誤差證明上式為什么維納濾波器比一般線性濾波器性能更好?MinimumPhaseSequence如果一個(gè)穩(wěn)定的因果序列具有有理Z變換并且其所有的零點(diǎn)和極點(diǎn)都位于單位圓內(nèi),則為最小相位序列。當(dāng)下式成立時(shí)為最小相位序列1.5一些預(yù)備知識1.最小相位序列例如:有限序列最小相位多項(xiàng)式,最小相位系統(tǒng)為什么一
5、個(gè)所有零點(diǎn)位于單位圓內(nèi)的序列具有最小相位滯后?假設(shè)A(Z)是一個(gè)M階多項(xiàng)式,僅有一個(gè)零點(diǎn)位于單位圓內(nèi):這里F(z)是一個(gè)M-1階多項(xiàng)式,其所有的零點(diǎn)都位于單位圓外Thezero:共軛倒序關(guān)系得到另一個(gè)序列B(z)Thezero:共軛系數(shù)倒序(1)幅度特性(單位圓上,即
6、z
7、=1)(2)相位滯后特性Phaselag:相位差:假設(shè)A(z)為M階多項(xiàng)式,其所有的零點(diǎn)都位于單位圓內(nèi)。最大相位序列最大相位多項(xiàng)式最大相位系統(tǒng)這里存在著個(gè)不同的序列(通過對各個(gè)因子作共軛倒序),這些序列具有相同的幅度特性但是不同的相位特性最小相位序列最小相位多項(xiàng)式最小相位系統(tǒng)任何實(shí)平穩(wěn)隨機(jī)
8、信號的有理功率譜都可以被唯一分解為下面形式:廣義平穩(wěn)隨機(jī)信號模型:(任何平穩(wěn)隨機(jī)信號可以看作是由白噪聲序列激勵(lì)一個(gè)因果和穩(wěn)定的線性時(shí)不變系統(tǒng)產(chǎn)生的輸出)信號合成信號白化一些預(yù)備知識:2.譜因子分解定理和廣義平穩(wěn)隨機(jī)信號模型任何平穩(wěn)隨機(jī)過程可以看作是一個(gè)穩(wěn)定的因果線性系統(tǒng)在白噪聲序列激勵(lì)下的輸出。(2)如果x(n)使用逆濾波器1/B(z)進(jìn)行濾波,則輸出為白噪聲序列。(3)由于信號x(n)和白噪聲由一個(gè)可逆的變換相互關(guān)聯(lián),可以相互求得,所以他們包含了同樣的信息。輸入為白噪聲,方差為因果IIR維納濾波器沖激響應(yīng)為:Stable,Causal.Thepolesi
9、nsidetheunitcircle.3.因果維納IIR(白噪聲作輸入)功率譜分解定理iscausalandoptimal4.因果維納IIR濾波器-----具有有理功率譜的輸入信號Given:信號模型測量模型其中,假設(shè)下面條件成立:白噪聲1.6維納因果IIR濾波器計(jì)算實(shí)例Solution:(1)輸入信號功率譜因式分解分子:RicattiEquation(2)分解為因果和非因果部分P為RicattiEquation的正解。(3)計(jì)算因果IIR濾波器系統(tǒng)函數(shù)f,G,P間的關(guān)系:methodforcomputingcausalIIRWienerfilter計(jì)算系統(tǒng)
10、函數(shù)步驟因式分解(譜分解定理)(2)分解為因果部分和