CN106997037A - 声矢量传感器阵列空间旋转解相干到达角估计方法 - Google Patents

声矢量传感器阵列空间旋转解相干到达角估计方法 Download PDF

Info

Publication number
CN106997037A
CN106997037A CN201710098048.4A CN201710098048A CN106997037A CN 106997037 A CN106997037 A CN 106997037A CN 201710098048 A CN201710098048 A CN 201710098048A CN 106997037 A CN106997037 A CN 106997037A
Authority
CN
China
Prior art keywords
array
axis
data
acoustic vector
covariance matrix
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.)
Pending
Application number
CN201710098048.4A
Other languages
English (en)
Inventor
王兰美
廖桂生
魏兵
孙长征
张艳艳
代少玉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201710098048.4A priority Critical patent/CN106997037A/zh
Publication of CN106997037A publication Critical patent/CN106997037A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Direction-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/80Direction-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/8003Diversity systems specially adapted for direction finding
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S11/00Systems for determining distance or velocity not using reflection or reradiation
    • G01S11/14Systems for determining distance or velocity not using reflection or reradiation using ultrasonic, sonic, or infrasonic waves

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)

Abstract

声矢量传感器阵列空间旋转解相干到达角估计方法,阵列接收K个相干窄带、平稳远场信号,利用接收阵列获取所有阵元的一次快拍输出数据,测量该阵列的N次快拍数据,充分利用声矢量传感器自身具有的正交矢量特性,将声矢量传感器阵列分成声压子阵、x轴的振速子阵、y轴的振速子阵和z轴的振速子阵四个子阵,通过子阵级处理恢复信号协方差矩阵的秩,得到解相干后的数据协方差矩阵;由解相干后的数据协方差矩阵获取噪声空间,利用噪声子空间构造MUSIC空间谱,通过角度域的一维MUSIC谱峰搜索得到到达角的估计值;相比于空间平滑解相干方法,本发明方法不仅适用于均匀线阵也适用于非均匀线阵,不存在孔径损失,具有更高的空间分辨力和测向精度,更广的适用范围。

Description

