CN1473403A - 用于基于回波消除器的快速收敛仿射投影的装置和方法 - Google Patents

用于基于回波消除器的快速收敛仿射投影的装置和方法 Download PDF

Info

Publication number
CN1473403A
CN1473403A CNA018185274A CN01818527A CN1473403A CN 1473403 A CN1473403 A CN 1473403A CN A018185274 A CNA018185274 A CN A018185274A CN 01818527 A CN01818527 A CN 01818527A CN 1473403 A CN1473403 A CN 1473403A
Authority
CN
China
Prior art keywords
weight vector
list entries
peak value
channel weight
channel
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
CNA018185274A
Other languages
English (en)
Other versions
CN1473403B (zh
Inventor
N��ɳ����
N·贝尔沙德
A·比斯特
S·希
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.)
Intel Corp
Original Assignee
Intel Corp
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 Intel Corp filed Critical Intel Corp
Publication of CN1473403A publication Critical patent/CN1473403A/zh
Application granted granted Critical
Publication of CN1473403B publication Critical patent/CN1473403B/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/02Details
    • H04B3/20Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other
    • H04B3/23Reducing 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

Abstract

一个自适应滤波器,使用仿射投影(affineprojection)(AP)修正估计回波信道的信道权向量。该回波信道接收一个发送输入序列和一个接收输入序列。当自适应滤波器分别运行在一个第一自适应模式和一个第二自适应模式时,该信道权向量具有第一和第二长度。一个延迟估计器,在第一自适应模式里使用自适应滤波器确定回波信道中的延迟。

Description

用于基于回波消除器的快速收敛仿射投影的装置和方法
相关申请
本申请要求2000年9月8日提交的美国临时专利申请No.60/231,420的权利(代理人登记号码004419.P015Z)。
技术领域
本发明涉及信号处理。特别地,本发明涉及回波消除。
背景技术
在电话通信里回波通常是不希望的。回波由许多原因引起。这里包括从扩音器回到麦克风的信号的多次反射,扩音器和电话之间的直接声耦合以及环境噪声。回波消除是减少回波的不希望影响的一种技术。回波消除器估计回声路径的脉冲响应并且产生一个该回波的估计值。估计的回波被从近端信号中减去。典型地,因为回声路径通常是未知的并且随机地随时间变化的,所以使用自适应滤波器用来估计该回波。
一种用于回波消除的现有方法是规格化的最小均方(NLMS)方法。该方法试图最小化平方误差的预期值。然而,NLMS方法具有许多缺点。一个显著的缺点是对于有色输入的缓慢收敛速度。
因此,需要一种高效的技术执行具有比NLMS方法快速的收敛速度的回波消除。
附图说明
本发明的特征和优点从下列发明的详细说明中便会变得更清楚,其中:
图1是说明一种系统的一个框图,其中本发明的一个实施例可以在其中实现。
图2是说明按照本发明的一个实施例的、用于图1中所示的基于AP的回波消除器的系统模型的一个框图。
图3是说明按照本发明的一个实施例的、图1中所示的基于AP的回波消除器的一个框图。
图4是说明按照本发明的一个实施例的、用于图3中所示的基于AP的回波消除器的自适应滤波器的一个框图。
图5是说明按照本发明的一个实施例的、用于图2中所示的基于AP的回波消除器的延迟估计器的一个框图
图6是说明按照本发明一个实施例的、用于回波消除的处理的流程图。
图7是说明按照本发明的一个实施例的、使用图6中所示的AP修正进行自适应滤波的处理的一个流程图。
图8是说明按照本发明一个实施例的、图6中所示的延迟估计的处理流程图。
图9A是按照本发明的一个实施例,说明使用M=256的AP(2)的回波消除器中衰减的分贝的框图。
图9B是按照本发明的一个实施例,说明使用M=256的NLMS的回波消除器中衰减的分贝的框图。
图10A是按照本发明的一个实施例,说明使用M=512的AP(2)的回波消除器中衰减的分贝的框图。
图10B是按照本发明的一个实施例,说明使用M=512的NLMS的回波消除器中衰减的分贝的框图。
图11A是按照本发明的一个实施例,说明使用M=1024的AP(2)的回波消除器中衰减的分贝的框图。
图11B是按照本发明的一个实施例,说明使用M=1024的NLMS的回波消除器中衰减的分贝的框图。
图12A是按照本发明的一个实施例,说明在第二自适应模式里使用M1=1024和M2=256的AP(2)的第一种方法的回波消除器中的分贝衰减的框图。
图12B是按照本发明的一个实施例,说明在第二自适应模式里使用M1=1024和M2=256的AP(2)的第二种方法的回波消除器中的分贝衰减的框图。
图13是说明按照本发明一个实施例,M+500迭代上的加权的框图。
具体实施方式
在下面描述中,为了解释起见,为了提供对本发明的彻底理解,提出许多的详述。然而,某些细节并不是实施本发明所必须的。在其它实例里,为了不使本发明模糊,公知的电结构与电路以方框图方式显示。
本发明可以通过硬件、软件、固件、微代码或者任何它们的任意组合实现。当以软件、固件或者微代码方式实现时,本发明的元件是执行必要任务的程序代码或者代码段。一个″代码段″可以表示一个过程、一个函数、一个子程序、一个程序、一个例行程序、一个子例行程序、一个模块、一个软件包、一个类程或者任何指令、数据结构或者程序语句的综合。一个代码段可以通过传递和/或接收信息、数据、自变量、参数或者存储内容连接到另一个代码段或者一个硬件电路。信息、自变量、参数、数据等等可以通过任何合适的方法包括存储器共享、报文传递、令牌传送、网络传输等等传递、转送或者发送。程序或者代码段可以保存在一个处理器可读媒介中或者通过一个包含在载波中的计算机数据信号或者由载波调制的信号在传输介质上传送。″机器可读的媒体″可以包括可以存储或者转送信息的任何媒介。处理器可读媒介的例子包括一种电子电路、一个半导体存储器装置、一个只读存储器(ROM)、一个快速[闪]存储器、一种可擦除ROM(EROM)、一个软盘磁盘、一个高密度磁盘(CD-ROM)、一种光盘、一种硬盘、一种光纤的媒介、一个无线电频率(RF)链路等等。计算机数据信号可以包括可以基于传输介质传播的任何信号,例如基于电子网络信道、光纤、空气、电磁、RF链路等等。代码段可以通过例如Internet、Intranet等等的计算机网络下载。
应当注意到本发明可以描述为一种处理过程,该处理过程通常被描述为一种流程图、一种程序框图、一种结构图或者一种方框图。虽然流程图可以以一种顺序进程描述该运行,许多运行可以并行执行或者同时执行。另外,该运行的顺序可以重新排列。当它的运行完成时处理过程停止。处理过程可以相应于一个方法、一个函数、一个过程、一个子例行程序、一个子程序等等。当一个处理过程相应于一个函数时,它的结束部分对应该函数返回到调用该函数或者该主函数。
本发明是提高回波消除性能的一个方法和装置。一种自适应滤波器使用仿射投影(AP)修正规则以便修正第一和第二自适应模式中的估计值。在第一自适应模式中,自适应滤波器使用一个长滤波长度估计信道权向量来确定容积延迟(bulk delay)。在容积延迟估算之后,自适应滤波器切换到第二自适应模式,基于估计容积延迟使用一个短滤波长度来估计来信道权向量。该技术显示出比NLMS方法快速的收敛速度。
图1是说明一种系统100的一个框图,其中本发明的一个实施例可以在其中实现。系统100包括一个发送输入解码器110、一个回波信道120、一个发送输出解码器130、一个接收输入解码器140、一个接收输出编码器150,一个网络145,一个发送输入解码器160,一个回波信道170,一个发送输出解码器180,一个接收输入解码器190和一个接收输出编码器195。
该发送输入解码器110从第一近端接收编码的语音并且将编码的语音解码成为线性语音数据Sin。在一个实施例里,发送输入解码器110是一个μ-Law/A-Law解码器。回波信道120包括一个使用仿射投影(AP)125的回波消除器。基于AP的回波消除器125从线性数据抽样Sin删除一个回波估计信号以便产生线性数据样本Sout。发送输出编码器130在分组化以前提供话音压缩。在一个实施例中,发送输出编码器130是一个G.7xx编码器,它使用用于低比特率话音(LBRV)的压缩标准中的任何一个压缩来自回波信道120的语音数据Sout,该压缩标准包括国际电信联盟(ITU)-T国际性地标准化的G.7xx系列。压缩的语音数据通过一个网络发给远端。接收输入解码器140解压缩来源于第一远端通过网络145接收的语音数据。解压缩技术与用于发送输出编码器130的压缩一致。回波信道120接收来自接收输入解码器140的Rin并且发送Rout线性数据抽样。接收输出编码器150编码线性数据抽样Rout成为将被发到第一近端的μ-Law和A-Iaw编码的语音。
网络145是具有发送且接收分组的数据性能的任何网络,其中分组的数据可以来自也可以是发送到发送输出解码器130,发送输入解码器160、接收输入解码器140、接收输出解码器195。网络145可以是Internet、intranet、extranet、局域网(LAN)或者广域网(WAN)。发送输入解码器160从网络145接收编码的语音并且解码编码的语音成为线性语音数据Sin。在一个实施例中,发送输入解码器160是一个μ-Law/A-Law解码器。回波信道170包括一个使用仿射投影(AP)175的回波消除器。基于AP的回波消除器175从线性数据抽样Sin删除一个回波估计信号以便产生线性数据样本Sout。发送输出编码器180在分组化以前提供语言压缩。在一个实施例中,发送输出编码器180是一个G.7xx编码器,它使用用于低比特率话音(LBRV)的压缩标准中的任何一个压缩来自回波信道170的语音数据Sout,该压缩标准包括国际电信联盟(ITU)-T国际性地标准化G.7xx系列。压缩的语音数据被发给第二远端上的一个接收设备。接收输入解码器190解压缩来源于第二远端的语音数据。解压缩技术与用于发送输出编码器180的压缩一致。回波信道170接收来自接收输入解码器190的Rin并且发送Rout线性数据抽样。接收输出编码器190编码线性数据抽样Rout成为将被发送到第二近端、到网络145的μ-Law和A-Iaw编码的语音。在一个实施例中,发送输入解码器160,回波信道170,发送输出解码器180,接收输入解码器190以及接收输出编码器195被集成到一个数字信号处理机165中。
图2是说明对于图1中所示的基于AP的回波消除器125的系统模型的一个框图,图1按照本发明的一个实施例。系统模型200包括一个设备(plant)210、一个加法器220、一个自适应滤波器230以及一个减法器240。
该系统模型200模仿回波消除处理过程。回波消除处理过程可以作为一个系统识别问题模仿,其中一个未知的线性系统将被确定。该系统的输入是一个输入序列u(k)。设备210表征系统功效。设备210可以作为具有权向量W=[W1,W2,...,WN]的有限脉冲响应(FIR)滤波器被模仿。设备输出由下式给出:
公式(1)
其中WT是权向量W的变换。权向量W是未知的。
加法器220添加一个随机噪声序列v(k)到设备输出y(k)以便产生期望输出值d(k)。自适应滤波器230作为具有权向量=[W1,W2...,Wn]的FIR滤波器被模仿。减法器240从自适应滤波器的输出中减去期望输出值d(k)以便产生一个误差e(k)。误差e(k)用来调整或者修正权向量因此误差e(k)在某些目标函数(例如最小均方)下被最小化。权向量被修正因此它接近权向量W。因此,权向量表示权向量W的一个估计值。
图3是说明按照本发明的一个实施例的、对于图1中所示的基于AP的回波消除器125的一个框图。回波消除器125包括一个自适应滤波器310以及一个延迟估计器320。
自适应滤波器310使用仿射投影(AP)修正估计回波信道120(图1)的信道权向量W。如图1所述,回波信道120接收一个发送输入序列Sin以及一个接收输入序列Rin。自适应滤波器310以两种模式运行:一个第一自适应模式312以及一个第二自适应模式314。第一自适应模式312是自适应滤波器310用来计算回波信道120的容积延迟的模式。第二自适应模式314是自适应滤波器310调整或者修正信道权向量W的模式。当自适应滤波器310分别地运行在第一自适应模式312和第二自适应模式314时,信道权向量W具有一个第一长度M1和一个第二长度M2。第一长度M1比第二长度M2长。例如,M1=1024而M2=256。可替换地,自适应滤波器310可以由具有滤波长度M1和M2的两个自适应的滤波器替代。
延迟估计器320使用由自适应滤波器310提供的估计的信道权向量W在回波信道120里确定容积延迟,自适应滤波器310运行在第一自适应模式312。延迟估计器320提供估算的延迟给自适应滤波器310。自适应滤波器310使用估算的延迟在第二自适应模式314里相应地定位该滤波器。
自适应滤波器310使用AP修正规则估计信道权向量。AP方法用来加速正常化的最小均方(NLMS)技术的收敛,尤其对于有色的输入。AP基于多个过去的输入向量修正加权,同时NLMS基于最近的输入向量修正加权。AP修正规则在下文中描述。
输入序列u(k)可以作为标记为AR(P)的阶(order)P的一个自动回归处理过程模仿:
公式(2)
其中z(k)是具有单元方差的白序列。
假设P是一个已知的先验,u(K)的取样可写为(Mx1)列矢量u(K),或者:
公式(3)
AR(P)过程可被写为:
公式(4)
其中U(K)是过去矢量P的集合
公式(5)
z(K)是白随机序列取样的(Mx1)列矢量:
公式(6)
a的参数的最小二乘估计量由下式给出
公式(7)
其中UT(K)*U(K)假定为P等级(rank P),*表示相乘。
AP循环修正规则对于μ=1由下列运行确定:
公式(8)
公式(9)
公式(10)
其中w(k)是由自适应滤波器310估计的信道权向量。
已经表明如果u(k)是一个阶P的AR处理过程,(k)是其AR系数的最小二乘估计量并且φ(k)≈Z(k)。换句话说,φ(k)是一个矢量,它的元素是白随机序列的估计值。
图4是说明按照本发明的一个实施例的、用于图3中所示的基于AP的回波消除器310的自适应滤波器310的一个框图,其中图3。基于AP的自适应滤波器310包括一个自动回归(AR)系数估计器410,一个随机估计器420,一个误差估计器430和一个加权修正器440。
AR系数估计器410使用过去接收输入序列和接收序列以第一修正速度估计AR系数向量a(k)。AR系数估计器410使用等式(7)产生AR系数向量(k)。
随机估计器420使用估计的AR系数向量以第二修正速度估计随机序列φ(k)。随机估计器420使用等式(8)确定随机序列φ(k)。
误差估计器430以第二修正速度使用发送序列、接收序列和估计的信道权向量估计误差。误差估计器430使用等式(9)计算误差e(k)。
权修改器440以第二修正速度使用估计误差和估计的随机序列修正信道权向量。加权修正器440使用等式(10)修正信道加权。
第一和第二修正速度可以相同或者不同。在一个实施例中,第一修正速度比第二修正速度慢。第二修正速度是修正环每个迭代上的速率。第一修正速度是每个K迭代上的速率,其中K是一个正整数。在一个实施例中,K等于100。换句话说,AR系数估计器410在每个K迭代上产生一个新的结果同时随机估计器420,误差估计器430和加权修正器440在每个迭代上产生新的结果。
图5是说明按照本发明的一个实施例的、图2中所示的基于AP的回波消除器的延迟估计器320的一个框图,其中图2。延迟估计器320包括一个峰值定位器510,一个峰值消除器520和一个上升沿定位器530。
延迟估计器320在回波信道120(图1)里从估计的信道权向量(k)中估计容积延迟。由延迟估计器320提供的估计延迟由自适应滤波器310(图3)使用用于自适应。延迟估计器320实质上在脉冲响应里设置许多峰值,脉冲响应作为由估计重量矢量(k)的分量提供。
峰值定位器510在搜索区域之内确定该峰值作为该加权的最大值。典型地L峰值被定位到L为从1到5的一个正整数。在一个实施例中,L等于5。L峰值位于彼此距离一个预定距离内。首先,最高峰值确定为信道权向量(k)中该加权的最高绝对值。因此,第二高峰值位于覆盖最高峰值区域的外部。该区域包括在最高峰值两个边上的LI抽样。在一个实施例中,LI=25。因此,第三个高峰值位于覆盖第二高峰值区域的外部。该区域包括在第二高峰值两个边上的L2抽样。在一个实施例中,L2=25。处理过程持续到所有的L峰值定位为止。
峰值消除器520消除一个由峰值定位器510定位的L峰值中的假峰值。当它的值小于阈值时,确定该假峰值。在一个实施例中,阈值等于β*最高峰值,其中β是0.4到0.8之间的一个数。在一个实施例中,β=0.6。
上升沿定位器530确定一个脉冲响应的上升沿。该上升沿是L峰值最小延迟上的峰值。在确定上升沿之后,该延迟作为上升沿位置和一个预定距离进行计算。在一个实施例中,预定距离近似等于25。上升沿定位器530的任务是配置回波响应的上升沿接近运行在第二自适应模式的自适应滤波器的第一抽头(tap)。
图6是说明按照本发明一个实施例的回波消除处理过程600的流程图。
在START,处理过程600预置迭代指数k=1,滤波长度M=M1,迭代数目Q=Q1,信道权向量(k),AR系数向量,随机序列φ(k),误差e(k),以及过去接收输入序列U(k)(方框610)。在一个实施例中,M1等于256、512或者1024。之后,处理过程600使用AP修正执行自适应滤波以便估计信道权向量(k)(方框620)。然后,处理过程600确定迭代指数k是否等于Q1(方框625)。如果是这样的话,处理过程600增加迭代指数k(方框627)并且返回到方框620。否则,处理过程600使用估计信道权向量(k)估计延迟(方框630)。
之后,处理过程600预置迭代指数k=Q1,滤波长度M=M2,迭代数目Q=Q2,信道权向量(k),AR系数向量,随机序列φ(k),误差e(k)以及过去接收输入序列U(k)(方框640)。M2典型地小于M1.。在一个实施例中,M1=1024,M2=256。然后,处理过程600使用AP修正执行自适应滤波以便估计信道权向量(k)(方框650)。在方框650里的计算和修正实质上与在方框620里的相同。这样,相同的滤波器实施例可被用于两个步骤。可替换地,可以使用两个相同的自适应滤波器。然后,处理过程600确定迭代指数k是否等于Q2(方框655)。如果不,处理过程600增加迭代指数k(方框657)并且返回到方框650。否则,处理过程600终止。
图7是说明按照本发明的一个实施例的、使用图6中所示的AP修正的自适应滤波处理过程620的一个流程图。处理过程620实质上与处理过程650相同。为了简便起见,仅使用参考号620。
在START,处理过程620计算短期平均功率,stavp,以及长期平均功率,Rinphat(方框710)。然后,处理过程620确定stavp是否小于Rinphat一个预定量(方框715)。在一个实施例中,该数量是20dB。如果是这样的话,处理过程620冻结估计的信道权向量(方框720)并且估计误差d(k),然后就终止。否则,处理过程620保存估计的接收矢量(k-1)(方框730)。
之后,处理过程620确定此时是否应该修正AR系数向量a(方框735)。在一个实施例中,第一修正速度相当于每个R迭代,其中R=100。如果不是修正的时机,估计的AR系数向量保持与保存在方框730中的以前的值相同。否则,AR系数向量按照等式(7)修正。之后,处理过程620使用等式(8)估计随机序列φ(k)(方框750)。处理过程620使用等式(9)估计误差e(k)。之后,处理过程620使用等式(10)修正信道权向量(k)。然后,处理过程620确定是否已经达到迭代的最大数(方框770)。如果没有,处理过程620增加迭代指数k(方框780)。否则,处理过程620终止。
图8是说明按照本发明一个实施例的图6中所示的延迟估计的处理过程630的流程图。
在START,处理过程630预置峰值指数i=1并确定最高峰值,peak(k)=max{(K)},其中(k)是由运行在第一自适应模式的自适应滤波器提供的信道权向量(方框810)。之后,处理过程630确定位于peak(k)周围的窗口外部的次峰值,next_peak=max{W(k)-R},其中R是围绕peak(i)的区域,最近定位的峰值(方框820)。在一个实施例中,R覆盖任一个peak(i)上的25个抽样,总共51抽样。然后,处理过程630确定next_peak是否大于β*最高峰值(方框830)。在一个实施例中,β等于0.6。如果next_peak不大于β*最高峰值,next_peak假定为一个假峰值然后就被删除(方框840)。否则,处理过程630确定peak(i)是next-peak(方框850)。然后,处理过程630确定迄今为止定位的峰值数目是否等于L,其中L是待定位的峰值的总额(方框860)。在一个实施例中,L等于5。否则,处理过程630增加峰值指数i=i+1(方框870),然后返回方框820。否则,处理过程630确定来自peak(i)的脉冲响应的上升沿,其中i=1,...,L(方框880)。脉冲响应的上升沿位于L峰值中最小峰值上。该上升沿提供估计的容积延迟。处理过程630则终止。
使用AP修正的该回波消除对于有色输入的收敛速度比NLMS方法快。基于AP的回波消除可以多种方式实现。在一个实施例中,自适应滤波器的长度M(也就是,信道权向量的大小)在第一和第二自适应模式里是相同的。M的三个值用于比较:M=256,512和1024。在另一个实施例里,在第一自适应模式中的自适应滤波器的长度M1比在第二自适应模式里的长度M2长。M1和M2的典型值分别是1024和256。还有两种方式启动第二自适应模式。在第一方法中,在Rin的第一期间估计容积延迟。容积延迟估计用来对准M=256抽头滤波器,因此回波信道脉冲响应属于短滤波器。256抽头滤波器在用于估计延迟的通道脉冲响应的估计上预置,并且相同的Rin数据重新运算。在第二方法里,自适应起始于延迟估计结束的时间。换句话说,它遵循与在第一方法中相同的程序,但是不重新启动该数据流。
用于回波消除的AP(2)方法的结果在图表I和II,以及图9A、9B、10A、10B、11A、11B、12A、12B、13中显示。标准AP(2)方法利用用于G.I68 test 2b的数据文件测试(Rin和Sin)。
图表I显示当Rin是非零时,用于AP(2)和NLMS技术的分贝上的衰减。AP(2)方法在第一和第二自适应模式里使用相同的滤波长度,M=256、512和1024。可见AP(2)方法在自适应的不同步骤上具有近乎10分贝的性能改善。
M=256  M=512  M=1024
迭代数目 AP(2) NLMS  AP(2) NLMS  AP(2) NLMS
 1000 -33.49 -26.31 -24.31 -22.58 -18.44 -18.4
 1500 -40.92 -30.22 -29.62 -25.51 -23.87 -23.18
 2000 -44.86 -36.61 -35.40 -29.23 -25.52 -23.86
 4800 -45.14 -44.57 -45.79 -36.92 -34.79 -27.89
 7600 -45.65 -47.07 -47.13 -42.47 -41.06 -31.04
 10400 -47.53 -49.96 -46.13 -45.54 -44.91 -33.24
 13200 -46.0 -49.06 -46.67 -47.80 -46.07 -35.38
 16000 -45.07 -47.89 -48.73 -49.99 -47.23 -37.52
 18800 -47.49 -48.75 -46.10 -48.67 -45.75 -38.57
 21600 -48.09 -50.47 -47.61 -50.27 -47.98 -40.66
 24400 -45.80 -48.22 -46.63 -48.7 -47.47 -41.78
