CN107071196A - 一种自适应回声消除方法 - Google Patents
一种自适应回声消除方法 Download PDFInfo
- Publication number
- CN107071196A CN107071196A CN201710306316.7A CN201710306316A CN107071196A CN 107071196 A CN107071196 A CN 107071196A CN 201710306316 A CN201710306316 A CN 201710306316A CN 107071196 A CN107071196 A CN 107071196A
- Authority
- CN
- China
- Prior art keywords
- mrow
- signal
- adaptive
- frequency domain
- sef
- 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
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 12
- 230000021615 conjugation Effects 0.000 claims abstract description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 description 5
- 230000002452 interceptive effect Effects 0.000 description 3
- 239000012636 effector Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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/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
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
- G10L21/0232—Processing in the frequency domain
-
- 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
- G10L21/0264—Noise filtering characterised by the type of parameter measurement, e.g. correlation techniques, zero crossing techniques or predictive techniques
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)
- 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)作为自适应滤波器的自适应步长因子,更新自适应滤波器系数。本发明的技术方案将残留信号与参考信号之间的相关系数作为自适应步长因子,对各话音状态均可进行有效的回声消除。
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)作为自适应滤波器的自适应步长因子μ(f),更新自适应滤波器系数;
其中,所述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
其中,
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)作为自适应滤波器的自适应步长因子μ(f),更新自适应滤波器系数;
本实施例中,
本实施例中,本步骤可以根据所述X(f)、所述E(f)和所述cohxe(f)直接更新自适应滤波器系数:
其中,所述k为时刻,k=1,2……M,M为所述自适应滤波器长度;所述f为频点;Wk+1(f)为频点f的更新后的自适滤波器系数;Wk(f)为频点f的更新前的自适应滤波器系数。
作为本实施例的一种优选实现方案,本步骤还可以包括:
S701、对所述频域残留信号E(f)进行归一化处理,获得归一化后的频域残留信号Enor(f);
S702、根据所述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);
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>h</mi>
<mi>x</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>x</mi>
<mi>e</mi>
<mi>P</mi>
<mi>o</mi>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>x</mi>
<mi>P</mi>
<mi>o</mi>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>e</mi>
<mi>P</mi>
<mi>o</mi>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
将cohxe(f)作为自适应滤波器的自适应步长因子μ(f),更新自适应滤波器系数;
其中,所述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)更新自适应滤波器系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710306316.7A CN107071196B (zh) | 2017-05-04 | 2017-05-04 | 一种自适应回声消除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710306316.7A CN107071196B (zh) | 2017-05-04 | 2017-05-04 | 一种自适应回声消除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107071196A true CN107071196A (zh) | 2017-08-18 |
CN107071196B CN107071196B (zh) | 2019-09-03 |
Family
ID=59596048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710306316.7A Active CN107071196B (zh) | 2017-05-04 | 2017-05-04 | 一种自适应回声消除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107071196B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109961798A (zh) * | 2017-12-26 | 2019-07-02 | 华平信息技术股份有限公司 | 回声消除系统、方法、可读计算机存储介质、及终端 |
CN110838300A (zh) * | 2019-11-18 | 2020-02-25 | 紫光展锐(重庆)科技有限公司 | 回声消除的处理方法及处理系统 |
CN112017679A (zh) * | 2020-08-05 | 2020-12-01 | 海尔优家智能科技(北京)有限公司 | 用于自适应滤波器系数更新的方法及装置、设备 |
CN112397080A (zh) * | 2020-10-30 | 2021-02-23 | 浙江大华技术股份有限公司 | 回声消除方法及装置、语音设备及计算机可读存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1806423A (zh) * | 2003-06-27 | 2006-07-19 | 诺基亚有限公司 | 使用残余回声滤波器增强声音回声消除系统的方法 |
CN101179294A (zh) * | 2006-11-09 | 2008-05-14 | 爱普拉斯通信技术(北京)有限公司 | 自适应回声消除器及其回声消除方法 |
CN101222555A (zh) * | 2008-01-25 | 2008-07-16 | 上海华平信息技术股份有限公司 | 改善音频通话质量的系统及方法 |
CN102387272A (zh) * | 2011-09-09 | 2012-03-21 | 南京大学 | 一种回声抵消系统中残留回声的抑制方法 |
JP4978352B2 (ja) * | 2007-07-11 | 2012-07-18 | ヤマハ株式会社 | エコーキャンセラ |
KR20150054424A (ko) * | 2013-11-12 | 2015-05-20 | 고려대학교 산학협력단 | 잔여 에코 제거 방법 및 그 장치 |
CN104954595A (zh) * | 2015-05-15 | 2015-09-30 | 百度在线网络技术(北京)有限公司 | 残留回声消除方法和装置 |
-
2017
- 2017-05-04 CN CN201710306316.7A patent/CN107071196B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1806423A (zh) * | 2003-06-27 | 2006-07-19 | 诺基亚有限公司 | 使用残余回声滤波器增强声音回声消除系统的方法 |
CN101179294A (zh) * | 2006-11-09 | 2008-05-14 | 爱普拉斯通信技术(北京)有限公司 | 自适应回声消除器及其回声消除方法 |
JP4978352B2 (ja) * | 2007-07-11 | 2012-07-18 | ヤマハ株式会社 | エコーキャンセラ |
CN101222555A (zh) * | 2008-01-25 | 2008-07-16 | 上海华平信息技术股份有限公司 | 改善音频通话质量的系统及方法 |
CN102387272A (zh) * | 2011-09-09 | 2012-03-21 | 南京大学 | 一种回声抵消系统中残留回声的抑制方法 |
KR20150054424A (ko) * | 2013-11-12 | 2015-05-20 | 고려대학교 산학협력단 | 잔여 에코 제거 방법 및 그 장치 |
CN104954595A (zh) * | 2015-05-15 | 2015-09-30 | 百度在线网络技术(北京)有限公司 | 残留回声消除方法和装置 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109961798A (zh) * | 2017-12-26 | 2019-07-02 | 华平信息技术股份有限公司 | 回声消除系统、方法、可读计算机存储介质、及终端 |
CN109961798B (zh) * | 2017-12-26 | 2021-06-11 | 华平信息技术股份有限公司 | 回声消除系统、方法、可读计算机存储介质、及终端 |
CN110838300A (zh) * | 2019-11-18 | 2020-02-25 | 紫光展锐(重庆)科技有限公司 | 回声消除的处理方法及处理系统 |
CN110838300B (zh) * | 2019-11-18 | 2022-03-25 | 紫光展锐(重庆)科技有限公司 | 回声消除的处理方法及处理系统 |
CN112017679A (zh) * | 2020-08-05 | 2020-12-01 | 海尔优家智能科技(北京)有限公司 | 用于自适应滤波器系数更新的方法及装置、设备 |
CN112017679B (zh) * | 2020-08-05 | 2024-01-26 | 海尔优家智能科技(北京)有限公司 | 用于自适应滤波器系数更新的方法及装置、设备 |
CN112397080A (zh) * | 2020-10-30 | 2021-02-23 | 浙江大华技术股份有限公司 | 回声消除方法及装置、语音设备及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107071196B (zh) | 2019-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109727604B (zh) | 用于语音识别前端的频域回声消除方法及计算机储存介质 | |
CN107134281A (zh) | 一种自适应回声消除中自适应滤波器系数更新方法 | |
CN107026950B (zh) | 一种频域自适应回声消除方法 | |
CN107123430B (zh) | 回声消除方法、装置、会议平板及计算机存储介质 | |
CN100524466C (zh) | 一种麦克风回声消除装置及回声消除方法 | |
CN104883462B (zh) | 一种用于消除声学回声的自适应滤波器及滤波方法 | |
CN111768796B (zh) | 一种声学回波消除与去混响方法及装置 | |
CN107071196A (zh) | 一种自适应回声消除方法 | |
CN111885275B (zh) | 语音信号的回声消除方法、装置、存储介质以及电子装置 | |
CN105810202B (zh) | 一种降低回声的方法、装置及通讯设备 | |
CN104395957A (zh) | 一种通用可重构回声消除系统 | |
US10880440B2 (en) | Echo canceller and method therefor | |
TW200948030A (en) | Apparatus and method for computing filter coefficients for echo suppression | |
CN106157964A (zh) | 一种确定回声消除中系统延时的方法 | |
CN110992923B (zh) | 回声消除方法、电子设备以及存储装置 | |
CN111755020B (zh) | 一种立体声回声消除方法 | |
Albu et al. | Proportionate algorithms for two-microphone active feedback cancellation | |
Ivry et al. | Nonlinear acoustic echo cancellation with deep learning | |
EP2930917B1 (en) | Method and apparatus for updating filter coefficients of an adaptive echo canceller | |
CN112259112A (zh) | 一种结合声纹识别和深度学习的回声消除方法 | |
CN106161820B (zh) | 一种用于立体声声学回声抵消的通道间去相关方法 | |
CN115457928A (zh) | 一种基于神经网络双讲检测的回声消除方法及系统 | |
CN105491256A (zh) | 一种声学回声消除器启动阶段稳健的步长调整方法 | |
CN204117590U (zh) | 语音采集降噪装置以及语音质量评价系统 | |
CN111370016B (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 |