CN106170049B - 一种有偏补偿的归一化子带自适应回声消除方法 - Google Patents

一种有偏补偿的归一化子带自适应回声消除方法 Download PDF

Info

Publication number
CN106170049B
CN106170049B CN201610312814.8A CN201610312814A CN106170049B CN 106170049 B CN106170049 B CN 106170049B CN 201610312814 A CN201610312814 A CN 201610312814A CN 106170049 B CN106170049 B CN 106170049B
Authority
CN
China
Prior art keywords
sampling instant
moment
filter
sef
subband
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610312814.8A
Other languages
English (en)
Other versions
CN106170049A (zh
Inventor
赵海全
郑宗生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201610312814.8A priority Critical patent/CN106170049B/zh
Publication of CN106170049A publication Critical patent/CN106170049A/zh
Application granted granted Critical
Publication of CN106170049B publication Critical patent/CN106170049B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04M—TELEPHONIC COMMUNICATION
    • H04M9/00—Arrangements for interconnection not involving centralised switching
    • H04M9/08—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic
    • H04M9/082—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic using echo cancellers

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

一种有偏补偿的归一化子带自适应回声消除方法，其步骤如下：A、获取子带信号；B、自适应滤波，将采样时刻n＝kM自适应滤波器的第m个输入向量Xm(k)通过自适应滤波器得到采样时刻n＝kM的第m个输出值ym(k)，C、回声抵消，将第m个近端子带信号的第k时段抽取值与采样时刻n＝kM的第m个滤波器输出值ym(k)相减得到采样时刻n＝kM的第m个子带残差em(k)，D、滤波器抽头权向量更新，使用有偏补偿的归一化子带的方法，得出自适应滤波器在n＝(k+1)M时刻的抽头权向量w(k+1)：E、令k＝k+1,重复B、C、D的步骤,直至通话结束。该方法对通信系统的声学回声的消除效果好，收敛速度快，稳态误差小。

Description

一种有偏补偿的归一化子带自适应回声消除方法
技术领域
本发明涉及一种有偏补偿的归一化子带自适应回声消除方法，属于通信的回声对消技术领域。
背景技术
回声存在于各种通信系统中，凡是在需要同时使用扬声器和麦克风的场合，例如免提电话、移动电话、网络电话以及电视电话会议系统中都不可避免地存在回声问题。以电视电话会议为例，因为扬声器和麦克风被置于同一空间，本地扬声器发出的远端语音被本地麦克风接收并传回远端，导致远端说话者听到自己的声音，这就是回声。在通信中，人耳又对于这种回声及其敏感，延迟10ms的回声就能被人耳捕捉感知到，超过32ms的回声就会对通信质量造成极大的干扰。因此，必须采取有效的措施来消除回声信号，提高语音通话质量。
目前，研究人员提出了多种方法来消除回声。在众多回声消除方法中，自适应回声消除技术具有逐步调节性能，应用成本低，收敛速度快，回声残差小，得到了研究学者们的一致认可，是目前国际上公认的最有前景的回声消除技术，也是回声消除目前采用的主流技术。自适应滤波器是能够根据输入信号自主调整性能进行数字信号处理的数字滤波器，作为对比，传统的非自适应滤波器(如有限脉冲响应数字滤波器，无限脉冲响应数字滤波器)的参数是静态的预先设定不可调整的。对于一些未知的环境来说，自适应滤波器能够在工作过程中逐步学习出所需的统计特性，并以此为依据自动调整滤波器系数，以达到最佳滤波的效果；一旦输入信号的统计特性发生变化，它又可以跟踪这种变化，自动调整滤波器参数，使滤波性能重新达到最佳。
自适应回声消除技术就是利用合适的自适应辨识策略来消除回声，通过一个自适应滤波器对实际未知的回声路径进行系统辨识是它的核心思想，并选取合适的自适应算法使滤波器的权系数与未知回声路径的脉冲响应相逼近，于是我们就可以获取回声估计信号，接着从近端接收到的信号中减去远端所获取的回声估计信号，从而把回声消除掉。自适应回声消除技术的核心是自适应回声消除算法。因此，如何完善和研究新的性能卓越的自适应回声消除算法是回声消除领域的主要研究方向。
在目前的自适应回声消除应用中主要采用归一化子带回声消除方法，见参考文献1“Improving convergence of the NLMS algorithm using constrained subbandupdates”(K.A.Lee和W.S.Gan，IEEE Signal Process.Lett.,vol.11,no.9,pp.736–739,Sep.2004.)与传统的归一化最小均方方法相比，该方法能够获得更快的收敛速度。然而，该方法在自适应滤波器输入信号带噪声的情况下表现出较大的稳态误差。
发明内容
本发明的发明目的就是提供一种有偏补偿的归一化子带自适应回声消除方法，该方法对通信系统的声学回声的消除效果好，收敛速度快，稳态误差小。
本发明实现其发明目的所采用的技术方案是，一种有偏补偿的归一化子带自适应回声消除方法，其步骤如下：
一种有偏补偿的归一化子带自适应回声消除方法，其步骤如下：
A、获取子带信号
将远端传来的远端信号采样得到在采样时刻n的远端信号离散值x(n)，通过分析滤波器组得到采样时刻n的M个远端子带信号x1(n),x2(n),...xm(n),...xM(n)，其中，m为子带序号，xm(n)为第m个远端子带信号，M为子带数目，其取值为2、4、8、16、32；
将近端麦克风拾取的带回声的近端信号采样得到在采样时刻n的近端信号离散值d(n)，通过分析滤波器组得到采样时刻n的M个近端子带信号d1(n),d2(n),...dm(n)...,dM(n)；其中，dm(n)为第m个近端子带信号；
在采样时刻n＝kM时，将该时刻的第m个近端子带信号dm(n)＝dm(kM)记为第m个近端子带信号的第k时段抽取值其中，k为抽取时段序号；
B、自适应滤波
将采样时刻n＝kM及其前L-1个时刻的第m个远端子带信号xm(n)，构成自适应滤波器的第m个输入向量Xm(k)，Xm(k)＝[xm(kM)，xm(kM-1),...,xm(kM-L+1)]T，其中L为自适应滤波器抽头数，其取值为32、64、128、256、512、1024，T代表转置运算；
将采样时刻n＝kM的第m个输入向量Xm(k)通过自适应滤波器得到采样时刻n＝kM的第m个输出值ym(k)，其中，w(k)为采样时刻n＝kM的自适应滤波器的抽头权向量，w(k)＝[w1(k),w2(k),...wl(k)...,wL(k)]T，wl(k)为采样时刻n＝kM的自适应滤波器的第l个抽头权系数，w(k)的初始值为零向量；
C、回声抵消
将第m个近端子带信号的第k时段抽取值与采样时刻n＝kM的第m个滤波器输出值ym(k)相减得到采样时刻n＝kM的第m个子带残差em(k)，
将时刻n＝kM的各个子带残差em(k)通过综合滤波器组得到时刻n＝kM的残差e(kM)；并将时刻n＝kM的残差e(kM)传送给远端；在时刻n＝kM+M'，且M'＝1,2,3,...M-1时，令时刻n＝kM+M'的残差e(kM+M')＝e(kM)，并将时刻n＝kM+M'的残差e(kM+M')传送给远端；
D、滤波器抽头权向量更新
D1、计算有偏补偿向量
由下式计算有偏补偿向量c(k)：
其中，||·||2表示2范数，μ为自适应滤波器的步长因子，其取值范围为0～2，δ为有偏补偿因子，其取值范围为1×10-8～1；
D2、滤波器抽头权向量更新
使用有偏补偿的归一化子带的方法，得出自适应滤波器在n＝(k+1)M时刻的抽头权向量w(k+1)：
E、令k＝k+1,重复B、C、D的步骤,直至通话结束。
与现有技术相比，本发明的有益效果是：
相比于传统的归一化子带回声消除方法，本发明在滤波器抽头权向量更新过程中添加了有偏补偿向量，能够有效地降低由输入噪声信号带来的滤波器抽头权向量与预期值的偏差，显著降低了算法稳态误差，表现出较好的回声消除效果。
下面结合附图和具体实施方式对本发明进行详细说明
附图说明
图1是参考文献1的方法和本发明方法的仿真实验归一化稳态失调曲线。
具体实施方式
实施例
本发明的一种具体实施方式是，一种有偏补偿的归一化子带自适应回声消除方法，其步骤如下：
A、获取子带信号
将远端传来的远端信号采样得到在采样时刻n的远端信号离散值x(n)，通过分析滤波器组得到采样时刻n的M个远端子带信号x1(n),x2(n),...xm(n),...xM(n)，其中，m为子带序号，xm(n)为第m个远端子带信号，M为子带数目，其取值为2、4、8、16、32；
将近端麦克风拾取的带回声的近端信号采样得到在采样时刻n的近端信号离散值d(n)，通过分析滤波器组得到采样时刻n的M个近端子带信号d1(n),d2(n),...dm(n)...,dM(n)；其中，dm(n)为第m个近端子带信号；
在采样时刻n＝kM时，将该时刻的第m个近端子带信号dm(n)＝dm(kM)记为第m个近端子带信号的第k时段抽取值其中，k为抽取时段序号；
B、自适应滤波
将采样时刻n＝kM及其前L-1个时刻的第m个远端子带信号xm(n)，构成自适应滤波器的第m个输入向量Xm(k)，Xm(k)＝[xm(kM)，xm(kM-1),...,xm(kM-L+1)]T，其中L为自适应滤波器抽头数，其取值为32、64、128、256、512、1024，T代表转置运算；
将采样时刻n＝kM的第m个输入向量Xm(k)通过自适应滤波器得到采样时刻n＝kM的第m个输出值ym(k)，其中，w(k)为采样时刻n＝kM的自适应滤波器的抽头权向量，w(k)＝[w1(k),w2(k),...wl(k)...,wL(k)]T，wl(k)为采样时刻n＝kM的自适应滤波器的第l个抽头权系数，w(k)的初始值为零向量；
C、回声抵消
将第m个近端子带信号的第k时段抽取值与采样时刻n＝kM的第m个滤波器输出值ym(k)相减得到采样时刻n＝kM的第m个子带残差em(k)，
将时刻n＝kM的各个子带残差em(k)通过综合滤波器组得到时刻n＝kM的残差e(kM)；并将时刻n＝kM的残差e(kM)传送给远端；在时刻n＝kM+M'，且M'＝1,2,3,...M-1时，令时刻n＝kM+M'的残差e(kM+M')＝e(kM)，并将时刻n＝kM+M'的残差e(kM+M')传送给远端；
D、滤波器抽头权向量更新
D1、计算有偏补偿向量
由下式计算有偏补偿向量c(k)：
其中，||·||2表示2范数，μ为自适应滤波器的步长因子，其取值范围为0～2，δ为有偏补偿因子，其取值范围为1×10-8～1；
D2、滤波器抽头权向量更新
使用有偏补偿的归一化子带的方法，得出自适应滤波器在n＝(k+1)M时刻的抽头权向量w(k+1)：
E、令k＝k+1,重复B、C、D的步骤,直至通话结束。
仿真实验：
为了验证本发明的有效性，进行了仿真实验，并与参考文献1的方法进行对比。
仿真实验的远端信号x(n)采用极点为0.8的一阶自回归(AR(1))信号加信噪比为10dB的高斯白噪声，采样频率为8000Hz，采样点个数80000。回声信道脉冲响应在高2.5m，宽3.75m，长6.25m，温度20℃，湿度50％的安静密闭房间内获得，脉冲响应长度即滤波器抽头数L＝128。实验的背景噪声为高斯白噪声，信噪比为30dB。
按照以上实验条件，用参考文献1的方法和本发明方法进行回声消除实验。各种方法的参数具体取值如表1。
表1实验各算法的最优参数近似取值
仿真实验结果是通过独立运行50次平均得到的。图1是参考文献1的方法和本发明方法的仿真实验归一化稳态失调曲线。
从图1可以看出，在收敛速度相同的情况下，本发明比参考文献1的方法具有更好的稳态误差，稳态误差降低约4dB。因此，本发明具有更小的稳态误差，表现出更好的回声消除效果。

