CN1503466B - 控制稀疏通道中回波消除的方法、电路和系统 - Google Patents
控制稀疏通道中回波消除的方法、电路和系统 Download PDFInfo
- Publication number
- CN1503466B CN1503466B CN2003101183656A CN200310118365A CN1503466B CN 1503466 B CN1503466 B CN 1503466B CN 2003101183656 A CN2003101183656 A CN 2003101183656A CN 200310118365 A CN200310118365 A CN 200310118365A CN 1503466 B CN1503466 B CN 1503466B
- Authority
- CN
- China
- Prior art keywords
- shade
- sparse
- weight
- subclass
- response
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B3/00—Line transmission systems
- H04B3/02—Details
- H04B3/20—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other
- H04B3/23—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other using a replica of transmitted signal in the time domain, e.g. echo cancellers
- H04B3/234—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other using a replica of transmitted signal in the time domain, e.g. echo cancellers using double talk detection
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B3/00—Line transmission systems
- H04B3/02—Details
- H04B3/20—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other
- H04B3/23—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other using a replica of transmitted signal in the time domain, e.g. echo cancellers
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Telephone Function (AREA)
Abstract
稀疏通道回波消除器中的控制逻辑,包括:使能发生器,响应于检测到回波通道中双通话状况和脉冲响应变化的至少一个确立稀疏更新使能,回波通道接收远端输入和近端输入,使能发生器包括:双通话检测器,检测双通话状况;变化检测器,检测脉冲响应变化;主权重更新器,根据该稀疏更新使能而将阴影滤波器的阴影权重传递给主滤波器的多个稀疏滤波器中的多个集合的稀疏权重,主滤波器过滤远端输入产生主输出,阴影滤波器过滤远端输入产生阴影输出,估算稀疏权重的峰值的延迟,该主权重更新器包括:峰值检测器,检测阴影权重中多个峰值的多个峰值位置;稀疏权重更新器,根据多个峰值位置基于阴影权重替换多个集合的稀疏权重。
Description
技术领域
本发明的实施例涉及通信领域,更具体而言涉及回波消除。
技术背景
在语音处理设备如电话网络和无线听筒的设计中,声回波是固有的问题。回波消除是一种技术,其通过估算在发送侧的话筒或线路上拾取的回波并将其从接收器侧的耳机或免提扬声器减去来减小声音传输过程中的回波量。由于环境的准确模型是未知且时间变化的,适应技术被典型地用于估算回波。
在回波消除中有两个重要的现实问题:双通话(double-talk)和由于扬声器封装话筒(LEM)变化导致的脉冲响应的变化。双通话是由来自四线通信系统的电子双工机的远端和近端端口两者的语音信号的同时产生而导致的。LEM变化是由各种局部变化如扬声器的移动、环境噪声和温度变化导致的。对于回波消除器,这两个问题需要相反的行为。双通话需要适应权重被保持不变,即没有适应性,而LEM变化需要快的适应性。另外,当通道脉冲响应稀疏时,计算可能是低效的并且可能是不精确的。
发明内容
本发明提供一种用于稀疏通道的回波消除器中的控制逻辑,包括:
使能发生器,其响应于检测到回波通道中的双通话的状况和脉冲响应变化的至少一个而确立稀疏更新使能,该回波通道接收一个远端输入和一个近端输入,所述使能发生器包括:
双通话检测器,用于检测双通话的状况,以及
变化检测器,用于检测脉冲响应变化;以及
主权重更新器,其根据所述稀疏更新使能而将阴影滤波器的阴影权重传递给主滤波器的多个稀疏滤波器中的多个集合的稀疏权重,主滤波器非适应地过滤远端输入以产生主输出,阴影滤波器适应地过滤远端输入以产生阴影输出,并且估算稀疏权重的峰值的延迟,所述主权重更新器包括:
峰值检测器,其检测阴影权重中多个峰值的多个峰值位置;以及
稀疏权重更新器,其根据多个峰值位置基于阴影权重来替换多个集合的稀疏权重。
本发明还提供一种控制用于稀疏通道的回波消除器的方法,包括:
响应于检测到回波通道中的双通话的状况和脉冲响应变化的至少一个而确立稀疏更新使能,该回波通道接收一个远端输入和一个近端输入,确立稀疏更新使能包括:
检测双通话的状况,以及
检测脉冲响应变化;以及
根据所述稀疏更新使能而将阴影滤波器的阴影权重传递给主滤波器的多个稀疏滤波器中的多个集合的稀疏权重,主滤波器非适应地过滤远端输入以产生主输出,阴影滤波器适应地过滤远端输入以产生阴影输出,并且估算稀疏权重的峰值的延迟,所述传递包括:
检测阴影权重中多个峰值的多个峰值位置;以及根据多个峰值位置基于阴影权重来替换多个集合的稀疏权重。
本发明还提供一种用于稀疏通道的回波消除器的系统,包括:
主滤波器,其具有有多个集合的稀疏权重的多个稀疏滤波器,以非适应地过滤回波通道中的远端输入,主滤波器产生主输出;
阴影滤波器,其具有阴影权重,用于在回波通道中过滤远端输入以适应脉冲响应变化,阴影滤波器产生阴影输出并估算稀疏权重中的峰值的延迟;以及
控制逻辑,其被耦合于主滤波器和阴影滤波器,检测双通话和脉冲响应变化之一,该控制逻辑包括:
使能发生器,响应于检测到双通话的状况和脉冲响应变化的至少一个而确立稀疏更新使能,所述使能发生器包括:
双通话检测器,用于检测双通话的状况,以及
变化检测器,用于检测脉冲响应变化;以及
主权重更新器,其被耦合于使能发生器,根据稀疏更新使能将阴影滤波器的阴影权重传递给多个集合的稀疏权重,所述主权重更新器包括:
峰值检测器,其检测阴影权重中多个峰值的多个峰值位置;以及稀疏权重更新器,其根据多个峰值位置基于阴影权重来替换多个集合的稀疏权重。
附图说明
可通过参考被用于说明本发明实施例的以下描述和附图而最好地了解本发明。在附图中:
图1是说明在其中可实践本发明的一个实施例的系统的图。
图2是说明依照本发明一个实施例的回波消除器的图。
图3是说明依照本发明一个实施例在图2中所示的使能发生器的图。
图4是说明依照本发明一个实施例在图2中所示的主权重更新器的图。
图5是说明依照本发明一个实施例执行回波消除的过程的流程图。
图6是说明依照本发明一个实施例更新图5中所示的稀疏权重集合的过程的流程图。
图7是说明依照本发明一个实施例在200K次迭代之后组合主权重的图。
图8是说明依照本发明一个实施例的平滑均方误差(MSE)的图。
图9是说明依照本发明一个实施例的使用延迟系数的双通话的状况的图。
图10是说明依照本发明一个实施例的平滑LEM变化状况的图。
图11是说明依照本发明一个实施例在迭代104800到106800处的平滑LEM变化状况的图。
具体实施方式
本发明是用于回波通道中的稀疏通道的回波消除器.回波通道接收远端输入和近端输入.响应于检测到双通话的状况和脉冲响应变化的至少一个,使能发生器确立(assert)稀疏更新使能.根据该稀疏更新使能,主权重更新器将阴影(shadow)滤波器的阴影权重传递给主滤波器的N个稀疏滤波器中的N个集合的稀疏权重.稀疏主滤波器非适应地过滤远端输入以产生主输出.长阴影滤波器适应地过滤远端输入以产生阴影输出并估算主滤波器中稀疏权重的峰值的延迟.
在以下说明书中,提出了诸多特定细节。然而应理解,本发明的实施例可无需这些特定细节而被实践。在其它实例中,众所周知的电路、结构和技术未被示出是为了不模糊对本说明书的理解。
本发明一个实施例的单元可由硬件、固件、软件或其任何组合来实施。当被实施以软件时,本发明实施例的单元基本上是执行必要任务的代码段。软件可包括实施在本发明一个实施例中描述的操作的实际代码,或者仿真(emulate)或模拟所述操作的代码。程序或代码段可被存储在处理器或机器可访问介质中,或者在传输介质上通过在载波中实施的计算机数据信号或由载波调制的信号来发送。“处理器可读或可访问介质”或“机器可读或可访问介质”可包括可存储、发送或传递信息的任何介质。机器可访问介质的实例包括电子电路、半导体存储器器件、只读存储器(ROM)、闪存、可擦除ROM(EROM)、软盘、高密度盘(CD)ROM、光盘、硬盘、光纤介质、射频(RF)链路等。计算机数据信号可包括可在诸如电子网络通道、光纤、无线电、电磁、RF链路等的传输介质上传播的任何信号。代码段可通过诸如互联网、内部网等的计算机网络来下载。机器可访问介质可被实施于一种产品中。机器可访问介质可包括这样的数据,即当被机器访问时,其使该机器执行在以下描述的操作。术语“数据”在这里指为机器可访问的目的而被编码的任何类型的信息。因此它可包括程序、代码、数据、文件等。
本发明实施例的全部或部分可由软件来实施。软件可具有被相互耦合的几个模块。一个软件模块被耦合于另一个模块以接收变量、参数、自变量、指针等,并且/或者产生或传递结果、更新的变量、指针等。软件模块亦可以是与在平台上运行的操作系统交互的软件驱动程序(driver)或接口。软件模块亦可以是配置、设置、初始化数据、发送数据给硬件装置并从其接收数据的硬件驱动程序。
本发明的一个实施例可被描述为一个过程,其通常被描述为流程图、程序图、结构图或方块图。尽管流程图可将操作描述成按顺序的过程,许多操作可被并行或同时执行。另外,操作的顺序可被重新安排。当其操作被完成时,过程被终止。过程可对应于方法、程序、步骤(procedure)等。
图1是说明在其中可实践本发明的一个实施例的系统100的图。系统100包括发送输入解码器110、回波通道120、发送输出编码器130、接收输入解码器140、接收输出编码器150、网络145、发送输入解码器160、回波通道170、发送输出编码器180、接收输入解码器190和接收输出编码器195。
发送输入解码器110从第一近端接收经编码的语音并将经编码的语音解码为线性语音数据Sin。在一个实施例中,发送输入解码器110是μ定律/A定律解码器。回波通道120包括回波消除器125。回波消除器125从线性数据采样Sin去除回波估算信号以产生线性数据采样Sout.发送输出编码器130在打包之前提供语音压缩.在一个实施例中,发送输出编码器130是G.7xx编码器,其使用用于低位率语音(LBRV)的任何一个压缩标准,包括国际电信联盟(ITU)-T国际标准化的G.7xx系列,以压缩来自回波通道120的语音数据Sout。经压缩的语音数据通过网络发送给远端。接收输入解码器140解压缩在网络145上从第一远端接收的语音数据。解压缩技术与在发送输出编码器130中使用的压缩兼容。回波通道120接收来自接收输入解码器140的Rin并发送出Rout线性数据采样。接收输出编码器150将线性数据采样Rout编码为待发送给第一近端的μ定律和A定律编码的语音。
网络145是具有以下能力的任何网络:从发送输出解码器130、发送输入解码器160、接收输入解码器140和接收输出解码器195发送打包数据并将打包数据发送给它们。网络145可以是互联网、内部网、外部网、局域网(LAN)或广域网(WAN)。发送输入解码器160从网络145接收经编码的语音并将经编码的语音解码为线性语音数据Sin。在一个实施例中,发送输入解码器160是μ定律/A定律解码器。回波通道170包括回波消除器175。回波消除器175从线性数据采样Sin去除回波估算信号以产生线性数据采样Sout。发送输出编码器180在打包之前提供语音压缩。在一个实施例中,发送输出编码器180是G.7xx编码器,其使用用于低位率语音(LBRV)的任何一个压缩标准,包括国际电信联盟(ITU)-T国际标准化的G.7xx系列,以压缩来自回波通道170的语音数据Sout。经压缩的语音数据被发送给第二远端处的接收装置。接收输入解码器190解压缩从第二远端接收的语音数据。解压缩技术与在发送输出编码器180中使用的压缩兼容。回波通道170接收来自接收输入解码器190的Rin并发送出Rout线性数据采样。接收输出编码器190将线性数据采样Rout编码为待发送给对网络145的第二近端的μ定律和A定律编码的语音。在一个实施例中,发送输入解码器160、回波通道170、发送输出解码器180、接收输入解码器190和接收输出编码器195被集成为数字信号处理器165。
在以下描述中,上标T表示矢量或矩阵的转置,符号*表示乘法,||.||表示绝对值。
图2是说明依照本发明一个实施例的回波消除器125的图。回波消除器125包括主滤波器210、减法器220、阴影滤波器230、减法器250和控制逻辑260。
稀疏主滤波器210是短或稀疏滤波器。它非适应地或没有适应性地过滤远端输入Rin以产生主输出ym(k)。主滤波器210包括矢量器212、乘法器214和权重组合器216。矢量器212将远端输入Rin延迟一个抽头延迟线(tapped delay line)以产生矢量u(k)。矢量器212可由缓冲存储器、移位寄存器或索引阵列来实施。乘法器214如下使u(k)乘以权重组合器216提供的组合权重Wc(k)以进行过滤过程:
ym(k)=u(k)*Wc(k)T (1)
权重组合器216如下组合或合并主滤波器210的N个集合的稀疏权重Wm1,Wm2,...,WmN以生成组合的权重Wc(k):
对于1≤k≤ND
否则
其中ND是耦合距离。耦合距离被用于检测双通话并将在以后被讨论。
减法器220如下从主输出减去近端输入Sin以产生主误差em(k):
em(k)=Sin(k)-ym(k) (4)
主误差em(n)由控制逻辑用来产生使能,用于更新稀疏权重Wmi(k),其中i=1,...,N.
阴影滤波器230是长适应滤波器。它适应地或有适应性地过滤远端输入Rin以产生阴影输出ys(k)。在一个实施例中,适应性规则是2阶仿射投影(AP),AP(2)适应性或更新规则。阴影滤波器230包括矢量器232、乘法器234、阴影权重Wsh 240和阴影权重更新器245。矢量器232类似于延迟线212并可被与主滤波器210共享。矢量器232将远端输入Rin延迟一个抽头延迟线以产生矢量u(k)。乘法器230如下使矢量u(k)乘以阴影权重Wsh 240以进行过滤过程:
ys(k)=u(k)*Wsh(k)T (5)
当被控制逻辑260使能时,阴影权重更新器245使用AP(2)适应性规则来更新阴影权重Wsh。
减法器250如下从阴影输出ys(k)减去近端输入Sin以产生阴影误差es(k):
es(k)=Sin(k)-ys(k) (6)
阴影误差es(k)由控制逻辑260用来产生稀疏更新使能。其亦被用于如下使用AP(2)规则来更新阴影权重Wsh:
Wsh(k+1)=Wsh(k)+[φ(k)/(φ(k)*φT(k)]*es(k) (9)
其中U(k)是P个过去矢量(past vector)的集合
U(k)=[u(k-1),u(k-2),...,u(k-P)]。
控制逻辑260执行许多功能以控制主滤波器和阴影滤波器210和230的过滤操作和权重更新。控制逻辑接收远端输入Rin以及主误差和阴影误差em(k)和es(n)。控制逻辑260包括使能发生器262和主权重更新器265。使能发生器262产生阴影更新使能以使能阴影滤波器230中的阴影权重更新器245。它亦产生稀疏更新使能以使能主权重更新器265。响应于检测到双通话的状况和脉冲响应变化的至少一个,如回波通道中的扬声器封装话筒(LEM)状况,它确立稀疏更新使能。当稀疏更新使能被确立时,主权重更新器265使用阴影权重Wsh 240的子集来更新主滤波器210的N个稀疏滤波器中的N个集合的稀疏权重。
图3是说明依照本发明一个实施例在图2中所示的使能发生器262的图。使能发生器262包括阴影更新使能发生器301和稀疏更新使能发生器302。
阴影更新使能发生器301产生阴影更新使能。当被确立时,阴影更新使能使阴影权重更新器245(图2)能更新阴影权重Wsh 240(图2)。它包括短期功率估算器310、长期功率估算器320和适应性使能器330。短期功率估算器310如下估算远端输入Rin的短期功率pS:
x1(k)=||Rin(k-M+1)||2 (10)
若x1(k)>pS(k-1)则pS(k)=(1-γr)x1(k)+γrpS(k-1) (11)
否则pS(k)=(1-γf)x1(k)+γfPS(k-1) (12)
其中:
γr和γf是有关过去采样数的时间常数,分别具有1-(1/24)和1-(1/240)的典型值;M是阴影滤波器的抽头数;M1是每个稀疏滤波器的抽头数。在一个实施例中,M=512+ND=544且M1=128。
长期功率估算器320如下估算远端输入Rin的长期功率pL:
若pS(k)>r1pL(k-1)则pL(k)=(1-α)x1(k)+αpL(k-1) (13)
否则pL(k)=pL(k-1) (14)
其中r1是正整数,具有从4到10的典型值。
适应性使能器330产生阴影更新使能。响应于短期功率pS超过长期功率pL的指示(indication),它确立阴影更新使能以使能阴影滤波器230的适应性。长期功率的指示是长期功率的倍数,例如r1pL,其中r1是正整数。换句话说,当Ps>r1pL时,阴影权重Wsh(k)被更新。在一个实施例中,r1的范围是从4到10。
稀疏更新使能发生器302包括延迟系数发生器340、双通话检测器350、主误差功率估算器360、阴影误差功率估算器370、脉冲响应变化检测器380和组合器390。延迟系数发生器340如下计算用于主滤波器210的主耦合因子βm和用于阴影滤波器230的阴影耦合因子βs:
βm=(M/ND)*[Wm(1:ND)*Wm(1:ND)T] (15)
βs=(M/ND)*[Wsh(1:ND)*Wsh(1:ND)T] (16)
其中Wm(1:ND)和Wsh(1:ND)指稀疏权重和阴影权重的第一ND权重;ND是耦合距离,其有依赖于滤波器尺寸M和M1的值。在一个实施例中,M=544,M1=128,且ND=32。
响应于主耦合因子βm超过阴影耦合因子βs的指示,双通话检测器检测双通话的状况。阴影耦合因子βs的指示是βs的倍数,即r2*βs,其中r2是正整数。在一个实施例中,r2=2。换句话说,状况是βm>2*βs。响应于检测到双通话的状况,双通话检测器350确立双通话的使能。
主误差功率估算器360如下估算主误差em(k)的主误差功率epm:
e1(k)=||em(k)||2 (17)
若e1(k)>epm(k-1)则epm(k)=(1-γr)e1(k)+γrepm(k-1) (18)
否则epm(k)=(1-γf)e1(k)+γfepm(k-1) (19)
其中γr和γf如以前被定义。
阴影误差功率估算器370如下估算阴影误差es(k)的阴影误差功率eps:
e2(k)=||es(k)||2 (20)
若e2(k)>eps(k-1)则eps(k)=(1-γr)e2(k)+γreps(k-1) (21)
否则eps(k)=(1-γf)e2(k)+γfeps(k-1) (22)
响应于主误差功率epm超过阴影误差功率eps的指示,变化检测器380检测脉冲响应变化。阴影误差功率eps的指示是eps的倍数,例如r3eps,其中r3是正整数。在一个实施例中,r3=100。换句话说,当epm>r3eps时,检测到脉冲响应变化。当检测到变化状况时或响应于检测到变化状况,变化检测器380确立变化使能。
组合器390组合如双通话检测器350和LEM检测器380分别确立的所检测的双通话和LEM变化,从而确立稀疏更新使能。如果这些状况的任何一个被确立,主权重更新器被使能以更新主滤波器210中的N个集合的稀疏权重。在一个实施例中,组合器350是逻辑OR运算器。
图4是说明依照本发明一个实施例在图2中所示的主权重更新器265的图。主权重更新器265包括峰值检测器410和稀疏权重更新器420。
峰值检测器410检测阴影权重Wsh 240(图2)中N个峰值的N个峰值位置。峰值检测器410包括阴影缓冲器412、最大值探测器(finder)414和局部零位器416。当过程第一次(first)开始时,阴影缓冲器412存储阴影权重Wsh 240。阴影缓冲器412在峰值检测器410的操作期间被更新。最大值探测器414识别或定位阴影权重Wsh的最大绝对值的位置。该位置对应于N个峰值位置之一。在峰值的位置被发现或识别之后,零位器416使位于该位置周围的阴影权重Wsh的子集为零。该子集以该位置为中心。子集的宽度W取决于稀疏滤波器的大小M。在一个实施例中,M的范围是从80-130,而W=50。被使得为零的子集然后被写回到阴影缓冲器412,并且最大值探测器414继续寻找被更新的阴影缓冲器412中的下一个峰值。每次检测到峰值位置时,峰值检测器410保存该位置以在稀疏权重更新中被使用。典型的是,待检测的峰值数等于稀疏权重Wmi的集合数。
根据N个峰值的位置,稀疏权重更新器420基于或使用阴影权重Wsh来替换N个集合的稀疏权重Wmi。稀疏权重更新器420包括阈值比较器430、第一更新器441、第二更新器442、第三更新器443和稀疏权重缓冲器450。
阈值比较器430比较峰值的位置j与低位置阈值TL和高位置阈值TH以确定如何替换N个集合的稀疏权重。基本上稀疏权重Wmi(k)是用阴影权重Wsh的子集或多个子集来替换的。子集由位置索引来识别如下:
第一子集:对于索引k在1和M1+ND之间,权重ND是耦合距离。
第二子集:对于索引k在1和ND之间
第三子集:对于索引k在M-N1+1≤k≤M之间
第四子集:对于索引k在j-M1/2和j+M1/2之间,即j是峰值位置并且是第四子集的中心。
稀疏滤波器Wmi(k)的替换或更新是根据以下方程来进行的:
如果j≤TL,则
对于1≤k≤M1+ND(第一子集),Wmi(k)=Wsh(k),且
否则Wmi(k)=0
如果j≥TH,则
对于1≤k≤ND(第二子集)和M-N1+1≤k≤M(第三子集),Wmi(k)=Wsh(k),且
否则Wmi(k)=0
如果TH<j<TL,则
对于1≤k≤ND(第二子集)和j-M1/2≤k≤j+M1/2(第四子集),Wmi(k)=Wsh(k),且
否则Wmi(k)=0。
响应于所述位置小于低位置阈值TL,第一更新器441用阴影权重的第一阴影子集来替换N个集合的稀疏权重之一的第一稀疏子集。响应于所述位置超过高位置阈值TH,第二更新器442分别用阴影权重的第二阴影子集和第三阴影子集来替换N个集合的稀疏权重之一的第二稀疏子集和第三稀疏子集。响应于所述位置在低和高位置阈值TL和TH之间,第三更新器443分别用阴影权重的第二阴影子集和第四阴影子集来替换N个集合的稀疏权重之一的第二稀疏子集和第四稀疏子集。第四阴影子集以所述位置为中心。
图5是说明依照本发明一个实施例执行回波消除的过程500的流程图。
一旦START,过程500分别使用主滤波器和阴影滤波器产生主输出和阴影输出(方块510)。这是通过延迟远端输入并使经延迟的远端输入乘以相应的滤波器权重而实现的。接下来,过程500计算主误差和阴影误差(方块515)。这是通过分别从主输出和阴影输出减去近端输入而实现的。
然后,过程500估算远端输入中的短期和长期功率(方块520)。接下来,过程500确定阴影更新是否被使能(方块525)。这是通过将短期功率比较于长期功率的指示而实现的。如果阴影更新被使能,则过程500使用AP适应性规则来更新阴影权重(方块535)。否则,过程500保持相同的阴影权重(方块530)。
接下来,过程500产生主耦合因子和阴影耦合因子(方块540)。然后,过程500确定是否已检测到双通话(方块545)。这是通过将主耦合因子比较于阴影耦合因子的指示而实现的。如果检测到双通话的状况,则过程500更新N个集合的稀疏权重(方块560)并且过程500被终止。否则,过程500估算主误差功率和阴影误差功率(方块550)。然后,过程500确定是否已检测到脉冲响应变化状况(方块555)。这是通过将主误差功率比较于阴影误差功率的指示而实现的。如果已检测到脉冲响应变化状况,则过程500更新N个集合的稀疏权重(方块560),然后被终止。否则,过程500保持N个集合的稀疏权重相同,然后被终止。
图6是说明依照本发明一个实施例更新图5中所示的稀疏权重集合的过程560的流程图。
一旦START,过程560识别或定位并保存阴影权重的最大绝对值的位置(方块610)。该位置对应于阴影权重中的峰值位置。接下来,过程560使以该位置为中心的该位置周围的局部子集为零(方块615)。然后,过程560确定是否所有的集合已被处理(方块620)。如果没有,过程560返回方块610以继续寻找下一个峰值位置。
如果所有的集合已被处理,则过程560从第一峰值位置开始(方块625)。接下来,过程560确定峰值位置是否小于低位置阈值(方块630)。如果是,则过程560更新稀疏权重的第一子集(方块635)并前进到方块660。否则,过程560确定峰值位置是否大于或超过高位置阈值(方块640)。如果是,过程560更新稀疏权重的第二和第三子集(方块645)并前进到方块660。否则,过程560更新稀疏权重的第二和第四子集(方块650)。
然后,过程560确定是否所有的峰值位置已被处理(方块660)。如果没有,过程560转到下一个峰值位置(方块665)并且返回方块630。如果所有的峰值位置已被处理,则过程560被终止。
图7-11示出仿真的结果。其参数为,对于首先的100,000次迭代,对在8、160和320的延迟的三个类型1滤波器,M=512+32=544和M1=140,然后是在88、240和400的延迟的滤波器。
图7是说明依照本发明一个实施例在200K次迭代之后组合主权重的图。图8是说明依照本发明一个实施例的平滑均方误差(MSE)的图。图9是说明依照本发明一个实施例的使用延迟系数的双通话的状况的图。在该图中,z=βm(k)-2βs(k)。图10是说明依照本发明一个实施例的平滑LEM变化状况的图。在该图中,z1=epm-r3eps。图11是说明依照本发明一个实施例在迭代104800到106800处的平滑LEM变化状况的图。
图8示出存在几个阶段。在从0到35,000次迭代的第一阶段内,阴影滤波器学习通道脉冲响应。如在从大约40dB减小到0dB的主滤波器MSE(esqta)中所反映的,控制逻辑将适当的权重传递给主滤波器。在从35,000到99,000次迭代的第二阶段内,如图9和10中所示,双通话出现并且被双通话检测器和LEM检测器检测到。主通道滤波器权重被冻结,而Wsh适应于(adapt on)双通话。在从100,000到125,000次迭代的第三阶段内,双通话的状况消失并且LEM变化出现。阴影滤波器学习LEM变化并将新的权重设置传递给主通道滤波器。该步骤对于双通话的检测(见图9)最初出现于104,000次迭代左右而对于LEM变化(见图11)最初出现于103K次迭代左右并且再次出现于105K次迭代左右。因此,初始的LEM变化的瞬态在大约5/8秒内消失。随着MSE继续减小(见图8),阴影滤波器继续适应。注意与第一学习阶段相比,消除现在是从40dB到10dB。然而,这可能是由于减小了学习时间。附加的训练采样可更多地缩减消除执行。
尽管已根据几个实施例描述了本发明,本领域的技术人员将认识到,本发明并不局限于所述的实施例,而是可在所附权利要求的精神和范围内实施以修改和变更。说明书由此被看作是说明性的而不是局限性的。
Claims (23)
1.一种用于稀疏通道的回波消除器中的控制逻辑电路,包括:
使能发生器,其响应于检测到回波通道中的双通话的状况和脉冲响应变化的至少一个而确立稀疏更新使能,该回波通道接收一个远端输入和一个近端输入,所述使能发生器包括:
双通话检测器,用于检测双通话的状况,以及
变化检测器,用于检测脉冲响应变化;以及
主权重更新器,其根据所述稀疏更新使能而将阴影滤波器的阴影权重传递给主滤波器的多个稀疏滤波器中的多个集合的稀疏权重,主滤波器非适应地过滤远端输入以产生主输出,阴影滤波器适应地过滤远端输入以产生阴影输出,并且估算稀疏权重的峰值的延迟,所述主权重更新器包括:
峰值检测器,其检测阴影权重中多个峰值的多个峰值位置;以及
稀疏权重更新器,其根据多个峰值位置基于阴影权重来替换多个集合的稀疏权重。
2.权利要求1的电路,其中峰值检测器包括:
阴影缓冲器,其存储阴影权重;
最大值探测器,其识别在阴影缓冲器中存储的阴影权重的最大绝对值的位置,该位置对应于多个峰值位置之一;以及
零位器,其使位于所述位置周围的阴影权重的子集为零,所述为零的子集被写回到阴影缓冲器。
3.权利要求2的电路,其中稀疏权重更新器包括:
阈值比较器,其比较所述位置与低位置阈值和高位置阈值;
第一更新器,其响应于所述位置小于低位置阈值而用阴影权重的第一阴影子集来替换多个集合的稀疏权重之一的第一稀疏子集;
第二更新器,其响应于所述位置超过高位置阈值而分别用阴影权重的第二阴影子集和第三阴影子集来替换多个集合的稀疏权重之一的第二稀疏子集和第三稀疏子集;以及
第三更新器,其响应于所述位置在低和高位置阈值之间而分别用阴影权重的第二阴影子集和第四阴影子集来替换多个集合的稀疏权重之一的第二稀疏子集和第四稀疏子集,所述位置在第四阴影子集的中心。
4.权利要求3的电路,其中稀疏权重更新器进一步包括:
稀疏权重缓冲器,其存储多个集合的稀疏权重。
5.权利要求1的电路,其中使能发生器包括:
短期功率估算器,其估算远端输入的短期功率;
长期功率估算器,其估算远端输入的长期功率;以及
适应性使能器,其响应于短期功率大于长期功率的指示而使能阴影滤波器的适应性。
6.权利要求1的电路,其中使能发生器还包括:
延迟系数发生器,其计算用于主滤波器的主耦合因子和用于阴影滤波器的阴影耦合因子;以及
其中双通话检测器响应于主耦合因子小于阴影耦合因子的指示而检测双通话的状况。
7.权利要求6的电路,其中使能发生器进一步包括:
主误差功率估算器,其估算从主输出和近端输入确定的主误差的主误差功率;
阴影误差功率估算器,其估算从阴影输出和近端输入确定的阴影误差的阴影误差功率;以及
其中变化检测器响应于主误差功率超过阴影误差功率的指示而检测脉冲响应变化。
8.权利要求7的电路,其中使能发生器进一步包括:
组合器,其组合所检测的双通话和脉冲响应变化以确立稀疏更新使能。
9.权利要求5的电路,其中适应性使能器使能阴影滤波器的仿射投影(AP)适应性。
10.一种控制用于稀疏通道的回波消除器的方法,包括:
响应于检测到回波通道中的双通话的状况和脉冲响应变化的至少一个而确立稀疏更新使能,该回波通道接收一个远端输入和一个近端输入,确立稀疏更新使能包括:
检测双通话的状况,以及
检测脉冲响应变化;以及
根据所述稀疏更新使能而将阴影滤波器的阴影权重传递给主滤波器的多个稀疏滤波器中的多个集合的稀疏权重,主滤波器非适应地过滤远端输入以产生主输出,阴影滤波器适应地过滤远端输入以产生阴影输出,并且估算稀疏权重的峰值的延迟,所述传递包括:
检测阴影权重中多个峰值的多个峰值位置;以及
根据多个峰值位置基于阴影权重来替换多个集合的稀疏权重。
11.权利要求10的方法,其中检测多个峰值位置包括:
在阴影缓冲器中存储阴影权重;
识别在阴影缓冲器中存储的阴影权重的最大绝对值的位置,该位置对应于多个峰值位置之一;以及
使位于所述位置周围的阴影权重的子集为零,所述为零的子集被写回到阴影缓冲器。
12.权利要求11的方法,其中替换多个集合的稀疏权重包括:
比较所述位置与低位置阈值和高位置阈值;
响应于所述位置小于低位置阈值而用阴影权重的第一阴影子集来替换多个集合的稀疏权重之一的第一稀疏子集;
响应于所述位置超过高位置阈值而分别用阴影权重的第二阴影子集和第三阴影子集来替换多个集合的稀疏权重之一的第二稀疏子集和第三稀疏子集;以及
响应于所述位置在低和高位置阈值之间而分别用阴影权重的第二阴影子集和第四阴影子集来替换多个集合的稀疏权重之一的第二稀疏子集和第四稀疏子集,所述位置在第四阴影子集的中心。
13.权利要求12的方法,其中替换多个集合的稀疏权重进一步包括:
存储多个集合的稀疏权重。
14.权利要求10的方法,其中确立稀疏更新使能包括:
估算远端输入的短期功率;
估算远端输入的长期功率;以及
响应于短期功率超过长期功率的指示而使能阴影滤波器的适应性。
15.权利要求10的方法,其中确立稀疏更新使能还包括:
计算用于主滤波器的主耦合因子和用于阴影滤波器的阴影耦合因子;以及
其中响应于主耦合因子小于阴影耦合因子的指示而检测双通话的状况。
16.权利要求15的方法,其中确立稀疏更新使能进一步包括:
估算从主输出和近端输入确定的主误差的主误差功率;
估算从阴影输出和近端输入确定的阴影误差的阴影误差功率;以及
其中响应于主误差功率超过阴影误差功率的指示而检测脉冲响应变化。
17.权利要求16的方法,其中确立稀疏更新使能进一步包括:
组合所检测的双通话和脉冲响应变化以确立稀疏更新使能。
18.权利要求14的方法,其中使能阴影滤波器的适应性包括使能阴影滤波器的仿射投影(AP)适应性。
19.一种用于稀疏通道的回波消除的系统,包括:
主滤波器,其具有有多个集合的稀疏权重的多个稀疏滤波器,以非适应地过滤回波通道中的远端输入,主滤波器产生主输出;
阴影滤波器,其具有阴影权重,用于在回波通道中过滤远端输入以适应脉冲响应变化,阴影滤波器产生阴影输出并估算稀疏权重中的峰值的延迟;以及
控制逻辑电路,其被耦合于主滤波器和阴影滤波器,检测双通话和脉冲响应变化之一,该控制逻辑电路包括:
使能发生器,响应于检测到双通话的状况和脉冲响应变化的至少一个而确立稀疏更新使能,所述使能发生器包括:
双通话检测器,用于检测双通话的状况,以及
变化检测器,用于检测脉冲响应变化;以及
主权重更新器,其被耦合于使能发生器,根据稀疏更新使能将阴影滤波器的阴影权重传递给多个集合的稀疏权重,所述主权重更新器包括:
峰值检测器,其检测阴影权重中多个峰值的多个峰值位置;以及
稀疏权重更新器,其根据多个峰值位置基于阴影权重来替换多个集合的稀疏权重。
20.权利要求19的系统,进一步包括:
发送输入解码器,其从近端接收经编码的语音;以及
发送输出编码器,其提供语音压缩。
21.权利要求20的系统,进一步包括:
接收器输入解码器,其解压缩从远端接收的语音数据;以及
接收器输出编码器,其编码来自回波通道的线性数据采样,经编码的线性数据采样被发送给近端。
22.权利要求21的系统,其中发送输入解码器是μ定律和A定律解码器之一。
23.权利要求22的系统,其中发送输出编码器压缩来自回波通道的语音数据,经压缩的语音数据被发送给远端。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/304,641 US7471788B2 (en) | 2002-11-25 | 2002-11-25 | Echo cancellers for sparse channels |
US10/304641 | 2002-11-25 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1503466A CN1503466A (zh) | 2004-06-09 |
CN1503466B true CN1503466B (zh) | 2010-05-12 |
Family
ID=32325268
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2003101183656A Expired - Fee Related CN1503466B (zh) | 2002-11-25 | 2003-11-25 | 控制稀疏通道中回波消除的方法、电路和系统 |
Country Status (7)
Country | Link |
---|---|
US (1) | US7471788B2 (zh) |
JP (1) | JP4663323B2 (zh) |
KR (1) | KR100738713B1 (zh) |
CN (1) | CN1503466B (zh) |
AU (1) | AU2003286814A1 (zh) |
TW (1) | TWI268093B (zh) |
WO (1) | WO2004049584A1 (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070189505A1 (en) * | 2006-01-31 | 2007-08-16 | Freescale Semiconductor, Inc. | Detecting reflections in a communication channel |
US8565416B2 (en) * | 2006-03-17 | 2013-10-22 | Cisco Technology, Inc. | Cache-based echo canceller |
US8320574B2 (en) | 2006-04-20 | 2012-11-27 | Hewlett-Packard Development Company, L.P. | Methods and systems for reducing acoustic echoes in communication systems |
US20080091415A1 (en) * | 2006-10-12 | 2008-04-17 | Schafer Ronald W | System and method for canceling acoustic echoes in audio-conference communication systems |
CN102859591B (zh) * | 2010-04-12 | 2015-02-18 | 瑞典爱立信有限公司 | 用于语音编码器中的噪声消除的方法和装置 |
CN104392727A (zh) * | 2014-11-12 | 2015-03-04 | 华为技术有限公司 | 音频信号处理方法和相关装置 |
KR20170032603A (ko) | 2015-09-15 | 2017-03-23 | 삼성전자주식회사 | 전자 장치, 그의 반향 신호 제거 방법 및 비일시적 컴퓨터 판독가능 기록매체 |
CN109766646B (zh) * | 2019-01-16 | 2021-06-04 | 北京大学 | 一种基于稀疏通道回波数据重建的超声成像方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5343522A (en) * | 1992-02-05 | 1994-08-30 | Northern Telecom Limited | Adaptive sparse echo canceller using a sub-rate filter for active tap selection |
WO1998027668A1 (en) * | 1996-12-19 | 1998-06-25 | Northern Telecom Limited | Method and apparatus for computing measures of echo |
US6028929A (en) * | 1997-11-14 | 2000-02-22 | Tellabs Operations, Inc. | Echo canceller employing dual-H architecture having improved non-linear echo path detection |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2947093B2 (ja) * | 1994-11-02 | 1999-09-13 | 日本電気株式会社 | 適応フィルタによるシステム同定の方法および装置 |
US5610909A (en) * | 1995-05-31 | 1997-03-11 | Lucent Technologies Inc. | Multistage echo canceler including time variation compensation |
KR100233463B1 (ko) * | 1997-03-07 | 1999-12-01 | 윤종용 | 반향제거장치 및 방법 |
US5857167A (en) * | 1997-07-10 | 1999-01-05 | Coherant Communications Systems Corp. | Combined speech coder and echo canceler |
US6240180B1 (en) * | 1997-11-14 | 2001-05-29 | Tellabs Operations, Inc. | Echo canceller employing dual-H architecture having split adaptive gain settings |
US6909782B2 (en) | 2000-09-08 | 2005-06-21 | Intel Corporation | Fast converging affine projection based echo canceller for sparse multi-path channels |
GB2379369B (en) * | 2001-08-29 | 2005-06-29 | Zarlink Semiconductor Inc | Subband echo location and double-talk detection in communication systems |
US20030219113A1 (en) | 2002-05-21 | 2003-11-27 | Bershad Neil J. | Echo canceller with double-talk and channel impulse response adaptation |
-
2002
- 2002-11-25 US US10/304,641 patent/US7471788B2/en not_active Expired - Fee Related
-
2003
- 2003-10-30 KR KR1020057009342A patent/KR100738713B1/ko not_active IP Right Cessation
- 2003-10-30 WO PCT/US2003/034688 patent/WO2004049584A1/en active Application Filing
- 2003-10-30 JP JP2004555361A patent/JP4663323B2/ja not_active Expired - Fee Related
- 2003-10-30 AU AU2003286814A patent/AU2003286814A1/en not_active Abandoned
- 2003-11-13 TW TW092131814A patent/TWI268093B/zh not_active IP Right Cessation
- 2003-11-25 CN CN2003101183656A patent/CN1503466B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5343522A (en) * | 1992-02-05 | 1994-08-30 | Northern Telecom Limited | Adaptive sparse echo canceller using a sub-rate filter for active tap selection |
WO1998027668A1 (en) * | 1996-12-19 | 1998-06-25 | Northern Telecom Limited | Method and apparatus for computing measures of echo |
US6028929A (en) * | 1997-11-14 | 2000-02-22 | Tellabs Operations, Inc. | Echo canceller employing dual-H architecture having improved non-linear echo path detection |
Also Published As
Publication number | Publication date |
---|---|
TW200423698A (en) | 2004-11-01 |
AU2003286814A1 (en) | 2004-06-18 |
KR100738713B1 (ko) | 2007-07-12 |
JP4663323B2 (ja) | 2011-04-06 |
WO2004049584A1 (en) | 2004-06-10 |
KR20050089009A (ko) | 2005-09-07 |
TWI268093B (en) | 2006-12-01 |
US20040101131A1 (en) | 2004-05-27 |
CN1503466A (zh) | 2004-06-09 |
JP2006507757A (ja) | 2006-03-02 |
US7471788B2 (en) | 2008-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4142715B2 (ja) | 自己回帰チャネルモデル化を用いて搬送周波数オフセットとフェージング・レートとを推定する方法および装置 | |
JPH04291522A (ja) | 少なくとも2つの受信ブランチからなる受信器 | |
CN1503466B (zh) | 控制稀疏通道中回波消除的方法、电路和系统 | |
EP0950305B1 (en) | Apparatus and method for sequential maximum likelihood estimating communications signals using whitening path metrics | |
CA2284989C (en) | Combined parallel adaptive equalizer/echo canceller | |
US6909782B2 (en) | Fast converging affine projection based echo canceller for sparse multi-path channels | |
WO1998038582A1 (en) | Adaptive dual filter echo cancellation | |
CN101711049B (zh) | 一种基于干扰消除的选径方法和装置 | |
CN1473403B (zh) | 用于基于回波消除器的快速收敛仿射投影的装置和方法 | |
KR100931118B1 (ko) | 등화 제어를 위한 방법 및 장치 | |
CN100512038C (zh) | 用于回波抵消器的噪声匹配 | |
KR20010006858A (ko) | 심볼 검출기에 대한 소프트 출력 측정치 발생 | |
Lawrence et al. | Multiprocessor implementation of adaptive digital filters | |
JPS62265824A (ja) | 双方向デイジタル伝送システム用エコ−・キヤンセラ | |
Adachi | Reduced state transition Viterbi differential detection of M-ary DPSK signals | |
CN103346822B (zh) | 用于解调器的可转换均衡器 | |
KR20050114735A (ko) | 데이터 신호의 효율적인 디인터리빙 및 다이버시티-결합을용이하게 하는 방법 및 장치 | |
Ganguly et al. | Performance evaluation of a packet voice system | |
Green et al. | Dynamic channel order estimation algorithm | |
Kafedziski | Joint source channel coding of images over frequency selective fading channels with feedback using dct and multicarrier block pulse amplitude modulation | |
Zhou et al. | Echo cancellation research of channel estimation based on PN sequence | |
KR20030010584A (ko) | 등화기 장치 및 등화기 장치에 대한 상태 감소의 방법 | |
Khatib et al. | Study of parameters effect on the performance of precoding and equalizer sharing systems | |
KR20000003240A (ko) | 무선통신 시스템의 송수신 장치 및 그 방법 | |
Agarwal et al. | Computationally efficient coherent correlator design for DVB-S2 receiver |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100512 Termination date: 20171125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |