CN107271989B - 一种具有最优旋转角的低波束旁瓣mimo声纳阵型 - Google Patents
一种具有最优旋转角的低波束旁瓣mimo声纳阵型 Download PDFInfo
- Publication number
- CN107271989B CN107271989B CN201710387452.3A CN201710387452A CN107271989B CN 107271989 B CN107271989 B CN 107271989B CN 201710387452 A CN201710387452 A CN 201710387452A CN 107271989 B CN107271989 B CN 107271989B
- Authority
- CN
- China
- Prior art keywords
- array
- transmitting
- rotation angle
- angle
- receiving
- 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
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/521—Constructional features
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
-
- 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/52003—Techniques for enhancing spatial resolution of targets
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种具有最优旋转角的低波束旁瓣MIMO声纳阵型,将发射和接收圆环阵的半径以及两者的阵元个数约束在一定范围内,根据已有的成本预算和平台能够提供的尺寸在此范围中选择合适的阵列半径和阵元个数,在此基础上计算出发射圆环阵和接收圆环阵之间的最优旋转角范围,从而获得具有最优旋转角的低波束旁瓣的MIMO声纳阵型。
Description
技术领域
本发明属于阵列信号处理领域,特别涉及一种低波束旁瓣MIMO声纳阵型及获得最优旋转角确定最优阵型的方法。
背景技术
主动声纳为了获得360°视角,通常采用单个阵元发射和多元圆环阵接收的阵型,如航空吊放声纳等(支绍龙,刘颖,刘莹,张红.FLASH吊放声纳声学性能及单指向性阵元验证.中国声学学会水声学分会2015年学术会议,2015)。使用该阵型的声纳可以简化为单输入多输出(single-input multiple-output:SIMO)声纳。由于使用均匀圆环阵接收目标场景回波,会出现阵列波束旁瓣处于较高的水平的问题。在阵列结构不变时,为了抑制圆环阵的波束图旁瓣,波束图主瓣会随之展宽,导致角度分辨率下降(H.L.Van Trees,Detection,Estimation,and Modulation Theory,Optimum Array Processing,Wiley Interscience,2002)。为了在保持一定角度分辨率的前提下降低波束图旁瓣,传统的方法一般通过增加接收圆环阵阵元个数(接收圆环阵按照半波长布阵)以达到目的。但是,这种方法不但会导致声纳系统成本过高,也会带来圆环阵阵列尺寸过大难以由小型平台搭载等问题。
不同于传统的SIMO声纳,多输入多输出(multiple-input multiple-output:MIMO)声纳通过各阵元发射相互正交信号,在接收端进行匹配滤波分离信号并实现波束形成,从而获得更大的虚拟阵列孔径(刘雄厚,孙超,卓颉,马茜,潘浩.一种用于高分辨扇扫成像的MIMO阵列.航空学报,2014,35(9):2540-2550.)。为了获得期望的阵列孔径和足够的自由度,MIMO声纳一般采用稀疏发射阵和密布接收阵进行布阵设计。但是,如此设计的MIMO声纳容易突破阵列尺寸约束,导致其难以应用到小型平台。因此,如何在一定的阵列尺寸约束下,设计合理的MIMO声纳阵型并结合相应的波束优化以获得期望的低旁瓣波束图,成为MIMO声纳在小型平台上应用过程中急需解决的问题之一。
发明内容
本发明解决的技术问题是:为了在一定的约束下获得低旁瓣的MIMO声纳阵型,本发明提出一种基于圆环阵结构的低旁瓣MIMO声纳阵列结构及通过寻找最优旋转角确定最优阵型的方法。本发明将发射和接收圆环阵的半径以及两者的阵元个数约束在一定范围内,根据已有的成本预算和平台能够提供的尺寸在此范围中选择合适的阵列半径和阵元个数,在此基础上计算出发射圆环阵和接收圆环阵之间的最优旋转角范围,从而获得具有低旁瓣波束图的MIMO声纳阵型。
本发明的技术方案是:一种低波束旁瓣MIMO声纳阵列结构,包括一个M元发射圆环阵和一个N元接收圆环阵,两个圆环阵同心平行排列,两个圆环阵的间距LD的取值范围是信号中心频率对应波长λ的0.5倍到2倍;发射圆环阵由M个发射阵元组成,接收圆环阵由N个接收阵元组成;接收圆环阵相邻两个接收阵元之间的弧长为发射信号中心频率波长的一半;旋转角大于等于0度且小于360/N度;其中旋转角定义为:首先在发射圆环阵所在的平面上建立一个平面直角坐标系xOy,将发射圆环阵圆心所在的位置定义为原点O,任选发射圆环阵的一个阵元作为发射阵的第一个阵元,以原点和发射阵的第一个阵元的几何中心的连线作为Ox轴,规定从原点O到发射阵的第一个阵元的几何中心的方向为x轴的正方向,垂直于Ox轴的直线为Oy轴。将xOy平面的法线作为Oz轴,建立空间直角坐标系xOyz,合理选取z轴的正方向和y轴的正方向,使空间直角坐标系xOyz是右手系。在接收圆环阵所在的平面上建立一个平面直角坐标系x1O1y1,将接收圆环阵圆心所在的位置定义为原点O1,与Ox轴平行的直线作为O1x1轴,x1轴的正方向与x轴的正方向同向,与Oy轴平行的直线作为O1y1轴,y1轴的正方向与y轴的正方向同向。在平面直角坐标系x1O1y1中,将N个接收阵元的几何中心与原点O1连线,形成N条直线。以x1轴为始边,N条直线为终边,沿逆时针方向构成N个角度。定义在这N个角度中,最小的角度对应的阵元为接收阵的第一个阵元。(规定逆时针方向为正方向)。定义以x1轴为始边,接收阵的第一个阵元的几何中心投影到xOy平面的投影点与原点O连线为终边,构成的角度为旋转角,用字母表示。
本发明的进一步技术方案是:所述M元发射圆环阵的阵型为M个发射阵元周向均布且圆心处不设有发射阵元,此时发射阵元个数M满足3≤M≤6,接收阵元个数N满足3≤N≤4M;发射圆环阵半径Rt和接收圆环阵半径Rr的关系为Rr≤Rt≤1.5Rr。
本发明的进一步技术方案是:所述M元发射圆环阵的阵型还可以为M-1个发射阵元周向均布且圆心处设有发射阵元,此时发射阵元个数M满足4≤M≤7,接收阵元个数N满足3≤N≤4(M-1);发射圆环阵半径Rt和接收圆环阵半径Rr的关系为Rr≤Rt≤1.5Rr。
本发明的进一步技术方案是:所述最优旋转角通过以下步骤获得:
步骤一:初始参数设置:
(1)在最初所取的阵列结构下,通过改变旋转角的取值,得到具有不同旋转角的阵型:令即可,其中<>表示用四舍五入的方法取整数的运算,N是接收阵元个数,Z是整数集;旋转角的变化间隔满足单位为角度;
(2)主瓣指向角θ0:针对一个具有任意旋转角的MIMO声纳阵型,在主瓣指向不同的角度θ0时进行波束形成,主瓣指向角θ0满足θ0,min≤θ0≤θ0,max;如果选择第一种发射阵型,最小的主瓣指向角θ0,min满足0°≤θ0,min≤360°-min{[360/M,360/N],360}°,最大的主瓣指向角θ0,max满足θ0,max=θ0,min+min{[360/M,360/N],360}°-Δθ0,如果选择第二种发射阵型,最小的主瓣指向角θ0,min满足0°≤θ0,min≤360°-min{[360/(M-1),360/N],360}°,最大的主瓣指向角θ0,max满足θ0,max=θ0,min+min{[360/(M-1),360/N],360}°-Δθ0;其中min{a,b}表示取括号内的数a,b的最小值,[a,b]表示求括号内的数a,b的最小公倍数,M是发射阵元个数,N是接收阵元个数,Δθ0是进行波束形成的主瓣指向角的变化的间隔,0°≤Δθ0≤5°;那么,针对一个阵型,进行波束形成的次数K=<(θ0,max-θ0,min)/Δθ0>+1,其中<>表示用四舍五入的方法取整数的运算;
步骤二:根据期望的平均旁瓣级的范围,确定最优旋转角的范围
子步骤一:针对步骤二种所选取的阵型,获得某个主瓣指向θ0下的波束图,得到该阵型波束图的旁瓣级改变主瓣指向θ0,获得某个旋转角下的一组多个不同主瓣指向下的波束图的旁瓣级;改变旋转角可以获得多组不同旋转角下的波束图的旁瓣级;
子步骤二:根据子步骤一中的不同旋转角下的波束图的旁瓣级,求具有不同旋转角的阵型的平均旁瓣级
子步骤三:由期望的平均旁瓣级的范围确定最优旋转角的范围;其中期望的平均旁瓣级的范围满足: 代表在子步骤二得到的所有的平均旁瓣级中最小的平均旁瓣级,代表在子步骤二得到的所有的平均旁瓣级中最大的平均旁瓣级;则最优旋转角范围可由下列集合确定:其中∈表示属于,表示使得平均旁瓣级小于等于期望旁瓣级的最优旋转角的集合;
步骤三:在步骤二中得到的最优旋转角的范围中,任取一个值,作为对应最初所取的阵列结构下最优旋转角,从而获得具有最优旋转角的低波束旁瓣MIMO声纳阵型。
发明效果
本发明的技术效果在于:本发明的基本原理和实施方案经过了计算机数值仿真的验证,其结果表明:与未经优化设计的阵型相比,使用本发明中设计的阵型可以在不增加阵列尺寸和阵元个数的前提下,使波束图的旁瓣级降低到期望的水平。
附图说明
表1是本发明实施实例中4个MIMO阵列结构所涉及到的参数。
表2是本发明实施实例中针对4个MIMO阵列结构,分别设计出的最优旋转角的范围和最优旋转角对应的旁瓣级的范围。
图1是本发明获得期望阵型的流程图。
图2是本发明确定旋转角的流程图。
图3是本发明MIMO声纳的三维视图,其中图(a)是3发6收发射阵采用第一种圆心没有阵元的方式布阵的MIMO声纳的三维视图,图(b)是4发6收发射阵采用第二种圆心有阵元的方式布阵的MIMO声纳的三维视图,图(c)是4发8收发射阵采用第一种圆心没有阵元的方式布阵的MIMO声纳的三维视图,图(d)是4发8收发射阵采用第二种圆心有阵元的方式布阵的MIMO声纳的三维视图。
图4是本发明MIMO声纳的二维俯视图,其中图(a)是3发6收发射阵采用第一种圆心没有阵元的方式布阵的MIMO声纳的二维俯视图,图(b)是4发6收发射阵采用第二种圆心有阵元的方式布阵的MIMO声纳的二维俯视图,图(c)是4发8收发射阵采用第一种圆心没有阵元的方式布阵的MIMO声纳的二维俯视图,图(d)是4发8收发射阵采用第二种圆心有阵元的方式布阵的MIMO声纳的二维俯视图。
图5是本发明3发6收的MIMO声纳,在不同的旋转角和不同的主瓣指向角下的波束图的旁瓣级。
图6是本发明3发6收的MIMO声纳,在不同的旋转角下的平均旁瓣级。
图7是本发明3发6收的MIMO声纳,以旋转角为29°的阵型和旋转角为0°的阵型为例,在不同的主瓣指向角下的波束图。其中图(a)是旋转角为29°的阵型在不同主瓣指向角下的波束图,图(b)是旋转角为0°的阵型在不同主瓣指向角下的波束图。
图8是本发明3发6收的MIMO声纳,以旋转角为29°的阵型和旋转角为0°的阵型为例,在不同的主瓣指向角下的旁瓣级。
图9是本发明4发6收的MIMO声纳,在不同的旋转角和不同的主瓣指向角下的波束图的旁瓣级。
图10是本发明4发6收的MIMO声纳,在不同的旋转角下的平均旁瓣级。
图11是本发明4发6收的MIMO声纳,以旋转角为30°的阵型和旋转角为0°的阵型为例,在不同的主瓣指向角下的波束图。其中图(a)是旋转角为30°的阵型在不同主瓣指向角下的波束图,图(b)是旋转角为0°的阵型在不同主瓣指向角下的波束图。
图12是本发明4发6收的MIMO声纳,以旋转角为30°的阵型和旋转角为0°的阵型为例,在不同的主瓣指向角下的旁瓣级。
图13是本发明4发8收的MIMO声纳,在不同的旋转角和不同的主瓣指向角下的波束图的旁瓣级。
图14是本发明4发8收的MIMO声纳,在不同的旋转角下的平均旁瓣级。
图15是本发明4发8收的MIMO声纳,以旋转角为34°的阵型和旋转角为0°的阵型为例,在不同的主瓣指向角下的波束图。其中图(a)是旋转角为34°的阵型在不同主瓣指向角下的波束图,图(b)是旋转角为0°的阵型在不同主瓣指向角下的波束图。
图16是本发明4发8收的MIMO声纳,以旋转角为34°的阵型和旋转角为0°的阵型为例,在不同的主瓣指向角下的旁瓣级。
图17是本发明5发8收的MIMO声纳,在不同的旋转角和不同的主瓣指向角下的波束图的旁瓣级。
图18是本发明5发8收的MIMO声纳,在不同的旋转角下的平均旁瓣级。
图19是本发明5发8收的MIMO声纳,以旋转角为30°的阵型和旋转角为0°的阵型为例,在不同的主瓣指向角下的波束图。其中图(a)是旋转角为30°的阵型在不同主瓣指向角下的波束图,图(b)是旋转角为0°的阵型在不同主瓣指向角下的波束图。
图20是本发明5发8收的MIMO声纳,以旋转角为30°的阵型和旋转角为0°的阵型为例,在不同的主瓣指向角下的旁瓣级。
图21是本发明3发6收的MIMO声纳,旋转角以旋转角为29°为例的最优阵型,其中,(a)是该阵型的三维视图,图(b)是该阵型的二维俯视图。
图22是本发明4发6收的MIMO声纳,旋转角以旋转角为30°为例的最优阵型,其中(a)是该阵型的三维视图,图(b)是该阵型的二维俯视图。
图23是本发明4发8收的MIMO声纳,旋转角以旋转角为34°为例的最优阵型,其中,(a)是该阵型的三维视图,图(b)是该阵型的二维俯视图。
图24是本发明5发8收的MIMO声纳,旋转角以旋转角为30°为例的最优阵型,其中,(a)是该阵型的三维视图,图(b)是该阵型的二维俯视图。
具体实施方式
为了方便描述本发明的内容,首先作以下术语定义:
定义1、旋转角:参见图4,首先在发射圆环阵所在的平面上建立一个平面直角坐标系xOy,将发射圆环阵圆心所在的位置定义为原点O,任选发射圆环阵的一个阵元作为发射阵的第一个阵元,以原点和发射阵的第一个阵元的几何中心的连线作为Ox轴,规定从原点O到发射阵的第一个阵元的几何中心的方向为x轴的正方向,垂直于Ox轴的直线为Oy轴。将xOy平面的法线作为Oz轴,建立空间直角坐标系xOyz,合理选取z轴的正方向和y轴的正方向,使空间直角坐标系xOyz是右手系。在接收圆环阵所在的平面上建立一个平面直角坐标系x1O1y1,将接收圆环阵圆心所在的位置定义为原点O1,与Ox轴平行的直线作为O1x1轴,x1轴的正方向与x轴的正方向同向,与Oy轴平行的直线作为O1y1轴,y1轴的正方向与y轴的正方向同向。在平面直角坐标系x1O1y1中,将N个接收阵元的几何中心与原点O1连线,形成N条直线。以x1轴为始边,N条直线为终边,沿逆时针方向构成N个角度。定义在这N个角度中,最小的角度对应的阵元为接收阵的第一个阵元。(规定逆时针方向为正方向)。定义以x1轴为始边,接收阵的第一个阵元的几何中心投影到xOy平面的投影点与原点O连线为终边,构成的角度为旋转角,用字母表示。
在四个不同的阵列结构下,MIMO声纳阵型的三维视图如图3所示,MIMO声纳阵型的旋转角如图4所示。图3(a)、图3(b)、图3(c)和图3(d)分别表示圆心位置没有发射阵元的3发6收声纳阵型的三维视图、圆心位置有发射阵元的4发6收声纳阵型的三维视图、圆心位置没有发射阵元的4发8收声纳阵型的三维视图、圆心位置有发射阵元的5发8收声纳阵型的三维视图。图4(a)、图4(b)、图4(c)和图4(d)分别表示圆心位置没有发射阵元的3发6收声纳阵型的二维俯视图、圆心位置有发射阵元的4发6收声纳阵型的二维俯视图、圆心位置没有发射阵元的4发8收声纳阵型的二维俯视图、圆心位置有发射阵元的5发8收声纳阵型的二维俯视图。其中字母表示旋转角,θ0表示波束形成时的主瓣指向角。
需要说明的是,发射圆环阵和接收圆环阵的相对上下位置可以和图4中相反。
定义2、平均旁瓣级
针对一个具有任意旋转角的MIMO声纳阵型,主瓣指向不同的角度θ0时进行波束形成,能得到一组多个波束图。搜索主瓣以外的区域,可以得到一组多个旁瓣级定义平均旁瓣级为:具有此旋转角的阵型在所有感兴趣的主瓣指向下旁瓣级的平均值
定义3、均匀圆环阵
如果有一个阵型满足:在由阵元组成的圆周上,任意两个相邻阵元之间的弧长相同,那么这个阵型是均匀圆环阵。
本发明的主要内容有:
1.提出了一种具有最优旋转角的低波束旁瓣MIMO声纳阵型。M发N收圆环阵MIMO声纳阵型的发射圆环阵和接收圆环阵位于两个相互平行的平面中,且这两个平面的间距LD的取值范围是信号中心频率所对应波长λ的0.5倍到2倍,即:0.5λ≤LD≤2λ。N元接收圆环阵为均匀圆环阵(其中相邻阵元的弧长为发射信号中心频率对应的半波长),M元发射圆环阵的阵型有两种:第一种是M元均匀圆环阵并且圆心位置没有发射阵元,第二种是M-1元均匀圆环阵并且圆心位置有1个发射阵元。在设计MIMO声纳阵型时,任意选择一种发射阵型即可。使用第一种发射圆环阵时,发射阵元个数M满足3≤M≤6,接收阵元个数N满足3≤N≤4M;使用第二种发射圆环阵时,发射阵元个数M满足4≤M≤7,接收阵元个数N满足3≤N≤4(M-1)。为约束阵列尺寸并获得期望波束性能,发射圆环阵半径Rt和接收圆环阵半径Rr的取值范围满足Rr≤Rt≤1.5Rr。在这些参数确定的基础上,利用波束图来确定发射阵和接收阵之间的最优旋转角范围,使得在满足阵元个数约束、阵列半径约束、主瓣宽度约束、阵列加权系数2-范数约束时,该旋转角范围下MIMO声纳具有期望的低旁瓣级。设期望的平均旁瓣级为其中期望的平均旁瓣级的范围满足: 代表不同旋转角对应的平均旁瓣级中最小的平均旁瓣级,代表不同旋转角对应的平均旁瓣级中最大的平均旁瓣级。则最优旋转角范围可由下列集合确定:其中∈表示属于,表示使得平均旁瓣级小于等于期望旁瓣级的最优旋转角的集合。在设计阵型时旋转角在最优旋转角范围内任意选取。
2.通过计算机数值仿真,给出了在一定约束条件下的最优MIMO声纳阵型。将具有最优旋转角阵型的MIMO声纳波束性能与未经优化设计的MIMO声纳波束性能进行对比,证明了在满足阵列尺寸约束、阵元个数约束、波束图主瓣宽度约束、加权系数二范数约束时,所设计的具有最优旋转角的MIMO声纳阵型能够获得期望的低旁瓣波束效果。
本发明的技术方案
本发明解决现存问题所采用的技术方案可分为以下4个步骤:
1)根据尺寸约束、阵元个数约束、信号频率确定阵列结构。本发明中MIMO声纳由一个M元发射圆环阵和一个N元接收圆环阵构成。发射圆环阵和接收圆环阵位于两个相互平行的平面中,且这两个平面的间距LD的取值范围是信号中心频率所对应波长λ的0.5倍到2倍,即:0.5λ≤LD≤2λ。在0.5λ≤LD≤2λ的范围内选择这两个平面的间距LD。N元接收圆环阵为按照相邻阵元的弧长为发射信号中心频率对应的半波长布阵的均匀圆环阵,M元发射圆环阵的阵型有两种:选择一种发射阵型:第一种是M元均匀圆环阵并且圆心位置没有发射阵元,第二种是M-1元均匀圆环阵并且圆心位置有1个发射阵元。在设计MIMO声纳阵型时,任意选择一种发射阵型即可。使用第一种发射圆环阵时,发射阵元个数M满足3≤M≤6,接收阵元个数N满足3≤N≤4M;使用第二种发射圆环阵时,发射阵元个数M满足4≤M≤7,接收阵元个数N满足3≤N≤4(M-1),在此约束的范围内选择发射阵元个数M和接收阵元个数N。发射圆环阵半径Rt和接收圆环阵半径Rr的关系为Rr≤Rt≤1.5Rr。在Rr≤Rt≤1.5Rr的范围内选择发射圆环阵半径Rt和接收圆环阵半径Rr。
2)在步骤1)设计的阵列结构的基础上,设计具有不同旋转角的阵型
在设计发射圆环阵和接收圆环阵之间的旋转角时,令即可。感兴趣的旋转角的变化间隔满足通过改变旋转角的取值,可以得到具有不同旋转角的阵型。感兴趣的旋转角其中<>表示用四舍五入的方法取整数的运算,N是接收阵元个数,Z是整数集,感兴趣的旋转角的变化间隔的单位是角度。
3)根据期望的平均旁瓣级的范围,确定最优旋转角的范围
针对某个旋转角下的阵型,获得某个主瓣指向θ0下的波束图,可以得到该阵型波束图的旁瓣级改变主瓣指向,可以获得某个旋转角下的一组多个不同主瓣指向下的波束图的旁瓣级。改变旋转角,可以获得多组不同旋转角下的波束图的旁瓣级。由公式(3)求出具有不同旋转角的阵型的平均旁瓣级由期望的平均旁瓣级的范围确定最优旋转角的范围。其中期望的平均旁瓣级的范围满足:代表在不同旋转角对应的平均旁瓣级中最小的平均旁瓣级,代表在不同旋转角对应的平均旁瓣级中最大的平均旁瓣级。则最优旋转角范围可由下列集合确定:其中∈表示属于,表示使得平均旁瓣级小于等于期望旁瓣级的最优旋转角的集合。
感兴趣的主瓣指向角满足θ0,min≤θ0≤θ0,max。如果选择第一种发射阵型,最小的主瓣指向角θ0,min满足0°≤θ0,min≤360°-min{[360/M,360/N],360}°,如果选择第二种发射阵型,最小的主瓣指向角θ0,min满足0°≤θ0,min≤360°-min{[360/(M-1),360/N],360}°。其中min{a,b}表示取括号内的数a,b的最小值,[a,b]表示求括号内的数a,b的最小公倍数,M是发射阵元个数,N是接收阵元个数。最大的主瓣指向角θ0,max由公式(1)确定,不同主瓣指向下的波束图个数K由公式(2)决定。
主瓣指向角的最大值θ0,max可表示为:
其中,min{a,b}表示取括号内的数a,b的最小值,[a,b]表示求括号内的数a,b的最小公倍数,Δθ0是感兴趣的主瓣指向角的变化的间隔,0°≤Δθ0≤5°。
在某个旋转角下所能获得的波束图的个数K可表示为:
其中,<>表示用四舍五入的方法取整数的运算。
平均旁瓣级可表示为:
以某个旋转角下的阵型为例,给出其在某个主瓣指向角下的波束图和波束旁瓣级的计算过程。
以坐标原点为参考点,M元发射圆环阵对θj方向的响应向量at(u)可表示为M×1维列向量:
其中,xt,i为第i(i=1,2,……,M)个发射阵元的位置坐标:
xt,i=[xt,i,yt,i] (5)
uj是方位角θj对应的方向向量:
uj=[cosθj,sinθj]T (6)
[]T代表转置,下标t表示发射,xt,i为第i(i=1,2,……,M)个发射阵元沿x轴方向的坐标,yt,i为第i(i=1,2,……,M)个发射阵元沿y轴方向的坐标,λ为信号中心频率所对应波长。
以坐标原点为参考点,N元接收圆环阵对θj方向的响应向量ar(uj)可表示为N×1维列向量:
其中,xr,i为第i(i=1,2,……,N)个接收阵元的位置坐标:
xr,i=[xr,i,yr,i] (8)
下标r表示接收,xr,i为第i(i=1,2,……,N)个接收阵元沿x轴方向的坐标,yr,i为第i(i=1,2,……,N)个接收阵元沿y轴方向的坐标。
MIMO声纳的发射与接收联合阵列对θj方向的响应向量atr(uj)可表示为(刘雄厚,孙超,杨益新,卓颉.单基地多输入多输出声呐的方位分辨力[J].声学学报,2016,第41卷(第2期):163-173):
其中,下标tr表示发射与接收,表示Hadamard积,即:MIMO声纳的响应向量为发射阵列的响应向量与接收阵列的响应向量的直积。
当主瓣指向角为θ0时,使用加权向量w(u0)进行波束形成,可以得到波束图:
BMIMO(u;u0)=w(u0)Hatr(u)
(10)
其中,[]H代表共轭转置,u是空间扫描角度θ的方向向量,u0是主瓣指向角θ0的方向向量,w(u0)是在主瓣指向角为θ0时,波束形成的加权系数,atr(u)是MIMO声纳对空间扫描角度θ的响应向量。
以波束旁瓣最低作为优化目标,加权系数w(u0)可以通过求解以下最优化问题得到
其中,atr(u0)是MIMO声纳对主瓣指向θ0的响应向量,θ是空间扫描角度,BWSL是主瓣宽度,||w(u0)||2代表加权系数的2-范数,T0是加权系数2-范数的约束系数。
搜索主瓣以外方向上波束形成的输出,其最大值就是此波束图的旁瓣级SLL,计算公式如下:
4)确定具有低波束旁瓣的MIMO声纳阵型
由步骤1)设计出发射阵、接收阵(包括基本阵型、阵元数、阵列半径)和两者之间的间距LD,由步骤3)设计出在此基础上的最优旋转角的范围。结合步骤1)和步骤3)的参数设计出具有最优旋转角的低波束旁瓣的MIMO声纳阵型。
以4个不同阵列结构的MIMO声纳(分别对应2种发射阵阵型,每个发射阵阵型下有2个阵列结构)为例,给出本发明的实施实例。
步骤1:根据尺寸约束、阵元个数约束、信号频率确定阵列结构。
设发射信号为声波,其在水下的传播速度为1500米/秒,发射信号中心频率f=1500Hz,对应的波长λ=1m。
设计以下4个MIMO声纳的阵列结构。发射圆环阵和接收圆环阵位于两个相互平行的平面中,且这两个平面的间距LD=1.5m。这4个阵列结构的发射阵元个数M、接收阵元个数N、发射阵半径Rt、接收阵半径Rr、所采用的发射阵型以及之后波束形成的主瓣宽度BWSL、加权系数w的2-范数约束T0如表1所示,并且前两个阵列结构的发射阵半径Rt均为接收阵半径Rr的1倍,即:Rt=Rr;后两个阵列结构的发射阵半径Rt均为接收阵半径Rr的1.5倍,即:Rt=1.5Rr。
表1
这4个阵列结构的三维视图分别如图3(a)、图3(b)、图3(c)和图3(d)所示,二维俯视图分别如图4(a)、图4(b)、图4(c)和图4(d)所示。
2)在步骤1)设计的阵列结构的基础上,设计具有不同旋转角的阵型
当N=6时,令即:设置感兴趣的旋转角的变化间隔那么,感兴趣的旋转角Z是整数集。通过改变旋转角的取值,得到具有不同旋转角的阵型。
当N=8时,令即:设置感兴趣的旋转角的变化间隔那么,感兴趣的旋转角Z是整数集。通过改变旋转角的取值,得到具有不同旋转角的阵型。
3)根据期望的平均旁瓣级的范围,确定最优旋转角的范围
当M=3、N=6时,取主瓣指向角的最小值θ0,min=60°,主瓣指向角的间隔Δθ0=2°,主瓣指向角的最大值θ0,max=θ0,min+min{[360/M,360/N],360}°-Δθ0=178°。波束形成的次数针对步骤2)设计的每个旋转角下的阵型,做60次不同主瓣指向下的波束形成,得到其在不同主瓣指向下旁瓣级 如图5所示。由公式(3)求出其平均旁瓣级如图6所示。根据期望的平均旁瓣级的范围确定最优旋转角的范围为
以旋转角为29°的阵型(在最优旋转角的范围内)和旋转角为0°的阵型(在最优旋转角的范围外)为例,在不同的主瓣指向角下的波束图如图7所示。其中图(a)是旋转角为29°的阵型在不同主瓣指向角下的波束图,图(b)是旋转角为0°的阵型在不同主瓣指向角下的波束图。
以旋转角为29°的阵型(在最优旋转角的范围内)和旋转角为0°的阵型(在最优旋转角的范围外)为例,这两种阵型在60°≤θ0≤178°的主瓣指向下,对应的旁瓣级如图8所示。由图7和图8可知,在相同的阵列结构下,本发明所设计的具有最优旋转角的阵型比未经最优化设计的阵型的旁瓣级更低。
当M=4、N=6时,取主瓣指向角的最小值θ0,min=60°,主瓣指向角的间隔Δθ0=2°,最大的主瓣指向θ0,max=θ0,min+min{[360/(M-1),360/N],360}°-Δθ0=178°。波束形成的次数针对步骤2)设计的每个旋转角下的阵型,做60次不同主瓣指向下的波束形成,得到其在不同主瓣指向下旁瓣级 如图9所示。由公式(3)求出其平均旁瓣级如图10所示。根据期望的平均旁瓣级的范围确定最优旋转角的范围为
以旋转角为30°的阵型(在最优旋转角的范围内)和旋转角为0°的阵型(在最优旋转角的范围外)为例,在不同的主瓣指向角下的波束图如图11所示。其中图(a)是旋转角为30°的阵型在不同主瓣指向角下的波束图,图(b)是旋转角为0°的阵型在不同主瓣指向角下的波束图。
以旋转角为30°的阵型(在最优旋转角的范围内)和旋转角为0°的阵型(在最优旋转角的范围外)为例,这两种阵型在60°≤θ0≤178°的主瓣指向下,对应的旁瓣级如图12所示。由图11和图12可知,在相同的阵列结构下,本发明所设计的具有最优旋转角的阵型比未经最优化设计的阵型的旁瓣级更低。
当M=4、N=8时,取主瓣指向角的最小值θ0,min=90°,主瓣指向角的间隔Δθ0=2°,主瓣指向角的最大值θ0,max=θ0,min+min{[360/M,360/N],360}°-Δθ0=178°。波束形成的次数针对步骤2)设计的每个旋转角下的阵型,做45次不同主瓣指向下的波束形成,得到其在不同主瓣指向下旁瓣级 如图13所示。由公式(3)求出其平均旁瓣级如图14所示。根据期望的平均旁瓣级的范围确定最优旋转角的范围为
以旋转角为34°的阵型(在最优旋转角的范围内)和旋转角为0°的阵型(在最优旋转角的范围外)为例,在不同的主瓣指向角下的波束图如图15所示。其中图(a)是旋转角为34°的阵型在不同主瓣指向角下的波束图,图(b)是旋转角为0°的阵型在不同主瓣指向角下的波束图。
以旋转角为34°的阵型(在最优旋转角的范围内)和旋转角为0°的阵型(在最优旋转角的范围外)为例,这两种阵型在90°≤θ0≤178°的主瓣指向下,对应的旁瓣级如图16所示。由图15和图16可知,在相同的阵列结构下,本发明所设计的具有最优旋转角的阵型比未经最优化设计的阵型的旁瓣级更低。
当M=5、N=8时,取主瓣指向角的最小值θ0,min=90°,主瓣指向角的间隔Δθ0=2°,最大的主瓣指向θ0,max=θ0,min+min{[360/(M-1),360/N],360}°-Δθ0=178°。波束形成的次数针对步骤2)设计的每个旋转角下的阵型,做45次不同主瓣指向下的波束形成,得到其在不同主瓣指向下旁瓣级 如图17所示。由公式(3)求出其平均旁瓣级如图18所示。根据期望的平均旁瓣级的范围确定最优旋转角的范围为
以旋转角为30°的阵型(在最优旋转角的范围内)和旋转角为0°的阵型(在最优旋转角的范围外)为例,在不同的主瓣指向角下的波束图如图19所示。其中图(a)是旋转角为30°的阵型在不同主瓣指向角下的波束图,图(b)是旋转角为0°的阵型在不同主瓣指向角下的波束图。
以旋转角为30°的阵型(在最优旋转角的范围内)和旋转角为0°的阵型(在最优旋转角的范围外)为例,这两种阵型在90°≤θ0≤178°的主瓣指向下,对应的旁瓣级如图20所示。由图19和图20可知,在相同的阵列结构下,本发明所设计的具有最优旋转角的阵型比未经最优化设计的阵型的旁瓣级更低。
步骤4:确定具有低波束旁瓣的MIMO声纳阵型
由步骤1)设计出4种发射阵、接收阵(包括基本阵型、阵元数、阵列半径)和两者之间的间距LD,由步骤3)设计出这4种阵列结构分别对应的最优旋转角的范围如表2所示。结合二者的参数设计出具有低波束旁瓣的MIMO声纳阵型。其中,[]代表闭区间,[a,b]代表大于等于a且小于等于b的集合。
表2
在第一种3发6收的阵列结构下,最优旋转角的范围为最优阵型的旋转角在[4°,56°]之间,以具有旋转角为29°的阵型为例,其阵型如图21所示,其中图(a)是该阵型的三维视图,图(b)是该阵型的二维俯视图。
在第二种4发6收的阵列结构下,最优旋转角的范围为最优阵型的旋转角在[7°,53°]之间,以具有旋转角为30°的阵型为例,其阵型如图22所示,其中图(a)是该阵型的三维视图,图(b)是该阵型的二维俯视图。
在第三种4发8收的阵列结构下,最优旋转角的范围为最优阵型的旋转角在[4°,41°]之间,以具有旋转角为34°的阵型为例,其阵型如图23所示,其中图(a)是该阵型的三维视图,图(b)是该阵型的二维俯视图。
在第四种5发8收的阵列结构下,最优旋转角的范围为最优阵型的旋转角在[5°,40°]之间,以具有旋转角为30°的阵型为例,其阵型如图24所示,其中图(a)是该阵型的三维视图,图(b)是该阵型的二维俯视图。
根据实施实例,可以认为:本发明中所提出的最优旋转角下的MIMO声纳阵型,可以在相同的阵元个数约束、阵列半径约束、主瓣宽度约束、阵列加权2-范数约束下,比未经优化设计的MIMO声纳阵型具有更低的波束旁瓣。
Claims (3)
1.一种具有最优旋转角的低波束旁瓣MIMO声纳阵型,其特征在于:所述声纳阵型包括一个M元发射圆环阵和一个N元接收圆环阵,两个圆环阵同心平行排列,两个圆环阵的间距LD的取值范围是信号中心频率对应波长λ的0.5倍到2倍;发射圆环阵由M个发射阵元组成,接收圆环阵由N个接收阵元组成;接收圆环阵相邻两个接收阵元之间的弧长为发射信号中心频率波长的一半;旋转角大于等于0度且小于360/N度;其中旋转角定义为:首先在发射圆环阵所在的平面上建立一个平面直角坐标系xOy,将发射圆环阵圆心所在的位置定义为原点O,任选发射圆环阵的一个阵元作为发射阵的第一个阵元,以原点和发射阵的第一个阵元的几何中心的连线作为Ox轴,规定从原点O到发射阵的第一个阵元的几何中心的方向为x轴的正方向,垂直于Ox轴的直线为Oy轴;将xOy平面的法线作为Oz轴,建立空间直角坐标系xOyz,合理选取z轴的正方向和y轴的正方向,使空间直角坐标系xOyz是右手系;在接收圆环阵所在的平面上建立一个平面直角坐标系x1O1y1,将接收圆环阵圆心所在的位置定义为原点O1,与Ox轴平行的直线作为O1x1轴,x1轴的正方向与x轴的正方向同向,与Oy轴平行的直线作为O1y1轴,y1轴的正方向与y轴的正方向同向;在平面直角坐标系x1O1y1中,将N个接收阵元的几何中心与原点O1连线,形成N条直线;以x1轴为始边,N条直线为终边,沿逆时针方向构成N个角度;定义在这N个角度中,最小的角度对应的阵元为接收阵的第一个阵元;规定逆时针方向为正方向;定义以x1轴为始边,接收阵的第一个阵元的几何中心投影到xOy平面的投影点与原点O连线为终边,构成的角度为旋转角,用字母表示;M元发射圆环阵的阵型有两种:第一种是M元均匀圆环阵并且圆心位置没有发射阵元,第二种是M-1元均匀圆环阵并且圆心位置有1个发射阵元;所述最优旋转角通过以下步骤获得:
步骤一:初始参数设置:
(1)在最初所取的阵列结构下,通过改变旋转角的取值,得到具有不同旋转角的阵型:令即可,其中<>表示用四舍五入的方法取整数的运算,N是接收阵元个数,Z是整数集;旋转角的变化间隔满足单位为角度;
(2)主瓣指向角θ0:针对一个具有任意旋转角的MIMO声纳阵型,在主瓣指向不同的角度θ0时进行波束形成,主瓣指向角θ0满足θ0,min≤θ0≤θ0,max;如果选择第一种发射阵型,最小的主瓣指向角θ0,min满足0°≤θ0,min≤360°-min{[360/M,360/N],360}°,最大的主瓣指向角θ0,max满足θ0,max=θ0,min+min{[360/M,360/N],360}°-Δθ0,如果选择第二种发射阵型,最小的主瓣指向角θ0,min满足0°≤θ0,min≤360°-min{[360/(M-1),360/N],360}°,最大的主瓣指向角θ0,max满足θ0,max=θ0,min+min{[360/(M-1),360/N],360}°-Δθ0;其中min{a,b}表示取括号内的数a,b的最小值,[a,b]表示求括号内的数a,b的最小公倍数,M是发射阵元个数,N是接收阵元个数,Δθ0是进行波束形成的主瓣指向角的变化的间隔,0°≤Δθ0≤5°;那么,针对一个阵型,进行波束形成的次数K=<(θ0,max-θ0,min)/Δθ0>+1,其中<>表示用四舍五入的方法取整数的运算;
步骤二:根据期望的平均旁瓣级的范围,确定最优旋转角的范围
子步骤一:针对步骤二种所选取的阵型,获得某个主瓣指向θ0下的波束图,得到该阵型波束图的旁瓣级改变主瓣指向θ0,获得某个旋转角下的一组多个不同主瓣指向下的波束图的旁瓣级;改变旋转角可以获得多组不同旋转角下的波束图的旁瓣级;
子步骤二:根据子步骤一中的不同旋转角下的波束图的旁瓣级,求具有不同旋转角的阵型的平均旁瓣级
子步骤三:由期望的平均旁瓣级的范围确定最优旋转角的范围;其中期望的平均旁瓣级的范围满足: 代表在子步骤二得到的所有的平均旁瓣级中最小的平均旁瓣级,代表在子步骤二得到的所有的平均旁瓣级中最大的平均旁瓣级;则最优旋转角范围可由下列集合确定:其中∈表示属于,表示使得平均旁瓣级小于等于期望旁瓣级的最优旋转角的集合;
步骤三:在步骤二中得到的最优旋转角的范围中,任取一个值,作为最初所取的阵列结构下最优旋转角,从而获得具有最优旋转角的低波束旁瓣MIMO声纳阵型。
2.如权利要求1所述的一种具有最优旋转角的低波束旁瓣MIMO声纳阵型,其特征在于:所述M元发射圆环阵的阵型为M个发射阵元周向均布且圆心处不设有发射阵元,此时发射阵元个数M满足3≤M≤6,接收阵元个数N满足3≤N≤4M;发射圆环阵半径Rt和接收圆环阵半径Rr的关系为Rr≤Rt≤1.5Rr。
3.如权利要求1所述的一种具有最优旋转角的低波束旁瓣MIMO声纳阵型,其特征在于:所述M元发射圆环阵的阵型还可以为M-1个发射阵元周向均布且圆心处设有发射阵元,此时发射阵元个数M满足4≤M≤7,接收阵元个数N满足3≤N≤4(M-1);发射圆环阵半径Rt和接收圆环阵半径Rr的关系为Rr≤Rt≤1.5Rr。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710387452.3A CN107271989B (zh) | 2017-05-27 | 2017-05-27 | 一种具有最优旋转角的低波束旁瓣mimo声纳阵型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710387452.3A CN107271989B (zh) | 2017-05-27 | 2017-05-27 | 一种具有最优旋转角的低波束旁瓣mimo声纳阵型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107271989A CN107271989A (zh) | 2017-10-20 |
CN107271989B true CN107271989B (zh) | 2019-12-20 |
Family
ID=60065449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710387452.3A Active CN107271989B (zh) | 2017-05-27 | 2017-05-27 | 一种具有最优旋转角的低波束旁瓣mimo声纳阵型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107271989B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108761433B (zh) * | 2018-08-02 | 2022-04-08 | 西北工业大学 | 一种使用mimo声纳差合阵处理的高分辨成像方法 |
CN110196428B (zh) * | 2019-06-10 | 2022-08-23 | 西北工业大学 | 一种mimo声纳高分辨三维前视成像方法 |
CN110736976B (zh) * | 2019-10-25 | 2022-05-10 | 海鹰企业集团有限责任公司 | 一种任意阵形的声纳波束形成器性能估计方法 |
CN112099017B (zh) * | 2020-08-29 | 2022-10-14 | 西北工业大学 | 一种基于过扫描的圆环阵低旁瓣波束优化方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5754140A (en) * | 1996-07-05 | 1998-05-19 | Sperry Marine Inc. | Beam sharpened, low sidelobe antenna system |
CN102707258A (zh) * | 2012-06-05 | 2012-10-03 | 西安交通大学苏州研究院 | 基于l型传感器阵列的信号方位角和仰角联合估计方法 |
CN103454626A (zh) * | 2013-08-31 | 2013-12-18 | 西北工业大学 | 一种基于mimo阵列和稀疏阵列的扇扫成像方法 |
CN105929384A (zh) * | 2016-07-06 | 2016-09-07 | 西北工业大学 | 基于十字阵发射圆环阵接收的联合低旁瓣多波束形成方法 |
CN106329153A (zh) * | 2016-08-31 | 2017-01-11 | 电子科技大学 | 一种用于大型异构四维天线阵综合的联合优化方法 |
-
2017
- 2017-05-27 CN CN201710387452.3A patent/CN107271989B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5754140A (en) * | 1996-07-05 | 1998-05-19 | Sperry Marine Inc. | Beam sharpened, low sidelobe antenna system |
CN102707258A (zh) * | 2012-06-05 | 2012-10-03 | 西安交通大学苏州研究院 | 基于l型传感器阵列的信号方位角和仰角联合估计方法 |
CN103454626A (zh) * | 2013-08-31 | 2013-12-18 | 西北工业大学 | 一种基于mimo阵列和稀疏阵列的扇扫成像方法 |
CN105929384A (zh) * | 2016-07-06 | 2016-09-07 | 西北工业大学 | 基于十字阵发射圆环阵接收的联合低旁瓣多波束形成方法 |
CN106329153A (zh) * | 2016-08-31 | 2017-01-11 | 电子科技大学 | 一种用于大型异构四维天线阵综合的联合优化方法 |
Non-Patent Citations (2)
Title |
---|
"Devising MIMO Arrays for Underwater 3-D Short-Range Imaging";Liu X H 等;《Oceans》;20121231;第41卷(第2期);第1-7页 * |
"单基地多输入多输出声呐的方位分辨力";刘雄厚 等;《声学学报》;20161231;第41卷(第2期);第163-173页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107271989A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107271989B (zh) | 一种具有最优旋转角的低波束旁瓣mimo声纳阵型 | |
CN108828603B (zh) | 一种基于十字型的三维成像声纳阵列的稀疏优化方法 | |
CN108037692B (zh) | 大规模共形数字相控阵的波束控制方法 | |
US6205224B1 (en) | Circularly symmetric, zero redundancy, planar array having broad frequency range applications | |
CN109765562B (zh) | 一种三维前视声像声纳系统和方法 | |
WO2022170783A1 (zh) | 基于频率分集的阵列天线的波束控制方法、系统及控制器 | |
CN112162266B (zh) | 一种基于凸优化理论的共形阵二维波束优化方法 | |
CN104849696A (zh) | 一种mimo雷达天线阵稀疏布阵方法 | |
Rasekh et al. | Design of an adaptive wideband beamforming algorithm for conformal arrays | |
CN105929384A (zh) | 基于十字阵发射圆环阵接收的联合低旁瓣多波束形成方法 | |
CN109188386B (zh) | 基于改进二维esprit算法的mimo雷达高分辨率参数估计方法 | |
CN113271270B (zh) | 基于深度卷积神经网络的太赫兹超大规模信道估计方法 | |
Khan et al. | Null placement and sidelobe suppression in failed array using symmetrical element failure technique and hybrid heuristic computation | |
CN110988868B (zh) | 一种等效分数阶模式涡旋电磁波产生及成像方法 | |
CN103630887A (zh) | 一种多波束形成方法及使用该方法的多波束声纳 | |
CN110995331B (zh) | 一种基于多点精确控制的波束赋形方法 | |
CN105974377B (zh) | 一种对数字阵列雷达自适应调零技术的干扰方法 | |
CN101149435A (zh) | 一种u型阵波束形成加权方法 | |
US10931021B2 (en) | Antenna lens array for tracking multiple devices | |
CN110554383B (zh) | 用于微波频段的mimo环形阵列方位向成像方法及装置 | |
CN109001754B (zh) | 用于太赫兹频段mimo弧形阵列方位向成像方法 | |
CN116679306A (zh) | 一种基于系数修正的扩展孔径声纳成像方法 | |
CN112134604A (zh) | 使用相控阵天线搜索信号的方法和相控阵天线系统 | |
CN108761433B (zh) | 一种使用mimo声纳差合阵处理的高分辨成像方法 | |
CN103777179A (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 |