CN109283492A - 多目标方位估计方法及水声垂直矢量阵列系统 - Google Patents

多目标方位估计方法及水声垂直矢量阵列系统 Download PDF

Info

Publication number
CN109283492A
CN109283492A CN201811264949.7A CN201811264949A CN109283492A CN 109283492 A CN109283492 A CN 109283492A CN 201811264949 A CN201811264949 A CN 201811264949A CN 109283492 A CN109283492 A CN 109283492A
Authority
CN
China
Prior art keywords
target
spectrum
vector
array
frequency point
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
Application number
CN201811264949.7A
Other languages
English (en)
Other versions
CN109283492B (zh
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.)
CETC 3 Research Institute
Original Assignee
CETC 3 Research Institute
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 CETC 3 Research Institute filed Critical CETC 3 Research Institute
Priority to CN201811264949.7A priority Critical patent/CN109283492B/zh
Publication of CN109283492A publication Critical patent/CN109283492A/zh
Application granted granted Critical
Publication of CN109283492B publication Critical patent/CN109283492B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
    • G01S5/20Position of source determined by a plurality of spaced direction-finders

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空间谱,对多目标对应的空间谱谱峰进行二维搜索判决,分别得到多目标的水平角和俯仰角;步骤六:对各线谱计算出的水平角和俯仰角进行整合,得到多目标方位估计结果。本发明通过三维垂直矢量水听器阵列实现了对目标信号水平角和俯仰角的全方位估计,提高了相干多目标信号的方位估计精确性。

Description

