CN109991567B - 一种水下滑翔机四面体阵三维被动测向方法 - Google Patents
一种水下滑翔机四面体阵三维被动测向方法 Download PDFInfo
- Publication number
- CN109991567B CN109991567B CN201910295799.4A CN201910295799A CN109991567B CN 109991567 B CN109991567 B CN 109991567B CN 201910295799 A CN201910295799 A CN 201910295799A CN 109991567 B CN109991567 B CN 109991567B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- coordinate system
- hydrophones
- signal
- signals
- 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
一种水下滑翔机四面体阵三维被动测向方法,属于探测信号处理领域。水下滑翔机平台无法对水下线谱声信号进行三维测向。本发明将四个水听器两两组合为各组合在直角坐标系中的向量;求取各向量与入射波方向向量所成余弦夹角对应的信号到达两个水听器的理论时延差;水听器的接收信号转换为频域信号,频域信号互谱计算结果取模平均得到互谱幅度谱平均值,之后经线谱检测得信号频率;由信号频率得到各组水听器的线谱信号的相位差;线谱频率及线谱信号的相位差求得实测时延差;两组时延差对应联立,得载体坐标系下的信号入射方向向量;求解二元方程组即得大地坐标系下的方位角和俯仰角。本发明能够测得目标在大地坐标系下的三维方位角。且测向过程简单。
Description
技术领域
本发明涉及水下滑翔机声学探测信号处理领域,特别涉及一种利用水下滑翔机装载的四面体阵对线谱声信号进行三维测向的方法。
背景技术
深海由于环境特殊性要实现对水下声目标被动探测和海洋环境噪声特性测量对探测系统平台有着严格的要求。深海水下滑翔机平台则是能够满足上述要求的良好载体。深海声传播存在明显的声汇聚和影区效应,系统的探测性能受工作深度的影响较大,探测系统要具有多个深度声信号采集记录能力,以保证探测设备可以对位于不同距离的目标进行有利探测。水下滑翔机作为新型水下探测平台,其具有水下工作时间长、可变深度采集声学信息、自控能力强、隐蔽性好、信息传输便捷等特点,在水声探测和声学海洋环境监测中可发挥重要的作用。水下滑翔机平台探测系统的缺点在于尺寸小,搭载能力差,电力有限,因此无法装载大型声学探测系统。这里给出了一种利用水下滑翔机平台搭载的小型四元阵构成的四面体阵对线谱声信号进行三维测向的方法。利用四面体阵的两两阵元做互谱,通过互谱法对线谱在各阵元间的相位差进行估计,从而得到目标相对于四元阵坐标系的方位信息,包括方位角信息和俯仰角信息。再利用水下滑翔机平台搭载的罗经测得的姿态信息通过三维坐标转换得到目标相对于滑翔机测量点的大地坐标系方位和俯仰角信息。实现对线谱声信号进行三维测向。
经过文献检索发现,有下列文献对互谱测向方法和罗经修正方法进行了研究:
李霞,李志舜.利用回波互谱进行方位估计的两种方法研究.船舶工程.2002,(3):55-57页.(以下简称文献1)
张新春,崔希民.轨道几何状态惯性检测系统零速姿态修正方法[J].测绘科学,2017,42(07):17-21.(以下简称文献2)
文献1研究了两种经典的互谱测向方法,一种是经典的互谱测向方法,一种是基于小波变换的互谱法。通过文献1的描述知道互谱法可以实现对目标的测向。
文献2给出了一种坐标修正的方法,表明在已知本地坐标系和大地坐标系之间的夹角情况下,可以通过坐标转换的方法将相对于本地坐标系下的测量结果转换为大地坐标结果。
本发明采用了文献1中的互谱相位测量方法。但是上述文献都未就双元线谱测向方位模糊问题给出解决办法,也未给出包含俯仰角测量的三维自动检测方法,仅仅是对方位等测量方法的研究,无法解决线谱目标的自动检测问题。
发明内容
本发明的目的是为了解决水下滑翔机平台无法对水下线谱声信号进行三维测向问题,提供了一种利用装载于滑翔机平台的小型四面体阵对线谱声信号进行三维测向方法。
一种水下滑翔机四面体阵三维被动测向方法,所述测向方法是通过以下步骤实现:
步骤2、求得步骤1中所得的各余弦夹角对应的信号到达两个水听器的理论时延差τ1、τ2、τ3、τ4、τ5、τ6;
步骤3、四个水听器的接收信号经傅里叶变换转换为频域信号X1(f)、X2(f)、X3(f)、X4(f),频域信号根据步骤1中的两两组合方式分别进行互谱计算,得到互谱计算结果X12(f)、X13(f)、X14(f)、X23(f)、X24(f)、X34(f);
步骤4、将步骤3测得的互谱结果取模平均,得到互谱幅度谱平均值,再对互谱幅度谱平均值进行线谱检测,测得信号的频率f0;
步骤6、根据步骤4、5中所得频率和各组水听器的线谱信号的相位差,求得各组水听器的线谱信号的实测时延差τ′1、τ′2、τ′3、τ′4、τ′5、τ′6;
步骤7、将步骤2与步骤6所求的时延差一一对应联立起来,得到一个由6个方程组成的二元超定方程组;
本发明的有益效果为:
1.文献1中给出的方法是传统的互谱测向方法,主要适用于两元阵阵元间距不大于测量频段上限半波间距的情况,只能对二维的波达方向角进行估计,且存在左右舷模糊;
2.文献2通过罗经数据对惯性导航系统下的速度进行修正,提高轨道几何状态检测仪长时间检测轨道的精度。本发明先在载体坐标系下进行来波方向估计,最后对来波方向向量进行姿态修正得到大地坐标系下的三维方位估计结果。
本发明利用四面体阵的布阵方式,能够克服传统线谱目标互谱测向方法只能测得二维方向角且存在测向左右舷模糊的弊端,本发明通过阵元两两互谱求时延,再联立解方程组,最后进行罗经姿态修正,能够测得目标在大地坐标系下的三维方位角。且测向过程简单。
附图说明
图1载体坐标系下的阵列模型示意图;
图2是四面体阵线谱三维测向方法处理框图;
图3是线谱检测流程框图;
图4是四面体阵互谱法方位角理论值与估计值对比结果图;
图5是四面体阵互谱法俯仰角理论值与估计值对比结果图。
具体实施方式
具体实施方式一:
本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述测向方法是通过以下步骤实现:
步骤3、四个水听器的接收信号经傅里叶变换(FFT)转换为频域信号X1(f)、X2(f)、X3(f)、X4(f),频域信号根据步骤1中的两两组合方式分别进行互谱计算,得到互谱计算结果X12(f)、X13(f)、X14(f)、X23(f)、X24(f)、X34(f);
步骤4、将步骤3测得的互谱结果取模平均,得到互谱幅度谱平均值,再对互谱幅度谱平均值进行线谱检测,测得信号的频率f0;
步骤6、根据步骤4、5中所得频率和各组水听器的线谱信号的相位差,求得各组水听器的线谱信号的实测时延差τ′1、τ′2、τ′3、τ′4、τ′5、τ′6;
步骤7、将步骤2与步骤6所求的时延差一一对应联立起来,得到一个由6个方程组成的二元超定方程组;
具体实施方式二:
步骤1-1、设以水下滑翔机平台为载体建立三维直角坐标系xyz即载体坐标系,示意图如图1所示,四个水听器所在位置分别为点A、B、C、D,L1为3号水听器到AB的垂直距离,L2为4号水听器到O的垂直距离,L3为4号水听器到中轴线的垂直距离,L4为两翼AB的距离,R为Glider的半径,其中A、B位于y轴且关于原点对称,坐标分别为A(0,L4/2,0),B(0,-L4/2,0),C、D点均位于平面xOz上,坐标分别为C(L1,0,R),D(-L2,0,L3);Glider表示滑翔机;
上式中,·表示向量的点乘,|·|为对向量取模;
同理,求得其余余弦夹角α2、α3、α4、α5和α6的值。
具体实施方式三:
与具体实施方式二不同的是,本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述的步骤2中,求得步骤1中所得的各向量与入射波方向向量所成的余弦夹角对应的信号到达两个水听器的理论时延差τ1、τ2、τ3、τ4、τ5、τ6的过程,具体为:
四个水听器A、B、C、D构成四面体阵,各水听器分别作为阵元,四面体阵与实验目标满足远场平面波条件,示意图如图2所示,阵元A和B之间,在已知入射信号与阵元所在直线的夹角α1的情况下,信号到达两个水听器A和B的时延差τ1为:
其中,c为水中声速;
同理,求得其余5组水听器的信号到达时延差τ2、τ3、τ4、τ5、τ6。
具体实施方式四:
与具体实施方式三不同的是,本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述的步骤3中,四个水听器的接收信号经傅里叶变换(FFT)转换为频域信号X1(f)、X2(f)、X3(f)、X4(f),频域信号根据步骤1中的两两组合方式分别进行互谱计算,得到互谱计算结果X12(f)、X13(f)、X14(f)、X23(f)、X24(f)、X34(f)的过程,具体为:
首先,对接收信号x1(t)、x2(t)、x3(t)、x4(t)分别做傅里叶变换,得到接收信号的频域信号X1(f)、X2(f)、X3(f)、X4(f),再根据步骤1中的两两组合方式将所得频域信号做互谱即得到X12(f)、X13(f)、X14(f)、X23(f)、X24(f)、X34(f)。
具体实施方式五:
与具体实施方式四不同的是,本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述的步骤4中,将步骤3测得的互谱结果取模平均,得到互谱幅度谱平均值,再对互谱幅度谱平均值进行线谱检测,测得信号的频率f0的过程,具体为:
如图3所示,首先,根据步骤3得到的互谱幅度谱取模平均得到|X(f)|,再对|X(f)|进行双向α滤波提取连续谱背景;
然后,在连续谱背景的基础上加上门限值DT得到检测背景曲线,判断|X(f)|中是否存在高于连续谱背景曲线的线谱,
若是,则判定为检测到线谱目标,并取最大值位置对应的频率为该目标线谱频率f0;
否则,认为未检测到目标。
具体实施方式六:
与具体实施方式五不同的是,本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述的步骤5中,根据步骤3、4所得信号频率f0和其所对应的相位谱进行相位差估计,得到各组水听器的线谱信号的相位差的过程,具体为:
已知互谱X(f)与信号中心频率f0,将所得的各互谱结果X12(f)、X13(f)、X14(f)、X23(f)、X24(f)、X34(f)分别带入下式所示的两水听器的接收信号的相位差的计算公式,即可求得到各组水听器的线谱信号的相位差其中,相位差的计算公式为:
具体实施方式七:
与具体实施方式六不同的是,本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述的步骤6中,根据步骤4、5中所得频率和各组水听器的线谱信号的相位差,求得各组水听器的线谱信号的实测时延差τ′1、τ′2、τ′3、τ′4、τ′5、τ′6的过程,具体为:
其中ω为信号的角频率且ω=2πf0;
具体实施方式八:
与具体实施方式七不同的是,本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述的步骤7中,将步骤2与步骤6所求的时延差一一对应联立起来,得到一个由6个方程组成的二元超定方程组的过程,具体为:
首先,联立步骤3与步骤7所得的时延差结果:τ1=τ′1、τ2=τ′2、τ3=τ′3、τ4=τ′4、τ5=τ′5、τ6=τ′6;
然后,整理成矩阵方程组Ax=b的形式,为:
其中,对于空间波数k有k=ω/c,则系数矩阵A与常数项矩阵b分别为:
具体实施方式九:
与具体实施方式八不同的是,本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述的步骤8中,用最小二乘法解步骤7所得的超定方程组,得载体坐标系下的信号入射方向向量并利用罗经值对信号入射方向向量进行修正,得到大地坐标系下的信号入射方向向量的过程,具体为:
按照最小二乘原理,方程组Ax=b的解满足:
式中,(·)-1表示对矩阵求逆;
式中,k′x、k′y、k′z分别为大地坐标系下k在x轴、y轴与z轴上的投影;
步骤8-3、将三维罗经值的Heading、Pitch、Roll分别用α、β、γ表示,通过旋转矩阵将信号入射方向向量转换到大地坐标系下,得:
具体实施方式十:
与具体实施方式九不同的是,本实施方式的一种水下滑翔机四面体阵三维被动测向方法,所述的步骤9中,根据步骤8的结果,将大地坐标系下的方位角θ和俯仰角联立构成二元方程组,求解二元方程组即得大地坐标系下的方位角θ和俯仰角的过程,具体为:
仿真实验:
图4、图5是对一个单频运动目标利用四面体阵线谱三维测向方法的测向结果,仿真条件四面体阵的尺寸为L1=2080mm,L2=36mm,L3=250mm,L4=1500mm,R=110.5mm,信号频率为100Hz,连续波信号,信噪比为20dB,图4为水平方位角测向结果,图5为垂直俯仰角的测向结果,图中“*”为仿真测向结果,“.”为理论角度。从图中结果可以看到,利用四面体阵线谱三维测向方法可以实现对目标的三维方向角的测量。
本发明还可有其它多种实施案例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明做出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种水下滑翔机四面体阵三维被动测向方法,其特征在于:所述测向方法是通过以下步骤实现:
步骤2、求得步骤1中所得的各余弦夹角对应的信号到达两个水听器的理论时延差τ1、τ2、τ3、τ4、τ5、τ6;
步骤3、四个水听器的接收信号经傅里叶变换转换为频域信号X1(f)、X2(f)、X3(f)、X4(f),频域信号根据步骤1中的两两组合方式分别进行互谱计算,得到互谱计算结果X12(f)、X13(f)、X14(f)、X23(f)、X24(f)、X34(f);
步骤4、将步骤3测得的互谱结果取模平均,得到互谱幅度谱平均值,再对互谱幅度谱平均值进行线谱检测,测得信号的频率f0;
步骤6、根据步骤4、5中所得频率和各组水听器的线谱信号的相位差,求得各组水听器的线谱信号的实测时延差τ′1、τ′2、τ′3、τ′4、τ′5、τ′6;
步骤7、将步骤2与步骤6所求的时延差一一对应联立起来,得到一个由6个方程组成的二元超定方程组;
步骤1-1、设以水下滑翔机平台为载体建立三维直角坐标系xyz即载体坐标系,四个水听器所在位置分别为点A、B、C、D,L1为3号水听器到AB的垂直距离,L2为4号水听器到O的垂直距离,L3为4号水听器到中轴线的垂直距离,L4为两翼AB的距离,R为Glider的半径,其中A、B位于y轴且关于原点对称,坐标分别为A(0,L4/2,0),B(0,-L4/2,0),C、D点均位于平面xOz上,坐标分别为C(L1,0,R),D(-L2,0,L3);Glider表示滑翔机;
同理,求得其余余弦夹角α2、α3、α4、α5和α6的值。
4.根据权利要求3所述的一种水下滑翔机四面体阵三维被动测向方法,其特征在于:所述的步骤3中,四个水听器的接收信号经傅里叶变换转换为频域信号X1(f)、X2(f)、X3(f)、X4(f),频域信号根据步骤1中的两两组合方式分别进行互谱计算,得到互谱计算结果X12(f)、X13(f)、X14(f)、X23(f)、X24(f)、X34(f)的过程,具体为:
首先,对接收信号x1(t)、x2(t)、x3(t)、x4(t)分别做傅里叶变换,得到接收信号的频域信号X1(f)、X2(f)、X3(f)、X4(f),再根据步骤1中的两两组合方式将所得频域信号做互谱即得到X12(f)、X13(f)、X14(f)、X23(f)、X24(f)、X34(f)。
5.根据权利要求4所述的一种水下滑翔机四面体阵三维被动测向方法,其特征在于:所述的步骤4中,将步骤3测得的互谱结果取模平均,得到互谱幅度谱平均值,再对互谱幅度谱平均值进行线谱检测,测得信号的频率f0的过程,具体为:
首先,根据步骤3得到的互谱幅度谱取模平均得到|X(f)|,再对|X(f)|进行双向α滤波提取连续谱背景;
然后,在连续谱背景的基础上加上门限值DT得到检测背景曲线,判断|X(f)|中是否存在高于连续谱背景曲线的线谱,
若是,则判定为检测到线谱目标,并取最大值位置对应的频率为该目标线谱频率f0;
否则,认为未检测到目标。
9.根据权利要求8所述的一种水下滑翔机四面体阵三维被动测向方法,其特征在于:所述的步骤8中,用最小二乘法解步骤7所得的超定方程组,得载体坐标系下的信号入射方向向量并利用罗经值对信号入射方向向量进行修正,得到大地坐标系下的信号入射方向向量的过程,具体为:
按照最小二乘原理,方程组Ax=b的解满足:
式中,(·)-1表示对矩阵求逆;
式中,k′x、k′y、k′z分别为大地坐标系下k在x轴、y轴与z轴上的投影;
步骤8-3、将三维罗经值的Heading、Pitch、Roll分别用α、β、γ表示,通过旋转矩阵将信号入射方向向量转换到大地坐标系下,得:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910295799.4A CN109991567B (zh) | 2019-04-12 | 2019-04-12 | 一种水下滑翔机四面体阵三维被动测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910295799.4A CN109991567B (zh) | 2019-04-12 | 2019-04-12 | 一种水下滑翔机四面体阵三维被动测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109991567A CN109991567A (zh) | 2019-07-09 |
CN109991567B true CN109991567B (zh) | 2021-03-09 |
Family
ID=67133508
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910295799.4A Active CN109991567B (zh) | 2019-04-12 | 2019-04-12 | 一种水下滑翔机四面体阵三维被动测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109991567B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110543179B (zh) * | 2019-08-21 | 2020-08-18 | 中国科学院声学研究所 | 一种基于三维稀疏阵列声源方位识别的水面目标躲避方法 |
CN110703199B (zh) * | 2019-10-22 | 2021-09-24 | 哈尔滨工程大学 | 基于罗经补偿的四元十字阵高精度方位估计方法 |
CN111398902B (zh) * | 2020-03-11 | 2022-03-22 | 哈尔滨工程大学 | 一种直角三角形三元水听器阵被动测距测向方法 |
CN111504253A (zh) * | 2020-04-14 | 2020-08-07 | 哈尔滨工程大学 | 一种波浪滑翔机的水下三维声学立体阵确定方法 |
CN111505568B (zh) * | 2020-04-14 | 2022-06-21 | 哈尔滨工程大学 | 基于差分优化的四面体阵目标方位估计方法 |
CN111610492A (zh) * | 2020-06-03 | 2020-09-01 | 电子科技大学 | 一种多声传感器阵列智能感知方法及系统 |
CN111679246B (zh) * | 2020-06-04 | 2022-11-15 | 哈尔滨工程大学 | 一种三维运动平台搭载的小尺度阵高分辨测向方法 |
CN113702960B (zh) * | 2021-06-29 | 2022-10-21 | 哈尔滨工程大学 | 一种基于时延和多普勒频移的水下机动平台高精度测速方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2191405C1 (ru) * | 2001-03-11 | 2002-10-20 | Федеральное государственное унитарное предприятие "Центральный научно-исследовательский институт "Морфизприбор" | Способ определения радиальной скорости объекта |
JP2005172760A (ja) * | 2003-12-15 | 2005-06-30 | Mitsubishi Electric Corp | 方位探知装置 |
CN104237871B (zh) * | 2013-06-08 | 2017-01-11 | 中国科学院声学研究所 | 一种基于相位补偿的时延差估计方法 |
CN104977565A (zh) * | 2014-04-02 | 2015-10-14 | 中国科学院声学研究所 | 一种基于分布式多阵列实现的水下目标三维定位方法 |
CN105182345B (zh) * | 2015-09-26 | 2017-10-03 | 哈尔滨工程大学 | 基于单矢量潜标的连续谱信号目标自动检测方法 |
CN106324559B (zh) * | 2016-08-30 | 2019-03-12 | 西安空间无线电技术研究所 | 一种大基线四元阵宽带信号测向系统及方法 |
CN109540172A (zh) * | 2018-12-27 | 2019-03-29 | 中国船舶重工集团公司第七0研究所 | 一种用于水雷平台的目标运动参数估计方法 |
-
2019
- 2019-04-12 CN CN201910295799.4A patent/CN109991567B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109991567A (zh) | 2019-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109991567B (zh) | 一种水下滑翔机四面体阵三维被动测向方法 | |
CN109737956B (zh) | 一种基于双应答器的sins/usbl相位差紧组合导航定位方法 | |
CN109884730B (zh) | 一种基于深海水下滑翔机平台的虚拟阵被动探测方法 | |
CN103438987A (zh) | 基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法 | |
CN108614268B (zh) | 低空高速飞行目标的声学跟踪方法 | |
CN103792513B (zh) | 一种雷声定位系统及方法 | |
CN108845325A (zh) | 拖曳线列阵声纳子阵误差失配估计方法 | |
CN108169731A (zh) | 基于单近场校正源的拖曳线列阵阵形估计方法和装置 | |
CN101915922A (zh) | 拖曳线列阵被动测距方法 | |
CN102645265A (zh) | 一种基于虚拟时间反转镜的舰船辐射噪声级测量方法 | |
CN111025273B (zh) | 一种畸变拖曳阵线谱特征增强方法及系统 | |
CN109855649B (zh) | 一种水下滑翔机平台运动轨迹被动确定方法 | |
CN105116372B (zh) | 一种浅海海底水平直线水听器阵朝向校准的方法 | |
CN106249196A (zh) | 三分量声矢量传感器稀疏阵列四元数解模糊方法 | |
CN110389318B (zh) | 一种基于立体六元阵的水下移动平台定位系统及方法 | |
CN104714235A (zh) | 一种低频双矢量水听器阵列的测距方法及系统 | |
CN109975813A (zh) | 一种基于水下滑翔机平台的十字阵线谱目标被动探测方法 | |
CN108562872A (zh) | 一种超短基线水下声学定位校准时异常值探测的方法 | |
CN104515974B (zh) | 微波着陆机载设备角度、测距数据处理方法 | |
CN102721966A (zh) | 相干类测深声纳正下方高精度测深方法及系统 | |
CN103513238A (zh) | 一种规整化最小二乘子空间相交的目标方位测向方法 | |
CN106908754A (zh) | L型声矢量传感器阵列esprit解相干参数估计方法 | |
CN111679246B (zh) | 一种三维运动平台搭载的小尺度阵高分辨测向方法 | |
CN111505568B (zh) | 基于差分优化的四面体阵目标方位估计方法 | |
CN202093166U (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 |