CN101917527A - 回声消除的方法和装置 - Google Patents

回声消除的方法和装置 Download PDF

Info

Publication number
CN101917527A
CN101917527A CN2010102694004A CN201010269400A CN101917527A CN 101917527 A CN101917527 A CN 101917527A CN 2010102694004 A CN2010102694004 A CN 2010102694004A CN 201010269400 A CN201010269400 A CN 201010269400A CN 101917527 A CN101917527 A CN 101917527A
Authority
CN
China
Prior art keywords
filtering processing
residual signals
energy
input signal
remote end
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.)
Granted
Application number
CN2010102694004A
Other languages
English (en)
Other versions
CN101917527B (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.)
New H3C Technologies Co Ltd
Original Assignee
Hangzhou H3C Technologies Co Ltd
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 Hangzhou H3C Technologies Co Ltd filed Critical Hangzhou H3C Technologies Co Ltd
Priority to CN 201010269400 priority Critical patent/CN101917527B/zh
Publication of CN101917527A publication Critical patent/CN101917527A/zh
Application granted granted Critical
Publication of CN101917527B publication Critical patent/CN101917527B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Telephone Function (AREA)

Abstract

本发明公开了一种回声消除的方法和装置,该方法包括:自适应滤波器根据上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号,确定所述远端输入信号和残差信号的能量;所述自适应滤波器根据上一次滤波处理使用的滤波器系数、所述远端输入信号和残差信号,以及所述远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数;所述自适应滤波器根据确定出的本次滤波器系数执行本次滤波处理。本发明在不进行双端发声检测的情况下实现语音通话中的回声消除。

Description

