CN108919176B - 一种单矢量传感器增秩music测向技术 - Google Patents
一种单矢量传感器增秩music测向技术 Download PDFInfo
- Publication number
- CN108919176B CN108919176B CN201810165467.XA CN201810165467A CN108919176B CN 108919176 B CN108919176 B CN 108919176B CN 201810165467 A CN201810165467 A CN 201810165467A CN 108919176 B CN108919176 B CN 108919176B
- Authority
- CN
- China
- Prior art keywords
- target
- vector
- sensor
- noise
- cos
- 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
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/802—Systems for determining direction or deviation from predetermined direction
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
本发明公开了一种单矢量传感器增秩MUSIC测向技术,属于传感器信号处理技术领域。MUSIC测向技术估计误差小,精度高,它要求声矢量传感器接收的海洋环境噪声的协方差矩阵为单位阵。在海洋环境噪声场中,声矢量传感器的声压通道和振速通道接收到的海洋环境噪声功率并不相等,导致MUSIC测向技术在水下目标探测中无法得到应有的高精度估计。本发明发现声矢量传感器的声压通道和振速通道接收到的环境噪声功率不一致性引起一个虚源。为保证目标导向矢量和噪声子空间的正交性,本发明把此虚源的导向矢量归于信号子空间而非噪声子空间。在较低信噪比下,本发明仍旧具有尖锐的空间谱峰和较小的估计误差。本发明可用于解决海洋环境噪声中对弱目标的被动测向问题。
Description
技术领域
本发明属于传感器信号处理技术领域,具体涉及一种单矢量传感器增秩MUSIC测向技术。
背景技术
声矢量传感器的声压通道和振速通道可共点同时获得声场的声压和振速信息,为水下目标测向提供了更有利的工具和更多信息。MUSIC测向技术是一种超分辨测向技术,MUSIC技术最初是由R.O.Schmidt(R.O.Schmidt.Multiple emitter location and signalparameter estimation[J].IEEE Transactions on Antennas and Propagation.1986,34(3),pp:276-280)在无线电测向应用中提出。近几年,曾雄风等(曾雄飞,孙贵青,李宇,黄海宁.单矢量水听器的几种DOA估计方法[J].仪器仪表学报.2012,33(3):499-507)把MUSIC测向技术拓展到单个矢量传感器水下目标的测向应用中。当声矢量传感器在水下应用时,海洋环境噪声是主要噪声源,海洋环境噪声在声压通道和振速通道的功率并不相等,孙贵青等对此做了细致的理论分析和实验验证(孙贵青,杨德森,时胜国.基于矢量水听器的声压和质点振速的空间相关系数[J].声学学报.2003,28(6):509-513)。已有的MUSIC测向技术在声矢量传感器测向中应用时,没有考虑声矢量传感器的声压通道和振速通道所接收的海洋环境噪声功率不一致性,导致MUSIC测向技术在水下目标探测中无法得到应有的超分辨能力。本发明发现声压通道和振速通道所接收的海洋环境噪声功率不一致性会引起虚源,破坏了接收数据协方差矩阵的噪声子空间和目标导向矢量之间的正交性,从而使得MUSIC测向技术在低信噪比时无法获得应有的性能。本发明考虑了此虚源的影响,提出一种增秩MUSIC测向技术,增秩MUSIC测向技术在低信噪比时仍旧能获得尖锐的谱峰,角度估计误差较小。
发明内容
本发明的目的在于提供消除虚源影响,角度估计误差较小,解决了旁瓣高、主瓣胖难题的一种单矢量传感器增秩MUSIC测向技术。
本发明的目的通过如下技术方案来实现:
声矢量传感器由声压传感器和三个空间轴向垂直的振速传感器空间共点组合而成。声矢量传感器空间共点同时测量声压以及x、y、z方向的三个振速分量。
一种单矢量传感器增秩MUSIC测向技术,包括以下步骤:
步骤一:声矢量传感器在海洋环境噪声场中接收一个远场信号,校准后的声矢量传感器的声压通道和振速通道之间不存在通道幅相误差。此时声矢量传感器输出为N个快拍数据r(n),r(n)是一个M×1的向量,M=4,n=1,…,N。
步骤四:在海洋环境噪声场中,声压通道和振速通道接收的海洋环境噪声功率不一致性导致虚源,虚拟源个数为1。此时,v2应归入信号子空间,而噪声子空间由向量组{v3,…,vM}扩展而成。
步骤六:通过搜索下面空间谱的谱峰位置,目标的角度估计如下:
其中,和分别是目标的方位角估计值和俯仰角估计值,S(θ,φ)=uH(θ,φ)Pnu(θ,φ)-1是空间谱,u(θ,φ)=[1,cos(θ)cos(φ),sin(θ)cos(φ),sin(φ)]T是声矢量传感器在角度(θ,φ)处的导向矢量,θ∈[-π,π]是搜索方位角,φ∈[-π/2,π/2]是搜索俯仰角。
本发明还可以包括:
步骤四还包括:
在各向同性噪声场条件下,声矢量传感器接收数据的协方差矩阵期望值为:
其中,u(θ1,φ1)=[1,cos(θ1)cos(φ1),sin(θ1)cos(φ1),sin(φ1)]T是目标的导向矢量,是目标功率,是声压传感器的海洋环境噪声功率。 为振速传感器接收到的海洋环境噪声功率,IM是M×M的单位阵,z1=[1,0,0,0]T是虚源的导向矢量。在非各向同性噪声场条件下,虚拟源个数仍旧可以取为1。
当声矢量传感器由声压传感器和在水平面内相互正交的二维振速传感器空间共点组合而成时,M=3,u(θ,φ)=[1,cos(θ)cos(φ),sin(θ)cos(φ)]T。
当目标的俯仰角已知时,通过搜索下面空间谱的谱峰位置,估计目标的方位角度如下:
此时,目标方位角的粗估计值如下:
当目标的方位角已知时,通过搜索下面空间谱的谱峰位置,估计目标的俯仰角度如下:
此时,目标俯仰角的粗估计值如下:
本发明的有益效果在于:
本发明发现声压通道和振速通道接收到的海洋环境噪声功率导致一个虚源,为消除此虚源的影响,提出了一种增秩MUSIC测向技术。增秩MUSIC测向技术把信号子空间的秩增加1,从而把虚源归入信号子空间,保证了噪声子空间和目标导向矢量的正交性。增秩MUSIC测向技术在低信噪比条件下,仍具有尖锐的谱峰和较高的估计精度。解决了已有MUSIC测向技术在海洋环境噪声场中应用时,旁瓣高、主瓣胖的难题。本发明主要应用于水下弱目标的被动探测。
附图说明
图1为增秩MUSIC测向技术流程图;
图2为增秩MUSIC测向技术的二维空间谱;
图3为MUSIC测向技术的二维空间谱;
图4为消声水池实验结果。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明:
实施例一:
结合图1,本发明的处理流程包含如下步骤:
声矢量传感器由声压传感器和三个空间轴向垂直的振速传感器空间共点组合而成。声矢量传感器空间共点同时测量声压以及x、y、z方向的三个振速分量。
步骤一:声矢量传感器在海洋环境噪声场中接收一个远场信号,校准后的声矢量传感器的声压通道和振速通道之间不存在通道幅相误差。声矢量传感器输出为N个快拍数据r(n),r(n)是M×1的向量,M=4,n=1,…,N。
步骤四:在海洋环境噪声场中,声压通道和振速通道接收的海洋环境噪声功率不一致性导致虚源,虚拟源个数为1。此时,v2应归入信号子空间,而噪声子空间由向量组{v3,…,vM}扩展而成。
步骤六:通过搜索下面空间谱的谱峰位置,估计目标的角度如下:
其中,和分别是射源的方位角估计值和俯仰角估计值,S(θ,φ)=uH(θ,φ)Pnu(θ,φ)-1是空间谱,u(θ,φ)=[1,cos(θ)cos(φ),sin(θ)cos(φ),sin(φ)]T是声矢量传感器在(θ,φ)角度的导向矢量,θ∈[-π,π]是搜索方位角,φ∈[-π/2,π/2]搜索俯仰角。
声矢量传感器增秩MUSIC测向技术,还包括:
当声矢量传感器由声压传感器和在水平面内相互正交的二维振速传感器空间共点组合而成时,M=3,u(θ,φ)=[1,cos(θ)cos(φ),sin(θ)cos(φ)]T。
当目标的俯仰角已知时,通过搜索下面空间谱的谱峰位置,估计目标的方位角度如下:
此时,目标方位角的粗估计值如下:
当目标的方位角已知时,通过搜索下面空间谱的谱峰位置,估计目标的俯仰角度如下:
此时,目标俯仰角的粗估计值如下:
下面通过仿真实例和试验实例对本发明具体实施做进一步描述。
仿真实例:
声矢量传感器由声压传感器和三个空间轴向垂直的振速传感器空间共点组合而成。声矢量传感器空间共点同时测量声压以及x、y、z方向的三个振速分量。海洋环境噪声为各向同性噪声,此时ηx=ηy=ηz=1/3,一个目标从(θ1,φ1)=(40°,10°)入射到声矢量传感器,信噪比为-3dB。图2为增秩MUSIC测向技术的空间谱,图3为传统MUSIC测向技术的空间谱。从图2可见,增秩MUSIC测向技术在(-180°+θ1,-φ1)附近存在角度模糊。根据图3所示,MUSIC测向技术所给出的峰值在(42.8°,10.6°),根据此角度,我们可判断在图2中,(40.1°,9.9°)是目标的角度估计值,而(-140.8°,-11.6°)是角度模糊。比较图2和图3,可知MUSIC测向技术的旁瓣较高,主瓣很宽,而本发明增秩MUSIC测向技术的旁瓣很低,主瓣很窄。另外,从图2和图3可知,本发明增秩MUSIC测向技术有较低的估计误差。
消声水池试验实例:
在哈尔滨工程大学消声水池开展声矢量目标方位角度估计实验。声矢量传感器由一个声压传感器和在水平面内相互正交的两个振速传感器组成。实验中目标声源发射单频信号,信号频率为2kHz,声源距声矢量传感器15m,声源和声矢量传感器位于同一深度,此时俯仰角φ1=0°。由于此声矢量传感器不含垂直于水平面的振速传感器,我们只给出方位角的空间谱,如图4所示。比较MUSIC测向技术和增秩MUSIC测向技术可知,增秩MUSIC测向技术在135.8度处的旁瓣对应真实的目标方位,而在-66.2度处的峰值是角度估计模糊。另外,MUSIC测向技术的旁瓣较高而且主瓣较宽,这是由于MUSIC测向技术没有考虑声矢量传感器声压通道和振速通道接收的海洋环境噪声功率不一致性造成的虚源。增秩MUSIC测向技术把此虚源归入信号子空间,从而保证了目标导向矢量和噪声子空间的正交性,得到较低旁瓣和较窄的主瓣。
实施例二:
一种单矢量传感器增秩MUSIC测向技术,包括以下步骤:
(1)声矢量传感器在海洋环境噪声场中接收一个远场信号,校准声矢量传感器使得声压通道和振速通道之间不存在通道幅相误差,声矢量传感器输出为N个快拍数据r(n),r(n)是一个M×1的向量,M=4,n=1,…,N,N为正整数;
(4)在海洋环境噪声场中,声压通道和振速通道接收的海洋环境噪声功率不一致性导致虚源,虚拟源个数为1,将v2归入信号子空间,噪声子空间由向量组{v3,…,vM}扩展而成;
(6)通过搜索空间谱的谱峰位置,目标的角度估计如下:
其中,和分别是目标的方位角估计值和俯仰角估计值,S(θ,φ)=uH(θ,φ)Pnu(θ,φ)-1是空间谱,u(θ,φ)=[1,cos(θ)cos(φ),sin(θ)cos(φ),sin(φ)]T是声矢量传感器在角度(θ,φ)处的导向矢量,θ∈[-π,π]是搜索方位角,φ∈[-π/2,π/2]是搜索俯仰角;
所述的步骤(4)具体包括:
(4.1)在各向同性噪声场条件下,声矢量传感器接收数据的协方差矩阵期望值为:
其中,u(θ1,φ1)=[1,cos(θ1)cos(φ1),sin(θ1)cos(φ1),sin(φ1)]T是目标的导向矢量,是目标功率,是声压传感器的海洋环境噪声功率, 为振速传感器接收到的海洋环境噪声功率,IM是M×M的单位阵,z1=[1,0,0,0]T是虚源的导向矢量,在非各向同性噪声场条件下,虚拟源个数取为1;
(4.2)当声矢量传感器由声压传感器和在水平面内相互正交的二维振速传感器空间共点组合而成时,M=3,u(θ,φ)=[1,cos(θ)cos(φ),sin(θ)cos(φ)]T,
(4.3)当目标的俯仰角已知时,设为φ1,通过搜索下面空间谱的谱峰位置,估计目标的方位角度如下:
此时,目标方位角的粗估计值如下:
(4.4)当目标的方位角已知时,设为θ1,通过搜索下面空间谱的谱峰位置,估计目标的俯仰角度如下:
此时,目标俯仰角的粗估计值如下:
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种单矢量传感器增秩MUSIC测向技术,其特征在于,包括以下步骤:
(1)声矢量传感器在海洋环境噪声场中接收一个远场信号,校准声矢量传感器使得声压通道和振速通道之间不存在通道幅相误差,声矢量传感器输出为N个快拍数据r(n),r(n)是一个M×1的向量,M=4,n=1,…,N,N为正整数;
(4)在海洋环境噪声场中,声压通道和振速通道接收的海洋环境噪声功率不一致性导致虚源,虚源个数为1,将v2归入信号子空间,噪声子空间由向量组{v3,…,vM}扩展而成;
(6)通过搜索空间谱的谱峰位置,目标的角度估计如下:
其中,和分别是目标的方位角估计值和俯仰角估计值,S(θ,φ)=uH(θ,φ)Pnu(θ,φ)-1是空间谱,u(θ,φ)=[1,cos(θ)cos(φ),sin(θ)cos(φ),sin(φ)]T是声矢量传感器在角度(θ,φ)处的导向矢量,θ∈[-π,π]是搜索方位角,φ∈[-π/2,π/2]是搜索俯仰角;
2.根据权利要求1所述的一种单矢量传感器增秩MUSIC测向技术,其特征在于,所述的步骤(4)具体包括:
(4.1)在各向同性噪声场条件下,声矢量传感器接收数据的协方差矩阵期望值为:
其中,u(θ1,φ1)=[1,cos(θ1)cos(φ1),sin(θ1)cos(φ1),sin(φ1)]T是目标的导向矢量,是目标功率,是声压传感器的海洋环境噪声功率, 为振速传感器接收到的海洋环境噪声功率,IM是M×M的单位阵,z1=[1,0,0,0]T是虚源的导向矢量,在非各向同性噪声场条件下,虚源个数取为1;
(4.2)当声矢量传感器由声压传感器和在水平面内相互正交的二维振速传感器空间共点组合而成时,M=3,u(θ,φ)=[1,cos(θ)cos(φ),sin(θ)cos(φ)]T,
(4.3)当目标的俯仰角已知时,设为φ1,通过搜索下面空间谱的谱峰位置,估计目标的方位角度如下:
此时,目标方位角的粗估计值如下:
(4.4)当目标的方位角已知时,设为θ1,通过搜索下面空间谱的谱峰位置,估计目标的俯仰角度如下:
此时,目标俯仰角的粗估计值如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810165467.XA CN108919176B (zh) | 2018-02-28 | 2018-02-28 | 一种单矢量传感器增秩music测向技术 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810165467.XA CN108919176B (zh) | 2018-02-28 | 2018-02-28 | 一种单矢量传感器增秩music测向技术 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108919176A CN108919176A (zh) | 2018-11-30 |
CN108919176B true CN108919176B (zh) | 2022-04-05 |
Family
ID=64403052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810165467.XA Active CN108919176B (zh) | 2018-02-28 | 2018-02-28 | 一种单矢量传感器增秩music测向技术 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108919176B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110837076A (zh) * | 2019-11-18 | 2020-02-25 | 哈尔滨工程大学 | 一种基于张量分解的矢量水听器阵列方位估计方法 |
CN113074810B (zh) * | 2021-04-01 | 2023-08-18 | 中国电子科技集团公司第三研究所 | 一种矢量传声器的校准***及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510874A (zh) * | 2015-12-15 | 2016-04-20 | 吉林大学 | 一种近场信源多参数联合估计的降维music方法 |
CN105676168A (zh) * | 2015-12-02 | 2016-06-15 | 江苏科技大学 | 一种声矢量阵方位估计方法 |
CN106249244A (zh) * | 2016-07-07 | 2016-12-21 | 哈尔滨工程大学 | 一种声矢量圆阵宽带相干信号源方位估计方法 |
CN106802403A (zh) * | 2017-02-22 | 2017-06-06 | 西安电子科技大学 | 声矢量传感器二维阵列music解相干参数估计方法 |
CN106997037A (zh) * | 2017-02-22 | 2017-08-01 | 西安电子科技大学 | 声矢量传感器阵列空间旋转解相干到达角估计方法 |
CN107092006A (zh) * | 2017-04-07 | 2017-08-25 | 广东精点数据科技股份有限公司 | 用于doa估计的阵列信号处理方法和*** |
-
2018
- 2018-02-28 CN CN201810165467.XA patent/CN108919176B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105676168A (zh) * | 2015-12-02 | 2016-06-15 | 江苏科技大学 | 一种声矢量阵方位估计方法 |
CN105510874A (zh) * | 2015-12-15 | 2016-04-20 | 吉林大学 | 一种近场信源多参数联合估计的降维music方法 |
CN106249244A (zh) * | 2016-07-07 | 2016-12-21 | 哈尔滨工程大学 | 一种声矢量圆阵宽带相干信号源方位估计方法 |
CN106802403A (zh) * | 2017-02-22 | 2017-06-06 | 西安电子科技大学 | 声矢量传感器二维阵列music解相干参数估计方法 |
CN106997037A (zh) * | 2017-02-22 | 2017-08-01 | 西安电子科技大学 | 声矢量传感器阵列空间旋转解相干到达角估计方法 |
CN107092006A (zh) * | 2017-04-07 | 2017-08-25 | 广东精点数据科技股份有限公司 | 用于doa估计的阵列信号处理方法和*** |
Non-Patent Citations (2)
Title |
---|
Research on acoustic vector sensor array MUSIC algorithm based on pretreatment of information in the field of array unit;Liang Guolong et al.;《2012 IEEE》;20121231;第248-251、257页 * |
基于矢量阵的自初始化MUSIC方位估计算法;张揽月等;《哈尔滨工程大学学报》;20060430;第27卷(第2期);第152-154页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108919176A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107015191B (zh) | 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法 | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN106526530B (zh) | 基于传播算子的2-l型阵列二维doa估计算法 | |
CN106950529B (zh) | 声矢量近场源esprit和music参数估计方法 | |
CN103267966B (zh) | 一种单个二维压差式矢量水听器方位估计的方法 | |
CN1664611A (zh) | 一种基于非直线天线阵列的无源通道校正方法 | |
CN108469599B (zh) | 一种声矢量传感器幅度加权music测向方法 | |
Zhao et al. | Open‐Lake Experimental Investigation of Azimuth Angle Estimation Using a Single Acoustic Vector Sensor | |
Liu et al. | Efficient 2-D DOA estimation for coherent sources with a sparse acoustic vector-sensor array | |
CN106249225B (zh) | 稀疏圆形声矢量传感器阵列四元数esprit参数估计方法 | |
CN106526531A (zh) | 基于三维天线阵列的改进传播算子二维doa估计算法 | |
CN108919176B (zh) | 一种单矢量传感器增秩music测向技术 | |
CN109696657A (zh) | 一种基于矢量水听器的相干声源定位方法 | |
CN107037398A (zh) | 一种二维music算法估计波达方向的并行计算方法 | |
CN104020440B (zh) | 基于l型干涉式线性阵列的二维波达角估计方法 | |
CN106249196A (zh) | 三分量声矢量传感器稀疏阵列四元数解模糊方法 | |
Zheng et al. | Fast DOA estimation of incoherently distributed sources by novel propagator | |
CN112285647A (zh) | 一种基于稀疏表示与重构的信号方位高分辨估计方法 | |
Lee et al. | Algebraic path-following algorithm for localising 3-D near-field sources in uniform circular array | |
CN108318855A (zh) | 基于均匀圆阵的近场和远场混合信号源定位方法 | |
CN103323810A (zh) | 一种l阵方位角和俯仰角配对的信号处理方法 | |
CN109061594B (zh) | 一种强杂波环境下基于子阵空时自适应处理的测角方法 | |
CN111157951B (zh) | 一种基于差分麦克风阵列的三维声源定位方法 | |
CN111431575B (zh) | 基于常规波束形成的来波方向稀疏重构方法 | |
CN114239251A (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 |