逆濾波器設(shè)計(jì)中若干基本問(wèn)題的探討*

曲萬(wàn)春，何國(guó)榮，梅鐵民

（沈陽(yáng)理工大學(xué)，遼寧 沈陽(yáng) 110159）

0 引言

在很多信號(hào)處理應(yīng)用中，如信道均衡、語(yǔ)音去混響[1-2]、能譜檢測(cè)[3]、軌道不平順檢測(cè)信號(hào)處理等[4]，都需要求解已知線性系統(tǒng)的逆系統(tǒng)，即逆濾波器設(shè)計(jì)。逆濾波器設(shè)計(jì)方法有很多，如最小均方誤差[5-7]、復(fù)倒譜法[8]、基于LMS或卡爾曼濾波的算法等[9]，其中最簡(jiǎn)單也最常用的方法是最小均方誤差設(shè)計(jì)方法。最小均方誤差法設(shè)計(jì)的逆系統(tǒng)與原系統(tǒng)的級(jí)聯(lián)在均方誤差最小化意義下逼近一個(gè)純時(shí)延全通系統(tǒng)。

在最小均方誤差逆濾波器設(shè)計(jì)中，一個(gè)經(jīng)常被忽視的問(wèn)題是原、逆濾波器級(jí)聯(lián)系統(tǒng)所逼近的全通系統(tǒng)的時(shí)延問(wèn)題（或群延遲問(wèn)題）[6]。級(jí)聯(lián)系統(tǒng)的時(shí)延大小取決于原系統(tǒng)的相位特性。最小均方誤差逆濾波器設(shè)計(jì)包括單輸入-單輸出系統(tǒng)（Single Input Single Output，SISO）、單輸入-多輸出（Single Input Multiple Output，SIMO）系統(tǒng)和多輸入-多輸出（Multiple Input Multiple Output，MIMO）系統(tǒng)的逆系統(tǒng)設(shè)計(jì)。本文中只討論前兩種情況。

1 SISO系統(tǒng)的逆系統(tǒng)設(shè)計(jì)

設(shè)單信道系統(tǒng)的輸入輸出關(guān)系為：

式中，s(n)和x(n)分別為系統(tǒng)的輸入、輸出信號(hào)；c(n)(n=0,1,2,…,Lc-1)為系統(tǒng)的單位取樣響應(yīng)。

式（1）的z域表示為：

在最小均方誤差意義下使g逼近理想系統(tǒng)ql：

式中，ql=[0,0,…,1(第l個(gè)分量),0,…,0]T是Lg×1維列矢量，代表一個(gè)純時(shí)延的理想全局響應(yīng)。

式（8）的最小均方解為：

式中，列矢Cl=CTql是卷積矩陣C的第l行的轉(zhuǎn)置。

最小均方誤差值為：

如上所述，SISO系統(tǒng)只能得到最小均方誤差意義下的最優(yōu)逆濾波器，不可能得到完全解卷積的逆濾波器，但SIMO系統(tǒng)可以得到完全解卷積逆系統(tǒng)。

2 SIMO逆系統(tǒng)設(shè)計(jì)

SIMO系統(tǒng)與SISO系統(tǒng)相比，具有本質(zhì)上的不同，有完全解卷積的有限沖激響應(yīng)逆系統(tǒng)。SIMO系統(tǒng)的輸入輸出關(guān)系為：

式中，i=1,2,…,N，xi(n)和ci(n)分別是第i路信道的輸出信號(hào)（觀測(cè)信號(hào)）和單位取樣響應(yīng)；s(n)是系統(tǒng)的輸入信號(hào)。在N路信道沖激響應(yīng)ci(k)已知且各個(gè)信道間沒(méi)有公共零點(diǎn)的情況下，多路解卷積系統(tǒng)由MINT定理給出[5-6]。

令N路解卷積濾波器為hi(n)(i=1,2,…,N;n=0,1,…,Lh)，則解卷積信號(hào)為：

式中，全局響應(yīng)g(k)可表示為：

則SIMO解卷積問(wèn)題等價(jià)于如下的優(yōu)化問(wèn)題：