声矢量传感器阵列空间旋转解相干到达角估计方法
技术领域
本发明属于信号处理技术领域,尤其涉及一种声矢量传感器阵列的相干源到达角估计方法。
背景技术
在实际中由于信号传输过程的多径现象和回波信号以及人为干扰,使得相关信号源非常普遍,所以基于相关源的检测和估计是信号处理的重要研究课题。当入射信号相互独立时子空间类方法有很好的参数估计性能,但是如果入射的信号是相干信号以MUSIC为代表的子空间类算法将失效。当入射信号是相干信号时,信号子空间的维数将小于信号源个数,信号子空间扩散到了噪声子空间,从而导致信号子空间和噪声子空间不完全垂直,从而导致子空间类算法失效,无法进行方向估计。
为使相干信源DOA估计中的上述难题得到有效的解决,国内外学者提出了一系列的解相干方法,这些方法的关键是如何恢复信号协方差矩阵的秩。目前解相干的方法大致可以分为两类,一类是以空间平滑和矩阵重构法为代表的降维方法,一类是以聚焦和虚拟变换法为代表的非降维方法,这些方法解决了信号相关问题但仍有一定的局限性和不足,以聚焦处理为代表的非降维算法的计算复杂度高;以空间平滑为代表的降维方法减小了阵列孔径,增大了阵列的波束宽度,降低了参数估计精度。且空间平滑一般只适用于均匀阵,严重限制了方法的应用范围。声矢量传感器是一种新型的声源信号测向设备,它是由三个相互正交的质点振速传感器和一个声压传感器构成,因而能够空间共点同步测量声场中某处的声压强度和质点振速。声矢量传感器阵列与标量传感器阵列相比较,声矢量传感器阵列不仅能够获取阵列孔径信息,而且蕴含矢量传感器各分量之间的正交信息,因而具有更高的空间分辨力和测向精度,近年来已成为国内外学者研究的热点问题。本发明针对现有方法的不足利用声矢量传感器阵列自身的矢量结构特性,将声矢量传感器阵列分成声压子阵、x轴振速子阵、y轴振速子阵和z轴振速子阵四个子阵,根据子阵的旋转不变特性解相干,该方法适用于均匀和非均匀声矢量传感器线阵。
发明内容
本发明的目的是提供一种远场窄带相干源声矢量传感器阵列解相干到达角估计方法。
为了实现上述目的,本发明采取如下的技术解决方案:
声矢量传感器阵列空间旋转解相干到达角估计方法,K个相干窄带、平稳远场声源信号从不同的方向θk入射到该接收阵列上,θk∈[0,π/2]是第k个信号的俯仰角,k=1,…,K,所述阵列由M个在x轴上任意分布的声矢量传感器阵元构成,所述阵元是具有空间共点同步测量声压以及x轴、y轴和z轴方向振速分量的声矢量传感器,所有传感器的对应通道相互平行:所有的声压传感器相互平行,所有的x轴方向振速传感器相互平行,所有的y轴方向振速传感器相互平行,以及所有的z轴方向振速传感器相互平行;任意两个相邻阵元间距小于等于λmin/2,λmin为入射声波信号的最小波长;
声矢量传感器阵列空间旋转解相干到达角估计方法步骤如下:
步骤一、M个任意分布在x轴的声矢量传感器阵元构成非均匀线阵接收阵列,阵列输出N次同步采样数据Z;
步骤二、通过子阵级处理恢复信号协方差矩阵的秩,得到解相干后的协方差数据矩阵RZ
根据阵列数据Z的排布规律将数据分成声压以及x轴、y轴和z轴方向振速子阵数据,计算四个子阵数据的协方差矩阵其中,分别是声压以及x轴、y轴和z轴方向振速子阵数据协方差矩阵,通过四个子阵数据协方差矩阵的算术平均得到解相干后的满秩数据协方差矩阵RZ
步骤三、由解相干后的数据协方差矩阵RZ获取噪声空间;
通过对数据协方差矩阵RZ进行特征分解,K个大特征对应的特征矢量构成信号子空间Us,M-K个小特征值对应的特征矢量构成噪声子空间Un
步骤四、利用噪声子空间Un构造MUSIC空间谱通过角度域θ∈[0,π/2]的一维MUSIC谱峰搜索得到到达角的估计值
本发明采用的x轴上任意分布的线阵,阵列的阵元为由声压传感器和x轴、y轴及z轴方向的振速传感器构成的声矢量传感器,并且所有的声压传感器相互平行,所有的x轴方向振速传感器相互平行,所有的y轴方向振速传感器相互平行,所有的z轴方向振速传感器相互平行。本发明利用声矢量传感器自身具有的正交矢量特性,给出了声矢量传感器空间旋转解相干方法,本发明方法没有减小阵列孔径信息,具有更高的空间分辨力和测向精度,打破了原有空间平滑方法仅适用于均匀阵列的限制,扩大了方法的适用范围。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中需要使用的附图做简单介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例声矢量传感器阵列的示意图;
图2为本发明方法的流程图;
图3为信噪比是0dB时三种方法在均匀线阵条件下到达角估计性能比较;
图4为信噪比是2dB时三种方法在均匀线阵条件下到达角估计性能比较;
图5为信噪比是4dB时三种方法在均匀线阵条件下到达角估计性能比较;
图6为信噪比是0dB时三种方法在非均匀线阵条件下到达角估计性能比较;
图7为信噪比是2dB时三种方法在非均匀线阵条件下到达角估计性能比较;
图8为信噪比是4dB时三种方法在非均匀线阵条件下到达角估计性能比较;
具体实施方式
为了让本发明的上述和其它目的、特征及优点能更明显,下文特举本发明实施例,并配合所附图示,做详细说明如下。
图1所示为本发明实施例的声矢量传感器阵列的示意图。本发明的声矢量传感器阵列由M个x轴上任意分布的阵元构成,所述阵元是具有空间共点同步测量声压以及x轴、y轴和z轴方向振速分量的声矢量传感器;任意两个相邻阵元间距小于等于λmin/2,λmin为入射声波信号的最小波长;
参照图2,本发明的解相干到达角估计方法的步骤如下:x轴上任意分布声矢量传感器线阵接收K个相干窄带、平稳远场声源信号,K为入射声源信号的数量,
步骤一、M个任意分布在x轴的声矢量传感器阵元构成非均匀线阵接收阵列,阵列输出N次同步采样数据Z;
步骤二、通过子阵级处理恢复信号协方差矩阵的秩,得到解相干后的协方差数据矩阵RZ
根据阵列数据Z的排布规律将数据分成声压以及x轴、y轴和z轴方向振速子阵数据,计算四个子阵数据的协方差矩阵其中, 分别是声压以及x轴、y轴和z轴方向振速子阵数据协方差矩阵,通过四个子阵数据协方差矩阵的算术平均得到解相干后的满秩数据协方差矩阵RZ
步骤三、由解相干后的数据协方差矩阵RZ获取噪声空间;
通过对数据协方差矩阵RZ进行特征分解,K个大特征对应的特征矢量构成信号子空间Us,M-K个小特征值对应的特征矢量构成噪声子空间Un
步骤四、利用噪声子空间Un构造MUSIC空间谱通过角度域θ∈[0,π/2]的一维MUSIC谱峰搜索得到到达角的估计值
本发明给出了x轴上任意分布声矢量传感器线阵解相干到达角估计方法,充分利用声矢量传感器阵列自身的矢量结构特性进行解相干处理,将声矢量传感器阵列分成声压子阵、x轴振速子阵、y轴振速子阵和z轴振速子阵四个子阵,然后通过子阵数据协方差矩阵算术平均解相干,然后利用子空间方法得到噪声子空间,构造MUSIC谱,通过一维角度搜索完成到达角的估计。相比于解相干前的传统MUSIC方法和空间平滑解相干方法,本发明方法突破了现有空间平滑解相干只适用于均匀阵局限,对于均匀线阵和非均匀线阵均成立,且不损失阵列孔径,相同条件下参数估计精度高于空间平滑方法。
本发明的效果可以通过以下的仿真结果进一步说明:
仿真实验1:
仿真实验条件如下:两个相干窄带、平稳远场声源信号入射到由8个在x轴上均匀分布的声矢量传感器阵元构成的线性阵列,如图1所示,阵元间隔等于0.5λmin,入射信号的参数为:(θ1,θ2)=(15°,50°),快拍数为512次;
仿真实验结果如图3至图5所示,图3-5为信噪比分别是0dB、2dB和4dB、时,x轴上均匀线阵下的解相干前传统MUSIC方法、空间平滑解相干方法和本发明方法的到达角估计结果,从图3-5可以看出在同样信噪比条件下,空间平滑解相干方法和本发明方法都能解相干,但相比于空间平滑解相干方法,本发明方法的空间谱更尖锐,有更高的到达角参数估计精度。
仿真实验2:
仿真实验条件如下:该仿真实验的两个相干窄带、平稳远场声源信号入射到由8个在x轴上任意分布的声矢量传感器阵元构成的非均匀线阵,阵元间隔小于等于0.5λmin且随机分布,入射信号的参数为:(θ1,θ2)=(15°,50°),快拍数为512次。
仿真实验结果如图6至图8所示,图6-8为信噪比分别是0、2dB和4dB时,x轴上非均匀线阵下的解相干前传统MUSIC方法、空间平滑解相干方法和本发明方法到达角空间谱估计结果,从图6-8可以看出在同样信噪比条件下,在非均匀线阵下空间平滑解相干方法没有形成谱峰,空间平滑解相干方法失效,而本发明方法能够形成尖锐谱峰,从而说明本发明方法对非均匀线阵也适用,从而验证了本发明方法的正确性。
以上所述,仅是本发明的较佳实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。

