[發明專利]一種重用權系數指數型誤差自適應回聲消除方法

有效

• 申請號: 201710362653.8
• 申請日: 2017-05-22
• 公開（公告）號: CN107333017B
• 公開（公告）日: 2019-08-20
• 發明（設計）人: 趙海全;羅正延
• 申請（專利權）人: 西南交通大學
• 主分類號: H04M9/08
• 分類號: H04M9/08;G10L21/0208
• 代理公司: 成都博通專利事務所 51208
• 代理人: 陳樹明

1.一種重用權系數指數型誤差自適應回聲消除方法，其步驟如下：

A、遠端信號采樣

將當前時刻n到前L-1個時刻的遠端采樣信號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為自適應濾波器的抽頭長度，其取值為16、32、128；

B、回聲信號估計

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

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

C、回聲信號消除

將當前時刻n的近端采樣信號d(n)，減去步驟B獲得的當前時刻n的回聲信號的估計值得到已消除回聲的當前時刻n的有用信號s(n)，

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

D1、p范數限制的濾波器抽頭權系數向量計算

由當前時刻n的自適應濾波器的第l個抽頭權系數w_l(n)，得到當前時刻n第l個抽頭權系數的p范數F_l(n)，F_l(n)＝(1+β)exp(-β|w_l(n)|)sgn[w_l(n)]；其中β為范數的冪次，其取值為5；exp(·)為自然對數的指數運算，|·|為絕對值運算，sgn(·)為符號運算；

將當前時刻n的所有抽頭權系數的p范數F₁(n),F₂(n),...F_l(n)...,F_L(n)，構成當前時刻n的抽頭權系數向量的p范數F(n)，F(n)＝[F₁(n),F₂(n),...F_l(n)...,F_L(n)]^T；

D2、重用權系數向量的計算

算出當前時刻n的第l個重用權系數S_l(n)，S_l(n)＝sgn[w_l(n)]/(1+ε|w_l(n)|)；其中，ε為重用權系數值的閾值，其取值為10；

將當前時刻n的所有的重用權系數S₁(n),S₂(n),...S_l(n)...,S_L(n)，構成當前時刻n的重用權系數向量S(n)，S(n)＝[S₁(n),S₂(n),...S_l(n)...,S_L(n)]^T；

D3、抽頭權系數向量的更新

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

其中，α為誤差陡度控制參數、其取值范圍為(0，1)，δ為正則化因子、其取值為0.0001，μ為步長因子，其取值范圍為(0，2)；γ是抽頭權系數的p范數的權值、其取值范圍為(0，1)，γ'是重用權系數向量的權值，取值范圍為(0，1)；

E、重復

令n＝n+1，重復步驟A、B、C、D的操作，直至通話結束。