CN102855409A - 近场径向干扰源抑制算法及其球麦克风阵列语音采集装置 - Google Patents
近场径向干扰源抑制算法及其球麦克风阵列语音采集装置 Download PDFInfo
- Publication number
- CN102855409A CN102855409A CN2012103501686A CN201210350168A CN102855409A CN 102855409 A CN102855409 A CN 102855409A CN 2012103501686 A CN2012103501686 A CN 2012103501686A CN 201210350168 A CN201210350168 A CN 201210350168A CN 102855409 A CN102855409 A CN 102855409A
- Authority
- CN
- China
- Prior art keywords
- array
- spherical
- microphone
- microphone array
- expression
- 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
- Obtaining Desirable Characteristics In Audible-Bandwidth Transducers (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
一种近场径向干扰源抑制算法及其球麦克风阵列语音采集装置,当期望源位于近场且期望源和干扰源的来波方向和频率相同时利用空间径向抵消滤波算法有效抑制干扰,该球麦克风阵列语音采集装置采用拓扑结构为正十二面体的球麦克风阵列,阵列冗余度小,设计、计算简单,同时可精确确定麦克风在球面上的位置,不会产生位置误差。
Description
技术领域
本发明涉及一种近场径向干扰源抑制算法及其球麦克风阵列语音采集装置。
背景技术
在实际复杂的声学环境下,麦克风拾取的语音信号不可避免地会受到房间混响、噪声以及其他说话人的干扰,即麦克风接收信号为带噪语音,这样不仅影响语音质量,还会使许多语音处理***的性能急剧恶化。因此,需要对麦克风拾取的语音信号进行有效的噪声抑制,减小干扰因素对语音信号的影响,以增强语音信号质量,有效地改善语音处理***的性能。与单麦克风相比,麦克风阵列融合了语音信号的空时信息,具有灵活的波束控制、较高的空间分辨率、较高的信号增益与较强的抗干扰能力等特点,已成为智能通信***中捕捉说话人语音、改善语音质量的重要手段。目前,麦克风阵列技术已广泛应用于车载语音***、语音识别***、助听器和视频会议等各种语音通信领域中。
常见麦克风阵列的拓扑结构主要包括线性阵列、平面阵列和球形阵列。在以上三种阵列形式中,球麦克风阵列(由多个麦克风按照不同的分布情况布置于同一个球面上)物理结构上具有全向旋转对称性,且当空间信号源进行任意方向的旋转移动时,不影响阵列的整体输出,所以能够更好地应用于三维空间波束形成、声源定位、噪声控制以及声场空间特性分析等领域。
实际应用球麦克风阵列进行声源信号采集时,不可能使用大量麦克风覆盖整个球面。因此,需要选择有限数目的麦克风以一定方式对声源信号进行采样。空间采样和时域采样相似,为避免空间信号发生混叠,需要满足有限的带宽(有限的谐波阶数N)。球面空间采样主要有三种形式,分别为等角度采样、高斯采样和均匀采样。其中,等角度采样需要不少于4(N+1)2个采样点,采样点具有规则的角度差异,这种采样需要大量的采样点,实现起来非常困难;高斯采样需要的采样点数至少为2(N+1)2个,同时存在很大的阵列冗余;均匀采样令每两个采样点之间的距离均相等,要求采样点至少为(N+1)2个。与前两种采样方式相比,当麦克风阵列采用均匀采样时,所需采样点麦克风最少,实现起来最方便,也比较广泛。
当采用均匀采样时,如果阵列阶数较小,采样点过少,则阵列***的分辨性能较差;若阵列阶数较大,阵列***的分辨性能虽然有所提高,但由于所需麦克风数目过多,带来设计方面的复杂性,麦克风在球面的位置难以确定,容易产生位置误差,同时对数据采集设备的硬件配置要求过高,难以实现。
在一些近场环境中,麦克风阵列采集空间源声学信号,若期望源与干扰源同向,则普通阵列及其波束形成方法无法抑制掉干扰源、提取期望源,进而实现空间滤波,无法取得较好的语音增强效果。
发明内容
本发明所要解决的技术问题是提供一种当期望源位于近场且期望源和干扰源的来波方向和频率相同时有效抑制干扰的近场径向干扰源抑制算法及其球麦克风阵列语音采集装置,该采集装置阵列冗余度小,设计、计算简单,同时可精确确定麦克风在球面上的位置,不会产生位置误差。
本发明的技术方案是:
一种近场径向干扰源抑制算法,其特殊之处是,具体步骤包括:
第一步,利用球麦克风阵列语音采集装置采集期望源的声压信号,所述球麦克风阵列语音采集装置的麦克风阵列的拓扑结构采用正十二面体,所述正十二面体的每个顶点处均设有一路麦克风,利用基于球阵列的MUSIC算法和基于低频不同模式差近场期望源距离估计算法估计出期望源的位置 ,r s 表示期望源距球心距离,θ s 表示入射俯仰角,φ s 表示入射方位角;
第二步,将球麦克风阵列接收到的声压信号进行球傅里叶变换,声压的球傅里叶变换公式如下:
上式中表示声压的球傅里叶变换系数,为与采样方案相关的权值参数, J表示球阵列所含的阵元数,表示第j个麦克风所采集的声压,为第j个麦克风在球阵表面所处的位置,和分别表示球坐标的俯仰角和方位角,表示n阶m度的球谐函数,为复数共轭运算符;
第三步,根据第一步所估计出的期望源位置,对径向陷波滤波器和多尔夫-切比雪夫径向滤波器进行设计,得到两个波束形成系数d n1(k)和d n2(k),相应滤波器的阵列权值分别为和,阵列权值按照以下公式计算为:
第四步,空间径向抵消滤波器的主支路利用多尔夫-切比雪夫径向滤波器对信号初步处理,抑制一部分噪声和干扰,辅支路用径向陷波滤波器对信号进行处理,完全抑制掉期望信号保留噪声和干扰,两支路输出结果分别为和,分别为:
该球麦克风阵列语音采集装置,包括承载球体、固定在承载球体上的麦克风阵列、支撑承载球体并内置麦克风连接线的支架,其特殊之处是:所述麦克风阵列的拓扑结构采用正十二面体,所述正十二面体的每个顶点处均设有一路麦克风。
所述承载球体的半径r=0.05m~0.5m。
所述支架由插装在一起的上、下支架组成,所述上、下支架的连接处设有控制阀,以自由调节承载球体的高度,方便语音信息采集。
本发明的有益效果是:
1、当期望源位于近场且期望源和干扰源的来波方向和频率相同时,利用空间径向抵消滤波算法有效抑制干扰,得到处理完善的语音信号。
2、由于选用具有20个阵元的球麦克风阵列,根据均匀采样方案可知,此时球麦克风阵列的阶数N可取1、2、3,其阵列冗余度小,计算量小,简单实用,具有良好的分辨性能,所需阵元数与阵列***性能均衡,实用性强。
附图说明
图1是本发明的结构示意图;
图2是本发明的麦克风阵元的分布示意图;
图3是本发明的应用算法的原理图;
图4是本发明的近场径向干扰源抑制算法所生成的径向方向的波束图。
图中:麦克风1、承载球体2、上支架3、控制阀4、下支架5。
具体实施方式
如图1所示,该球麦克风阵列语音采集装置,包括承载球体2、固定在承载球体2上的麦克风阵列、支撑承载球体2并布置麦克风连接线的支架。所述麦克风阵列的拓扑结构采用正十二面体,所述正十二面体的每个顶点处均设有一路麦克风1,即麦克风1数目为M=20,完成接收空间语音信号的任务。承载球体2由钢质材料制成,主要用来安放麦克风1并形成球面的刚性反射。所述支架由插装在一起的上、下支架3、5组成,所述上、下支架3、5的连接处设有控制阀4,以自由调节承载球体2的高度,进而自由调节麦克风阵列的高度,方便语音信息采集。
该球麦克风阵列的设计方法是:
1、根据近场声学环境要求,阵列半径r(承载球体的半径)需满足
r ≈ N / k max ,
式中N为阵列阶数,k = ω/c为来波语音源信号波数,ω为来波语音源角频率,c为声音在空气中的传播速度,k max 为最大波数值,由语音信号频率决定。根据均匀分布采样定理,这里最大N值取3,根据语音信号的频率可知,k max 的取值范围可以为6~60,由此可得承载球体的半径范围为0.05m~0.5m。本实施例中,r取值0.1m。
2、以承载球体2的球心为原点,建立空间直角坐标系,由正十二面体特性计算各顶点的精确位置,各顶点坐标(x\y\z)如下:
其中。
如图3所示,利用了上述球麦克风阵列语音采集装置的近场径向干扰源抑制算法(球麦克风阵列波束形成方法),其特殊之处是,采用空间径向抵消滤波算法。
具体步骤包括:
第一步,利用基于球阵列的MUSIC算法和基于低频不同模式差近场期望源距离估计算法估计出期望源的位置,r s 表示期望源距球心距离,表示入射俯仰角,φs表示入射方位角。本实施例中期望源位于距球心处,入射俯仰角θ s 为 0,入射方位角φ s 为0,则期望源的位置为(0.5,0,0),与期望源同方向的干扰源的空间位置为(1.5, 0, 0)。
第二步,对拓扑结构为正十二面体的球麦克风阵列接收到的声压信号进行球傅里叶变换。声压的球傅里叶变换公式如下:
上式中表示声压的球傅里叶变换系数,为与采样方案相关的权值参数,, J表示球阵列所含的阵元数,表示第j个麦克风所采集的声压,为第j个麦克风在球阵表面所处的位置,和分别表示球坐标的俯仰角和方位角,表示n阶m度的球谐函数,‘’为复数共轭运算符。
第三步,根据第一步所估计的参数,对径向陷波滤波器和多尔夫-切比雪夫径向滤波器进行设计,本实施例中多尔夫-契比雪夫径向滤波器的主旁瓣抑制比选取R=1000,径向陷波滤波器零陷位置选取,得到两个波束形成系数d n1(k)和d n2(k),相应滤波器的阵列权值分别为和,阵列权值所采用的形式为:
第四步,空间径向抵消滤波器的主支路用多尔夫-切比雪夫径向滤波器对信号初步处理,它可以抑制一部分噪声和干扰;辅支路用径向陷波滤波器对信号进行处理,它基本上能够完全抑制掉期望信号保留噪声和干扰,两支路所得的输出结果分别为和,分别为:
。
式中的为利用自适应算法得出的主从支路加权的最优权值,本实施例中确定最优权值为。
利用上述方法得到如图4所示的波束图,如波束图所示,期望源(0.5,0,0)的语音得到最大限度的提取,而与期望源(0.5,0,0)同向的干扰源(1.5, 0, 0)被完全抵消或抑制掉。
表1
Claims (4)
1.一种近场径向干扰源抑制算法,其特征在于,包括如下具体步骤:
第一步,利用球麦克风阵列语音采集装置采集期望源的声压信号,所述球麦克风阵列语音采集装置的麦克风阵列的拓扑结构采用正十二面体,所述正十二面体的每个顶点处均设有一路麦克风,利用基于球阵列的MUSIC算法和基于低频不同模式差近场期望源距离估计算法估计出期望源的位置 ,r s 表示期望源距球心距离,θ s 表示入射俯仰角,φ s 表示入射方位角;
第二步,将球麦克风阵列接收到的声压信号进行球傅里叶变换,声压的球傅里叶变换公式如下:
上式中表示声压的球傅里叶变换系数,为与采样方案相关的权值参数, J表示球阵列所含的阵元数,表示第j个麦克风所采集的声压,为第j个麦克风在球阵表面所处的位置,和分别表示球坐标的俯仰角和方位角,表示n阶m度的球谐函数,为复数共轭运算符;
第三步,根据第一步所估计出的期望源位置,对径向陷波滤波器和多尔夫-切比雪夫径向滤波器进行设计,得到两个波束形成系数d n1(k)和d n2(k),相应滤波器的阵列权值分别为和,阵列权值按照以下公式计算为:
第四步,空间径向抵消滤波器的主支路利用多尔夫-切比雪夫径向滤波器对信号初步处理,抑制一部分噪声和干扰,辅支路用径向陷波滤波器对信号进行处理,完全抑制掉期望信号保留噪声和干扰,两支路输出结果分别为和,分别为:
2.一种用于如权利要求1所述的近场径向干扰源抑制算法的球麦克风阵列语音采集装置,包括承载球体、固定在承载球体上的麦克风阵列、支撑承载球体并内置麦克风连接线的支架,其特征在于:所述麦克风阵列的拓扑结构采用正十二面体,所述正十二面体的每个顶点处均设有一路麦克风。
3.根据权利要求1所述的球麦克风阵列语音采集装置,其特征在于:所述承载球体的半径r=0.05m~0.5m。
4.根据权利要求1所述的球麦克风阵列语音采集装置,其特征在于:所述支架由插装在一起的上、下支架组成,所述上、下支架的连接处设有控制阀。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210350168.6A CN102855409B (zh) | 2012-09-20 | 2012-09-20 | 近场径向干扰源抑制算法及其球麦克风阵列语音采集装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210350168.6A CN102855409B (zh) | 2012-09-20 | 2012-09-20 | 近场径向干扰源抑制算法及其球麦克风阵列语音采集装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102855409A true CN102855409A (zh) | 2013-01-02 |
CN102855409B CN102855409B (zh) | 2014-02-26 |
Family
ID=47401992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210350168.6A Expired - Fee Related CN102855409B (zh) | 2012-09-20 | 2012-09-20 | 近场径向干扰源抑制算法及其球麦克风阵列语音采集装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102855409B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103196546A (zh) * | 2013-01-24 | 2013-07-10 | 中国船舶重工集团公司第七一五研究所 | 一种水声传感器声学参数近场测量方法 |
CN103308889A (zh) * | 2013-05-13 | 2013-09-18 | 辽宁工业大学 | 复杂环境下被动声源二维doa估计方法 |
CN105933820A (zh) * | 2016-04-28 | 2016-09-07 | 冠捷显示科技(中国)有限公司 | 一种外置无线音箱自动定位方法 |
CN107884741A (zh) * | 2017-10-30 | 2018-04-06 | 北京理工大学 | 一种多球阵列多宽带声源快速定向方法 |
JP2018151306A (ja) * | 2017-03-14 | 2018-09-27 | 大成建設株式会社 | 音源探査装置およびマイクロホンアレイ |
CN111627456A (zh) * | 2020-05-13 | 2020-09-04 | 广州国音智能科技有限公司 | 噪音排除方法、装置、设备及可读存储介质 |
CN114878854A (zh) * | 2022-03-03 | 2022-08-09 | 北京航空航天大学 | 一种可用于全流向测量的测速球 |
CN115278449A (zh) * | 2022-09-26 | 2022-11-01 | 中国飞行试验研究院 | 确定麦克风阵列单元坐标的方法、装置、设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9012090D0 (en) * | 1990-05-31 | 1990-07-18 | Wealthy Peter B | Telephone for the disabled |
CN101548554A (zh) * | 2006-10-06 | 2009-09-30 | 彼得·G·克拉文 | 麦克风阵列 |
CN101674508A (zh) * | 2009-09-27 | 2010-03-17 | 上海大学 | 三经线交点球麦克风阵列及其设计方法 |
-
2012
- 2012-09-20 CN CN201210350168.6A patent/CN102855409B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9012090D0 (en) * | 1990-05-31 | 1990-07-18 | Wealthy Peter B | Telephone for the disabled |
CN101548554A (zh) * | 2006-10-06 | 2009-09-30 | 彼得·G·克拉文 | 麦克风阵列 |
CN101674508A (zh) * | 2009-09-27 | 2010-03-17 | 上海大学 | 三经线交点球麦克风阵列及其设计方法 |
Non-Patent Citations (1)
Title |
---|
吴涛,王冬霞,郭庆龙,程宗明: "基于球麦克风阵列的语音增加方法", 《辽宁工业大学学报(自然科学报)》, no. 31, 31 December 2011 (2011-12-31), pages 351 - 357 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103196546A (zh) * | 2013-01-24 | 2013-07-10 | 中国船舶重工集团公司第七一五研究所 | 一种水声传感器声学参数近场测量方法 |
CN103308889A (zh) * | 2013-05-13 | 2013-09-18 | 辽宁工业大学 | 复杂环境下被动声源二维doa估计方法 |
CN103308889B (zh) * | 2013-05-13 | 2014-07-02 | 辽宁工业大学 | 复杂环境下被动声源二维doa估计方法 |
CN105933820A (zh) * | 2016-04-28 | 2016-09-07 | 冠捷显示科技(中国)有限公司 | 一种外置无线音箱自动定位方法 |
JP2018151306A (ja) * | 2017-03-14 | 2018-09-27 | 大成建設株式会社 | 音源探査装置およびマイクロホンアレイ |
CN107884741A (zh) * | 2017-10-30 | 2018-04-06 | 北京理工大学 | 一种多球阵列多宽带声源快速定向方法 |
CN107884741B (zh) * | 2017-10-30 | 2021-01-19 | 北京理工大学 | 一种多球阵列多宽带声源快速定向方法 |
CN111627456A (zh) * | 2020-05-13 | 2020-09-04 | 广州国音智能科技有限公司 | 噪音排除方法、装置、设备及可读存储介质 |
CN114878854A (zh) * | 2022-03-03 | 2022-08-09 | 北京航空航天大学 | 一种可用于全流向测量的测速球 |
CN115278449A (zh) * | 2022-09-26 | 2022-11-01 | 中国飞行试验研究院 | 确定麦克风阵列单元坐标的方法、装置、设备及存储介质 |
CN115278449B (zh) * | 2022-09-26 | 2023-03-10 | 中国飞行试验研究院 | 确定麦克风阵列单元坐标的方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102855409B (zh) | 2014-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102855409B (zh) | 近场径向干扰源抑制算法及其球麦克风阵列语音采集装置 | |
CN102804809B (zh) | 音频源定位 | |
US9291697B2 (en) | Systems, methods, and apparatus for spatially directive filtering | |
AU2011322560B2 (en) | Apparatus and method for deriving a directional information and computer program product | |
CN108549052B (zh) | 一种时频-空域联合加权的圆谐域伪声强声源定位方法 | |
CN109932689A (zh) | 一种适用于特定定位场景的任意阵列优化方法 | |
CN112180329B (zh) | 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法 | |
CN109991567B (zh) | 一种水下滑翔机四面体阵三维被动测向方法 | |
CN102147458B (zh) | 一种针对宽带声源的波达方向估计方法及其装置 | |
CN107167770B (zh) | 一种混响条件下的麦克风阵列声源定位装置 | |
CN106371057B (zh) | 语音声源测向方法及装置 | |
CN206057554U (zh) | 一种多球麦克风阵列声场声压采集装置 | |
CN206114888U (zh) | 语音声源测向装置 | |
CN202799044U (zh) | 球麦克风阵列语音采集装置 | |
Nannuru et al. | Sparse planar arrays for azimuth and elevation using experimental data | |
CN108595758A (zh) | 一种任意形式传感器阵列最优宽带波束图综合方法 | |
CN104008287B (zh) | 基于pso‑mvdr的声场重构与鬼影抑制方法 | |
CN115061089B (zh) | 一种声源定位方法、***、介质、设备及装置 | |
CN111060867A (zh) | 一种指向性麦克风微阵列波达方向估计方法 | |
CN115547354A (zh) | 波束形成方法、装置及设备 | |
CN111157952B (zh) | 一种基于移动麦克风阵列的房间边界估计方法 | |
JP6821836B2 (ja) | 収音装置 | |
CN110876100B (zh) | 一种音源定向方法与*** | |
CN112666521A (zh) | 一种基于改进的自适应陷波器的室内声源定位方法 | |
Li et al. | Conformal cylindrical array sound source localization at the presence of shadowed elements |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140226 Termination date: 20140920 |
|
EXPY | Termination of patent right or utility model |