Claims (1)

1.声矢量传感器阵列空间旋转解相干到达角估计方法,其特征在于:
所述声矢量传感器线阵由x轴上任意分布的阵元构成,所述阵元是具有空间共点同步测量声压以及x轴、y轴和z轴方向振速分量的声矢量传感器;任意两个相邻阵元间距小于等于λmin/2,λmin为入射声波信号的最小波长;
空间旋转解相干到达角估计方法的步骤如下:阵列接收K个相关窄带、平稳远场声源信号,
步骤一、M个任意分布在x轴的声矢量传感器阵元构成非均匀线阵接收阵列,阵列输出N次同步采样数据Z;
步骤二、通过子阵级处理恢复信号协方差矩阵的秩,得到解相干后的协方差数据矩阵RZ
根据阵列数据Z的排布规律将数据分成声压以及x轴、y轴和z轴方向振速子阵数据,计算四个子阵数据的协方差矩阵其中,分别是声压以及x轴、y轴和z轴方向振速子阵数据协方差矩阵,通过四个子阵数据协方差矩阵的算术平均得到解相干后的满秩数据协方差矩阵RZ
步骤三、由解相干后的数据协方差矩阵RZ获取噪声空间;
通过对数据协方差矩阵RZ进行特征分解,K个大特征对应的特征矢量构成信号子空间Us,M-K个小特征值对应的特征矢量构成噪声子空间Un
步骤四、利用噪声子空间Un构造MUSIC空间谱通过角度域θ∈[0,π/2]的一维MUSIC谱峰搜索得到到达角的估计值
CN201710098048.4A 2017-02-22 2017-02-22 声矢量传感器阵列空间旋转解相干到达角估计方法 Pending CN106997037A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710098048.4A CN106997037A (zh) 2017-02-22 2017-02-22 声矢量传感器阵列空间旋转解相干到达角估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710098048.4A CN106997037A (zh) 2017-02-22 2017-02-22 声矢量传感器阵列空间旋转解相干到达角估计方法