Claims (1)

1.一种有偏补偿的归一化子带自适应回声消除方法，其步骤如下：
A、获取子带信号
将远端传来的远端信号采样得到在采样时刻n的远端信号离散值x(n)，通过分析滤波器组得到采样时刻n的M个远端子带信号x1(n),x2(n),...xm(n),...xM(n)，其中，m为子带序号，xm(n)为第m个远端子带信号，M为子带数目，其取值为2、4、8、16、32；
将近端麦克风拾取的带回声的近端信号采样得到在采样时刻n的近端信号离散值d(n)，通过分析滤波器组得到采样时刻n的M个近端子带信号d1(n),d2(n),...dm(n)...,dM(n)；其中，dm(n)为第m个近端子带信号；
在采样时刻n＝kM时，将该时刻的第m个近端子带信号dm(n)＝dm(kM)记为第m个近端子带信号的第k时段抽取值其中，k为抽取时段序号；
B、自适应滤波
将采样时刻n＝kM及其前L-1个时刻的第m个远端子带信号xm(n)，构成自适应滤波器的第m个输入向量Xm(k)，Xm(k)＝[xm(kM)，xm(kM-1),...,xm(kM-L+1)]T，其中L为自适应滤波器抽头数，其取值为32、64、128、256、512、1024，T代表转置运算；
将采样时刻n＝kM的第m个输入向量Xm(k)通过自适应滤波器得到采样时刻n＝kM的第m个输出值其中，w(k)为采样时刻n＝kM的自适应滤波器的抽头权向量，w(k)＝[w1(k),w2(k),...w l (k)...,wL(k)]T，w l (k)为采样时刻n＝kM的自适应滤波器的第l个抽头权系数，w(k)的初始值为零向量；
C、回声抵消
将第m个近端子带信号的第k时段抽取值与采样时刻n＝kM的第m个滤波器输出值ym(k)相减得到采样时刻n＝kM的第m个子带残差em(k)，
将时刻n＝kM的各个子带残差em(k)通过综合滤波器组得到时刻n＝kM的残差e(kM)；并将时刻n＝kM的残差e(kM)传送给远端；在时刻n＝kM+M'，且M'＝1,2,3,...M-1时，令时刻n＝kM+M'的残差e(kM+M')＝e(kM)，并将时刻n＝kM+M'的残差e(kM+M')传送给远端；
D、滤波器抽头权向量更新
D1、计算有偏补偿向量
由下式计算有偏补偿向量c(k)：
其中，||·||2表示2范数，μ为自适应滤波器的步长因子，其取值范围为0～2，δ为有偏补偿因子，其取值范围为1×10-8～1；
D2、滤波器抽头权向量更新
使用有偏补偿的归一化子带的方法，得出自适应滤波器在n＝(k+1)M时刻的抽头权向量w(k+1)：
E、令k＝k+1,重复A、B、C、D的步骤,直至通话结束。
CN201610312814.8A 2016-05-12 2016-05-12 一种有偏补偿的归一化子带自适应回声消除方法 Expired - Fee Related CN106170049B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610312814.8A CN106170049B (zh) 2016-05-12 2016-05-12 一种有偏补偿的归一化子带自适应回声消除方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610312814.8A CN106170049B (zh) 2016-05-12 2016-05-12 一种有偏补偿的归一化子带自适应回声消除方法

