技術(shù)特征：

1.一種基于權(quán)向量重用的零范數(shù)集員仿射投影自適應(yīng)回聲消除方法，其步驟如下：

A、遠(yuǎn)端信號(hào)采樣

將遠(yuǎn)端傳來(lái)的遠(yuǎn)端信號(hào)在當(dāng)前時(shí)刻n及之前L-1個(gè)時(shí)刻采樣得到的遠(yuǎn)端信號(hào)x(n)，x(n-1),...,x(n-L+1)構(gòu)成當(dāng)前時(shí)刻n的輸入向量X(n)，X(n)＝[x(n)，x(n-1),...,x(n-L+1)]^T，其中L＝512是濾波器抽頭數(shù)，T代表轉(zhuǎn)置運(yùn)算；

B、回聲估計(jì)

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

C、回聲抵消

將近端麥克風(fēng)拾取的帶回聲的當(dāng)前時(shí)刻n的近端信號(hào)d(n)，與自適應(yīng)濾波器的當(dāng)前時(shí)刻n的輸出值y(n)相減得到當(dāng)前時(shí)刻n的殘差信號(hào)e(n)，e(n)＝d(n)-y(n)；再將當(dāng)前時(shí)刻n的殘差信號(hào)e(n)回送給遠(yuǎn)端；

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

D1、計(jì)算相干權(quán)向量

由下式計(jì)算相干權(quán)向量

$\stackrel{\‾}{w}\left(n\right)=\frac{1}{M}\underset{m=0}{\overset{M-1}{\Σ}}w(n-m)$

其中，m為權(quán)向量相干時(shí)刻數(shù)，M為權(quán)向量最大相干時(shí)刻數(shù)，其取值范圍為2～5；

D2、計(jì)算輸入加權(quán)矩陣

將當(dāng)前時(shí)刻n及之前P-1個(gè)時(shí)刻的輸入向量X(n),X(n-1),...,X(n-P+1)構(gòu)成當(dāng)前時(shí)刻n的仿射投影輸入矩陣U(n)，U(n)＝[X(n),X(n-1),...,X(n-P+1)]；其中，P為仿射投影階數(shù)，其取值范圍為2～8；

由下式計(jì)算輸入加權(quán)矩陣Φ(n)：

Φ(n)＝[U^T(n)U(n)+δI_P]^-1

其中，δ為防止矩陣求逆計(jì)算困難的正常數(shù)，其取值為0.001，I_P為P×P的單位矩陣；

D3、計(jì)算相干殘差

由下式計(jì)算當(dāng)前時(shí)刻n的相干殘差

$\stackrel{\‾}{e}\left(n\right)=\mathrm{\θ}\stackrel{\‾}{e}(n-1)+(1-\mathrm{\θ})|d\left(n\right)-X^T\left(n\right)\stackrel{\‾}{w}\left(n\right)|$

其中，θ為平滑因子，其取值范圍為0.800～0.999，相干殘差初始值為零；

D4、計(jì)算步長(zhǎng)因子

由下式計(jì)算步長(zhǎng)因子μ(n)：

$\mathrm{\μ}\left(n\right)=\left\{\begin{array}{cc}1-\frac{\mathrm{\γ}}{\stackrel{\‾}{e}},& |d\left(n\right)-X^T\left(n\right)\stackrel{\‾}{w}\left(n\right)|\>\mathrm{\γ}\\ 0,& |d\left(n\right)-X^T\left(n\right)\stackrel{\‾}{w}\left(n\right)|\≤\mathrm{\γ}\end{array}\right.$

其中，γ為閾值參數(shù)，其取值范圍為0.0001～1；

D5、獲取零吸引向量

計(jì)算當(dāng)前時(shí)刻n的自適應(yīng)濾波器的第l個(gè)抽頭權(quán)系數(shù)的零吸引因子f_β(w_l(n))：

$f_{\mathrm{\β}}\left(w_l\right(n\left)\right)=\left\{\begin{array}{cc}-{\mathrm{\β}}^2{w}_l\left(n\right)-\mathrm{\β},& -\frac{1}{\mathrm{\β}}\≤w_l\left(n\right)\≤0\\ -{\mathrm{\β}}^2{w}_l\left(n\right)+\mathrm{\β},& 0\≤w_l\left(n\right)\≤\frac{1}{\mathrm{\β}}\\ 0,& \left|w_l\left(n\right)\right|\>\frac{1}{\mathrm{\β}}\end{array}\right.,$

其中，β為控制因子，其取值范圍為1～20；

將所有的抽頭權(quán)系數(shù)的零吸引因子f_β(w_l(n))構(gòu)成自適應(yīng)濾波器當(dāng)前時(shí)刻n的零吸引向量F_β(w(n))，F(xiàn)_β(w(n))＝[f_β(w₁(n)),f_β(w₂(n)),...,f_β(w_L(n))]^T，

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

將當(dāng)前時(shí)刻n及之前P-1個(gè)時(shí)刻的近端信號(hào)d(n),d(n-1),...,d(n-P+1)構(gòu)成當(dāng)前時(shí)刻n的近端信號(hào)向量D(n)，D(n)＝[d(n),d(n-1),...,d(n-P+1)]^T；

使用基于權(quán)向量重用的零范數(shù)集員仿射投影的方法，得出下一時(shí)刻n+1的自適應(yīng)濾波器的抽頭權(quán)向量w(n+1)：

$w(n+1)=\stackrel{\‾}{w}\left(n\right)+\mathrm{\μ}\left(n\right)U\left(n\right)\mathrm{\Φ}\left(n\right)\[D\left(n\right)-U^T\left(n\right)\stackrel{\‾}{w}\left(n\right)\]+\frac{\mathrm{\α}}{2M}\[U\left(n\right)\mathrm{\Φ}\left(n\right)U^T\left(n\right)-I_L\]F_{\mathrm{\β}}\left(w\right(n\left)\right)$

其中，α為零吸引參數(shù)，其取值范圍為1×10^-8～1，I_L為L(zhǎng)×L的單位矩陣；

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