表I
图表II显示在第二自适应模式的第一和第二方式中对于AP(2)的分贝上的衰减。在第一方法中延迟估计所需要的时间被增加给迭数字用于一个精确的对比。估计发生在大约1500次迭代并且使用Rin的1500个抽样。当滤波器充满数据时,这相当于AP(2)技术的500次迭代。因此,修正值由500以及1500次迭代限制。然而,上界迭代数落于用于没有取消特性结果的Rin的死区。对于第一方法,下界产生一个显著的优点。对于第一方法,收敛发生在2500以及3500迭代之间的某处,并且对于第二方法,收敛发生在小于4800迭代的某处。
    迭代     第一种方式     第二种方式
    1000     -38.48     -18.44
    1250     -41.91     -23.01
    1500     -43.03     -23.87
    1750     -47.05     -29.24
    2000     -45.99     -33.9
    3800     -     -37.53
    4300 -     -42.77
    4800 -     -44.89
表II
图9A是按照本发明的一个实施例,说明使用M=256的AP(2)的回波消除器中衰减的分贝的框图。
图9B是按照本发明的一个实施例,说明使用M=256的NLMS的回波消除器中衰减的分贝的框图。
图10A是按照本发明的一个实施例,说明使用M=512的AP(2)的回波消除器中衰减的分贝的框图。
图10B是按照本发明的一个实施例,说明使用M=512的NLMS的回波消除器中衰减的分贝的框图。
图11A是按照本发明的一个实施例,说明使用M=1024的AP(2)的回波消除器中衰减的分贝的框图。
图11B是按照本发明的一个实施例,说明使用M=1024的NLMS的回波消除器中衰减的分贝的框图。
图12A是按照本发明的一个实施例,说明在第二自适应模式里使用M1=1024和M2=256的AP(2)的第一种方法的回波消除器中的分贝衰减的框图。
图12B是按照本发明的一个实施例,说明在第二自适应模式里使用M1=1024和M2=256的AP(2)的第二种方法的回波消除器中的分贝衰减的框图。
图13是说明按照本发明一个实施例,M+500迭代上的加权的框图。可以清楚地看到回波信道脉冲响应的三个峰值。
尽管本发明已经参考说明性的实施例进行了描述,该描述不用于解释为限制的含义。说明性实施例的不同改进以及该发明其它实施例,对于所属技术领域的专业人员来说是显而易见的,被认为是归于本发明实质和范围之内。