多目标方位估计方法及水声垂直矢量阵列系统
技术领域
本发明属于海洋工程技术领域,尤其涉及一种多目标方位估计方法及水声垂直矢量阵列系统。
背景技术
一个矢量水听器可以同时测量声场中的声压标量和三个正交方向的质点振速矢量。矢量信息的获得,加深了人们对声场特性的全面认识。与传统的声压水听器阵列相比,矢量水听器阵列具有更多的处理手段,可以获得更大的处理增益。现有水声目标方位估计方法通常利用二维矢量水听器阵列,所采用的算法仅对水平角度起到了较为理想的效果,不能保证多目标信号的方位估计精确性。
发明内容
本发明的目的在于提供一种多目标方位估计方法及水声垂直矢量阵列系统,通过三维垂直矢量水听器阵列实现对目标信号水平角和俯仰角的全方位估计,利用矢量正则相关算法估计目标声源个数,并利用矢量加权空间平滑算法来提高相干多目标信号的方位估计精确性。
本发明提供了一种多目标方位估计方法,包括如下步骤:
步骤一:基于垂直矢量水听器阵列,接收多目标信号;
步骤二:对接收的多目标信号作变换得到频谱和功率谱,并提取线谱;
步骤三:利用矢量正则相关算法,对目标声源个数进行估计;
步骤四:根据提取的线谱,利用矢量加权空间平滑算法,获得MUSIC空间谱;
步骤五:基于所述目标声源个数的估计结果及所述MUSIC空间谱,对多目标对应的空间谱谱峰进行二维搜索判决,分别得到多目标的水平角和俯仰角;
步骤六:对各线谱计算出的水平角和俯仰角进行整合,得到多目标方位估计结果。
进一步地,所述垂直矢量水听器阵列包括M个三维矢量阵元,M为不小于2的自然数,各个三维矢量阵元均匀排列于垂直线列上。
进一步地,步骤二包括:
对所有阵元功率谱整合并提取噪声背景的趋势项,通过设定功率谱门限提取反映目标情况的功率谱频点值和功率谱峰值;
获取M个三维矢量阵元反映多目标情况的功率谱频点值,并获取其相邻间隔对应的复数频谱值,将上述多个复数频谱值进一步整合为4M通道的接收频域数据。
进一步地,步骤三包括:
将功率谱频点值处构造的接收频域数据拆分为阵元数nx和ny的上下两个矢量阵列X和Y形成的数据,nx≦ny,nx+ny=M;
构造满足自由度为m=2(4nx-k)(4ny-k),k=1,2,…,4nx-1的近似x2分布的序列,设定判决门限,序列与对应的门限比较,第1个小于门限的位置数减1即为判决出的目标声源个数。
进一步地,步骤四包括:
将功率谱频点值处构造的接收频域数据拆分为1个声压和3个振速数据,再分别计算4个数据协方差矩阵;
对4个数据协方差矩阵分别进行前后向平滑处理,并引入矢量权重得到矢量加权平滑矩阵;
对矢量加权平滑矩阵进行特征值分解,结合步骤三获取的目标声源个数,计算得到噪声子空间矩阵;
计算垂直矢量阵列的阵列流型,根据阵列的目标信号子空间导向矢量和环境噪声子空间所具有的正交特性得到基于频率点功率谱频点值的空间谱。
本发明还提供了一种水声垂直矢量阵列系统,包括:
垂直矢量水听器阵列,用于接收多目标信号;
滤波放大模块,用于对接收的多目标信号进行带通滤波和信号放大初步处理;
信号处理模块,用于对经滤波及信号放大处理后的多目标信号作变换得到频谱和功率谱,并提取线谱;利用矢量正则相关算法,对目标声源个数进行估计;根据提取的线谱,利用矢量加权空间平滑算法,获得MUSIC空间谱;基于所述目标声源个数的估计结果及所述空间谱,对多目标对应的空间谱谱峰进行二维搜索判决,分别得到多目标的水平角和俯仰角;对各线谱计算出的水平角和俯仰角进行整合,得到多目标方位估计结果;
数据存储模块,用于存储所述多目标方位估计结果;
显示控制终端,用于显示多目标方位估计结果;
进一步地,垂直矢量水听器阵列包括M个三维矢量阵元,M为不小于2 的自然数,各个三维矢量阵元均匀排列于垂直线列上。
进一步地,所述信号处理模块对所有阵元功率谱整合并提取噪声背景的趋势项,通过设定功率谱门限提取反映目标情况的功率谱频点值和功率谱峰值;以及
获取M个三维矢量阵元反映多目标情况的功率谱频点值,并获取其相邻间隔对应的复数频谱值,将上述多个复数频谱值进一步整合为4M通道的接收频域数据。
进一步地,所述信号处理模块将功率谱频点值处构造的接收频域数据拆分为阵元数nx和ny的上下两个矢量阵列X和Y形成的数据,nx≦ny,nx+ny=M;以及
构造满足自由度为m=2(4nx-k)(4ny-k),k=1,2,…,4nx-1的近似x2分布的序列,设定判决门限,序列与对应的门限比较,第1个小于门限的位置数减 1即为判决出的目标声源个数。
进一步地,所述信号处理模块将功率谱频点值处构造的接收频域数据拆分为1个声压和3个振速数据,再分别计算4个数据协方差矩阵;
对4个数据协方差矩阵分别进行前后向平滑处理,并引入矢量权重得到矢量加权平滑矩阵;
对矢量加权平滑矩阵进行特征值分解,结合所述目标声源个数,计算得到噪声子空间矩阵;以及
计算垂直矢量阵列的阵列流型,根据阵列的目标信号子空间导向矢量和环境噪声子空间所具有的正交特性得到基于频率点功率谱频点值的空间谱。
与现有技术相比本发明的有益效果是:
通过三维垂直矢量水听器阵列实现了对目标信号水平角和俯仰角的全方位估计,利用矢量正则相关算法估计目标声源个数,并利用矢量加权空间平滑算法来提高了相干多目标信号的方位估计精确性。
附图说明
图1是本发明垂直矢量阵列示意图;
图2是本发明水声垂直矢量阵列系统的结构框图;
图3是本发明前后向空间平滑算法原理图;其中,图3(a)为前向平滑,图3(b)为后向平滑;
图4是本发明垂直矢量阵列系统信号处理流程图。
具体实施方式
下面结合附图所示的各实施方式对本发明进行详细说明,但应当说明的是,这些实施方式并非对本发明的限制,本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代,均属于本发明的保护范围之内。
本实施例提供了一种多目标方位估计方法,包括如下步骤:
步骤一:基于垂直矢量水听器阵列,接收多目标信号;
步骤二:对接收的多目标信号作变换得到频谱和功率谱,并提取线谱;
步骤三:利用矢量正则相关算法,对目标声源个数进行估计;
步骤四:根据提取的线谱,利用矢量加权空间平滑算法,获得MUSIC空间谱;
步骤五:基于所述目标声源个数的估计结果及所述MUSIC空间谱,对多目标对应的空间谱谱峰进行二维搜索判决,分别得到多目标的水平角和俯仰角;
步骤六:对各线谱计算出的水平角和俯仰角进行整合,得到多目标方位估计结果。
该多目标方位估计方法通过三维垂直矢量水听器阵列实现了对目标信号水平角和俯仰角的全方位估计,利用矢量正则相关算法估计目标声源个数,并利用矢量加权空间平滑算法来提高了相干多目标信号的方位估计精确性。
在本实施例中,垂直矢量水听器阵列包括M个三维矢量阵元,M为不小于2的自然数,各个三维矢量阵元均匀排列于垂直线列上。
在本实施例中,步骤二包括:
对所有阵元功率谱整合并提取噪声背景的趋势项,通过设定功率谱门限提取反映目标情况的功率谱频点值和功率谱峰值;
获取M个三维矢量阵元反映多目标情况的功率谱频点值,并获取其相邻间隔对应的复数频谱值,将上述多个复数频谱值进一步整合为4M通道的接收频域数据。
在本实施例中,步骤三包括:
将功率谱频点值处构造的接收频域数据拆分为阵元数nx和ny的上下两个矢量阵列X和Y形成的数据,nx≦ny,nx+ny=M;
构造满足自由度为m=2(4nx-k)(4ny-k),k=1,2,…,4nx-1的近似x2分布的序列,设定判决门限,序列与对应的门限比较,第1个小于门限的位置数减 1即为判决出的目标声源个数。
在本实施例中,步骤四包括:
将功率谱频点值处构造的接收频域数据拆分为1个声压和3个振速数据,再分别计算4个数据协方差矩阵;
对4个数据协方差矩阵分别进行前后向平滑处理,并引入矢量权重得到矢量加权平滑矩阵;
对矢量加权平滑矩阵进行特征值分解,结合步骤三获取的目标声源个数,计算得到噪声子空间矩阵;
计算垂直矢量阵列的阵列流型,根据阵列的目标信号子空间导向矢量和环境噪声子空间所具有的正交特性得到基于频率点功率谱频点值的空间谱。
本实施例还提供了一种水声垂直矢量阵列系统,包括:
垂直矢量水听器阵列,用于接收多目标信号;
滤波放大模块,用于对接收的多目标信号进行带通滤波和信号放大初步处理;
信号处理模块,用于对经滤波及信号放大处理后的多目标信号作变换得到频谱和功率谱,并提取线谱;利用矢量正则相关算法,对目标声源个数进行估计;根据提取的线谱,利用矢量加权空间平滑算法,获得MUSIC空间谱;基于所述目标声源个数的估计结果及所述空间谱,对多目标对应的空间谱谱峰进行二维搜索判决,分别得到多目标的水平角和俯仰角;对各线谱计算出的水平角和俯仰角进行整合,得到多目标方位估计结果。;
数据存储模块,用于存储多目标方位估计结果;
显示控制终端,用于显示多目标方位估计结果。
该水声垂直矢量阵列系统通过三维垂直矢量水听器阵列实现了对目标信号水平角和俯仰角的全方位估计,利用矢量正则相关算法估计目标声源个数,并利用矢量加权空间平滑算法来提高了相干多目标信号的方位估计精确性。
在本实施例中,垂直矢量水听器阵列包括M个三维矢量阵元,M为不小于2的自然数,各个三维矢量阵元均匀排列于垂直线列上。
在本实施例中,信号处理模块对所有阵元功率谱整合并提取噪声背景的趋势项,通过设定功率谱门限提取反映目标情况的功率谱频点值和功率谱峰值;以及
获取M个三维矢量阵元反映多目标情况的功率谱频点值,并获取其相邻间隔对应的复数频谱值,将上述多个复数频谱值进一步整合为4M通道的接收频域数据。
在本实施例中,信号处理模块将功率谱频点值处构造的接收频域数据拆分为阵元数nx和ny的上下两个矢量阵列X和Y形成的数据,nx≦ny,nx+ny=M;以及
构造满足自由度为m=2(4nx-k)(4ny-k),k=1,2,…,4nx-1的近似x2分布的序列,设定判决门限,序列与对应的门限比较,第1个小于门限的位置数减 1即为判决出的目标声源个数。
在本实施例中,信号处理模块将功率谱频点值处构造的接收频域数据拆分为1个声压和3个振速数据,再分别计算4个数据协方差矩阵;
对4个数据协方差矩阵分别进行前后向平滑处理,并引入矢量权重得到矢量加权平滑矩阵;
对矢量加权平滑矩阵进行特征值分解,结合所述目标声源个数,计算得到噪声子空间矩阵;以及
计算垂直矢量阵列的阵列流型,根据阵列的目标信号子空间导向矢量和环境噪声子空间所具有的正交特性得到基于频率点功率谱频点值的空间谱。
下面对本发明作进一步详细说明。
参图2所示,水声垂直矢量阵列系统,包括:
(1)M个三维矢量阵元,各三维矢量阵元均匀排列于间距为d垂直线列上。其中,M为不小于2的自然数。如图1所示。
(2)滤波放大和数字采集模块,每个三维矢量阵元均配有声压振速4通道的滤波放大模块,用于对接收的三维矢量信号进行带通滤波和信号放大初步处理。经放大滤波初步处理后的信号,再利用数字采集模块进行多通道数字采集成为数字信号,为后续信号处理做好准备。
(3)信号处理模块
负责对数字信号进行谱分析,线谱提取,空间谱估计和目标方位估计等运算,运算结果为数据存储模块和显示控制终端提供输入。
(4)数据存储模块
与信号处理模块和显示控制终端相连,负责将信号处理模块输出结果进行存储,用于结果分析和回放。
(5)显示控制终端
分别与数据存储模块和信号处理模块,负责显示和回放多目标方位估计结果。
基于上述水声垂直矢量阵列系统的高精度多目标方位估计方法,包括如下步骤:
(1)以第1个三维矢量阵元为坐标原点,该阵元的三个振速方向为坐标轴,建立空间参考坐标系。如图1所示。考虑K个目标,每个目标发出某一窄带信号,所有信号均被M元矢量阵列接收。于是矢量阵列的接收信号为:
S(t)=A(θ)·E(t)+N(t) (1)
在此,E(t)=[e1(t),e2(t),...eK(t)]T为K个目标的窄带信号矢量, A(θ)=[a(θ1),a(θ2),...a(θK)]为目标信号方向矢量矩阵;为A(θ)的方向矢量矩阵元素,为克罗内克直积。间距为dm的各阵元的矢量矩阵ack)在远场等效为:
ack)=[1,exp(-j2πfkdmsin(pitk)/c),... ,exp(-j2πfkdmsin(pitk)/c)]T (2)
目标单位方向矢量为:
u(θk)=[1,cos(pitk)cos(azik),cos(pitk)sin(azik),sin(pitk)]T (3)
其中,azik和pitk分别表示第k个目标的水平角和俯仰角;
公式(1)中代表了海洋环境有色噪声的矢量矩阵为:
N(t)=[np(t),nvx(t),nvy(t),nvz(t)]T (4)
其中,N(t)的元素分别为:
np(t)=[np1(t),np2(t),...,npM(t)]T
nvx(t)=[nvx1(t),nvx2(t),...,nvxM(t)]T
nvy(t)=[nvy1(t),nvy2(t),...,nvyM(t)]T
nvz(t)=[nvz1(t),nvz2(t),...,nvzM(t)]T (5)
进一步考虑K个目标,每个目标最多发出I个窄带谐波信号,所有信号均被M元矢量阵列接收。于是矢量阵列的接收信号为:
(2)对M元矢量阵列接收信号S1(t)作频谱F(f)和功率谱P(f):
F(f)=fft(S1(t)),P(f)=psd(S1(t)) (7)
对所有阵元P(f)整合并提取噪声背景的趋势项,然后通过设定功率谱门限提取反映目标情况的num个功率谱频点值fnum和功率谱峰值pnum
为保证后续算法数据的可靠性,获取M个矢量阵元的频谱F(f)上每一个频率点fnum附近的复数频谱值。设定finterval为频率点fnum相邻间隔的频率点,则可以将获取的复数频谱值进一步整合为4M通道的接收频域数据:
data=F(fnum-finterval)+F(fnum)+F(fnum+finterval) (8)
(3)将频率点fnum处构造的接收数据data拆分为阵元数nx和ny的上下两个矢量阵列X和Y形成的数据,nx≦ny,nx+ny=M。由于在色噪声背景下声矢量阵声压和振速通道噪声有一定程度的相关性,在此采用矢量阵列的正则相关算法对目标声源个数进行估计,此时:
计算两者的自相关矩阵Rxx,Ryy和互相关矩阵Rxy,构造矩阵Rsvd并进行奇异值分解:
Rsvd=Rxx -1/2Rxy(Ryy -1/2)H=UTV (10)
上式中U、V分别是左右奇异矩阵,即为奇异值对应由大到小排列后组成的对角阵,且有如下关系:
上式中的奇异值也称为正则相关因子。正则相关因子在特定条件下满足Bartletts近似,即当接收数据满足高斯分布,且满足时,构造的如下序列Seq(k)满足自由度为m=2(4nx-k)(4ny-k)的近似χ2分布:
设定判决门限Tk,数列与对应的门限比较,第1个小于门限的位置数减1 即为判决出的目标声源个数。该判决门限Tk通过设定下式的虚警概率Pf获得:
(4)将频率点fnum处构造的接收数据data拆分为1个声压和3个振速数据,再分别计算4个数据协方差矩阵。
将Rp,Rvx,Rvy,Rvz对应的等间距垂直阵列分成P个重叠的子阵列,每个子阵列的阵元数为Q,K为入射信源数,子阵列从上逐步下移,如图3所示。很明显, 子阵阵元数Q和子阵列总个数P满足:
P+Q-1=M (16)
将Rp,Rvx,Rvy,Rvz分别进行上述的前后向平滑处理,得到以Rp为例,R′p为元素共轭,则有:
基于上述处理,考虑引入矢量权重wp,wvx,wvy,wvz后的矢量加权空间平滑矩阵Rxx为:
其中,wp,wvx,wvy,wvz∈(0,1),且和为1。
(5)采用MUSIC方法计算空间谱,从而实现基于频率点fnum的多目标空间谱计算。对式(19)进行特征值分解,再结合步骤(3)获取的目标源数,可以计算得到噪声子空间矩阵Un
计算垂直矢量阵列的阵列流型A(θ),根据阵列的目标信号子空间导向矢量 A(θ)和环境噪声子空间Un所具有的正交特性得到基于频率点fnum的空间谱PMUSIC
(6)基于步骤(3)和(5)的结果,对多目标对应的空间谱谱峰进行二维搜索判决,分别得到多目标的水平角azis和俯仰角pits。
[P1,P2,...,Pk]=max(PMUSIC),1≤k≤K
(azisk,pitsk)=position(Pk) (22)
整合基于所有频率点fnum的多目标的水平角azis和俯仰角pits,最终实现垂直矢量阵列的相干高精度全方位估计。
基于上述多目标方位估计方法的垂直矢量阵列系统信号处理流程参图4 所示。
本实施例提供的水声多方位目标信号估计方法及水声垂直矢量阵列系统具有如下技术效果:
(1)实现了三维垂直矢量水听器阵列对多目标信号水平角和俯仰角的全方位估计。
(2)实现了高精度方位估计方法在垂直矢量水听器阵列系统上的应用。
(3)水声垂直矢量阵列系统和高精度方位估计方法对色噪声背景下的相干多目标方位估计给出了解决方案,为技术研究和更深入的工程应用提供了基础。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。

