楊瑞麗，郭 瑩

（沈陽工業(yè)大學(xué)信息科學(xué)與工程學(xué)院,沈陽110870）

1 引言

聲學(xué)回聲一直是影響通信質(zhì)量的關(guān)鍵問題[1-2]。目前，典型的回聲消除器設(shè)計都在假設(shè)回聲路徑是線性的，可以采用有限長脈沖響應(yīng)（FIR）濾波器進行建模。然而，近幾年，隨著低成本小型設(shè)備在免提手機和視頻會議等應(yīng)用中的廣泛使用，放大器和揚聲器引入非線性回聲的問題日益嚴重[3-4]。這種非線性回聲難以由傳統(tǒng)的線性回聲消除器消除，回聲路徑中存在的非線性也會限制線性濾波器的性能，干擾聲學(xué)脈沖響應(yīng)的正確估計，降低語音信號的感知質(zhì)量[5-6]。因此，如何消除非線性回聲日益成為信號處理領(lǐng)域的研究熱點[7-10]。

目前，雖然已有部分學(xué)者研究了很多非線性回聲消除法，但其都缺少相應(yīng)的收斂性分析。為此，新提出的一種非線性回聲消除方法，在此基礎(chǔ)上，進一步對非線性濾波器系數(shù)和線性濾波器系數(shù)進行詳細的收斂性分析，并進行了大量的MATLAB實驗仿真，驗證對非線性回聲消除的收斂性分析的合理和可行性。

2 非線性回波消除

本研究采用基于Hammerstein模型的無記憶非線性回聲消除方案，以非線性濾波器建模非線性器件，以線性濾波器建模線性回聲路徑。系統(tǒng)整體構(gòu)成原理如圖1所示。

圖1 非線性回聲消除系統(tǒng)

對于無記憶非線性濾波器，采取一種級數(shù)展開的截斷模型來描述，如下式所示：

其中，

這里，

由式(2)、(3)、(4)可以得到：

此處所用的兩個濾波器的自適應(yīng)算法均采用NLMS算法，非線性濾波器系數(shù)的級聯(lián)NLMS迭代方程如下式：

相似地，線性濾波器系數(shù)的迭代方程可以表示為下式：

3 完美假設(shè)下的收斂性分析

3.1 線性系數(shù)誤差的收斂性分析

完美假設(shè)指的是濾波器系數(shù)對未知系統(tǒng)的估計是完全匹配的，這里是假設(shè)非線性系數(shù)估計完美時即，來對線性系數(shù)誤差進行收斂性分析。

線性系數(shù)誤差和非線性系數(shù)誤差分別表示如下：

公式(7)又可以表示為：

這里假設(shè)系統(tǒng)中的近端噪聲v(n)是零均值且獨立的，即這兩項不存在。

對其取跡可得：

通過遞歸，公式(14)可以重新表示如下：

3.2 非線性系數(shù)誤差的收斂性分析

4 不完美假設(shè)下的收斂性分析

4.1 線性系數(shù)誤差的收斂性分析

在此不再假設(shè)非線性系數(shù)是完美估計，而是考慮到了非線性系數(shù)的估計誤差所產(chǎn)生的影響。公式(9)將會成為線性參數(shù)的更新方程，有以下關(guān)系，其中

將線性系數(shù)估計誤差的相關(guān)矩陣劃分為兩個部分，可以將其重寫為：

由于假設(shè)非線性系數(shù)的估計誤差和線性系數(shù)的估計誤差的期望值均為零，以及當(dāng)?shù)卦肼暤木狄矠榱?，在此條件下，公式(22)中只保留了第五項，即如下式表示：

將公式(18)和線性濾波器的輸出代入公式(22)中，其可以被重寫為：

綜合以上分析，公式(19)可以重寫表示為：

接下來，對公式(24)進行求跡，

當(dāng)Q=3時，

4.2 非線性系數(shù)誤差的收斂性分析

在不完美假設(shè)下對非線性系數(shù)誤差進行分析，其過程與在4.1中類似，

接下來，對此式求跡，具體表示如下：

當(dāng)Q=3時，

5 仿真實驗

5.1 仿真條件

在仿真中，采用零均值的高斯白噪聲作為系統(tǒng)的輸入信號。無記憶非線性濾波器的階數(shù)取3,FIR濾波器的長度和線性聲學(xué)回聲路徑的長度一致，均為128。背景噪聲采用20dB的高斯噪聲。無記憶非線性系數(shù)向量定義為a=[1.0,0.1,0.33]T，線性聲學(xué)回聲路徑可以由下式產(chǎn)生[11]：

其中，n,i=0,1,...,L-1,L為脈沖響應(yīng)長度，r(n)為[-0.2,0.2]的隨機數(shù)，δ為Dirac函數(shù)，L取128。如圖2所示。

圖2 聲學(xué)回聲路徑脈沖響應(yīng)h

在初始化時，a與h不同，h初始化為零向量，而a初始化為[1,0,0]。

5.2 性能指標(biāo)

非線性系數(shù)的權(quán)誤差向量范數(shù)WEVN表達式為：

線性系數(shù)的WEVN表達式為：

WEVN的值越小則說明自適應(yīng)濾波器越逼近于所跟蹤的未知系統(tǒng)。

5.3 仿真結(jié)果分析

仿真生成曲線與理論曲線的對比情況如圖3、圖4、圖5和圖6所示，分別對應(yīng)為公式(15)、公式(16)、公式(25)和公式(28)。

從圖3～圖6可以得出，無論是在完美假設(shè)條件下，還是在不完美假設(shè)條件下，線性系數(shù)和非線性系數(shù)誤差的仿真結(jié)果與理論結(jié)果之間的誤差保持在3dB左右，可以說明文中對線性濾波器系數(shù)和非線性濾波器系數(shù)的收斂性分析是合理的，這對非線性回波消除的研究具有很重要的意義。而且，在實際生活中有很多場合都會存在非線性回聲，比如，免提通話系統(tǒng)，視頻會議系統(tǒng)等。對于有關(guān)提高電信設(shè)備通話質(zhì)量的研究，此仿真結(jié)果可提供有力的參考。

圖3 完美假設(shè)下線性系數(shù)誤差收斂曲線

圖4 完美假設(shè)下非線性系數(shù)誤差收斂曲線

圖5 不完美假設(shè)下線性系數(shù)誤差收斂曲線

圖6 不完美假設(shè)下非線性系數(shù)誤差收斂曲線

6 結(jié)束語

本研究在提出一種新的非線性回聲消除方法的基礎(chǔ)上，在完美假設(shè)和不完美假設(shè)的兩個方面分別對線性系數(shù)誤差和非線性系數(shù)誤差進行了詳細的收斂性分析，并進行了相應(yīng)的MATLAB實驗仿真，與理論值進行對比，驗證了對線性系數(shù)誤差和非線性系數(shù)誤差的收斂性分析的合理性，對于非線性回聲消除的進一步研究具有非常重要的意義以及實用價值。