回声消除的方法和装置
技术领域
本发明涉及通信领域,尤其涉及一种回声消除的方法和装置。
背景技术
在语音通信系统中,远端输入信号到达本地信号接收设备(例如电话机)后,经过本地信号接收设备的音箱、房间等器件到达听筒,在此过程中经常会由于音箱、房间等器件内的声音反射产生回声。为了消除回声,现有技术中在语音信号输出之前使用自适应滤波器进行滤波,目前广泛使用的滤波算法为自适应滤波归一化最小均方算法(NLMS)。
基于NLMS算法的回声消除过程如图1所示,其中x表示远端输入信号,y表示x经过房间形成的实际回声信号,v表示本地说话者的声音及背景噪声,d表示自适应滤波器的近端输入信号,
Figure BSA00000252357400011
表示经过自适应滤波器运算得到的估计回声,e表示滤波输出的残差信号,h表示实际房间脉冲响应,表示自适应滤波器系数,即h的估计值。
图1所示的回声消除过程分为两部分,即
滤波过程:
e ( n ) = d ( n ) - w ^ T ( n ) x ( n ) , …………公式(1)
滤波器系数更新过程:
w ^ ( n + 1 ) = w ^ ( n ) + &mu; | | x ( n ) | | 2 e ( n ) x ( n ) , 0 < &mu; < 2 …………公式(2)
(上式中用斜体表示的是标量(即u、e(n)),粗体表示的是矢量(即
Figure BSA00000252357400015
Figure BSA00000252357400016
)
其中:
w ^ ( n ) = [ w ^ 0 ( n ) , w ^ 1 ( n ) , L , w K - 1 ^ ( n ) ] T
x(n)=[x0(n),x1(n),L,xK-1(n)]T=[x(n),x(n-1),L,x(n-K+1)]T
其中,K为滤波器阶数。
当本地没有说话声音及背景噪声v时,
Figure BSA00000252357400021
取值正常,此时滤波器可以正常工作,但是当滤波器的本地有说话声音及背景噪声v时,自适应滤波器会受到v的干扰,滤波过程中d的量包括了v,导致e的取值变大,滤波器系数
Figure BSA00000252357400022
产生较大的误差,甚至引起
Figure BSA00000252357400023
的发散。这就是回声消除中双端发声问题的原因,双端发声问题已成为回声消除的技术瓶颈。
现有技术中为了克服双端发声问题,采用了双端发声检测技术,在检测到双端发声时,根据双端发声检测的结果调整滤波器系数更新过程中的步长μ,即检测到双端发声后将迭代步长μ直接赋0,停止自适应滤波器的迭代更新。现有技术中双端发声检测技术主要包括以下两种:
(1)能量比较法
能量比较法主要是通过比较回声消除后的残差能量和消除前的近端信号能量来进行双端发声检测。该方法中定义了回音往返损耗增强(ERLE):
ERLE = 10 log 10 ( &Sigma; n d 2 ( n ) &Sigma; n e 2 ( n ) )
当ERLE满足式下式时认为检测到双端发声:
ERLE≤ERLEopt
式中的ERLEopt为预先设定的阈值。该方法的物理意义为:当v不存在时,自适应滤波器应具有较好的稳态收敛性能,所以ERLE应该超过这个阈值;当v存在时,自适应滤波器不可能消除近端信号v,因此输出的残差信号e也包含了v(其中d=y+v),ERLE应该小于这个阈值。
(2)相关计算法
相关计算法利用近端和残差信号的相关性来进行双端发声检测。当不存在双端发声时,由于自适应滤波器能够很好地逼近回声信道,d(n)和
Figure BSA00000252357400025
相关性较大,d(n)和e(n)相关性较小;而当下列两式的条件同时满足时,则认为检测到了双端发声,(式中T1、T2为设定的阈值)
E [ d ( n ) y ^ ( n ) ] &sigma; d 2 &sigma; y ^ 2 &le; T 1
E [ d ( n ) e ( n ) ] &sigma; d 2 &sigma; e 2 &GreaterEqual; T 2
其中d(n)和
Figure BSA00000252357400033
e(n)能量的计算由下式给出:
&sigma; d 2 = ( 1 - &beta; ) &sigma; d 2 + &beta; | d ( n ) | 2
&sigma; y ^ 2 = ( 1 - &beta; ) &sigma; y ^ 2 + &beta; | y ^ ( n ) | 2
&sigma; e 2 = ( 1 - &beta; ) &sigma; e 2 + &beta; | e ( n ) | 2
但是,上述能量法和相关计算法均要求自适应滤波器已经处于收敛或接近收敛状态,而在回声消除的实际应用中往往很难满足这一条件,如在一次通信中对话开始时,产生双端发声的概率比通话过程中要高,此时,自适应滤波器还在收敛过程中,所以利用能量法和相关计算法进行双端发声检测,并控制自适应滤波器停止迭代,将导致回声消除性能恶化和检测错误的恶性循环。同时,双端发声检测所存在的误检将导致自适应滤波器系数不能有效更新,同样也可能造成回声消除性能恶化和检测错误的恶性循环。因此,提供一种既可以消除回声又不需要进行双端发声检测的方案成为亟待解决的问题。
发明内容
本发明实施例提供了一种回声消除的方法和装置,在不进行双端发声检测的情况下实现语音通话中的回声消除。
本发明实施例提供了一种回声消除的方法,包括:
自适应滤波器根据上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号,确定所述远端输入信号和残差信号的能量;
所述自适应滤波器根据上一次滤波处理使用的滤波器系数、所述远端输入信号和残差信号,以及所述远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数;
所述自适应滤波器根据确定出的本次滤波器系数执行本次滤波处理。
所述自适应滤波器通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + 1 &sigma; x 2 + &sigma; e 2 e ( n ) x ( n )
其中,
Figure BSA00000252357400042
表示本次滤波处理使用的滤波器系数,
Figure BSA00000252357400043
表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量。
所述自适应滤波器通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + u N ( &sigma; x 2 + &sigma; e 2 ) e ( n ) x ( n )
其中,
Figure BSA00000252357400045
表示本次滤波处理使用的滤波器系数,
Figure BSA00000252357400046
表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量,N为所述σX 2e 2的系数。
所述自适应滤波器根据所述远端输入信号和残差信号的权重确定所述远端输入信号和残差信号的能量和。
在确定所述远端输入信号和残差信号的能量和之前,还将所述残差信号的能量乘以参数γ,所述参数γ具体为:
Figure BSA00000252357400047
其中,d(n)为当前近端信号,σd 2表示当前近端信号的能量,表示所述滤波得到的残差信号的能量。
本发明实施例提供了一种回声消除的装置,包括:
第一确定单元,用于根据上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号,确定所述远端输入信号和残差信号的能量;
第二确定单元,与所述第一确定单元连接,用于根据上一次滤波处理使用的滤波器系数、所述远端输入信号和残差信号、以及所述第一确定单元确定的所述远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数;
滤波单元,与所述第二确定单元连接,用于根据所述第二确定单元确定出的滤波器系数执行滤波处理。
所述第二确定单元具体用于:
通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + 1 &sigma; x 2 + &sigma; e 2 e ( n ) x ( n )
其中,表示本次滤波处理使用的滤波器系数,
Figure BSA00000252357400053
表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量。
所述第二确定单元具体用于:
通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + u N ( &sigma; x 2 + &sigma; e 2 ) e ( n ) x ( n )
其中,
Figure BSA00000252357400055
表示本次滤波处理使用的滤波器系数,
Figure BSA00000252357400056
表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量,N为所述σX 2e 2的系数。
所述第二确定单元还用于:根据所述远端输入信号和残差信号的权重确定所述远端输入信号和残差信号的能量和。
所述第二确定单元还用于:在确定所述远端输入信号和残差信号的能量和之前,将所述残差信号的能量乘以参数γ,所述参数γ具体为:
其中,d(n)为当前近端信号,σd 2表示当前近端信号的能量,表示所述滤波得到的残差信号的能量。
与现有技术相比,本发明实施例至少具有以下优点:
本发明实施例中,自适应滤波器根据上一次滤波处理使用的滤波器系数、上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号、以及该远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数,当存在近端信号时,由于残差信号的能量增大,自适应滤波器可以通过残差信号的能量调节滤波器系数的变化,当没有近端信号时,由于残差信号的能量很小,可以忽略对滤波器系数变化的影响,保证自适应滤波器的正常工作。
附图说明
图1是现有技术中回声消除过程示意图;
图2是本发明实施例提供的回声消除的方法的流程示意图;
图3是本发明实施例提供的回声消除的装置的结构示意图。
具体实施方式
本发明实施例中,对现有的滤波器系数更新过程进行调整,自适应滤波器根据上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号、该远端输入信号和残差信号的能量和、以及上一次滤波处理使用的滤波器系数确定本次滤波处理使用的滤波器系数,使得滤波器系数的变化与近端信号的能量变化成反比,起到在近端信号出现时减缓滤波器系数迭代更新的作用。
本发明实施例中,对现有技术中的公式(1)和(2)进行调整,根据现有技术中信号能量的表达方式,远端信号x(n)的能量可以表达为:
Figure BSA00000252357400061
其中N为信号能量系数,
Figure BSA00000252357400062
则公式(2)可以调整为如下公式(3)
w ^ ( n + 1 ) = w ^ ( n ) + u N ( &sigma; x 2 + c&sigma; e 2 ) e ( n ) x ( n ) ……公式(3)
公式(3)中,u的取值优选地为:0<μ<2,参数c用来调节σ2 e与σX 2的权重,其取值一般由实验确定,优选地取2~4。
需要说明的是,系数N、u以及c的取值可以根据实际需要灵活调整,例如取1,也可以根据实际需要去除公式(3)中的上述一种或多种系数。当系数N、u以及c都取1,或者说系数N、u以及c都从公式(3)去除后,得到下述公式(4):
w ^ ( n + 1 ) = w ^ ( n ) + 1 &sigma; x 2 + &sigma; e 2 e ( n ) x ( n ) ……公式(4)
结合上述公式,图2示出了本发明实施例提供的回声消除流程,该流程可包括以下步骤:
步骤201,自适应滤波器获取上一次滤波处理时记录的远端输入信号x(n)、滤波得到的残差信号e(n)和上一次滤波器系数
Figure BSA00000252357400071
具体的,自适应滤波器内缓存上一次或几次滤波处理时的远端输入信号、滤波得到的残差信号和上一次滤波器系数,优选的,缓存上一次滤波处理时的远端输入信号、滤波得到的残差信号和上一次滤波器系数即可。当缓存多次远端输入信号、残差信号和滤波器系数时,可以设置标识以表明各次缓存数据的次序关系,使得自适应滤波器可以识别出上一次滤波处理时的远端输入信号、滤波得到的残差信号和上一次滤波器系数。
步骤202,自适应滤波器根据记录的远端输入信号x(n)和残差信号e(n)确定远端输入信号x(n)和残差信号e(n)的能量。
本发明是实例中,自适应滤波器采用现有的能量确定方式确定信号的能量,具体的,x(n)的能量
Figure BSA00000252357400072
e(n)的能量
Figure BSA00000252357400073
步骤203,自适应滤波器根据滤波器系数计算公式(公式(3)或(4))计算本次滤波器系数
Figure BSA00000252357400074
步骤204,自适应滤波器根据本次滤波器系数执行滤波处理。
自适应滤波器确定出
Figure BSA00000252357400076
后,即可根据下述公式进行本次滤波处理:
e ( n ) = d ( n ) - w ^ T ( n ) x ( n ) .
滤波处理结束后,自适应滤波器可记录滤波处理时的远端输入信号x(n+1)、滤波得到的残差信号e(n+1)和该次滤波器系数
Figure BSA00000252357400078
以作为计算下一次滤波处理过程所使用的滤波器系数的依据。
采用本发明实施例提供的方式更新滤波器系数时,在单端说话时(即只有远端输入信号x、不存在本地说话者的声音及背景噪声v时),由于σ2 e的取值很小,因此采用公式(3)或(4)进行计算的结果与采用现有技术中的公式(2)进行计算的结果近似相等,滤波器系数的更新过程与经典算法基本保持一致。在双端说话时(即同时有远端输入信号x和本地说话者的声音及背景噪声v),误差信号中包含了近端信号,且σ2 e与σ2x相当,因此分母会因为误差信号能量的增大而增大,就相当于减小步长μ值,减缓滤波器系数的更新,且减缓程度与误差信号能量成正比,因此可以在不同的状态下进行平滑的过渡。而且,采用本发明实施例提供的方法,在回声路径发生变化时也可以保证滤波器的正常工作:在回声路径刚发生变化的时刻,误差信号也比较大,这时滤波器更新速度相对较慢,但随着滤波器逐渐收敛,误差信号不断减少,从而滤波器更新速度也逐渐加快,如此形成正反馈,滤波器最终的收敛速度没有太大的影响。
需要说明的是,本发明实施例中提供的上述公式(3)或公式(4)还可以根据本发明实施例提供的方法的思想进行其他调整,例如在公式(3)中cσ2 e之前再乘以一个参数γ,将公式(3)调整为:
w ^ ( n + 1 ) = w ^ ( n ) + &mu; N ( &sigma; x 2 + &gamma;c&sigma; e 2 ) e ( n ) x ( n ) …………公式(5)
其中, &gamma; = E [ d ( n ) e ( n ) ] &sigma; d 2 &sigma; e 2
其中,参数γ是近端信号d(n)与误差信号e(n)的互相关系数,在回声消除不同的状态时存在不同的取值范围,其主要作用在于在单端说话时进一步降低误差能量项
Figure BSA00000252357400083
对滤波器正常更新的影响。具体的,在只有远端输入信号x、不存在本地说话者的声音及背景噪声v时,近端信号d(n)与误差信号e(n)相关性很低,参数γ取值趋向于0,进一步减少滤波器系数更新过程中误差能量项对滤波器正常更新的影响。同时,在同时有远端输入信号x和本地说话者的声音及背景噪声v时,近端信号与误差信号都包含近端信号,相关度比较高,因此参数γ取值趋向于1,进而保证滤波器系数更新表达式分母项中误差能量项
Figure BSA00000252357400085
的作用,减缓滤波器系数的更新。在回声路径发生变化时,近端信号d(n)只含有回声信号,残差信号e(n)包含因为回声路径变化没有正常消的回声信号,近端信号与残差信号存在一定的相关度,但可能存在较大的波动范围,因此参数γ取值一般在0.5±0.2,部分减小误差能量项
Figure BSA00000252357400086
对滤波器更新过程的影响;并且随着滤波器逐渐收敛,参数γ取值逐渐减小并趋向于0,进一步减小
Figure BSA00000252357400091
对滤波器更新过程的影响;同时误差信号不断减少,从而滤波器更新速度也逐渐加快,如此形成的正反馈能很快的使滤波器系数恢复到正常的更新过程,即没有影响误差能量项
Figure BSA00000252357400092
对滤波器更新过程的作用。
本发明实施例中,自适应滤波器根据上一次滤波处理使用的滤波器系数、上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号、以及该远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数,当存在近端信号时,由于残差信号的能量增大,自适应滤波器可以通过残差信号的能量调节滤波器系数的变化,当没有近端信号时,由于残差信号的能量很小,可以忽略对滤波器系数变化的影响,保证自适应滤波器的正常工作。
本发明实施例提供了一种回声消除的装置,如图3所示,包括:
第一确定单元11,用于根据上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号,确定所述远端输入信号和残差信号的能量;
第二确定单元12,与所述第一确定单元11连接,用于根据上一次滤波处理使用的滤波器系数、所述远端输入信号和残差信号、以及所述第一确定单元确定的所述远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数;
滤波单元13,与所述第二确定单元12连接,用于根据所述第二确定单元12确定出的滤波器系数执行滤波处理。
所述第二确定单元12具体用于:通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + 1 &sigma; x 2 + &sigma; e 2 e ( n ) x ( n )
其中,表示本次滤波处理使用的滤波器系数,
Figure BSA00000252357400095
表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量。
所述第二确定单元12还可以用于:通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + u N ( &sigma; x 2 + &sigma; e 2 ) e ( n ) x ( n )
其中,
Figure BSA00000252357400102
表示本次滤波处理使用的滤波器系数,
Figure BSA00000252357400103
表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量,N为所述σX 2e 2的系数。
所述第二确定单元12使用上述公式(1)或(2)确定滤波器系数时,还使根据所述远端输入信号和残差信号的权重确定所述远端输入信号和残差信号的能量和,即使用所述σ2 e乘以调节所述σ2 e与σX 2的权重的参数c。
所述第二确定单元12使用上述公式(1)或(2)确定滤波器系数时,还用于:在确定所述远端输入信号和残差信号的能量和之前,将所述残差信号的能量乘以参数γ,所述参数γ具体为:
其中,d(n)为当前近端信号,σd 2表示当前近端信号的能量,
Figure BSA00000252357400105
表示所述滤波得到的残差信号的能量。
本发明实施例中,自适应滤波器根据上一次滤波处理使用的滤波器系数、上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号、以及该远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数,当存在近端信号时,由于残差信号的能量增大,自适应滤波器可以通过残差信号的能量调节滤波器系数的变化,当没有近端信号时,由于残差信号的能量很小,可以忽略对滤波器系数变化的影响,保证自适应滤波器的正常工作。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
本领域技术人员可以理解附图只是一个优选实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
本领域技术人员可以理解实施例中的装置中的模块可以按照实施例描述进行分布于实施例的装置中,也可以进行相应变化位于不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块,也可以进一步拆分成多个子模块。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上公开的仅为本发明的几个具体实施例,但是,本发明并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。

Claims (10)

1.一种回声消除的方法,其特征在于,包括:
自适应滤波器根据上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号,确定所述远端输入信号和残差信号的能量;
所述自适应滤波器根据上一次滤波处理使用的滤波器系数、所述远端输入信号和残差信号,以及所述远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数;
所述自适应滤波器根据确定出的本次滤波器系数执行本次滤波处理。
2.如权利要求1所述的方法,其特征在于,所述自适应滤波器通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + 1 &sigma; x 2 + &sigma; e 2 e ( n ) x ( n )
其中,表示本次滤波处理使用的滤波器系数,表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量。
3.如权利要求1所述的方法,其特征在于,所述自适应滤波器通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + u N ( &sigma; x 2 + &sigma; e 2 ) e ( n ) x ( n )
其中,
Figure FSA00000252357300015
表示本次滤波处理使用的滤波器系数,表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量,N为所述σX 2e 2的系数。
4.如权利要求3所述的方法,其特征在于,所述自适应滤波器根据所述远端输入信号和残差信号的权重确定所述远端输入信号和残差信号的能量和。
5.如权利要求2-4任一项所述的方法,其特征在于,在确定所述远端输入信号和残差信号的能量和之前,还将所述残差信号的能量乘以参数γ,所述参数γ为:
&gamma; = E [ d ( n ) e ( n ) ] &sigma; d 2 &sigma; e 2
其中,d(n)为当前近端信号,σd 2表示当前近端信号的能量,
Figure FSA00000252357300022
表示所述滤波得到的残差信号的能量。
6.一种回声消除的装置,其特征在于,包括:
第一确定单元,用于根据上一次滤波处理时的远端输入信号、上一次滤波处理得到的残差信号,确定所述远端输入信号和残差信号的能量;
第二确定单元,与所述第一确定单元连接,用于根据上一次滤波处理使用的滤波器系数、所述远端输入信号和残差信号、以及所述第一确定单元确定的所述远端输入信号和残差信号的能量和,确定本次滤波处理使用的滤波器系数;
滤波单元,与所述第二确定单元连接,用于根据所述第二确定单元确定出的滤波器系数执行滤波处理。
7.如权利要求6所述的装置,其特征在于,所述第二确定单元具体用于:
通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + 1 &sigma; x 2 + &sigma; e 2 e ( n ) x ( n )
其中,
Figure FSA00000252357300024
表示本次滤波处理使用的滤波器系数,
Figure FSA00000252357300025
表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量。
8.如权利要求6所述的装置,其特征在于,所述第二确定单元具体用于:
通过下述公式确定本次滤波处理使用的滤波器系数,
w ^ ( n + 1 ) = w ^ ( n ) + u N ( &sigma; x 2 + &sigma; e 2 ) e ( n ) x ( n )
其中,
Figure FSA00000252357300027
表示本次滤波处理使用的滤波器系数,
Figure FSA00000252357300028
表示上一次滤波处理使用的滤波器系数,e(n)表示上一次滤波处理得到的残差信号,x(n)表示上一次滤波处理时的远端输入信号,σX 2表示所述远端输入信号的能量,σe 2表示所述残差信号的能量,N为所述σX 2e 2的系数。
9.如权利要求8所述的装置,其特征在于,所述第二确定单元还用于:
根据所述远端输入信号和残差信号的权重确定所述远端输入信号和残差信号的能量和。
10.如权利要求6-8任一项所述的装置,其特征在于,所述第二确定单元还用于:在确定所述远端输入信号和残差信号的能量和之前,将所述残差信号的能量乘以参数γ,所述参数γ具体为:
Figure FSA00000252357300031
其中,d(n)为当前近端信号,σd 2表示当前近端信号的能量,
Figure FSA00000252357300032
表示所述滤波得到的残差信号的能量。
CN 201010269400 2010-09-02 2010-09-02 回声消除的方法和装置 Active CN101917527B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010269400 CN101917527B (zh) 2010-09-02 2010-09-02 回声消除的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010269400 CN101917527B (zh) 2010-09-02 2010-09-02 回声消除的方法和装置

Publications (2)

Publication Number Publication Date
CN101917527A true CN101917527A (zh) 2010-12-15
CN101917527B CN101917527B (zh) 2013-07-03

Family

ID=43324905

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010269400 Active CN101917527B (zh) 2010-09-02 2010-09-02 回声消除的方法和装置

Country Status (1)

Country Link
CN (1) CN101917527B (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102185991A (zh) * 2011-03-01 2011-09-14 杭州华三通信技术有限公司 回声消除方法、系统和装置
CN104036784A (zh) * 2014-06-06 2014-09-10 华为技术有限公司 一种回声消除方法及装置
CN106157965A (zh) * 2016-05-12 2016-11-23 西南交通大学 一种基于权向量重用的零范数集员仿射投影自适应回声消除方法
CN106796779A (zh) * 2014-06-13 2017-05-31 美国思睿逻辑有限公司 用于选择性地启用和停用自适应消噪系统的调整的系统及方法
CN106898359A (zh) * 2017-03-24 2017-06-27 上海智臻智能网络科技股份有限公司 音频信号处理方法、系统、音频交互装置及计算机设备
CN107888792A (zh) * 2017-10-19 2018-04-06 浙江大华技术股份有限公司 一种回声消除方法、装置及系统
CN108550371A (zh) * 2018-03-30 2018-09-18 北京云知声信息技术有限公司 智能语音交互设备快速稳定的回声消除方法
CN108831497A (zh) * 2018-05-22 2018-11-16 出门问问信息科技有限公司 一种回声压缩方法及装置、存储介质、电子设备
CN108962276A (zh) * 2018-07-24 2018-12-07 北京三听科技有限公司 一种语音分离方法及装置
CN109559756A (zh) * 2018-10-26 2019-04-02 北京佳讯飞鸿电气股份有限公司 滤波系数确定方法、回声消除方法、相应装置及设备
CN109785853A (zh) * 2019-03-11 2019-05-21 出门问问信息科技有限公司 一种回声消除方法、装置、系统及存储介质
CN109961798A (zh) * 2017-12-26 2019-07-02 华平信息技术股份有限公司 回声消除系统、方法、可读计算机存储介质、及终端
CN110021289A (zh) * 2019-03-28 2019-07-16 腾讯科技(深圳)有限公司 一种声音信号处理方法、装置及存储介质
CN110211599A (zh) * 2019-06-03 2019-09-06 Oppo广东移动通信有限公司 应用唤醒方法、装置、存储介质及电子设备
CN110326307A (zh) * 2017-03-16 2019-10-11 松下知识产权经营株式会社 声学回声抑制装置以及声学回声抑制方法
CN110995951A (zh) * 2019-12-13 2020-04-10 展讯通信(上海)有限公司 基于双端发声检测的回声消除方法、装置及系统
CN111355855A (zh) * 2020-03-12 2020-06-30 紫光展锐(重庆)科技有限公司 回声处理方法、装置、设备及存储介质
CN112492112A (zh) * 2020-11-19 2021-03-12 睿云联(厦门)网络通讯技术有限公司 一种基于对讲系统的回音消除方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000016497A1 (en) * 1998-09-17 2000-03-23 Ericsson Inc. Echo canceler adaptive filter optimization
CN1937432A (zh) * 2006-09-30 2007-03-28 南京大学 基于最优参数预测的声回声抵消处理方法
CN101043560A (zh) * 2006-03-22 2007-09-26 北京大学深圳研究生院 回声消除器及回声消除方法
CN101179294A (zh) * 2006-11-09 2008-05-14 爱普拉斯通信技术(北京)有限公司 自适应回声消除器及其回声消除方法
CN101179296A (zh) * 2006-11-10 2008-05-14 索尼株式会社 回音消除器和话音处理装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000016497A1 (en) * 1998-09-17 2000-03-23 Ericsson Inc. Echo canceler adaptive filter optimization
CN101043560A (zh) * 2006-03-22 2007-09-26 北京大学深圳研究生院 回声消除器及回声消除方法
CN1937432A (zh) * 2006-09-30 2007-03-28 南京大学 基于最优参数预测的声回声抵消处理方法
CN101179294A (zh) * 2006-11-09 2008-05-14 爱普拉斯通信技术(北京)有限公司 自适应回声消除器及其回声消除方法
CN101179296A (zh) * 2006-11-10 2008-05-14 索尼株式会社 回音消除器和话音处理装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JULIE E. GREENBERG ET AL: "Modified LMS algorithms for speech processing with an adaptive noise canceller", 《IEEE TRANSACTIONS ON SPEECH AND AUDIO PROCESSING》 *
JULIE E. GREENBERG ET AL: "Modified LMS algorithms for speech processing with an adaptive noise canceller", 《IEEE TRANSACTIONS ON SPEECH AND AUDIO PROCESSING》, vol. 6, no. 4, 31 July 1998 (1998-07-31), pages 338 - 351, XP011054314 *

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102185991A (zh) * 2011-03-01 2011-09-14 杭州华三通信技术有限公司 回声消除方法、系统和装置
CN104036784A (zh) * 2014-06-06 2014-09-10 华为技术有限公司 一种回声消除方法及装置
CN106796779A (zh) * 2014-06-13 2017-05-31 美国思睿逻辑有限公司 用于选择性地启用和停用自适应消噪系统的调整的系统及方法
CN106796779B (zh) * 2014-06-13 2020-12-22 美国思睿逻辑有限公司 用于选择性地启用和停用自适应消噪系统的调整的系统及方法
CN106157965A (zh) * 2016-05-12 2016-11-23 西南交通大学 一种基于权向量重用的零范数集员仿射投影自适应回声消除方法
CN106157965B (zh) * 2016-05-12 2019-05-17 西南交通大学 一种基于权向量重用的零范数集员仿射投影自适应回声消除方法
CN110326307A (zh) * 2017-03-16 2019-10-11 松下知识产权经营株式会社 声学回声抑制装置以及声学回声抑制方法
CN106898359A (zh) * 2017-03-24 2017-06-27 上海智臻智能网络科技股份有限公司 音频信号处理方法、系统、音频交互装置及计算机设备
CN106898359B (zh) * 2017-03-24 2020-03-17 上海智臻智能网络科技股份有限公司 音频信号处理方法、系统、音频交互装置及计算机设备
US11151976B2 (en) 2017-10-19 2021-10-19 Zhejiang Dahua Technology Co., Ltd. Methods and systems for operating a signal filter device
CN107888792A (zh) * 2017-10-19 2018-04-06 浙江大华技术股份有限公司 一种回声消除方法、装置及系统
CN107888792B (zh) * 2017-10-19 2019-09-17 浙江大华技术股份有限公司 一种回声消除方法、装置及系统
CN109961798B (zh) * 2017-12-26 2021-06-11 华平信息技术股份有限公司 回声消除系统、方法、可读计算机存储介质、及终端
CN109961798A (zh) * 2017-12-26 2019-07-02 华平信息技术股份有限公司 回声消除系统、方法、可读计算机存储介质、及终端
CN108550371A (zh) * 2018-03-30 2018-09-18 北京云知声信息技术有限公司 智能语音交互设备快速稳定的回声消除方法
CN108550371B (zh) * 2018-03-30 2021-06-01 云知声智能科技股份有限公司 智能语音交互设备快速稳定的回声消除方法
CN108831497A (zh) * 2018-05-22 2018-11-16 出门问问信息科技有限公司 一种回声压缩方法及装置、存储介质、电子设备
CN108831497B (zh) * 2018-05-22 2020-06-09 出门问问信息科技有限公司 一种回声压缩方法及装置、存储介质、电子设备
CN108962276A (zh) * 2018-07-24 2018-12-07 北京三听科技有限公司 一种语音分离方法及装置
CN108962276B (zh) * 2018-07-24 2020-11-17 杭州听测科技有限公司 一种语音分离方法及装置
CN109559756B (zh) * 2018-10-26 2021-05-14 北京佳讯飞鸿电气股份有限公司 滤波系数确定方法、回声消除方法、相应装置及设备
CN109559756A (zh) * 2018-10-26 2019-04-02 北京佳讯飞鸿电气股份有限公司 滤波系数确定方法、回声消除方法、相应装置及设备
CN109785853B (zh) * 2019-03-11 2020-06-16 出门问问信息科技有限公司 一种回声消除方法、装置、系统及存储介质
CN109785853A (zh) * 2019-03-11 2019-05-21 出门问问信息科技有限公司 一种回声消除方法、装置、系统及存储介质
CN110021289A (zh) * 2019-03-28 2019-07-16 腾讯科技(深圳)有限公司 一种声音信号处理方法、装置及存储介质
CN110021289B (zh) * 2019-03-28 2021-08-31 腾讯科技(深圳)有限公司 一种声音信号处理方法、装置及存储介质
CN110211599A (zh) * 2019-06-03 2019-09-06 Oppo广东移动通信有限公司 应用唤醒方法、装置、存储介质及电子设备
CN110211599B (zh) * 2019-06-03 2021-07-16 Oppo广东移动通信有限公司 应用唤醒方法、装置、存储介质及电子设备
CN110995951A (zh) * 2019-12-13 2020-04-10 展讯通信(上海)有限公司 基于双端发声检测的回声消除方法、装置及系统
CN111355855A (zh) * 2020-03-12 2020-06-30 紫光展锐(重庆)科技有限公司 回声处理方法、装置、设备及存储介质
CN111355855B (zh) * 2020-03-12 2021-06-15 紫光展锐(重庆)科技有限公司 回声处理方法、装置、设备及存储介质
CN112492112A (zh) * 2020-11-19 2021-03-12 睿云联(厦门)网络通讯技术有限公司 一种基于对讲系统的回音消除方法和装置
CN112492112B (zh) * 2020-11-19 2022-03-18 睿云联(厦门)网络通讯技术有限公司 一种基于对讲系统的回音消除方法和装置

Also Published As

Publication number Publication date
CN101917527B (zh) 2013-07-03

Similar Documents

Publication Publication Date Title
CN101917527B (zh) 回声消除的方法和装置
CN102065190B (zh) 一种回声消除方法及其装置
CN101461257B (zh) 一种自适应声学回声消除方法和系统
US9653092B2 (en) Method for controlling acoustic echo cancellation and audio processing apparatus
CN103718538B (zh) 可实现声学回声消除的音频信号非线性后期处理方法和系统
EP2804177A2 (en) Method for processing an audio signal and audio receiving circuit
CN109686381A (zh) 用于信号增强的信号处理器和相关方法
US11450335B2 (en) Method and device for updating coefficient vector of finite impulse response filter
US10880440B2 (en) Echo canceller and method therefor
US11349525B2 (en) Double talk detection method, double talk detection apparatus and echo cancellation system
CN105391879A (zh) 一种无回声残留双端通话鲁棒的声学回声消除方法
CN107134281A (zh) 一种自适应回声消除中自适应滤波器系数更新方法
CN104270539B (zh) 基于系数差的成比例仿射投影回声消除方法
SE533956C2 (sv) Enhet och metod för styrning av restekodämpning
CN110992923B (zh) 回声消除方法、电子设备以及存储装置
CN106533500A (zh) 一种优化回声消除器收敛特性的方法
CN105721729B (zh) 基于块稀疏成比例重用权系数仿射投影的回声消除方法
CN109040497B (zh) 一种基于m估计的成比例类仿射投影自适应回声消除方法
US8345860B1 (en) Method and system for detection of onset of near-end signal in an echo cancellation system
CN110191245B (zh) 一种基于时变参数的自适应回声消除方法
CN105491256A (zh) 一种声学回声消除器启动阶段稳健的步长调整方法
CN107071196B (zh) 一种自适应回声消除方法
CN110148421B (zh) 一种残余回声检测方法、终端和装置
CN106033673A (zh) 一种近端语音信号检测方法及装置
CN102655558B (zh) 一种双端发音鲁棒结构及其消除声学回声的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 310052 Binjiang District Changhe Road, Zhejiang, China, No. 466, No.

Patentee after: Xinhua three Technology Co., Ltd.

Address before: 310053 Hangzhou hi tech Industrial Development Zone, Zhejiang province science and Technology Industrial Park, No. 310 and No. six road, HUAWEI, Hangzhou production base

Patentee before: Huasan Communication Technology Co., Ltd.