Claims (10)

1.一种多目标方位估计方法,其特征在于,包括如下步骤:
步骤一:基于垂直矢量水听器阵列,接收多目标信号;
步骤二:对接收的多目标信号作变换得到频谱和功率谱,并提取线谱;
步骤三:利用矢量正则相关算法,对目标声源个数进行估计;
步骤四:根据提取的线谱,利用矢量加权空间平滑算法,获得MUSIC空间谱;
步骤五:基于所述目标声源个数的估计结果及所述MUSIC空间谱,对多目标对应的空间谱谱峰进行二维搜索判决,分别得到多目标的水平角和俯仰角;
步骤六:对各线谱计算出的水平角和俯仰角进行整合,得到多目标方位估计结果。
2.根据权利要求1所述的多目标方位估计方法,其特征在于,所述垂直矢量水听器阵列包括M个三维矢量阵元,M为不小于2的自然数,各个三维矢量阵元均匀排列于垂直线列上。
3.根据权利要求2所述的多目标方位估计方法,其特征在于,所述步骤二包括:
对所有阵元功率谱整合并提取噪声背景的趋势项,通过设定功率谱门限提取反映目标情况的功率谱频点值和功率谱峰值;
获取M个三维矢量阵元上反映多目标情况的功率谱频点值,并获取其相邻间隔对应的复数频谱值,将上述多个复数频谱值进一步整合为4M通道的接收频域数据。
4.根据权利要求3所述的多目标方位估计方法,其特征在于,所述步骤三包括:
将功率谱频点值处构造的接收频域数据拆分为阵元数nx和ny的上下两个矢量阵列X和Y形成的数据,nx≦ny,nx+ny=M;
构造满足自由度为m=2(4nx-k)(4ny-k),k=1,2,…,4nx-1的近似x2分布的序列,设定判决门限,序列与对应的门限比较,第1个小于门限的位置数减1即为判决出的目标声源个数。
5.根据权利要求4所述的多目标方位估计方法,其特征在于,所述步骤四包括:
将功率谱频点值处构造的接收频域数据拆分为1个声压和3个振速数据,再分别计算4个数据协方差矩阵;
对4个数据协方差矩阵分别进行前后向平滑处理,并引入矢量权重得到矢量加权平滑矩阵;
对矢量加权平滑矩阵进行特征值分解,结合步骤三获取的目标声源个数,计算得到噪声子空间矩阵;
计算垂直矢量阵列的阵列流型,根据阵列的目标信号子空间导向矢量和环境噪声子空间所具有的正交特性得到基于频率点功率谱频点值的空间谱。
6.一种水声垂直矢量阵列系统,其特征在于,包括:
垂直矢量水听器阵列,用于接收多目标信号;
滤波放大模块,用于对接收的多目标信号进行带通滤波和信号放大初步处理;
信号处理模块,用于对经滤波及信号放大处理后的多目标信号作变换得到频谱和功率谱,并提取线谱;利用矢量正则相关算法,对目标声源个数进行估计;根据提取的线谱,利用矢量加权空间平滑算法,获得MUSIC空间谱;基于所述目标声源个数的估计结果及所述空间谱,对多目标对应的空间谱谱峰进行二维搜索判决,分别得到多目标的水平角和俯仰角;对各线谱计算出的水平角和俯仰角进行整合,得到多目标方位估计结果;
数据存储模块,用于存储所述多目标方位估计结果;
显示控制终端,用于显示多目标方位估计结果。
7.根据权利要求6所述的水声垂直矢量阵列系统,其特征在于,垂直矢量水听器阵列包括M个三维矢量阵元,M为不小于2的自然数,各个三维矢量阵元均匀排列于垂直线列上。
8.根据权利要求7所述的水声垂直矢量阵列系统,其特征在于,所述信号处理模块对所有阵元功率谱整合并提取噪声背景的趋势项,通过设定功率谱门限提取反映目标情况的功率谱频点值和功率谱峰值;以及
获取M个三维矢量阵元上反映多目标情况的功率谱频点值,并获取其相邻间隔对应的复数频谱值,将上述多个复数频谱值进一步整合为4M通道的接收频域数据。
9.根据权利要求8所述的水声垂直矢量阵列系统,其特征在于,所述信号处理模块将功率谱频点值处构造的接收频域数据拆分为阵元数nx和ny的上下两个矢量阵列X和Y形成的数据,nx≦ny,nx+ny=M;以及
构造满足自由度为m=2(4nx-k)(4ny-k),k=1,2,…,4nx-1的近似x2分布的序列,设定判决门限,序列与对应的门限比较,第1个小于门限的位置数减1即为判决出的目标声源个数。
10.根据权利要求9所述的水声垂直矢量阵列系统,其特征在于,所述信号处理模块将功率谱频点值处构造的接收频域数据拆分为1个声压和3个振速数据,再分别计算4个数据协方差矩阵;
对4个数据协方差矩阵分别进行前后向平滑处理,并引入矢量权重得到矢量加权平滑矩阵;
对矢量加权平滑矩阵进行特征值分解,结合所述目标声源个数,计算得到噪声子空间矩阵;以及
计算垂直矢量阵列的阵列流型,根据阵列的目标信号子空间导向矢量和环境噪声子空间所具有的正交特性得到基于频率点功率谱频点值的空间谱。
CN201811264949.7A 2018-10-29 2018-10-29 多目标方位估计方法及水声垂直矢量阵列系统 Active CN109283492B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811264949.7A CN109283492B (zh) 2018-10-29 2018-10-29 多目标方位估计方法及水声垂直矢量阵列系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811264949.7A CN109283492B (zh) 2018-10-29 2018-10-29 多目标方位估计方法及水声垂直矢量阵列系统

