[发明专利]一种重用权系数指数型误差自适应回声消除方法

摘要

本发明公开了一种重用权系数指数型误差自适应回声消除方法，其步骤主要是：A、远端信号采样，B、回声信号估计，C、回声信号消除，D、滤波器抽头权系数更新：先算出当前时刻n第l个抽头权系数的p范数F_l(n)，F_l(n)＝(1+β)exp(‑β|w_l(n)|)sgn[w_l(n)]；进而得到当前时刻n的抽头权系数向量的p范数F(n)，F(n)＝[F₁(n),F₂(n),...F_l(n)...,F_L(n)]^T；再算出当前时刻n的第l个重用权系数S_l(n)，S_l(n)＝sgn[w_l(n)]/(1+ε|w_l(n)|)；最后，得出下一时刻n+1的自适应滤波器抽头权系数向量W(n+1)，E、重复。该方法能够获得更快的收敛性，回声消除效果更好，且抗冲激能力强。

主权项

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的操作，直至通话结束。