Claims (40)

1.一个装置包括:
一个自适应滤波器,使用仿射投影(AP)修正估计回波信道的信道权向量,该回波信道接收一个发送输入序列和一个接收输入序列,当自适应滤波器分别运行在一个第一自适应模式和一个第二自适应模式时,该信道权向量具有第一和第二长度;和
一个连接到自适应滤波器的延迟估计器,使用在第一自适应模式里的自适应滤波器确定回波信道中的延迟。
2.权利要求1的装置,其中该自适应滤波器基于由延迟估计器确定的延迟运行在第二自适应模式中。
3.权利要求1的装置,其中第一长度比第二长度长。
4.权利要求1的装置,其中接收输入序列由随机序列和由一个自回归(AR)系数向量加权的过去接收输入序列之和表示。
5.权利要求4的装置,其中该自适应滤波器包括:
一个AR系数估计器,在第一修正速度上使用过去接收输入序列和接收输入序列估计AR系数向量;
一个连接到AR系数估计器的随机估计器,在第二修正速度上使用估计的AR系数向量估计该随机序列;
一个误差估计器,在第二修正速度上使用发送输入序列、接收输入序列以及估计的信道权向量估计误差;和
一个连接到该误差估计器和随机估计器的加权修正器,在第二修正速度上使用估计误差和估计的随机序列修正该信道权向量。
6.权利要求1的装置,其中延迟估计器包括:
一个峰值定位器,使用估计的信道权向量来定位回波信道脉冲响应的L峰值,L峰值位于彼此一预定距离以内。
7.权利要求6的装置,其中延迟估计器进一步地包括:
一个峰值消除器,在L峰值中的一个里删除一个假峰值;和
一个上升沿定位器,用于定位脉冲响应的一个上升沿,该上升沿是在L峰值里具有最短延迟的峰值并且距离该延迟一预定距离以内。
8.权利要求5的装置,其中当运行在第二自适应模式时,该自适应滤波器通过在第一自适应模式中修正的估计的信道权向量来预置信道权向量。
9.权利要求5的装置,其中第一修正速度比第二修正速度慢。
10.权利要求5的装置,其中当短期平均功率小于长期平均功率一个阈值功率数量时,加权修正器冻结信道权向量。
11、一个方法包括:
通过使用仿射投影(AP)修正的自适应滤波器估计回波信道的信道权向量,该回波信道接收一个发送输入序列和一个接收输入序列,当自适应滤波器运行在第一自适应模式时,信道权向量具有一个第一长度,当自适应滤波器运行在一个第二自适应模式时,信道权向量具有一个第二长度;以及
通过延迟估计器使用在第一自适应模式里的自适应滤波器来确定回波信道中的延迟。
12.权利要求11的方法,其中估计信道权向量包括基于由延迟估计器确定的延迟来在第二自适应模式运行。
13.权利要求11的方法,其中第一长度比第二长度长。
14.权利要求11的方法,其中接收输入序列由随机序列和由一个自回归(AR)系数向量加权的过去接收输入序列之和表示。
15.权利要求14的方法,其中估计信道权向量包括:
在第一修正速度上使用过去接收输入序列和接收输入序列估计AR系数向量;
在第二修正速度上使用估计AR系数向量估计该随机序列;
在第二修正速度上使用发送输入序列、接收输入序列和估计的信道权向量估计误差;以及
在第二修正速度上使用估计误差和估计的随机序列修正信道权向量。
16.权利要求11的方法,其中确定该延迟包括:
使用估计的信道权向量利用峰值定位器来定位回波信道脉冲响应的L峰值,L峰值位于彼此一预定距离以内。
17.权利要求16的方法其中确定该延迟进一步地包括:
在L峰值中的一个里删除一个假峰值;和
定位脉冲响应的一个上升沿,该上升沿是在L峰值里具有最短延迟的峰值并且距离该延迟一预定距离以内。
18.权利要求15的方法,其中运行在第二自适应模式包括通过在第一自适应模式中修正的估计的信道权向量来预置信道权向量。
19.权利要求15的方法,其中第一修正速度比第二修正速度慢。
20.权利要求15的方法,其中修正信道权向量包括当短期平均功率小于长期平均功率一个阈值功率数量时,冻结信道权向量。
21.一个系统包括:
一个连接到声道远端的第一解码器,解码远端信号,第一解码器产生一个接收输入序列;
一个连接到声道近端的第二解码器,解码一个近端信号,第二解码器产生一个发送输入序列;以及
连接到第一和第二解码器的回波信道中的一个回波消除器,执行回波消除,该回波信道接收接收和发送输入序列,该回波消除器包括:
一个自适应滤波器,使用仿射投影(AP)修正估计回波信道的信道权向量,当自适应滤波器分别运行在一个第一自适应模式和一个第二自适应模式时,该信道权向量具有第一和第二长度,以及
一个连接到该自适应滤波器的延迟估计器,使用在第一自适应模式里的自适应滤波器确定一个延迟。
22、权利要求21的系统,其中自适应滤波器基于由延迟估计器确定的延迟运行在第二自适应模式中。
23.权利要求21的系统,其中第一长度比第二长度长。
24.权利要求21的系统,其中接收输入序列由随机序列和由一个自回归(AR)系数向量加权的过去接收输入序列之和表示。
25.权利要求24的系统,其中该自适应滤波器包括:
一个AR系数估计器,在第一修正速度上使用过去接收输入序列和接收输入序列估计AR系数向量;
一个连接到该AR系数估计器的随机估计器,在第二修正速度上使用估计AR系数向量估计该随机序列;
一个误差估计器,在第二修正速度上使用发送输入序列、接收输入序列和估计的信道权向量估计误差;以及
一个连接到该误差估计器和随机估计器的加权修正器,在第二修正速度上使用估计误差和估计的随机序列修正该信道权向量。
26.权利要求21的系统,其中延迟估计器包括:
一个峰值定位器,使用估计的信道权向量来确定回波信道脉冲响应的L峰值,L峰值位于彼此距离一个预定距离以内。
27.权利要求26的系统,其中延迟估计器进一步包括:
一个峰值消除器,在L峰值中的一个里删除一个假峰值;
一个上升沿定位器,定位脉冲响应的一个上升沿,该上升沿是在L峰值里具有最短延迟的峰值并且距离该延迟一预定距离以内。
28.权利要求25的系统,其中当运行在第二自适应模式时,该自适应滤波器通过在第一自适应模式中修正的估计的信道权向量来预置信道权向量。
29.权利要求25的系统,其中第一修正速度比第二修正速度慢。
30.权利要求25的系统,其中当短期平均功率小于长期平均功率一个阈值功率数量时,加权修正器冻结信道权向量。
31.一种计算机程序产品包括:
一个在其中嵌入程序代码的机器可读的媒体,该程序代码包括:
计算机可读程序代码,通过使用一个仿射投影(AP)修正的自适应滤波器估计回波信道的一个信道权向量,该回波信道接收一个发送输入序列和一个接收输入序列,当自适应滤波器运行在第一自适应模式时,信道权向量具有一个第一长度,当自适应滤波器运行在一个第二自适应模式时,信道权向量具有一个第二长度;以及
计算机可读程序代码,通过延迟估计器使用在第一自适应模式里的自适应滤波器来确定回波信道中的延迟。
32.权利要求31的计算机程序产品,其中估计信道权向量的计算机可读程序代码包括基于由延迟估计器确定的延迟运行在第二自适应模式中的计算机可读程序代码。
33.权利要求31的计算机程序产品,其中第一长度比第二长度长。
34.权利要求31的计算机程序产品,其中接收输入序列由随机序列和由一个自回归(AR)系数向量加权的过去接收输入序列之和表示。
35.权利要求34的计算机程序产品,其中估计信道权向量的计算机可读程序代码包括:
在第一修正速度上使用过去接收输入序列和接收输入序列估计AR系数向量的计算机可读程序代码;
在第二修正速度上使用估计的AR系数向量估计随机序列的计算机可读程序代码;
在第二修正速度上使用发送输入序列、接收输入序列和估计的信道权向量估计误差的计算机可读程序代码;和
在第二修正速度上使用估计误差和估计的随机序列修正信道权向量的计算机可读程序代码。
36.权利要求31的计算机程序产品,其中确定延迟的计算机可读程序代码包括:
使用估计的信道权向量由峰值定位器来定位回波信道脉冲响应的L峰值的计算机可读程序代码,该L峰值位于彼此一个预定距离以内。
37.权利要求36的计算机程序产品,其中确定延迟的计算机可读程序代码进一步地包括:
在L峰值中的一个里删除一个假峰值的计算机可读程序代码;和
定位脉冲响应的上升沿的计算机可读程序代码,该上升沿是L峰值中具有最短延迟的峰值并且位于距离该延迟一预定距离以内。
38.权利要求35的计算机程序产品,其中运行在第二自适应模式的计算机可读程序代码包括通过在第一自适应模式中修正的估计的信道权向量来预置信道权向量的计算机可读程序代码。
39.权利要求35的计算机程序产品,其中第一修正速度比第二修正速度慢。
40.权利要求35的计算机程序产品,其中修正的计算机可读程序代码包括:当短期平均功率小于长期平均功率一阈值功率数量时,冻结该信道权向量的计算机可读程序代码。
CN018185274A 2000-09-08 2001-09-07 用于基于回波消除器的快速收敛仿射投影的装置和方法 Expired - Fee Related CN1473403B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US23142000P 2000-09-08 2000-09-08
US60/231,420 2000-09-08
US09/947,804 2001-09-06
US09/947,804 US6700978B2 (en) 2000-09-08 2001-09-06 Method and apparatus for fast converging affine projection based echo canceller
PCT/US2001/042088 WO2002021717A2 (en) 2000-09-08 2001-09-07 Digital echo cancellor

