CN116299156B - 一种水听器垂直阵阵元位置估计方法及其优化策略 - Google Patents
一种水听器垂直阵阵元位置估计方法及其优化策略 Download PDFInfo
- Publication number
- CN116299156B CN116299156B CN202310601868.6A CN202310601868A CN116299156B CN 116299156 B CN116299156 B CN 116299156B CN 202310601868 A CN202310601868 A CN 202310601868A CN 116299156 B CN116299156 B CN 116299156B
- Authority
- CN
- China
- Prior art keywords
- far
- hydrophone
- vertical array
- array
- array element
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005457 optimization Methods 0.000 title claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000005315 distribution function Methods 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 230000005236 sound signal Effects 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000009304 pastoral farming Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 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
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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
本申请属于水声信号处理技术领域,提供一种水听器垂直阵阵元位置估计方法及其优化策略,所述估计方法包括以下步骤:通过至少两个远场声源分别向水听器垂直阵发射标定信号,对于每个远场声源发射的标定信号,分别通过各个阵元进行接收并对接收信号进行波束形成处理,基于波束形成处理的结果估计水听器垂直阵的实际布设轨迹的俯仰角在该远场声源所在的远场平面上对应的投影俯仰角;基于至少两个所述投影俯仰角获取的估计值以及水听器垂直阵的实际布设轨迹的方位角的估计值;基于所述、估计各个阵元的实际位置。本申请利用简正波的方向特性,通过两个非共面的远场声源发射的标定信号进行波束形成,能够同时获得水听器垂直阵的完整阵型信息。
Description
技术领域
本申请属于水声信号处理技术领域,具体提供一种水听器垂直阵阵元位置估计方法及其优化策略。
背景技术
随着声纳设备的不断发展和波束形成技术的成熟,多基元的水听器阵列被越来越多地应用于目标探测和定位的领域中,水听器垂直阵列作为水中传感器的一种常用分布方式,在其中扮演了重要的角色。对于水听器垂直阵,其各个阵元位置的准确获取对与垂直阵增益处理、探测定位性能有着显著的影响,尤其对于高频信号,对阵型和阵元位置的估计即使存在较小的偏差也会使得阵列信号处理算法的性能大幅度下降,特别是当垂直阵的阵列尺度较大时,在海流及重力作用下会产生明显的倾斜、悬垂等现象,因此,在使用水听器垂直阵进行水声信号处理前,必须对阵元实际位置进行准确地估计以用于信号校准。
目前关于垂直阵阵型的估计方法相关研究较少,主要分为两类,一类是利用垂直阵布设压力传感器,通过测量垂直阵首尾和中间位置所处深度来估计垂直阵的倾斜程度,但该方法只能得到垂直阵的俯仰角,不能估计垂直阵的倾斜方位角;另一类是在垂直阵正上方布设声源发射信号或者利用环境噪声互相关处理等声学的手段估计各阵元间的相对到达时间,结合水体内的声速分布来估计垂直阵的倾斜程度,该方法同样也只能估计垂直阵的俯仰角,如要准确估计垂直阵在三维坐标中的阵型,仍需要结合其他测量手段。
为此,亟需一种算法简单,计算量小,并能够同时估计水听器垂直阵俯仰角与方位角的方法。本专利所提的算法基于简正波的特性,利用两个正交平面或者已知夹角的两个非正交平面内的远场声源发射标定信号,并通过水听器垂直阵接收信号做常规波束形成处理,就可以估计垂直阵的倾斜俯仰角和方位角,算法简单,计算量小,对比其他方法有着较为明显的优势。
发明内容
本申请的目的在于解决上述现有技术中存在的问题,提供一种水听器垂直阵阵元位置估计方法,以及对该估计方法进一步进行优化的策略,该估计方法及其优化策略利用简正波的方向特性,通过两个正交平面或者已知夹角的两个非正交平面内的远场声源辐射信号的波束形成结果同时获得水听器垂直阵的完整阵型信息。
本申请的第一方面提供一种水听器垂直阵阵元位置估计方法,用于对构成水听器垂直阵的各个阵元的真实位置进行估计,包括以下步骤:
S1,通过至少两个远场声源分别向水听器垂直阵发射标定信号,其中,任意两个远场声源所在位置与水听器垂直阵的理想布设轨迹所成的平面互相之间不共面;
S2,对于每个远场声源发射的标定信号,分别通过各个阵元进行接收并对接收信号进行波束形成处理,基于波束形成处理的结果估计水听器垂直阵的实际布设轨迹的俯仰角在该远场声源所在的远场平面上对应的投影俯仰角,其中,每个远场声源所在的远场平面为该远场声源与水听器垂直阵的理想布设轨迹构成的平面;
S3,基于至少两个所述投影俯仰角获取的估计值/>以及水听器垂直阵的实际布设轨迹的方位角/>的估计值/>;
S4,基于所述、/>估计各个阵元的实际位置。
优选地,所述标定信号为宽带短脉冲声信号;以及,所述远场声源与水听器垂直阵的距离大于等于水听器垂直阵所在位置的水深的10倍。
优选地,至少两个远场声源对应的远场平面互相垂直。
优选地,所述水听器垂直阵的阵元数量大于等于8,且各个阵元等间距设置。
进一步地,在步骤S2中对于每个远场声源发射的标定信号,基于以下步骤进行所述波束形成处理:
S21,通过各个阵元接收所述标定信号并进行下式的傅里叶变换:
,
其中,为水听器垂直阵的阵元数量,/>为阵元的序号,/>为第/>个阵元的深度,/>为积分时间,/>为所述标定信号传播至第/>个阵元所产生的时域接收信号,为/>的频谱;
S22,将变换为简正波表示:
,
其中为水体密度,/>为标定信号的频谱,/>为远场声源的深度,/>为远场声源和水听器垂直阵的水平间距,/>为简正波的序号,/>为有效简正波的最大序号,/>、/>分别为第/>号简正波的本征波数和本征函数;
S23,进一步将进行变换为下式:
,
其中,为第/>号简正波模态射线的掠射角,/>,/>为水体平均声速;
S24,基于下式对各个阵元的进行波束形成,得到接收信号的频率-掠射角二维分布函数:
,
其中,为接收信号关于频率/>与掠射角/>的二维分布函数,/>为阵元的间距;
S25,基于下式对所述进行逆傅里叶变换,得到接收信号的时间-掠射角二维分布函数:
,
其中,为接收信号关于接收时间/>与掠射角/>的二维分布函数,/>、/>分别为频率积分的上下限;
S26,基于下式估计所述在每个远场声源所在的远场平面上对应的投影俯仰角:
,
其中,为远场声源的序号,/>为所述/>在第/>个远场声源所在的远场平面上对应的投影俯仰角,/>分别为第/>号简正波的上行波的掠射角及下行波的掠射角。
进一步地,基于下式确定所述、/>:
,
其中,分别为任意两个不同的远场声源的序号。
本申请的第二方面还提供一种水听器垂直阵阵元位置估计方法的优化策略,用于对前述的水听器垂直阵阵元位置估计方法进行优化,包括以下步骤:
A1,使用前述的水听器垂直阵阵元位置估计方法估计各个阵元的实际位置;
A2,基于水听器垂直阵的前个阵元构造第一垂直阵子阵,以及根据水听器垂直阵的第/>个阵元至最后一个阵元构造第二垂直子阵,其中,第一垂直子阵与第二垂直子阵的阵元数量均大于等于8;
A3,分别对第一垂直阵子及第二垂直子阵执行步骤S1至S3,获取第一垂直阵子阵的实际布设轨迹的俯仰角、方位角/>的估计值/>、/>,以及获取第二垂直子阵的实际布设轨迹的俯仰角/>、方位角/>的估计值/>、/>;
A4,如果与/>的差值或/>与/>的差值大于预设的阈值,则执行步骤A5,否则改变/>的取值并返回执行步骤A1,直到遍历/>的取值范围;
A5,基于所述、/>、/>、/>重新估计各个阵元的实际位置。
优选地,的取值范围为水听器垂直阵阵元数量的1/4至3/4。
本申请的实施例提供的一种一种水听器垂直阵阵元位置估计方法及其优化策略,利用倾斜状态下的水听器垂直阵对多个非共面的远场声源的发射的标定信号进行波数形成,利用简正波上行波与下行波相对于入射方向的对称性方便地确定水听器垂直阵实际布设轨迹在各个远场平面上的投影信息,从而能够同时获取对水听器垂直阵布设轨迹的俯仰角及方位角的估计结果。
附图说明
图1为一个具体的布设于海面以下的水听器垂直阵的轨迹示意图;
图2为根据本申请实施例提供的水听器垂直阵阵元位置估计方法的流程图;
图3为根据本申请优选的实施例,在两个正交的平面作为设置远场声源的远场平面的示意图;
图4为根据本申请优选的实施例,以XZ平面为远场平面布置远场声源的示意图;
图5为水听器垂直阵实际布设轨迹在XZ平面上的投影相对于OZ轴具有逆时针方向的偏移时,远场声信号入射的示意图;
图6为水听器垂直阵实际布设轨迹在XZ平面上的投影相对于OZ轴具有顺时针方向的偏移时,远场声信号入射的示意图;
图7为在一些具体的实施例中,简正波分解的原理示意图;
图8为在一些具体的实施例中,水听器垂直阵在远场平面内的投影俯仰角为0°时,时间-掠射角二维分布函数的图像;
图9为在一些具体的实施例中,水听器垂直阵在远场平面内的投影俯仰角为5°时,时间-掠射角二维分布函数的图像;
图10为在一些实施例中,水听器垂直阵的阵元在海流冲击及线缆拖拽的共同作用下形成的三维轨迹的示意图。
具体实施方式
以下,基于优选的实施方式并参照附图对本申请进行进一步说明。
在本申请实施例中的描述中,需要说明的是,若出现术语“上”、“下”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,或者是本申请实施例的产品使用时惯常摆放的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。此外,本申请的描述中,为了区分不同的单元,本说明书上用了第一、第二等词汇,但这些不会受到制造的顺序限制,也不能理解为指示或暗示相对重要性,其在本申请的详细说明与权利要求书上,其名称可能会不同。此外,为了方便理解,放大或者缩小了图纸上的各种构件,但这种做法不是为了限制本申请的保护范围。
图1示出了一个具体的布设于海面以下的水听器垂直阵的轨迹示意图,在使用水听器垂直阵进行水声信号处理的过程中,为了便于对位置、方向进行表示,一般以位于最下端的阵元为原点O建立正交坐标系,以Z轴作为水听器垂直阵的理想布设轨迹。
在水听器垂直阵的实际布设及使用过程中,其实际布设轨迹一般受海流的影响产生倾斜,从而与理想布设轨迹并不重合。如图1所示,原点O与最上端的阵元A的连线为水听器垂直阵的实际布设轨迹,其与理想布设轨迹(即OZ轴)的姿态偏离可以通过俯仰角及方位角/>表示,其中,/>为Z轴正方向与OA连线之间的夹角,/>为OA连线在水平面(即XY平面)上的投影与X轴正方向之间的夹角。上述水听器垂直阵的实际布设轨迹与理想布设轨迹的偏离对后续垂直阵数据的处理造成了很大的影响,尤其是高频声信号处理对垂直阵的阵型更加敏感,因此有必要获取水听器垂直阵实际布设轨迹的精确数据以保证对后续的水声信号处理提供精确校正。
由于海洋信道中存在多途声传播特性,多途干涉使得海洋信道形成天然的角度滤波器,从而使得水听器垂直阵各阵元接收的声信号在不同到达时刻强弱不一致,声信号的强弱分布跟声源的深度和距离有关,这对利用各阵元脉冲相对到达时间估计垂直阵倾斜角度造成困难;此外,还有利用压力传感器海试测试声信号的方法对垂直阵布设轨迹进行估计,然而其只能估计垂直阵的俯仰角,而不能同时有效地估计垂直阵的方位角。
为此,本申请提供水听器垂直阵阵元位置估计方法,该方法能有同时获取对水听器垂直阵布设轨迹的俯仰角及方位角的估计结果,图2示出了在一些优选的实施例中,该估计方法的实施流程示意图,如图2所示,该估计方法包括以下步骤:
S1,通过至少两个远场声源分别向水听器垂直阵发射标定信号,其中,任意两个远场声源所在位置与水听器垂直阵的理想布设轨迹所成的平面互相之间不共面;
S2,对于每个远场声源发射的标定信号,分别通过各个阵元进行接收并对接收信号进行波束形成处理,基于波束形成处理的结果估计水听器垂直阵的实际布设轨迹的俯仰角在该远场声源所在的远场平面上对应的投影俯仰角,其中,每个远场声源所在的远场平面为该远场声源与水听器垂直阵的理想布设轨迹构成的平面;
S3,基于至少两个所述投影俯仰角获取的估计值/>以及水听器垂直阵的实际布设轨迹的方位角/>的估计值/>;
S4,基于所述、/>估计各个阵元的实际位置。
本申请提供的水听器垂直阵阵元位置估计方法,首先通过步骤S1在至少两个与垂直阵的理想布设轨迹不共面的远场平面发射标定信号,然后通过步骤S2分别对垂直阵的实际俯仰角在各个远场平面的投影进行估计,最后利用至少两个远场平面的估计结果,通过步骤S3及S4得到对垂直阵俯仰角和方位角的估计,以下结合附图及具体实施例,对步骤S1至步骤S4的具体实施方式进行详细说明。
图3示出了在一个优选的实施例中,分别以两个正交的平面,即XZ平面与YZ平面为作为设置远场声源的远场平面的示意图,水听器垂直阵的实际布设轨迹OA在XZ平面与YZ平面的投影分别与OZ成夹角与/>,在本申请中,/>与/>被称为OA的俯仰角/>在远场平面的投影俯仰角,当获取了/>、/>的估计值/>、/>后,利用几何关系即可得到俯仰角/>及方位角的估计值/>、/>,进而利用阵元之间的间距信息估计各个阵元的位置。
需要说明的是,图3所示的以两个正交平面作为远场平面的实施例,仅为本申请的一种优选的实施方式,在本申请的其他实施例中,远场平面的数量还可以大于2,且任意两个远场平面之间的角度可以成不等于90°的夹角,只要保证远场平面的数量不小于2,且任意两个远场平面相对于水听器垂直阵的理想布设轨迹不共面,均不脱离本申请的发明构思。
图4以XZ平面为例,示意性地说明在远场平面中布置远场声源的实施方式,本领域技术人员通过图4可以想象在其他由远场声源与水听器垂直阵的理想布设轨迹所形成的远场平面(如YZ平面或其他远场平面)内,远场声源的具体布置方式。
如图4所示,声源1与水听器垂直阵的距离远大于水听器垂直阵所在位置的水深,例如,在一些优选的实施例中,两者之比大于等于10,因此,可以将从该声源1作为XZ平面内的远场声源,从其发射的声信号可以视为从水平方向入射的远场声信号;此外,为保证估计精度,在本申请的一些优选的实施例中,水听器垂直阵的阵元数量大于等于8(例如,水听器垂直阵可以是8元阵、16元阵、32元阵等),且各个阵元等间距设置。
图5、图6进一步示出了水听器垂直阵实际布设轨迹在XZ平面上的投影相对于OZ轴具有逆时针方向及顺时针方向的偏移时,远场声信号入射的情景,如图5、图6所示,水听器垂直阵在XZ平面上的投影具有投影俯仰角,并且上述倾斜,使得垂直阵进行水声定位时,其用于标定水平方向的X’轴相对于真实的X轴产生了夹角/>,且该夹角/>=/>,且在此种情况下,在XZ平面上水平入射的远场声信号(即相对于X轴以0°掠射角入射)将被倾斜的水听器垂直阵“认为”是沿掠射角/>入射;由此可见,可以利用水平入射的远场声信号作为标定信号,倾斜的水听器垂直阵将对标定信号相对于X’轴所具有的掠射角/>进行估计,且该/>的估计值也就是/>的估计值/>。
通过步骤S1由多个远场声源发射标定信号后,步骤S2分别对各个远场声源发射的标定信号进行处理,得到水听器垂直阵在各个远场平面内对应的投影俯仰角的估计值,具体地,在本申请的实施例中,S2包括以下步骤:
S21,通过各个阵元接收所述标定信号并进行下式的傅里叶变换:
,
其中,为水听器垂直阵的阵元数量,/>为阵元的序号,/>为第/>个阵元的深度,为积分时间,/>为所述标定信号传播至第/>个阵元所产生的时域接收信号,为/>的频谱。
S22,将变换为简正波表示:
,
其中为水体密度,/>为标定信号的频谱,/>为远场声源的深度,/>为远场声源和水听器垂直阵的水平间距,/>为简正波的序号,/>为有效简正波的最大序号,/>、/>分别为第/>号简正波的本征波数和本征函数。
S23,进一步将进行变换为下式:
,
其中,为第/>号简正波模态射线的掠射角,/>,/>为水体平均声速;
S24,基于下式对各个阵元的进行波束形成,得到接收信号的频率-掠射角二维分布函数:
,
其中,为接收信号关于频率/>与掠射角/>的二维分布函数,/>为阵元的间距。
S25,基于下式对所述进行逆傅里叶变换,得到接收信号的时间-掠射角二维分布函数:
,
其中,为接收信号关于接收时间/>与掠射角/>的二维分布函数,/>、/>分别为频率积分的上下限;
步骤S22至S25通过简正波波束形成得到简正波上行波与下行波掠射角信息,图7以OX平面为例说明简正波分解的原理,如图7所示,对于垂直阵接收远场声信号的情况,可以将海洋环境简化为水平不变的分层介质波导,在此模型下将各个阵元接收到的频域信号分解为多个简正波,其模态射线可以由上行波与下行波合成得到。
基于上述原理,通过对垂直阵各个阵元的频域信号进行波束形成,则可以得到接收信号的频率-掠射角二维分布函数,其中,当阵元数量远大于1时,各号简正波上行波与下行波的掠射角相对于其自身模态射线的掠射角近似地成对称关系,因此,可以通过频率-掠射角二维分布函数获取各上行波与下行波的掠射角,进而得到各号简正波模态射线的掠射角信息。
对于宽带短脉冲信号,当信号频率远大于简正波截至频率时,简正波的本征函数、群速度和模态射线俯仰角随频率几乎不变,因此可以进一步对进行逆傅里叶变换处理,以得到接收信号的时间-掠射角二维分布函数/>,图8、图9分别示出了水听器垂直阵在某一远场平面内的投影俯仰角为0°和5°时,通过上述步骤得到的时间-掠射角二维分布函数/>的图像,图中声强较高的区域分别代表了多号简正波的上行波与下行波的到达时刻与掠射角。
通过图8、图9可以看出,由于各有效简正波群速度不同,使得各号简正波到达时间不同,各号简正波上行波和下行波不仅在掠射角上是可分离的,在相对到达时间上也是可分离的,从而提高了简正波上行波掠射角和下行波掠射角的估计精度。
进一步地,可以通过人工识别,或者现有的各种图像识别算法,从图8、图9中提取前个上行波与下行波的掠射角,并在步骤S26中基于下式估计所述/>在各个远场声源所在的远场平面上对应的投影俯仰角:
,
其中,为远场声源的序号,/>为所述/>在第/>个远场声源所在的远场平面上对应的投影俯仰角,/>分别为第/>号简正波的上行波的掠射角及下行波的掠射角。
进一步地,如在获取了至少两个远场平面的投影俯仰角的估计值后,通过步骤S3即可水听器垂直阵实际布设轨迹的俯仰角及方位角进行估计。
具体地,在一些优选的实施例中,对于两个互相垂直的远场平面,可以通过下式进行估计:
,
其中,分别为任意两个不同的远场声源的序号,以图3为例,图中水听器垂直阵的实际布设轨迹OA在XZ平面和YZ平面的投影俯仰角分别为/>、/>,其估计值分别为/>、,带入上式得到:
。
此外,对于两个夹角不为90°的远场平面,进行俯仰角及方位角估计时还需要考虑两个远场平面之间的夹角关系,上述处理均为本领域技术人员所知晓,在此不再赘述。
通过前述水听器垂直阵阵元位置优化方法虽然能够获得对整体阵型的俯仰角及方位角的估计,然而,在一些存在复杂海流的区域,如图10所示,由于海水流速随深度变化的程度较为剧烈,此时垂直阵的阵型很难保持一致,在部分流速较大的深度,各个阵元在海流冲击及线缆拖拽的共同作用下形成三维轨迹上的如图中A’点所示的“拐点”,该拐点将使得其上下两部分阵型轨迹呈现明显的差异,甚至无法使用单一的俯仰角及方位角对其进行表征,为此,本申请的一些优选的实施例还提供了一种对前述水听器垂直阵阵元位置估计方法的优化策略,包括以下步骤:
A1,使用前述的水听器垂直阵阵元位置估计方法估计各个阵元的实际位置;
A2,基于水听器垂直阵的前个阵元构造第一垂直阵子阵,以及根据水听器垂直阵的第/>个阵元至最后一个阵元构造第二垂直子阵,其中,第一垂直子阵与第二垂直子阵的阵元数量均大于等于8;
A3,分别对第一垂直阵子及第二垂直子阵执行步骤S1至S3,获取第一垂直阵子阵的实际布设轨迹的俯仰角、方位角/>的估计值/>、/>,以及获取第二垂直子阵的实际布设轨迹的俯仰角/>、方位角/>的估计值/>、/>;
A4,如果与/>的差值或/>与/>的差值大于预设的阈值,则执行步骤A5,否则改变/>的取值并返回执行步骤A1,直到遍历/>的取值范围;
A5,基于所述、/>、/>、/>重新估计各个阵元的实际位置。
上述各步骤中,首先获取水听器垂直阵的整体阵型信息,然后分别在一定遍历范围内将垂直阵分为两个子阵,分别获取其阵型信息并判断是否存在拐点,如果存在,则利用该拐点分别估计上下两个子阵的阵元位置,否则仍以整个垂直阵的阵型信息估计阵元位置。
通过对垂直阵在海流冲击下的阵型变化进行分析发现,上述拐点往往出现在整个垂直阵的中间区域,为此,在一些优选的实施例中,的取值范围为水听器垂直阵阵元数量的1/4至3/4。
以上对本申请的具体实施方式作了详细介绍,对于本技术领域的技术人员来说,在不脱离本申请原理的前提下,还可以对本申请进行若干改进和修饰,这些改进和修饰也属于本申请权利要求的保护范围。
Claims (8)
1.一种水听器垂直阵阵元位置估计方法,用于对构成水听器垂直阵的各个阵元的真实位置进行估计,其特征在于,包括以下步骤:
S1,通过至少两个远场声源分别向水听器垂直阵发射标定信号,其中,任意两个远场声源所在位置与水听器垂直阵的理想布设轨迹所成的平面互相之间不共面;
S2,对于每个远场声源发射的标定信号,分别通过各个阵元进行接收并对接收信号进行波束形成处理,基于波束形成处理的结果估计水听器垂直阵的实际布设轨迹的俯仰角在该远场声源所在的远场平面上对应的投影俯仰角,其中,每个远场声源所在的远场平面为该远场声源与水听器垂直阵的理想布设轨迹构成的平面;
S3,基于至少两个所述投影俯仰角获取的估计值/>以及水听器垂直阵的实际布设轨迹的方位角/>的估计值/>;
S4,基于所述、/>估计各个阵元的实际位置。
2.根据权利要求1所述的水听器垂直阵阵元位置估计方法,其特征在于:
所述标定信号为宽带短脉冲声信号;以及,
所述远场声源与水听器垂直阵的距离大于等于水听器垂直阵所在位置的水深的10倍。
3.根据权利要求1所述的水听器垂直阵阵元位置估计方法,其特征在于:
至少两个远场声源对应的远场平面互相垂直。
4.根据权利要求1所述的水听器垂直阵阵元位置估计方法,其特征在于:
所述水听器垂直阵的阵元数量大于等于8,且各个阵元等间距设置。
5.根据权利要求4所述的水听器垂直阵阵元位置估计方法,其特征在于,在步骤S2中对于每个远场声源发射的标定信号,基于以下步骤进行所述波束形成处理:
S21,通过各个阵元接收所述标定信号并进行下式的傅里叶变换:
,
其中,为水听器垂直阵的阵元数量,/>为阵元的序号,/>为第/>个阵元的深度,/>为积分时间,/>为所述标定信号传播至第/>个阵元所产生的时域接收信号,/>为的频谱;
S22,将变换为简正波表示:
,
其中为水体密度,/>为标定信号的频谱,/>为远场声源的深度,/>为远场声源和水听器垂直阵的水平间距,/>为简正波的序号,/>为有效简正波的最大序号,/>、/>分别为第/>号简正波的本征波数和本征函数;
S23,进一步将进行变换为下式:
,
其中,为第/>号简正波模态射线的掠射角,/>,/>为水体平均声速;
S24,基于下式对各个阵元的进行波束形成,得到接收信号的频率-掠射角二维分布函数:
,
其中,为接收信号关于频率/>与掠射角/>的二维分布函数,/>为阵元的间距;
S25,基于下式对所述进行逆傅里叶变换,得到接收信号的时间-掠射角二维分布函数:
,
其中,为接收信号关于接收时间/>与掠射角/>的二维分布函数,/>、/>分别为频率积分的上下限;
S26,基于下式估计所述在每个远场声源所在的远场平面上对应的投影俯仰角:
,
其中,为远场声源的序号,/>为所述/>在第/>个远场声源所在的远场平面上对应的投影俯仰角,/>分别为第/>号简正波的上行波的掠射角及下行波的掠射角。
6.根据权利要求5所述的水听器垂直阵阵元位置估计方法,其特征在于,基于下式确定所述、/>:
,
其中,分别为任意两个不同的远场声源的序号。
7.一种水听器垂直阵阵元位置估计方法的优化策略,用于对权利要求1所述的水听器垂直阵阵元位置估计方法进行优化,其特征在于,包括以下步骤:
A1,使用如权利要求1所述的水听器垂直阵阵元位置估计方法估计各个阵元的实际位置;
A2,基于水听器垂直阵的前个阵元构造第一垂直阵子阵,以及根据水听器垂直阵的第/>个阵元至最后一个阵元构造第二垂直子阵,其中,第一垂直子阵与第二垂直子阵的阵元数量均大于等于8;
A3,分别对第一垂直阵子及第二垂直子阵执行步骤S1至S3,获取第一垂直阵子阵的实际布设轨迹的俯仰角、方位角/>的估计值/>、/>,以及获取第二垂直子阵的实际布设轨迹的俯仰角/>、方位角/>的估计值/>、/>;
A4,如果与/>的差值或/>与/>的差值大于预设的阈值,则执行步骤A5,否则改变/>的取值并返回执行步骤A1,直到遍历/>的取值范围;
A5,基于所述、/>、/>、/>重新估计各个阵元的实际位置。
8.根据权利要求7所述的水听器垂直阵阵元位置估计方法的优化策略,其特征在于:
的取值范围为水听器垂直阵阵元数量的1/4至3/4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310601868.6A CN116299156B (zh) | 2023-05-26 | 2023-05-26 | 一种水听器垂直阵阵元位置估计方法及其优化策略 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310601868.6A CN116299156B (zh) | 2023-05-26 | 2023-05-26 | 一种水听器垂直阵阵元位置估计方法及其优化策略 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116299156A CN116299156A (zh) | 2023-06-23 |
CN116299156B true CN116299156B (zh) | 2023-08-15 |
Family
ID=86824342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310601868.6A Active CN116299156B (zh) | 2023-05-26 | 2023-05-26 | 一种水听器垂直阵阵元位置估计方法及其优化策略 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116299156B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116593965B (zh) * | 2023-07-19 | 2023-09-29 | 中国海洋大学 | 一种基于远场机会声源的水平阵阵型估计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109270518A (zh) * | 2018-10-29 | 2019-01-25 | 中国电子科技集团公司第三研究所 | 一种水声垂直矢量阵列声呐的多目标模拟方法、装置及声呐 |
CN109283492A (zh) * | 2018-10-29 | 2019-01-29 | 中国电子科技集团公司第三研究所 | 多目标方位估计方法及水声垂直矢量阵列*** |
CN111398902A (zh) * | 2020-03-11 | 2020-07-10 | 哈尔滨工程大学 | 一种直角三角形三元水听器阵被动测距测向方法 |
CN112098938A (zh) * | 2020-08-31 | 2020-12-18 | 黑龙江工程学院 | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 |
CN113495275A (zh) * | 2021-05-24 | 2021-10-12 | 中国海洋大学 | 一种单水听器垂直合成孔径被动定位方法、***及应用 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9835764B2 (en) * | 2014-09-25 | 2017-12-05 | The United States Of America As Represented By The Secretary Of The Navy | System and method for the calibration of a hydrophone line array |
-
2023
- 2023-05-26 CN CN202310601868.6A patent/CN116299156B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109270518A (zh) * | 2018-10-29 | 2019-01-25 | 中国电子科技集团公司第三研究所 | 一种水声垂直矢量阵列声呐的多目标模拟方法、装置及声呐 |
CN109283492A (zh) * | 2018-10-29 | 2019-01-29 | 中国电子科技集团公司第三研究所 | 多目标方位估计方法及水声垂直矢量阵列*** |
CN111398902A (zh) * | 2020-03-11 | 2020-07-10 | 哈尔滨工程大学 | 一种直角三角形三元水听器阵被动测距测向方法 |
CN112098938A (zh) * | 2020-08-31 | 2020-12-18 | 黑龙江工程学院 | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 |
CN113495275A (zh) * | 2021-05-24 | 2021-10-12 | 中国海洋大学 | 一种单水听器垂直合成孔径被动定位方法、***及应用 |
Non-Patent Citations (1)
Title |
---|
史文佳.基于垂直阵的差频匹配场声源距离估计方法.《声学技术》.2019,第第38卷卷(第第2期期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN116299156A (zh) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6924221B2 (ja) | 合成開口ソナーのためのシステムおよび方法 | |
CN113702909B (zh) | 一种基于声音信号到达时间差的声源定位解析解计算方法及装置 | |
CN109239712B (zh) | 基于水下声场和声能流的噪声探测方法 | |
CN116299156B (zh) | 一种水听器垂直阵阵元位置估计方法及其优化策略 | |
CN108845325A (zh) | 拖曳线列阵声纳子阵误差失配估计方法 | |
WO2009048683A2 (en) | Controlling seismic source elements based on determining a three-dimensional geometry of the seismic source elements | |
CN113884986B (zh) | 波束聚焦增强的强冲击信号空时域联合检测方法及*** | |
Battle et al. | Geoacoustic inversion of tow-ship noise via near-field-matched-field processing | |
Zhao et al. | Open‐Lake Experimental Investigation of Azimuth Angle Estimation Using a Single Acoustic Vector Sensor | |
CN112098938B (zh) | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 | |
Gerstoft et al. | Adaptive beamforming of a towed array during a turn | |
CN109001297B (zh) | 基于单矢量水听器的大样本水声材料声反射系数测量方法 | |
CN116879901B (zh) | 一种矢量水听器垂直阵列深海宽带声源深度估计方法 | |
CN111551942B (zh) | 一种基于解卷积算法的水下自主航行器回坞方法 | |
US9897712B2 (en) | Method of acquiring ghost-free signatures for broadband source calibration | |
CN116068493A (zh) | 一种深海大深度垂直分布式水听器被动声源定位方法 | |
Santos et al. | Dual accelerometer vector sensor mounted on an autonomous underwater vehicle (AUV)-Experimental results | |
CN107202975A (zh) | 一种二维矢量阵阵元姿态误差校正方法 | |
Kerstens et al. | An optimized planar MIMO array approach to in-air synthetic aperture sonar | |
Touzé et al. | Double-Capon and double-MUSICAL for arrival separation and observable estimation in an acoustic waveguide | |
CN115390052A (zh) | 一种成像声呐实物模拟装置及方法 | |
KR101135456B1 (ko) | 수동 소나의 센서 신호 모의 장치 | |
CN113126029A (zh) | 适用于深海可靠声路径环境的多传感器脉冲声源定位方法 | |
CN113126030A (zh) | 基于宽带声场干涉结构的深海直达声区目标深度估计方法 | |
Sinha et al. | Study of acoustic vector sensor based direction of arrival estimation of in-air maneuvering tonal source |
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 |