Publications (2)

Publication Number Publication Date
CN109283492A true CN109283492A (zh) 2019-01-29
CN109283492B CN109283492B (zh) 2021-02-19

Family

ID=65178055

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811264949.7A Active CN109283492B (zh) 2018-10-29 2018-10-29 多目标方位估计方法及水声垂直矢量阵列系统

Country Status (1)

Country Link
CN (1) CN109283492B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110390073A (zh) * 2019-08-19 2019-10-29 西北工业大学 一种矢量传感的多通道空间合成方位滤波方法
CN110389319A (zh) * 2019-07-22 2019-10-29 北京工业大学 一种基于低空多径情况下的mimo雷达doa估计方法
CN110716177A (zh) * 2019-10-22 2020-01-21 哈尔滨工程大学 分布式声学定位网络中多目标分类方法
CN110764053A (zh) * 2019-10-22 2020-02-07 浙江大学 一种基于水下传感器网络的多目标被动定位方法
CN110837076A (zh) * 2019-11-18 2020-02-25 哈尔滨工程大学 一种基于张量分解的矢量水听器阵列方位估计方法
CN111273263A (zh) * 2019-06-25 2020-06-12 哈尔滨工程大学 基于信息融合的自主探测声纳多目标doa估计方法
CN111323750A (zh) * 2020-03-19 2020-06-23 哈尔滨工程大学 一种基于声矢量阵列网络的直接定位方法
CN111665489A (zh) * 2019-03-08 2020-09-15 中国科学院声学研究所 一种基于目标特性的线谱提取方法
CN112033520A (zh) * 2020-09-07 2020-12-04 中国电子科技集团公司第三研究所 一种水声矢量探测系统的校正系统及空气声学校正方法
CN112071297A (zh) * 2020-09-07 2020-12-11 西北工业大学 一种矢量声的自适应滤波方法
CN112528216A (zh) * 2020-12-03 2021-03-19 昆明理工大学 基于声强矢量分解的多点声源声功率频率谱的分离方法
CN112924926A (zh) * 2021-01-27 2021-06-08 东南大学 一种多途信道下的doa估计方法
CN113376579A (zh) * 2021-06-23 2021-09-10 华中科技大学 基于分布式声波传感技术的水下声源定位及信号获取方法
CN113391257A (zh) * 2020-03-13 2021-09-14 光宝科技新加坡私人有限公司 用于物件的角度估测的运算装置及物件的角度估测方法
CN113624330A (zh) * 2021-07-12 2021-11-09 哈尔滨工程大学 一种水下目标辐射噪声测量组合体积阵及测量方法
CN116299156A (zh) * 2023-05-26 2023-06-23 中国海洋大学 一种水听器垂直阵阵元位置估计方法及其优化策略

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6901362B1 (en) * 2000-04-19 2005-05-31 Microsoft Corporation Audio segmentation and classification
CN101034158A (zh) * 2007-02-25 2007-09-12 四川川大智胜软件股份有限公司 基于麦克风阵联网的低空目标监视方法
CN102663423A (zh) * 2012-03-28 2012-09-12 北京航空航天大学 一种简谱图像的自动识别和演奏的方法
CN104749563A (zh) * 2015-03-26 2015-07-01 武汉大学 从高频地波雷达海洋回波一阶Bragg谐振峰中提取浪高的方法
CN105093070A (zh) * 2014-05-06 2015-11-25 长沙理工大学 一种大型变压器多放电源超声定位方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6901362B1 (en) * 2000-04-19 2005-05-31 Microsoft Corporation Audio segmentation and classification
CN101034158A (zh) * 2007-02-25 2007-09-12 四川川大智胜软件股份有限公司 基于麦克风阵联网的低空目标监视方法
CN102663423A (zh) * 2012-03-28 2012-09-12 北京航空航天大学 一种简谱图像的自动识别和演奏的方法
CN105093070A (zh) * 2014-05-06 2015-11-25 长沙理工大学 一种大型变压器多放电源超声定位方法
CN104749563A (zh) * 2015-03-26 2015-07-01 武汉大学 从高频地波雷达海洋回波一阶Bragg谐振峰中提取浪高的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
井岩: "基于压缩感知的水声信号波达方向估计方法研究", 《中国博士学位论文全文数据库 信息科技辑》 *
宫君乐,等: "基于VS2010的水声垂直矢量阵列声呐程序设计与实现", 《电声技术》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111665489A (zh) * 2019-03-08 2020-09-15 中国科学院声学研究所 一种基于目标特性的线谱提取方法
CN111273263A (zh) * 2019-06-25 2020-06-12 哈尔滨工程大学 基于信息融合的自主探测声纳多目标doa估计方法
CN111273263B (zh) * 2019-06-25 2022-07-12 哈尔滨工程大学 基于信息融合的自主探测声纳多目标doa估计方法
CN110389319A (zh) * 2019-07-22 2019-10-29 北京工业大学 一种基于低空多径情况下的mimo雷达doa估计方法
CN110389319B (zh) * 2019-07-22 2021-04-27 北京工业大学 一种基于低空多径情况下的mimo雷达doa估计方法
CN110390073A (zh) * 2019-08-19 2019-10-29 西北工业大学 一种矢量传感的多通道空间合成方位滤波方法
CN110390073B (zh) * 2019-08-19 2023-03-24 西北工业大学 一种矢量传感的多通道空间合成方位滤波方法
CN110716177B (zh) * 2019-10-22 2022-09-09 哈尔滨工程大学 分布式声学定位网络中多目标分类方法
CN110764053B (zh) * 2019-10-22 2021-08-17 浙江大学 一种基于水下传感器网络的多目标被动定位方法
CN110764053A (zh) * 2019-10-22 2020-02-07 浙江大学 一种基于水下传感器网络的多目标被动定位方法
CN110716177A (zh) * 2019-10-22 2020-01-21 哈尔滨工程大学 分布式声学定位网络中多目标分类方法
CN110837076A (zh) * 2019-11-18 2020-02-25 哈尔滨工程大学 一种基于张量分解的矢量水听器阵列方位估计方法
CN113391257A (zh) * 2020-03-13 2021-09-14 光宝科技新加坡私人有限公司 用于物件的角度估测的运算装置及物件的角度估测方法
CN111323750B (zh) * 2020-03-19 2022-08-02 哈尔滨工程大学 一种基于声矢量阵列网络的直接定位方法
CN111323750A (zh) * 2020-03-19 2020-06-23 哈尔滨工程大学 一种基于声矢量阵列网络的直接定位方法
CN112071297A (zh) * 2020-09-07 2020-12-11 西北工业大学 一种矢量声的自适应滤波方法
CN112033520A (zh) * 2020-09-07 2020-12-04 中国电子科技集团公司第三研究所 一种水声矢量探测系统的校正系统及空气声学校正方法
CN112071297B (zh) * 2020-09-07 2023-11-10 西北工业大学 一种矢量声的自适应滤波方法
CN112528216A (zh) * 2020-12-03 2021-03-19 昆明理工大学 基于声强矢量分解的多点声源声功率频率谱的分离方法
CN112924926A (zh) * 2021-01-27 2021-06-08 东南大学 一种多途信道下的doa估计方法
CN113376579A (zh) * 2021-06-23 2021-09-10 华中科技大学 基于分布式声波传感技术的水下声源定位及信号获取方法
CN113624330A (zh) * 2021-07-12 2021-11-09 哈尔滨工程大学 一种水下目标辐射噪声测量组合体积阵及测量方法
CN113624330B (zh) * 2021-07-12 2023-11-17 哈尔滨工程大学 一种水下目标辐射噪声测量组合体积阵及测量方法
CN116299156A (zh) * 2023-05-26 2023-06-23 中国海洋大学 一种水听器垂直阵阵元位置估计方法及其优化策略
CN116299156B (zh) * 2023-05-26 2023-08-15 中国海洋大学 一种水听器垂直阵阵元位置估计方法及其优化策略

