CN109961798A - 回声消除系统、方法、可读计算机存储介质、及终端 - Google Patents
回声消除系统、方法、可读计算机存储介质、及终端 Download PDFInfo
- Publication number
- CN109961798A CN109961798A CN201711437294.4A CN201711437294A CN109961798A CN 109961798 A CN109961798 A CN 109961798A CN 201711437294 A CN201711437294 A CN 201711437294A CN 109961798 A CN109961798 A CN 109961798A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- time
- echo
- peak
- varying
- 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
Links
- 238000000034 method Methods 0.000 title claims description 26
- 230000005540 biological transmission Effects 0.000 claims abstract description 18
- 239000000872 buffer Substances 0.000 claims description 45
- 238000001228 spectrum Methods 0.000 claims description 31
- 230000003044 adaptive effect Effects 0.000 claims description 27
- 238000004422 calculation algorithm Methods 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000005055 memory storage Effects 0.000 claims description 3
- 230000001934 delay Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 7
- 230000025925 convergent extension Effects 0.000 abstract description 2
- 238000001514 detection method Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
-
- 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
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L2021/02082—Noise filtering the noise being echo, reverberation of the speech
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明提供一种回声消除系统,包括:自适应滤波器,用于跟踪回声传输信道;辅助滤波器,用于与所述自适应滤波器一同消除回声;回声消除控制器,用于控制所述自适应滤波器和辅助滤波器,从而实现消除回声。本发明的回声消除系统可有效检测是否存在近端信号;在自适应滤波器收敛较好时利用自适应滤波器消除回声并保存其时变系数;在自适应滤波器发散时,可利用收敛的辅助滤波器消除回声;自适应滤波器作为主滤波器能够及时跟踪回声传输信道的突变情况,故本发明的回声消除系统的适应性非常强,从而实现了存在近端信号时依然能够有效地消除回声的技术效果。
Description
技术领域
本发明涉及语音通信领域,特别是涉及回声消除系统、方法、可读计算机存储介质、及终端。
背景技术
在现代多媒体通信系统中,语音业务占有着举足轻重的位置。无论在传统的固定电话业务还是在VOIP网络电话业务或者各种手持设备和免提电话系统中,为提高通信质量,必须在通信设备中集成回声消除器以消除或抑制回声。
一般而言,回声消除算法通过自适应滤波来完成,其基本原理如图1所示。其中,远端声音x(k)通过回声传输信道h产生回声y(k),近端信号d(k)是由回声y(k)和近端声音v(k)(可包含噪声信号)混合而成。通过使用自适应滤波器来模拟回声传输信道h,可以使所得逼近回声信号,进而达到回声消除的目的。由此可见,回声消除的关键是自适应地调整所述自适应滤波器使其逼近回声传输信道h,以实现回声的消除。
但是,上述方法只在不存在近端声音v(k)的情况下才能得到满意的回声消除效果,当存在明显的近端声音v(k)时,自适应滤波器的性能将会恶化,甚至不能保证算法的收敛性。因此,如何有效地检测出是否存在近端声音,成为回声消除算法在实际应用中必须解决的关键问题。
检测是否存在近端声音也被称为双端发声检测,是回声消除的技术瓶颈。现有双端发声检测算法大致可以分为:能量比较法,相关计算法及其他方法。但是,这两种双端发声检测方法都会有一定的误检概率。若在误检之后停止自适应滤波器的迭代,则可能会使检测和自适应滤波陷入不能正常工作的状态,而无法恢复。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供回声消除系统和方法,用于解决现有技术中双端发声情况下因近端声的存在而导致自适应滤波器无法收敛回声无法消除的问题。
为实现上述目的及其他相关目的,本发明提供一种回声消除系统,包括:自适应滤波器,用于跟踪回声传输信道;辅助滤波器,用于与所述自适应滤波器一同消除回声;回声消除控制器,用于分析自适应滤波器输出信号与近端信号之间的关联程度,据以计算所述自适应滤波器的时变系数;将所述时变系数与预设阈值范围做比较,以判断所述时变系数是否满足收敛条件;若满足收敛条件,则通过所述自适应滤波器消除回声;若不满足收敛条件,则将缓存器中前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声。
于本发明的一实施例中,所述分析自适应滤波器输出信号与近端信号之间的关联程度,具体包括:判断k是否为M的整数倍;其中,k为时间,M为预设向量长度;若k不为M的整数倍,则令k=k+1直至k为M的倍数;若k为M的整数倍,则基于自适应算法对近端信号d(k)以及自适应滤波器输出信号进行2M次自适应计算,以分别得到第一信号d′(k)和第二信号将所述第一信号d′(k)和第二信号分别依次进行频域变换、功率谱迭代值计算、以及互功率谱迭代值计算,据以得到关联系数其中,所述关联系数代表所述自适应滤波器输出信号与近端信号d(k)之间的关联程度。
于本发明的一实施例中,所述自适应算法包括NLMS算法。
于本发明的一实施例中,基于自适应算法对所述近端信号d(k)以及自适应滤波器输出信号进行2M次自适应计算,以分别得到第一信号d′(k)和第二信号具体包括:d′(k)=[d(k-2M),d(k-2M+1),...,d(k-M),...,d(k-1),d(k)]T;2M,yk-2M+1,...,yk-M,...,yk-1,ykT。
于本发明的一实施例中,所述将所述第一信号d′(k)和第二信号分别依次进行频域变换、功率谱迭代值计算、以及互功率谱迭代值计算,据以得到关联系数具体包括:将所述第一信号d′(k)和第二信号进行频域变换,以分别得到对应的第一频域信号Dk和第二频域信号Yk;其中,Dk=FFT[d′(k)]T,将所述第一频域信号Dk和第二频域信号Yk进行功率谱计算,以分别得到对应的近端信号d(k)的功率谱迭代值以及自适应滤波器输出信号的功率谱迭代值 其中,U为控制迭代估计时间常数的常量,D* k为所述第一频域信号Dk的共轭值,Y* k为所述第二频域信号Yk的共轭值;计算互功率谱迭代值 根据所述功率谱迭代值功率谱迭代值Syk,以及互功率谱迭代值计算所述关联系数其中,为所述互功率谱迭代值的共轭值。
于本发明的一实施例中,所述根据所述输出信号与近端信号的关联程度计算所述自适应滤波器的时变系数,具体包括:根据所述关联系数计算所述自适应滤波器的时变系数r,
于本发明的一实施例中,根据所述时变系数分析所述自适应滤波器与回声传输信道的逼近程度,据以判断所述时变系数是否满足收敛条件;若满足收敛条件,则将所述时变系数存入缓存器中以形成缓存系数;若不满足收敛条件,则将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声,具体包括:分别比较所述时变系数与第一阈值和第二阈值之间的大小关系;其中,所述第一阈值小于第二阈值;若所述时变系数大于所述第二阈值,则所述时变系数收敛到较好程度,回声消除误差并且:判断所述时变系数与所述缓存器前一次保存的缓存系数的收敛程度;若所述时变系数的收敛程度高于所述缓存器前一次保存的缓存系数的收敛程度,则将所述时变系数存入所述缓存器中以作为此次保存的缓存系数;若所述时变系数小于所述第一阈值,则变所述时系数发散;将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声,所述回声消除误差其中,为所述辅助滤波器的输出信号;若所述时变系数大于等于第一阈值,且小于等于第二阈值,则判定为无所述近端信号,所述回声消除误差
于本发明的一实施例中,所述判断所述时变系数与所述缓存器前一次保存的缓存系数的收敛程度,具体包括:比较所述时变系数r的峰值rpeak与所述缓存器前一次保存的缓存系数r′的峰值r′peak,且比较所述时变系数r的峰值位置kpeak与缓存系数r′的峰值位置k′peak;若rpeak>r′peak或者kpeak<k′peak+to,则所述时变系数r的收敛程度甚于所述缓存系数,令r′peak=rpeak,k'peak=kpeak;其中,to表示所述时变系数r最多维持峰值的有效时间。
为实现上述目的及其他相关目的,本发明提供一种回声消除方法,应用于回声消除控制器,其包括:分析自适应滤波器输出信号与近端信号之间的关联程度;根据所述输出信号与近端信号之间的关联程度计算所述自适应滤波器的时变系数;根据所述时变系数分析所述自适应滤波器与回声传输信道的逼近程度,据以判断所述时变系数是否满足收敛条件;若满足收敛条件,则将所述时变系数存入缓存器中以形成缓存系数;若不满足收敛条件,则将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声。
于本发明的一实施例中,根据所述时变系数分析所述自适应滤波器与回声传输信道的逼近程度,据以判断所述时变系数是否满足收敛条件;若满足收敛条件,则将所述时变系数存入缓存器中以形成缓存系数;若不满足收敛条件,则将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声,具体包括:分别比较所述时变系数与第一阈值和第二阈值之间的大小关系;其中,所述第一阈值小于第二阈值;若所述时变系数大于所述第二阈值,则所述时变系数收敛到较好程度,回声消除误差并且:判断所述时变系数与所述缓存器前一次保存的缓存系数的收敛程度;若所述时变系数的收敛程度高于所述缓存器前一次保存的缓存系数的收敛程度,则将所述时变系数存入所述缓存器中以作为此次保存的缓存系数;若所述时变系数小于所述第一阈值,则变所述时系数发散;将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声,所述回声消除误差其中,为所述辅助滤波器的输出信号;若所述时变系数大于等于第一阈值,且小于等于第二阈值,则判定为无所述近端信号,所述回声消除误差
于本发明的一实施例中,所述判断所述时变系数与所述缓存器前一次保存的缓存系数的收敛程度,具体包括:比较所述时变系数r的峰值rpeak与所述缓存器前一次保存的缓存系数r′的峰值r′peak,且比较所述时变系数r的峰值位置kpeak与缓存系数r′的峰值位置k′peak;若rpeak>r'peak或者kpeak<k'peak+to,则所述时变系数r的收敛程度甚于所述缓存系数,令r'peak=rpeak,k′peak=kpeak;其中,to表示所述时变系数r最多维持峰值的有效时间。
为实现上述目的及其他相关目的,本发明提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现所述回声消除方法。
为实现上述目的及其他相关目的,本发明提供一种回声消除终端,包括:处理器及存储器;所述存储器用于存储计算机程序,所述处理器用于执行所述存储器存储的计算机程序,以使所述终端执行所述回声消除方法。
如上所述,本发明的回声消除系统、方法、计算机可读存储介质、以及终端,具有以下有益效果:本发明的回声消除系统可有效检测是否存在近端信号;在自适应滤波器收敛较好时利用自适应滤波器消除回声并保存其时变系数;在自适应滤波器发散时,可利用收敛的辅助滤波器消除回声;自适应滤波器作为主滤波器能够及时跟踪回声传输信道的突变情况,故本发明的回声消除系统的适应性非常强,从而实现了存在近端信号时依然能够有效地消除回声的技术效果。
附图说明
图1显示为现有技术中回声消除系统的示意图。
图2显示为本发明一实施例中回声消除系统的示意图。
图3显示为本发明一实施例中回声消除系统具体实施步骤的流程图。
图4显示为本发明一实施例中回声消除方法的流程图。
元件标号说明
21 自适应滤波器
22 辅助滤波器
23 回声消除控制器
24 回声传输信道
25 缓存器
S301~S316 步骤
S401~S406 步骤
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
本发明提供回声消除系统及方法,用于在双端发声的情况下消除回声。
如图2所示,展示本发明一实施例中的回声消除系统,其包括:自适应滤波器21、辅助滤波器22、以及回声消除控制器23。其中,所述自适应滤波器用于跟踪回声传输信道24;所述辅助滤波器22,用于辅助所述自适应滤波器21一同消除回声;所述回声消除控制器23,用于分析自适应滤波器21输出信号与近端信号之间的关联程度,根据所述输出信号与近端信号之间的关联程度计算所述自适应滤波器21的时变系数,并根据所述时变系数分析所述自适应滤波器21与回声传输信道24的逼近程度,据以判断所述时变系数是否满足收敛条件;若满足收敛条件,则将所述时变系数存入缓存器25中以形成缓存系数;若不满足收敛条件,则将所述缓存器25前一次保存的缓存系数载入辅助滤波器22中,以通过所述辅助滤波器22消除回声。
所述自适应滤波器21,是一种指根据环境的改变,使用自适应算法改变其参数和结构的滤波器;具体的,所述自适应滤波器21设有时变系数,所述时变系数由自适应算法更新,可自动连续地适应于给定信号,以获得期望响应。所述自适应滤波器21例如可选用FIR(Finite Impulse Response)非递归型滤波器,也可选用IIR(Infinite ImpulseResponse)递归型滤波器等等,下面以以FIR自适应滤波器以及FIR辅助滤波器为例,说明所述回声消除系统的工作原理。
令所述FIR自适应滤波器的时变系数为r,则用rpeak代表所述时变系数r的峰值,用kprak代表峰值的位置,且用to代表峰值最多维持的有效时间。于该实施例中,设定所述时变系数r的第一阈值T1=0.2,第二阈值T2=0.8,所述时变系数r与T1、T2的大小关系即代表该FIR自适应滤波器收敛或发散的程度。需要说明的是,T1取0.2,T2取0.8只是本发明的其中一个实施例,而非限定T1和T2的取值。
如图3所示,展示本发明一实施例中所述回声消除系统具体的执行步骤,实现消除回声。
S301:所述FIR自适应滤波器和FIR辅助滤波器接收输入信号x(k),并分别计算N阶滤波的结果。其中,N阶滤波意为滤波器过滤谐波的次数;令所述FIR自适应滤波器的输出信号为FIR辅助滤波器的输出信号为k则代表时间。回声消除误差信号e(k)由回声消除控制器23控制,可选接所述FIR自适应滤波器或者FIR辅助滤波器。也即,回声消除误差信号或者图3中,用虚线代表所述回声消除控制器控制所述回声误差信号。
需要说明的是,为保证所述FIR自适应滤波器能够稳定且快速地跟踪所述回声传输信道,优选的,时变系数的更新不受所述回声消除控制器23的控制而持续进行更新。具体的,所述FIR自适应滤波器可基于NLMS自适应算法更新时变系数;所述NLMS自适应算法,全称为Normalized Least Mean Square,意为归一化最小均方算法;因NLMS自适应算法为现有技术,故不再赘述。
S302:判断k是否为M的整数倍。其中,M意为预设向量长度,于本实施例中M取值64。需要说明的是,M可取64但并不限定为64,在其他的实施例中也可取其他值作为预设向量长度,本发明对此不做限定。
S303:若k为M的整数倍,则基于NLMS自适应算法对近端信号d(k)以及自适应滤波器21输出信号进行2M次自适应计算,以分别得到第一信号d′(k)和第二信号其中,2M次自适应计算是指由M次自适应计算外加上一次迭代的M次自适应计算,共计2M次计算所得。
本发明选用2M次自适应算法,相比于现有技术的M次算法,不仅大大减少了算法的计算量还因自适应次数的加倍而大幅提升了计算精度。简单来说,M次自适应算法的基本原理是把一个M长的向量分解为两个M/2长的向量,再把M/2长向量分解为M/4长向量,再依次分解直至分解到长度为2的向量位置。在选取2M次算法的情况下,以M次历史信号作为输入,每迭代M次即计算一次相关值,不仅可大大减少计算复杂度,使傅里叶变换的计算量大为减少,还可大幅提升计算所得的相关值的精确度,减少误判的可能性。具体而言:
d′(k)=[d(k-2M),d(k-2M+1),...,d(k-M),...,d(k-1),d(k)]T;
S304:将所述第一信号d′(k)和第二信号进行频域变换。所谓频域变换,是指将所述第一信号d′(k)和第二信号进行FFT快速傅氏变换,以分别得到第一频域信号Dk和第二频域信号Yk。具体的,Dk=FFT[d′(k)]T,
S305:将所述第一频域信号Dk和第二频域信号Yk进行功率谱计算,以分别得到和所述为近端信号d(k)的功率谱迭代值,Syk为自适应滤波器21输出信号的功率谱迭代值。具体的,Syk=(1-U)Syk-1+UYkY* k。
S306:计算互功率谱迭代值所述互功率谱迭代值用于表示所述第一频域信号Dk和第二频域信号Yk之间的关联性。具体的,其中,U为控制迭代估计时间常数的常量,例如可取值0.9,D* k为所述第一频域信号Dk的共轭值,Y* k为所述第二频域信号Yk的共轭值。
S307:根据所述功率谱迭代值功率谱迭代值Syk,以及互功率谱迭代值计算所述关联系数所述关联系数代表所述自适应滤波器21输出信号与近端信号d(k)之间的关联程度。具体的,为所述互功率谱迭代值的共轭值。
S308:根据所述关联系数计算所述FIR自适应滤波器的时变系数r,其中,所述时变系数r为(M+1)个关联系数的平均值,所述时变系数r用于与第一阈值T1和第二阈值T2做比较,以判断该时变系数的收敛或者发散程度。
S309:判断时变系数r是否满足r>T2。
S310:若满足r>T2,则判定时变系数r收敛到较好的程度,则无需FIR辅助滤波器的辅助,所述FIR自适应滤波器自身即能完成消除回声,故回声消除误差
S311:判断时变系数r的收敛程度是否甚于缓存器25前一次保存的缓存系数的收敛程度。具体而言,比较所述时变系数r的峰值rpeak与所述缓存器25前一次保存的缓存系数r′的峰值r′peak,且比较所述时变系数r的峰值位置kpeak与缓存系数r′的峰值位置k′peak;判断时变系数r是否满足:rpeak>r′peak或者kpeak<k′peak+to,to表示所述时变系数r最多维持峰值的有效时间。
S312:若满足rpeak>r′peak或者kpeak<k′peak+to,则所述时变系数r的收敛程度甚于所述缓存系数的收敛程度,将所述时变系数r存入所述缓存器25中,以待后续时刻中出现FIR自适应滤波器发散时供FIR辅助滤波器使用。其中,将所述时变系数r存入所述缓存器25中,具体是指:令r′peak=rpeak,k′peak=kpeak。
S313:若时变系数r不满足r>T2,则判断时变系数r是否满足r<T1。
S314:若时变系数r满足r<T1,则说明时变系数r已经发散,此时所述FIR自适应滤波器已经发散,需要所述FIR辅助滤波器的辅助才能实现消除回声。具体而言,所述回声消除控制器23将缓存器25前一次保存的缓存系数载入所述FIR辅助滤波器中,且控制回声消除误差接通所述辅助滤波器22,以使所述回声消除误差
由此可知,本发明保证了在自适应滤波器发生发散的情况下能够利用保存有收敛系数的辅助滤波器进行回声的消除,解决了自适应滤波器在周围环境中存在近端信号时性能恶化而无法保证算法收敛性的技术难题。
S315:若时变系数r不满足r<T1,也即T1<r<T2,则判定回声误差系统无近端信号,所述自适应滤波器的误差
S316:若k不为M的整数倍,则令k=k+1,返回至步骤S302。
如图4所示,展示本发明提供的一种回声消除方法,应用于回声消除控制器。所述回声消除方法具体包括:
S401:分析自适应滤波器输出信号与近端信号之间的关联程度;
S402:根据所述输出信号与近端信号之间的关联程度计算所述自适应滤波器的时变系数;
S403:根据所述时变系数分析所述自适应滤波器与回声传输信道的逼近程度;
S404:判断所述时变系数是否满足收敛条件;
S405:若满足收敛条件,则将所述时变系数存入缓存器中以形成缓存系数;
S406:若不满足收敛条件,则将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声。所述回声消除方法的具体实施方式与回声消除系统的实时方式雷同,故不再赘述。
本领域普通技术人员可以理解:实现上述各方法实施例的全部或部分步骤可以通过计算机程序相关的硬件来完成。前述的计算机程序可以存储于一计算机可读存储介质中。该程序在执行时,执行包括上述各方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
本发明还提供一种回声消除终端,包括处理器和存储器,所述存储器用于存储计算机程序,所述处理器用于执行所述存储器存储的计算机程序,以使所述终端执行所述回声消除方法。上述的处理器可以是通用处理器,包括中央处理器(CentralProcessingUnit,简称CPU)、网络处理器(NetworkProcessor,简称NP)等;还可以是数字信号处理器(DigitalSignalProcessing,简称DSP)、专用集成电路(ApplicationSpecificIntegratedCircuit,简称ASIC)、现场可编程门阵列(Field-ProgrammableGateArray,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
综上所述,本发明提供的回声消除系统、方法、可读计算机存储介质、以及终端,可有效检测是否存在近端信号;在自适应滤波器收敛较好时利用自适应滤波器消除回声并保存其时变系数;在自适应滤波器发散时,可利用收敛的辅助滤波器消除回声;自适应滤波器作为主滤波器能够及时跟踪回声传输信道的突变情况,故本发明的回声消除系统的适应性非常强,从而实现了存在近端信号时依然能够有效地消除回声的技术效果。所以,本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (13)
1.一种回声消除系统,其特征在于,包括:
自适应滤波器,用于跟踪回声传输信道;
辅助滤波器,用于与所述自适应滤波器一同消除回声;
回声消除控制器,用于分析自适应滤波器输出信号与近端信号之间的关联程度,据以计算所述自适应滤波器的时变系数;将所述时变系数与预设阈值范围做比较,以判断所述时变系数是否满足收敛条件;若满足收敛条件,则通过所述自适应滤波器消除回声;若不满足收敛条件,则将缓存器中前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声。
2.根据权利要求1所述的回声消除系统,其特征在于,所述分析自适应滤波器输出信号与近端信号之间的关联程度,具体包括:
判断k是否为M的整数倍;其中,k为时间,M为预设向量长度;
若k不为M的整数倍,则令k=k+1直至k为M的倍数;
若k为M的整数倍,则基于自适应算法对近端信号d(k)以及自适应滤波器输出信号分别进行2M次自适应计算,以得到对应的第一信号d'(k)和第二信号
将所述第一信号d'(k)和第二信号分别依次进行频域变换、功率谱迭代值计算、以及互功率谱迭代值计算,据以得到关联系数其中,所述关联系数代表所述自适应滤波器输出信号与近端信号d(k)之间的关联程度。
3.根据权利要求2所述的回声消除系统,其特征在于,所述自适应算法包括NLMS算法。
4.根据权利要求2所述的回声消除系统,其特征在于,基于自适应算法对所述近端信号d(k)以及自适应滤波器输出信号分别进行2M次自适应计算,以分别得到第一信号d'(k)和第二信号具体包括:
d'(k)=[d(k-2M),d(k-2M+1),...,d(k-M),...,d(k-1),d(k)]T;
5.根据权利要求2所述的回声消除系统,其特征在于,将所述第一信号d'(k)和第二信号分别依次进行频域变换、功率谱迭代值计算、以及互功率谱迭代值计算,据以得到关联系数具体包括:
将所述第一信号d'(k)和第二信号进行FFT卷积计算,以分别得到对应的第一频域信号Dk和第二频域信号Yk;其中,Dk=FFT[d′(k)]T,
将所述第一频域信号Dk和第二频域信号Yk进行功率谱迭代值计算,以分别得到对应的近端信号d(k)的功率谱迭代值以及自适应滤波器输出信号的功率谱迭代值Syk,Syk=(1-U)Syk-1+UYkY* k;其中,U为控制迭代估计时间常数的常量,D* k为所述第一频域信号Dk的共轭值,Y* k为所述第二频域信号Yk的共轭值;
计算互功率谱迭代值
根据所述功率谱迭代值功率谱迭代值Syk,以及互功率谱迭代值计算所述关联系数 其中,为所述互功率谱迭代值的共轭值。
6.根据权利要求5所述的回声消除系统,其特征在于,所述根据所述输出信号与近端信号的关联程度计算所述自适应滤波器的时变系数,具体包括:根据所述关联系数计算所述自适应滤波器的时变系数r,
7.根据权利要求1所述的回声消除系统,其特征在于,所述将所述时变系数与预设阈值范围做比较,以判断所述时变系数的收敛程度,具体包括:
分别比较所述时变系数与第一阈值和第二阈值之间的大小关系;其中,所述第一阈值小于第二阈值;
若所述时变系数大于所述第二阈值,则所述时变系数收敛到较好程度,回声消除误差并且:
判断所述时变系数与所述缓存器前一次保存的缓存系数的收敛程度;
若所述时变系数的收敛程度高于所述缓存器前一次保存的缓存系数的收敛程度,
则将所述时变系数存入所述缓存器中替代所述前一次保存的缓存系数;
若所述时变系数小于所述第一阈值,则所述时系数发散;将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声;所述回声消除误差其中,为所述辅助滤波器的输出信号;
若所述时变系数大于等于第一阈值,且小于等于第二阈值,则判定为无所述近端信号,所述回声消除误差
8.根据权利要求7所述的回声消除系统,其特征在于,所述判断所述时变系数与所述缓存器前一次保存的缓存系数的收敛程度,具体包括:
令所述时变系数为r;比较所述时变系数r的峰值rpeak与所述缓存器前一次保存的缓存系数r′的峰值r′peak,且比较所述时变系数r的峰值位置kpeak与缓存系数r′的峰值位置k′peak;
若rpeak>r′peak或者kpeak<k′peak+to,则所述时变系数r的收敛程度甚于所述缓存系数,令r′peak=rpeak,k′peak=kpeak;其中,to表示所述时变系数r最多维持峰值的有效时间。
9.一种回声消除方法,其特征在于,应用于回声消除控制器,所述方法包括:
分析自适应滤波器输出信号与近端信号之间的关联程度;
根据所述输出信号与近端信号之间的关联程度计算所述自适应滤波器的时变系数;
根据所述时变系数分析所述自适应滤波器与回声传输信道的逼近程度,据以判断所述时变系数是否满足收敛条件;
若满足收敛条件,则将所述时变系数存入缓存器中以形成缓存系数;
若不满足收敛条件,则将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声。
10.根据权利要求9所述的回声消除方法,其特征在于,根据所述时变系数分析所述自适应滤波器与回声传输信道的逼近程度,据以判断所述时变系数是否满足收敛条件;若满足收敛条件,则将所述时变系数存入缓存器中以形成缓存系数;若不满足收敛条件,则将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声,具体包括:
分别比较所述时变系数与第一阈值和第二阈值之间的大小关系;其中,所述第一阈值大于第二阈值;
若所述时变系数大于所述第一阈值,则所述时变系数收敛到较好程度,所述回声消除误差并且:
判断所述时变系数与所述缓存器前一次保存的缓存系数的收敛程度;
若所述时变系数的收敛程度高于所述缓存器前一次保存的缓存系数的收敛程度,
则将所述时变系数存入所述缓存器中以作为此次保存的缓存系数;
若所述时变系数小于所述第二阈值,则变所述时系数发散;将所述缓存器前一次保存的缓存系数载入辅助滤波器中,以通过所述辅助滤波器消除回声,所述回声消除误差其中,为所述辅助滤波器的输出信号;
若所述时变系数大于等于第二阈值,且小于等于第一阈值,则判定为无所述近端信号,所述回声消除误差
11.根据权利要求10所述的回声消除方法,其特征在于,所述判断所述时变系数与所述缓存器前一次保存的缓存系数的收敛程度,具体包括:
比较所述时变系数r的峰值rpeak与所述缓存器前一次保存的缓存系数r′的峰值r′peak,且比较所述时变系数r的峰值位置kpeak与缓存系数r′的峰值位置k′peak;
若rpeak>r′peak或者kpeak<k′peak+to,则所述时变系数r的收敛程度甚于所述缓存系数,且令r′peak=rpeak,k′peak=kpeak;其中,to表示所述时变系数r最多维持峰值的有效时间。
12.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现权利要求9至11中任一项所述的回声消除方法。
13.一种回声消除终端,其特征在于,包括:处理器及存储器;
所述存储器用于存储计算机程序,所述处理器用于执行所述存储器存储的计算机程序,以使所述终端执行权利要求9至11中任一项所述的回声消除方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711437294.4A CN109961798B (zh) | 2017-12-26 | 2017-12-26 | 回声消除系统、方法、可读计算机存储介质、及终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711437294.4A CN109961798B (zh) | 2017-12-26 | 2017-12-26 | 回声消除系统、方法、可读计算机存储介质、及终端 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109961798A true CN109961798A (zh) | 2019-07-02 |
CN109961798B CN109961798B (zh) | 2021-06-11 |
Family
ID=67022593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711437294.4A Active CN109961798B (zh) | 2017-12-26 | 2017-12-26 | 回声消除系统、方法、可读计算机存储介质、及终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109961798B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111031448A (zh) * | 2019-11-12 | 2020-04-17 | 西安讯飞超脑信息科技有限公司 | 回声消除方法、装置、电子设备和存储介质 |
CN111277718A (zh) * | 2020-01-21 | 2020-06-12 | 上海推乐信息技术服务有限公司 | 一种回声消除系统及其方法 |
CN111370015A (zh) * | 2020-02-28 | 2020-07-03 | 北京字节跳动网络技术有限公司 | 回声消除方法、装置、电子设备及存储介质 |
CN111445918A (zh) * | 2020-03-23 | 2020-07-24 | 深圳市友杰智新科技有限公司 | 降低智能语音设备误唤醒的方法、装置和计算机设备 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0053202A1 (en) * | 1980-12-06 | 1982-06-09 | Communications Satellite Corporation | Improved double talk detector for echo cancellers |
US6766019B1 (en) * | 2000-07-21 | 2004-07-20 | Agere Systems Inc. | Method and apparatus for performing double-talk detection in acoustic echo cancellation |
CN101461257A (zh) * | 2006-05-30 | 2009-06-17 | 微软公司 | 自适应声学回声消除 |
CN101917527A (zh) * | 2010-09-02 | 2010-12-15 | 杭州华三通信技术有限公司 | 回声消除的方法和装置 |
CN102065190A (zh) * | 2010-12-31 | 2011-05-18 | 杭州华三通信技术有限公司 | 一种回声消除方法及其装置 |
CN102117620A (zh) * | 2010-01-06 | 2011-07-06 | 杭州华三通信技术有限公司 | 一种双滤波器传递滤波器系数的方法及装置 |
CN104202018A (zh) * | 2014-08-22 | 2014-12-10 | 华为技术有限公司 | 确定自适应滤波器的稳定因子的方法和装置 |
CN107026950A (zh) * | 2017-05-04 | 2017-08-08 | 重庆第二师范学院 | 一种频域自适应回声消除方法 |
CN107071196A (zh) * | 2017-05-04 | 2017-08-18 | 重庆第二师范学院 | 一种自适应回声消除方法 |
CN107123430A (zh) * | 2017-04-12 | 2017-09-01 | 广州视源电子科技股份有限公司 | 回声消除方法、装置、会议平板及计算机存储介质 |
CN107333017A (zh) * | 2017-05-22 | 2017-11-07 | 西南交通大学 | 一种重用权系数指数型误差自适应回声消除方法 |
-
2017
- 2017-12-26 CN CN201711437294.4A patent/CN109961798B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0053202A1 (en) * | 1980-12-06 | 1982-06-09 | Communications Satellite Corporation | Improved double talk detector for echo cancellers |
US6766019B1 (en) * | 2000-07-21 | 2004-07-20 | Agere Systems Inc. | Method and apparatus for performing double-talk detection in acoustic echo cancellation |
CN101461257A (zh) * | 2006-05-30 | 2009-06-17 | 微软公司 | 自适应声学回声消除 |
CN102117620A (zh) * | 2010-01-06 | 2011-07-06 | 杭州华三通信技术有限公司 | 一种双滤波器传递滤波器系数的方法及装置 |
CN101917527A (zh) * | 2010-09-02 | 2010-12-15 | 杭州华三通信技术有限公司 | 回声消除的方法和装置 |
CN102065190A (zh) * | 2010-12-31 | 2011-05-18 | 杭州华三通信技术有限公司 | 一种回声消除方法及其装置 |
CN104202018A (zh) * | 2014-08-22 | 2014-12-10 | 华为技术有限公司 | 确定自适应滤波器的稳定因子的方法和装置 |
CN107123430A (zh) * | 2017-04-12 | 2017-09-01 | 广州视源电子科技股份有限公司 | 回声消除方法、装置、会议平板及计算机存储介质 |
CN107026950A (zh) * | 2017-05-04 | 2017-08-08 | 重庆第二师范学院 | 一种频域自适应回声消除方法 |
CN107071196A (zh) * | 2017-05-04 | 2017-08-18 | 重庆第二师范学院 | 一种自适应回声消除方法 |
CN107333017A (zh) * | 2017-05-22 | 2017-11-07 | 西南交通大学 | 一种重用权系数指数型误差自适应回声消除方法 |
Non-Patent Citations (3)
Title |
---|
JACOB BENESTY: ""a new class of doubletalk detectors based on cross-correlation"", 《IEEE TRANSACTIONS ON SPEECH AND AUDIO PROCESSING》 * |
宋彦: ""带双端发声检测的回声消除改进算法"", 《清华大学学报自然科学版》 * |
汪萌: ""一种使用双滤波器的回声消除算法"", 《数据采集与处理》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111031448A (zh) * | 2019-11-12 | 2020-04-17 | 西安讯飞超脑信息科技有限公司 | 回声消除方法、装置、电子设备和存储介质 |
CN111277718A (zh) * | 2020-01-21 | 2020-06-12 | 上海推乐信息技术服务有限公司 | 一种回声消除系统及其方法 |
CN111277718B (zh) * | 2020-01-21 | 2021-10-08 | 上海推乐信息技术服务有限公司 | 一种回声消除系统及其方法 |
CN111370015A (zh) * | 2020-02-28 | 2020-07-03 | 北京字节跳动网络技术有限公司 | 回声消除方法、装置、电子设备及存储介质 |
CN111370015B (zh) * | 2020-02-28 | 2021-03-19 | 北京字节跳动网络技术有限公司 | 回声消除方法、装置、电子设备及存储介质 |
CN111445918A (zh) * | 2020-03-23 | 2020-07-24 | 深圳市友杰智新科技有限公司 | 降低智能语音设备误唤醒的方法、装置和计算机设备 |
CN111445918B (zh) * | 2020-03-23 | 2020-12-15 | 深圳市友杰智新科技有限公司 | 降低智能语音设备误唤醒的方法、装置和计算机设备 |
Also Published As
Publication number | Publication date |
---|---|
CN109961798B (zh) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109961798A (zh) | 回声消除系统、方法、可读计算机存储介质、及终端 | |
CN103718538B (zh) | 可实现声学回声消除的音频信号非线性后期处理方法和系统 | |
Shan et al. | Adaptive algorithms with an automatic gain control feature | |
CN100524466C (zh) | 一种麦克风回声消除装置及回声消除方法 | |
CN108711433B (zh) | 一种回声消除方法和装置 | |
US20100135483A1 (en) | Systems and methods for double-talk detection in acoustically harsh environments | |
CN109102821B (zh) | 时延估计方法、系统、存储介质及电子设备 | |
WO2017099728A1 (en) | System and method for suppression of non-linear acoustic echoes | |
US20210013927A1 (en) | Double talk detection method, double talk detection apparatus and echo cancellation system | |
CN110402540A (zh) | 主动降噪方法、装置、芯片、主动控制系统和存储介质 | |
US8831210B2 (en) | Method and system for detection of onset of near-end signal in an echo cancellation system | |
CN110138990A (zh) | 一种消除移动设备VoIP电话回声的方法 | |
CN110431624A (zh) | 残余回声检测方法、残余回声检测装置、语音处理芯片及电子设备 | |
CN103688522B (zh) | 时钟漂移补偿方法和装置 | |
CN109727605A (zh) | 处理声音信号的方法及系统 | |
CN109785853A (zh) | 一种回声消除方法、装置、系统及存储介质 | |
CN102739286B (zh) | 一种用于通信系统中的回声抵消方法 | |
CN107071196A (zh) | 一种自适应回声消除方法 | |
CN103327201A (zh) | 残留回声消除方法及系统 | |
CN113345459A (zh) | 一种双讲状态的检测方法、装置、计算机设备及存储介质 | |
CN111355855B (zh) | 回声处理方法、装置、设备及存储介质 | |
AU4869200A (en) | Apparatus and method for estimating an echo path delay | |
Hafsi et al. | Identification of Wiener-Hammerstein model with multisegment piecewise-linear characteristic | |
CN107071197A (zh) | 一种基于全相位多延迟分块频域的回音消除方法及系统 | |
US20120140939A1 (en) | Method and device for cancelling acoustic echo |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |