CN114563778B - 基于时延冗余测量的阵形校准方法 - Google Patents
基于时延冗余测量的阵形校准方法 Download PDFInfo
- Publication number
- CN114563778B CN114563778B CN202210168090.XA CN202210168090A CN114563778B CN 114563778 B CN114563778 B CN 114563778B CN 202210168090 A CN202210168090 A CN 202210168090A CN 114563778 B CN114563778 B CN 114563778B
- Authority
- CN
- China
- Prior art keywords
- array
- sound source
- auxiliary sound
- array element
- time delay
- 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 46
- 238000005259 measurement Methods 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 50
- 239000013598 vector Substances 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000001514 detection method Methods 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 7
- 238000009825 accumulation Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000000342 Monte Carlo simulation Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- 238000004880 explosion Methods 0.000 description 2
- 241001386813 Kraken Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000012736 aqueous medium Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000013519 translation 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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/78—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 electromagnetic waves other than radio waves
- G01S3/781—Details
-
- 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/82—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 with means for adjusting phase or compensating for time-lag errors
-
- 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/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了基于时延冗余测量的阵形校准方法,该方法包括:步骤1)分别针对第一辅助声源S1和第二辅助声源S2,依次进行步骤2)‑步骤3):步骤2)对待校准阵列A的所有阵元接收的一个辅助声源的信号进行带通滤波处理;步骤3)对经带通滤波处理后的待校准阵列A的两两阵元组合信号进行相关时延估计,得到相对参考阵元的时延向量;步骤4)根据分别针对第一辅助声源S1和第二辅助声源S2得到的时延向量,建立待校准阵列A每个阵元的位置反演目标函数;步骤5)通过优化搜索得到每个阵元的坐标,从而完成阵列A的校准。本发明通过充分利用阵元空间信息,提高了阵形校准精度与稳定性,采用该方法校准阵形后的阵列探测目标,得到的声源方位角精度更高。
Description
技术领域
本发明属于声纳信号处理、阵列信号处理技术领域,尤其涉及基于时延冗余测量的阵形校准方法。
背景技术
阵形校准精度对目标测向性能有着显著影响,布放于海底的柔性水平阵,由于受到布放、内波、浪涌等因素的影响,实际阵形往往与预设阵形有较大差异。进行阵列形状校准工作,可以提高后续信号处理算法的性能,尤其是各种自适应算法的性能。已有的阵形校准方法大致分为两类:一类为无源校准,无需额外的校准声源,利用噪声源就可以实现,但是反演出的阵列只知道相对形状,需要根据其他信息对初始结果进行平移、旋转、翻转等操作才能获取各阵元的真实位置;另一类为有源校准,通过利用校准声源,结合声源与参考阵元的GPS位置信息,测量各相邻阵元间接收信号的时间延迟,从而获得各阵元相对参考阵元的真实位置;但由于浅海信道存在纵向相关振荡现象、信道随机起伏,以及受到海洋环境噪声等因素影响,尤其在校准信号信噪比较低时,传统方法会因为时延估计误差累积造成校准阵形的偏差。
发明内容
本发明的目的在于克服现有技术缺陷,提出了基于时延冗余测量的阵形校准方法。
为了实现上述目的,本发明提出了一种基于时延冗余测量的阵形校准方法,所述方法包括:
步骤1)分别针对第一辅助声源S1和第二辅助声源S2,依次进行步骤2)-步骤3):
步骤2)对待校准阵列A的所有阵元接收的一个辅助声源的信号进行带通滤波处理;
步骤3)对经带通滤波处理后的待校准阵列A的两两阵元组合信号进行相关时延估计,得到相对参考阵元的时延向量;
步骤4)根据分别针对第一辅助声源S1和第二辅助声源S2得到的时延向量,建立待校准阵列A每个阵元的位置反演目标函数;
步骤5)通过优化搜索得到每个阵元的坐标,从而完成阵列A的校准。
作为上述方法的一种改进,所述待校准阵列为I元的坐底直线阵,其中第1个阵元为参考阵元O,其位置坐标为O(x1,y1)。
作为上述方法的一种改进,所述方法还包括先后布设第一辅助声源S1和第二辅助声源S2的步骤,具体包括:
第一辅助声源S1的位置为S1(X1,Y1),第二辅助声源S2的位置为S2(X2,Y2),第一辅助声源S1与参考阵元O的连线为OS1,第二辅助声源S2与参考阵元O的连线为OS2,OS1与OS2垂直。
作为上述方法的一种改进,所述步骤3)具体包括:
步骤301)根据下式,通过相关峰最大值的位置,得到第i个阵元Ai较第j个阵元Aj接收到第k个辅助声源的信号时延估计值为:
其中,为第j个阵元Aj接收到第k个辅助声源的时域信号,/>为第i个阵元Ai接收到第k个辅助声源时间延迟/>后的时域信号,上标*表示取共轭,t表示时间;
步骤302)根据步骤301)得到的全部阵元组合,由下式得到接收第k个辅助声源时,两两阵元间的时延估计矩阵Gk为:
其中, 为第j个阵元Aj较第i个阵元Ai接收到第k个辅助声源的信号时延估计值,当j=i时,/>
步骤303)剔除主对角线元素,取两两阵元间的时延估计矩阵Gk的上三角元素得到矩阵上标T表示转置;
基于最小二乘准则得到矩阵Q为:
根据下式得到待校准阵列A接收第k个辅助声源时相对参考阵元O的时延向量P为:
P=Q+B
其中, 为第i个阵元Ai接收第k个辅助声源时相对参考阵元O的时延量,上标+表示取矩阵伪逆。
作为上述方法的一种改进,所述步骤4)具体包括:
建立待校准阵列A各个阵元的位置反演目标函数,从时延向量P中获取k=1,2,则第i个阵元Ai的坐标J(xi,yi)为:
其中,c0为阵形校准参考声速。
作为上述方法的一种改进,所述步骤5)具体包括:
通过优化搜索得第i个阵元Ai的坐标为:
与现有技术相比,本发明的优势在于:
1、在信噪比偏低时,所提的基于时延冗余测量的阵形校准方法,更充分地利用了阵元空间信息,从而提高了阵形校准精度与稳定性;
2、用校准阵形后的阵列探测目标,得到的声源方位角精度更高。
附图说明
图1是本发明的校准方法位置示意图;
图2是声速剖面;
图3是信噪比10dB时的仿真阵形校准结果,其中图3(a)是阵形真实位置,图3(b)是采用常规相邻累积的时延估计阵形校准结果,图3(c)是采用所提冗余测量的时延估计阵形校准结果;
图4是不同信噪比下阵形校准误差的均值与方差;
图5是实验数据,其中图5(a)是各阵元接收的第一辅助声源S1带限信号时域波形的实验数据,图5(b)是各阵元接收的第二辅助声源S2带限信号时域波形;
图6是实验数据的阵形校准结果,其中图6(a)是阵形预设的布设位置,图6(b)是采用常规相邻累积的时延估计阵形校准结果,图6(c)是采用所提冗余测量的时延估计阵形校准结果;
图7是不同阵形数据对测试声源信号的CBF方位估计结果,其中图7(a)是预设阵形CBF波束扫描图,图7(b)是所提时延冗余测量后校准阵形的CBF波束扫描图,图7(c)是常规时延相邻累积测量后校准阵形的CBF波束扫描图;
图8(a)是加噪声后各阵元接收第一辅助声源S1信号时域波形,图8(b)是加噪声后各阵元接收第二辅助声源S2信号时域波形;
图9是50次蒙特卡洛仿真实验后的阵形校准,其中图9(a)是时延冗余测量,图9(b)是时延累积测量。
具体实施方式
根据简正波理论,在水平不变分层的海洋波导环境中,点声源激发的单频声场可以表示为M阶简正波的线性叠加
其中,krm为第m阶模态的水平波数,Ψm为模态深度函数,ρ(zs)为声源处水介质密度,zs和z分别为声源深度和接收深度,r为阵元与声源的水平距离。假定频带[fL,fH]内声源信号频谱为S(f),则第i号阵元接收该频段信号的时域波形为:
其中,ri为第i号阵元与声源的水平距离。
如图1所示,S1、S2分别为校准声源的位置,O为海底水平阵参考阵元的位置,Ai为某待校准阵元的位置。
已知声源S1(X1,Y1)、S2(X2,Y2)与参考阵元O(x1,y1)的GPS位置信息,通过估计各阵元较参考阵元的信号到达时延,以求解Ai点坐标(xi,yi)。
其中为第i阵元较参考阵元的信号到达时延,c0为阵形校准参考声速。选取Ai点坐标反演目标函数如下
通过优化搜索可得Ai点坐标:
任意两个阵元间的信号时延可以根据相关峰的位置来估计:
其中,上标*表示取共轭,最大相关系数对应的即为阵元Ai较阵元Aj接收到声源Sk的信号时延估计值。测量全部阵元组合后,可得两两阵元间的时延估计矩阵Gk:
其中,Gk需要估计的阵元间时延变量总数为/>
当信噪比偏低、或者声场信道存在起伏的情况下,仅仅测量(I-1)组相邻阵元的时延并没有充分利用信号的空间信息。
通过全部两两阵元组合的信号相关,开展全部时延估计,进而由最小二乘方法计算第i阵元相对参考阵元的时延量,可实现对时延的更好估计,从而更准确地估计阵形。
在全部组两两阵元组合的信号相关时延估计完成后,建立相对参考阵元的时延向量P和两两阵元组合的信号时延向量B的线性关系如下
P=Q+B (8)
其中,上标+表示取矩阵伪逆;
对Gk采用最小二乘方法得到
由最小二乘方法计算第i阵元相对参考阵元的时延量,可实现对时延的更好估计。
下面结合附图和实施例对本发明的技术方案进行详细的说明。
实施例1
本发明的实施例提出了一种基于时延冗余测量的阵形校准方法。实施的步骤如下:
1、于阵列适当距离布设两信号到达方向近似垂直的辅助声源发射信号。记录声源S1(X1,Y1)、S2(X2,Y2)与参考阵元O(x1,y1)的位置坐标;
2、对待校准的阵列所有通道接收信号进行带通滤波处理;
3、对全部I(I-1)/2组两两阵元组合的信号相关时延开展估计,
4、根据两两阵元组合的信号时延向量B,求取建立相对参考阵元的时延向量P:P=Q+B
其中,
其中上标+表示取矩阵伪逆;
5、定义各个阵元的位置反演目标函数如下:
通过优化搜索可得第i阵元的坐标:
定义阵形校准误差
其中,为校准得到的位置坐标,(x0i,y0i)为位置坐标真值。
以下分别通过数值仿真和实验数据分析来对比验证本发明的技术效果。
●数值仿真
通过数值仿真来对比上述本发明和传统方法用于阵形校准的差异。声场仿真所用的信道参数如图2所示,接收阵位置处实测值声速剖面为温跃层结构,海深设为96m。海底假设为平坦半无限基底,声速1630m/s、密度1.76g/cm3、衰减系数0.33dB/λ。声源深度设为50m,接收阵列布放于海底。信号频段为20-30Hz,频域仿真分辨率为0.2Hz,使用Kraken程序计算声场,阵形估计中的参考声速选取为1568m/s。
阵列的坐标位置如图3中符号“o”所示,为10米间隔的11阵元直线阵。布设校准声源位置:S1在参考阵元O的0°方向,距离为20公里;S2在参考阵元O的270°方向,距离为20公里,其中以x轴正半轴为0°方向,逆时针方向旋转角度递增。需要说明的是S1和S2距离参考阵元O的距离可以不相等。
图3为信噪比10dB时的阵形校准结果:图3(a)阵形真实位置,用“o”表示,图3(b)用常规相邻累积的时延估计阵形校准结果,用“^”表示,图3(c)用所提冗余测量的时延估计阵形校准结果,用“.”表示。对比可以看出冗余测量的方法明显更好。
为考察不同信噪比下,两种时延估计方法的阵形校准精度,给接收信号添加不同大小的高斯白噪声,并进行100次蒙特卡洛仿真实验,统计分析两种时延估计方法的阵形校准结果的误差。将两种方法的上述误差均值和误差标准差绘图对比,如图4所示,“*”为冗余测量结果,“o”为累积测量结果。
从图4可以看出,在信噪比大于30dB时,时延冗余测量方法跟相邻累积求解方法都可以准确地实现阵形校准,而信噪比降低时,冗余测量的优势性则更为明显,不仅误差均值低,而且方差也更小。如以误差校准精度0.1为门限,冗余测量方法对信号的信噪比要求比累积测量方法要低约十分贝。
●实验数据分析
2020年9月某近海实验,水平接收阵布放于海底,阵列的信号采样率为16kHz,灵敏度为-170dB,实验船沿着测线投掷爆炸声源,标称爆炸深度为50m。实验期间的声速剖面、波导环境等参数与上文数值仿真所用参数一致。两辅助声源的布设位置与数值仿真接近,S1在参考阵元O的358.6°方向、距离19.67公里处;S2在参考阵元O的269.4°方向、距离20.78公里处。
对接收信号进行20-30Hz带通滤波处理,各阵元接收到的带内信号波形如图5所示,图5(a)为S1的接收信号,图5(b)为源S2的接收信号。
对滤波后信号分别用两种方法进行时延估计,得到阵形校准结果如图6所示,图6(a)阵形预设的布设位置,用“o”表示,图6(b)用常规相邻累积的时延估计阵形校准结果,用“^”表示,图6(c)用所提冗余测量的时延估计阵形校准结果,用“.”表示。可以看出校准后阵形与实验预设的直线阵形差别较大。为了验证校准阵形的有效性,对测试声源Stest[参考阵元O的269.4°方向、距离为21.98公里]的信号进行波束形成,并与GPS测量的结果进行比较。
图7给出了三种阵形数据对测试声源信号的CBF结果,图7(a)根据实验预设阵形CBF测量得到的方位角为272.4度,图7(b)冗余测量时延估计阵形校准后CBF测量得到的方位角为269.3度,图7(c)相邻累积时延估计阵形校准后CBF测量得到的方位角为269.2度,而根据CPS得到的方位角为269.4,可以看出校准后阵形CBF测向更准确,验证了所提方法的有效性。
从图5可以看出,由于各阵元接收信号的信噪比很高,两种阵形校准方法的差别并不明显。如图8所示,给接收信号添加10dB的高斯白噪声,并进行50次蒙特卡洛仿真实验,以分析两种方法的阵形校准结果的稳定性。结果如图9所示,当接收信号信噪比降低后,充分利用信号的空间信息进行冗余时延估计如图9(a),相较时延累积方法图如9(b)的阵形校准结果更稳定。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种基于时延冗余测量的阵形校准方法,所述方法包括:
步骤1)分别针对第一辅助声源S1和第二辅助声源S2,依次进行步骤2)-步骤3):
步骤2)对待校准阵列A的所有阵元接收的一个辅助声源的信号进行带通滤波处理;
步骤3)对经带通滤波处理后的待校准阵列A的两两阵元组合信号进行相关时延估计,得到相对参考阵元的时延向量;
步骤4)根据分别针对第一辅助声源S1和第二辅助声源S2得到的时延向量,建立待校准阵列A每个阵元的位置反演目标函数;
步骤5)通过优化搜索得到每个阵元的坐标,从而完成阵列A的校准;
所述待校准阵列为I元的坐底直线阵,其中第1个阵元为参考阵元O,其位置坐标为O(x1,y1);
所述方法还包括先后布设第一辅助声源S1和第二辅助声源S2的步骤,具体包括:
第一辅助声源S1的位置为S1(X1,Y1),第二辅助声源S2的位置为S2(X2,Y2),第一辅助声源S1与参考阵元O的连线为OS1,第二辅助声源S2与参考阵元O的连线为OS2,OS1与OS2垂直;
所述步骤3)具体包括:
步骤301)根据下式,通过相关峰最大值的位置,得到第i个阵元Ai较第j个阵元Aj接收到第k个辅助声源的信号时延估计值为:
其中,为第j个阵元Aj接收到第k个辅助声源的时域信号,/>为第i个阵元Ai接收到第k个辅助声源时间延迟/>后的时域信号,上标*表示取共轭,t表示时间;
步骤302)根据步骤301)得到的全部阵元组合,由下式得到接收第k个辅助声源时,两两阵元间的时延估计矩阵Gk为:
其中,为第j个阵元Aj较第i个阵元Ai接收到第k个辅助声源的信号时延估计值,当j=i时,/>
步骤303)剔除主对角线元素,取两两阵元间的时延估计矩阵Gk的上三角元素得到矩阵上标T表示转置;
基于最小二乘准则得到矩阵Q为:
根据下式得到待校准阵列A接收第k个辅助声源时相对参考阵元O的时延向量P为:
P=Q+B
其中,为第i个阵元Ai接收第k个辅助声源时相对参考阵元O的时延量,上标+表示取矩阵伪逆。
2.根据权利要求1所述的基于时延冗余测量的阵形校准方法,其特征在于,所述步骤4)具体包括:
建立待校准阵列A各个阵元的位置反演目标函数,从时延向量P中获取k=1,2,则第i个阵元Ai的坐标J(xi,yi)为:
其中,c0为阵形校准参考声速。
3.根据权利要求2所述的基于时延冗余测量的阵形校准方法,其特征在于,所述步骤5)具体包括:
通过优化搜索得第i个阵元Ai的坐标为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210168090.XA CN114563778B (zh) | 2022-02-23 | 2022-02-23 | 基于时延冗余测量的阵形校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210168090.XA CN114563778B (zh) | 2022-02-23 | 2022-02-23 | 基于时延冗余测量的阵形校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114563778A CN114563778A (zh) | 2022-05-31 |
CN114563778B true CN114563778B (zh) | 2023-08-29 |
Family
ID=81713113
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210168090.XA Active CN114563778B (zh) | 2022-02-23 | 2022-02-23 | 基于时延冗余测量的阵形校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114563778B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115130362B (zh) * | 2022-07-08 | 2023-04-07 | 中国科学院声学研究所 | 一种基于信道匹配的阵列阵元位置校准方法及系统 |
CN115015895B (zh) * | 2022-07-08 | 2024-08-27 | 中国科学院声学研究所 | 一种水平阵列阵元位置校准方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018063200A (ja) * | 2016-10-14 | 2018-04-19 | 日本電信電話株式会社 | 音源位置推定装置、音源位置推定方法、及びプログラム |
CN112834986A (zh) * | 2020-12-31 | 2021-05-25 | 海鹰企业集团有限责任公司 | 一种海底阵形测量方法 |
-
2022
- 2022-02-23 CN CN202210168090.XA patent/CN114563778B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018063200A (ja) * | 2016-10-14 | 2018-04-19 | 日本電信電話株式会社 | 音源位置推定装置、音源位置推定方法、及びプログラム |
CN112834986A (zh) * | 2020-12-31 | 2021-05-25 | 海鹰企业集团有限责任公司 | 一种海底阵形测量方法 |
Non-Patent Citations (1)
Title |
---|
海底布放阵形位置的精确测量方法;吕冰;葛辉良;;舰船电子工程(第09期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114563778A (zh) | 2022-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114563778B (zh) | 基于时延冗余测量的阵形校准方法 | |
EP2030041B1 (en) | Methods and systems for passive range and depth localization | |
CN108828522A (zh) | 一种利用垂直阵lcmv波束形成的水下目标辐射噪声测量方法 | |
CN109870695B (zh) | 基于深海海底反射声的非合作目标多节点水声定位方法 | |
CN109541546A (zh) | 一种基于tdoa的水下长基线声学定位方法 | |
CN113591011B (zh) | 一种快速的多波束声速改正方法及实时处理终端 | |
CN112782704B (zh) | 一种海面声源被动定位声波覆盖特性及定位精度分析方法 | |
CN103076590A (zh) | 一种基于频率预估的水声脉冲信号的定位方法 | |
Mohammadloo et al. | Correcting multibeam echosounder bathymetric measurements for errors induced by inaccurate water column sound speeds | |
CN110146924A (zh) | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 | |
CN112904316B (zh) | 一种针对海面波动的机载激光测深数据折射误差改正方法 | |
CN109975760A (zh) | 基于双矢量重构的基阵安装误差角高精度标定方法和装置 | |
CN113534161A (zh) | 一种用于远距离定位水下声源的波束镜像聚焦方法 | |
CN110471032A (zh) | 一种水下目标被动定位的方法 | |
CN116243240A (zh) | 一种基于空间多元信息融合的地面nlos环境精确定位方法 | |
RU2356060C1 (ru) | Способ измерения характеристики направленности гидроакустической антенны в условиях распространения многолучевого сигнала | |
CN111427044A (zh) | 水域目标的定位系统及定位方法 | |
Yang et al. | Acoustic inversion of bottom reflectivity and bottom sound-speed profile | |
CN116524954B (zh) | 基于反卷积理论的深海垂直阵水下目标高精度定位方法 | |
Yayu et al. | Research on location of underwater sound source target in deep sea sound field based on bellhop model | |
CN112881981B (zh) | 多种失配情况下无线传感器阵列空间处理增益损失的方法 | |
Barlee et al. | Geoacoustic model parameter estimation using a bottom-moored hydrophone array | |
RU2820807C1 (ru) | Способ определения координат шумящих объектов с использованием вертикально развитых бортовых антенн гидроакустических комплексов | |
CN112114287B (zh) | 一种方位观测数据的野值实时剔除方法 | |
CN117630894B (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 |