CN111679248B - 一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 - Google Patents
一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 Download PDFInfo
- Publication number
- CN111679248B CN111679248B CN202010415286.5A CN202010415286A CN111679248B CN 111679248 B CN111679248 B CN 111679248B CN 202010415286 A CN202010415286 A CN 202010415286A CN 111679248 B CN111679248 B CN 111679248B
- Authority
- CN
- China
- Prior art keywords
- target
- distance
- sparse
- azimuth
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000011159 matrix material Substances 0.000 claims abstract description 54
- 238000001228 spectrum Methods 0.000 claims abstract description 13
- 238000005457 optimization Methods 0.000 claims abstract description 10
- 230000005236 sound signal Effects 0.000 claims abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000013535 sea water Substances 0.000 claims description 4
- 235000007164 Oryza sativa Nutrition 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 235000009566 rice Nutrition 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims 1
- 238000007796 conventional method Methods 0.000 description 13
- 230000001427 coherent effect Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- WABPQHHGFIMREM-UHFFFAOYSA-N lead(0) Chemical compound [Pb] WABPQHHGFIMREM-UHFFFAOYSA-N 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 235000020046 sherry Nutrition 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
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
- 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/20—Position of source determined by a plurality of spaced direction-finders
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
本发明是一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法。所述方法在海底布放水平L型阵列,接收临近海面附近多目标发出的声信号,获得多个目标的源矩阵;在目标的方位和距离二维平面中,获取关于目标的稀疏源矩阵;获得关于N元L型阵列接收信号的稀疏阵列矩阵;利用得到的稀疏阵列矩阵和归一化后的稀疏源矩阵,重构稀疏信号;采用范数正则法对重构的稀疏信号进行最优化求解;通过空间谱的谱峰位置和强度,确定目标的方位和距离信息,并判断目标功率的相对大小。本发明可以有效剥离减少多途干扰对目标方位‑距离联合估计精度的影响,获得无空间模糊的目标定位结果,并可同时准确获得多目标的位置信息和强度相对大小。
Description
技术领域
本发明涉及海底信号定位检测技术领域,是一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法。
背景技术
在浅海环境下,声信道包括海面、海底,声传播存在明显的多途效应,多途径声信号对于浅海目标方位估计和定位等问题均是最为重要的干扰之一。对于浅海中近程目标的声学定位问题而言,多途效应的存在会使定位谱图中结果出现大量虚假的“伪峰”,无法准确判断真实目标的位置,同时“伪峰”的存在还可能影响到对真实“弱目标”所在位置的判断,导致多目标分辨结果出现严重错误。
此外,常规CBF和MVDR等信号处理器在对同频相干目标的分辨上能力不足,常规CBF处理器可适用于相干目标,但空间分辨率不高,无法对较近空间距离的目标进行分辨,而MVDR处理器虽然可提高分辨率,但无法直接适用于相干声源(罗方方,生雪莉,梅继丹,郭咏.基于MVDR高分辨算法的时反定位技术研究.哈尔滨工程大学学报,2010,31(7):945-950),使得以上两种方法在实际适用中受到较大限制。
专利ZL201210162659.8适用于浅海多途环境的水下噪声源高稳健性聚焦定位方法,针对浅海环境下的信号处理问题,基于模型匹配的声纳被动定位思想,对多途信号获得了一定的抑制效果,但由于采用的是线列阵,其缺点是无法保证对目标方位上全空间的“无模糊”定位。在浅海噪声源聚焦定位方面,对浅海多途效应的影响以及可行性方法的研究普遍集中于时间反转镜方法(惠娟,胡丹,惠俊英等.聚焦波束形成声图测量原理研究.声学学报,2007,32(4):356-361),仍未出现利用稀疏信号重构以获得高精度方位-距离联合估计的定位方法。
发明内容
本发明为有效的抑制“伪峰”的出现,便于对“弱目标”的分辨,保证了多目标定位和分辨的准确性和有效性,本发明提供了一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,本发明提供了以下技术方案:
一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,所述方法包括以下步骤:
步骤1:在海底布放水平L型阵列,接收临近海面附近多目标发出的声信号,获得多个目标的源矩阵;
步骤2:在目标的方位和距离二维平面中,获取关于目标的稀疏源矩阵;
步骤3:获得关于N元L型阵列接收信号的稀疏阵列矩阵;
步骤4:利用得到的稀疏阵列矩阵和归一化后的稀疏源矩阵,重构稀疏信号;
步骤5:采用范数正则法对重构的稀疏信号进行最优化求解;
步骤6:获得在目标的方位和距离二维平面上包含目标所在方位和距离信息的空间稀疏信号的波形估计结果,对波形估计结果进行的功率计算,获得方位和距离联合空间谱估计结果;
步骤7:通过空间谱的谱峰位置和强度,确定目标的方位和距离信息,并判断目标功率的相对大小。
优选地,所述步骤1具体为:
步骤1.1:将N元L型阵列的参考阵元位于坐标系原点,第i个阵元的坐标为(xi,yi),i=1,2,…,N,采用浅海声场的Pekeries模型,海水深度为H米,阵列置于水下hr米处;
步骤1.2:目标的方位和距离平面S中共存在Q个目标,第q个目标的深度为hsq,q=1,2,…,Q,且与L型阵列参考阵元的水平距离为rsq,方位角为θsq,各目标对基阵接收信号的虚源个数为M,第q个目标的第p个虚源的方位角和俯仰角分别为θpq和p=0,1,…,M,且θpq=θsq;
步骤1.3:确定第q个声源对应的源矢量Bq,通过下式表示Bq:
其中,αpq和Rpq分别表示第q个声源的第p个虚源的反射系数和到参考阵元的声程距离;Aq为第q个源的镜像空间矩阵,Wq为第q个源的幅度加权系数矢量,φpqi表示第q个声源的第p个虚源到第i个阵元与参考阵元间的相位差;
得到Q个目标的源矩阵B,B=[B1 B2…Bq…BQ]。
优选地,所述步骤2具体为:
目标的方位和距离二维平面S中,目标水平方位扫描范围为[0 2π],方位扫描步长为θstep,方位扫描区域个数为Kθ,目标距离参考阵元的水平距离扫描范围为[0 rmax],距离扫描步长为rstep,水平距离扫描区域个数为Kr,将声源所在方位和距离平面划分为{Θ1,Θ2,...,Θk,...,ΘK}个区域,其中,k对应于第kθi个方位和第kri个距离的区域序号,kθi=1,2,…,Kθ,kri=1,2,…,Kr,通过下式表示k:
k=(kri-1)Kr+krθ (4)
确定对应于第kθi方位和第kri个距离所在区域的源矢量Ak,通过下式表示Ak:
其中,αpk和Rpk分别表示第k个存在声源的空间区域点对应的第p个虚源的反射系数和到参考阵元的声程距离,φpki表示第k个存在声源的空间点对应的第p个虚源到第i个阵元与参考阵元间的相位差;
得到目标的方位和距离二维平面S中的稀疏源矩阵A,通过下式表示A:
A=[A1 A2…Ak…AK] (6)
对目标的方位和距离二维平面S中的稀疏源矩阵A进行2范数下的归一化处理,得到归一化后的稀疏源矩阵,通过下式表示归一化后的稀疏源矩阵:
优选地,所述步骤3具体为:
获得关于N元L型阵列接收信号的稀疏阵列矩阵,通过下式表示N元L型阵列接收信号的稀疏阵列矩阵:
y=As+n (8)
其中,y为N×Ls维阵列接收信号矩阵稀疏阵列矩阵,s表示包含声源位置信息的空间稀疏信号,Ls为信号长度,n为阵列接收到的噪声信号。
优选地,步骤4具体为:从L型阵列接收数据中恢复声源信号波形,将稀疏信号处理最优化求解,通过下式表示重构后的稀疏信号:
min||s||l s.t.||y-As||2≤δ (9)
其中,δ为噪声约束参数,s.t.表示约束条件,下标l表示某范数。
优选地,采用l1范数正则法对重构稀疏信号进行最优化求解,通过下式表示最优化解:
min||s||1 s.t.||y-As||2≤δ (10)
其中,||.||1表示l1范数,||.||2表示l2范数,s.t.为满足右侧式子约束条件下使左侧式子的l1范数最小。
优选地,通过线性规划理论进行求解,包括使用l1-MAGIC数据包,SeDuMi软件以及CVX工具箱进行求解。
本发明具有以下有益效果:
本发明构造了新颖的浅海海底L型阵列中近程目标方位和距离联合的声源信号稀疏表示形式;利用L阵型和信道多途匹配的源矩阵构造方法,保证了目标方位和距离估计结果在全空间上“无模糊”;利用信道多途特性,针对不同的方位和距离扫描区域生成相应的源矩阵,保证了对目标信息的有效利用,并尽可能保证对多途干扰的剥离,定位结果也显示该方法可有效抑制多途效应带来的“伪峰”问题;利用稀疏重构思想,解决了强相干多目标的方位和距离联合估计问题,可在同时获得目标定位结果的同时,获得目标强度的相对大小。该方法保留了对强弱目标的分辨能力,可在获得“强目标”位置的同时,兼顾对“弱目标”的定位,定位结果更为有效和可靠。
本发明提出的一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,克服了传统方法抗多途影响能力差、空间模糊、相干目标不适用等问题,可在浅海环境下对相干多目标进行准确的方位和距离估计,获得强弱目标准确的定位结果和强度大小分辨结果,并保证了目标方位和距离估计结果在全空间上“无模糊”,保证了多目标定位和分辨的准确性和有效性。该方法可以满足在浅海多途条件下,海底监测用阵列对水面和水中多目标定位的迫切需求。
本发明涉及一种新的浅海中目标方位-距离联合估计方法,利用坐底L型阵列实现对水面船只、水中UUV等目标位置的定位。海底水平L型阵列可布放于水下油气管道、通信光纤等重要设备附近,对设备上方作业船只等进行声学监测和预警防护。
该方法充分利用多目标在方位-距离二维空间上的稀疏性,通过构造与L型阵列和浅海信道匹配的源矩阵稀疏表示,可以有效剥离减少多途干扰对目标方位-距离联合估计精度的影响,获得无空间模糊的目标定位结果,并可同时准确获得多目标的位置信息和强度相对大小。
附图说明
图1为基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法流程图;
图2为浅海L型阵列接收声源信号和声线结构示意图;
图3为本专利方法和传统方法对比图,Kθ=37,Kr=50;图3-(a)为本专利,图3-(b)为传统方法;
图4为本专利方法和传统方法对比图,Kθ=37,Kr=99;图4-(a)为本专利,图4-(b)为传统方法。
具体实施方式
以下结合具体实施例,对本发明进行了详细说明。
具体实施例一:
根据图1所示,本发明提供一种基于电机运行状态的磁悬浮轴承控制***,包括以下步骤:
步骤1:在海底布放水平L型阵列,接收临近海面附近多目标发出的声信号,获得Q个目标的源矩阵B。
图2所示为浅海海底布放的L型阵列接收目标声信号及其声线的示意图。将N元L型阵列的参考阵元位于坐标系原点,第i(i=1,2,…,N)个阵元的坐标为(xi,yi)。考虑浅海声场的Pekeries模型,海水深度为H米,阵列置于水下hr米处,靠近海底。
考虑在方位和距离平面S中共存在Q个目标,均位于水面附近,第q(q=1,2,…,Q)个目标的深度为hsq,且与L型阵列参考阵元的水平距离为rsq,方位角为θsq。假设各目标对基阵接收信号有意义的虚源个数为M,第q(q=1,2,…,Q)个目标的第p(p=0,1,…,M)个虚源的方位角和俯仰角分别为θpq,有θpq=θsq。
得到第q个声源对应的源矢量Bq为:
得到Q个目标的源矩阵B可表示为B=[B1 B2…Bq…BQ]。
B为考虑全部声源的源矩阵,由各个声源考虑了多途信道影响下的方位和距离联合导向矢量构成,Bq为第q个源的源矢量,Aq的列张成了第q个源的镜像空间矩阵,Wq为第q个源的幅度加权系数矢量。
步骤2:在目标方位和距离二维平面S中,获得关于目标的稀疏源矩阵A。
设置目标的方位和距离二维平面S中,目标水平方位扫描范围为[0 2π],方位扫描步长为θstep,方位扫描区域个数为Kθ。目标距离参考阵元的水平距离扫描范围为[0 rmax],距离扫描步长为rstep,水平距离扫描区域个数为Kr。则声源所在方位和距离平面S可划分为{Θ1,Θ2,...,Θk,...,ΘK}个区域,其中,k对应于第kθi(kθi=1,2,…,Kθ)个方位和第kri(kri=1,2,…,Kr)个距离的区域序号,有
k=(kri-1)Kr+krθ K=Kθ·Kr (4)
得到对应于第kθi方位和第kri个距离所在区域的源矢量Ak,有
同样,αpk和Rpk分别表示第k个可能存在声源的空间区域点对应的第p个虚源的反射系数和到参考阵元的声程距离。φpki表示第k个可能存在声源的空间点对应的第p个虚源到第i个阵元与参考阵元间的相位差,计算方法参考步骤1,不再赘述。
得到感兴趣的方位和距离二维平面S中的稀疏源矩阵可表示为:
A=[A1 A2…Ak…AK] (6)
进行2范数下的归一化处理,得到归一化后的稀疏源矩阵
步骤3:获得关于N元L型阵列接收信号的稀疏表示。
y=As+n (8)
y为N×Ls维阵列接收信号矩阵,A为N×K维源导向矩阵,s为K×Ls维源信号矩阵,Ls为信号长度。考虑噪声影响,n为阵列接收到的噪声信号。
将方位和距离平面S划分为{Θ1,Θ2,...,Θk,...,ΘK}离散区域,因为空间中只存在Q个真实声源(Q<<K),则在S上,只有Q个空间区域存在信号,即构造的K×Ls维信号矩阵s中只有Q行非零元素的波形数据。据此(8)式的物理意义已经十分清晰,A为N×K维空间源导向矩阵,y表示阵列接收的信号波形,s表示包含声源位置信息的空间稀疏信号,y本质上就是s的稀疏表示。
步骤4:利用得到的阵列接收数据y矩阵和完备空间源导向矩阵A,重构稀疏信号s。
利用得到的阵列接收数据y矩阵和完备空间源导向矩阵A,重构稀疏信号s。就是从L型阵列接收数据中恢复声源信号波形。将稀疏信号处理过程表示为以下最优化求解:
min||s||l s.t.||y-As||2≤δ (9)
其中,δ为噪声约束参数,s.t.表示约束条件,下标l表示某范数。
步骤5:利用l1范数正则法对重构稀疏信号进行最优化求解。
min||s||1 s.t.||y-As||2≤δ (10)
其中,||.||1表示l1范数,||.||2表示l2范数,min表示取最小值。s.t.含义为满足右侧式子约束条件下使左侧式子的l1范数最小。上式为一个欠定方程,数学上是一个凸优化问题,可通过线性规划理论进行求解,例如可使用l1-MAGIC数据包,SeDuMi软件以及CVX工具箱等进行求解。
步骤6:获得平面S上{Θ1,Θ2,...,Θk,...,ΘK}范围内包含目标所在方位和距离信息的空间稀疏信号s的波形估计结果后,对波形数据进行简单的功率计算,可获得方位和距离联合空间谱估计结果。
步骤7:通过空间谱的谱峰位置和强度,同时确定Q个目标的方位和距离信息,并可判断Q个目标功率的相对大小。
具体实施例二:
本发明提出了一种基于海底水平L型阵列的目标方位-距离联合稀疏重构定位方法,该方法利用L型阵列,其阵列流程为获得对目标方位上全空间的“无模糊”定位提供了条件,通过构建新颖的与L型阵列和浅海信道匹配的稀疏表示源矩阵,充分利用多目标在空间上的稀疏性,通过稀疏重构算法获得了目标信号的重构结果,可同时准确获得多目标的方位-距离联合定位结果,并可获得多个目标强度的相对大小,较一般的定位方法仅可获得定位结果而无法获得强度大小具有较大的优越性,同时,稀疏重构算法可适用于相干信号,且不受目标数量的限制,对于解决水下强相干目标的定位问题具有较大优势。通过与L型阵列和浅海信道匹配的信源稀疏重构,可以有效剥离多途干扰对目标方位-距离联合估计精度的影响,与常规方法相比,可以有效的抑制“伪峰”的出现,便于对“弱目标”的分辨,保证了多目标定位和分辨的准确性和有效性。
实例参数设置如下:L型阵列的阵元个数为29,阵元间距为12m,浅海水深为250m,基阵布放于距离海底1m深度处。两个单频目标声源的频率均为250Hz,目标1的水平方位角为120°,距离参考阵元的水平距离为6倍水深,即1500m,目标1的水平方位角为250°,距离参考阵元的水平距离为5倍水深,即1250m,目标2声源强度是目标1声源强度的0.75倍。双目标均位于水下2m深度处。***采样率为4kHz。水中声速取为1480m/s。设置水平方位角的扫描范围为[0 2π],方位扫描步长为θstep=10π/180,方位扫描区域个数为Kθ=37。目标距离参考阵元的水平距离扫描范围为[50 2500]m,距离扫描步长为rstep=50m,水平距离扫描区域个数为Kr=50。信噪比为5dB,选取为δ=0.15。
两种方法的定位对比结果如图3所示,图中(a)表示本专利方法的方位-距离定位结果,(b)表示常规方法的方位-距离定位结果,图中的横坐标表示距离,单位m,纵坐标表示方位角,单位度。伪彩图的色棒颜色代表估计目标功率的归一化值。
对比可以看出:
(1)利用L型阵列和多途信道匹配的方式,可以保证在0度到360度方位内,方位-距离估计结果不存在空间模糊问题,这与线列阵相比具有较大优势;
(2)本专利方法和常规方法相比,常规方法虽然利用了信道的匹配性,但在整个“方位-距离”平面上存在大量的起伏“伪峰”,当弱目标存在时,基本无法准确判断“弱目标”存在,导致存在弱目标丢失;而本专利方法可以准确给出双相干目标的定位结果,并可给出目标功率的强弱相对大小,且具有非常强的背景压制能力,除了真实目标出现的问题,完全不存在虚假“伪峰”的影响,可十分准确的对强弱目标的定位和强度信息进行判断;
(3)从常规方法的谱峰宽度可以看出,常规方法在距离方向上的分辨率较差,常规方法在500m左右的距离方位内,均显示有强目标出现,说明常规方法无法给出准确的定位结果;而本专利方法的分辨率非常优秀,可以直接给出准确而谱峰尖锐的定位结果,且与仿真预设的方位-距离参数完全一致,定位精度非常高。
进一步调整仿真参数,距离扫描步长为rstep=25m,水平距离扫描区域个数为Kr=99。其它参数不变。两种方法的定位对比结果如图4(a)和图4(b)所示。可以看出:
常规方法由于扫描细致程度的改进,谱图的显示也更为精细,但在主峰处出现了大量强度接近的“伪峰”群,导致无法准确判断目标位置,同样大量的背景起伏对弱目标产生了掩盖的影响,仍然无法对弱目标进行分辨。而本专利方法在扫描精细度增大的同时,仍然保持准确的定位结果和强弱目标分辨能力,性能优异。
以上所述仅是一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法的优选实施方式,一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法的保护范围并不仅局限于上述实施例,凡属于该思路下的技术方案均属于本发明的保护范围。应当指出,对于本领域的技术人员来说,在不脱离本发明原理前提下的若干改进和变化,这些改进和变化也应视为本发明的保护范围。
Claims (6)
1.一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,其特征是:所述方法包括以下步骤:
步骤1:在海底布放水平L型阵列,接收临近海面附近多目标发出的声信号,获得多个目标的源矩阵;
所述步骤1具体为:
步骤1.1:将N元L型阵列的参考阵元位于坐标系原点,第i个阵元的坐标为(xi,yi),i=1,2,…,N,采用浅海声场的Pekeries模型,海水深度为H米,阵列置于水下hr米处;
步骤1.2:目标的方位和距离平面S中共存在Q个目标,第q个目标的深度为hsq,q=1,2,…,Q,且与L型阵列参考阵元的水平距离为rsq,方位角为θsq,各目标对基阵接收信号的虚源个数为M,第q个目标的第p个虚源的方位角和俯仰角分别为θpq和p=0,1,…,M,且θpq=θsq;
步骤1.3:确定第q个声源对应的源矢量Bq,通过下式表示Bq:
其中,αpq和Rpq分别表示第q个声源的第p个虚源的反射系数和到参考阵元的声程距离;Aq为第q个源的镜像空间矩阵,Wq为第q个源的幅度加权系数矢量,φpqi表示第q个声源的第p个虚源到第i个阵元与参考阵元间的相位差;
步骤2:在目标的方位和距离二维平面中,获取关于目标的稀疏源矩阵;
步骤3:获得关于N元L型阵列接收信号的稀疏阵列矩阵;
步骤4:利用得到的稀疏阵列矩阵和归一化后的稀疏源矩阵,重构稀疏信号;
步骤5:采用范数正则法对重构的稀疏信号进行最优化求解;
步骤6:获得在目标的方位和距离二维平面上包含目标所在方位和距离信息的空间稀疏信号的波形估计结果,对波形估计结果进行的功率计算,获得方位和距离联合空间谱估计结果;
步骤7:通过空间谱的谱峰位置和强度,确定目标的方位和距离信息,并判断目标功率的相对大小。
2.根据权利要求1所述的一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,其特征是:所述步骤2具体为:
目标的方位和距离二维平面S中,目标水平方位扫描范围为[02π],方位扫描步长为θstep,方位扫描区域个数为Kθ,目标距离参考阵元的水平距离扫描范围为[0rmax],距离扫描步长为rstep,水平距离扫描区域个数为Kr,将声源所在方位和距离平面划分为{Θ1,Θ2,...,Θk,...,ΘK}个区域,其中,k对应于第kθi个方位和第kri个距离的区域序号,kθi=1,2,…,Kθ,kri=1,2,…,Kr,通过下式表示k:
k=(kri-1)Kr+krθ (4)
确定对应于第kθi方位和第kri个距离所在区域的源矢量Ak,通过下式表示Ak:
其中,αpk和Rpk分别表示第k个存在声源的空间区域点对应的第p个虚源的反射系数和到参考阵元的声程距离,φpki表示第k个存在声源的空间点对应的第p个虚源到第i个阵元与参考阵元间的相位差;
得到目标的方位和距离二维平面S中的稀疏源矩阵A,通过下式表示A:
A=[A1 A2 … Ak … AK] (6)
对目标的方位和距离二维平面S中的稀疏源矩阵A进行2范数下的归一化处理,得到归一化后的稀疏源矩阵,通过下式表示归一化后的稀疏源矩阵:
3.根据权利要求1所述的一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,其特征是:所述步骤3具体为:
获得关于N元L型阵列接收信号的稀疏阵列矩阵,通过下式表示N元L型阵列接收信号的稀疏阵列矩阵:
y=As+n (8)
其中,y为N×Ls维阵列接收信号矩阵稀疏阵列矩阵,s表示包含声源位置信息的空间稀疏信号,Ls为信号长度,n为阵列接收到的噪声信号。
4.根据权利要求1所述的一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,其特征是:步骤4具体为:从L型阵列接收数据中恢复声源信号波形,将稀疏信号处理最优化求解,通过下式表示重构后的稀疏信号:
min||s||ls.t.||y-As||2≤δ (9)
其中,δ为噪声约束参数,s.t.表示约束条件,下标l表示某范数。
5.根据权利要求1所述的一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,其特征是:采用l1范数正则法对重构稀疏信号进行最优化求解,通过下式表示最优化解:
min||s||1s.t.||y-As||2≤δ (10)
其中,||·||1表示l1范数,||·||2表示l2范数,s.t.为满足右侧式子约束条件下使左侧式子的l1范数最小。
6.根据权利要求5所述的一种基于海底水平L型阵列的目标方位和距离联合稀疏重构定位方法,其特征是:通过线性规划理论进行求解,包括使用l1-MAGIC数据包、SeDuMi软件以及CVX工具箱进行求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010415286.5A CN111679248B (zh) | 2020-05-15 | 2020-05-15 | 一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010415286.5A CN111679248B (zh) | 2020-05-15 | 2020-05-15 | 一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111679248A CN111679248A (zh) | 2020-09-18 |
CN111679248B true CN111679248B (zh) | 2023-04-21 |
Family
ID=72451924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010415286.5A Active CN111679248B (zh) | 2020-05-15 | 2020-05-15 | 一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111679248B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112269163B (zh) * | 2020-09-30 | 2023-04-25 | 黑龙江工程学院 | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 |
CN113671493B (zh) * | 2021-08-09 | 2023-08-11 | 黑龙江工程学院 | 一种基于特征融合的海面小目标检测方法及*** |
Citations (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6424596B1 (en) * | 2001-10-01 | 2002-07-23 | The United States Of America As Represented By The Secretary Of The Navy | Method and apparatus for reducing noise from near ocean surface sources |
JP2006125986A (ja) * | 2004-10-28 | 2006-05-18 | Nippon Telegr & Teleph Corp <Ntt> | 波源探索方法および波源探索システム |
JP2007047427A (ja) * | 2005-08-10 | 2007-02-22 | Hitachi Ltd | 音声処理装置 |
WO2007116398A2 (en) * | 2006-04-10 | 2007-10-18 | Noam Amram | A method device and system for receiving a communication signal |
CN101644774A (zh) * | 2009-07-07 | 2010-02-10 | 中国科学院声学研究所 | 一种基于线性规划的近场目标定位方法及其*** |
JP2014098568A (ja) * | 2012-11-13 | 2014-05-29 | Advanced Telecommunication Research Institute International | 音源位置推定装置、音源位置推定方法および音源位置推定プログラム |
CN103983948A (zh) * | 2014-05-29 | 2014-08-13 | 西安电子科技大学 | 基于稀疏表示的目标角度距离联合估计方法 |
CN104678372A (zh) * | 2015-03-14 | 2015-06-03 | 西安电子科技大学 | 正交频分复用雷达超分辨距离与角度值联合估计方法 |
CN105372623A (zh) * | 2015-12-15 | 2016-03-02 | 吉林大学 | 一种基于l型阵列的信源仰角和方位角估计方法 |
CN105553894A (zh) * | 2016-01-19 | 2016-05-04 | 南京航空航天大学 | 声矢量阵中一种基于平行因子压缩感知的盲源分离方法 |
CN105828266A (zh) * | 2016-03-11 | 2016-08-03 | 苏州奇梦者网络科技有限公司 | 一种麦克风阵列的信号处理方法与*** |
CN106680762A (zh) * | 2016-12-15 | 2017-05-17 | 哈尔滨工程大学 | 一种基于互协方差稀疏重构的声矢量阵方位估计方法 |
CN106872936A (zh) * | 2017-03-21 | 2017-06-20 | 陕西理工大学 | 近场源l型声矢量传感器阵列解模糊多参数估计方法 |
CN107966677A (zh) * | 2017-11-16 | 2018-04-27 | 黑龙江工程学院 | 一种基于空间稀疏约束的圆阵模态域方位估计方法 |
JP2018077139A (ja) * | 2016-11-09 | 2018-05-17 | 日本電信電話株式会社 | 音場推定装置、音場推定方法、プログラム |
US10206036B1 (en) * | 2018-08-06 | 2019-02-12 | Alibaba Group Holding Limited | Method and apparatus for sound source location detection |
US10310074B1 (en) * | 2014-03-27 | 2019-06-04 | Hrl Laboratories, Llc | System and method for denoising synthetic aperture radar (SAR) images via sparse and low-rank (SLR) decomposition and using SAR images to image a complex scene |
EP3520437A1 (en) * | 2016-09-29 | 2019-08-07 | Dolby Laboratories Licensing Corporation | Method, systems and apparatus for determining audio representation(s) of one or more audio sources |
CN110196427A (zh) * | 2019-05-29 | 2019-09-03 | 西北工业大学 | 一种基于距离方位字典的目标定位算法 |
CN110244272A (zh) * | 2019-06-14 | 2019-09-17 | 西安电子科技大学 | 基于秩一去噪模型的波达方向估计方法 |
CN110927662A (zh) * | 2019-05-16 | 2020-03-27 | 陕西理工大学 | 基于声矢量传感器的近场多参数估计四阶累积量方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110317522A1 (en) * | 2010-06-28 | 2011-12-29 | Microsoft Corporation | Sound source localization based on reflections and room estimation |
WO2012077287A1 (ja) * | 2010-12-09 | 2012-06-14 | パナソニック株式会社 | 姿勢状態推定装置および姿勢状態推定方法 |
US9562968B2 (en) * | 2012-10-22 | 2017-02-07 | Saab-Sensis Corporation | Sensor system and method for determining target location using sparsity-based processing |
-
2020
- 2020-05-15 CN CN202010415286.5A patent/CN111679248B/zh active Active
Patent Citations (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6424596B1 (en) * | 2001-10-01 | 2002-07-23 | The United States Of America As Represented By The Secretary Of The Navy | Method and apparatus for reducing noise from near ocean surface sources |
JP2006125986A (ja) * | 2004-10-28 | 2006-05-18 | Nippon Telegr & Teleph Corp <Ntt> | 波源探索方法および波源探索システム |
JP2007047427A (ja) * | 2005-08-10 | 2007-02-22 | Hitachi Ltd | 音声処理装置 |
WO2007116398A2 (en) * | 2006-04-10 | 2007-10-18 | Noam Amram | A method device and system for receiving a communication signal |
CN101644774A (zh) * | 2009-07-07 | 2010-02-10 | 中国科学院声学研究所 | 一种基于线性规划的近场目标定位方法及其*** |
JP2014098568A (ja) * | 2012-11-13 | 2014-05-29 | Advanced Telecommunication Research Institute International | 音源位置推定装置、音源位置推定方法および音源位置推定プログラム |
US10310074B1 (en) * | 2014-03-27 | 2019-06-04 | Hrl Laboratories, Llc | System and method for denoising synthetic aperture radar (SAR) images via sparse and low-rank (SLR) decomposition and using SAR images to image a complex scene |
CN103983948A (zh) * | 2014-05-29 | 2014-08-13 | 西安电子科技大学 | 基于稀疏表示的目标角度距离联合估计方法 |
CN104678372A (zh) * | 2015-03-14 | 2015-06-03 | 西安电子科技大学 | 正交频分复用雷达超分辨距离与角度值联合估计方法 |
CN105372623A (zh) * | 2015-12-15 | 2016-03-02 | 吉林大学 | 一种基于l型阵列的信源仰角和方位角估计方法 |
CN105553894A (zh) * | 2016-01-19 | 2016-05-04 | 南京航空航天大学 | 声矢量阵中一种基于平行因子压缩感知的盲源分离方法 |
CN105828266A (zh) * | 2016-03-11 | 2016-08-03 | 苏州奇梦者网络科技有限公司 | 一种麦克风阵列的信号处理方法与*** |
EP3520437A1 (en) * | 2016-09-29 | 2019-08-07 | Dolby Laboratories Licensing Corporation | Method, systems and apparatus for determining audio representation(s) of one or more audio sources |
JP2018077139A (ja) * | 2016-11-09 | 2018-05-17 | 日本電信電話株式会社 | 音場推定装置、音場推定方法、プログラム |
CN106680762A (zh) * | 2016-12-15 | 2017-05-17 | 哈尔滨工程大学 | 一种基于互协方差稀疏重构的声矢量阵方位估计方法 |
CN106872936A (zh) * | 2017-03-21 | 2017-06-20 | 陕西理工大学 | 近场源l型声矢量传感器阵列解模糊多参数估计方法 |
CN107966677A (zh) * | 2017-11-16 | 2018-04-27 | 黑龙江工程学院 | 一种基于空间稀疏约束的圆阵模态域方位估计方法 |
US10206036B1 (en) * | 2018-08-06 | 2019-02-12 | Alibaba Group Holding Limited | Method and apparatus for sound source location detection |
CN110927662A (zh) * | 2019-05-16 | 2020-03-27 | 陕西理工大学 | 基于声矢量传感器的近场多参数估计四阶累积量方法 |
CN110196427A (zh) * | 2019-05-29 | 2019-09-03 | 西北工业大学 | 一种基于距离方位字典的目标定位算法 |
CN110244272A (zh) * | 2019-06-14 | 2019-09-17 | 西安电子科技大学 | 基于秩一去噪模型的波达方向估计方法 |
Non-Patent Citations (23)
Title |
---|
H. Song.Vector-Sensor Array DOA Estimation Based on Spatial Time-Frequency Distribution.2018 Eighth International Conference on Instrumentation & Measurement, Computer, Communication and Control (IMCCC).2018,1351-1356. * |
H. Song.Vector-Sensor Array DOA Estimation Based on Spatio-Temporal Correlation Matrices Joint Diagonalization Using Jacobi Rotation.2018 Eighth International Conference on Instrumentation & Measurement, Computer, Communication and Control (IMCCC).2018,1357-1362. * |
Hu, K.Near-Field Source Localization: Sparse Recovery Techniques and Grid Matching.2014 IEEE 8TH SENSOR ARRAY AND MULTICHANNEL SIGNAL PROCESSING WORKSHOP (SAM).2014,369-372. * |
Li, S.Near-Field Localization Algorithm Based on Sparse Reconstruction of the Fractional Lower Order Correlation Vector.12th International Conference on Wireless Algorithms, Systems, and Applications (WASA).2018,903-908. * |
Ulrich, M.Multi-Carrier MIMO Radar: A Concept of Sparse Array for Improved DOA Estimation.2016 IEEE RADAR CONFERENCE .2016,494-498. * |
刘自成.基于稀疏表示的雷达目标角度与距离估计.中国优秀硕士学位论文全文数据库.2014,I136-833. * |
吴世龙.星载稀布L阵对地面辐射源的高精度测向研究.上海航天.2006,24-27. * |
宋海岩.基于虚拟阵列的压缩波束形成方位估计方法.黑龙江工程学院学报.2018,32-36. * |
巩朋成;刘刚;黄禾;王文钦.频控阵MIMO雷达中基于稀疏迭代的多维信息联合估计方法.雷达学报.2014,(第02期),46-53. * |
干练.基于压缩感知的非规则阵DOA估计方法.声学与电子工程.2018,9-13. * |
张群飞,保铮,***.一种水下多目标方位、频率、距离联合估计新方法.电子学报.2014,(第09期),2-6. * |
王伟东.基于矢量水听器阵的迭代稀疏协方差矩阵拟合波达方向估计方法.西北工业大学学报.2020,14-23. * |
王秀.基于超声波传感器陈列DOA估计算法的研究与应用.中国优秀硕士学位论文全文数据库.2019,全文. * |
王秀红.基于稀疏表示的波达方向估计方法研究.中国优秀硕士学位论文全文数据库.2016,全文. * |
罗争.一种稀疏表示的二维DOA降维估计新算法.中国电子科学研究院学报.2013,501-506. * |
蒋佳佳.水声阵列探测算法及***关键技术研究.中国优秀硕士学位论文全文数据库.2014,全文. * |
许彬.基于块稀疏恢复的高频雷达空间角谱和极化谱联合估计方法.电波科学学报.2019,741-750. * |
许萌.基于压缩感知的空间虚拟阵列波束形成技术.光学与光电技术.2020,47-54. * |
许萌;杨阳;徐磊;***明.基于压缩感知的空间虚拟阵列波束形成技术.光学与光电技术.2010,(第02期),51-58. * |
赵春雷.基于压缩感知的高频地波雷达二维DOA估计.***工程与电子技术.2017,733-741. * |
陈峰.基于矢量水听器的水声目标方位估计方法研究.中国优秀硕士学位论文全文数据库.2018,全文. * |
顾飞飞.基于稀疏重构的极化敏感阵列信号参数估计.中国优秀硕士学位论文全文数据库.2016,I136-10. * |
黄麟舒;察豪;叶慧娟;徐慨.均匀线阵目标到达角估计的压缩感知方法研究.通信学报.2016,(第02期),172-178. * |
Also Published As
Publication number | Publication date |
---|---|
CN111679248A (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108828522B (zh) | 一种利用垂直阵lcmv波束形成的水下目标辐射噪声测量方法 | |
CN108226933B (zh) | 一种基于条纹干涉结构的深海宽带目标深度估计方法 | |
CN112269164B (zh) | 深海可靠声路径下基于干涉结构匹配处理弱目标定位方法 | |
CN112083404B (zh) | 一种基于多途特征匹配的单矢量水听器声源深度估计方法 | |
CN112987004B (zh) | 一种浅海环境下基于水平阵列的水面水下目标分类方法 | |
CN109444864B (zh) | 一种深海微弱多目标深度长时累积估计方法 | |
CN111679248B (zh) | 一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 | |
CN113640808B (zh) | 浅水区海底电缆埋深探测方法及装置 | |
CN107085216A (zh) | 一种基于单水听器的深海水声被动测距测深方法 | |
CN113011006B (zh) | 一种基于互相关函数脉冲波形匹配的目标深度估计方法 | |
CN105005026A (zh) | 一种近场目标声源三维被动定位方法 | |
CN104793212A (zh) | 利用声波海底反射实现主动声纳远程探测的方法 | |
CN104656073A (zh) | 三维成像声纳波束形成方法及在多核处理器上的实现方法 | |
CN112098938B (zh) | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 | |
CN112904316B (zh) | 一种针对海面波动的机载激光测深数据折射误差改正方法 | |
CN113534161B (zh) | 一种用于远距离定位水下声源的波束镜像聚焦方法 | |
Zhang et al. | Bistatic localization of objects in very shallow water | |
CN114909610B (zh) | 一种水下油气管道泄漏检测与定位的方法及控制*** | |
CN113009417B (zh) | 利用声场干涉特性的海底声学阵列阵形估计方法 | |
CN115902849A (zh) | 一种基于波束输出强度重采样的深海声源深度估计方法 | |
Tao et al. | Single hydrophone passive localization of transiting acoustic sources | |
RU2510045C2 (ru) | Фазовый гидролокатор бокового обзора | |
Shi et al. | An improved generalized inverse beamforming-noise source localization method using acoustic vector sensor arrays | |
CN117590369B (zh) | 一种深海目标深度估计方法、装置、设备及存储介质 | |
Guo et al. | Review of research on sonar imaging technology in China |
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 |