Publications (2)

Publication Number Publication Date
CN106170049A CN106170049A (zh) 2016-11-30
CN106170049B true CN106170049B (zh) 2019-05-17

Family

ID=57358967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610312814.8A Expired - Fee Related CN106170049B (zh) 2016-05-12 2016-05-12 一种有偏补偿的归一化子带自适应回声消除方法

Country Status (1)

Country Link
CN (1) CN106170049B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109040498B (zh) * 2018-08-12 2022-01-07 瑞声科技（南京）有限公司 一种提升回声抵消效果的方法及其系统
CN109697986B (zh) * 2018-09-19 2020-12-18 四川大学 一种基于最小三次方绝对值的自适应偏差补偿回声消除方法
CN111739548B (zh) * 2020-06-22 2022-03-22 西南交通大学 抗输入噪声的最小误差熵自适应回声消除方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102195616A (zh) * 2010-03-02 2011-09-21 哈曼贝克自动系统股份有限公司 高效子带自适应有限脉冲响应滤波
CN102724152A (zh) * 2012-07-12 2012-10-10 西南交通大学 基于Laguerre结构的多项式自适应有源噪声对消方法
CN104185866A (zh) * 2012-03-30 2014-12-03 苹果公司 用于有源噪声消除自适应滤波器的预成形串联滤波器
CN105070295A (zh) * 2015-07-10 2015-11-18 西南交通大学 一种活性因子成比例子带的自适应方法应用于回声消除

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102195616A (zh) * 2010-03-02 2011-09-21 哈曼贝克自动系统股份有限公司 高效子带自适应有限脉冲响应滤波
CN104185866A (zh) * 2012-03-30 2014-12-03 苹果公司 用于有源噪声消除自适应滤波器的预成形串联滤波器
CN102724152A (zh) * 2012-07-12 2012-10-10 西南交通大学 基于Laguerre结构的多项式自适应有源噪声对消方法
CN105070295A (zh) * 2015-07-10 2015-11-18 西南交通大学 一种活性因子成比例子带的自适应方法应用于回声消除

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Affine projection M-estimate subband adaptive filters for robust adaptive filtering in impulsive noise;Zongsheng Zheng et al.;《SIGNAL PROCESSING》;20160331;64-70

