基于块稀疏成比例重用权系数仿射投影的回声消除方法
【技术领域】
[0001 ]本发明属于通信的自适应回声对消技术领域。
【背景技术】
[0002] 随着科技的发展，多功能多种类的通信应用一直层出不穷，用户在享受着便利的 同时也会遇到通信中的回声问题所带来的困扰。根据回声产生的机理，回声可以分为声学 回声和电学回声。电学回声通常比较小，因此，当前回声消除研究的重点已转向了声学回声 消除。通信过程中，声学回声产生的机理是远端说话人的声音在近端扬声器播放出来后又 被麦克风接受又传播到远端的说话人，使之又听到了自己的声音，从而形成了声学上的回 声。声学回声会让通信用户觉得非常不适应，严重时还会影响通话质量。因此，如何采取有 效的方法来消除声学回声成为了信号处理领域热门课题之一。回声消除技术的难点主要是 由于要消除的电声反馈信号中有大量通过室内反射声途径产生的电声反馈，它们是无限多 的、连续不断的、无规则延时的、来自各个方向的反射声信号。一般说来，在已投入使用的回 声抑制器中或仍在研究、开发的相关技术中，均采取了适当的自适应辨识方案来消除回声。
[0003] 在处理语音信号时，仿射投影技术(APA)算法可看作归一化最小均方算法(NLMS) 算法的时域扩展，仿射投影技术(APA)比传统的NLMS算法具有更快的收敛速度，但APA的稳 态误差性能变差，在高噪声下，APA的估计波动极大。针对上述问题，Clark和Shynk均提出对 输入数据序列通过串-并变换器被分成L点的块，而且这样产生的输入数据块被一次一块地 加到长度为M的有限冲击响应(FIR)滤波器。在收集到每一块数据样值后，进行滤波器抽头 权值的更新，使得滤波器的自适应一块一块地进行，而不是像传统LMS滤波器那样一个样值 的进行，减少计算复杂性。基于此基础上，之后又有许多新的有效算法被提出。但是，通信中 的大多数回声信道都是稀疏信道，这种稀疏信道的脉冲响应大多很小，而阶数又比较长。传 统自适应滤波器在这种稀疏信道下的回声消除收敛慢，稳态误差大。针对这种情况引入成 比例滤波器，可以提高算法的收敛速度。目前，回声消除应用中，性能较好的算法有：
[0004] (1)基于重用权系数的仿射投影回声消除方法
[0005] 参考文献"A noise-resilient affine projection algorithm and its convergence analysis'，（Seong-Eun Kim，Signal Processing,vol .121 ,pp.94-101,2016) 该方法将重用权系数方法引入仿射投影算法，提高了自适应滤波器的收敛速度和跟踪性 能。但是，由于在稀疏系统下该算法没有利用稀疏系统的特性，因此仍表现出收敛速度慢的 缺陷。
[0006] (2)基于块稀疏的成比例仿射投影回声消除方法
[0007] 参考文南犬"Proportionate Adaptive Filtering for Block-sparse System Identification"（IEEE/ACM Transactions on Audio,Speech,and Language Processing，Vol .24,pp.623-630,April .2016)该算法将块稀疏算法引入仿射投影算法。在 块稀疏系统中，该方法具有较快的收敛速度和较小的稳态误差。然而，在低信噪比系统下， 由于没有使用重用权系数方法去消除高噪声对系统的影响，该方法的收敛速度仍然较慢。
【发明内容】
[0008] 本发明的目的就是提供一种基于块稀疏成比例的重用权系数仿射投影方法。该方 法对回声消除效果好，收敛速度快，稳态误差小。
[0009] 本发明实现其发明目的所采用的技术方案是，一种基于块稀疏成比例的重用权系 数仿射投影方法，其步骤如下：
[0010] A、远端信号滤波
[0011] A1、将远端传来的远端信号采样得到当前时刻n的远端信号离散值x(n)，远端信号 离散值x(n)在n到n-L+1时刻的值x(n)，x(n-l)，. . .，x(n-L+l)构成当前时刻n的滤波器输入 向量父(1〇，父(11)=[以11)，以11-1)，...，奴11-1+1)]7，其中1 = 512是滤波器抽头数，1'代表转置 运算；
[0012] A2、将滤波器当前时刻n的输入向量X(n)通过自适应滤波器得到当前时刻n的滤波 器输出值y(n)，y(n) =XT(n)w(n);其中w(n)为自适应回声消除滤波器当前时刻n的抽头权 系数向量，w(n) = [wi(n)，W2(n)，. . .wi(n)，. . .，WL(n)]T，其长度为L,初始值为零向量，wi(n) 为抽头权系数向量w(n)中的第1个元素；
[0013]将n到n-M+l时刻的滤波器输出值y(n)进行平滑，得到当前时刻n的滤波器输出平 滑值
为重用权向量的个数，其取值为2,4,6;
[0014] B、回声抵消
[0015] 将近端麦克风拾取到的当前时刻n的带回声的近端信号d(n)与当前时刻n的滤波 器输出平滑值灭《)的差值，作为当前时刻n的残差信号，^(?)=外《)-只》)，并送回给 远端；
[0016] C、滤波器抽头权向量更新
[0017] C1、残差向量的组成:将时刻n到时刻n-K+1的残差信号一),一--1)，..., 一-[ + 1) 组成当前时刻n的残差向量豆(")它(/2) = f⑷；(/2- 1).....~e{n-K + \)\ ;
[0018] C2、欧几里德2，1范数的计算:将A2步的当前时刻n的抽头权系数向量w(n)均分成I = L/P个分组，每组的抽头权系数个数为P，得到当前时刻n第i组抽头权系数向量W[1](n)， w[i](n) = [w((i-i)p+i)，w((i-i)p+2), ? ? .w((i-i)p+P), ? ? ?，w((i-i)p+p)]，进而得到当前时刻n第i组抽 头权系数向量w[i](n)的欧几里德2，1范数| |w[i](n) | h
其中，p 为第i组抽头权系数向量中的元素 w((i-：〇p+P)的序号，p=l，2,. . .，P;
[0019] C3、欧几里德范数的修定:求出各组抽头权系数向量W[1](n)的欧几里德2，1范数| wwU) | |2与抽头权系数设定阈值q间的最大值，将该最大值与设定加权系数P的积作为当 前时刻n的范数下限值Mn)，人(n)=pmax{q, | |w[i](n) | 丨2, | |w[2](n) | 丨2,…，| |w[i](n) ，...，||?[1](11)||2}，其中，11^表示取最大值运算，加权系数0的取值为〇.〇1-〇.〇5;
[0020] 再由下式得出当前时刻n第i组抽头权系数向量W[1](n)的欧几里德2，1范数修定值 T i(n),
[0021] y i(n)=max{A(n)，| |w[i](n) | I2}，
[0022] C4、成比例矩阵的计算：
[0023] 由下式算出当前时亥Ijn第i组抽头权系数向量w[i](n)的范数归一化值y 'i(n)，
[0025] 将当前时刻n各组抽头权系数向量W[1](n)的范数归一化值由下式扩展组成抽头权 系数向量w[i](n)的范数归一化值矩阵r (n)，
[0026] r (n) = [y 'i(n) ? lp，y，2(n) ? lp，.. y 'i(n) ? lp.，y 'i(n) ? lp]T，式中，lp 表示 长度为P的元素全为一的行向量，?表示向量的点乘；
[0027] 再由当前时刻n的抽头权系数向量w[i](n)的范数归一化值矩阵r (n)和当前时刻n 的滤波器输入向量X(n)，算出当前时刻n的成比例矩阵G(n)，
[0028] G(n) = [ r (n) 0 X(n) ,G-i(n)]
[0029] 其中，〇表示两个矩阵的各个元素相乘，G-Kn)表示n-1时刻成比例矩阵的前K-l 列，初始时为零矩阵；
[0030] C5、滤波器抽头权向量更新
[0031]将A1步的当前时刻n和前K-1个时刻的滤波器输入向量X(n)，X(n-l)，. . .，X(n-K+ 1)组合，得到滤波器仿射投影输入向量1](1〇，1](11)=[父(11)，父(11-1)...父(11-1(+1)]，其中1(代 表仿射投影阶数，K = 4,8，16;
[0032]计算下一时刻n+1的自适应回声消除滤波器的抽头权系数向量w(n+l);
[0034]其中，y为自适应滤波器的步长，其取值范围为0<y<2，S为防止矩阵求逆计算困 难的常数，其取值为0.001~0.01，Ik为K XK的单位矩阵；
[0035] D、令n = n+1，重复A、B、C的步骤，直至通话结束。
[0036] 与现有技术相比，本发明的有益效果是：
[0037] -、本发明引入重用权系数仿射投影方法，通过对前M个时刻的滤波器输出值进行 平滑处理作为当前时刻滤波器输出值，对相关信号的解相关作用强，增强了对信号的抗干 扰能力，其收敛速度快，稳态误差小，回声消除效果好。
[0038] 二、权系数更新时，通过对当前滤波器抽头权系数向量进行分组，即分块投影计算 各块抽头权系数向量的当前时刻成比例值，再和当前时刻滤波器输出值相乘，修正得到下 一时刻滤波器抽头权系数向量更新值。通过这样块稀疏成比例仿射投影算法，使滤波器能 够更快速、准确的跟踪非零信号集中成块(块稀疏)的声音通信信号，也使其收敛速度快，稳 态误差小，回声消除效果好。
[0039]下面结合附图和【具体实施方式】对本发明进行详细说明
【附图说明】
[0040] 图1是本发明仿真实验的块稀疏信道图。
[0041] 图2和图3分别是文献1、文献2和本发明在一阶相关信号为输入信号时仿真实验的 归一化稳态失调曲线。 实施例
[0042]本发明的一种【具体实施方式】是，一种基于块稀疏成比例重用权系数仿射投影的回 声消除方法，其步骤如下：
[0043] A、远端信号滤波
[0044] A1、将远端传来的远端信号采样得到当前时刻n的远端信号离散值x(n)，远端信号 离散值x(n)在n到n-L+1时刻的值x(n)，x(n-l)，. . .，x(n-L+l)构成当前时刻n的滤波器输入 向量父(1〇，父(11)=[以11)，以11-1)，...，奴11-1+1)]7，其中1 = 512是滤波器抽头数，1'代表转置 运算；
[0045] A2、将滤波器当前时刻n的输入向量X(n)通过自适应滤波器得到当前时刻n的滤波 器输出值y(n)，y(n) =XT(n)w(n);其中w(n)为自适应回声消除滤波器当前时刻n的抽头权 系数向量，w(n) =