[發明專利]一種基于熵吸引的最大熵子帶回聲消除方法

有效

• 申請號: 201710551228.3
• 申請日: 2017-07-07
• 公開（公告）號: CN107369458B
• 公開（公告）日: 2019-09-24
• 發明（設計）人: 趙海全;劉倩倩
• 申請（專利權）人: 西南交通大學
• 主分類號: G10L21/0208
• 分類號: G10L21/0208;G10L21/0264;H04M9/08
• 代理公司: 成都博通專利事務所 51208
• 代理人: 陳樹明

1.一種基于熵吸引的最大熵子帶回聲消除方法，其步驟如下：

A、信號的采樣與處理

將當前時刻n到時刻n-L+1之間的采樣遠端信號u(n),u(n-1),...,u(n-L+1),構成當前時刻n分析濾波器一的輸入向量U(n)，U(n)＝[u(n),u(n-1),...,u(n-L+1)]^T；L＝512是濾波器抽頭數，上標T表示轉置運算；

分析濾波器一將輸入向量U(n)按頻帶分割成I個遠端子帶向量U_i(n)，U_i(n)＝[u_i(n),u_i(n-1),...,u_i(n-L+1)]^T；

同時，分析濾波器二將近端麥克風拾取到的當前時刻n的帶回聲的近端信號d(n)按頻帶分割分割成I個近端子帶信號d_i(n)；

其中，i為遠端子帶向量或近端子帶信號的序號，i＝1,2,...,I，I為遠端子帶向量和近端子帶信號的總個數,其取值為2、4、6、8；

B、信號抽取

將遠端子帶向量U_i(n)經抽取器進行I抽取，即將n＝k＝KI時刻的遠端子帶向量U_i(n)抽出，得到抽取時刻k的遠端子帶抽取向量U_i(k),U_i(k)＝[u_i(k),u_i(k-1),...,u_i(k-L+1)]^T；其中，K為抽取的序號；

對近端子帶信號d_i(n)也經抽取器進行I抽取，即將n＝k＝KI時刻的近端子帶信號d_i(n)抽出，得到抽取時刻k的近端子帶抽取信號d_i(k)，d_i(k)＝d_i(KI)；

C、濾波器的輸出

將抽取時刻k的遠端子帶抽取向量U_i(k)通過自適應回聲消除濾波器中的子帶濾波器得到抽取時刻k的輸出子帶信號y_i(k)，其中W(k)為子帶濾波器在抽取時刻k的權系數向量，W(k)＝[w₁(k),w₂(k),..w_l(k).,w_L(k)]^T；w_l(k)為權系數向量W(k)中的第l個權系數，l＝1,2,...,L為權系數w_l(k)的序號；W(k)的初始值為零，即W(1)＝0；

D、回聲抵消

將抽取時刻k的近端子帶抽取信號d_i(k)與抽取時刻k的輸出子帶信號y_i(k)相減得到抽取時刻k的子帶誤差信號e_i(k)，即抽取時刻k消除回聲后的子帶回送信號e_i(k)，e_i(k)＝d_i(k)-y_i(k)，并將抽取時刻k的子帶誤差信號e_i(k)送回給遠端；

E、權系數向量的更新

E1、根據抽取時刻k的權向量W(k)，計算抽取時刻k的熵吸引量P(k)，

其中表示熵吸引的核寬度，取值為0.2；exp(·)表示以自然常數為底的指數運算；表示向量或矩陣的對應元素相乘；

E2、更新得到下一個抽取時刻k+1的權系數向量W(k+1)，

其中μ為步長參數，取值為0.1；γ為零吸引因子，取值為0.000001；

F、重復

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