1.一种基于相关熵诱导的集员的自适应回声消除方法,其步骤是:A、远端信号采集
对远端传来的信号进行采样,获得当前时刻n的远端输入信号的离散值x(n);将当前时刻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=512代表滤波器抽头数;
B、回声信号估计
将当前时刻n的输入信号向量x(n)通过自适应滤波器,得到自适应滤波器的当前时刻n的输出值,即回声信号的估计值y(n)y(n)=xT(n)w(n)
其中w(n)为自适应滤波器的当前时刻n的抽头权向量,w(n)=[w1(n),w2(n),...,wL-1T(n)],w(n)的初始值为零向量;
C、回声消除
对近端麦克风采样得到带回声的当前时刻n的近端信号d(n),将其减去回声信号的估计值y(n),得到当前时刻n的误差信号e(n),再送回给远端,e(n)=d(n)-y(n);
D、滤波器抽头权向量更新
D1、计算输入信号仿射投影矩阵
将当前时刻n与前P-1个时刻的输入向量x(n),x(n-1),...,x(n-P+1)构成当前时刻n的输入信号仿射投影矩阵U(n),U(n)=[x(n),x(n-1),...,x(n-P+1)];其中P表示仿射投影阶数,取值范围为2~9;
D2、计算误差信号向量
将当前时刻n与前P-1个时刻的误差信号e(n),e(n-1),...,e(n-P+1)构成当前时刻n误差信号向量E(n),E(n)=[e(n),e(n-1),...,e(n-P+1)]T;
D3、计算步长因子
由下式算出当前时刻n的步长因子μ(n):其中γ表示误差阈值参数,取值范围为0.0001~1;
D4、计算相关熵诱导因子
由自适应滤波器的当前时刻n的抽头权向量w(n),计算出当前时刻n的相关熵因子C(n):其中ρ为控制因数,取值范围为10-8~1;σ为核宽度,范围为0.001~2;exp(·)表示指数运算;
D5、滤波器抽头权向量更新
由下式得出下一时刻n+1的滤波器抽头权向量w(n+1):w(n+1)=w(n)+μ(n)U(n)(UT(n)U(n))-1E(n)-C(n)E、令n=n+1,重复上述A、B、C、D的过程,直至通话结束。