Publications (2)

Publication Number Publication Date
CN1473403A true CN1473403A (zh) 2004-02-04
CN1473403B CN1473403B (zh) 2010-05-05

Family

ID=26925112

Family Applications (1)

Application Number Title Priority Date Filing Date
CN018185274A Expired - Fee Related CN1473403B (zh) 2000-09-08 2001-09-07 用于基于回波消除器的快速收敛仿射投影的装置和方法

Country Status (9)

Country Link
US (1) US6700978B2 (zh)
EP (1) EP1320941B1 (zh)
CN (1) CN1473403B (zh)
AT (1) ATE332040T1 (zh)
AU (1) AU2001293264A1 (zh)
CA (1) CA2421759C (zh)
DE (1) DE60121186D1 (zh)
HK (1) HK1054136A1 (zh)
WO (1) WO2002021717A2 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100394579C (zh) * 2005-03-03 2008-06-11 中国科学院微电子研究所 超大规模集成电路设计中保持时间快速收敛的方法
CN103107838A (zh) * 2013-01-17 2013-05-15 杭州电子科技大学 用于上行多用户mimo系统中基于用户选择的仿射投影解码法
CN103152137A (zh) * 2013-01-17 2013-06-12 杭州电子科技大学 上行多用户多输入多输出系统中的仿射投影解码法
CN104392727A (zh) * 2014-11-12 2015-03-04 华为技术有限公司 音频信号处理方法和相关装置
CN109089004A (zh) * 2018-08-14 2018-12-25 西南交通大学 一种基于相关熵诱导的集员自适应回声消除方法

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6792106B1 (en) * 1999-09-17 2004-09-14 Agere Systems Inc. Echo canceller and method of echo cancellation using an NLMS algorithm
DE10012868C1 (de) * 2000-03-16 2001-07-26 Siemens Ag Rechnergerät mit Mehrwege-Akustikaus- und -eingabesystem
US7627111B2 (en) * 2002-11-25 2009-12-01 Intel Corporation Noise matching for echo cancellers
CN1859519B (zh) * 2005-11-19 2010-09-29 华为技术有限公司 一种自适应滤波器及回波抵消器
JP4432916B2 (ja) * 2006-03-03 2010-03-17 ソニー株式会社 音声処理装置
CN101536342B (zh) 2006-11-15 2012-10-10 西门子公司 用于信号的自适应滤波的方法和装置
US8254561B1 (en) * 2007-04-17 2012-08-28 Plantronics, Inc. Headset adapter with host phone detection and characterization
KR101975251B1 (ko) 2012-07-09 2019-05-07 삼성전자주식회사 오디오 신호 처리 시스템 및 이의 에코 신호 제거 방법
CN104144269B (zh) * 2014-08-08 2016-03-02 西南交通大学 一种基于解相关的成比例自适应电话回声消除方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04123621A (ja) 1990-09-14 1992-04-23 Nippon Telegr & Teleph Corp <Ntt> 反響消去装置
FI935834A (fi) * 1993-12-23 1995-06-24 Nokia Telecommunications Oy Menetelmä kaikukohtaan sovittautumiseksi kaiunpoistajassa
JP3654470B2 (ja) 1996-09-13 2005-06-02 日本電信電話株式会社 サブバンド多チャネル音声通信会議用反響消去方法
US6137881A (en) * 1997-02-28 2000-10-24 Texas Instruments Incorporated Adaptive filtering method and apparatus employing modified fast affine projection algorithm
JP3236242B2 (ja) 1997-05-12 2001-12-10 沖電気工業株式会社 エコーキャンセラ装置

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100394579C (zh) * 2005-03-03 2008-06-11 中国科学院微电子研究所 超大规模集成电路设计中保持时间快速收敛的方法
CN103107838A (zh) * 2013-01-17 2013-05-15 杭州电子科技大学 用于上行多用户mimo系统中基于用户选择的仿射投影解码法
CN103152137A (zh) * 2013-01-17 2013-06-12 杭州电子科技大学 上行多用户多输入多输出系统中的仿射投影解码法
CN104392727A (zh) * 2014-11-12 2015-03-04 华为技术有限公司 音频信号处理方法和相关装置
CN109089004A (zh) * 2018-08-14 2018-12-25 西南交通大学 一种基于相关熵诱导的集员自适应回声消除方法
CN109089004B (zh) * 2018-08-14 2020-10-09 西南交通大学 一种基于相关熵诱导的集员自适应回声消除方法

