CN107026950A - 一种频域自适应回声消除方法 - Google Patents

一种频域自适应回声消除方法 Download PDF

Info

Publication number
CN107026950A
CN107026950A CN201710306327.5A CN201710306327A CN107026950A CN 107026950 A CN107026950 A CN 107026950A CN 201710306327 A CN201710306327 A CN 201710306327A CN 107026950 A CN107026950 A CN 107026950A
Authority
CN
China
Prior art keywords
frequency domain
signal
calculate
adaptive
cohxe
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
CN201710306327.5A
Other languages
English (en)
Other versions
CN107026950B (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.)
Chongqing University of Education
Original Assignee
Chongqing University of Education
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 Chongqing University of Education filed Critical Chongqing University of Education
Priority to CN201710306327.5A priority Critical patent/CN107026950B/zh
Publication of CN107026950A publication Critical patent/CN107026950A/zh
Application granted granted Critical
Publication of CN107026950B publication Critical patent/CN107026950B/zh
Active 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/085—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 digital techniques

Landscapes

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

Abstract

本发明公开了一种频域自适应回声消除方法，包括：根据发音单元的参考信号x(k)通过自适应滤波器估计回声信号y(k)；根据所述y(k)和麦克风拾取的信号d(k)计算残留信号e(k)；计算频域参考信号X(f)的功率谱xPow(f)；计算频域残留信号E(f)的功率谱ePow(f)；计算所述E(f)的共轭残留信号E*(f)与所述X(f)的相关功率谱xePow(f)；计算X(f)与E*(f)的相关系数cohxe(f)；计算所有频点的cohxe(f)的平均值cohxeaver；将cohxeaver作为自适应滤波器在各频点的自适应步长因子，更新自适应滤波器系数。本发明的技术方案将残留信号与参考信号之间的相关系数的平均值作为自适应步长因子，对各话音状态均可进行有效的回声消除。

Description

