CN103165137A - 一种非平稳噪声环境下传声器阵列的语音增强方法 - Google Patents
一种非平稳噪声环境下传声器阵列的语音增强方法 Download PDFInfo
- Publication number
- CN103165137A CN103165137A CN2011104273292A CN201110427329A CN103165137A CN 103165137 A CN103165137 A CN 103165137A CN 2011104273292 A CN2011104273292 A CN 2011104273292A CN 201110427329 A CN201110427329 A CN 201110427329A CN 103165137 A CN103165137 A CN 103165137A
- Authority
- CN
- China
- Prior art keywords
- omega
- signal
- microphone
- probability
- noise
- 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
Images
Landscapes
- Circuit For Audible Band Transducer (AREA)
Abstract
本发明涉及一种在非平稳噪声环境下传声器阵列的语音增强方法,通过利用通道之间的相位信息,准确判断有效信号的概率,以提高抗非平稳噪声干扰的性能。该方法估计出有效信号存在概率具体包括如下步骤:(51)根据所需波束宽度计算门限T1=NL sin(θ)/Fs,其中N表示帧长,L表示传声器单元间距,θ表示波束宽度,Fs为采样频率;(52)计算各通道输入信号在各频率点相位差ω的和值pd(n,ω);(53)当pd(n,ω)>ω(M-1)T1时,该频率点存在有效信号概率p1(n,ω)=1,否则p1(n,ω)=0,其中,M表示传声器个数。本发明由于利用了通道之间的相位信息,使得对有效信号的概率判断更加准确,提高了拾取语音的信噪比,提高了整体性能,特别是抗非平稳噪声干扰性能。
Description
技术领域
本发明涉及传声器阵列语音增强领域,特别涉及一种在非平稳噪声环境下传声器阵列的语音增强方法。
背景技术
在视频会议系统、车载免提系统、语音识别系统等应用场合,通常需要在较远距离拾取目标语音。然而,由于背景噪声、干扰噪声及回声的存在,拾取的语音信噪比通常较低,而基于传声器阵列的语音增强技术是提高接收语音信噪比的常用手段之一。
目前,传声器阵列语音增强技术主要有两类,一类是基于常规波束形成的方法,这种方法的波束形状是固定的,不会随着干扰声信号的变化而变化,其对非方向性噪声如混响声有较好的抑制效果,语音增强性能随着传声器数量增多而提高,但是,对方向性干扰噪声,其抑制效果有限;另一类是基于自适应波束形成算法,其波束形状会随着干扰语音的变化而发生改变,自适应波束形成算法能够有效的抑制方向性干扰噪声但在弱相关噪声场或在散射噪声场中,性能不如传统波束形成算法。上述两种算法通常均需要干扰噪声是平稳过程,然而在实际使用环境下,干扰噪声通常是突发非平稳的,如语音干扰,在这种情况下,目前大部分的语音增强技术性能有限。
发明内容
本发明的目的在于,针对现有技术中存在的不足局限性,从而提供一种在非平稳噪声环境下传声器阵列的语音增强方法。
为实现上述发明目的,本发明提供的一种在非平稳噪声环境下传声器阵列的语音增强方法,通过利用通道之间的相位信息,准确判断有效信号的概率,以提高抗非平稳噪声干扰的性能。
该方法的步骤包括:
(1)根据传声器阵列的输入信号估计出噪声参考信号r(n);
(2)采用自适应算法初步滤除各通道的输入信号中所含的噪声信号;
(3)对各通道中自适应滤波器的输出信号进行分帧加窗,再采用快速傅立叶变换变换至频域;
(4)估计背景噪声功率谱n(ω),并由此得到当前帧的信噪比SNR(n,ω);
(5)利用各通道之间的相位信息,再进一步根据和值与预定义门限T1的关系得到有效语音信号存在概率p1(n,ω);
(6)估计出当前帧是有效语音帧的概率pf(n);
(7)计算出输出信号在各个频率点增益G(n,ω);
(8)把处理后的频谱信号通过反傅里叶变换变换到时间域,作为输出信号;
所述步骤(5)中的估计出有效信号存在概率具体包括如下步骤:
(51)根据所需波束宽度计算门限T1=NL sin(θ)/Fs,其中N表示帧长,L表示传声器单元间距,θ表示波束宽度,Fs为采样频率;
(52)计算各通道输入信号在各频率点相位差ω的和值pd(n,ω);
(53)当pd(n,ω)>ω(M-1)T1时,该频率点存在有效信号概率p1(n,ω)=1,否则p1(n,ω)=0,其中,M表示传声器个数。
作为上述技术方案的一种改进,所述的步骤(6)中的估计当前帧是有效语音帧的概率具体包括如下步骤:
(62)计算满足pd(n,ω)>ω(M-1)T1频率点个数C;
作为上述技术方案的一种选择,所述的传声器阵列为两个传声器单元,该方法的具体步骤包括:
(1)根据两个传声器输入信号计算参考背景噪声信号:
r(n)=x1(n)-x2(n)
其中,xi(n)表示第i个传声器在时刻n的输入信号;
(2)以背景噪声信号作为参考信号对传声器信号进行自适应滤波初步去除传声器输入信号中的噪声,所述的自适应滤波算法采用NLMS算法,具体描述如下:
wi(k,n)=wi(k,n-1)-μ(n)ei(n)n(n) k=0...K-1;
式中,K=20为自适应滤波器阶数,wi(k,n)为自适应滤波器系数,μ为归一化收敛系数,μ(n)为收敛系数。
(4)采用MCRA方法估计出背景噪声功率谱n(ω),并由此得到信号信噪比:
(5)根据两个传声器单元在各频率点的相位差:
其中,∠表示求取复数相位,
再进一步根据和值与预定义门限T1的关系得到有效语音信号存在概率:
其中,预定义门限T1为:
T1=NL sin(θ)/Fs,
式中,L表示传声器单元间距,θ表示波束宽度,Fs为采样频率;
(6)首先,计算满足pd(n,ω)>ω(M-1)T1频率点个数C;并计算当前帧总体能量:
最后,得到当前帧为有效语音帧概率:
(7)计算总体增益:
增强后信号为:z(n,ω)=G(n,ω)y(n,ω);
(8)通过反傅里叶变换转为时域信号。
本发明的优点在于,本发明的传声器阵列语音增强具有如下优点:相比于其他方法,由于利用了通道之间的相位信息,使得对有效信号的概率判断更加准确,提高了拾取语音的信噪比,提高了整体性能,特别是抗非平稳噪声干扰性能。
附图说明
图1是为本发明传声器阵列的语音增强方法的流程示意图。
具体实施方式
以一个具体实施例对本发明的传声器阵列语音增强方法进行说明。实施例中传声器阵列单元个数为2,间距为8cm,波束宽度设为10度,采样频率为16000Hz。
如图1所示,首先,根据两个传声器输入信号计算参考背景噪声信号:
r(n)=x1(n)-x2(n);
其中,xi(n)表示第i个传声器在时刻n的输入信号。以背景噪声信号作为参考信号对传声器信号进行自适应滤波,初步去除传声器输入信号中的噪声,自适应滤波算法可采用常用的NLMS算法,具体描述如下:
wi(k,n)=wi(k,n-1)-μ(n)ei(n)n(n) k=0...K-1
采用Israel Cohen提出的MCRA方法估计出背景噪声率谱n(ω),从而信号信噪比
计算两个传声器单元在各频率点的相位差
其中∠表示求取复数相位。根据和值与预定义门限T1的关系得到有效语音信号存在概率
预定义门限T1通常选取为
T1=NL sin(θ)/Fs,
其中L=0.08表示传声器单元间距,θ=10°表示波束宽度,Fs=16000为采样频率。
计算满足pd(n,ω)>ω(M-1)T1频率点个数C;计算当前帧总体能量
计算当前帧为有效语音帧概率
计算总体增益
增强后信号
z(n,ω)=G(n,ω)y(n,ω)
最后通过反傅里叶变换转为时域信号。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种非平稳噪声环境下传声器阵列的语音增强方法,该方法包括的步骤如下:
(1)根据传声器阵列的输入信号估计出噪声参考信号r(n);
(2)采用自适应算法初步滤除各通道的输入信号中所含的噪声信号;
(3)对各通道中自适应滤波器的输出信号进行分帧加窗,再采用快速傅立叶变换变换至频域;
(4)估计背景噪声功率谱n(ω),并由此得到当前帧的信噪比SNR(n,ω);
(5)利用各通道之间的相位信息,再进一步根据和值与预定义门限T1的关系得到有效语音信号存在概率p1(n,ω);
(6)估计出当前帧是有效语音帧的概率pf(n);
(7)计算出输出信号在各个频率点增益G(n,ω);
(8)把处理后的频谱信号通过反傅里叶变换变换到时间域,作为输出信号;
所述步骤(5)中的估计出有效信号存在概率具体包括如下步骤:
(51)根据所需波束宽度计算门限T1=NL sin(θ)/Fs,其中N表示帧长,L表示传声器单元间距,θ表示波束宽度,Fs为采样频率;
(52)计算各通道输入信号在各频率点相位差ω的和值pd(n,ω);
(53)当pd(n,ω)>ω(M-1)T1时,该频率点存在有效信号概率p1(n,ω)=1,否则p1(n,ω)=0,其中,M表示传声器个数。
3.根据权利要求1或2所述的非平稳噪声环境下传声器阵列的语音增强方法,其特征在于,所述的传声器阵列为两个传声器单元,该方法的具体步骤包括:
(1)根据两个传声器输入信号计算参考背景噪声信号:
r(n)=x1(n)-x2(n);
其中,xi(n)表示第i个传声器在时刻n的输入信号;
(2)以背景噪声信号作为参考信号对传声器信号进行自适应滤波初步去除传声器输入信号中的噪声,所述的自适应滤波算法采用NLMS算法,具体描述如下:
wi(k,n)=wi(k,n-1)-μ(n)ei(n)n(n) k=0...K-1;
式中,K=20为自适应滤波器阶数,wi(k,n)为自适应滤波器系数;μ为归一化收敛系数,μ(n)为收敛系数;
(4)采用MCRA方法估计出背景噪声功率谱n(ω),并由此得到信号信噪比:
(5)根据两个传声器单元在各频率点的相位差:
其中,∠表示求取复数相位,
再进一步根据和值与预定义门限T1的关系得到有效语音信号存在概率:
其中,预定义门限T1为:
T1=NL sin(θ)/Fs,
式中,L表示传声器单元间距,θ表示波束宽度,Fs为采样频率;
(6)首先,计算满足pd(n,ω)>ω(M-1)T1频率点个数C;并计算当前帧总体能量:
最后,得到当前帧为有效语音帧概率:
(7)计算总体增益:
增强后信号为:z(n,ω)=G(n,ω)y(n,ω);
(8)通过反傅里叶变换转为时域信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110427329.2A CN103165137B (zh) | 2011-12-19 | 2011-12-19 | 一种非平稳噪声环境下传声器阵列的语音增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110427329.2A CN103165137B (zh) | 2011-12-19 | 2011-12-19 | 一种非平稳噪声环境下传声器阵列的语音增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103165137A true CN103165137A (zh) | 2013-06-19 |
CN103165137B CN103165137B (zh) | 2015-05-06 |
Family
ID=48588161
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110427329.2A Active CN103165137B (zh) | 2011-12-19 | 2011-12-19 | 一种非平稳噪声环境下传声器阵列的语音增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103165137B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105225673A (zh) * | 2014-06-09 | 2016-01-06 | 杜比实验室特许公司 | 噪声水平估计 |
CN105321528A (zh) * | 2014-06-27 | 2016-02-10 | 中兴通讯股份有限公司 | 一种麦克风阵列语音检测方法及装置 |
CN105590630A (zh) * | 2016-02-18 | 2016-05-18 | 南京奇音石信息技术有限公司 | 基于指定带宽的定向噪音抑制方法 |
CN105635453A (zh) * | 2015-12-28 | 2016-06-01 | 上海博泰悦臻网络技术服务有限公司 | 一种通话音量自动调节方法、系统、车载设备及汽车 |
CN106024018A (zh) * | 2015-03-27 | 2016-10-12 | 大陆汽车系统公司 | 实时风冲击噪声检测 |
CN106052852A (zh) * | 2016-06-01 | 2016-10-26 | 中国电子科技集团公司第三研究所 | 一种脉冲声信号的检测方法及装置 |
CN107248412A (zh) * | 2017-05-18 | 2017-10-13 | 歌尔科技有限公司 | 音频降噪方法和装置 |
CN107369460A (zh) * | 2017-07-31 | 2017-11-21 | 深圳海岸语音技术有限公司 | 基于声学矢量传感器空间锐化技术的语音增强装置及方法 |
CN108831499A (zh) * | 2018-05-25 | 2018-11-16 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 利用语音存在概率的语音增强方法 |
CN109150238A (zh) * | 2018-10-25 | 2019-01-04 | 深圳星禾设计有限公司 | 无线音频一体化测试系统和方法 |
CN109473118A (zh) * | 2018-12-24 | 2019-03-15 | 苏州思必驰信息科技有限公司 | 双通道语音增强方法及装置 |
CN111933165A (zh) * | 2020-07-30 | 2020-11-13 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 突变噪声快速估计方法 |
CN113689875A (zh) * | 2021-08-25 | 2021-11-23 | 湖南芯海聆半导体有限公司 | 一种面向数字助听器的双麦克风语音增强方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079267A (zh) * | 2006-05-26 | 2007-11-28 | 富士通株式会社 | 定向集音装置、定向集音方法以及存储器产品 |
CN101236250A (zh) * | 2007-01-30 | 2008-08-06 | 富士通株式会社 | 声音判定方法和声音判定装置 |
US20080235013A1 (en) * | 2007-03-22 | 2008-09-25 | Samsung Electronics Co., Ltd. | Method and apparatus for estimating noise by using harmonics of voice signal |
WO2010130928A1 (fr) * | 2009-05-14 | 2010-11-18 | Parrot | Procede de selection d'un microphone parmi deux microphones ou plus, pour un systeme de traitement de la parole tel qu'un dispositif telephonique "mains libres" operant dans un environnement bruite |
US20110022361A1 (en) * | 2009-07-22 | 2011-01-27 | Toshiyuki Sekiya | Sound processing device, sound processing method, and program |
-
2011
- 2011-12-19 CN CN201110427329.2A patent/CN103165137B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079267A (zh) * | 2006-05-26 | 2007-11-28 | 富士通株式会社 | 定向集音装置、定向集音方法以及存储器产品 |
CN101236250A (zh) * | 2007-01-30 | 2008-08-06 | 富士通株式会社 | 声音判定方法和声音判定装置 |
US20080235013A1 (en) * | 2007-03-22 | 2008-09-25 | Samsung Electronics Co., Ltd. | Method and apparatus for estimating noise by using harmonics of voice signal |
WO2010130928A1 (fr) * | 2009-05-14 | 2010-11-18 | Parrot | Procede de selection d'un microphone parmi deux microphones ou plus, pour un systeme de traitement de la parole tel qu'un dispositif telephonique "mains libres" operant dans un environnement bruite |
US20110022361A1 (en) * | 2009-07-22 | 2011-01-27 | Toshiyuki Sekiya | Sound processing device, sound processing method, and program |
Non-Patent Citations (1)
Title |
---|
袁榕嵘; 吴鸣; 杨军;: "结合相位差和谱减法的双传声器语音增强技术", 《电声技术》, 31 May 2011 (2011-05-31), pages 50 - 54 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105225673A (zh) * | 2014-06-09 | 2016-01-06 | 杜比实验室特许公司 | 噪声水平估计 |
CN105225673B (zh) * | 2014-06-09 | 2020-12-04 | 杜比实验室特许公司 | 用于噪声水平估计的方法、系统和介质 |
CN105321528A (zh) * | 2014-06-27 | 2016-02-10 | 中兴通讯股份有限公司 | 一种麦克风阵列语音检测方法及装置 |
CN105321528B (zh) * | 2014-06-27 | 2019-11-05 | 中兴通讯股份有限公司 | 一种麦克风阵列语音检测方法及装置 |
CN106024018A (zh) * | 2015-03-27 | 2016-10-12 | 大陆汽车系统公司 | 实时风冲击噪声检测 |
CN106024018B (zh) * | 2015-03-27 | 2022-06-03 | 大陆汽车系统公司 | 实时风冲击噪声检测 |
CN105635453A (zh) * | 2015-12-28 | 2016-06-01 | 上海博泰悦臻网络技术服务有限公司 | 一种通话音量自动调节方法、系统、车载设备及汽车 |
CN105590630A (zh) * | 2016-02-18 | 2016-05-18 | 南京奇音石信息技术有限公司 | 基于指定带宽的定向噪音抑制方法 |
CN105590630B (zh) * | 2016-02-18 | 2019-06-07 | 深圳永顺智信息科技有限公司 | 基于指定带宽的定向噪音抑制方法 |
CN106052852B (zh) * | 2016-06-01 | 2019-03-08 | 中国电子科技集团公司第三研究所 | 一种脉冲声信号的检测方法及装置 |
CN106052852A (zh) * | 2016-06-01 | 2016-10-26 | 中国电子科技集团公司第三研究所 | 一种脉冲声信号的检测方法及装置 |
CN107248412A (zh) * | 2017-05-18 | 2017-10-13 | 歌尔科技有限公司 | 音频降噪方法和装置 |
CN107369460A (zh) * | 2017-07-31 | 2017-11-21 | 深圳海岸语音技术有限公司 | 基于声学矢量传感器空间锐化技术的语音增强装置及方法 |
CN107369460B (zh) * | 2017-07-31 | 2020-08-21 | 深圳海岸语音技术有限公司 | 基于声学矢量传感器空间锐化技术的语音增强装置及方法 |
CN108831499A (zh) * | 2018-05-25 | 2018-11-16 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 利用语音存在概率的语音增强方法 |
CN108831499B (zh) * | 2018-05-25 | 2020-07-21 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 利用语音存在概率的语音增强方法 |
CN109150238A (zh) * | 2018-10-25 | 2019-01-04 | 深圳星禾设计有限公司 | 无线音频一体化测试系统和方法 |
CN109473118A (zh) * | 2018-12-24 | 2019-03-15 | 苏州思必驰信息科技有限公司 | 双通道语音增强方法及装置 |
CN109473118B (zh) * | 2018-12-24 | 2021-07-20 | 思必驰科技股份有限公司 | 双通道语音增强方法及装置 |
CN111933165A (zh) * | 2020-07-30 | 2020-11-13 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 突变噪声快速估计方法 |
CN113689875A (zh) * | 2021-08-25 | 2021-11-23 | 湖南芯海聆半导体有限公司 | 一种面向数字助听器的双麦克风语音增强方法和装置 |
CN113689875B (zh) * | 2021-08-25 | 2024-02-06 | 湖南芯海聆半导体有限公司 | 一种面向数字助听器的双麦克风语音增强方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103165137B (zh) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103165137B (zh) | 一种非平稳噪声环境下传声器阵列的语音增强方法 | |
CN101510426B (zh) | 一种噪声消除方法及系统 | |
CN102855880B (zh) | 用于多话筒音频设备的去噪方法 | |
EP3047483B1 (en) | Adaptive phase difference based noise reduction for automatic speech recognition (asr) | |
CN101192411B (zh) | 大距离麦克风阵列噪声消除的方法和噪声消除系统 | |
CN102498709B (zh) | 从两个或多个麦克风中为语音处理系统,例如工作在嘈杂环境中的免提电话设备,选择一个麦克风的方法 | |
CN101369427B (zh) | 用于音频信号处理的方法和装置 | |
US8891785B2 (en) | Processing signals | |
CN105244036A (zh) | 一种麦克风语音增强方法及装置 | |
CN101325061A (zh) | 用于音频信号处理的方法和装置 | |
Niwa et al. | Post-filter design for speech enhancement in various noisy environments | |
JP6225245B2 (ja) | 信号処理装置、方法及びプログラム | |
CN108538306B (zh) | 提高语音设备doa估计的方法及装置 | |
CN104835503A (zh) | 一种改进gsc自适应语音增强方法 | |
CN110085247A (zh) | 一种针对复杂噪声环境的双麦克风降噪方法 | |
KR100940629B1 (ko) | 잡음 제거 장치 및 방법 | |
CN104575513A (zh) | 突发噪声的处理系统、突发噪声的检测及抑制方法与装置 | |
CN105353340A (zh) | 一种双层圆柱阵水下被动目标检测方法 | |
CN106782618B (zh) | 基于二阶锥规划的目标方向语音检测方法 | |
CN110992977B (zh) | 一种目标声源的提取方法及装置 | |
CN109243476A (zh) | 混响语音信号中后混响功率谱的自适应估计方法及装置 | |
KR20090098552A (ko) | 위상정보를 이용한 자동 이득 조절 장치 및 방법 | |
CN111210836B (zh) | 一种麦克风阵列波束形成动态调整方法 | |
Chen et al. | Filtering techniques for noise reduction and speech enhancement | |
Panda | A fast approach to psychoacoustic model compensation for robust speaker recognition in additive noise. |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20180313 Address after: 215000 401 Room 401, No. 158, Jinfeng Road, science and Technology City, high tech Zone, Jiangsu Province Patentee after: In the new science and Technology Co. Ltd. (Suzhou) Address before: 100190 Beijing City, North Fourth Ring Road West, No. 21, No. Patentee before: Institute of acoustics, Chinese Academy of Sciences |
|
TR01 | Transfer of patent right |