技術(shù)特征：

1.一種基于M估計(jì)的歸一化子帶自適應(yīng)回聲消除方法，其步驟如下：

A、獲取子帶信號(hào)

將遠(yuǎn)端傳來的遠(yuǎn)端信號(hào)采樣得到在采樣時(shí)刻n的遠(yuǎn)端信號(hào)離散值x(n)，通過分析濾波器組得到采樣時(shí)刻n的M個(gè)遠(yuǎn)端子帶信號(hào)x₁(n),x₂(n),...x_m(n),...x_M(n)，其中，m為子帶序號(hào)，x_m(n)為第m個(gè)遠(yuǎn)端子帶信號(hào)，M為子帶數(shù)目，其取值為2、4、8、16、32；

將近端麥克風(fēng)拾取的帶回聲的近端信號(hào)采樣得到在采樣時(shí)刻n的近端信號(hào)離散值d(n)，通過分析濾波器組得到采樣時(shí)刻n的M個(gè)近端子帶信號(hào)d₁(n),d₂(n),...d_m(n)...,d_M(n)；其中，d_m(n)為第m個(gè)近端子帶信號(hào)；

在采樣時(shí)刻n＝kM時(shí)，將該時(shí)刻的第m個(gè)近端子帶信號(hào)d_m(n)＝d_m(kM)記為第m個(gè)近端子帶信號(hào)的第k時(shí)段抽取值其中，k為抽取時(shí)段序號(hào)；

B、自適應(yīng)濾波

將采樣時(shí)刻n＝kM及其前L-1個(gè)時(shí)刻的第m個(gè)遠(yuǎn)端子帶信號(hào)x_m(n)，構(gòu)成自適應(yīng)濾波器的第m個(gè)輸入向量X_m(k)，X_m(k)＝[x_m(kM)，x_m(kM-1),...,x_m(kM-L+1)]^T，其中L為自適應(yīng)濾波器抽頭數(shù)，其取值為32、64、128、256、512、1024，T代表轉(zhuǎn)置運(yùn)算；

將采樣時(shí)刻n＝kM的第m個(gè)輸入向量X_m(k)通過自適應(yīng)濾波器得到采樣時(shí)刻n＝kM的第m個(gè)輸出值y_m(k)，其中，w(k)為采樣時(shí)刻n＝kM的自適應(yīng)濾波器的抽頭權(quán)向量，w(k)＝[w₁(k),w₂(k),…w_l(k)...,w_L(k)]^T，w_l(k)為采樣時(shí)刻n＝kM的自適應(yīng)濾波器的第l個(gè)抽頭權(quán)系數(shù)，w(k)的初始值為零向量；

C、回聲抵消

將第m個(gè)近端子帶信號(hào)的第k時(shí)段抽取值與采樣時(shí)刻n＝kM的第m個(gè)濾波器輸出值y_m(k)相減得到采樣時(shí)刻n＝kM的第m個(gè)子帶殘差e_m(k)，

將時(shí)刻n＝kM的各個(gè)子帶殘差e_m(k)通過綜合濾波器組得到時(shí)刻n＝kM的殘差e(kM)；并將時(shí)刻n＝kM的殘差e(kM)傳送給遠(yuǎn)端；在時(shí)刻n＝kM+M'；時(shí)，令時(shí)刻n＝kM+M'的殘差e(kM+M')＝e(kM)，并將時(shí)刻n＝kM+M'的殘差e(kM+M')傳送給遠(yuǎn)端；其中M'＝1,2,3,…M-1；

D、濾波器抽頭權(quán)向量更新

D1、子帶殘差方差的計(jì)算

由下式算出第m個(gè)子帶殘差的方差

${\mathrm{\σ}}_m^2\left(k\right)={\mathrm{\λ\σ}}_m^2(k-1)+2.13(1-\mathrm{\λ})med(\[e_m^2(k),e_m^2(k-1),...,e_m^2(k-K_w+1)\])$

其中，λ為遺忘因子，其取值范圍為0.800～0.999，函數(shù)med(·)表示取向量中間值，K_w為時(shí)間段估計(jì)窗的值，其取值范圍為2～20，子帶殘差的初始值為零；

D2、計(jì)算M估計(jì)核函數(shù)

由下式計(jì)算出第m個(gè)子帶殘差的M估計(jì)核函數(shù)值ψ(e_m(k))：

$\mathrm{\ψ}\left(e_m\right(k\left)\right)=\left\{\begin{array}{cc}{e}_m\left(k\right),& \left|e_m\left(k\right)\right|\<2.576{\mathrm{\σ}}_m\left(k\right)\\ 0,& \left|e_m\left(k\right)\right|\≥2.576{\mathrm{\σ}}_m\left(k\right)\end{array}\right.$

D3、濾波器抽頭權(quán)向量更新

使用基于M估計(jì)的歸一化子帶的方法，得出自適應(yīng)濾波器在n＝(k+1)M時(shí)刻的抽頭權(quán)向量w(k+1)：

$w(k+1)=w\left(k\right)+\mathrm{\μ}\underset{m=1}{\overset{M}{\Σ}}\frac{X_m\left(k\right)}{\left|\right|X_m\left(k\right)|{|}_2^2}\mathrm{\ψ}\left(e_m\right(k\left)\right)$

其中，||·||₂表示2范數(shù)，μ為自適應(yīng)濾波器的步長因子，其取值范圍為(0，2]；

E、令k＝k+1,重復(fù)A、B、C、D步的操作,直至通話結(jié)束。