CN103513238A - 一种规整化最小二乘子空间相交的目标方位测向方法 - Google Patents
一种规整化最小二乘子空间相交的目标方位测向方法 Download PDFInfo
- Publication number
- CN103513238A CN103513238A CN201210202013.8A CN201210202013A CN103513238A CN 103513238 A CN103513238 A CN 103513238A CN 201210202013 A CN201210202013 A CN 201210202013A CN 103513238 A CN103513238 A CN 103513238A
- Authority
- CN
- China
- Prior art keywords
- regularization
- square
- subspace
- eigenwert
- array
- 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
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
- G01S3/803—Systems for determining direction or deviation from predetermined direction using amplitude comparison of signals derived from receiving transducers or transducer systems having differently-oriented directivity characteristics
- G01S3/8034—Systems for determining direction or deviation from predetermined direction using amplitude comparison of signals derived from receiving transducers or transducer systems having differently-oriented directivity characteristics wherein the signals are derived simultaneously
- G01S3/8036—Systems for determining direction or deviation from predetermined direction using amplitude comparison of signals derived from receiving transducers or transducer systems having differently-oriented directivity characteristics wherein the signals are derived simultaneously derived directly from separate directional systems
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
本发明涉及一种规整化最小二乘子空间相交的目标方位测向方法,该方法将吉洪诺夫规整化方法与最小二乘子空间相交的算法结合起来,处理浅海声纳装置的接收信号,在利用子空间相交的原理构造成最小二乘问题的基础上,再通过规整化方法确定最优的规整化因子,最终,实现对浅海目标的检测和精确定位。本发明将Tikhonov规整化方法与最小二乘子空间相交的算法结合起来,形成新的测向算法,该方法比常规波束形成算法具有更高的方位分辨力和准确度,并且数值计算稳定。
Description
技术领域
本发明属于声纳数字信号处理领域,特别涉及一种适用于浅海环境中的规整化最小二乘子空间相交测向方法。
背景技术
在浅海环境中,声传播受海洋波导界面的影响显著,特别是海底的影响,在这种情况下,传统的以平面波模型为基础的目标方位估计算法已不再适用,简正波模型可以更准确地描述声场,根据简正波理论,声源在海洋波导中会激发出若干号简正波,阵列的各阵元接收到的信号是各号简正波迭加的结果,各号简正波的相速度是不同的,特别是对于低频信号,相速度的差别尤为明显,这种差别在时域上反映为多途效应,采用平面波假设必然会带来目标方位估计偏差。匹配场处理(MFP)测向方法对声场用简正波模型表示,消除了测向算法中多途效应的影响,与常规平面波波束形成法相比,有明显的优点。
MFP测向的思想出现较早,因巨大的建模和信号处理计算量,未能在工程中使用。早在上世纪90年代,为了减小浅海低频测向的误差,就有人研究了相干特性对大尺度阵列的匹配场测向的性能影响,如文献1“P.S.Naidu,“On subspace method forsource localization,”JASA,90(5),1991:2489-2491”;后来子空间的思想被用于匹配测向,如文献2“J.V.Candy,“Ocean acoustic signal processing:a model-basedapproach”,JASA,1992”;同一频率的声波,各模态相速度不同,参考声速的选取不当导致阵列在端射附近出现明显的测向误差,文献3“宫在晓,林京,郭良浩,“浅海声传播相速度对测向精度的影响”,声学学报,27(6),2002:492-496”提出用声场匹配的方法做测向,同时指出匹配场测向方法需要知道海洋环境参数,而且需要对三维空间进行极大值搜索,计算量很大,计算量和环境失配是要面临的问题。
为解决MFP测向计算量大的问题,文献4“S.Lakshmipathi,and G.V.Anand,“Subspace intersection method of high-resolution bearing estimation in shallow ocean”,Signal Processing,84,2004:1367-1384”提出了一种子空间相交(Subspace intersection,SI)的算法,在已知各号简正波的波数的情况下,只需要对角度进行搜索,可以以较小的计算量实现匹配测向,较好地解决了建模计算量的问题。SI算法使用QR分解,在要分解的矩阵条件数很大时,数值计算不稳定;文献5“张爱民,基于子空间的目标方位估计的若干问题研究,中科院声学所博士学位论文”将子空间相交问题描述为最小二乘问题,提出了最小二乘SI的算法,并提出了总体最小二乘的算法,使算法在存在系数的扰动时稳健性能提高。
然而,正如所有的最小二乘问题一样,上述最小二乘SI算法中的方程的求解面临着病态问题,通常情况下,不同号数简正波的波数相差很小,在有限阵长条件下,很难满足在每个扫描向量线性无关的条件,矩阵的条件数很大,即该最小二乘问题的方程是病态的,数值计算的稳定性很差。
总体而言,最小二乘SI算法比传统波束形成算法具有无可比拟的性能优势,但数值计算稳定性差,测向精度低,可行性不高。实际应用中需要一种能够在模型参数扰动的情形下依然可用的处理方法。
发明内容
本发明目的在于,为克服现有最小二乘SI测向方法的病态问题,提出了一种规整化最小二乘子空间相交测向方法,从而使在浅海环境中的测向精度和数值计算稳定度大大提高。
为实现上述发明目的,本发明的规整化最小二乘子空间相交的测向方法,该方法将吉洪诺夫(Tikhonov)规整化方法与最小二乘子空间相交的算法结合起来,形成新的测向算法,其中Tikhonov规整化方法在文献6“Chne.H.Golub,et al.Tikhnovregularization and totle least square”中有详细论述。该方法比常规波束形成算法具有更高的方位分辨力和准确度,并且数值计算稳定。
本发明的一种规整化最小二乘子空间相交的目标方位测向方法,该方法将吉洪诺夫规整化方法与最小二乘子空间相交的算法结合起来,处理浅海声纳装置的接收信号,在利用子空间相交的原理构造成最小二乘问题的基础上,再通过规整化方法确定最优的规整化因子,最终,实现对浅海目标的检测和精确定位。
所述的浅海声纳装置由多个水听器组成,是一条海底直线阵或是拖曳阵。
该方法的步骤具体包括:
1)所述的浅海声纳装置由多个水听器组成,数目为N,阵元间距为d,目标入射方向θ,快拍长度为L;用线阵接收空间信号,得到N个阵元的时域信号x(t);
2)对L个时域快拍的数据进行快速傅里叶变换:
3)根据傅里叶变换后的快拍数据估计阵列协方差矩阵R:
R=E[XXH]
上式中,X是阵列接收数据矩阵,其中元素为步骤2)中的Xk,m(fi),H表示共轭转置,E表示对不同的快拍求统计平均;
4)对阵列协方差矩阵R进行特征值分解:
R=U∑U
式中,U为特征矢量矩阵;
其中,由特征值组成的对角阵Σ为:
如果有J个声源,则上式中的特征值满足如下关系:
λ1≥λ2≥…λJ≥λJ+1=…λN=σ2
上式中,λ表示阵列协方差的特征值,σ表示噪声的标准差;
由此,定义两个对角阵:
其中,ΣS为大特征值组成的对角阵,ΣN为小特征值组成的对角阵;
再将特征矢量矩阵分为与特征值对应的两部分:
大特征值对应的信号子空间:
US=[u1,u2,…,uJ]
和小特征值对应的噪声子空间:
UN=[uJ+1,uJ+2,…,uN]
5)根据海洋环境建模得到简正波响应矩阵:
A(θ)=[a(θ,k1),...,a(θ,kM)]∈CN×M
上式中,A表示简正波响应矩阵,ki表示第i号简正波的波数,M表示简正波号数,a(θ,ki)表示当目标方位为θ时,接收阵列相对于第i号简正波的导向矢量;
6)利用子空间相交的原理,若某一方位存在目标,利用信号子空间和简正波响应矩阵构造成最小二乘问题:
A(θ)x=Usy
上式中,A表示简正波响应矩阵,US表示大特征值对应的特征向量,x,y代表最小二乘问题的未知向量;
7)通过规整化方法,遍历所有可行的规整化因子,将规整化解向量的范数与对应的余差作出问题的L曲线,确定最优的规整化因子λ;
8)最后,得到问题的规整化最小二乘解:
x=(ATA+λI)-1ATUsy
继而得到目标方位的判别量:
上式中,P(θ)表示θ方位的判别量,将θ从1~360遍历,可得到全方位的判别量,其中最大值对应的方位即为目标方位。
本发明的优点在于:本发明将Tikhonov规整化方法与最小二乘子空间相交的算法结合起来,形成新的测向算法,该方法比常规波束形成算法具有更高的方位分辨力和准确度,并且数值计算稳定。
附图说明
图1是本发明的算法流程图;
图2是本发明实施例中的解向量示意图;
图3是本发明实施例中的海试数据处理得到的L曲线示意图;
图4是采用本发明、未采用本发明处理得到噪声目标波束输出对比示意图;
图5是采用本发明、未采用本发明处理得到宽带声源目标波束输出对比示意图。
具体实施方式
下面结合附图对本发明进行进一步说明。
本发明的基本构思是:实际应用系统中,所有的最小二乘问题的求解都面临着病态问题,最小二乘SI算法也不例外。通常情况下,不同号数简正波的波数相差很小,在有限阵长条件下,很难满足在每个扫描向量线性无关的条件,矩阵的条件数很大,即该最小二乘问题的方程是病态的,数值计算的稳定性很差。
为实现精确测向的目的,设实际阵元数目为N,阵元间距d,目标入射方向θ,阵元接收信号表示为x(t);快拍长度为L。本发明提出一种规整化最小二乘子空间相交的测向方法包含如下步骤:
1)用线阵接收空间信号,得到N个阵元的时域信号;
2)对L个时域快拍的数据做快速傅里叶变换;
3)使用傅里叶变换后的快拍数据估计阵列协方差矩阵R
R=E[XXH]
4)对R进行特征值分解如下
R=U∑U
式中,U为特征矢量矩阵,其中由特征值组成的对角阵Σ如下:
上式中的特征值满足如下关系:
λ1≥λ2≥…λJ≥λJ+1=…λN=σ2
定义如下两个对角阵:
前者为大特征值组成的对角阵,后者为小特征值组成的对角阵。
将特征矢量矩阵分为与特征值对应的两部分:一是与大特征值对应的信号子空间
US=[u1,u2,…,uJ]
二是小特征值对应的噪声子空间
UN=[uJ+1,uJ+2,…,uN]
5)根据海洋环境建模得到简正波响应矩阵
A(θ)=[a(θ,k1),...,a(θ,kM)]∈CN×M
6)利用子空间相交的原理,若某一方位存在目标,则信号子空间式及简正波响应矩阵式可以构造成最小二乘问题
A(θ)x=Usy
7)作出问题的L曲线,确定最优的规整化因子λ;
L曲线是遍历所有可行的规整化因子,将规整化解向量的范数与对应的余差画成一条曲线,将该曲线画在对数坐标轴上,它往往显示为L型,这是它名字的由来。是规整化方法中的常规手段。
8)得到问题的规整化最小二乘解:
x=(ATA+λI)-1ATUsy
继而得到目标方位的判别量:
上述步骤5)中的建模操作,在文献“M.Porter,A numerical method for computingocean acoustic normal modes.PHD 1984”中有详细的阐述。对本领域的技术人员来说,理解和实现是没问题的。
下面结合某次海试数据和附图对本发明的具体实施方式做进一步的详细描述。
试验参数:海深30m,46个阵元,2m、4m不等间隔水平直线阵,目标位于声阵的20°方位左右,信号采样率fs=2048Hz。试验中采用了两种声源信号,分别是试验船的辐射噪声信号和由信号发生器发出频率700-800Hz的宽带信号,其中,试验船辐射噪声的处理频带为150-300Hz,声速1512m/s,快拍长度60。
本发明的规整化最小二乘子空间相交测向方法,如图1所示,包括以下步骤:
步骤1:对应图1中的101,用线阵接收空间信号,得到46个阵元的时域信号;
步骤2:对应图1中的102和103,对60个时域快拍的数据做快速傅里叶变换,如下式所示,行表示时间采样,列表示阵元。
步骤3:对应图1中的104,使用傅里叶变换后的快拍数据估计阵列协方差矩阵R
R=E[XXH]
步骤4:对应图1中的105,对R进行特征值分解如下
R=U∑U
大特征值对应的空间为信号子空间,在本例中,只有1个声源,故J=1,得到
US=[u1]
步骤5:对应图1中的106,根据海洋环境建模得到简正波响应矩阵,根据本例中的海洋环境,使用Kraken建模工具,得到波导所能承载的简正波,共有10号简正波,得到阵列的简正波响应矩阵:
A(θ)=[a(θ,k1),...,a(θ,k10)]∈C46×10
步骤6:对应图1中的107,利用信号子空间式及简正波响应矩阵式构造最小二乘问题,在本例中,方程简化为:
A(θ)x=u1
图2是本方程采用原算法的解向量示意图,可以看到,该方程呈现明显的病态,系数矩阵的大小特征值相差很大(图中的三角),约差5个数量级;而系数并不随特征值减小(图中的星),这样造成的结果是在最终解当中小特征值所对应的分量占主导地位(图中的圈),所以该方程的解受噪声和计算取舍误差的影响很大。
步骤7:对应图1中的108,作出问题的L曲线,如图3所示,根据此图确定最优的规整化因子λ=0.31168;
步骤8:对应图1中的109,得到问题最终的的规整化最小二乘解:
x=(ATA+λI)-1ATu1
继而得到目标方位的判别量:
图4是应用本发明前后的船噪声目标的测向结果比较,菱形连线是原方法的处理结果,三角连线是本发明的处理结果,可以看到规整化之后波束主瓣宽度明显改善,平台效应消除,另外旁瓣也较规整化之前低,表明本发明对于原方法是由明显优势的。
图5是700-800Hz的宽带声源数据的处理结果,菱形连线是原方法的处理结果,三角连线是本发明的处理结果。同样,应用本方法后测向性能得到改善。
总之,本发明能够准确稳定地得到浅海中目标方位估计值。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种规整化最小二乘子空间相交的目标方位测向方法,该方法将吉洪诺夫规整化方法与最小二乘子空间相交的算法结合起来,处理浅海声纳装置的接收信号,在利用子空间相交的原理构造成最小二乘问题的基础上,再通过规整化方法确定最优的规整化因子,最终,实现对浅海目标的检测和精确定位。
2.根据权利要求1所述的规整化最小二乘子空间相交的目标方位测向方法,其特征在于,所述的浅海声纳装置由多个水听器组成,是一条海底直线阵或是拖曳阵。
3.根据权利要求1或2所述的规整化最小二乘子空间相交的目标方位测向方法,其特征在于,该方法的步骤包括:
1)所述的浅海声纳装置由多个水听器组成,数目为N,阵元间距为d,目标入射方向θ,快拍长度为L;用线阵接收空间信号,得到N个阵元的时域信号x(t);
2)对L个时域快拍的数据进行快速傅里叶变换:
3)根据傅里叶变换后的快拍数据估计阵列协方差矩阵R:
R=E[XXH]
上式中,X是阵列接收数据矩阵,其中元素为步骤2)中的Xk,m(fi),H表示共轭转置,E表示对不同的快拍求统计平均;
4)对阵列协方差矩阵R进行特征值分解:
R=U∑U
式中,U为特征矢量矩阵;
其中,由特征值组成的对角阵Σ为:
如果有J个声源,则上式中的特征值满足如下关系:
λ1≥λ2≥…λJ≥λJ+1=…λN=σ2
上式中,λ表示阵列协方差的特征值,σ表示噪声的标准差;
由此,定义两个对角阵:
其中,ΣS为大特征值组成的对角阵,ΣN为小特征值组成的对角阵;
再将特征矢量矩阵分为与特征值对应的两部分:
大特征值对应的信号子空间:
US=[u1,u2,…,uJ]
和小特征值对应的噪声子空间:
UN=[uJ+1,uJ+2,…,uN]
5)根据海洋环境建模得到简正波响应矩阵:
A(θ)=[a(θ,k1),...,a(θ,kM)]∈CN×M
上式中,A表示简正波响应矩阵,ki表示第i号简正波的波数,M表示简正波号数,a(θ,ki)表示当目标方位为θ时,接收阵列相对于第i号简正波的导向矢量;
6)利用子空间相交的原理,若某一方位存在目标,利用信号子空间和简正波响应矩阵构造成最小二乘问题:
A(θ)x=Usy
上式中,A表示简正波响应矩阵,US表示大特征值对应的特征向量,x,y代表最小二乘问题的未知向量;
7)通过规整化方法,遍历所有可行的规整化因子,将规整化解向量的范数与对应的余差作出问题的L曲线,确定最优的规整化因子λ;
8)最后,得到问题的规整化最小二乘解:
x=(ATA+λI)-1ATUsy
继而得到目标方位的判别量:
上式中,P(θ)表示θ方位的判别量,将θ从1~360遍历,可得到全方位的判别量,其中最大值对应的方位即为目标方位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210202013.8A CN103513238B (zh) | 2012-06-15 | 2012-06-15 | 一种规整化最小二乘子空间相交的目标方位测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210202013.8A CN103513238B (zh) | 2012-06-15 | 2012-06-15 | 一种规整化最小二乘子空间相交的目标方位测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103513238A true CN103513238A (zh) | 2014-01-15 |
CN103513238B CN103513238B (zh) | 2015-11-18 |
Family
ID=49896256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210202013.8A Expired - Fee Related CN103513238B (zh) | 2012-06-15 | 2012-06-15 | 一种规整化最小二乘子空间相交的目标方位测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103513238B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104679976A (zh) * | 2014-10-31 | 2015-06-03 | 哈尔滨工业大学深圳研究生院 | 用于信号处理的收缩线性和收缩广义线性复最小二乘算法 |
CN105954709A (zh) * | 2016-05-06 | 2016-09-21 | 哈尔滨工程大学 | 一种基于特征值多阈值修正的声矢量圆阵信源数检测方法 |
CN108226851A (zh) * | 2016-12-21 | 2018-06-29 | 罗德施瓦兹两合股份有限公司 | 用于测向的方法和测向器 |
CN108549069A (zh) * | 2018-03-08 | 2018-09-18 | 哈尔滨工程大学 | 一种非理想波导下的warping变换修正方法 |
CN108802683A (zh) * | 2018-05-30 | 2018-11-13 | 东南大学 | 一种基于稀疏贝叶斯学习的源定位方法 |
CN108845325A (zh) * | 2018-05-25 | 2018-11-20 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 拖曳线列阵声纳子阵误差失配估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1800874A (zh) * | 2004-12-31 | 2006-07-12 | 中国科学院声学研究所 | 一种用声纳测量目标的信号处理方法 |
CN101793965A (zh) * | 2010-02-08 | 2010-08-04 | 中国石化集团胜利石油管理局地球物理勘探开发公司 | 浅海检波器声学定位系统 |
EP2369799A1 (en) * | 2010-03-24 | 2011-09-28 | Pace Plc | Receiver with high performence channel estimation |
CN102695970A (zh) * | 2009-06-04 | 2012-09-26 | 道达尔公司 | 用于表征石油或者天然气储集层随时间演变的改进方法 |
-
2012
- 2012-06-15 CN CN201210202013.8A patent/CN103513238B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1800874A (zh) * | 2004-12-31 | 2006-07-12 | 中国科学院声学研究所 | 一种用声纳测量目标的信号处理方法 |
CN102695970A (zh) * | 2009-06-04 | 2012-09-26 | 道达尔公司 | 用于表征石油或者天然气储集层随时间演变的改进方法 |
CN101793965A (zh) * | 2010-02-08 | 2010-08-04 | 中国石化集团胜利石油管理局地球物理勘探开发公司 | 浅海检波器声学定位系统 |
EP2369799A1 (en) * | 2010-03-24 | 2011-09-28 | Pace Plc | Receiver with high performence channel estimation |
Non-Patent Citations (2)
Title |
---|
宋雷鸣: "基于相控麦克风阵列的逆向噪声源识别原理与技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
张爱民等: "最小二乘子空间相交方法用于浅海目标方位估计", 《哈尔滨工程大学学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104679976A (zh) * | 2014-10-31 | 2015-06-03 | 哈尔滨工业大学深圳研究生院 | 用于信号处理的收缩线性和收缩广义线性复最小二乘算法 |
CN104679976B (zh) * | 2014-10-31 | 2017-11-28 | 哈尔滨工业大学深圳研究生院 | 用于信号处理的收缩线性和收缩广义线性复最小二乘算法 |
CN105954709A (zh) * | 2016-05-06 | 2016-09-21 | 哈尔滨工程大学 | 一种基于特征值多阈值修正的声矢量圆阵信源数检测方法 |
CN105954709B (zh) * | 2016-05-06 | 2018-12-07 | 哈尔滨工程大学 | 一种基于特征值多阈值修正的声矢量圆阵信源数检测方法 |
CN108226851A (zh) * | 2016-12-21 | 2018-06-29 | 罗德施瓦兹两合股份有限公司 | 用于测向的方法和测向器 |
CN108226851B (zh) * | 2016-12-21 | 2023-03-24 | 罗德施瓦兹两合股份有限公司 | 用于测向的方法和测向器 |
CN108549069A (zh) * | 2018-03-08 | 2018-09-18 | 哈尔滨工程大学 | 一种非理想波导下的warping变换修正方法 |
CN108845325A (zh) * | 2018-05-25 | 2018-11-20 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 拖曳线列阵声纳子阵误差失配估计方法 |
CN108845325B (zh) * | 2018-05-25 | 2022-07-05 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 拖曳线列阵声纳子阵误差失配估计方法 |
CN108802683A (zh) * | 2018-05-30 | 2018-11-13 | 东南大学 | 一种基于稀疏贝叶斯学习的源定位方法 |
CN108802683B (zh) * | 2018-05-30 | 2021-04-27 | 东南大学 | 一种基于稀疏贝叶斯学习的源定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103513238B (zh) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103513238B (zh) | 一种规整化最小二乘子空间相交的目标方位测向方法 | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN101813772B (zh) | 一种快速宽带频域扩展拖曳阵波束形成方法 | |
CN103513250B (zh) | 一种基于鲁棒自适应波束形成原理的模基定位方法及系统 | |
CN105301580B (zh) | 一种基于分裂阵互谱相位差方差加权的被动探测方法 | |
CN104777453A (zh) | 舰船线谱噪声源定位的波束域时频分析方法 | |
CN101581785B (zh) | 用于小型水下载体的三维前视声像声纳系统和方法 | |
CN103267966A (zh) | 一种单个二维压差式矢量水听器方位估计的方法 | |
CN112269164A (zh) | 深海可靠声路径下基于干涉结构匹配处理弱目标定位方法 | |
CN101915922A (zh) | 拖曳线列阵被动测距方法 | |
CN104730513A (zh) | 一种分级子阵聚焦mvdr波束形成方法 | |
CN102692620A (zh) | 适用于浅海多途环境的水下噪声源高稳健性聚焦定位方法 | |
CN109581388A (zh) | 一种实时三维成像声纳的近场宽视角波束形成方法 | |
CN104656073A (zh) | 三维成像声纳波束形成方法及在多核处理器上的实现方法 | |
CN110687538A (zh) | 一种基于近场聚焦的超波束形成方法 | |
CN103454616A (zh) | 一种十字型振速梯度水听器的方位估计方法 | |
CN112987004A (zh) | 一种浅海环境下基于水平阵列的水面水下目标分类方法 | |
CN104502904A (zh) | 一种鱼雷自导波束锐化方法 | |
CN103323832A (zh) | 一种相控阵三维摄像声纳系统换能器阵列的幅相误差校正方法 | |
CN112612027B (zh) | 一种浅海环境下利用声能量起伏的海洋内波监测方法 | |
CN103048641A (zh) | 具有联合约束优化形式的矢量阵稳健聚焦处理方法 | |
CN109597034B (zh) | 一种基于欧几里得距离的空时自适应处理方法 | |
CN115201821B (zh) | 基于强目标成像对消的小目标检测方法 | |
CN103513249A (zh) | 一种宽带相干模基信号处理方法及系统 | |
CN102183755A (zh) | 一种基于柯西-高斯模型的新型高分辨方位估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151118 Termination date: 20180615 |
|
CF01 | Termination of patent right due to non-payment of annual fee |