CN106066472A - 一种二维振速梯度水听器的被动目标相关检测方法 - Google Patents

一种二维振速梯度水听器的被动目标相关检测方法 Download PDF

Info

Publication number
CN106066472A
CN106066472A CN201510606683.XA CN201510606683A CN106066472A CN 106066472 A CN106066472 A CN 106066472A CN 201510606683 A CN201510606683 A CN 201510606683A CN 106066472 A CN106066472 A CN 106066472A
Authority
CN
China
Prior art keywords
vibration velocity
velocity gradient
alpha
gradient
target
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
CN201510606683.XA
Other languages
English (en)
Other versions
CN106066472B (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.)
Shandong Beiming Technology Co ltd
Original Assignee
Qindao University Of Technology
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 Qindao University Of Technology filed Critical Qindao University Of Technology
Priority to CN201510606683.XA priority Critical patent/CN106066472B/zh
Publication of CN106066472A publication Critical patent/CN106066472A/zh
Application granted granted Critical
Publication of CN106066472B publication Critical patent/CN106066472B/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

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)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种二维振速梯度水听器的被动目标相关检测方法。该方法利用振速水听器的接收数据,先计算出中心点处的振速和振速梯度,后对二者进行互相关处理,利用互相关处理器的输出来检测是否存在被动目标,并根据不同角度上检测概率的分布粗略估计目标的方位。有益效果是:本发明充分利用了振速和振速梯度的指向性,以及噪声及目标信号振速与振速梯度信息相关性的差异性特点,降低了被动目标的可检测信噪比,从而有效的提高被动目标的探测距离。本发明方法对水下弱目标的远程探测与参数估计均有重要意义,可以在浮标、潜标、鱼探仪、水下航行器、水下定位与导航等产品中应用。

Description

一种二维振速梯度水听器的被动目标相关检测方法
技术领域
本发明涉及的是一种二维振速梯度水听器的被动目标相关检测方法,具体的说,是一种利用二维振速梯度水听器中心点处的振速和振速梯度信息,通过互相关处理来检测是否存在被动目标的方法。
背景技术
随着科学的发展、技术的进步,军用、民用海洋装备对小孔径、远程、高精度的水下探测设备需求增加,这使得水听器检测技术不断发展进步,也促使新的水下传感技术发展,如近十几年来快速发展的矢量水听器技术。矢量水听器可同时共点的获取声场的声压(标量)、振速(矢量)信息。振速具有与频率无关的偶极子指向性,且声压和振速信息联合处理可以有效地抑制噪声,提高被动目标的探测能力。在低频段,矢量水听器相比于声压水听器阵具有体积小,重量轻,便于布放等优点,广泛应用于军事探潜、海洋生物研究、地震监测和环境噪声测量等系统中。
为了得到指向性更加尖锐的传感器,很多学者对振速梯度水听器开展了系列研究,已经通过理论分析和实场实验证明了振速梯度水听器可以获得更好的指向性。如何利用振速梯度水听器实现被动目标的远程检测和有效参数估计,还有待开展进一步的研究。各向同性噪声的振速和振速梯度是不相关的,而目标信号的振速与振速梯度是相干的。如果有一种处理方法或技术,能够将振速与振速梯度信息联合处理,则可提高振速梯度水听器的检测性能,有效增强其探测距离,这种技术必然会扩展振速梯度水听器的应用范围,特别是在水下小孔径测量和监视系统中的应用。
发明内容
本发明的目的是:利用目标信号和环境噪声信号振速与振速梯度相关信息的差异性,提供一种简单实用、稳健可靠、可检测信噪比低的二维振速梯度水听器被动目标的相关检测方法。
本发明的技术方案是:一种二维振速梯度水听器被动目标的相关检测方法,包括如下的步骤:
步骤一:将二维振速梯度水听器接收的八个通道实数信号数据vxi、vyi(i=1,…,4),通过希尔伯特变换器转换成复数信号
步骤二:对四个通道的x轴复振速信号求平均,求取振速梯度水听器中心点处x轴方向振速
步骤三:对四个通道的y轴复振速信号求平均,求取振速梯度水听器中心点处y轴方向振速
步骤四:利用八通道振速复信号求振速梯度水听器中心点处横向振速梯度和纵向振速梯度
步骤五:利用引导方向α,对振速梯度水听器中心点处振速分量进行加权处理,得到合成振速 v ~ c ( α ) = v ~ o x c o s α + v ~ o y s i n α ;
步骤六:利用引导方向α,对振速梯度水听器中心点处横向振速梯度 和纵向振速梯度进行加权处理,得到合成振速梯度 d ~ c ( α ) = ∂ v ~ x / ∂ xcos 2 α + ∂ v ~ y / ∂ ysin 2 α + ( ∂ v ~ x / ∂ y + ∂ v ~ y / ∂ x ) cos α sin α ;
步骤七:对M个采样点的合成振速与合成振速梯度进行相关处理,得到被动目标的检测统计量 D ( α , M ) = Σ n = 1 M | | v ~ c ( α , n ) · d ~ c ( α , n ) | | ;
步骤八:根据系统要求的最大虚警概率和上述检验统计量D(α,M),在噪声环境中计算出不同引导方位上的检测门限Dth(α,M);
步骤九:实时计算D(α,M),将D(α,M)与Dth(α,M)进行比较,如果D(α,M)>Dth(α,M),则目标存在,反之目标不存在,并在不同的α上进行扫描,根据成功检测概率的分布粗略估计目标的方位区域。
本发明的有益效果是:本发明提供了一种稳健、低信噪比的振速梯度水听器被动目标检测方法,根据水听器的接收数据,计算出中心点处组合振速和组合振速梯度,对二者进行互相关处理,利用互相关处理器的输出来检测被动弱目标有无,并粗略估计目标的方位区域,对水声探测与通信均有重要意义。
本发明充分利用了振速和振速梯度的指向性,以及噪声和目标信号振速与振速梯度相关信息的差异性特点,降低了被动目标的可检测信噪比,从而有效的提高被动目标的探测距离。本发明方法对水下弱目标检测与参数估计具有重要意义,可以在浮标、潜标、鱼探仪、水下航行器、水下定位与导航等产品中应用。
附图说明
图1是本发明提出的二维振速梯度水听器被动目标检测的信号处理流程图。
图2是二维振速梯度水听器的结构示意图。
图3为实施例中目标信号中心频率f0=1kHz,系统采样频率为fs=8f0,噪声为数值模拟的各向同性噪声,系统地虚警概率Pf=0.05,目标在90°方位时,振速梯度水听器振速与振速梯度互相关检测器检测成功概率统计结果图。
图4为实施例中目标信号中心频率f0=1kHz,系统采样频率为fs=8f0,噪声为数值模拟的各向同性噪声,系统地虚警概率Pf=0.05,目标在45°方位时,振速梯度水听器振速与振速梯度互相关检测器检测成功概率统计结果图。
具体实施方式
现结合实施例、附图对本发明做进一步描述:
如图1所示,本发明利用四个对角振速传感器的平均和差分运算,计算出中心点处的振速和振速梯度,并通过两者的互相关处理来检测是否存在目标。振速梯度水听器是由对角布放的四个振速传感器组成,如图2所示。
本实施例中振速梯度水听器通过振速与振速梯度相关处理进行被动目标检测的方法,包括如下的步骤:
步骤一:二维振速梯度水听器八个通道的接收信号用vx1、vy1、vx2、vy2、vx3、vy3、vx4和vy4来表示,通过希尔伯特变换器将接收的实信号变换成复信号
步骤二:由四通道x轴复振速信号求振速梯度水听器中心点处x轴方向振速 v ~ o x = [ v ~ x 1 + v ~ x 2 + v ~ x 3 + v ~ x 4 ] / 4 ;
步骤三:由四通道y轴复振速信号求振速梯度水听器中心点处y轴方向振速 v ~ o y = [ v ~ y 1 + v ~ y 2 + v ~ y 3 + v ~ y 4 ] / 4 ;
步骤四:利用求解x轴横向振速梯度 ∂ v ~ x / ∂ x = j ( v ~ x 2 - v ~ x 1 ) / 2 , 利用求解y轴横向振速梯度 ∂ v ~ y / ∂ y = j ( v ~ y 4 - v ~ y 3 ) / 2
步骤五:利用求解纵向振速梯度 ∂ v ~ y / ∂ x = j ( v ~ y 2 - v ~ y 1 ) / 2 , 利用求解纵向振速梯度 ∂ v ~ x / ∂ y = j ( v ~ x 4 - v ~ x 3 ) / 2 ;
步骤六:利用引导方向α,对振速梯度水听器中心点处振速分量进行加权处理,得到合成振速 v ~ c ( α ) = v ~ o x c o s α + v ~ o y s i n α ;
步骤七:利用引导方向α,对振速梯度水听器中心点处横向振速梯度 和纵向振速梯度进行加权处理,获得合成振速梯度 d ~ c ( α ) = ∂ v ~ x / ∂ x · cos 2 α + ∂ v ~ y / ∂ y · sin 2 α + ( ∂ v ~ x / ∂ y + ∂ v ~ y / ∂ x ) · cos α sin α ;
步骤八:对M个采样点的合成振速与合成振速梯度进行相关处理,得到被动目标的检测统计量 D ( α , M ) = Σ n = 1 M | | v ~ c ( α , n ) · d ~ c ( α , n ) | | ;
步骤九:根据系统限定最大虚警概率Pf和上述检测统计量D(α,M),在噪声环境中统计计算出不同引导方向上的检测门限Dth(α,M);
步骤十:在引导方向α上,实时计算D(α,M),将D(α,M)与Dth(α,M)进行比较,如果D(α,M)>Dth(α,M),则目标存在,反之目标不存在;在不同的α上进行扫描,根据检测成功概率的大小可粗略估计目标存在的方位区域。
图3是实施例振速梯度水听器半径为0.01米,目标信号为中心频率为f0=1kHz的CW脉冲信号,采样频率为fs=8f0,噪声为数值模拟的各向同性噪声,系统地虚警概率Pf=0.05,目标在90°方位时,单个振速梯度水听器振速与振速梯度互相关检测器检测成功概率统计结果;图4是相同条件下,目标从45°方向入射时,单个振速梯度水听器振速与振速梯度互相关检测器检测成功概率统计结果。为了比较,同时给出了声压能量检测器和矢量水听器声能流检测器的成功检测概率随信噪比的变化曲线。从图3和图4可以看出,在给定的仿真条件下,本发明方法(振速与振速梯度相关检测)的可检测信噪比低于矢量水听器的声强检测法和声压水听器的能量检测法。这说明振速梯度水听器的振速与振速梯度相关检测法对被动目标探测距离更远,在军用、民用水下探测设备中有广泛的应用价值。
当然,本发明还可以有其他的实施例,熟悉本领域的技术人员可根据本发明做相应的改变,但这些改变都应该属于本发明所附权利要求的保护范围。

Claims (1)

1.一种二维振速梯度水听器的被动目标相关检测方法,包括如下步骤:
步骤一:将二维振速梯度水听器八个通道的接收数据vxi、vyi(i=1,2,3,4),通过希尔伯特变换器转换成复数信号(i=1,2,3,4);
步骤二:对四个通道的x轴复振速信号(i=1,2,3,4)求平均,求取振速梯度水听器中心点处x轴方向振速
步骤三:对四个通道的y轴复振速信号(i=1,2,3,4)求平均,求取振速梯度水听器中心点处y轴方向振速
步骤四:利用八通道振速复信号(i=1,2,3,4)求振速梯度水听器中心点处横向振速梯度和纵向振速梯度
步骤五:利用引导方向α,对振速梯度水听器中心点处振速分量进行加权处理,得到合成振速 v ~ c ( α ) = v ~ o x c o s α + v ~ o y s i n α ;
步骤六:利用引导方向α,对振速梯度水听器中心点处横向振速梯度 和纵向振速梯度进行加权处理,得到合成振速梯度 d ~ c ( α ) = ∂ v ~ x / ∂ xcos 2 α + ∂ v ~ y / ∂ y sin 2 α + ( ∂ v ~ x / ∂ y + ∂ v ~ y / ∂ x ) c o s α s i n α ;
步骤七:对M个采样点的合成振速与合成振速梯度进行相关处理,得到被动目标的检测统计量 D ( α , M ) = Σ n = 1 M | | v ~ c ( α , n ) · d ~ c ( α , n ) | | ;
步骤八:根据系统要求的最大虚警概率Pf和上述检验统计量D(α,M),在噪声环境中计算出不同引导方位上的检测门限Dth(α,M);
步骤九:实时计算D(α,M),将D(α,M)与Dth(α,M)进行比较,如果D(α,M)>Dth(α,M),则目标存在,反之目标不存在,并在不同的α上进行扫描,根据成功检测概率的分布粗略估计目标的方位。
CN201510606683.XA 2015-09-22 2015-09-22 一种二维振速梯度水听器的被动目标相关检测方法 Active CN106066472B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510606683.XA CN106066472B (zh) 2015-09-22 2015-09-22 一种二维振速梯度水听器的被动目标相关检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510606683.XA CN106066472B (zh) 2015-09-22 2015-09-22 一种二维振速梯度水听器的被动目标相关检测方法

Publications (2)

Publication Number Publication Date
CN106066472A true CN106066472A (zh) 2016-11-02
CN106066472B CN106066472B (zh) 2018-06-19

Family

ID=57418974

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510606683.XA Active CN106066472B (zh) 2015-09-22 2015-09-22 一种二维振速梯度水听器的被动目标相关检测方法

Country Status (1)

