CN106842111B - 基于传声器镜像的室内声源定位方法 - Google Patents
基于传声器镜像的室内声源定位方法 Download PDFInfo
- Publication number
- CN106842111B CN106842111B CN201611230608.9A CN201611230608A CN106842111B CN 106842111 B CN106842111 B CN 106842111B CN 201611230608 A CN201611230608 A CN 201611230608A CN 106842111 B CN106842111 B CN 106842111B
- Authority
- CN
- China
- Prior art keywords
- formula
- mirror image
- microphone
- frequency
- room
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
- G01S1/72—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using ultrasonic, sonic or infrasonic waves
- G01S1/76—Systems for determining direction or position line
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
本发明公开了一种基于传声器镜像的室内声源定位方法,用于解决现有室内声源定位方法复杂的技术问题。技术方案是以室内多传声器的采集数据为输入,将室内各接收点位置、房间尺寸信息与虚声源法结合建立信号模型,并借助传声器镜像构建求解字典,最后采取Lasso优化算法进行求解,求出的位置向量可直接获取声源的位置信息;在定位计算时,使用基于传声器镜像的字典构建字典,可以显著提升算法计算效率,提高室内声源定位的实时效率,方法简单易行。
Description
技术领域
本发明涉及一种室内声源定位方法,特别涉及一种基于传声器镜像的室内声源定位方法。
背景技术
文献“Sparse sound field decomposition using group sparse Bayesianlearning,in 2015Asia-Pacific Signal and Information Processing AssociationAnnual Summit and Conference(APSIPA),2015,pp.850-855.”公开了一种室内声源定位方法。该方法将直达声和混响声分解为单独的声源信息以及多个虚源的混响信息,采集多个接收信号,在室内和室外空间划分网格,构建字典求解来进行室内声源定位。文献所述方法基于虚声源法,将室内的混响信息分解为由声源以及房间外多个虚源辐射声波的叠加,基于这些信息可将房间内部的直达声和混响声分解。然而,当室内混响较强时,为满足声能量衰减为10%以下的要求,空间内设定网格点数目会极大地增加,这将导致求解时所构造字典的相应的维数产生相应的扩张。而稀疏求解时,字典原子个数的增多会增加求解的计算成本,降低计算效率,这将严重影响声源定位的实时输出。
发明内容
为了克服现有室内声源定位方法复杂的不足,本发明提供一种基于传声器镜像的室内声源定位方法。该方法以室内多传声器的采集数据为输入,将室内各接收点位置、房间尺寸信息与虚声源法结合建立信号模型,并借助传声器镜像构建求解字典,最后采取Lasso优化算法进行求解,求出的位置向量可直接获取声源的位置信息;在定位计算时,使用基于传声器镜像的字典构建字典,可以显著提升算法计算效率,提高室内声源定位的实时效率,方法简单易行。
本发明解决其技术问题所采用的技术方案:一种基于传声器镜像的室内声源定位方法,其特点是包括以下步骤:
步骤一、在室内布置多个传声器,所述传声器为无指向性型。对于1000~2000频率段,传声器与声源的距离大于0.2m。采集信号工作实时同步完成,数据预处理的采样频率为16kHz以上,录制时间不超过30s。
步骤二、对于尺寸为lx×ly×lz的房间,任意接收点镜像的空间坐标表示为
xi=2l·lx±xm
yi=2m·ly±ym (1)
zi=2n·lz±zm
式中,接收点位置为rm=(xm,ym,zm)),某级接收点镜像的位置为ri=(xi,yi,zi),l、m、n是界于-∞和+∞之间的整数,镜像对应的反射级数由(2)式给出
N=|2l-Δl|+|2m-Δm|+|2n-Δn| (2)
Δl=1表示求xi式中xm前取负;Δl=0表示求xi式中xm前取正;Δm=1表示求yi式中ym前取负;Δm=0表示求yi式中ym前取正;Δn=1表示求zi式中zm前取负;Δn=0表示求zi式中zm前取正;N=0时表示传声器本身,所以无反射级数。
步骤三、对于房间内某一二维平面,平均吸声系数大于0.1,房间内部空间沿各方向按照距离间隔Δx、Δy和Δz划分网格,每个网格点gj的位置信息取网格点中心位置,设总网格点数为G。
步骤四、计算各传声器位置及其一定阶数镜像的位置rni与每个网格点位置gj之间的距离|rni-gj|,将所述距离|rni-gj|与对应的镜像级数和频率联合运算,求解字典矩阵Df中的各元素,Df的维度为M×G。
式中,f为频率,c为空气中的声速,N0为传声器镜像的最高阶数。
步骤五、将传声器信号的频域表示yf和各频率的字典矩阵Df作为输入,采取Lasso优化算法求解xf。
yf=Dfxf (5)
步骤六、基于频域表示yf和字典矩阵Df,根据公式(5),选择l1范数作为求解方法获取各单独频率下以及联合频率下的位置向量xf和x,二者具有相同的维度。式中,y∈CM×1,x,xf∈CG×1,Df∈CM×G。在进行室内声定位时,房间内的网格点数目N远远大于 传声器个数M。相对于整个房间,x具有稀疏性,给x加稀疏约束求解方程。
步骤七、选取基于Lasso优化算法,构建如下的单频率及多频率的目标函数:
其中,λ是需要调节的参数,用来平衡估计误差和信号稀疏度。
步骤八、根据式(6)和式(7)求解结果xf和x,其最大值对应的网格点位置为该频率点下定位得到的声源位置。
本发明的有益效果是:该方法以室内多传声器的采集数据为输入,将室内各接收点位置、房间尺寸信息与虚声源法结合建立信号模型,并借助传声器镜像构建求解字典,最后采取Lasso优化算法进行求解,求出的位置向量可直接获取声源的位置信息;在定位计算时,使用基于传声器镜像的字典构建字典,可以显著提升算法计算效率,提高室内声源定位的实时效率,方法简单易行。
假设室内划分为a2个网格点,则对于背景技术方法字典对应的网格点数目为(2n+1)2*a2个,n为与房间混响相关的虚源层数,房间混响越强烈,层数n数值也越大。而对于本发明方法,在对房间混响或平均吸声有预估的基础上,最终所构造字典对应的网格点数目为a2个,且网格点数目不会因房间混响的强弱而改变。假设n=2,房间内部网格点为10×10,背景技术方法计算时的网格点为2500,求解需要1小时,然而本发明方法在预先不考虑接收点虚源计算时间时,只需要100个网格点,求解只需要2.4分钟。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明基于传声器镜像的室内声源定位方法的流程图。
图2是本发明方法用测试系统连接框图。
图3是不同信噪比下的定位对比结果曲线。
具体实施方式
参照图1-3。本发明基于传声器镜像的室内声源定位方法具体步骤如下:
步骤一、在室内布置多个传声器采集数据,要求传声器为无指向性型,此外,根据所选的频频范围,传声器应选择尽可能小尺寸的,这有利于提高在测点附近的所采 集声场数据的精度。对于1000Hz,对应的半波长为0.172m。在本方法中,一般选取1000~2000的频率段,对于这一频段,传声器与声源的距离大于0.2m便可满足定位要求。在采集信号时,应保证采集工作是实时同步完成的,对于数据预处理,保证16kHz以上的采样频率,录制时间不超过30s。
步骤二、根据虚声源法,对于规则的矩形房间,对于给定点声源位置可直接计算获得该点声源不同反射级数的虚声源。同样的,对于内部的每个传声器位置,也可以直接计算获得该位置对称于各壁面的镜像及其相应的反射级数。对于尺寸为lx×ly×lz的房间,任意接收点镜像的空间坐标可表示为
xi=2l·lx±xm
yi=2m·ly±ym (1)
zi=2n·lz±zm
式中,接收点位置为rm=(xm,ym,zm)),某级接收点镜像的位置为ri=(xi,yi,zi),l、m、n是界于-∞和+∞之间的整数,镜像对应的反射级数可以由下式给出
N=|2l-Δl|+|2m-Δm|+|2n-Δn| (2)
Δl=1表示求xi式中xm前取负;Δl=0表示求xi式中xm前取正;Δm=1表示求yi式中ym前取负;Δm=0表示求yi式中ym前取正;Δn=1表示求zi式中zm前取负;Δn=0表示求zi式中zm前取正;N=0时表示传声器本身,所以无反射级数。
步骤三、对于定位环境,由于声源的虚源的求解受房间形状的限制,此方法只可用于矩形或其他根据形状可直接获取虚源的房间。对于房间的吸声条件,在平均吸声系数大于0.1时,均有高效的定位效果。理论上对于二维、三维房间问题,本方法都可以实现,然而二维到三维的变化会使字典维度急剧增加,限于实际计算机的计算能力,传统方法与改进方法均是对房间内某一二维平面内的进行声源定位。
针对目标房间,将其内部空间沿各方向按照一定的距离间隔Δx、Δy和Δz划分网格,每个网格点gj的位置信息取网格点中心位置,设总网格点数为G。
步骤四、计算各传声器位置及其一定阶数镜像的位置rni与每个网格点位置gj之间的距离,将该距离|rni-gj|与对应的镜像级数、频率联合运算,求解字典矩阵Df中的各元素,Df的维度为M×G。
式中,f为频率,c为空气中的声速,N0为传声器镜像的最高阶数。
步骤五、将传声器信号的频域表示yf和各频率的字典矩阵Df作为输入,采取Lasso优化算法求解xf。
yf=Dfxf (5)
步骤六、基于yf和Df,根据公式(5),选择l1范数作为求解方法来获取各单独频率下以及联合频率下的位置向量xf和x,二者具有相同的维度。式中,y∈CM×1,x,xf∈CG×1,Df∈CM ×G。在进行室内声定位时,房间内的网格点数目N远远大于传声器个数M。实际环境中,室内一般仅存在少数几个发声声源,因此相对于整个房间,x具有稀疏性,可以给x加稀疏约束来求解方程。
步骤七、选取基于Lasso优化算法,构建如下的单频率及多频率的目标函数:
其中,λ是需要调节的参数,用来平衡估计误差和信号稀疏度。
步骤八、根据式(6)和式(7)求解结果xf和x,其最大值对应的网格点位置为该频率点下使用本方法定位得到的声源位置。
应用实施例:
房间为一封闭矩形空间。其长宽高分别为lx=3m,ly=3m,lz=3m。
步骤1、在封闭房间内部选定同一高度z=1m,布置15阵元的线性传声器阵列,阵列位置为x=1.44m,y方向间隔为0.2m,起始端传声器位置为y=0.1m。
步骤2、根据房间尺寸,选定z=1m的平面为主要参考区域,x、y方向各划分10段来进行网格划分,即房间内部的网格点为100个,网格间距为10cm。网格坐标取网格中心位置坐标。
步骤3、根据房间壁面材料,通过查表或测量室内混响时间确定房间壁面的平均声压反射系数β0。
步骤4、随机选取一个预划分的网格点作为声源位置,按照附图2连接电脑,以及相应的采集设备。
步骤5、开启所有设备和声源,持续播放预定的白噪声信号并设定所采集时域信号的采样频率并进行录制,保存数据。对各保存数据做短时傅立叶变换,获取各传声器信号的频域表示yf。
步骤6、选定坐标参考点,测量各传声器在房间内的相对位置,根据房间内部网格点划分,求解各传声器镜像的坐标信息以及各镜像的反射级数。
步骤7、根据各传声器镜像、镜像反射级数以及各壁面吸声系数,根据公式(3)求解房间内各网格点对应的字典原子,形成各频率点下的字典Df。
步骤8、基于yf和Df,根据式(5)、(6)、(7),选择l1范数作为超定方程组的求解方法来获取各单独频率下的位置向量xf以及联合各频率点数据求解的位置向量x。
步骤9、根据求解结果,xf、x中最大值对应的网格点即为该频率点下使用本方法定位得到的声源位置。
在本实施例中,分别对得到的采集信号做不同的信噪比处理并进行定位,选定频率范围为2kHz~2.3kHz。分别统计两种方法的计算时间,结果见表1。
表1两种方法计算时间对比统计
根据计算效率结果,背景技术方法的运行速度受房间混响程度的影响,在房间混响较大时,需要划分较多的网格点来进行求解,因而耗费时间较多,而本发明方法的计算时间不受房间混响程度的影响,它只于房间划分的网格点数有直接的关系,网格点数目越多,求解向量维度越高,耗费时间越多。
最后,验证所提方法在不同信噪比下的定位精度性能,参照图3,单一频率点的信号定位准确率达到85%以上需要信噪比高于15dB,而对于多频点的联合求解,在信噪比高于-5dB时,就可以达到100%的定位精度,这证明了本发明方法的有效性。
Claims (1)
1.一种基于传声器镜像的室内声源定位方法,其特征在于包括以下步骤:
步骤一、在室内布置多个传声器,所述传声器为无指向性型;对于1000~2000频率段,传声器与声源的距离大于0.2m;采集信号工作实时同步完成,数据预处理的采样频率为16kHz以上,录制时间不超过30s;
步骤二、对于尺寸为lx×ly×lz的房间,任意接收点镜像的空间坐标表示为
式中,接收点位置为rm=(xm,ym,zm)),某级接收点镜像的位置为ri=(xi,yi,zi),l、m、n是界于-∞和+∞之间的整数,镜像对应的反射级数由(2)式给出
N=|2l-Δl|+|2m-Δm|+|2n-Δn| (2)
Δl=1表示求xi式中xm前取负;Δl=0表示求xi式中xm前取正;Δm=1表示求yi式中ym前取负;Δm=0表示求yi式中ym前取正;Δn=1表示求zi式中zm前取负;Δn=0表示求zi式中zm前取正;N=0时表示传声器本身,所以无反射级数;
步骤三、对于房间内某一二维平面,平均吸声系数大于0.1,房间内部空间沿各方向按照距离间隔Δx、Δy和Δz划分网格,每个网格点位置gj的位置信息取网格点中心位置,设总网格点数为G;
步骤四、计算各传声器位置及其一定阶数镜像的位置与每个网格点位置gj之间的距离将所述距离与对应的镜像级数和频率联合运算,求解字典矩阵Df中的各元素,Df的维度为M×G;
式中,f为频率,c为空气中的声速,N0为传声器镜像的最高阶数;
步骤五、将传声器信号的频域表示yf和各频率的字典矩阵Df作为输入,采取Lasso优化算法求解xf;
yf=Dfxf (5)
步骤六、基于频域表示yf和字典矩阵Df,根据公式(5),选择l1范数作为求解方法获取各单独频率下以及联合频率下的位置向量xf和x,二者具有相同的维度;式(5)中,yf∈CM×1,x,xf∈CG×1,Df∈CM×G;在进行室内声定位时,房间内的网格点数目G远远大于传声器个数M;相对于整个房间,x具有稀疏性,给x加稀疏约束求解方程;
步骤七、选取基于Lasso优化算法,构建如下的单频率及多频率的目标函数:
其中,λ是需要调节的参数,用来平衡估计误差和信号稀疏度;
步骤八、根据式(6)和式(7)求解结果xf和x,其最大值对应的网格点位置为该频率点下定位得到的声源位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611230608.9A CN106842111B (zh) | 2016-12-28 | 2016-12-28 | 基于传声器镜像的室内声源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611230608.9A CN106842111B (zh) | 2016-12-28 | 2016-12-28 | 基于传声器镜像的室内声源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106842111A CN106842111A (zh) | 2017-06-13 |
CN106842111B true CN106842111B (zh) | 2019-03-29 |
Family
ID=59113234
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611230608.9A Active CN106842111B (zh) | 2016-12-28 | 2016-12-28 | 基于传声器镜像的室内声源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106842111B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107942314B (zh) * | 2017-11-22 | 2021-06-04 | 中南大学 | 基于lasso特征提取的多普勒穿墙雷达定位方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201399A (zh) * | 2007-12-18 | 2008-06-18 | 北京中星微电子有限公司 | 一种声源定位方法及系统 |
CN102901950A (zh) * | 2012-09-20 | 2013-01-30 | 浙江工业大学 | 平面阵列识别声源三维坐标的方法 |
JP2013545382A (ja) * | 2010-10-28 | 2013-12-19 | フラウンホーファー−ゲゼルシャフト・ツール・フェルデルング・デル・アンゲヴァンテン・フォルシュング・アインゲトラーゲネル・フェライン | 方向情報を取得する装置および方法、ならびにシステムおよびコンピュータプログラム |
CN103995252A (zh) * | 2014-05-13 | 2014-08-20 | 南京信息工程大学 | 一种三维空间声源定位方法 |
CN104407328A (zh) * | 2014-11-20 | 2015-03-11 | 西北工业大学 | 基于空间脉冲响应匹配的封闭空间声源定位方法及系统 |
CN106030331A (zh) * | 2013-10-01 | 2016-10-12 | 奥尔德巴伦机器人公司 | 声源定位方法和使用这种方法的仿人机器人 |
-
2016
- 2016-12-28 CN CN201611230608.9A patent/CN106842111B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201399A (zh) * | 2007-12-18 | 2008-06-18 | 北京中星微电子有限公司 | 一种声源定位方法及系统 |
JP2013545382A (ja) * | 2010-10-28 | 2013-12-19 | フラウンホーファー−ゲゼルシャフト・ツール・フェルデルング・デル・アンゲヴァンテン・フォルシュング・アインゲトラーゲネル・フェライン | 方向情報を取得する装置および方法、ならびにシステムおよびコンピュータプログラム |
CN102901950A (zh) * | 2012-09-20 | 2013-01-30 | 浙江工业大学 | 平面阵列识别声源三维坐标的方法 |
CN106030331A (zh) * | 2013-10-01 | 2016-10-12 | 奥尔德巴伦机器人公司 | 声源定位方法和使用这种方法的仿人机器人 |
CN103995252A (zh) * | 2014-05-13 | 2014-08-20 | 南京信息工程大学 | 一种三维空间声源定位方法 |
CN104407328A (zh) * | 2014-11-20 | 2015-03-11 | 西北工业大学 | 基于空间脉冲响应匹配的封闭空间声源定位方法及系统 |
Non-Patent Citations (1)
Title |
---|
一种用于舱室声源定位的双传声器方法及其性能研究;曾向阳等;《西北工业大学学报》;20160430;第34卷(第2期);189-193 |
Also Published As
Publication number | Publication date |
---|---|
CN106842111A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103583054B (zh) | 用于产生音频输出信号的装置和方法 | |
JP5814476B2 (ja) | 空間パワー密度に基づくマイクロフォン位置決め装置および方法 | |
Dmochowski et al. | On spatial aliasing in microphone arrays | |
CN106842112B (zh) | 强混响环境下基于参数化贝叶斯字典学习的声源定位方法 | |
JPH09512676A (ja) | 適応性ビーム形成方法及び装置 | |
Tervo et al. | Acoustic reflection localization from room impulse responses | |
RU2005105753A (ru) | Система и способ для формирования луча с использованием микрофонной решетки | |
KR20130116299A (ko) | 음향 삼각 측량에 의한 공간 선택적 사운드 취득 장치 및 방법 | |
Richter et al. | On the influence of continuous subject rotation during high-resolution head-related transfer function measurements | |
CN109188362A (zh) | 一种麦克风阵列声源定位信号处理方法 | |
CN107566969A (zh) | 一种封闭环境内部低频声场重构方法 | |
CN110221249A (zh) | 基于压缩感知的宽带声源定位方法 | |
Pezzoli et al. | Sparsity-based sound field separation in the spherical harmonics domain | |
Ajdlee et al. | Plenacoustic function on the circle with application to HRTF interpolation | |
CN106842111B (zh) | 基于传声器镜像的室内声源定位方法 | |
CN104703111B (zh) | 一种房间混响合成方法 | |
US10375501B2 (en) | Method and device for quickly determining location-dependent pulse responses in signal transmission from or into a spatial volume | |
US20160249153A1 (en) | Real-time loudspeaker distance estimation with stereo audio | |
KR20090128221A (ko) | 음원 위치 추정 방법 및 그 방법에 따른 시스템 | |
Szwajcowski | Approximating head-related transfer functions in the domain of common basis functions | |
Pertilä | Acoustic source localization in a room environment and at moderate distances | |
JP5826502B2 (ja) | 音響処理装置 | |
ÇATALBAŞ et al. | 3D moving sound source localization via conventional microphones | |
Firoozabadi et al. | Combination of nested microphone array and subband processing for multiple simultaneous speaker localization | |
Mazur et al. | Robust room equalization using sparse sound-field reconstruction |
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 |