《維納濾波》word版

《《維納濾波》word版》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫。

1、第3講：Wiener濾波Wiener濾波器是從統(tǒng)計意義上的最優(yōu)濾波,它要求輸入信號是寬平穩(wěn)隨機序列,本章主要集中在FIR結(jié)構(gòu)的Wiener濾波器的討論。由信號當(dāng)前值與它的各階延遲，估計一個期望信號，輸入信號是寬平穩(wěn)的，和是聯(lián)合寬平穩(wěn)的,要求這個估計的均方誤差最小.。Wiener濾波器的幾個實際應(yīng)用實例如下：①通信的信道均衡器。圖1.信道均衡器的結(jié)構(gòu)示意②系統(tǒng)辨識：圖2.線性系統(tǒng)辨識的結(jié)構(gòu)③一般結(jié)構(gòu):12圖3.Wiener濾波器的一般結(jié)構(gòu)Wiener濾波器的目的是求最優(yōu)濾波器系數(shù)，使最小?！?.1從

2、估計理論觀點導(dǎo)出Wiener濾波FIR結(jié)構(gòu)(也稱為橫向)的Wiener濾波器的核心結(jié)構(gòu)如圖4所示.圖4.橫向Wiener濾波器為了與第2講中估計理論一致，假設(shè)信號，濾波器權(quán)值均為實數(shù)由輸入和它的1至（M-1）階延遲，估計期望信號，確定權(quán)系數(shù)使估計誤差均方值最小，均方誤差定義為：這里估計寫為：除了現(xiàn)在是波形估計外，與線性Bayesian估計一一對應(yīng)。12R（零均值假設(shè)）這里,Wiener濾波與線性Bayesian估計變量之間具有一一對應(yīng)關(guān)系,設(shè)最優(yōu)濾波器系數(shù)為，由線性Bayesian估計得到Wien

3、er濾波器系數(shù)對應(yīng)式：上式后一個方程稱為Wiener-Hopf方程,或結(jié)論：1)Wiener濾波器是線性FIR濾波器中的最優(yōu)濾波器，但非線性濾波可能會達到更好結(jié)果。2)在聯(lián)合高斯條件下，Wiener濾波也是總體最優(yōu)的（①從Bayesian估計意義上講是這樣，②要滿足平穩(wěn)條件）3)從線性貝葉斯估計推導(dǎo)過程知，在濾波器系數(shù)取非最優(yōu)的w時，其誤差性能表示：它是w的二次曲面，只有一個最小點，時，§3.2維納濾波：從正交原理和線性濾波觀點分析Wiener濾波器Wiener濾波器是一個最優(yōu)線性濾波器，圖3是一

4、個一般表示框圖，濾波器核是IIR或FIR的，為了導(dǎo)出后續(xù)常用的一些工具,我們導(dǎo)出最優(yōu)濾波器的正交原理,并從正交原理出發(fā)重新導(dǎo)出一般的Wiener濾波器方程,先推導(dǎo)適應(yīng)于IIR和FIR的一般結(jié)論，然后重點討論FIR。討論一般的復(fù)數(shù)形式?！ぽ斎脒^程?！V波器系數(shù)，（權(quán)系數(shù)）·希望的響應(yīng)d[n]·輸出誤差：12·正交性原理對復(fù)數(shù)據(jù)情況，推導(dǎo)一般結(jié)論，實數(shù)據(jù)是特例。均方誤差是：設(shè)權(quán)系數(shù)定義遞度算子其中符號是遞度算子作用于J，其中第k項為：要求的值,使得J最小,即或等價：將代入遞度算子得到：由得到：12代入

5、表達式整理得：當(dāng)時，J達到最小。設(shè)J達最小時，用表示權(quán)系數(shù)和誤差e[n]，且則有：，和，以上兩式為正交性原理，達到最優(yōu)濾波時，誤差和輸入正交。推論：·最小均方誤差：在達最優(yōu)時，也寫成：，表示由張成的空間對d[n]的估計（最優(yōu)線性估計）。也可以寫成：由和正交性得：12即：·維納－霍夫方程（Wiener-Hopf）由正交性原理即：得：由于：得：這就是Wiener-Hopf方程，解此方程，可得到最優(yōu)權(quán)。對于M階FIR濾波器，（橫向濾波器）wiener-Hopf方程變?yōu)椋?，·矩陣形式：令?12Winer

6、-Hopf方程寫為：這里解方程求得：·再看最小均方誤差由得則分析：由正交原理導(dǎo)出Wiener濾波器和從線性貝葉斯估計得到wiener濾波器,分別得到一些有益的啟示,從估計理論看,Wiener濾波器只是一個線性Beyesian估計,它是最優(yōu)估計的一個線性逼近,只有在高斯情況下,它才是真正的最優(yōu)濾波器,在其它分布情況下,非線性濾波器可以達到比線性最優(yōu)濾波器更優(yōu)的結(jié)果.從一般線性最優(yōu)濾波器的正交原理出發(fā),我們?nèi)菀缀鲆曔@些限制?！ふ`差性能表面12由直接代入整理得：由上式，可以看出，J是Wk的二次曲面，是碗

7、狀曲面，碗口向上，Jmin在碗底，其實，由上式直接對wk求導(dǎo)，得到一組方程，正是wiener-Hopf方程。上式也可以直接寫成矩陣形式它可以整理成如下形式：上式在時，達最小，性能表面也可以寫成：由于故令通過坐標(biāo)變換，得到如上規(guī)范形式，對于一個給定,有：12這是超相圓，為其一個軸。數(shù)值例子：①希望響應(yīng)是一個AR(1)過程，，是白噪聲，,由白噪聲驅(qū)動的產(chǎn)生該過程的傳輸函數(shù)為:②經(jīng)過了一個通信信通，信道的傳輸函數(shù)為，并加入了白噪聲即：通道模型如圖5所示:圖5.通道模型③求解：一個二階FIR結(jié)構(gòu)Wiene

8、r濾波器，目的是由u[n]恢復(fù)d[n]解：①是一個過程，②在中,是一個二階過程,相當(dāng)于12由二階參數(shù)，確定,由Yule-walker方程：反解.得由上確定x[n]的自相關(guān)矩陣為:但：③求由：,和代入上式12得：故④性能表面⑤最優(yōu)系數(shù)⑥最小均方誤差：①規(guī)范誤差性能表面解12這是一個隨圓，主軸，副軸性能表面等高線如圖6.圖6.性能表面的等高線12