[發明專利]基于系數差的成比例仿射投影回聲消除方法

有效

• 申請號: 201410465261.0
• 申請日: 2014-09-12
• 公開（公告）號: CN104270539A
• 公開（公告）日: 2015-01-07
• 發明（設計）人: 趙海全;鄭宗生
• 申請（專利權）人: 西南交通大學
• 主分類號: H04M9/08
• 分類號: H04M9/08
• 代理公司: 成都博通專利事務所 51208
• 代理人: 陳樹明

1.一種基于系數差的成比例仿射投影回聲消除方法，其步驟如下：

A、遠端信號采樣

通過對遠端信號進行采樣獲得當前時刻n的輸入信號x(n)，將最近的L個輸入信號x(n),x(n-1),...,x(n-L+1)組成當前時刻n的輸入信號向量X(n)，X(n)＝[x(n),x(n-1),...,x(n-L+1)]^T；上標T表示向量轉置，L為自適應濾波器抽頭長度，L＝64、128、256、512；

B、回聲信號估計

將當前時刻n的輸入信號向量X(n)通過自適應濾波器得到回聲信號的估計值即

y^(n)=WT(n)X(n)]]>

其中，W(n)＝[w₁(n),w₂(n),K,w_i(n),Kw_L(n)]^T為當前時刻的自適應濾波器抽頭權系數向量，其初始值為零，w_i(n)為當前時刻的自適應濾波器的第i個抽頭權系數；

C、回聲信號消除

對近端信號采樣得到帶回聲的當前時刻的近端信號d(n)，將其減去步驟B獲得的回聲信號的估計值得到消除回聲的當前時刻的有用信號e(n)，e(n)=d(n)-y^(n);]]>

D、濾波器抽頭權系數更新

將最近P個時刻的輸入信號向量X(n),X(n-1),...,X(n-P+1)組成輸入信號矩陣U(n)，U(n)＝[X(n),X(n-1),...,X(n-P+1)]；其中，P為仿射投影階數，P＝2～8；

將最近P個時刻的有用信號e(n),e(n-1),K,e(n-P+1)構成當前時刻的有用信號向量E(n)，E(n)＝[e(n),e(n-1),K,e(n-P+1)]^T；

計算出當前時刻的自適應濾波器的第i個抽頭權系數w_i(n)的步長g_i(n)，g_i(n)＝max{ρ,|w_i(n)-w_i(kM)|}；其中，max{}表示取最大值，|·|表示取絕對值，ρ為一個大于零的常數，其取值范圍為0.001～0.1；w_i(kM)為kM時刻的自適應濾波器的第i個抽頭權系數，M為一常數，其取值范圍為100～300，表示向下取整數；

將當前時刻的自適應濾波器的所有抽頭權系數w_i(n)的步長g_i(n)構成對角化的步長矩陣G(n)，G(n)＝diag{g₁(n),g₂(n),...,g_i(n)...,g_L(n)}；其中，diag{}表示形成對角化矩陣；

最后，按下式得出下一時刻的自適應濾波器抽頭權系數向量W(n+1)，

W(n+1)＝W(n)+μG(n)U(n)[U^T(n)G(n)U(n)+δI_P]^-1E(n)

其中，μ為自適應濾波器的全局步長參數，其取值范圍為0<μ<2，δ為防止矩陣求逆計算困難的正常數，其取值通常為0.01，I_P為P×P的單位矩陣；

E、令n＝n+1，重復步驟B、C、D，即可實現實時的回聲消除。