CN105635500B - 双麦克风回声及噪声的抑制系统及其方法 - Google Patents
双麦克风回声及噪声的抑制系统及其方法 Download PDFInfo
- Publication number
- CN105635500B CN105635500B CN201410597191.4A CN201410597191A CN105635500B CN 105635500 B CN105635500 B CN 105635500B CN 201410597191 A CN201410597191 A CN 201410597191A CN 105635500 B CN105635500 B CN 105635500B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- echo
- noise
- region
- 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.)
- Active
Links
- 230000009977 dual effect Effects 0.000 title claims abstract description 78
- 230000005764 inhibitory process Effects 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000001629 suppression Effects 0.000 claims abstract description 26
- 238000009499 grossing Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 241000209140 Triticum Species 0.000 claims description 4
- 235000021307 Triticum Nutrition 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000008030 elimination Effects 0.000 description 4
- 238000003379 elimination reaction Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Telephone Function (AREA)
Abstract
本发明提供了一种双麦克风回声及噪声的抑制系统及其方法，短时傅里叶变换模块接收双麦克风终端采集的语音信号并将其转换为频域信号，所述回声估计模块及所述噪声估计模块接收所述频域信号并分别计算所述频域信号中的回声信号的方差及噪声信号的方差，所述抑制模块根据所述频域信号中的回声信号的方差及噪声信号的方差将所述频域信号中的回声信号及噪声信号抑制消除。将抑制噪声信号的方式融合于抑制回声信号的方法中，利用抑制模块将回声信号及噪声信号一同抑制消除，整个过程无需采用自适应滤波器，避免现有技术消除回声所存在的问题，提高了语音质量。
Description
技术领域
本发明涉及通信技术领域，特别涉及一种双麦克风回声及噪声的抑制系统及其方法。
背景技术
随着手机的智能化发展，对多媒体技术方面的要求也越来越高。一方面，手机要求在各种复杂环境和不同场合下使用，比如车站、机场、嘈杂的市区，为了提高手机的通话质量，就要求手机具有很好的抑制声源背景噪声的功能。另一方面，为了适合特殊人群的需要，希望有语音识别功能，其前提就是必须有纯净的语音提供，这也需要强大的噪声消除功能。
针对噪声消除的问题，各种噪声消除装置及噪声消除的方法已逐渐成熟，但是通话的质量还受到回声的影响。针对回声问题，本领域技术人员通常是利用自适应滤波器将语音信号的回声消除，再经过后处理模块将自适应滤波器未消除的回声进一步抑制消除，但是这种方法也存在不足，基于自适应滤波器是是采用残差反馈的方式进行滤波器系数更新，消除回声时需要采用至少一个自适应滤波器，当多个自适应滤波器串联使用时就会存在相互干扰的现象，消除回声的效果不佳。此外由于回声消除的自适应滤波器所采用的抑制回声的算法和双麦克风内的自适应算法会形成相互干扰，因此回声和噪声的消除较难达到平衡，致使手机的通信质量不高，无法达到最优化。
发明内容
本发明的目的在于提供一种双麦克风回声及噪声的抑制系统及其方法，以克服传统消除回声的方法所存在的不足。
为解决上述技术问题，本发明提供一种双麦克风回声及噪声的抑制系统及其方法，所述双麦克风回声及噪声的抑制系统包括：相互连接的短时傅里叶变换模块、回声估计模块、噪声估计模块及抑制模块；所述短时傅里叶变换模块接收双麦克风终端采集的语音信号并将其转换为频域信号，所述回声估计模块及所述噪声估计模块接收所述频域信号并分别计算所述频域信号中的回声信号的方差及噪声信号的方差，所述抑制模块根据所述频域信号中的回声信号的方差及噪声信号的方差将所述频域信号中的回声信号及噪声信号抑制消除。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述回声估计模块计算所述回声信号的方差包括如下步骤：
对所述频域信号进行平滑处理；
利用平滑处理后的频域信号计算比值；
根据所述频域信号及所述比值计算所述回声信号的方差。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述双麦克风终端上设置有第一麦克风及第二麦克风，所述第一麦克风与所述第二麦克风分别分布在所述双麦克风终端的底部和顶部；所述频域信号包括第一频域信号Y1(ω,m)及第二频域信号Y2(ω,m)，所述第一麦克风采集的语音信号转换为频域信号为所述第一频域信号Y1(ω,m)，所述第二麦克风采集的语音信号转换为所述频域信号为第二频域信号Y2(ω,m)，其中，
Y1(ω,m)＝a1X(ω,m)+D1(ω,m)+b1E(ω,m)；
Y2(ω,m)＝a2X(ω,m)+D2(ω,m)+b2E(ω,m)；
m为帧序号，X(ω,m)为所述频域信号中的近端语音信号，E(ω,m)为所述频域信号中的回声信号，D1(ω,m)为所述第一频域信号Y1(ω,m)中的噪声信号，D2(ω,m)为所述第二频域信号Y2(ω,m)中的噪声信号，a1，a2，b1及b2为由于近场信号的传播引起的幅值差。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述对所述频域信号进行平滑处理采用的公式如下：
Φ11(ω,m)＝α1Φ11(ω,m-1)+(1-α1)|Y1(ω,m)|2；
Φ22(ω,m)＝α1Φ22(ω,m-1)+(1-α1)|Y2(ω,m)|2；
其中，α1为平滑因子，Φ11(ω,m)为第一频域信号Y1(ω,m)进行平滑处理后的信号,Φ22(ω,m)为第二频域信号Y2(ω,m)进行平滑处理后的信号。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述利用平滑处理后的频域信号计算比值采用的公式如下：
其中，r(ω,m)为比值，r1(ω,m)为利用第一频域信号Y1(ω,m)进行平滑处理后的信号Φ11(ω,m)计算的第一比值，r2(ω,m)为利用第二频域信号Y2(ω,m)进行平滑处理后的信号Φ22(ω,m)计算的第二比值。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述第一比值r1(ω,m)及第二比值r2(ω,m)的计算公式如下：
可选的，在所述的双麦克风回声及噪声的抑制系统中，根据所述频域信号及所述比值计算所述回声信号的方差所采用的公式如下：
其中，λe(ω,m)为回声信号的方差，β为平滑因子，η为阈值。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述噪声估计模块计算所述噪声信号的方差采用的方法为：最小统计量方式、语音激活检测方式或者最小控制递归平均方式。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述抑制模块抑制消除所述回声信号及所述噪声信号所采用的算法为MMSE-LSA算法。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述抑制模块包括信噪比估计模块及与所述信噪比模块相连的抑制准则模块；其中，所述信噪比估计模块根据所述回声信号的方差及所述噪声信号的方差计算先验信噪比ξ及后验信噪比γ；所述抑制准则模块通过所述先验信噪比ξ及所述后验信噪比γ计算增益因子G(ξ,γ)；其中，
ξ＝λx(ω,m)/λD(ω,m)；
γ＝|Y1(ω,m)+Y2(ω,m)|2/λD(ω,m)；
λD(ω,m)＝λe(ω,m)+λd(ω,m)；
v＝ξγ/(1+ξ)；
λx(ω,m)为所述语音信号的方差，λd(ω,m)为噪声信号的方差。
可选的，在所述的双麦克风回声及噪声的抑制系统中，还包括与所述短时傅里叶变换模块和所述抑制模块连接的短时傅里叶反变换模块，所述短时傅里叶反变换模块将已抑制消除掉所述回声信号及所述噪声信号的频域信号转化为时域信号输出。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述短时傅里叶变换模块采用窗函数为Hamming窗或者平顶窗。
可选的，在所述的双麦克风回声及噪声的抑制系统中，所述双麦克风终端为手机。
本发明还提供一种双麦克风回声及噪声的抑制方法，所述双麦克风回声及噪声的抑制方法，包括如下步骤：短时傅里叶变换模块接收双麦克风终端采集的语音信号并将其转换为频域信号；
回声估计模块及噪声估计模块接收所述频域信号并分别计算所述频域信号中的回声信号的方差及噪声信号的方差；以及
抑制模块根据所述频域信号中的回声信号的方差及噪声信号的方差将所述频域信号中的回声信号及噪声信号抑制消除。
经过发明人不断的实验发现，由于双麦克风终端上所设置的两个麦克风的位置不同，不同时刻两个麦克风所采集的语音信号强弱有所不同，发明人正是利用两个麦克风采集的语音信号的差异来判断不同时刻双麦克风终端所采集语音信号中存在信号的类型，从而有效的、有针对性的去除两个麦克风采集语音信号中需要抑制消除的信号。。
在本发明所提供的双麦克风回声及噪声的抑制系统及其方法中，短时傅里叶变换模块接收双麦克风终端采集的语音信号并将其转换为频域信号，所述回声估计模块及所述噪声估计模块接收所述频域信号并分别计算所述频域信号中的回声信号的方差及噪声信号的方差，所述抑制模块根据所述频域信号中的回声信号的方差及噪声信号的方差将所述频域信号中的回声信号及噪声信号抑制消除。将抑制噪声信号的方式融合于抑制回声信号的方法中，利用抑制模块将回声信号及噪声信号一同抑制消除，整个过程无需采用自适应滤波器，避免现有技术消除回声所存在的问题，提高了语音质量。
附图说明
图1是本发明一实施例中双麦克风回声及噪声的抑制系统抑制回声的框图；
图2是本发明一实施例中回声估计模块的工作原理图；
图3是本发明中回声估计模块计算回声信号的方差的流程图；
图4是使用本发明的双麦克风回声及噪声的抑制系统之前双麦克风终端采集到的声音的语谱图；
图5是使用本发明的双麦克风回声及噪声的抑制系统之后第一麦克风采集到的声音的语谱图。
具体实施方式
以下结合附图和具体实施例对本发明提出的双麦克风回声及噪声的抑制系统及其方法作进一步详细说明。根据下面说明和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
请参考图1，其为本发明一实施例中双麦克风回声及噪声的抑制系统抑制回声的框图，如图1所示，所述的双麦克风回声及噪声的抑制系统包括：相互连接的短时傅里叶变换模块20、回声估计模块40、噪声估计模块30及抑制模块50；所述短时傅里叶变换模块20接收双麦克风终端10采集的语音信号并将其转换为频域信号，所述回声估计模块40及所述噪声估计模块30接收所述频域信号并分别计算所述频域信号中的回声信号的方差λe(ω,m)及噪声信号的方差λd(ω,m)，所述抑制模块50根据所述频域信号中的回声信号的方差λe(ω,m)及噪声信号的方差λd(ω,m)将所述频域信号中的回声信号及噪声信号抑制消除。
较佳的，所述抑制模块50抑制消除所述回声信号及所述噪声信号所采用的算法为MMSE-LSA算法。所述MMSE-LSA算法为本领域计算人员所公知的用于抑制噪声信号的算法，申请人只是利用回声估计模块40将频域信号中的回声信号的方差计算出来，采用抑制噪声信号的方法将回声信号当成噪声信号抑制消除掉，换句话说，也就是采用抑制噪声信号的方法将语音信号所包括的回声信号及噪声信号一起抑制消除，有效的去除了语音信号中回声信号，提高了语音质量。
进一步地，所述双麦克风回声及噪声的抑制系统还包括与所述短时傅里叶变换模块20和所述抑制模块50连接的短时傅里叶反变换模块60，如图1所示，第二短时傅里叶变换模块20b与短时傅里叶反变换模块60相连，所述短时傅里叶反变换模块60将已抑制消除掉所述回声信号及所述噪声信号的频域信号转化为时域信号输出，进而获得高质量的语音。
本实施例中，请参考图2及图3，其中，图2为本发明一实施例中回声估计模块的工作原理图；图3为本发明中回声估计模块计算回声信号的方差的流程图。如图2所示，所述回声估计模块40计算所述回声信号的方差λe(ω,m)包括如下步骤：
首先，执行步骤S1，对所述频域信号进行平滑处理；
具体的，步骤S1中是采用频域信号，原因是频域信号可以忽略时域中两个麦克风由于位置不同所采集到的信号之间的延时差异，为后续估计噪声及回声能够得到较为精准的数值奠定了基础。本实施例中的所述双麦克风终端上设置有第一麦克风及第二麦克风，所述第一麦克风与所述第二麦克风分别分布在所述双麦克风终端的底部和顶部；所述频域信号包括第一频域信号Y1(ω,m)及第二频域信号Y2(ω,m)，所述第一麦克风采集的语音信号转换为频域信号为所述第一频域信号Y1(ω,m)，所述第二麦克风采集的语音信号转换为所述频域信号为第二频域信号Y2(ω,m)，
Y1(ω,m)＝a1X(ω,m)+D1(ω,m)+b1E(ω,m) (1)
Y2(ω,m)＝a2X(ω,m)+D2(ω,m)+b2E(ω,m) (2)
公式(1)及(2)表征了第一麦克风及第二麦克风的时域信号经短时傅里叶变换模块(即20a及20b)后变为频域信号时的式子，其中，m为帧序号，X(ω,m)为所述频域信号中的近端语音信号，E(ω,m)为所述频域信号中的回声信号，D1(ω,m)为所述第一频域信号Y1(ω,m)中的噪声信号，D2(ω,m)为所述第二频域信号Y2(ω,m)中的噪声信号，a1，a2，b1及b2为由于近场信号的传播引起的幅值差。
较佳的，短时傅里叶变换模块20包括第一短时傅里叶变换模块20a及第二短时傅里叶变换模块20b，即双麦克风终端上设置的第一麦克风及第二麦克风的输出端分别与第一短时傅里叶变换模块20a及第二短时傅里叶变换模块20b相连，对双麦克风终端的两个麦克风所输出的语音信号单独进行频域变换。
其中，所述双麦克风终端10采集的语音信号包括：近端语音信号、噪声信号及回声信号，故所述语音信号转换成频域信号后仍然包括这三种信号，只是变成频域的形式，本申请文件中对于双麦克风回声及噪声的抑制系统涉及计算的过程均在频域下进行的，避免时域时外界的信号干扰，提高了计算回声信号的方差及噪声信号的方差等数值的精准度。
为了较好的理解语音信号所包括的近端语音信号、回声信号及噪声信号，以双麦克风终端10为手机为例进行分析解释。假设将通话时靠近人嘴的麦克风定义为第一麦克风，将设置于手机背面的顶部的麦克风定义为第二麦克风，第一麦克风与第二麦克风之间的垂直距离为10cm左右。近端语音信号的判定标准为：手机在手持模式下，近端信号(本地使用者的语音信号)离第一麦克风大概在3-4cm，近端信号离第二麦克风大概在13cm左右，由此可判定近端语音信号为近场信号。此时第一麦克风采集到的信号比第二麦克风采集到的近端信号大10dB。回声信号的判定标准为：手机在手持模式下，回声离第二麦克风大概2cm左右，离第一麦克风大概在10cm左右，此时第二麦克风采集到的回声比第一麦克风的大10dB。对于噪声信号通常指，指距离手机大于50cm的远场信号(下行链路的信号)，此时第一麦克风与第二麦克风所采集到的远场信号幅值差异较小。
其中，所述对所述频域信号进行平滑处理采用的公式如下：
Φ11(ω,m)＝α1Φ11(ω,m-1)+(1-α1)|Y1(ω,m)|2 (3)
Φ22(ω,m)＝α1Φ22(ω,m-1)+(1-α1)|Y2(ω,m)|2 (4)
公式(3)及(4)中，α1为平滑因子，Φ11(ω,m)为第一频域信号Y1(ω,m)进行平滑处理后的信号,Φ22(ω,m)为第二频域信号Y2(ω,m)进行平滑处理后的信号。
接着，执行步骤S2，利用平滑处理后的频域信号计算比值r(ω,m)；
具体的，所述利用平滑处理后的频域信号计算比值r(ω,m)采用的公式如下：
公式(5)中，r(ω,m)为比值，r1(ω,m)为利用第一频域信号Y1(ω,m)进行平滑处理后的信号Φ11(ω,m)计算的第一比值，r2(ω,m)为利用第二频域信号Y2(ω,m)进行平滑处理后的信号Φ22(ω,m)计算的第二比值。其中，所述第一比值r1(ω,m)及第二比值r2(ω,m)的计算公式如下：
结合前面列举的有关手机的例子、公式(3)至公式(7)可知，当双麦克风终端10采集的语音信号中有回声信号时(也指所述频域信号中有回声信号时)，第一比值r1(ω,m)将变得接近0，第二比值r2(ω,m)变得接近1。当双麦克风终端10采集的语音信号中有近端语音信号时(也指所述频域信号中有近端语音信号时)，第一比值r1(ω,m)和第二比值r2(ω,m)均接近1。当双麦克风终端10采集的语音信号中有噪声信号时(也指所述频域信号中有噪声信号时)，第一比值r1(ω,m)和第二比值r2(ω,m)相近并均处于0～1之间。之后再计算公式(5)可以进一步的区分双麦克风终端10采集的语音信号中包括的信号类型，即在存在回声信号的时候，r(ω,m)将接近1，而存在近端语音信号的时候，r(ω,m)将接近0。理论上，存在噪声信号的时候，r(ω,m)也趋近于0，但在实际信号内，由于噪声信号的随意性，往往和语音时存在较大的台阶。
接着，执行步骤S3，根据所述频域信号及所述比值r(ω,m)计算所述回声信号的方差λe(ω,m)。
进一步地，根据所述频域信号及所述比值r(ω,m)计算所述回声信号的方差λe(ω,m)所采用的公式如下：
公式(8)中β为平滑因子，η为阈值。本文中涉及平滑因子及阈值的都是根据开发者的不断试验所确定的经验值。
另外，考虑到单个频点的信号可能造成误判，则可以进一步采用全局的判决方式来降低这一缺陷，采用的计算公式是将公式(5)演变而成的，具体如下，
此时，公式(8)变为：
此外，还可利用(9)计算的全局比值rg(ω,m)去修正比值r(ω,m)以保证每个频点单独更新，以便控制不同频点的更新步长，采用如下式子实现全局比值rg(ω,m)修正比值r(ω,m)：
公式(11)中，η1及η2均为阈值，阈值η1>η2。
最后将估计到回声信号的方差λe(ω,m)和噪声信号的方差λd(ω,m)直接相加，作为整体的噪声信号的方差进行传统的后处理，当然对于回声信号和噪声信号的结合方式不局限于上述方式，其他方式也可以。
此处需要注意的是，文中涉及角标有g的所表示频域信号的所有频点，不带角标g的表示单独的某个频点，文中若式子中没有写m时，此时表示当前帧，例如r(ω)相当于r(ω,m)。
较佳的，所述噪声估计模块计算所述噪声信号的方差采用的方法为：最小统计量方式、语音激活检测方式或者最小控制递归平均方式。当然计算所述噪声信号的方差所采用的方法包括但不局限于上述这几种，只要能将噪声信号的方差计算出来即可。
进一步地，所述抑制模块包括信噪比估计模块及与所述信噪比模块相连的抑制准则模块；其中，所述信噪比估计模块根据所述回声信号的方差及所述噪声信号的方差计算先验信噪比ξ及后验信噪比γ；所述抑制准则模块通过所述先验信噪比ξ及所述后验信噪比γ计算增益因子G(ξ,γ)；其中，
ξ＝λx(ω,m)/λD(ω,m)；
γ＝|Y1(ω,m)+Y2(ω,m)|2/λD(ω,m)；
λD(ω,m)＝λe(ω,m)+λd(ω,m)；
v＝ξγ/(1+ξ)；
λx(ω,m)为所述语音信号的方差，λd(ω,m)为噪声信号的方差。
较佳的，所述先验信噪比ξ还可以采用DD方式计算。当然计算先验信噪比ξ的方式包括但不局限于这里记载的方式，也可以采用其他能够计算出先验信噪比ξ的方式。
进一步地，所述短时傅里叶变换模块(20a及20b)采用窗函数为Hamming窗或者平顶窗。
本实施例中，所述双麦克风终端10为手机，当然也可以为其它具有双麦克风的终端。
请参考图4，其为使用本发明的双麦克风回声及噪声的抑制系统之前双麦克风终端采集到的声音的语谱图。如图4所示，噪声信号和近端语音信号均使用箭头进行了标识，其它无标识的条纹状信号表征回声信号。从图中可以清楚的观察到，对于近端语音信号，第一麦克风明显比第二麦克风采集的多；对于噪声信号，第一麦克风与第二麦克风采集的相接近；对于回声信号，第二麦克风明显比第一麦克风采集的多。进一步通过图4验证了通过公式(5)、(6)及(7)可以区分出双麦克风终端10所采集信号的类型的方法，符合实际情况。
请参考图5，其为使用本发明的双麦克风回声及噪声的抑制系统之后第一麦克风采集到的声音的语谱图。结合图4及图5所示，可以明显的看出回声信号及噪声信号已被抑制干净，而近端语音信号则无损失，克服了现有技术中消除回声所存在的问题。对于第二麦克风经过本发明的双麦克风回声及噪声的抑制系统之后的语谱图同样可以有回声信号明显被抑制消除的现象，这里就不具体阐述了。
本发明还提供一种双麦克风回声及噪声的抑制方法，包括如下步骤：
短时傅里叶变换模块接收双麦克风终端采集的语音信号并将其转换为频域信号；
回声估计模块及噪声估计模块接收所述频域信号并分别计算所述频域信号中的回声信号的方差及噪声信号的方差；以及
抑制模块根据所述频域信号中的回声信号的方差及噪声信号的方差将所述频域信号中的回声信号及噪声信号抑制消除。
综上，在本发明所提供的双麦克风回声及噪声的抑制系统及其方法中，短时傅里叶变换模块接收双麦克风终端采集的语音信号并将其转换为频域信号，所述回声估计模块及所述噪声估计模块接收所述频域信号并分别计算所述频域信号中的回声信号的方差及噪声信号的方差，所述抑制模块根据所述频域信号中的回声信号的方差及噪声信号的方差将所述频域信号中的回声信号及噪声信号抑制消除。将抑制噪声信号的方式融合于抑制回声信号的方法中，利用抑制模块将回声信号及噪声信号一同抑制消除，整个过程无需采用自适应滤波器，避免现有技术消除回声所存在的问题，提高了语音质量。
上述描述仅是对本发明较佳实施例的描述，并非对本发明范围的任何限定，本发明领域的普通技术人员根据上述揭示内容做的任何变更、修饰，均属于权利要求书的保护范围。
Claims (9)
1.一种双麦克风回声及噪声的抑制系统，包括：短时傅里叶变换模块分别与回声估计模块和噪声估计模块连接，所述回声估计模块和所述噪声估计模块分别与抑制模块连接；所述短时傅里叶变换模块接收双麦克风终端采集的语音信号并将其转换为频域信号，所述回声估计模块及所述噪声估计模块接收所述频域信号并分别计算所述频域信号中的回声信号的方差及噪声信号的方差，所述抑制模块根据所述频域信号中的回声信号的方差及噪声信号的方差将所述频域信号中的回声信号及噪声信号抑制消除；
所述回声估计模块计算所述回声信号的方差包括如下步骤：对所述频域信号进行平滑处理；利用平滑处理后的频域信号计算比值；根据所述频域信号及所述比值计算所述回声信号的方差；
所述双麦克风终端上设置有第一麦克风及第二麦克风，所述第一麦克风与所述第二麦克风分别分布在所述双麦克风终端的底部和顶部；所述频域信号包括第一频域信号Y1(ω,m)及第二频域信号Y2(ω,m)，所述第一麦克风采集的语音信号转换为频域信号为所述第一频域信号Y1(ω,m)，所述第二麦克风采集的语音信号转换为所述频域信号为第二频域信号Y2(ω,m)，其中，
Y1(ω,m)＝a1X(ω,m)+D1(ω,m)+b1E(ω,m)；
Y2(ω,m)＝a2X(ω,m)+D2(ω,m)+b2E(ω,m)；
m为帧序号，X(ω,m)为所述频域信号中的近端语音信号，E(ω,m)为所述频域信号中的回声信号，D1(ω,m)为所述第一频域信号Y1(ω,m)中的噪声信号，D2(ω,m)为所述第二频域信号Y2(ω,m)中的噪声信号，a1，a2，b1及b2为由于近场信号的传播引起的幅值差；
所述对所述频域信号进行平滑处理采用的公式如下：
Φ11(ω,m)＝α1Φ11(ω,m-1)+(1-α1)|Y1(ω,m)|2；
Φ22(ω,m)＝α1Φ22(ω,m-1)+(1-α1)|Y2(ω,m)|2；
其中，α1为平滑因子，Φ11(ω,m)为第一频域信号Y1(ω,m)进行平滑处理后的信号,Φ22(ω,m)为第二频域信号Y2(ω,m)进行平滑处理后的信号；
所述利用平滑处理后的频域信号计算比值采用的公式如下：
其中，r(ω,m)为比值，r1(ω,m)为利用第一频域信号Y1(ω,m)进行平滑处理后的信号Φ11(ω,m)计算的第一比值，r2(ω,m)为利用第二频域信号Y2(ω,m)进行平滑处理后的信号Φ22(ω,m)计算的第二比值；
所述第一比值r1(ω,m)及第二比值r2(ω,m)的计算公式如下：
2.如权利要求1所述的双麦克风回声及噪声的抑制系统，其特征在于，根据所述频域信号及所述比值计算所述回声信号的方差所采用的公式如下：
其中，λe(ω,m)为回声信号的方差，β为平滑因子，η为阈值。
3.如权利要求1所述的双麦克风回声及噪声的抑制系统，其特征在于，所述噪声估计模块计算所述噪声信号的方差采用的方法为：最小统计量方式、语音激活检测方式或者最小控制递归平均方式。
4.如权利要求1所述的双麦克风回声及噪声的抑制系统，其特征在于，所述抑制模块抑制消除所述回声信号及所述噪声信号所采用的算法为MMSE-LSA算法。
5.如权利要求1所述的双麦克风回声及噪声的抑制系统，其特征在于，所述抑制模块包括信噪比估计模块及与所述信噪比模块相连的抑制准则模块；其中，所述信噪比估计模块根据所述回声信号的方差及所述噪声信号的方差计算先验信噪比ξ及后验信噪比γ；所述抑制准则模块通过所述先验信噪比ξ及所述后验信噪比γ计算增益因子G(ξ,γ)；其中，
ξ＝λx(ω,m)/λD(ω,m)；
γ＝|Y1(ω,m)+Y2(ω,m)|2/λD(ω,m)；
λD(ω,m)＝λe(ω,m)+λd(ω,m)；
v＝ξγ/(1+ξ)；
λx(ω,m)为所述语音信号的方差，λd(ω,m)为噪声信号的方差，λe(ω,m)为回声信号的方差。
6.如权利要求1所述的双麦克风回声及噪声的抑制系统，其特征在于，还包括与所述短时傅里叶变换模块和所述抑制模块连接的短时傅里叶反变换模块，所述短时傅里叶反变换模块将已抑制消除掉所述回声信号及所述噪声信号的频域信号转化为时域信号输出。
7.如权利要求1-6中任一项所述的双麦克风回声及噪声的抑制系统，其特征在于，所述短时傅里叶变换模块采用窗函数为Hamming窗或者平顶窗。
8.如权利要求1-6中任一项所述的双麦克风回声及噪声的抑制系统，其特征在于，所述双麦克风终端为手机。
9.一种双麦克风回声及噪声的抑制方法，其特征在于，包括如下步骤：
短时傅里叶变换模块接收双麦克风终端采集的语音信号并将其转换为频域信号；
回声估计模块及噪声估计模块接收所述频域信号并分别计算所述频域信号中的回声信号的方差及噪声信号的方差；以及
抑制模块根据所述频域信号中的回声信号的方差及噪声信号的方差将所述频域信号中的回声信号及噪声信号抑制消除；
其中，所述双麦克风终端上设置有第一麦克风及第二麦克风，所述第一麦克风与所述第二麦克风分别分布在所述双麦克风终端的底部和顶部；所述频域信号包括第一频域信号Y1(ω,m)及第二频域信号Y2(ω,m)，所述第一麦克风采集的语音信号转换为频域信号为所述第一频域信号Y1(ω,m)，所述第二麦克风采集的语音信号转换为所述频域信号为第二频域信号Y2(ω,m)，其中，
Y1(ω,m)＝a1X(ω,m)+D1(ω,m)+b1E(ω,m)；
Y2(ω,m)＝a2X(ω,m)+D2(ω,m)+b2E(ω,m)；
m为帧序号，X(ω,m)为所述频域信号中的近端语音信号，E(ω,m)为所述频域信号中的回声信号，D1(ω,m)为所述第一频域信号Y1(ω,m)中的噪声信号，D2(ω,m)为所述第二频域信号Y2(ω,m)中的噪声信号，a1，a2，b1及b2为由于近场信号的传播引起的幅值差；
所述对所述频域信号进行平滑处理采用的公式如下：
Φ11(ω,m)＝α1Φ11(ω,m-1)+(1-α1)|Y1(ω,m)|2；
Φ22(ω,m)＝α1Φ22(ω,m-1)+(1-α1)|Y2(ω,m)|2；
其中，α1为平滑因子，Φ11(ω,m)为第一频域信号Y1(ω,m)进行平滑处理后的信号,Φ22(ω,m)为第二频域信号Y2(ω,m)进行平滑处理后的信号。
所述利用平滑处理后的频域信号计算比值采用的公式如下：
其中，r(ω,m)为比值，r1(ω,m)为利用第一频域信号Y1(ω,m)进行平滑处理后的信号Φ11(ω,m)计算的第一比值，r2(ω,m)为利用第二频域信号Y2(ω,m)进行平滑处理后的信号Φ22(ω,m)计算的第二比值；
所述第一比值r1(ω,m)及第二比值r2(ω,m)的计算公式如下：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410597191.4A CN105635500B (zh) | 2014-10-29 | 2014-10-29 | 双麦克风回声及噪声的抑制系统及其方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410597191.4A CN105635500B (zh) | 2014-10-29 | 2014-10-29 | 双麦克风回声及噪声的抑制系统及其方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105635500A CN105635500A (zh) | 2016-06-01 |
| CN105635500B true CN105635500B (zh) | 2019-01-25 |
Family
ID=56049961
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410597191.4A Active CN105635500B (zh) | 2014-10-29 | 2014-10-29 | 双麦克风回声及噪声的抑制系统及其方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105635500B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106297818B (zh) * | 2016-09-12 | 2019-09-13 | 广州酷狗计算机科技有限公司 | 一种获取去噪语音信号的方法和装置 |
| CN106302906A (zh) * | 2016-09-15 | 2017-01-04 | 上海斐讯数据通信技术有限公司 | 一种双麦克降噪方法及其装置、一种移动终端 |
| CN108447496B (zh) * | 2018-06-22 | 2020-07-21 | 重庆瑞可利科技有限公司 | 一种基于麦克风阵列的语音增强方法及装置 |
| CN111445918B (zh) * | 2020-03-23 | 2020-12-15 | 深圳市友杰智新科技有限公司 | 降低智能语音设备误唤醒的方法、装置和计算机设备 |
| CN113362843B (zh) * | 2021-06-30 | 2023-02-17 | 北京小米移动软件有限公司 | 音频信号处理方法及装置 |
| CN113489854B (zh) * | 2021-06-30 | 2024-03-01 | 北京小米移动软件有限公司 | 声音处理方法、装置、电子设备和存储介质 |
| CN113823314B (zh) * | 2021-08-12 | 2022-10-28 | 北京荣耀终端有限公司 | 语音处理方法和电子设备 |
| CN113707167A (zh) * | 2021-08-31 | 2021-11-26 | 北京地平线信息技术有限公司 | 残留回声抑制模型的训练方法和训练装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1875611A (zh) * | 2003-11-20 | 2006-12-06 | 摩托罗拉公司(在特拉华州注册的公司) | 用于自适应回声和噪声控制的方法与装置 |
| CN101222555A (zh) * | 2008-01-25 | 2008-07-16 | 上海华平信息技术股份有限公司 | 改善音频通话质量的系统及方法 |
| CN101458943A (zh) * | 2008-12-31 | 2009-06-17 | 北京中星微电子有限公司 | 一种录音控制方法和录音设备 |
| CN101778183A (zh) * | 2009-01-13 | 2010-07-14 | 华为终端有限公司 | 一种残留回声抑制方法及设备 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101719969B (zh) * | 2009-11-26 | 2013-10-02 | 美商威睿电通公司 | 判断双端对话的方法、系统以及消除回声的方法和系统 |
| US9100466B2 (en) * | 2013-05-13 | 2015-08-04 | Intel IP Corporation | Method for processing an audio signal and audio receiving circuit |
-
2014
- 2014-10-29 CN CN201410597191.4A patent/CN105635500B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1875611A (zh) * | 2003-11-20 | 2006-12-06 | 摩托罗拉公司(在特拉华州注册的公司) | 用于自适应回声和噪声控制的方法与装置 |
| CN101222555A (zh) * | 2008-01-25 | 2008-07-16 | 上海华平信息技术股份有限公司 | 改善音频通话质量的系统及方法 |
| CN101458943A (zh) * | 2008-12-31 | 2009-06-17 | 北京中星微电子有限公司 | 一种录音控制方法和录音设备 |
| CN101778183A (zh) * | 2009-01-13 | 2010-07-14 | 华为终端有限公司 | 一种残留回声抑制方法及设备 |
Non-Patent Citations (1)
| Title |
|---|
| 基于最小统计与MMSE-LSA的语音增强;陈立平等;《电声技术》;20090817;第33卷(第08期);第55-59页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105635500A (zh) | 2016-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105635500B (zh) | 双麦克风回声及噪声的抑制系统及其方法 | |
| CN108735213B (zh) | 一种基于相位补偿的语音增强方法及系统 | |
| CN105575397B (zh) | 语音降噪方法及语音采集设备 | |
| JP5952434B2 (ja) | 携帯電話に適用する音声強調方法及び装置 | |
| CN106875938B (zh) | 一种改进的非线性自适应语音端点检测方法 | |
| TWI398855B (zh) | 多重麥克風聲音活動偵測器 | |
| CN102859591B (zh) | 用于语音编码器中的噪声消除的方法和装置 | |
| KR20130108063A (ko) | 다중 마이크로폰의 견고한 잡음 억제 | |
| KR20090127709A (ko) | 적응 빔포밍을 위한 사용자 방향의 소리 검출 기반의 적응모드 제어 장치 및 방법 | |
| CA2458428A1 (en) | System for suppressing wind noise | |
| CN104427143B (zh) | 残留回声检测方法及系统 | |
| CN110211602B (zh) | 智能语音增强通信方法及装置 | |
| CN104754430A (zh) | 终端麦克风降噪装置和方法 | |
| CN102739886A (zh) | 基于回声频谱估计和语音存在概率的立体声回声抵消方法 | |
| CN103730126B (zh) | 噪声抑制方法和噪声抑制器 | |
| CN106448692A (zh) | 应用语音存在概率优化的retf混响消除方法及系统 | |
| WO2016034915A1 (en) | Audio processing circuit and method for reducing noise in an audio signal | |
| CN109068012A (zh) | 一种用于音频会议系统的双端通话检测方法 | |
| CN108010536A (zh) | 回声消除方法、装置、系统及存储介质 | |
| WO2022218254A1 (zh) | 语音信号增强方法、装置及电子设备 | |
| CN103824563A (zh) | 一种基于模块复用的助听器去噪装置和方法 | |
| CN106328160B (zh) | 一种基于双麦克的降噪方法 | |
| CN109817235A (zh) | 一种VoIP设备的回声消除方法 | |
| CN102739286B (zh) | 一种用于通信系统中的回声抵消方法 | |
| Zhang et al. | Noise estimation based on an adaptive smoothing factor for improving speech quality in a dual-microphone noise suppression system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20160601 Assignee: Shanghai Li Ke Semiconductor Technology Co.,Ltd. Assignor: LEADCORE TECHNOLOGY Co.,Ltd. Contract record no.: 2018990000159 Denomination of invention: System and method for inhibiting echo and noise of double microphones License type: Common License Record date: 20180615 |
|
| TA01 | Transfer of patent application right |
Effective date of registration: 20180827 Address after: 201206 China (Shanghai) free trade pilot area, 1258 moon 3, fourth floor, A406 room. Applicant after: Chen core technology Co.,Ltd. Address before: 200233 4 building, No. 333, No. 41, Qinjiang Road, Shanghai, Xuhui District Applicant before: LEADCORE TECHNOLOGY Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20221121 Address after: 266500 No. 687, Binhai street, Huangdao District, Qingdao, Shandong Patentee after: Chenxin Technology Co.,Ltd. Patentee after: Qingdao Weixuan Technology Co.,Ltd. Address before: 201206 China (Shanghai) free trade pilot area, 1258 moon 3, fourth floor, A406 room. Patentee before: Chen core technology Co.,Ltd. |
|
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: Room 102, Building 16, No. 1699, the Pearl River Road, Huangdao District, Qingdao, Shandong 266500 Patentee after: Chenxin Technology Co.,Ltd. Patentee after: Qingdao Weixuan Technology Co.,Ltd. Address before: 266500 No. 687, Binhai street, Huangdao District, Qingdao, Shandong Patentee before: Chenxin Technology Co.,Ltd. Patentee before: Qingdao Weixuan Technology Co.,Ltd. |