一种频域自适应回声消除方法
技术领域
本发明涉及到实时语音通信、人机交互过程中的回声处理技术，特别涉及到一种频域自适应回声消除方法。
背景技术
在实时语音通信过程中，通信终端扬声器发出的声音，总会被该终端的麦克风拾取到，若是不处理就发送出去，对方总能听到自己说话的声音；在人机交互领域，由于交互终端发出的声音又被麦克风拾取回去，同时拾取了控制者的说话声，若在麦克风拾取信号中不消除交互终端发出的声音，那么交互终端在识别控制者说话声音时将引入很强的干扰，降低了识别的成功率，最终造成交互困难。
针对实时语音通信过程中的回声干扰问题，现有技术也提出了一些技术方案，回声消除可以在时域或频域进行，由于在时域进行回声消除的复杂度较高，通常都是在频域进行回声消除，典型的频域回声消除方法如中国专利申请“一种用于消除声学回声的自适应滤波器及滤波方法”(申请号:201410073711.1)中所提出一种用于回声消除的自适应滤波方案。该技术方案是一种频域自适应滤波方法，利用远端信号x(n)作为参考信号，通过自适应滤波对传播路径h(n)进行系统估计，并利用估计的系统参数w(n)对参考信号进行滤波，进而消除终端麦克风采集到的信号d(n)中的回声得到残差信号e(n)。一方面将e(n)作为滤波结果直接输出，另一方面利用e(n)对自适应滤波器系数进行迭代控制。该方案中在对滤波器系数进行迭代控制时，自适应滤波器每迭代R次对滤波器系数限制参数进行一次更新，其中R是预先设定的固定值。
现有方案的问题在于：一方面，自适应滤波器系数的更新与滤波效果(残留信号)没有关联起来，且没有根据回声强度实时调整自适应步长因子，不能很好的保证在各种不同的话音状态下回声消除的精确度。
发明内容
为了解决现有技术中存在的上述问题，本发明提出了一种频域自适应回声消除方法，以实现在各种话音状态下有效的自适应回声消除。
为了实现上述目的，本发明采用了以下技术方案：
一种频域自适应回声消除方法，包括：
根据发音单元的参考信号x(k)通过自适应滤波器估计回声信号y(k)；
根据所述y(k)和麦克风拾取的信号d(k)计算残留信号e(k):
e(k)＝d(k)-y(k)
计算频域参考信号X(f)的功率谱xPow(f)；
计算频域残留信号E(f)的功率谱ePow(f)；
其中，E(f)＝FFT[0M个0,e(k)]；
计算所述E(f)的共轭残留信号E*(f)与所述X(f)的相关功率谱xePow(f)；
xePow(f)＝||X(f)·E*(f)||2
计算X(f)与E*(f)的相关系数cohxe(f)；
计算所有频点的cohxe(f)的平均值cohxeaver；
其中，所述F为频点数量；
将cohxeaver作为自适应滤波器在各频点的自适应步长因子，更新自适应滤波器系数
其中，所述k为时刻，k＝1,2……M,M为所述自适应滤波器长度；所述f为频点；Wk+1(f)为频点f的更新后的自适滤波器系数；Wk(f)为频点f的更新前的自适应滤波器系数。
进一步的，所述估计回声信号y(k)包括：
将扬声器单元的参考信号x(k)转换为频域参考信号X(f)；
X(f)＝FFT[x(k-M),...,x(k),...,x(k+M-1)]
估计频域回声信号Y(f)，
对所述Y(f)进行快速傅立叶逆变换获得y(k)；
y(k)＝IFFT[Y(f)]
其中，所述IFFT为快速傅立叶逆变换，为频域卷积操作；W(f)为自适应滤波器当前在频点f所使用的自适应滤波器系数。
进一步，所述更新自适应滤波器系数包括：
对所述频域残留信号E(f)进行归一化处理，获得归一化后的频域残留信号Enor(f)；
根据所述Enor(f)更新自适应滤波器系数。
本发明的技术方案将残留信号与参考信号之间的相关系数作为自适应步长因子，对各话音状态均可进行极大限度的回声消除。通话状态通常分为四种，远端单讲、近端单讲、双端讲话、双端静默。对远端单讲状态，麦克风拾取的是回声信号，在自适应滤波器的初始阶段，自适应滤波器有一个收敛过程，残留信号与麦克风拾取的回声信号强度相当，那么参考信号与残留信号直接的相关性很强，则对应的自适应步长因子也较大，当自适应滤波器进入稳定状态时，残留信号几乎为0，此时回声信号与残留信号的相关性非常弱，仅对自适应滤波器系数进行微调。对近端单讲状态，残留信号为麦克风拾取的本端话音，参考信号为0，此时计算出的相关系数为0，无需对自适应滤波器进行更新，而实际上此时并没有回声，确实不需进行回声消除。双端静默状态，残留信号与参考信号都很弱，相关系数也会很小。同时，本发明技术方案中针自适应步长因子，处理后频点之间仍保持原有的数据关系。本发明的技术方案能有效的提高回声消除的精确度。
附图说明
图1本发明频域自适应回声消除方法流程图。
具体实施方式
为了更好的说明本发明的技术方案，下面结合附图对本发明的具体实施方式进行详细描述。
具体实施例
本实施例为本发明频域自适应回声消除方法的一种优选实施方式。
参见图1，本实施例的频域自适应回声消除方法流程如图1所示，包括：
S1、根据发音单元的参考信号x(k)通过自适应滤波器估计回声信号y(k)；
作为本实施例的一种优选实现方案，本步骤还可以经一部包括：
S101、将扬声器单元的参考信号x(k)转换为频域参考信号X(f)。
X(f)＝FFT[x(k-M),...,x(k),...,x(k+M-1)]
S102、估计频域回声信号Y(f)，
其中，W(f)为当前的自适应滤波器系数；
S103、根据所述Y(f)获得所述y(k)：
y(k)＝IFFT[Y(f)]
其中，所述FFT为傅立叶变换，所述IFFT为傅立叶逆变换；
本是实例中，本步骤也可以通过其他方式对回声信号进行估计，本发明对此没有限制。
S2、根据所述y(k)和麦克风拾取的信号d(k)计算残留信号e(k):
e(k)＝d(k)-y(k)
残留信号e(k)作为滤波消除回声干扰后的结果输出。
S3、计算频域参考信号X(f)的功率谱xPow(f)：
xPow(f)＝||X(f)||2
S4、计算频域残留信号E(f)的功率谱ePow(f)：
ePow(f)＝||E(f)||2
其中，E(f)＝FFT[0M个0,e(k)]；
S5、计算所述E(f)的共轭残留信号E*(f)与所述X(f)的相关功率谱
xePow(f)；
xePow(f)＝||X(f)·E*(f)||2
S6、计算X(f)与E*(f)的相关系数cohxe(f)；
S7、计算所有频点的cohxe(f)的平均值cohxeaver；
其中，所述F为频点数量；
S8、将cohxeaver作为自适应滤波器的自适应步长因子μ(f)，更新自适应滤波器系数；
本实施例中，
本实施例中，本步骤可以根据所述X(f)、所述E(f)和所述cohxe(f)直接更新自适应滤波器系数：
其中，所述k为时刻，k＝1,2……M,M为所述自适应滤波器长度；所述f为频点；Wk+1(f)为频点f的更新后的自适滤波器系数；Wk(f)为频点f的更新前的自适应滤波器系数。
作为本实施例的一种优选实现方案，本步骤还可以包括：
S801、对所述频域残留信号E(f)进行归一化处理，获得归一化后的频域残留信号Enor(f)；
S802、根据所述Enor(f)更新自适应滤波器系数。
其中，所述k为时刻，k＝1,2……M,M为所述自适应滤波器长度；所述f为频点；Wk+1 (f)为频点f的更新后的自适滤波器系数；Wk(f)为频点f的更新前的自适应滤波器系数。
更新后的自适应滤波器系数Wk+1(f)将用于消除下一个时刻麦克风拾取信号d(k+1)中的回声信号y(k+1)以得到下一时刻的残留信号e(k+1)。
双端讲话状态时，自适应滤波器是需要更新的，且参考信号与残留信号的相关性很强，但作为反馈信号的残留信号包含有本端话音，若直接反馈到自适应滤波器系数，势必会造成滤波器发散，本优选实现方案首先将频域残留信号进行归一化处理，再利用归一化后的频域残留信号来更新自适应滤波器系数，能很好地解决滤波器发散的问题。
需要说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (3)