當(dāng)Lc、N不同值時(shí)，由得到的Lh值不一定是整數(shù)。如果Lh取大于的整數(shù)，則LgNLh，式（23）為超定問(wèn)題，有唯一最小二乘解：

方法描述如下：設(shè)N路解卷積濾波器的長(zhǎng)度分別為L(zhǎng)h1,Lh2,…,LhN，不失一般性地假設(shè)Lh1最大，且則全局響應(yīng)的長(zhǎng)度Lg=Lh1+Lc-1。分配各個(gè)解卷積濾波器的長(zhǎng)度值，使得解卷積濾波器組構(gòu)成的矢量總長(zhǎng)度滿足條件：

考慮到Lg=Lh1+Lc-1，式（26）也等價(jià)于：

式（28）同時(shí)隱含了對(duì)每一個(gè)解卷積濾波器長(zhǎng)度的約束條件：

當(dāng)Lh1,Lh2,…,LhN給定后，即可構(gòu)造矩陣其中Ci矩陣仍如式（22）定義，但只取其前Lhi列，這樣得到的濾波器組可以實(shí)現(xiàn)完全解卷積。由于Lh1,Lh2,…,LhN的取值有不同的組合，再加上不同的取值，因此完全解卷積濾波器組并不唯一。此外，如果假設(shè)Lh2=Lc-1，則Lhi=0(i≠1,2)，即只有兩路觀測(cè)信號(hào)即可實(shí)現(xiàn)完全解卷積，也說(shuō)明多路信號(hào)解卷積中至少要有兩路觀測(cè)信號(hào)才能實(shí)現(xiàn)完全解卷積。

3 仿真實(shí)驗(yàn)

3.1 SISO系統(tǒng)仿真實(shí)驗(yàn)

在這個(gè)仿真實(shí)驗(yàn)中，原系統(tǒng)沖激響應(yīng)c0(n)是仿真得到的Lc=800的房間沖激響應(yīng)序列，如圖1所示。

圖1 原系統(tǒng)的沖激響應(yīng)序列c0(n)

為了模擬原系統(tǒng)沖激響應(yīng)c0(n)估計(jì)誤差對(duì)全局響應(yīng)g(n)的影響，在原系統(tǒng)沖激響應(yīng)上增加不同強(qiáng)度的零均值高斯白噪聲，即c(n)=c0(n)+v(n)，用估計(jì)誤差信噪比來(lái)表示c(n)估計(jì)精度。在最佳全局響應(yīng)時(shí)延下得到最小均方誤差意義下的解卷積濾波器h(n)后，計(jì)算全局響應(yīng)g(n)=c0(n)*h(n)相對(duì)理想全局響應(yīng)q(n)的均方誤差為：

圖2中給出了ERR(SIR,Lh)與SIR、Lh之間的關(guān)系曲線?？梢钥闯觯瑢?duì)于給定的原系統(tǒng)估計(jì)誤差信噪比SIR，解卷積濾波器越長(zhǎng)（圖2中的曲線由下向上對(duì)應(yīng)Lh由100增加到1 500），解卷積效果越好，但隨著解卷積濾波器長(zhǎng)度的增加，解卷積效果的改善越來(lái)越?。粚?duì)于給定長(zhǎng)度Lh解卷積濾波器，原系統(tǒng)估計(jì)誤差的SIR達(dá)到某個(gè)值（隨Lh增大而增大）以上時(shí)，對(duì)解卷積結(jié)果影響很小，這時(shí)卷積結(jié)果的誤差主要取決于解卷積系統(tǒng)h(n)的逼近誤差。因此，對(duì)于給定長(zhǎng)度的解卷積濾波器，一味增加原系統(tǒng)估計(jì)精度是沒(méi)有意義的。

3.2 SIMO系統(tǒng)仿真實(shí)驗(yàn)

4路實(shí)測(cè)房間沖激響應(yīng)用于SIMO系統(tǒng)解卷積仿真實(shí)驗(yàn)Lc=818，如圖3所示。

在滿足MINT的定理?xiàng)l件下，首先定義不同原系統(tǒng)估計(jì)誤差信噪比SIR、不同全局響應(yīng)時(shí)延l對(duì)應(yīng)的全局響應(yīng)均方誤差函數(shù)ERRMINT：

4個(gè)解卷積濾波器的長(zhǎng)度分別為L(zhǎng)h1=Lh2=273，Lh3=Lh4=272。原系統(tǒng)c0(n)估計(jì)誤差信噪比SIR對(duì)全局響應(yīng)均方誤差ERR的影響如圖4所示。圖4中每條曲線對(duì)應(yīng)不同的全局響應(yīng)時(shí)延l=0,20,40,…,1 080。從圖4可見(jiàn)，只有在高估計(jì)誤差信噪比下，ERR近似地隨SIR的增加而線性增加；當(dāng)SIR較低時(shí)ERR很小，甚至為負(fù)值，無(wú)法實(shí)現(xiàn)解卷積。

圖2 原系統(tǒng)c0(n)估計(jì)誤差SIR（橫軸）、解卷積濾波器長(zhǎng)度Lh（對(duì)應(yīng)圖中不同曲線）對(duì)全局響應(yīng)均方誤差（縱軸）的影響（曲線的排列順序與圖例的順序一致）

圖3 房間4通道沖激響應(yīng)序列

對(duì)ERRMINT(SIR,l)在SIR方向取平均，得到，如圖5所示。找到圖5中最大值對(duì)應(yīng)的l值，即為最優(yōu)全局響應(yīng)時(shí)延，在本例中l(wèi)=100。

對(duì)于SIMO系統(tǒng)，在不滿足MINT定理時(shí)得到的是在最佳全局響應(yīng)時(shí)延下的最小均方誤差解。解卷積濾波器等長(zhǎng)并滿足關(guān)系（向下取整，在本例中Lh≤272），仍用式（30）中的ERR來(lái)衡量算法性能，結(jié)果如圖6所示。與SISO系統(tǒng)的結(jié)果類似，即越逼近其上限，解卷積效果越好，但噪聲穩(wěn)定性越差。圖6中每條曲線對(duì)應(yīng)不同解卷積濾波器長(zhǎng)度Lh=27、54、81、108、135、162、189、216、243、270。

圖4 在滿足MINT定理?xiàng)l件下，全局響應(yīng)均方誤差ERR與原系統(tǒng)c0(n)估計(jì)誤差信噪比SIR關(guān)系曲線（每條曲線對(duì)應(yīng)不同的全局響應(yīng)時(shí)延0,20,40,…,1 080；粗黑線對(duì)應(yīng)最佳全局響應(yīng)時(shí)延下的ERR）

圖5 不同SIR下的平均ERR曲線（曲線最大值位置對(duì)應(yīng)最優(yōu)全局響應(yīng)時(shí)延l=100）

圖6 在不滿足MINT定理?xiàng)l件下，全局響應(yīng)均方誤差與原系統(tǒng)估計(jì)誤差信噪比SIR關(guān)系曲線

4 結(jié)語(yǔ)

理論上，SISO系統(tǒng)的逆濾波器長(zhǎng)度越長(zhǎng)越好，但由于有限沖激響應(yīng)逆系統(tǒng)只是理想逆系統(tǒng)的最小均方誤差逼近，因此在原系統(tǒng)存在估計(jì)誤差時(shí)，過(guò)長(zhǎng)的逆濾波器并不能提高解卷積性能，反而增加計(jì)算負(fù)擔(dān)。對(duì)于SIMO系統(tǒng)，正如文獻(xiàn)所述，滿足MINT定理的逆系統(tǒng)對(duì)噪聲過(guò)于敏感，而不滿足MINT定理的最小均方誤差逆系統(tǒng)相對(duì)來(lái)說(shuō)噪聲穩(wěn)定性更好，且該最小均方誤差逆系統(tǒng)的沖激響應(yīng)的長(zhǎng)度是有上限的。它的上限由MINT定理給出，這與SISO系統(tǒng)截然不同。此外，無(wú)論是SISO系統(tǒng)還是SIMO系統(tǒng)，都可以通過(guò)對(duì)全局響應(yīng)時(shí)延優(yōu)化來(lái)改善逆系統(tǒng)性能。