Also Published As

Publication number Publication date
WO2002021717A2 (en) 2002-03-14
CA2421759C (en) 2008-09-02
CN1473403B (zh) 2010-05-05
EP1320941A2 (en) 2003-06-25
ATE332040T1 (de) 2006-07-15
EP1320941B1 (en) 2006-06-28
US6700978B2 (en) 2004-03-02
DE60121186D1 (de) 2006-08-10
US20020071547A1 (en) 2002-06-13
AU2001293264A1 (en) 2002-03-22
HK1054136A1 (en) 2003-11-14
CA2421759A1 (en) 2002-03-14
WO2002021717A3 (en) 2003-02-27

Similar Documents

Publication Publication Date Title
CN1473403A (zh) 用于基于回波消除器的快速收敛仿射投影的装置和方法
CN1075692C (zh) 通信系统中噪声抑制方法及装置
CN1097360C (zh) 数字移动通信系统中的回声消除
CN1218506C (zh) 采用状态判断来控制数字电话系统中的功能元件的方法和装置
CN1143265C (zh) 具有改进的语音编码器的传输系统
CN1041374C (zh) 通信网回波抵消器
CN1223108C (zh) 回声消除方法和回声消除器
CN1129263C (zh) 分组编码信号的方法和装置
CN1118956C (zh) 回声路径延迟估计
CN1111973C (zh) 改进了的数字蜂窝应用的回波消除器
CN1114122A (zh) 用于回声抑制器的话音作用检测器和回声抑制器
CN1671161A (zh) 回波消除器电路和方法
CN1269926A (zh) 用于回波估算和抑制的方法和设备
CN1706174A (zh) 综合噪声消除及残余回波抑制
CN1529882A (zh) 用于扩展窄带滤波的语音信号、特别是由通信设备发送的语音信号的带宽的方法
CN102855881B (zh) 一种回声抑制方法和装置
KR100486732B1 (ko) 블럭제한된 트렐리스 부호화 양자화방법과 음성부호화시스템에있어서 이를 채용한 라인스펙트럼주파수 계수양자화방법 및 장치
CN1505894A (zh) 移动站内的网络回波抑制
US20020093919A1 (en) Fast converging affine projection based echo canceller for sparse multi-path channels
CN1350727A (zh) 纯延迟估计
CN101163178B (zh) 一种无线通信系统中回声抵消方法及其系统
JP2003514264A (ja) 雑音抑圧装置
CN100346577C (zh) 信号编码装置和信号解码装置,以及信号编码方法和信号解码方法
JP2005506581A (ja) 正弦波モデルパラメータの周波数差分符号化
CN1503467A (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
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: 20100505

Termination date: 20170907