1.一种零吸引变核宽度的最大熵回声消除方法,其步骤如下:A、获取信号
将远端传来的语音信号采样得到当前时刻n的远端信号离散值x(n);同时,对近端麦克风收集到的回声信号采样得到当前时刻n的期望信号d(n)。
B、计算自适应滤波器的输出y(n)
将远端信号离散值x(n)在当前时刻n到时刻n-L+1的值构成当前时刻n的自适应滤波器输入向量x(n),x(n)=[x(n),x(n-1)...,x(n-L+1)]T,其中,L表示自适应滤波器抽头数,L=512或1024,上标T代表转置;
计算当前时刻n的自适应滤波器的输出信号y(n),y(n)=WT(n)x(n),其中,W(n)为当前时刻n的自适应滤波器的抽头权向量,其长度等于L,初始值为零向量即W(1)=0;
C、回声抵消
将当前时刻n的期望信号d(n)减去当前时刻n的自适应滤波器的输出信号y(n),得到当前时刻n的误差信号e(n)作为当前时刻n消除回声后的近端信号再传送给远端,即e(n)=d(n)-y(n);
D、权系数向量的更新
D1、根据当前时刻n的误差信号e(n),算出当前时刻n去冲激干扰的误差信号其中τ1表示误差信号估计量的平滑参数,取值为0.2;min(·)表示取最小值运算;
的初始值为零,即 Nw为误差信号平方平滑估计的平滑窗口的大小,取值为10~20;
D2、计算当前时刻n的核宽度σn:
如果 则当前时刻n的核宽度σn等于默认核宽度σ0,即σn=σ0,其中kσ为限定常数,取值为20,默认核宽度σ0的取值为5~20;
如果 则当前时刻n的核宽度σn,由下式得出,
D3、更新得到下一个时刻n+1的抽头权向量W(n+1),其中,μ为步长参数,取值为0.1;exp[·]表示指数运算;sgn[·]表示取符号运算;ρ为零吸引因子,取值为0.000001;ε为正参数,取值为1~20;
E、令n=n+1,重复A、B、C、D的步骤,直至通话结束。