Publications (1)

Publication Number Publication Date
CN106997037A true CN106997037A (zh) 2017-08-01

Family

ID=59431318

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710098048.4A Pending CN106997037A (zh) 2017-02-22 2017-02-22 声矢量传感器阵列空间旋转解相干到达角估计方法

Country Status (1)

Country Link
CN (1) CN106997037A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108469599A (zh) * 2018-02-28 2018-08-31 哈尔滨工程大学 一种声矢量传感器幅度加权music测向方法
CN108919176A (zh) * 2018-02-28 2018-11-30 哈尔滨工程大学 一种单矢量传感器增秩music测向技术
CN109116296A (zh) * 2018-08-27 2019-01-01 陕西理工大学 存在阵列位置误差的多输出支持向量回归机参数估计方法
CN109696657A (zh) * 2018-06-06 2019-04-30 南京信息工程大学 一种基于矢量水听器的相干声源定位方法
CN109870671A (zh) * 2017-12-05 2019-06-11 常熟海量声学设备科技有限公司 一种稳健性高分辨的高效doa估计算法
WO2021064468A1 (en) * 2019-10-04 2021-04-08 Soundskrit Inc. Sound source localization with co-located sensor elements
CN113901657A (zh) * 2021-10-08 2022-01-07 金陵科技学院 一种复合材料分布式损伤源定位补偿方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103941221A (zh) * 2014-03-12 2014-07-23 西安电子科技大学 空间拉伸电磁矢量传感器阵列参数估计方法
CN105005038A (zh) * 2015-06-18 2015-10-28 许昌学院 一种改进的声矢量阵相干源doa估计算法
CN105510874A (zh) * 2015-12-15 2016-04-20 吉林大学 一种近场信源多参数联合估计的降维music方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103941221A (zh) * 2014-03-12 2014-07-23 西安电子科技大学 空间拉伸电磁矢量传感器阵列参数估计方法
CN105005038A (zh) * 2015-06-18 2015-10-28 许昌学院 一种改进的声矢量阵相干源doa估计算法
CN105510874A (zh) * 2015-12-15 2016-04-20 吉林大学 一种近场信源多参数联合估计的降维music方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
XIN YUAN: "Coherent Source Direction-Finding using a Sparsely-Distributed Acoustic Vector-Sensor Array", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 *
马伯乐等: "改进声矢量阵相干信号源方位估计算法", 《系统工程与电子技术》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109870671A (zh) * 2017-12-05 2019-06-11 常熟海量声学设备科技有限公司 一种稳健性高分辨的高效doa估计算法
CN108469599A (zh) * 2018-02-28 2018-08-31 哈尔滨工程大学 一种声矢量传感器幅度加权music测向方法
CN108919176A (zh) * 2018-02-28 2018-11-30 哈尔滨工程大学 一种单矢量传感器增秩music测向技术
CN108469599B (zh) * 2018-02-28 2021-11-23 哈尔滨工程大学 一种声矢量传感器幅度加权music测向方法
CN108919176B (zh) * 2018-02-28 2022-04-05 哈尔滨工程大学 一种单矢量传感器增秩music测向技术
CN109696657A (zh) * 2018-06-06 2019-04-30 南京信息工程大学 一种基于矢量水听器的相干声源定位方法
CN109116296A (zh) * 2018-08-27 2019-01-01 陕西理工大学 存在阵列位置误差的多输出支持向量回归机参数估计方法
CN109116296B (zh) * 2018-08-27 2023-07-18 陕西理工大学 存在阵列位置误差的多输出支持向量回归机参数估计方法
WO2021064468A1 (en) * 2019-10-04 2021-04-08 Soundskrit Inc. Sound source localization with co-located sensor elements
US12075210B2 (en) * 2019-10-04 2024-08-27 Soundskrit Inc. Sound source localization with co-located sensor elements
CN113901657A (zh) * 2021-10-08 2022-01-07 金陵科技学院 一种复合材料分布式损伤源定位补偿方法及系统

