CN112305497A - 一种近场麦克风阵列doa估计测向模糊消除方法 - Google Patents
一种近场麦克风阵列doa估计测向模糊消除方法 Download PDFInfo
- Publication number
- CN112305497A CN112305497A CN202011029596.XA CN202011029596A CN112305497A CN 112305497 A CN112305497 A CN 112305497A CN 202011029596 A CN202011029596 A CN 202011029596A CN 112305497 A CN112305497 A CN 112305497A
- Authority
- CN
- China
- Prior art keywords
- array
- doa estimation
- representing
- signal
- direction finding
- 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
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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/86—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves with means for eliminating undesired waves, e.g. disturbing noises
-
- 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
- G01S3/803—Systems for determining direction or deviation from predetermined direction using amplitude comparison of signals derived from receiving transducers or transducer systems having differently-oriented directivity characteristics
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (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
本发明公开了近场麦克风阵列DOA估计测向模糊消除方法,包括以下步骤:S1、将阵元数为M的阵列进行分布式阵列布阵,将M个阵元按照接收信号功率从大到小标记并构建集合A;S2、选取集合A中与发声源直线距离最短的麦克风阵元M1作为参考阵元;在剩下的M‑1个阵元中选取接收信号功率最大阵元M2、M3与参考阵元一起构成V1阵列;从剩下的M‑3个阵元中选择接收信号功率最大阵元M4、M5与参考阵元一起构成V2阵列;S3、将宽带信号划分为L个中心角频率分别为ω1、ω2,……,ωL的窄带信号,通过MUSIC算法分别对V1和V2阵列接收信号进行信号DOA估计,得到两个空间谱,然后对两个空间谱进行归一化处理;S4、将V1阵列和V2阵列DOA估计结果的空间谱进行谱峰相减,滤除虚假峰。
Description
技术领域
本发明涉及近场阵列解模糊技术领域,具体涉及一种近场麦克风阵列DOA估计测向模糊消除方法。
背景技术
阵列信号处理技术广泛应用于雷达、声纳、地震勘探、无线通信等领域。波达方向(DOA)估计是阵列信号处理技术的重要分支,DOA估计可以通过阵列接收数据准确地获得信源的空间位置和方位信息,与传统规则阵列相比,分布式阵列具有布阵灵活,分辨力高等优点,在舰载雷达目标探测、机场声源定位等场合具有广阔的应用前景。但是,当分布式阵列阵元间距大于半波长时,不满足空域奈奎斯特采样定律,会产生测向模糊问题。
发明内容
本发明的目的在于解决上述问题,提供一种近场麦克风阵列DOA估计测向模糊消除方法。
本发明的目的之一是这样实现的,提供一种近场麦克风阵列DOA估计测向模糊消除方法,包括以下步骤:
S1、将阵元数为M的阵列进行分布式阵列布阵,将所述M个阵元按照接收信号功率从大到小标记并构建集合A;
S2、选取集合A中与发声源直线距离最短的麦克风阵元M1作为参考阵元;在剩下的M-1个阵元中选取接收信号功率最大阵元M2、M3与参考阵元一起构成V1阵列;从剩下的M-3个阵元中选择接收信号功率最大阵元M4、M5与参考阵元一起构成V2阵列;
S3、将宽带信号划分为L个中心角频率分别为ω1、ω2,……,ωL的窄带信号,通过MUSIC算法分别对V1和V2阵列接收信号进行信号DOA估计,得到两个空间谱,然后对所述两个空间谱进行归一化处理;
S4、将V1阵列和V2阵列DOA估计结果的空间谱进行谱峰相减,滤除虚假峰。
相对现有技术,本发明的显著进步性至少体现在:以同一参考阵元为基础,构建共点双V形麦克风阵列组合进行解模糊,根据不同麦克风阵列的DOA虚假峰位置不同,而真实峰的位置保持不变的特性,将量V形阵列DOA估计的谱峰差值与设定阈值比较,真实谱峰保持不,滤除虚假峰,有效的解决了分布式麦克风阵列选阵后存在的测向模糊问题。
附图说明
图1为本发明实施例的分布式麦克风阵列示意图;
图2为本发明实施例方法步骤示意图;
图3为本发明实施例的方法流程图。
具体实施方式
下面结合附图和具体实施例,对本发明作详细说明。应该说明的是,本发明的具体实施方式并不限于所提供的具体实施例。
参阅图1-3所示,本发明提供一种近场麦克风阵列DOA估计测向模糊消除方法,包括以下步骤:
S1、将阵元数为M的阵列进行分布式阵列布阵,将所述M个阵元按照接收信号功率从大到小标记并构建集合A;
S2、选取集合A中与发声源直线距离最短的麦克风阵元M1作为参考阵元;在剩下的M-1个阵元中选取接收信号功率最大阵元M2、M3与参考阵元一起构成V1阵列;从剩下的M-3个阵元中选择接收信号功率最大阵元M4、M5与参考阵元一起构成V2阵列;
S3、将宽带信号划分为L个中心角频率分别为ω1、ω2,……,ωL的窄带信号,通过MUSIC算法分别对V1和V2阵列接收信号进行信号DOA估计,得到两个空间谱,然后对所述两个空间谱进行归一化处理;
S4、将V1阵列和V2阵列DOA估计结果的空间谱进行谱峰相减,滤除虚假峰。
V形阵列的建立较其他多元阵列而言简单,并且其构建难度接近于线阵,能够保证计算的效率,并且分布式阵列具有布阵灵活,分辨力高等优点。
作为优选的方案,以参考阵元为坐标原点构建空间坐标系,计算声源Sk到第i个阵元的距离,表征为:
如图1所述,以参考阵元为坐标原点,假设某一声源Sk入射至如图1所示的分布式阵列,则第i个阵元的位置矢量为:
通过第i个阵元的位置矢量计算得到上述声源相对于参考阵元的位置矢量,进而计算出声源到第i个阵元的距离,以此方式验证阵元接收信号功率以及参考阵元与声源距离关系的合理性,确保集合A中对于阵元与功率标记的准确性。
为保证具有较高的信噪比,选择与发射声源距离最近的麦克风阵元作为参考阵元,其接收信号的功率最大。参考阵元确定后,按照同样的选阵策略,根据功率大小选择其他阵元,与参考阵元一起构造新的阵列;
设单个入射声源S入射至分布式麦克风阵列,选用V形阵列进行DOA估计。再根据功率的大小,选取第二个阵元。假设第二个阵元坐标为(x2,y2,0),它与参考阵元距离为d。为了满足近场模型以及阵元间距小于半波长的条件,理想情况下,d应满足:
其中,r是声源到参考阵元的距离,λmin和λmax分别表示宽带信号最高频率对应的波长和最低频率对应的波长。
通常情况下,建立L形阵列即线阵进行DOA估计,实际上L阵的特殊情形,在选择第三个阵元时,可以优先选择能与前面两个阵元构成L阵的阵元。假设第三个阵元坐标为(x3,y3,0),若第三个阵元与前两个阵元能构成L阵,则其坐标应满足余弦定理:
若找不到坐标满足上式关系的第三个阵元,则寻找使得上式中cosα尽量小的阵元,使其构造的V形阵列接近L阵。
作为优选的方案,所述步骤S2中,将观测时间长度为T的宽带信号划分为J段长度为Td的时间子序列,同时将所述宽带信号划分为Q个窄带,再进行Q点离散傅里叶变换,得到宽带信号频域模型,结合所述宽带信号频域模型进行V形阵列构建,所述宽带信号频域模型表征为:
Xj(fq)=Aj(fq)Sj(fq)+Nj(fq)
j=1,2,L,J
q=1,2,L,Q (3)
其中Xj(fq),Sj(fq),Nj(fq)分别表示在频率fq下的阵列接收信号、源信号及噪声信号的傅里叶变换。
作为优选的方案,将选取的M2、M3、M4、M5放入所述空间坐标系中进行V形阵列构建,M2、M3、M4、M5是除去擦考阵元外按照集合A中接收功率从大到小选取出前4个阵元。
作为优选的方案,所述V1阵列表征为:
其中a1表示V1阵列,r2表示声源到阵元M2的距离,r3表示声源到阵元M3的距离,ω表示入射信号的角频率,c表示声音传播的速度,j表示宽带信号段数,e表示自然对数函数的底数。
麦克风阵列接收到信号的时域矢量形式为:
MUSIC算法是一种基于特征分解的子空间算法,在使用二维MUSIC算法前,先对接收信号加窗分帧,再进行短时傅里叶变换,得到接收信号的频域表达式:
其中,ω为角频率。由于信源是宽带信号,需进行窄带划分,将宽带信号划分为中心角频率分别是ω1,ω2,L,ω3的L个窄带信号,设第l个窄带信号的频率为ω1,则有:
X(ωl)的协方差矩阵为:
其中,RS(ωl)为S(ωl)的自相关矩阵,为噪声方差,I为单位矩阵。与经典MUSIC算法相同,首先对RS(ωl)进行特征值分解,按从大到小的顺序排列,其中,大于的特征值对应的特征向量与信号相关,构成信号子空间US,等于的特征值对应的特征向量与噪声相关,构成噪声子空间UN。
作为优选的方案,所述V2阵列表征为:
其中a2表示V2阵列,r4表示声源到阵元M4的距离,r5表示声源到阵元M5的距离。
作为优选的方案,所述空间图谱表征为:
信号子空间US和噪声子空间UN具有正交性,进行谱峰搜索便可以得到第l个窄带的MUSIC空间谱,其中表示第l个窄带信号的MUSIC空间谱,ωl表示第l个窄带信号的中心角频率,r,θ分别表示第l个窄带信号与参考阵元的距离、俯仰角、方位角,表示导向矢量的共轭转置,UN表示噪声子空间,表示噪声子空间的共轭转置,表示第l个窄带信号对应的导向矢量。
整个宽带信号的平均空间谱为:
根据平均空间谱,找到谱峰峰值点所对应的位置,该位置即为所要估计的目标信源的位置。
作为优选的方案,所述步骤S4中,预设真假峰判断阈值,将相减得到的差值与所述阈值比较,差值大于所述阈值对应为虚假峰,小于所述阈值对应为真实峰。
作为优选的方案,所述判断阈值为0到0.8之间的任意数值,可以理解的是,归一化处理后得到的各个峰的数值为0到1之间的任意数,由于真实峰的位置保持不变,两个空间谱相减后真实峰相减后的数值为0,而虚假峰由于位置随机,很难重叠相减,并且噪声的功率相对主要声源的功率而言较小,选取的参考阵元为距离声源最近且接收功率最大,对于阈值的选择需要参考当前声源所处的环境,容易理解的是数值越接近于0的越接近真实峰,根据多次试验的经验选取合理的阈值能够一定程度上保证模糊消除的准确性。
Claims (9)
1.一种近场麦克风阵列DOA估计测向模糊消除方法,其特征在于,包括以下步骤:
S1、将阵元数为M的阵列进行分布式阵列布阵,将所述M个阵元按照接收信号功率从大到小标记并构建集合A;
S2、选取集合A中与发声源直线距离最短的麦克风阵元M1作为参考阵元;在剩下的M-1个阵元中选取接收信号功率最大阵元M2、M3与参考阵元一起构成V1阵列;从剩下的M-3个阵元中选择接收信号功率最大阵元M4、M5与参考阵元一起构成V2阵列;
S3、将宽带信号划分为L个中心角频率分别为ω1、ω2,……,ωL的窄带信号,通过MUSIC算法分别对V1和V2阵列接收信号进行信号DOA估计,得到两个空间谱,然后对所述两个空间谱进行归一化处理;
S4、将V1阵列和V2阵列DOA估计结果的空间谱进行谱峰相减,滤除虚假峰。
3.根据权利要求1所述的近场麦克风阵列DOA估计测向模糊消除方法,其特征在于,所述步骤S2中,将观测时间长度为T的宽带信号划分为J段长度为Td的时间子序列,同时将所述宽带信号划分为Q个窄带,再进行Q点离散傅里叶变换,得到宽带信号频域模型,结合所述宽带信号频域模型进行V形阵列构建,所述宽带信号频域模型表征为:
Xj(fq)=Aj(fq)Sj(fq)+Nj(fq)
j=1,2,L,J
q=1,2,L,Q (3)
其中Xj(fq),Sj(fq),Nj(fq)分别表示在频率fq下的阵列接收信号、源信号及噪声信号的傅里叶变换。
4.根据权利要求3所述的近场麦克风阵列DOA估计测向模糊消除方法,其特征在于,将选取的M2、M3、M4、M5放入所述空间坐标系中进行V形阵列构建。
8.根据权利要求1所述的近场麦克风阵列DOA估计测向模糊消除方法,其特征在于,所述步骤S4中,预设真假峰判断阈值,将相减得到的差值与所述阈值比较,差值大于所述阈值对应为虚假峰,小于所述阈值对应为真实峰。
9.根据权利要求8所述的近场麦克风阵列DOA估计测向模糊消除方法,其特征在于,所述判断阈值为0到0.8之间的任意数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011029596.XA CN112305497B (zh) | 2020-09-27 | 2020-09-27 | 一种近场麦克风阵列doa估计测向模糊消除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011029596.XA CN112305497B (zh) | 2020-09-27 | 2020-09-27 | 一种近场麦克风阵列doa估计测向模糊消除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112305497A true CN112305497A (zh) | 2021-02-02 |
CN112305497B CN112305497B (zh) | 2023-09-01 |
Family
ID=74488731
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011029596.XA Active CN112305497B (zh) | 2020-09-27 | 2020-09-27 | 一种近场麦克风阵列doa估计测向模糊消除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112305497B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965026A (zh) * | 2021-02-08 | 2021-06-15 | 华南理工大学 | 一种具有先验定位角度范围的doa阵元间距设置方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5525996A (en) * | 1995-02-10 | 1996-06-11 | Applied Concepts, Inc. | Police traffic radar for calculating and simultaneously displaying fastest target speed |
TW200824485A (en) * | 2006-11-20 | 2008-06-01 | Matsushita Electric Ind Co Ltd | Device and method for detecting sound |
CN101193460A (zh) * | 2006-11-20 | 2008-06-04 | 松下电器产业株式会社 | 检测声音的装置及方法 |
CN101231385A (zh) * | 2008-02-26 | 2008-07-30 | 四川大学 | 用于啁啾脉冲放大的光谱整形调制方法 |
CN102830386A (zh) * | 2012-09-03 | 2012-12-19 | 西安建筑科技大学 | 一种强干扰背景下的任意阵型微弱信号源角度估计方法 |
US20140241549A1 (en) * | 2013-02-22 | 2014-08-28 | Texas Instruments Incorporated | Robust Estimation of Sound Source Localization |
CN104392114A (zh) * | 2014-11-11 | 2015-03-04 | 西北大学 | 一种基于空时数据的高分辨目标方位估计方法 |
US20200182995A1 (en) * | 2015-07-17 | 2020-06-11 | Origin Wireless, Inc. | Method, apparatus, and system for outdoor target tracking |
EP3695783A1 (en) * | 2019-02-15 | 2020-08-19 | Origin Wireless, Inc. | Method, apparatus, and system for wireless gait recognition |
-
2020
- 2020-09-27 CN CN202011029596.XA patent/CN112305497B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5525996A (en) * | 1995-02-10 | 1996-06-11 | Applied Concepts, Inc. | Police traffic radar for calculating and simultaneously displaying fastest target speed |
TW200824485A (en) * | 2006-11-20 | 2008-06-01 | Matsushita Electric Ind Co Ltd | Device and method for detecting sound |
CN101193460A (zh) * | 2006-11-20 | 2008-06-04 | 松下电器产业株式会社 | 检测声音的装置及方法 |
CN101231385A (zh) * | 2008-02-26 | 2008-07-30 | 四川大学 | 用于啁啾脉冲放大的光谱整形调制方法 |
CN102830386A (zh) * | 2012-09-03 | 2012-12-19 | 西安建筑科技大学 | 一种强干扰背景下的任意阵型微弱信号源角度估计方法 |
US20140241549A1 (en) * | 2013-02-22 | 2014-08-28 | Texas Instruments Incorporated | Robust Estimation of Sound Source Localization |
CN104392114A (zh) * | 2014-11-11 | 2015-03-04 | 西北大学 | 一种基于空时数据的高分辨目标方位估计方法 |
US20200182995A1 (en) * | 2015-07-17 | 2020-06-11 | Origin Wireless, Inc. | Method, apparatus, and system for outdoor target tracking |
EP3695783A1 (en) * | 2019-02-15 | 2020-08-19 | Origin Wireless, Inc. | Method, apparatus, and system for wireless gait recognition |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965026A (zh) * | 2021-02-08 | 2021-06-15 | 华南理工大学 | 一种具有先验定位角度范围的doa阵元间距设置方法 |
CN112965026B (zh) * | 2021-02-08 | 2023-07-18 | 华南理工大学 | 一种具有先验定位角度范围的doa阵元间距设置方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112305497B (zh) | 2023-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN108828522B (zh) | 一种利用垂直阵lcmv波束形成的水下目标辐射噪声测量方法 | |
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN105301563B (zh) | 一种基于一致聚焦变换最小二乘法的双声源定位方法 | |
CN106249244B (zh) | 一种声矢量圆阵宽带相干信号源方位估计方法 | |
CN106950529B (zh) | 声矢量近场源esprit和music参数估计方法 | |
WO2009114578A1 (en) | Autonomous sonar system and method | |
CN103353588B (zh) | 基于天线均匀平面阵的二维波达方向角估计方法 | |
CN104360310A (zh) | 一种多目标近场源定位方法和装置 | |
WO2003096587A9 (en) | System and method for localizing targets using multiple arrays | |
CN103713276A (zh) | 基于最小互熵谱分析的波达方向估计方法 | |
CN109696657A (zh) | 一种基于矢量水听器的相干声源定位方法 | |
CN111650556A (zh) | 一种宽带辐射源参数估计方法 | |
CN108089146B (zh) | 一种对预估角误差鲁棒的高分辨宽带波达方向估计方法 | |
CN115436896A (zh) | 快速的雷达单快拍music测角方法 | |
CN106526563A (zh) | 一种基于互相关虚拟阵的五元体积阵多目标方位估计方法 | |
CN109491009B (zh) | 一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法 | |
CN112305497A (zh) | 一种近场麦克风阵列doa估计测向模糊消除方法 | |
CN116559778B (zh) | 一种基于深度学习的车辆鸣笛定位方法及系统 | |
CN112098938B (zh) | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 | |
Chen et al. | Insight into split beam cross-correlator detector with the prewhitening technique | |
Rao et al. | DOA estimation using sparse vector sensor arrays | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
Jiang et al. | Beam-time delay domain deconvolved scheme for high-resolution active localization of underwater targets | |
CN114114163A (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 |