1.一种频域自适应回声消除方法，包括：
根据发音单元的参考信号x(k)通过自适应滤波器估计回声信号y(k)；
根据所述y(k)和麦克风拾取的信号d(k)计算残留信号e(k):
e(k)＝d(k)-y(k)
计算频域参考信号X(f)的功率谱xPow(f)；
计算频域残留信号E(f)的功率谱ePow(f)；
计算所述E(f)的共轭残留信号E*(f)与所述X(f)的相关功率谱xePow(f)；
计算X(f)与E*(f)的相关系数cohxe(f)；
计算所有频点的cohxe(f)的平均值cohxeaver；
将cohxeaver作为自适应滤波器在各频点的自适应步长因子，更新自适应滤波器系数；
其中，所述k为时刻，所述f为频点,k＝1,2……M,M为所述自适应滤波器长度。
2.根据权利要求1所述的方法，其特征在于，所述估计回声信号y(k)包括：
将扬声器单元的参考信号x(k)转换为频域参考信号X(f)；
估计频域回声信号Y(f)，
对所述Y(f)进行快速傅立叶逆变换获得y(k)；
其中，为频域卷积操作；W(f)为自适应滤波器当前在频点f所使用的自适应滤波器系数。
3.根据权利要求1或2所述的方法，其特征在于，所述更新自适应滤波器系数包括：
对所述频域残留信号E(f)进行归一化处理，获得归一化后的频域残留信号Enor(f)；
根据所述Enor(f)更新自适应滤波器系数。
CN201710306327.5A 2017-05-04 2017-05-04 一种频域自适应回声消除方法 Active CN107026950B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710306327.5A CN107026950B (zh) 2017-05-04 2017-05-04 一种频域自适应回声消除方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710306327.5A CN107026950B (zh) 2017-05-04 2017-05-04 一种频域自适应回声消除方法

Publications (2)

Publication Number Publication Date
CN107026950A true CN107026950A (zh) 2017-08-08
CN107026950B CN107026950B (zh) 2019-09-17

Family

ID=59527198

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710306327.5A Active CN107026950B (zh) 2017-05-04 2017-05-04 一种频域自适应回声消除方法

Country Status (1)

