CN105158734A - 一种基于阵不变量的单矢量水听器被动定位方法 - Google Patents
一种基于阵不变量的单矢量水听器被动定位方法 Download PDFInfo
- Publication number
- CN105158734A CN105158734A CN201510400349.9A CN201510400349A CN105158734A CN 105158734 A CN105158734 A CN 105158734A CN 201510400349 A CN201510400349 A CN 201510400349A CN 105158734 A CN105158734 A CN 105158734A
- Authority
- CN
- China
- Prior art keywords
- mrow
- tau
- time
- mtd
- frequency distribution
- 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
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000012417 linear regression Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 4
- 238000011160 research Methods 0.000 description 7
- 230000001052 transient effect Effects 0.000 description 4
- 230000004807 localization Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 101100206075 Escherichia coli (strain K12) tauA gene Proteins 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-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/22—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
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
本发明公开了一种基于阵不变量的单矢量水听器被动定位方法,包括以下步骤,对声压信号p(t)和水平振速信号vr(t)进行短时傅立叶变换,得到声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f);利用步骤一中的声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f)得到简正波俯仰角与简正波到达时间τ的模糊度平面根据步骤二中的简正波俯仰角与简正波到达时间τ的模糊度平面利用阵不变量的方法实现目标定位。本发明降低了定位方法对海洋环境先验知识的依赖性。
Description
技术领域
本发明属于水声定位领域,尤其涉及一种基于阵不变量的单矢量水听器被动定位方法。
背景技术
随着水声定位方法的发展,被动定位技术由于其良好的隐蔽性,在军事领域中有着良好的应用前景,多年来一直备受世界各国水声学者的关注和重视。一方面,随着声呐系统工作频段不断向低频发展,常规的声压水听器构成的水听器阵列在保持一定增益和束宽的条件下,其阵列孔径越来越大,这在一定程度上限制了其在实际工程中的应用。另一方面,当前成熟的被动定位技术对于较为复杂的浅海环境,还面临着诸多问题,如匹配场定位技术的环境失配问题等。矢量水听器迅速发展为被动定位开辟了新的视野,矢量水听器作为一种新型的水声传感器,在低频和甚低频段依然可以获得良好的空间增益,应用在低频和甚低频领域可以有效解决设备庞大的问题;矢量水听器可以同步共点测量声场的声压和振速信息,为水声信号的处理提供了全面的声场信息。
文献“Thearrayinvariant”(J.Acoust.Soc.Am.Vol.119,336-351,2006)在水平分层海洋波导环境下提出了一种基于浅海基阵不变量理论的瞬态声源定位方法,该方法利用了基阵常规波束形成后的被动波束-时间强度数据偏移不变性原理。
文献“Anewinvariantmethodforinstantaneoussourcerangeestimationinanoceanwaveguidefrompassivebeam-timeintensitydata”(J.Acoust.Soc.Am.116(4),2646,2004)推导出了一种基于阵不变量利用常规基阵估计瞬态源距离的方法。
文献“Rangeestimationofbroadbandnoisesourcesinanoceanwaveguideusingthearrayinvariant”(J.Acoust.Soc.Am.117(4),2577,2005)提出了在环境的先验知识未知情况下用于宽带噪声源距离估计的方法,这种方法适用于在水平分层海洋波导定位多个不相关噪声源。
文献“浅海波导中宽带声源被动测距研究”(中国海洋大学硕士论文,2011)对阵不变量测距进行了推广发展,提出了广义阵不变量,并对水平阵和垂直阵的阵不变量测距进行了试验验证。
文献“基于基阵不变量理论的水平阵被动估距方法研究”(声学与电子工程,104(4):14-17,2011)针对浅海水平阵测距问题,研究了一种基于浅海基阵不变量理论的瞬态声源定位方法。
以上文献的共同特点是基于阵不变量利用常规声压阵对目标进行定位,该方法对海洋环境先验知识依赖性较弱。
文献“基于单矢量水听器的水面目标运动分析”(声学技术,29(4):361-364,2010)介绍了一种基于单个三维矢量水听器对水面目标进行跟踪定位和运动速度估计的方法。
文献“基于干涉谱分析的单水听器被动定位技术研究”(哈尔滨工程大学硕士学位论文,2006)和文献“浅海声场干涉结构与宽带声源测距研究”(中国海洋大学硕士学位论文,2010)通过对声场干涉特性的分析,利用单矢量水听器对目标进行定位。
文献“时间反转镜被动定位技术研究”(哈尔滨工程大学博士论文,2008)和文献“单水听器被动时反定位研究”(声学与电子工程,84(4):15-17,2006)利用基于时间反转镜的单矢量水听器进行定位。
以上文献的共同特点是定位方法需要事先确知海洋环境信息。
发明内容
本发明的目的是为了解决现有方法对海洋环境先验知识依赖性弱以及定位方法需要事先确知海洋环境信息的问题,而提出了一种基于阵不变量的单矢量水听器被动定位方法。
一种基于阵不变量的单矢量水听器被动定位方法,包括以下步骤,
步骤一:对声压信号p(t)和水平振速信号vr(t)进行短时傅立叶变换,得到声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f);
步骤二:利用步骤一中的声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f)得到简正波俯仰角与简正波到达时间τ的模糊度平面
步骤三:根据步骤二中的简正波俯仰角与简正波到达时间τ的模糊度平面利用阵不变量的方法实现目标定位。
一种基于阵不变量的单矢量水听器被动定位方法,还可以包括:
1、声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f)分别为:
其中,数据窗函数h(t)定义为:
其中,T为窗长,e为自然常数,j为虚数单位,f为频率,t为时间,τ为时移。
2、简正波俯仰角与简正波到达时间τ的模糊度平面为:
其中,Rr(τ,f)为计算模糊度平面的统计量,为水平振速的时频分布Vr(τ,f)的平方,P2(τ,f)为声压的时频分布P(τ,f)的平方。
3、利用阵不变量的方法实现目标定位的方法为:
对简正波俯仰角与简正波到达时间τ的模糊度平面进行峰值检测得到线性回归数据点
MR为模糊度平面最大值,D为设定阈值,
目标定位结果为:
其中,r为目标与接收间的相对距离,为峰值检测到的线性回归数据点的俯仰角余割值,c为平均声速。
有益效果:
本发明的一种基于阵不变量的单矢量水听器被动定位方法,结合阵不变量理论和矢量水听器的特点,利用矢量水听器代替声压阵,降低阵列孔径;利用阵不变量理论,降低定位方法对海洋环境先验知识的依赖性。需要目标信号形式为宽带源。解决了现有方法对海洋环境先验知识依赖性弱以及定位方法需要事先确知海洋环境信息的问题。
附图说明
图1为本发明流程图;
图2为具体实施方式三中俯仰角与到达时间模糊度平面;
图3为具体实施方式四中峰值检测后的线性回归数据点;
图4为具体实施方式四中τ和之间线性关系斜率(阵不变量)估计。
具体实施方式
下面将结合附图对本发明做进一步详细说明。
具体实施方式一:结合图1说明本实施方式,一种基于阵不变量的单矢量水听器被动定位方法具体是按照以下步骤进行的:
步骤一、对矢量水听器接收的声压信号p(t)和水平振速信号vr(t)进行短时傅立叶变换,得到声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f);
步骤二、利用步骤一中的声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f)得到简正波俯仰角与简正波到达时间τ的模糊度平面
步骤三、根据步骤二中的简正波俯仰角与简正波到达时间τ的模糊度平面利用阵不变量的方法实现目标定位。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中对矢量水听器接收的声压信号p(t)和水平振速信号vr(t)进行短时傅立叶变换,得到声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f);具体过程为:
式中,数据窗函数h(t)定义为:
式中,T为窗长,e为自然常数,j为虚数单位,f为频率,t为时间,τ为时移。.
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中利用步骤一中的声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f)得到简正波俯仰角与简正波到达时间τ的模糊度平面具体过程为:
根据简正波理论,浅海波导中矢量声场中声压P(r,z,f)和水平振速Vr(r,z,f)的表述形式如下:
式中,为第n号简正波的俯仰角,r为目标与接收点间的水平相对距离,z为接收点深度,zs为目标深度,Zn为简正模态深度分布函数,μn为水平波数,e为自然常数,j为虚数单位,f为频率,n为简正波计数,为正整数;
其中,
联立(1)、(2)、(4)和(5)得到:
记俯仰角统计间隔为度,则空间角度允许间隔总数为:
式中,N为正整数;
由式(7)可得到空间角度允许间隔分别为: k为正整数,0<k<N;
若时间俯仰角落入各统计角度间隔的频点数量分别为m1、m2、m3、……mk……、mN,则可得到简正波俯仰角与简正波到达时间τ的模糊度平面如图2所示,
其中,Rr(τ,f)满足如下关系式:
式中,Rr(τ,f)为计算模糊度平面的统计量,为水平振速的时频分布Vr(τ,f)的平方,P2(τ,f)为声压的时频分布P(τ,f)的平方。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一、二或三不同的是:所述步骤三中根据步骤二中的模糊度平面利用阵不变量的方法实现目标定位;具体过程为:
对简正波俯仰角与简正波到达时间τ的模糊度平面进行峰值检测得到线性回归数据点如图3所示,并结合阵不变量的方法进行目标定位:
记MR为模糊度平面最大值,峰值检测按照式(10)进行:
式中,为线性回归数据点,D为设定阈值;
阵不变量的方法为:目标信号在平均声速c条件下从目标到接收器的绝对传播时间通过简正波到达时间τ和简正波俯仰角余割值的一阶导数进行计算,即:
式中,为俯仰角余割值,r为目标与接收间的相对距离;
根据式(11),将峰值检测到的线性回归数据点的俯仰角对应坐标转化为余割函数然后按照最小二乘原则估计τ和间线性关系的斜率为:如图4所示,从而由式(11)可实现目标定位,目标定位方程为:
其它步骤及参数与具体实施方式一、二或三相同。
一种基于阵不变量的单矢量水听器被动定位方法,本发明涉及基于阵不变量的单矢量水听器被动定位方法。本发明的目的是为了解决现有方法对海洋环境先验知识依赖性弱以及定位方法需要事先确知海洋环境信息的问题。通过以下技术方案实现的:步骤一、对声压信号p(t)和水平振速信号vr(t)进行短时傅立叶变换,得到声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f);步骤二、利用步骤一中的声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f)得到简正波俯仰角与简正波到达时间τ的模糊度平面步骤三、根据步骤二中的简正波俯仰角与简正波到达时间τ的模糊度平面利用阵不变量的方法实现目标定位。本发明属于水声定位领域。
Claims (4)
1.一种基于阵不变量的单矢量水听器被动定位方法,其特征在于:包括以下步骤,
步骤一:对声压信号p(t)和水平振速信号vr(t)进行短时傅立叶变换,得到声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f);
步骤二:利用步骤一中的声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f)得到简正波俯仰角与简正波到达时间τ的模糊度平面
步骤三:根据步骤二中的简正波俯仰角与简正波到达时间τ的模糊度平面利用阵不变量的方法实现目标定位。
2.根据权利要求1所述的一种基于阵不变量的单矢量水听器被动定位方法,其特征在于:
所述的声压的时频分布P(τ,f)和水平振速的时频分布Vr(τ,f)分别为:
其中,数据窗函数h(t)定义为:
其中,T为窗长,e为自然常数,j为虚数单位,f为频率,t为时间,τ为时移。
3.根据权利要求1所述的一种基于阵不变量的单矢量水听器被动定位方法,其特征在于:
所述的简正波俯仰角与简正波到达时间τ的模糊度平面为:
其中,Rr(τ,f)为计算模糊度平面的统计量,为水平振速的时频分布Vr(τ,f)的平方,P2(τ,f)为声压的时频分布P(τ,f)的平方。
4.根据权利要求1所述的一种基于阵不变量的单矢量水听器被动定位方法,其特征在于:
所述的利用阵不变量的方法实现目标定位的方法为:
对简正波俯仰角与简正波到达时间τ的模糊度平面进行峰值检测得到线性回归数据点
MR为模糊度平面最大值,D为设定阈值,
目标定位结果为:
其中,r为目标与接收间的相对距离,为峰值检测到的线性回归数据点的俯仰角余割值,c为平均声速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510400349.9A CN105158734B (zh) | 2015-07-09 | 2015-07-09 | 一种基于阵不变量的单矢量水听器被动定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510400349.9A CN105158734B (zh) | 2015-07-09 | 2015-07-09 | 一种基于阵不变量的单矢量水听器被动定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105158734A true CN105158734A (zh) | 2015-12-16 |
CN105158734B CN105158734B (zh) | 2018-04-24 |
Family
ID=54799652
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510400349.9A Active CN105158734B (zh) | 2015-07-09 | 2015-07-09 | 一种基于阵不变量的单矢量水听器被动定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105158734B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107783135A (zh) * | 2016-08-25 | 2018-03-09 | 中国科学院声学研究所 | 一种三元矢量阵被动测距方法 |
CN109632258A (zh) * | 2019-01-30 | 2019-04-16 | 杭州电子科技大学 | 一种基于矢量传感器的收发分离的海洋内波声学检测方法 |
CN110824429A (zh) * | 2019-10-28 | 2020-02-21 | 西北工业大学 | 深海环境下利用非同步垂直阵的宽带声源被动定位方法 |
CN111273265A (zh) * | 2020-03-02 | 2020-06-12 | 中国人民解放军国防科技大学 | 基于demon谱的单矢量水听器自主测距装置及方法 |
CN112269163A (zh) * | 2020-09-30 | 2021-01-26 | 黑龙江工程学院 | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 |
CN112684437A (zh) * | 2021-01-12 | 2021-04-20 | 浙江大学 | 一种基于时域warping变换的被动测距方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101271155A (zh) * | 2008-04-25 | 2008-09-24 | 哈尔滨工程大学 | 用浅海低频相干矢量声场干涉结构探测并判别目标的方法 |
CN102981146A (zh) * | 2012-11-19 | 2013-03-20 | 哈尔滨工程大学 | 一种单矢量水听器被动定位方法 |
CN104091048A (zh) * | 2014-06-19 | 2014-10-08 | 哈尔滨工程大学 | 被动垂直双矢量水听器的水下目标深度测定方法 |
-
2015
- 2015-07-09 CN CN201510400349.9A patent/CN105158734B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101271155A (zh) * | 2008-04-25 | 2008-09-24 | 哈尔滨工程大学 | 用浅海低频相干矢量声场干涉结构探测并判别目标的方法 |
CN102981146A (zh) * | 2012-11-19 | 2013-03-20 | 哈尔滨工程大学 | 一种单矢量水听器被动定位方法 |
CN104091048A (zh) * | 2014-06-19 | 2014-10-08 | 哈尔滨工程大学 | 被动垂直双矢量水听器的水下目标深度测定方法 |
Non-Patent Citations (3)
Title |
---|
姚玉玲 等: "水声信号时频分析方法比较及应用研究", 《中国海洋大学学报》 * |
师俊杰 等: "基于阵不变量的垂直阵被动定位方法研究", 《声学技术》 * |
张宇 等: "应用匹配场实现单矢量水听器被动定位", 《传感器与微系统》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107783135A (zh) * | 2016-08-25 | 2018-03-09 | 中国科学院声学研究所 | 一种三元矢量阵被动测距方法 |
CN107783135B (zh) * | 2016-08-25 | 2019-10-22 | 中国科学院声学研究所 | 一种三元矢量阵被动测距方法 |
CN109632258A (zh) * | 2019-01-30 | 2019-04-16 | 杭州电子科技大学 | 一种基于矢量传感器的收发分离的海洋内波声学检测方法 |
CN110824429A (zh) * | 2019-10-28 | 2020-02-21 | 西北工业大学 | 深海环境下利用非同步垂直阵的宽带声源被动定位方法 |
CN110824429B (zh) * | 2019-10-28 | 2022-09-13 | 西北工业大学 | 深海环境下利用非同步垂直阵的宽带声源被动定位方法 |
CN111273265A (zh) * | 2020-03-02 | 2020-06-12 | 中国人民解放军国防科技大学 | 基于demon谱的单矢量水听器自主测距装置及方法 |
CN111273265B (zh) * | 2020-03-02 | 2020-09-29 | 中国人民解放军国防科技大学 | 基于demon谱的单矢量水听器自主测距装置及方法 |
CN112269163A (zh) * | 2020-09-30 | 2021-01-26 | 黑龙江工程学院 | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 |
CN112269163B (zh) * | 2020-09-30 | 2023-04-25 | 黑龙江工程学院 | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 |
CN112684437A (zh) * | 2021-01-12 | 2021-04-20 | 浙江大学 | 一种基于时域warping变换的被动测距方法 |
CN112684437B (zh) * | 2021-01-12 | 2023-08-11 | 浙江大学 | 一种基于时域warping变换的被动测距方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105158734B (zh) | 2018-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105158734B (zh) | 一种基于阵不变量的单矢量水听器被动定位方法 | |
CN103076594B (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN101957443B (zh) | 声源定位方法 | |
CN102183435A (zh) | 一种基于多路径反射理论的海底密度和声速测量方法 | |
CN101915922A (zh) | 拖曳线列阵被动测距方法 | |
CN104765037A (zh) | 基于短垂直阵的水下目标定位稳健方法 | |
CN103713276B (zh) | 基于最小互熵谱分析的波达方向估计方法 | |
CN103076590A (zh) | 一种基于频率预估的水声脉冲信号的定位方法 | |
CN112285647B (zh) | 一种基于稀疏表示与重构的信号方位高分辨估计方法 | |
CN113011006B (zh) | 一种基于互相关函数脉冲波形匹配的目标深度估计方法 | |
CN112098938B (zh) | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 | |
CN116609725B (zh) | 一种利用深海垂直阵的窄带线谱目标深度估计方法及系统 | |
CN102333052A (zh) | 一种适用于浅海低频条件的水声信号盲解卷方法 | |
CN109597021A (zh) | 一种波达方向估计方法及装置 | |
Pham et al. | Real-time implementation of MUSIC for wideband acoustic detection and tracking | |
CN113866718B (zh) | 一种基于互质阵的匹配场被动定位方法 | |
CN115236592B (zh) | 一种基于单阵元时频曲线匹配的冰声定位方法 | |
CN110780340A (zh) | 一种浅海弹性软沉积层横波声速或厚度的反演方法 | |
CN106066472A (zh) | 一种二维振速梯度水听器的被动目标相关检测方法 | |
CN103513249A (zh) | 一种宽带相干模基信号处理方法及系统 | |
CN104215964A (zh) | 一种多列等差频率原波相互作用形成参量阵的声场获取方法 | |
CN106595835B (zh) | 一种基于错位叠加算法的声信标信号提取方法及系统 | |
Zamanizadeh et al. | Source localization from time-differences of arrival using high-frequency communication signals | |
Soares et al. | Environmental inversion using high-resolution matched-field processing | |
Kumar et al. | Doa estimation using compressive beamforming in shallow ocean using acoustic vector sensors |
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 |