Also Published As

Publication number Publication date
CN109283492B (zh) 2021-02-19

Similar Documents

Publication Publication Date Title
CN109283492A (zh) 多目标方位估计方法及水声垂直矢量阵列系统
CN110458055B (zh) 一种障碍物检测方法及系统
CN108226892B (zh) 一种基于深度学习的复杂噪声环境下的雷达信号恢复方法
CN104237879B (zh) 一种雷达系统中的多目标跟踪方法
CN101661104B (zh) 基于雷达/红外量测数据坐标转换的目标跟踪方法
CN106019254B (zh) 一种uwb冲击生物雷达多人体目标距离向分离辨识方法
CN107202989B (zh) 一种适用于被动拖曳线列阵声呐的复杂弱目标检测和跟踪方法
CN110018461B (zh) 基于高分辨距离像和单脉冲测角的群目标识别方法
CN109444896B (zh) 水声阵列定位系统及其定位方法
CN110058222B (zh) 一种基于传感器选择的双层粒子滤波检测前跟踪方法
CN102213759A (zh) 基于功率谱的水声目标特征匹配方法
CN106500671A (zh) 一种基于lm算法分解激光雷达波形确定海水深度的方法
CN108318876A (zh) 一种利用单个水听器估计水下目标深度和距离的方法
CN107390229A (zh) 一种测风激光雷达信号的处理方法及其处理装置
CN109884591A (zh) 一种基于麦克风阵列的多旋翼无人机声信号增强方法
CN106546949A (zh) 一种基于频率预估计的双阵元正弦信号来波方向估计方法
CN105022046A (zh) 一种基于图像特征的雷达弱目标检测方法
CN106546947A (zh) 一种联合波导不变量和线谱的单水听器被动定位方法
Janakaraj et al. STAR: Simultaneous tracking and recognition through millimeter waves and deep learning
CN116338684A (zh) 基于毫米波雷达与深度学习的人体跌倒检测方法及系统
CN105676169A (zh) 一种脉冲声目标的定向方法及装置
CN105445724B (zh) 单水听器自由场被动测距方法
US20240012083A1 (en) Method and apparatus for measuring directions of arrival of multiple sound sources
Zhang et al. Underwater target detection based on machine learning
CN110334674A (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