Similar Documents

Publication Publication Date Title
CN106997037A (zh) 声矢量传感器阵列空间旋转解相干到达角估计方法
CN104730491B (zh) 一种基于l型阵的虚拟阵列doa估计方法
CN106802403B (zh) 声矢量传感器二维阵列music解相干参数估计方法
CN111123192B (zh) 一种基于圆形阵列和虚拟扩展的二维doa定位方法
CN108375751B (zh) 多信源波达方向估计方法
CN106950529B (zh) 声矢量近场源esprit和music参数估计方法
CN102411136A (zh) 一种扩展基线解模糊的相位干涉仪测向方法
CN106932087B (zh) 圆形声矢量传感器阵列近场源多参数估计方法
CN103353588B (zh) 基于天线均匀平面阵的二维波达方向角估计方法
CN108872926A (zh) 一种基于凸优化的幅相误差校正及doa估计方法
CN111352063B (zh) 一种均匀面阵中基于多项式求根的二维测向估计方法
CN106970348B (zh) 电磁矢量传感器阵列解相干二维music参数估计方法
CN106872934B (zh) L型电磁矢量传感器阵列解相干esprit参数估计方法
CN110837076A (zh) 一种基于张量分解的矢量水听器阵列方位估计方法
CN106249196A (zh) 三分量声矢量传感器稀疏阵列四元数解模糊方法
CN106980105B (zh) 电磁矢量传感器阵列空间旋转解相干测向方法
CN111812581B (zh) 基于原子范数的球面阵列声源波达方向估计方法
CN106908754B (zh) L型声矢量传感器阵列esprit解相干参数估计方法
CN104793177A (zh) 基于最小二乘法的麦克风阵列测向方法
CN108089146A (zh) 一种对预估角误差鲁棒的高分辨宽带波达方向估计方法
CN106154220B (zh) L型简化声矢量传感器阵列多参数联合估计四元数方法
CN114460531A (zh) 一种均匀线阵music空间谱估计方法
CN112763972A (zh) 基于稀疏表示的双平行线阵二维doa估计方法及计算设备
CN107064862B (zh) 电磁矢量传感器线阵解相干music参数估计方法
CN113238184B (zh) 一种基于非圆信号的二维doa估计方法

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170801

WD01 Invention patent application deemed withdrawn after publication