Country Link
CN (1) CN106066472B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108088547A (zh) * 2017-12-20 2018-05-29 青岛理工大学 一种基于小孔径二维矢量水听器阵的微弱目标被动检测方法
CN108195459A (zh) * 2017-12-20 2018-06-22 青岛理工大学 一种小孔径声压水听器阵的被动目标空间谱检测方法
CN112816940A (zh) * 2020-12-23 2021-05-18 中国船舶重工集团有限公司第七一0研究所 一种基于声压及质点振速的目标距离估计方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3976969A (en) * 1974-03-21 1976-08-24 Fried. Krupp Gesellschaft Mit Beschrankter Haftung Underwater sound wave direction finding system
CN102043145A (zh) * 2010-11-03 2011-05-04 中国科学院声学研究所 基于声矢量传感器均匀直线阵的快速宽带频域波束形成方法
CN102981146A (zh) * 2012-11-19 2013-03-20 哈尔滨工程大学 一种单矢量水听器被动定位方法
CN103267966A (zh) * 2013-05-20 2013-08-28 西北工业大学 一种单个二维压差式矢量水听器方位估计的方法
CN103454616A (zh) * 2013-08-27 2013-12-18 西北工业大学 一种十字型振速梯度水听器的方位估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3976969A (en) * 1974-03-21 1976-08-24 Fried. Krupp Gesellschaft Mit Beschrankter Haftung Underwater sound wave direction finding system
CN102043145A (zh) * 2010-11-03 2011-05-04 中国科学院声学研究所 基于声矢量传感器均匀直线阵的快速宽带频域波束形成方法
CN102981146A (zh) * 2012-11-19 2013-03-20 哈尔滨工程大学 一种单矢量水听器被动定位方法
CN103267966A (zh) * 2013-05-20 2013-08-28 西北工业大学 一种单个二维压差式矢量水听器方位估计的方法
CN103454616A (zh) * 2013-08-27 2013-12-18 西北工业大学 一种十字型振速梯度水听器的方位估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王绪虎: "单个二阶组合水听器的性能分析", 《振动与冲击》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108088547A (zh) * 2017-12-20 2018-05-29 青岛理工大学 一种基于小孔径二维矢量水听器阵的微弱目标被动检测方法
CN108195459A (zh) * 2017-12-20 2018-06-22 青岛理工大学 一种小孔径声压水听器阵的被动目标空间谱检测方法
CN112816940A (zh) * 2020-12-23 2021-05-18 中国船舶重工集团有限公司第七一0研究所 一种基于声压及质点振速的目标距离估计方法及装置
CN112816940B (zh) * 2020-12-23 2023-06-06 中国船舶重工集团有限公司第七一0研究所 一种基于声压及质点振速的目标距离估计方法及装置

Also Published As

Publication number Publication date
CN106066472B (zh) 2018-06-19

Similar Documents

Publication Publication Date Title
CN102183435B (zh) 一种基于多路径反射理论的海底密度和声速测量方法
CN103076594B (zh) 一种基于互相关的水声脉冲信号双阵元定位的方法
CN105629220B (zh) 一种基于单水听器的深海水声被动测距方法
Najeem et al. Open lake experiment for direction of arrival estimation using acoustic vector sensor array
Hui et al. Underwater acoustic channel
CN104678384B (zh) 一种波束域的声压差互相关谱分析水下目标速度估计方法
CN103076604B (zh) 一种基于频散特征的低频水声脉冲信号距离的测量方法
CN103969639B (zh) 五波束探鱼仪的信号处理系统及其信号处理方法
CN101915922A (zh) 拖曳线列阵被动测距方法
Morgado et al. Design and experimental evaluation of an integrated USBL/INS system for AUVs
CN103076590A (zh) 一种基于频率预估的水声脉冲信号的定位方法
CN109444864A (zh) 一种深海微弱多目标深度长时累积估计方法
Liang et al. The generalized cross-correlation method for time delay estimation of infrasound signal
CN106066472B (zh) 一种二维振速梯度水听器的被动目标相关检测方法
CN108612519A (zh) 油气井出砂的监测方法及装置
Zhang et al. Real-time localization for underwater equipment using an extremely low frequency electric field
CN109632258A (zh) 一种基于矢量传感器的收发分离的海洋内波声学检测方法
Zou et al. Vector hydrophone array development and its associated DOA estimation algorithms
Stafsudd et al. Localization of short-range acoustic and seismic wideband sources: Algorithms and experiments
RU2572085C1 (ru) Способ поиска, обнаружения и слежения за местоположением подвижных подводных объектов технической природы в акваториях морей и океанов
Gebbie et al. Passive acoustic array harbor security applications
RU2300781C1 (ru) Устройство гидрометеорологоакустических наблюдений за акваторией морского полигона
Alam et al. Optimal maneuvering of seismic sensors for localization of subsurface targets
Welton Cross correlation of omnidirectional, broadband signals scattered by a random pressure-release surface
Li et al. Acoustics vector sensor linear array passive ranging based on waveguide invariant

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210825

Address after: 266000 502, Hisense entrepreneurship center, No. 17, Shandong Road, Shinan District, Qingdao, Shandong Province

Patentee after: Shandong Beiming Technology Co.,Ltd.

Address before: 266520 No. 777 Jialingjiang Road, Huangdao District, Qingdao City, Shandong Province

Patentee before: Qindao University of Technology