CN103969633B - 杂波中检测目标用mimo雷达发射波形的分级设计方法 - Google Patents
杂波中检测目标用mimo雷达发射波形的分级设计方法 Download PDFInfo
- Publication number
- CN103969633B CN103969633B CN201410122330.8A CN201410122330A CN103969633B CN 103969633 B CN103969633 B CN 103969633B CN 201410122330 A CN201410122330 A CN 201410122330A CN 103969633 B CN103969633 B CN 103969633B
- Authority
- CN
- China
- Prior art keywords
- clutter
- matrix
- signal
- waveform
- represent
- 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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
-
- 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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/282—Transmitters
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达技术领域,涉及集中式MIMO(Multiple?Input?Multiple?Output)雷达的发射波形设计,公开了一种杂波中检测目标用MIMO雷达发射波形的分级设计方法,能够提高杂波的抑制性能,可用于强杂波中的目标检测;其包括以下步骤:步骤1,建立MIMO雷达信号模型,获取杂波在空域上的强度以及在时间域上的主瓣响应先验知识;步骤2,基于杂波在空域上的强度先验知识,设计发射信号的自相关矩阵Rx;步骤3,联合优化主瓣合成信号(Mainlobe?Synthesized?Signal,MSS)和失配滤波器来抑制主瓣杂波信号;步骤4,在恒模条件约束下,设计发射波形矩阵。
Description
技术领域
本发明属于雷达技术领域,涉及集中式MIMO(MultipleInputMultipleOutput)雷达的发射波形设计,特别涉及一种杂波中检测目标用MIMO雷达发射波形的分级设计方法,能够提高杂波的抑制性能,可用于强杂波中的目标检测。
背景技术
在传统雷达中,发射波形通常是固定的。因此,接收到的信号是适应电磁环境的唯一来源,这会限制系统在复杂的传播和干扰环境下的性能。然而基于先验知识的雷达发射波形设计方法可以提高整个系统的灵活性,从而可以在复杂的传播和干扰环境下得到满意的性能。在杂波区域中检测目标,杂波相对于系统热噪声占主要地位。在输出信号能量与杂波能量比(SignaltoClutterRatio,SCR)低的情况下,我们很难得到满意的检测性能,特别是针对低速运动目标的检测。此时,增加发射信号的能量对提高检测性能是没有效果的,唯一的方法就是设计发射的波形。
受通信领域中MIMO技术及综合脉冲孔径雷达SIAR的启发,2003年Rabideau和Parker提出了MIMO雷达概念,见[RabideauD.J.andParkerP..UbiquitousMIMOMultifunctionDigitalArrayRadar[C].ConferenceRecordofthe37thAsilomarConferenceonSignals,SystemsandComputers,2003,vol.1,pp.1057-1064]。之后,这一概念在雷达领域引起了人们的广泛关注。根据发射天线和接收天线的间距大小,可以将MIMO雷达分为分布式MIMO雷达和集中式MIMO雷达两类。对于分布式MIMO雷达来说,由于各个天线对目标有不同的观测视角以及目标回波的独立性,在统计意义下,这类MIMO雷达可以克服目标的闪烁效应从而提高雷达对目标的探测性能。对于集中式MIMO雷达来说,其特点是阵元间距较小,具有自由地设计每副天线波形的能力。与相控阵雷达相比,集中式MIMO雷达的自由度提高了,从而集中式MIMO雷达呈现出更多的优越性,如参量分辨力的提高、更自由的发射方向图设计能力等,见[LiJ.andStoicaP..MIMORadarWithColocatedAntennas[J].IEEESignalProcessingMagazine,Sep.2007,vol.24,pp.106-114]。因此,可以通过设计集中式MIMO雷达每个阵元的发射波形,从而提高它在目标检测、跟踪与识别上的系统性能。
由于MIMO雷达可以通过各阵元发射不同的信号,从而获得某种实际需要的方向图,为此根据实际需要设计发射波形成为必需。目前发射波形设计的主要方法有正交波形设计、发射方向图的波形设计以及基于先验知识的波形设计。
正交波形设计是MIMO雷达波形设计的基础,它主要考虑波形的自相关和互相关性质。然而,正交波形的发射方向图是各向同性的。与相控阵雷达相比,主瓣增益损失可能会降低正交波形的MIMO雷达性能,这限制了正交波形的应用。
在给定期望的发射方向图下,PetreStoica和JianLi提出了在恒模约束下的波形优化设计算法。具体设计方法是先设计发射信号的自相关矩阵R,再根据已设计的自相关矩阵R合成发射波形矩阵X。对于自相关矩阵R的设计,PetreStoica和JianLi提出了方向图匹配设计和最小化旁瓣设计,如[StoicaP.,LiJ.,XieY..OnprobingsignaldesignforMIMOradar[J].IEEETrans.onSignalProcessing.2007,vol.55(8).4151-4161]文中所述对设计中所建立的凸规划模型可采用凸优化工具包cvx求解,cvx的具体用法见[M.GrantandS.Boyd.CVX:Matlabsoftwarefordisciplinedconvexprogramming.http://stanford.edu/~boyd/cvx,Dec.2008]。该设计方法是基于最小二乘准则的,使由自相关矩阵R得到的方向图尽可能的逼近期望方向图。对于根据自相关矩阵R合成发射波形矩阵X,PetreStoica和JianLi提出了cyclicalgorithm(CA)算法,见[StoicaP.,LiJ.,ZhuX..WaveformSynthesisforDiversity-BasedTransmitBeampatternDesign.IEEETrans.onSignalProcessing.2008,vol.56(6).2593-2598]。该设计方法是基于最小二乘准则的循环算法,逐步使得恒模信号发射波形矩阵X的相关矩阵逼近已优化的自相关矩阵R。
但是,这些算法适用于均匀杂波中。在实际中,环境杂波通常是异构的,这将会降低已存在的发射方向图设计算法的性能。基于环境杂波的先验知识,BenjaminFriendlander提出了一种MIMO雷达的发射波形设计算法,其具体方法是基于已经得到的杂波响应,通过最大化系统输出信号能量与杂波能量比SCR来获得用于目标检测的最优发射波形,见[Benjamin.,WaveformDesignforMIMORadars.IEEETrans.onAerospaceandElectronicsystems.2007,vol.43(3).1227-1238]。但是,该算法的代价函数非常复杂并且没有考虑到发射波形的恒模约束条件。在恒模约束条件下,代价函数是一个NP(Non-deterministicPolynomial)难问题,这将很难找到一个满意的解。
发明内容
本发明的目的是为了克服上述基于杂波先验知识的波形设计方法的不足,提出了一种杂波中检测目标用MIMO雷达发射波形的分级设计方法,能够提高杂波的抑制性能,可用于强杂波中的目标检测,提高MIMO雷达在杂波区域中检测目标的性能。
为达到上述目的,本发明采用以下技术方案予以实现。
1、一种杂波中检测目标用MIMO雷达发射波形的分级设计方法,其特征在于,包括以下步骤:
步骤1,建立MIMO雷达信号模型,获取杂波在空域上的强度以及在时间域上的主瓣响应;
步骤2,根据空域上杂波强度以及MIMO雷达信号模型,分别计算出目标回波能量Pt和θk方向上杂波的回波能量Pck,得到输出信号能量与杂波能量比SCR为其中,k=1,2,…,K,K表示杂波区域在空域上的等分个数;在自相关矩阵半正定以及发射天线发射能量恒定条件约束下,通过最大化输出信号能量与杂波能量比SCR得到最优的发射信号的自相关矩阵Rx;
步骤3,根据MIMO雷达的信号模型分别得到主瓣合成信号和主瓣杂波响应分别为st=XaT(θt)和qt=[qt(1),qt(2),…,qt(Nc)]T,其中θt表示目标所在的方位角,aT(θt)表示θt方向上的发射导向矢量,NP表示信号码长,N表示阵元个数,C表示复数域;发射信号波形矩阵X表示为xm=[xm(1),xm(2),…,xm(Np)]T,其中xm(n)(m=1,…,N)表示第m个阵元天线发射的恒模信号, 其中d表示阵元间距,λ表示信号波长;利用联合优化主瓣合成信号和失配滤波器的设计得到最优的主瓣合成信号st以及失配滤波器h;
步骤4,在恒模条件约束下,设计发射波形矩阵:对已得到的发射信号自相关矩阵Rx,采用cyclicalgorithm算法设计发射信号波形矩阵XCA;利用得到的主瓣合成信号st,最终得到发射波形矩阵为:
Xopt=diag(ejφ)XCA
其中,φ=arg(st)-arg(XCAaT(θt)),arg(a)表示矢量a的相位角度。
上述技术方案的特点和改进在于:
(1)步骤2的子步骤为:
2.0)根据空域上杂波强度以及MIMO雷达信号模型,目标回波的能量表示为:
其中Rx=XHX表示发射信号的自相关矩阵,tr(·)表示计算矩阵迹的运算,(·)H表示矩阵的共轭转置运算,||·||2表示Euclidean范数,ζ表示目标的复反射系数;
方位角θk方向上的杂波回波能量可近似表示为:
其中表示方位角θk方向上的杂波强度,w是一个N×1维向量,ωk=wHaR(θk),
发射信号的自相关矩阵Rx的设计如下:
2.1)写出目标函数:
输出信号能量与杂波能量比SCR如<3>式所示:
其中
2.2)最大化<3>式,建立如下凸规划数学模型:
s.t.Rx≥0<4>
Rx(i,i)=E,i=1,2,…,N
其中Rx≥0表示矩阵Rx为半正定矩阵,[M]ij表示矩阵M的第(i,j)个元素,E表示每根阵元天线的发射能量,N表示阵元个数;
将<4>式中的代价函数转化为最小化旁瓣杂波能量,数学模型等价于<5>式:
s.t.tr(RxRat)=N2E(1-β)<5>
Rx≥0
[Rx]ii=Efori=1,2,…,N
其中β是用来控制主瓣增益损失的一个小的正数,N表示阵元个数;
2.3)通过凸优化工具包cvx可以求得<5>式中最优的发射信号的自相关矩阵Rx。
(2)步骤3的子步骤为:
3.0)首先,根据MIMO雷达信号模型,经过滤波器后输出的信号分量和杂波分量为<6>式:
zt=ζHXaT(θt)
其中矩阵H是(Np+Nh-1)×Np维的矩阵,Qk为(Np+NC+Nh-2)×(Np+Nh-1)维矩阵,Nh表示失配滤波器的长度,分别如下所示:
然后,根据<6>式中zt的中心元素表示主瓣自相关函数,定义自相关主瓣增益损失为:
其中[M]i,:表示矩阵第i行所有元素;
自相关峰值旁瓣电平表示为:
其中max{·}表示选取所有元素中最大值的运算;
杂波峰值电平表示为:
CPL=max{|zc(i)|2},i=1,…,Np+Nh+Nc-2<11>
得到主瓣合成信号波形和主瓣杂波响应分别为:
st=XaT(θt)<12>
qt=[qt(1),qt(2),…,qt(Nc)]T<13>
得到主瓣杂波的回波表达式为:
ct=qt*st<14>
经接收失配滤波器h后,输出的主瓣杂波回波信号和目标回波信号分别为:
zct=h*qt*st=qt*h*st<15>
zt=h*ζst<16>
分别写成矩阵形式,输出的主瓣杂波回波信号和目标回波信号可以重新写为:
zct=QtHst=QtSth<17>
zt=Hst=ζSth<18>
其中矩阵St为(Np+Nh-1)×Nh维矩阵,表示为:
3.1)在控制自相关主瓣增益损失情况下,最小化自相关旁瓣电平APSL和杂波峰值电平CPL,得到数学模型为:
|zct(i)|2≤p1,i=1,2,…,NC+NP+Nh-2<20>
其中λ1和λ2都是小的正数,分别用来控制自相关主瓣增益损失以及目标函数中杂波峰值电平和自相关旁瓣电平比例的正尺度因子,p1表示自相关峰值旁瓣电平,p2表示杂波峰值电平。
不失一般性,将主瓣合成信号的能量归一为单位能量,数学模型重新写为:
|zc(i)|2≤p1,i=1,2,…,Nc+Np+Nh-2<21>
3.2)对<21>式按如下步骤求解:
3.2a)按式<12>设置初始的主瓣合成信号矩阵st,按式<13>计算主瓣杂波响应向量qt;
3.2b)令i=0,设置Nh×1维初始失配滤波器hi为
其中表示将st中对第n个元素取共轭。同时设置小的正数λ1、λ2和ξ,利用<8>式计算矩阵Qt;
3.2c)令i值增1,利用<7>式计算矩阵H,同时固定hi-1,利用cvx优化工具包寻找<22>式的最优解:
s.t.max|QtHst|≤p1<22>
max|Hst-d|≤p2
||st||2≤1
其中表示期望的自相关函数。求解得到最优的主瓣合成信号为 其中arg(a)表示矢量a的相位角度;
3.2d)固定sti,通过<23>式优化接收滤波器hi:
max|QtSthi|≤p1
max|Sthi-d|≤p2
其中矩阵St是步骤3.2c)中取得的sti通过<19>式计算得到的;
3.2e)重复步骤3.2c)和步骤3.2d),直到两次相邻两次循环得到的主瓣合成信号sti与满足终止条件||sti-st(i-1)≤ξ,停止并输出最优解st=sti,h=hi;否则,转入步骤3.2c)中继续循环。
(3)步骤4的子步骤为:
4.1)在恒模条件约束下,发射波形矩阵设计的数学模型写为:
s.t.XaT(θt)=st<24>
在每次发射过程中,如果给发射波形加任一相位,MIMO雷达的发射方向图将保持不变,即:
其中X表示改变相位后的波形矩阵,表示增加的相位向量,K表示将整个空域等分的个数,diag(·)表示对角矩阵;
4.2)利用这一性质,<24>式简化成在满足给定发射方向图条件下的发射波形矩阵设计问题,等价的数学模型为:
利用CA算法求解<26>式,可以得到初始发射波形矩阵XCA;
4.3)利用步骤3中得到的最优主瓣合成信号st,最终的发射波形矩阵通过改变初始发射波形矩阵XCA的相位得到,如<27>式所示:
Xopt=diag(ejφ)XCA<27>
其中φ=arg(st)-arg(XCAaT(θt)),arg(a)表示矢量a的相位角度。
本发明相对现有技术具有以下显著优点:(1)本发明中引入了分级的思想,它可以将复杂的波形设计问题转化为一些简单的优化问题。在分级思想中,分别在空间域和时间域上对发射方向图和发射波形进优化设计;实现了空域上的旁瓣杂波的抑制,时域的主瓣杂波的抑制。(2)基于空间域上杂波强度的先验知识,本发明提出了基于发射方向图优化的凸规划模型,它可以提高接收波束形成后在空间域上的输出信号能量与杂波能量比SCR,用来抑制旁瓣杂波。(3)基于主瓣杂波响应的先验知识,本发明提出了联合优化主瓣合成信号和失配滤波器的优化算法,用来抑制主瓣杂波。(4)给定最优的发射方向图和主瓣合成信号条件下,本发明采用基于CA循环算法的发射波形矩阵设计算法,提高了雷达在杂波中检测目标的性能。
附图说明
下面结合附图和具体实施方式对本发明做进一步详细说明。
图1是本发明的分级波形设计方案的设计框图;
图2是本发明在均匀线性阵列条件下一种简化的收发共置的集中式MIMO雷达示意图(图中d表示阵元间距,θ表示目标的方位角度);
图3是本实验中的杂波信号响应图(图中横坐标和纵坐标分别表示方位角和时间延迟序号,横坐标的单位为度);
图4是本实验中杂波信号在空域上的强度分布图(图中横坐标和纵坐标分别表示方位角和杂波的强度,单位分别为角以及分贝(dB));
图5是本实验中位于0°目标方向上的主瓣杂波响应(其中横坐标和纵坐标分别表示时间延迟序号以及幅度);
图6是用本发明得到的波形1的发射方向图同波形2和波形3的发射方向图(图中横纵坐标分别表示方位角和发射方向图,单位分别为角和分贝(dB));
图7是用本发明得到波形1的主瓣合成信号相位图(图中横坐标和纵坐标分别表示时间序号和相位,纵坐标的单位为弧度);
图8是用本发明得到波形1的主瓣合成信号功率谱密度(PSD)(图中横坐标和纵坐标分别表示频率序号和幅度);
图9是用本发明得到波形1的主瓣合成信号经过失配滤波器后自相关图(图中横坐标和纵坐标分别表示时间序号和自相关幅度,纵坐标单位为分贝(dB));
图10是本发明中用于比较的波形2的主瓣合成信号经过失配滤波器后自相关图(图中横坐标和纵坐标分别表示时间序号和自相关幅度,纵坐标单位为分贝(dB));
图11是利用本发明得到的波形1的主瓣杂波输出图(图中横坐标和纵坐标分别表示时间序号和幅度);
图12是利用用于比较的波形2的主瓣杂波输出图(图中横坐标和纵坐标分别表示时间序号和幅度);
图13是利用用于比较的波形3的主瓣杂波输出图(图中横坐标和纵坐标分别表示时间序号和幅度);
图14是利用本发明得到的波形1的杂波输出图(图中横坐标和纵坐标分别表示时间序号和幅度);
图15是利用用于比较的波形2的杂波输出图(图中横坐标和纵坐标分别表示时间序号和幅度);
图16是利用用于比较的波形3的杂波输出图(图中横坐标和纵坐标分别表示时间序号和幅度);
图17是利用本发明得到的波形1与用于比较的波形2和波形3在不同主瓣增益损失情况下对输出信号能量与杂波能量比的影响(图中横坐标和纵坐标分别表示主瓣增益损失和提升大小,单位均为分贝(dB))。
具体实施方式
参照图1,说明本发明的杂波中检测目标用MIMO雷达发射波形的分级设计方法,其包括以下步骤:
步骤1,建立MIMO雷达信号模型,获取杂波在空域上的强度以及在时间域上的主瓣响应先验知识。
其中,MIMO雷达系统由N个间距为d的各向同性的发射和接收共置的均匀线阵组成。
其中,发射信号模型中,发射信号波形矩阵X表示为X=[x1,x2,…,xN],xm=[xm(1),xm(2),…,xm(Np)]T,xm(n)(m=1,…,N)表示第m个阵元天线发射的恒模信号,N表示阵元个数;远场方位向θ处的合成信号表示为s=Xax(θ),其中aT(θ)=[1,e-j2πdsin(θ)/λ,…,e-j2π(N-1)dsin(θ)/λ]T,d表示阵元间距,λ表示信号波长。
其中,在回波信号模型中,将杂波所在的空域和时域分别用Δθ和Ts间隔大小进行离散化,得到杂波响应矩阵为Q=[q1,q2,…,qK],其中qk=[qk(1),qk(2),…,qk(Nc)]T表示在θk方向上的杂波响应,qk(n)表示θk方向上时间延迟为nTs的路径反射系数,Nc表示杂波响应向量的长度,k=1,2,…,K,K表示杂波区域在空域上的等分个数。
由杂波响应矩阵可以得到杂波的回波信号模型,目标回波信号和杂波回波信号在分别经过波束形成器以及失配滤波器后得到最终输出信号模型。根据杂波响应矩阵利用统计学方法得到空域上的杂波强度知识,利用目标的方位角信息可以得到杂波在时域上的主瓣响应知识。
具体地,
1)建立MIMO雷达的信号模型,即MIMO雷达的发射信号和输出信号模型:
1.1a)发射信号模型分析:
如图2所示,MIMO雷达由N个间距为d的各向同性的发射和接收共置的均匀线阵组成。xm(t)表示第m个阵元发射的有限能量为E的基带波形,它的频谱范围为[-B/2,B/2],其中B表示基带波形的带宽。Tp和Ts分别表示脉冲的持续时间和采样周期,需要Np=Tp/Ts个基带波形来完整表示发射波形。第m个阵元的发射波形矩阵为:
xm=[xm(1),xm(2),…,xm(Np)]T<28>
其中
用大小为Δθ的角度间隔将空域离散后,方位角θk上的远场合成信号可以写为:
sk=XaT(θk)<29>
其中X=[x1,x2,…,xN]表示Np×N维离散的发射波形矩阵, 表示θk方向上的发射导向矢量,d表示阵元间距,λ表示信号波长,N表示发射阵元个数,符号(·)T表示矩阵转置运算。
1.1b)接收回波信号模型分析:
将杂波所在的空域和时域分别用Δθ和Ts离散化,杂波响应矩阵可以表示为:
Q=[q1,q2,…,qK]<30>
其中qk=[qk(1),qk(2),…,qk(Nc)]T表示在θk方向上的杂波响应向量,qk(n)表示θk方向上时间延迟为nTs的路径的反射系数,Nc表示杂波响应向量的长度,K表示杂波区域在空域上的等分个数。
方位角θk方向上的杂波回波信号为:
ck=qk*sk<31>
其中符号*表示卷积运算,sk表示在θk方向上的合成信号。
<31>式写成矩阵形式,等价于<32>式:
ck=Skqk<32>
其中Sk是一个(Np+Nc-1)×Nc维矩阵,形式为:
Nc表示杂波响应向量的长度,NP表示发射信号的码长。
在MIMO雷达中,设未知目标为点目标,方位角为θt,目标的回波信号为:
rt=ζst<33>
其中ζ表示目标的复反射系数,st表示在θt方向上的合成信号。
MIMO雷达接收到的目标回波信号以及杂波回波信号分别为:
其中aR(θ)=[1,e-j2πdsin(θ)/λ,…,e-j2π(N-1)dsin(θ)/λ]T表示θ方向上的接收导向矢量,K表示将杂波所在区域在空域上的等分个数。
在杂波区域中,杂波能量通常大于噪声能量,因此可以忽略系统的噪声。总的接收信号矩阵可以近似表示为:
其中,V=[v1,v2,…,vN]表示(Np+Nc-1)×N维接收回波信号矩阵,vm是第m个阵元接收到的(Np+Nc-1)×1维回波信号矩阵。
1.1c)输出信号模型分析:
接收回波信号经过接收波束形成向量w后,最终得到的输出信号为:
其中w是一个N×1维向量,为了发射方向图设计的方便,用w=μaR(θt)代替w,其中是归一化因子,||·||2表示Euclidean范数,输出信号中的信号分量yt和杂波分量yc可以分别写为:
yt=ζXaT(θt)
其中
选择失配滤波器h作为接收滤波器,经过滤波器后输出的信号分量zt以及杂波分量zc可以分别表示为:
zt=h*(ζXaT(θt))
其中h是一个Nh×1维的向量,Nh表示失配滤波器的长度。
为了数学推导的方便,经过滤波器后输出的信号分量和杂波分量重新写为<39>式:
zt=ζHXaT(θt)
其中矩阵H是(Np+Nh-1)×Np维的矩阵,Qk为(Np+NC+Nh-2)×(Np+Nh-1)维矩阵,Nh表示失配滤波器的长度,分别如下所示:
2)获取杂波的先验知识
在空域上,各个方向上的杂波包括一组不同时延的杂波路径。从统计意义上说,一个方向上的杂波反射特性可以近似认为是零均值的复高斯分布。从而,角度θk上的杂波反射特性的方差可以近似地认为是空域上的杂波强度,例如:
其中τ表示初始时刻,Δτ表示时延采样的间隔,K表示将杂波所在空域等分的个数。
假设目标的方位角为θt,根据<30>式可以得到主瓣杂波响应为:
qt=[qt(1),qt(2),…,qt(Nc)]T<43>
其中Nc表示杂波响应向量的长度。
步骤2,基于杂波在空域上的强度先验知识,设计发射信号的自相关矩阵Rx。根据空域上杂波强度先验知识以及MIMO雷达信号模型,分别计算出目标回波能量Pt和θk方向上杂波的回波能量Pck,得到输出信号能量与杂波能量比SCR为其中,k=1,2,…,K,K表示杂波区域在空域上的等分个数。在自相关矩阵半正定以及发射天线发射能量恒定条件约束下,通过最大化输出信号能量与杂波能量比SCR得到最优的发射信号的自相关矩阵Rx。
具体地,根据步骤1中经过波束形成后得到的输出信号模型,目标回波的能量可以表示为:
其中Rx=XHX表示发射信号的自相关矩阵,tr(·)表示计算矩阵迹的运算,(·)H表示矩阵的共轭转置运算,||·||2表示Euclidean范数。
已知杂波在空间上强度的先验知识条件下,方位角θk方向上的杂波回波能量可近似表示为:
其中ωk=wHaR(θk),
根据以上分析,已知杂波在空间上强度的先验条件下,发射信号的自相关矩阵Rx的设计可按如下步骤进行:
2.1)写出目标函数:
输出信号能量与杂波能量比SCR如<46>式所示:
其中
2.2)不失一般性,去掉式<46>中的常数ζ不会影响最终的优化结果。最大化<46>式,建立如下凸规划数学模型:
s.t.Rx≥0<47>
Rx(i,i)=E,i=1,2,…,N
其中Rx≥0表示矩阵Rx为半正定矩阵,[M]ij表示矩阵M的第(i,j)个元素,E表示每根阵元天线的发射能量(即每个阵元发射的有限能量)。
实际上,最小化旁瓣杂波能量可以通过优化发射方向图来实现,但会带来一定的主瓣增益的损失,因此<47>式中的代价函数可以转化为最小化旁瓣杂波能量,数学模型等价于<48>式:
s.t.tr(RxRat)=N2E(1-β)<48>
Rx≥0
[Rx]ii=Efori=1,2,…,N
其中β是用来控制主瓣增益损失的一个小的正数,N表示阵元个数。
2.3)通过凸优化工具包cvx可以求得<48>式中最优的发射信号的自相关矩阵Rx。
步骤3,联合优化主瓣合成信号(MainlobeSynthesizedSignal,MSS)和失配滤波器来抑制主瓣杂波信号。根据MIMO雷达的信号模型分别得到主瓣合成信号和主瓣杂波响应分别为st=XaT(θt)和qt=[qt(1),qt(2),…,qt(Nc)]T,其中θt表示目标所在的方位角,aT(θt)表示θt方向上的发射导向矢量,NP表示信号码长,N表示阵元个数,C表示复数域;发射信号波形矩阵X表示为xm=[xm(1),xm(2),…,xm(Np)]T,其中xm(n)(m=1,…,N)表示第m个阵元天线发射的恒模信号, 其中d表示阵元间距,λ表示信号波长;利用联合优化主瓣合成信号和失配滤波器的设计得到最优的主瓣合成信号st以及失配滤波器h。
具体地,根据<39>式中zt的中心元素表示主瓣自相关函数,定义自相关主瓣增益损失(Auto-correlationMainlobeGainLoss,AMGL)为:
其中[M]i,:表示矩阵第i行所有元素。
自相关峰值旁瓣电平(Auto-correlationPeakSidelobeLevel,APSL)可以表示为:
其中max{·}表示选取所有元素中最大值的运算。
杂波峰值电平(ClutterPeakLevel,CPL)可以表示为:
CPL=max{|zc(i)|2},i=1,…,Np+Nh+Nc-2<51>
根据步骤1中的信号模型可以得到主瓣合成信号波形和主瓣杂波响应分别为:
st=XaT(θt)<52>
qt=[qt(1),qt(2),…,qt(Nc)]T<53>
得到主瓣杂波的回波表达式为:
ct=qt*st<54>
经接收失配滤波器h后,输出的主瓣杂波回波信号和目标回波信号分别为:
zct=h*qt*st=qt*h*st<55>
zt=h*ζst<56>
分别写成矩阵形式,输出的主瓣杂波回波信号和目标回波信号可以重新写为:
zct=QtHst=QtSth<57>
zt=Hst=ζSth<58>
其中矩阵St为(Np+Nh-1)×Nh维矩阵,表示为:
根据以上知识,联合优化主瓣合成信号MSS和失配滤波器来抑制主瓣杂波步骤为:
3.1)在控制自相关主瓣增益损失情况下,最小化自相关旁瓣电平APSL和杂波峰值电平CPL,得到数学模型为:
|zct(i)|2≤p1,i=1,2,…,NC+NP+Nh-2<60>
其中λ1和λ2都是小的正数,分别用来控制自相关主瓣增益损失以及目标函数中杂波峰值电平和自相关旁瓣电平比例的正尺度因子,p1表示自相关峰值旁瓣电平,p2表示杂波峰值电平。
不失一般性,将主瓣合成信号的能量归一为单位能量,数学模型重新写为:
|zc(i)|2≤p1,i=1,2,…,Nc+Np+Nh-2<61>
3.2)对<61>式按如下步骤求解:
3.2a)按式<52>设置初始的主瓣合成信号矩阵st,按式<53>计算主瓣杂波响应向量qt;
3.2b)令i=0,设置Nh×1维初始失配滤波器hi为
其中表示将st中对第n个元素取共轭。同时设置小的正数λ1、λ2和ξ,利用<41>式计算矩阵Qt;
3.2c)令i值增1,利用<40>式计算矩阵H,同时固定hi-1,利用cvx优化工具包寻找<62>式的最优解
s.t.max|QtHst|≤p1<62>
max|Hst-d|≤p2
||st||2≤1
其中表示期望的自相关函数。求解得到最优的主瓣合成信号为其中arg(a)表示矢量a的相位角度(弧度)。
3.2d)固定sti,通过<63>式优化接收滤波器hi:
max|QtSthi|≤p1
max|Sthi-d|≤p2
其中矩阵St是步骤3.2c)中取得的sti通过<59>式计算得到的。
3.2e)重复步骤3.2c)和步骤3.2d),直到两次相邻两次循环得到的主瓣合成信号sti与st(i-1)满足终止条件||sti-st(i-1)||≤ξ,停止并输出最优解st=sti,h=hi;否则,转入步骤3.2c)中继续循环。
步骤4,在恒模条件约束下,设计发射波形矩阵。对已得到的发射信号自相关矩阵Rx,采用cyclicalgorithm算法设计发射信号波形矩阵XCA;利用得到的主瓣合成信号st,最终得到发射波形矩阵为:
Xopt=diag(ejφ)XCA
其中,φ=arg(st)-arg(XCAaT(θt)),arg(a)表示矢量a的相位角度(弧度)。
具体地,在恒模条件约束下,发射波形矩阵设计的数学模型写为:
s.t.XaT(θt)=st<64>
在每次发射过程中,如果给发射波形加任一相位,MIMO雷达的发射方向图将保持不变,即:
其中X表示改变相位后的波形矩阵,表示增加的相位向量,K表示将整个空域等分的个数,diag(·)表示对角矩阵。
利用这一性质,<64>式可以简化成在满足给定发射方向图条件下的发射波形矩阵设计问题,等价的数学模型为:
利用CA算法求解<66>式,可以得到初始发射波形矩阵XCA。
利用步骤3中得到的最优主瓣合成信号st,最终的发射波形矩阵可以通过改变初始发射波形矩阵XCA的相位得到,如<67>式所示:
Xopt=diag(ejφ)XCA<67>
其中φ=arg(st)-arg(XCAaT(θt)),arg(a)表示矢量a的相位角度(弧度)。
本发明的效果通过以下仿真对比试验进一步说明:
(1)实验场景:考虑收发共置的集中式MIMO雷达,其阵元数为16,阵元间距为半波长,发射波形的中心频率、带宽、脉冲持续时间分别为10GHz、4MHz和16μs。发射波形的采样频率为fs=4MHz,采样周期为Ts=0.25μs,码长为64,目标所在的方位角为0°。杂波存在的空间范围为[-60°,60°],且主要分布在两个方向上,并用7°角度对杂波空域进行离散采样,杂波响应向量长度为128。
(2)仿真内容:
本实验中的杂波信号响应如图3所示。通过<15>式计算空域上不同延迟路径的杂波反射系数的方差,得到杂波信号在空域上的强度分布如图4所示,0°目标方向上的主瓣杂波响应如图5所示。
实验中比较了以下三种不同波形对杂波的抑制性能:
波形1:利用本发明得到的在空间域和时间域上分别进行优化的发射波形矩阵(β=0.5dB,λ1=0.03);
波形2:利用CA算法从<39>式得到的发射波形矩阵(β=0.5dB);
波形3:未优化发射方向图的线性调频波形。
三种波形得到的发射方向图如图6所示,分别为通过本发明得到的波形1的发射方向图(最优方向图)、通过CA方法得到的波形2的发射方向图(CA方向图)和传统线性调频波形3的发射方向图(未优化的方向图)。
经过联合优化主瓣合成信号和失配滤波器后,波形1的主瓣合成信号的相位和功率谱密度分别如图7和图8所示。经过失配滤波器滤波后,波形1的主瓣合成信号的自相关图如图9所示,波形2的主瓣合成信号的自相关图如图10所示。
分别仿真了波形1、波形2和波形3的主瓣输出杂波,结果分别如图11-13所示。分别仿真了波形1、波形2和波形3的输出杂波,结果分别如图14-16所示。分别仿真了波形1与用于比较的波形2和波形3在不同主瓣增益损失情况下对输出信号能量与杂波能量比SCR的影响,仿真结果如图17所示。
(3)仿真结果分析:
图6中圆圈加横线表示本发明得到的波形1的发射方向图,星号加横线表示利用CA算法得到波形2的发射方向图,虚线表示传统发射波形3的方向图。图17中圆圈加横线表示波形1与波形2相比,主瓣增益损失对输出信号能量与杂波能量比SCR的影响,星号加横线表示波形1与波形3相比,主瓣增益损失对输出信号能量与杂波能量比SCR的影响。
从图6中可以看出,本发明得到的发射方向图在杂波区域非常低,可以看成是在空域上的杂波白化,从而有效地抑制杂波。这是由于本发明在联合优化主瓣合成信号和接收滤波器的设计中,考虑了降低自相关峰值旁瓣电平和杂波峰值电平。
从图9-12中可以看出,杂波峰值电平的降低是以提高自相关旁瓣电平为代价的。
从图11-16中可以看出,最终输出的主瓣杂波强于旁瓣杂波。定义输出的主瓣杂波峰值与信号之比(peakcluttertosignalratio,PCSR)为输出主瓣杂波的峰值能量与输出目标信号能量的比值,通过计算得到波形1的PCSR比波形2和波形3的PCSR低了15.11dB和16.57dB,并且波形1的输出杂波峰值功率与输出目标信号功率的比值比波形2和波形3的输出杂波峰值功率与输出目标信号功率的比值要低11.24dB和12.84dB。通过计算输出信号能量与杂波能量比SCR,得到波形1的输出SCR比波形2和波形3的SCR分别高4.49dB和7.25dB。
从图17中可以看出,输出信号能量与杂波能量比SCR的变化并不是随着主瓣增益损失的增加而增加的,从实验中得到一个合适的主瓣增益损失系数应选择在[0.2dB,0.5dB]之间。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (4)
1.一种杂波中检测目标用MIMO雷达发射波形的分级设计方法,其特征在于,包括以下步骤:
步骤1,建立MIMO雷达信号模型,获取杂波在空域上的强度以及在时间域上的主瓣响应;
步骤2,根据空域上杂波强度以及MIMO雷达信号模型,分别计算出目标回波能量Pt和θk方向上杂波的回波能量Pck,得到输出信号能量与杂波能量比SCR为其中,k=1,2,…,K,K表示杂波区域在空域上的等分个数;在自相关矩阵半正定以及发射天线发射能量恒定条件约束下,通过最大化输出信号能量与杂波能量比SCR得到最优的发射信号的自相关矩阵Rx;
步骤3,根据MIMO雷达的信号模型分别得到主瓣合成信号和主瓣杂波响应分别为st=XaT(θt)和qt=[qt(1),qt(2),…,qt(Nc)]T,其中θt表示目标所在的方位角,aT(θt)表示θt方向上的发射导向矢量,Np表示信号码长,N表示阵元个数,C表示复数域;发射信号波形矩阵X表示为xm=[xm(1),xm(2),…,xm(Np)]T,其中xm(n)(m=1,…,N)表示第m个阵元天线发射的恒模信号,其中d表示阵元间距,λ表示信号波长;利用联合优化主瓣合成信号和失配滤波器的设计得到最优的主瓣合成信号st以及失配滤波器h;
步骤4,在恒模条件约束下,设计发射波形矩阵:对已得到的发射信号自相关矩阵Rx,采用cyclicalgorithm算法设计发射信号波形矩阵XCA;利用得到的主瓣合成信号st,最终得到发射波形矩阵为:
Xopt=diag(ejφ)XCA
其中,φ=arg(st)-arg(XCAaT(θt)),arg(a)表示矢量a的相位角度。
2.根据权利要求1所述的杂波中检测目标用MIMO雷达发射波形的分级设计方法,其特征在于,步骤2的子步骤为:
2.0)根据空域上杂波强度以及MIMO雷达信号模型,目标回波的能量表示为:
其中Rx=XHX表示发射信号的自相关矩阵,tr(·)表示计算矩阵迹的运算,(·)H表示矩阵的共轭转置运算,||·||2表示Euclidean范数,表示目标的复反射系数,yt表示输出信号的信号分量;
方位角θk方向上的杂波回波能量可近似表示为:
其中表示方位角θk方向上的杂波强度,w是一个N×1维向量,ωk=wHaR(θk),yck表示方位角θk方向上的杂波回波信号;其中aR(θ)表示θ方向上的接收导向矢量,
发射信号的自相关矩阵Rx的设计如下:
2.1)写出目标函数:
输出信号能量与杂波能量比SCR如<3>式所示:
其中
2.2)最大化<3>式,建立如下凸规划数学模型:
s.t.Rx≥0<4>
Rx(i,i)=E,i=1,2,…,N
其中Rx≥0表示矩阵Rx为半正定矩阵,[M]ij表示矩阵M的第(i,j)个元素,E表示每根阵元天线的发射能量,N表示阵元个数;
将<4>式中的代价函数转化为最小化旁瓣杂波能量,数学模型等价于<5>式:
s.t.tr(RxRat)=N2E(1-β)<5>
Rx≥0
[Rx]ii=Efori=1,2,…,N
其中β是用来控制主瓣增益损失的一个小的正数,N表示阵元个数;
2.3)通过凸优化工具包cvx可以求得<5>式中最优的发射信号的自相关矩阵Rx。
3.根据权利要求1所述的杂波中检测目标用MIMO雷达发射波形的分级设计方法,其特征在于,步骤3的子步骤为:
3.0)首先,根据MIMO雷达信号模型,经过滤波器后输出的信号分量和杂波分量为<6>式:
其中矩阵H是(Np+Nh-1)×Np维的矩阵,Qk为(Np+NC+Nh-2)×(Np+Nh-1)维矩阵,Nh表示失配滤波器的长度,表示目标的复反射系数,分别如下所示:
其中,qk表示方位角θk上对应的杂波响应,然后,根据<6>式中zt的中心元素表示主瓣自相关函数,定义自相关主瓣增益损失为:
其中[M]i,:表示矩阵第i行所有元素;
自相关峰值旁瓣电平表示为:
其中max{·}表示选取所有元素中最大值的运算;
杂波峰值电平表示为:
CPL=max{|zc(i)|2},i=1,…,Np+Nh+Nc-2<11>
得到主瓣合成信号波形和主瓣杂波响应分别为:
st=XaT(θt)<12>
qt=「qt(1),qt(2),…,qt(Nc)]T<13>
得到主瓣杂波的回波表达式为:
ct=qt*st<14>
经接收失配滤波器h后,输出的主瓣杂波回波信号和目标回波信号分别为:
zct=h*qt*st=qt*h*st<15>
分别写成矩阵形式,输出的主瓣杂波回波信号和目标回波信号可以重新写为:
zct=QtHst=QtSth<17>
其中矩阵St为(Np+Nh-1)×Nh维矩阵,表示为:
Qt为由qt组成的矩阵,利用<8>式根据qt计算矩阵Qt;
3.1)在控制自相关主瓣增益损失情况下,最小化自相关旁瓣电平APSL和杂波峰值电平CPL,得到数学模型为:
|zct(i)|2≤p1,i=1,2,…,NC+NP+Nh-2<20>
其中λ1和λ2都是小的正数,分别用来控制自相关主瓣增益损失以及目标函数中杂波峰值电平和自相关旁瓣电平比例的正尺度因子,p1表示自相关峰值旁瓣电平,p2表示杂波峰值电平,β是用来控制主瓣增益损失的一个小的正数,
不失一般性,将主瓣合成信号的能量归一为单位能量,数学模型重新写为:
|zc(i)|2≤p1,i=1,2,…,Nc+Np+Nh-2<21>
3.2)对<21>式按如下步骤求解:
3.2a)按式<12>设置初始的主瓣合成信号矩阵st,按式<13>计算主瓣杂波响应向量qt;
3.2b)令i=0,设置Nh×1维初始失配滤波器hi为
其中表示将st中对第n个元素取共轭,同时设置小的正数λ1、λ2和ξ;
3.2c)令i值增1,利用<13>式计算矩阵H,同时固定hi-1,利用cvx优化工具包寻找<22>式的最优解
s.t.max|QtHst|≤p1<22>
max|Hst-d|≤p2
||st||2≤1
其中表示期望的自相关函数,求解得到最优的主瓣合成信号为其中arg(a)表示矢量a的相位角度;
3.2d)固定sti,通过<23>式优化接收滤波器hi:
max|QtSthi|≤p1
maxSthi-d|≤p2
其中矩阵St是步骤3.2c)中取得的sti通过<19>式计算得到的;
3.2e)重复步骤3.2c)和步骤3.2d),直到两次相邻两次循环得到的主瓣合成信号sti与st(i-1)满足终止条件||sti-st(i-1)||≤ξ,停止并输出最优解st=sti,h=hi;否则,转入步骤3.2c)中继续循环。
4.根据权利要求1所述的杂波中检测目标用MIMO雷达发射波形的分级设计方法,其特征在于,步骤4的子步骤为:
4.1)在恒模条件约束下,发射波形矩阵设计的数学模型写为:
s.t.XaT(θt)=st<24>
其中,Rx表示发射信号自相关矩阵,E表示每根阵元天线的发射能量;在每次发射过程中,如果给发射波形加任一相位,MIMO雷达的发射方向图将保持不变,即:
其中X表示改变相位后的波形矩阵,表示增加的相位向量,K表示将整个空域等分的个数,diag(·)表示对角矩阵;
4.2)利用这一性质,<24>式简化成在满足给定发射方向图条件下的发射波形矩阵设计问题,等价的数学模型为:
利用CA算法求解<26>式,可以得到初始发射波形矩阵XCA;
4.3)利用步骤3中得到的最优主瓣合成信号st,最终的发射波形矩阵可以通过初始发射波形矩阵XCA的相位得到,如<27>式所示:
Xopt=diag(ejφ)XCA<27>
其中φ=arg(st)-arg(XCAaT(θt)),arg(a)表示矢量a的相位角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410122330.8A CN103969633B (zh) | 2014-03-28 | 2014-03-28 | 杂波中检测目标用mimo雷达发射波形的分级设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410122330.8A CN103969633B (zh) | 2014-03-28 | 2014-03-28 | 杂波中检测目标用mimo雷达发射波形的分级设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103969633A CN103969633A (zh) | 2014-08-06 |
CN103969633B true CN103969633B (zh) | 2016-05-18 |
Family
ID=51239363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410122330.8A Active CN103969633B (zh) | 2014-03-28 | 2014-03-28 | 杂波中检测目标用mimo雷达发射波形的分级设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103969633B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102015203454A1 (de) * | 2014-10-07 | 2016-04-07 | Robert Bosch Gmbh | Verfahren und MIMO-Radarvorrichtung zum Bestimmen eines Lagewinkels eines Objekts |
CN105044681B (zh) * | 2015-06-05 | 2017-06-20 | 西安电子科技大学 | 高码率低距离分辨率相位编码信号的失配滤波器优化方法 |
CN105044682B (zh) * | 2015-06-05 | 2017-05-24 | 西安电子科技大学 | 宽主瓣相位编码信号和失配滤波器的联合优化方法 |
CN104865558B (zh) * | 2015-06-05 | 2017-09-01 | 西安电子科技大学 | 基于p范数的相位编码信号和失配滤波器的联合优化方法 |
CN105655727B (zh) * | 2015-09-28 | 2018-07-24 | 中国电子科技集团公司第二十九研究所 | 一种基于凸优化的增益损失恒定波束的形成方法及装置 |
CN106093877B (zh) * | 2016-07-19 | 2018-09-21 | 西安电子科技大学 | 正交宽主瓣相位编码信号和失配滤波器联合优化方法 |
CN106646420B (zh) * | 2016-09-28 | 2019-08-06 | 西安电子科技大学 | 基于lfm信号的mimo雷达发射方向图设计方法 |
CN106682405B (zh) * | 2016-12-14 | 2019-10-18 | 西北工业大学 | 基于凸优化的低旁瓣波束图综合设计方法 |
CN107622036B (zh) * | 2017-09-30 | 2020-07-21 | 中国人民解放军战略支援部队航天工程大学 | 一种基于蚁群优化的多项式相位信号自适应时频变换方法 |
EP3762742B1 (en) * | 2018-03-08 | 2023-11-15 | IEE International Electronics & Engineering S.A. | Method and system for target detection using mimo radar |
CN108490401B (zh) * | 2018-03-27 | 2020-11-13 | 北京理工大学 | 一种频率步进相位编码信号的旁瓣抑制方法 |
CN110531326B (zh) * | 2018-05-24 | 2023-06-30 | 中安锐达(南京)电子科技有限公司 | 低慢小雷达对地杂波抑制的发射波束控制算法 |
CN109061599B (zh) * | 2018-08-28 | 2022-06-03 | 电子科技大学 | 一种基于循环平稳和对称先验知识的stap方法 |
CN111025258B (zh) * | 2019-12-04 | 2022-02-08 | 北京理工大学 | 一种针对雷达波形分集的联合失配滤波器及其设计方法 |
CN112558021B (zh) * | 2020-11-18 | 2022-03-04 | 中国人民解放军国防科技大学 | 用于脉冲压缩雷达系统的波形及滤波器设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4580139A (en) * | 1983-06-22 | 1986-04-01 | The United States Of America As Represented By The Secretary Of The Army | Waveform design for optimized ambiguity response |
CN103076596A (zh) * | 2013-01-17 | 2013-05-01 | 西安电子科技大学 | 基于先验信息的mimo雷达发射方向图设计方法 |
-
2014
- 2014-03-28 CN CN201410122330.8A patent/CN103969633B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4580139A (en) * | 1983-06-22 | 1986-04-01 | The United States Of America As Represented By The Secretary Of The Army | Waveform design for optimized ambiguity response |
CN103076596A (zh) * | 2013-01-17 | 2013-05-01 | 西安电子科技大学 | 基于先验信息的mimo雷达发射方向图设计方法 |
Non-Patent Citations (4)
Title |
---|
MIMO雷达稳健的发射波束形成算法;王勇 等;《电子与信息学报》;20120229;第34卷(第2期);全文 * |
Minimax Robust MIMO Radar Waveform Design;Yang Yang 等;《IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING》;20070630;第1卷(第1期);全文 * |
基于凸优化的宽带雷达波形优化方法;纠博 等;《电波科学学报》;20090430;第24卷(第2期);全文 * |
针对方向图综合的MIMO雷达双边自适应矩阵算法;王勇 等;《电子与信息学报》;20120430;第34卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103969633A (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103969633B (zh) | 杂波中检测目标用mimo雷达发射波形的分级设计方法 | |
CN103076596B (zh) | 基于先验信息的mimo雷达发射方向图设计方法 | |
CN102156279B (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
CN101251597B (zh) | 一种多输入多输出雷达系统阵列误差自校正的方法 | |
CN103901417B (zh) | L型阵列mimo雷达低复杂度空间目标二维角度估计方法 | |
CN103353591B (zh) | 基于mimo的双基地雷达局域化降维杂波抑制方法 | |
CN105699953B (zh) | 频率分集mimo雷达距离‑角度解耦合波束形成方法 | |
CN103592642B (zh) | Mimo雷达波形的设计方法 | |
CN105467365A (zh) | 一种改善mimo雷达doa估计性能的低旁瓣发射方向图设计方法 | |
CN103353592B (zh) | 基于mimo的双基地雷达多通道联合降维杂波抑制方法 | |
CN103592647B (zh) | 阵列三维sar数据获取方法 | |
CN107390208A (zh) | 一种频率分集阵列雷达通信一体化方法 | |
CN103983952A (zh) | 一种非圆信号双基地mimo雷达低复杂度收发角度联合估计方法 | |
CN106646394A (zh) | 一种fda‑mimo雷达中基于循环优化的发射波束矩阵设计方法 | |
CN106646387A (zh) | 基于发射波束域的mimo雷达抗有源干扰方法 | |
CN102707264A (zh) | 基于圆形阵列双基地mimo雷达的波达方向估计方法 | |
CN103018732B (zh) | 基于空时联合优化的mimo雷达波形合成方法 | |
CN104833959A (zh) | 基于目标先验信息的mimo雷达波形优化方法 | |
CN102841335B (zh) | 基于迭代fft的mimo雷达快速波形合成方法 | |
CN104808179A (zh) | 杂波环境下基于克拉美罗界的mimo雷达波形优化方法 | |
CN101887117A (zh) | 基于三迭代的机载mimo雷达空时降维自适应处理方法 | |
CN106772337A (zh) | 一种基于music算法的频率分集阵列两层介质目标定位方法 | |
CN103364762B (zh) | 任意阵列流形的单基地mimo雷达波达方向估计方法 | |
CN109828252A (zh) | 一种mimo雷达参数估计方法 | |
CN110471033A (zh) | 基于相干fda的二维空时编码干扰抑制方法 |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Jiu Bo Inventor after: Chen Yanfei Inventor after: Liu Hongwei Inventor after: Wang Yinghua Inventor after: Dai Fengzhou Inventor before: Jiu Bo Inventor before: Chen Yanfei Inventor before: Liu Hongwei Inventor before: Wang Yinghua Inventor before: Dai Fengzhou |