Country Link
CN (1) CN107026950B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109727604A (zh) * 2018-12-14 2019-05-07 上海蔚来汽车有限公司 用于语音识别前端的频域回声消除方法及计算机储存介质
CN109961798A (zh) * 2017-12-26 2019-07-02 华平信息技术股份有限公司 回声消除系统、方法、可读计算机存储介质、及终端
CN110646769A (zh) * 2019-09-03 2020-01-03 武汉大学深圳研究院 一种适用于lte外辐射源雷达的时域杂波抑制方法
CN111355855A (zh) * 2020-03-12 2020-06-30 紫光展锐(重庆)科技有限公司 回声处理方法、装置、设备及存储介质
CN111742541A (zh) * 2017-12-08 2020-10-02 华为技术有限公司 声学回波抵消方法和装置
CN112017679A (zh) * 2020-08-05 2020-12-01 海尔优家智能科技（北京）有限公司 用于自适应滤波器系数更新的方法及装置、设备
CN113362844A (zh) * 2021-07-26 2021-09-07 西南交通大学 一种低复杂度分离去相关自适应声学回声消除方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102387272A (zh) * 2011-09-09 2012-03-21 南京大学 一种回声抵消系统中残留回声的抑制方法
EP2701145A1 (en) * 2012-08-24 2014-02-26 Retune DSP ApS Noise estimation for use with noise reduction and echo cancellation in personal communication
CN105794190A (zh) * 2013-12-12 2016-07-20 皇家飞利浦有限公司 回声消除

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102387272A (zh) * 2011-09-09 2012-03-21 南京大学 一种回声抵消系统中残留回声的抑制方法
EP2701145A1 (en) * 2012-08-24 2014-02-26 Retune DSP ApS Noise estimation for use with noise reduction and echo cancellation in personal communication
CN105794190A (zh) * 2013-12-12 2016-07-20 皇家飞利浦有限公司 回声消除

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111742541A (zh) * 2017-12-08 2020-10-02 华为技术有限公司 声学回波抵消方法和装置
CN111742541B (zh) * 2017-12-08 2021-11-30 华为技术有限公司 声学回波抵消方法、装置、存储介质
CN109961798A (zh) * 2017-12-26 2019-07-02 华平信息技术股份有限公司 回声消除系统、方法、可读计算机存储介质、及终端
CN109961798B (zh) * 2017-12-26 2021-06-11 华平信息技术股份有限公司 回声消除系统、方法、可读计算机存储介质、及终端
CN109727604A (zh) * 2018-12-14 2019-05-07 上海蔚来汽车有限公司 用于语音识别前端的频域回声消除方法及计算机储存介质
CN109727604B (zh) * 2018-12-14 2023-11-10 上海蔚来汽车有限公司 用于语音识别前端的频域回声消除方法及计算机储存介质
CN110646769B (zh) * 2019-09-03 2021-07-20 武汉大学深圳研究院 一种适用于lte外辐射源雷达的时域杂波抑制方法
CN110646769A (zh) * 2019-09-03 2020-01-03 武汉大学深圳研究院 一种适用于lte外辐射源雷达的时域杂波抑制方法
CN111355855A (zh) * 2020-03-12 2020-06-30 紫光展锐(重庆)科技有限公司 回声处理方法、装置、设备及存储介质
CN111355855B (zh) * 2020-03-12 2021-06-15 紫光展锐(重庆)科技有限公司 回声处理方法、装置、设备及存储介质
CN112017679A (zh) * 2020-08-05 2020-12-01 海尔优家智能科技（北京）有限公司 用于自适应滤波器系数更新的方法及装置、设备
CN112017679B (zh) * 2020-08-05 2024-01-26 海尔优家智能科技（北京）有限公司 用于自适应滤波器系数更新的方法及装置、设备
CN113362844A (zh) * 2021-07-26 2021-09-07 西南交通大学 一种低复杂度分离去相关自适应声学回声消除方法及装置
CN113362844B (zh) * 2021-07-26 2022-05-10 西南交通大学 一种低复杂度分离去相关自适应声学回声消除方法及装置

Also Published As

Publication number Publication date
CN107026950B (zh) 2019-09-17

Similar Documents

Publication Publication Date Title
CN107026950B (zh) 一种频域自适应回声消除方法
CN107134281A (zh) 一种自适应回声消除中自适应滤波器系数更新方法
CN107123430B (zh) 回声消除方法、装置、会议平板及计算机存储介质
CN109727604B (zh) 用于语音识别前端的频域回声消除方法及计算机储存介质
US9697846B2 (en) Method and apparatus for reducing echo
CN110838300B (zh) 回声消除的处理方法及处理系统
KR101469739B1 (ko) 오디오 신호들을 처리하는 디바이스 및 방법
CN111768796B (zh) 一种声学回波消除与去混响方法及装置
CN104883462B (zh) 一种用于消除声学回声的自适应滤波器及滤波方法
CN107071196B (zh) 一种自适应回声消除方法
CN109754813B (zh) 基于快速收敛特性的变步长回声消除方法
CN104395957A (zh) 一种通用可重构回声消除系统
CN105810202B (zh) 一种降低回声的方法、装置及通讯设备
CN110992923B (zh) 回声消除方法、电子设备以及存储装置
CN111755020B (zh) 一种立体声回声消除方法
CN110956975A (zh) 回声消除方法及装置
CN112259112A (zh) 一种结合声纹识别和深度学习的回声消除方法
EP2930917A1 (en) Method and apparatus for updating filter coefficients of an adaptive echo canceller
CN114242100B (zh) 音频信号处理方法、训练方法及其装置、设备、存储介质
WO2014000658A1 (zh) 消除噪音的方法和装置、以及移动终端
CN115457928A (zh) 一种基于神经网络双讲检测的回声消除方法及系统
CN107005268A (zh) 回波消除装置和回波消除方法
WO2015044915A1 (en) Acoustic feedback cancellation based on cesptral analysis
CN111370016B (zh) 一种回声消除方法及电子设备
CN111355855A (zh) 回声处理方法、装置、设备及存储介质

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