Also Published As

Publication number Publication date
CN106170049A (zh) 2016-11-30

Similar Documents

Publication Publication Date Title
CN109040499B (zh) 一种抗冲击干扰的自适应回声消除方法
US8934620B2 (en) Acoustic echo cancellation for high noise and excessive double talk
Halimeh et al. Combining adaptive filtering and complex-valued deep postfiltering for acoustic echo cancellation
CN104506746B (zh) 一种改进的凸组合解相关成比例自适应回声消除方法
CN106210370B (zh) 一种仿射投影最大熵子带自适应回声消除方法
CN103369162B (zh) 一种低复杂度的电话回声自适应消除方法
CN104410761B (zh) 一种仿射投影符号子带凸组合自适应回声消除方法
CN110138990A (zh) 一种消除移动设备VoIP电话回声的方法
CN106170049B (zh) 一种有偏补偿的归一化子带自适应回声消除方法
US9020144B1 (en) Cross-domain processing for noise and echo suppression
EP2101480A3 (en) Echo canceller and echo cancelling method
CN107369458B (zh) 一种基于熵吸引的最大熵子带回声消除方法
CN110956975A (zh) 回声消除方法及装置
Cheng et al. Semi-blind source separation for nonlinear acoustic echo cancellation
CN106412352B (zh) 一种时频记忆子带比例自适应回声消除方法
CN106128471B (zh) 一种收缩变步长子带回声消除方法
CN106157965B (zh) 一种基于权向量重用的零范数集员仿射投影自适应回声消除方法
CN109697986A (zh) 一种基于最小三次方绝对值的自适应偏差补偿回声消除方法
CN106161822B (zh) 一种基于m估计的归一化子带自适应回声消除方法
CN111756943B (zh) 一种增广复值成比例子带自适应立体声学回声消除方法
CN107527622B (zh) 一种基于最大熵的变步长符号子带回声消除方法
CN106060295B (zh) 一种凸组合系数差的成比例仿射投影回声消除方法
CN107578782B (zh) 一种变核宽度的最大熵子带回声消除方法
CN113409806A (zh) 一种基于反正切函数的零吸引回声消除方法
Tiwari et al. A Review Paper on Adaptive Noise Cancellation Implementation using TMS320C6713 DSP Board

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190517