CN111856410B - 恒模约束下机载mimo雷达的稳健发射接收联合优化方法 - Google Patents
恒模约束下机载mimo雷达的稳健发射接收联合优化方法 Download PDFInfo
- Publication number
- CN111856410B CN111856410B CN202010705704.4A CN202010705704A CN111856410B CN 111856410 B CN111856410 B CN 111856410B CN 202010705704 A CN202010705704 A CN 202010705704A CN 111856410 B CN111856410 B CN 111856410B
- Authority
- CN
- China
- Prior art keywords
- waveform
- transmitting
- receiving
- target
- constant modulus
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000005457 optimization Methods 0.000 title claims abstract description 74
- 238000013461 design Methods 0.000 claims abstract description 44
- 238000012545 processing Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 48
- 230000001427 coherent effect Effects 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000002592 echocardiography Methods 0.000 claims description 5
- 241000764238 Isis Species 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000004891 communication Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 15
- 230000004044 response Effects 0.000 description 11
- 230000001419 dependent effect Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000005428 wave function 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
- 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
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
-
- 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/285—Receivers
- G01S7/292—Extracting wanted echo-signals
-
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
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雷达的稳健发射接收联合优化方法,包括以下步骤:机载MIMO雷达快时间波形相关信号模型的建立;稳健的发射接收联合设计模型的建立;接收滤波器的优化;基于分式规划和Power method‑like迭代,或者基于半正定规划和随机化的稳健恒模发射波形优化;利用优化结果进行STAP处理。本发明考虑目标多普勒频率和空间锥角存在着不确定性,以平均输出SCNR作为优化目标,即最大化平均输出SCNR,同时保持发射波形的恒模约束,提出了基于分式规划和Power method‑like迭代,以及基于半正定规划和随机化的稳健恒模发射波形优化方法,这两种方法都能够获得迭代递增的输出SCNR性能,而且对目标参数的不确定性具有稳健性,可用于雷达和通信等领域。
Description
技术领域
本发明涉及MIMO雷达技术领域,尤其涉及一种恒模约束下机载 MIMO雷达的稳健发射接收联合优化方法。
背景技术
多输入多输出(Multiple-input-multiple-output,MIMO)雷达可以通过天线发射可自由控制的多个脉冲信号,使得MIMO雷达比一般的相控阵雷达具有更好的性能,从而提升了雷达系统设计的灵活性。根据收发天线配置的方式,MIMO雷达可以分为统计MIMO雷达和相干MIMO雷达。 MIMO雷达由于其能够发射灵活的波形得到了广泛的应用,如统计MIMO雷达利用空间分集特性来提升目标的检测性,相干MIMO雷达具有更好的参数辨识能力和更灵活的发射方向图设计能力。
由于MIMO雷达系统具有上述优点,学者们将其应用到了机载雷达 STAP系统中。与传统的单输入多输出(single-input-multiple-output,SIMO) STAP雷达相比,MIMO-STAP具有更高的多普勒分辨性,能够获得更窄的杂波凹口。
近年来,MIMO雷达波形设计得到了学者的广泛关注,这些工作可以分为两类:第一类主要利用一些先验信息来解决快时间波形设计问题。文献《Waveform design for MIMOradars》(B.Friedlander,IEEE Trans.Aerosp. Electron.Syst.,vol.43,no.3,pp.1227–1238,Jul.2007.)采用梯度算法来设计能量约束的针对扩展目标的波形优化问题。文献《MIMO radar space time adaptive processing using prolate spheroidal wavefunctions》(C.-Y.Chen and P.P.Vaidyanathan,IEEE Trans.Signal Process.,vol.56,no.2,pp.623–635, 2008.)设计了能量约束的发射波形和接收滤波器优化问题来解决信号依赖杂波下的扩展目标检测问题。文献《多输入多输出雷达恒模稳健波形联合优化算法[J]》(崔琛,孙从易,郝天铎,等。国防科技大学学报,2018,40 (4):166-172.)针对目标先验信息不准确的问题,建立了极大极小化的恒模约束下的波形和接收滤波器联合优化模型,并采用了迭代优化步骤进行求解。文献《MIMO radar waveform design with constantmodulus and similarity constraints》(G.Cui,H.Li,and M.Rangaswamy,IEEETrans.Signal Process., vol.62,no.2,pp.343–353,Jan.2014.)提出了两种序贯优化算法来解决发射波形和接收滤波器联合优化的问题,考虑了恒模和相似度约束。
第二类主要是MIMO雷达慢时间编码设计来进行运动目标检测。文献《Beamspaceslow-time MIMO radar for multipath clutter mitigation》(V.F. Mecca,J.L.Krolik,and F.C.Robey,“Beamspace slow-time MIMO radar for multipath cluttermitigation,”in Proc.IEEE Int.Conf.Acoust.Speech,Signal Process.,Mar.2008,pp.2313–2316.)提出了空时编码框架来解决信号依赖干扰下基于波束空间的STAP问题。文献《Time-division beamforming for MIMO radar waveform design》(A.J.Duly andJ.V.Krogmeier, “Time-division beamforming for MIMO radar waveform design,”IEEE Trans. Aerosp.Electron.Syst.,vol.49,no.2,pp.1210–1223,Apr.2015.)采用最大化输出信干噪比(SINR)来设计了针对多目标场景的时分波束形成的信号。文献《信号相关杂波背景下稳健的恒模序列与接收滤波器设计方法[J]》 (付月,崔国龙,余显祥.雷达学报,2017,6(3):292-299.)研究了一种稳健的慢时间发射波形和接收滤波器联合优化设计的方法,该方法具有对抗不确定杂波信息的能力,而且计算复杂度较低。针对MIMO-STAP雷达发射波形和接收滤波器联合优化问题,文献《Joint design of transmit waveformsand receive filter for MIMO radar space-time adaptive processing [J]》(TangB.,Tang J.IEEE Trans.Signal Process.,2016,64(18):4707-4722.) 针对能量约束、恒模约束和相似度约束提出了多种优化算法,然而并没有考虑目标先验信息的不确定性。现有技术对机载MIMO雷达的研究均是对理想条件下的机载MIMO雷达发射波形和接收滤波器联合设计,无法解决先验信息不准确下的机载MIMO雷达稳健联合设计问题。
发明内容
针对上述存在的问题,本发明旨在提供一种恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,考虑目标多普勒频率和空间锥角存在的不确定性,以平均输出SCNR作为优化目标,即以最大化平均输出SCNR 作为联合设计的目标函数,考虑发射波形的恒模约束,利用迭代优化方案来迭代求解发射波形和接收滤波器,从而实现目标参数不确定性时的稳健性设计。
为了实现上述目的,本发明所采用的技术方案如下:
恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,其特征在于,包括以下步骤,
S1:构建机载MIMO雷达快时间波形相关信号模型;
S2:基于步骤S1中的快时间波形相关信号模型,构建具有恒模约束的稳健发射接收联合设计模型;
S3:对稳健发射接收联合设计模型中的接收滤波器进行优化;
S4:对稳健发射接收联合设计模型中的发射波形进行优化。
进一步的,步骤S1中构建机载MIMO雷达快时间波形相关信号模型的方法包括以下步骤,
S11:对快时间变化的机载MIMO雷达系统,假设机载雷达的发射阵元和接收阵元个数分别为NT和NR,dT和dR分别表示发射阵元之间和接收阵元之间的间距,它们之间满足如下关系dT=γdR,每个发射阵元在一个相干处理间隔内发射K个慢时间编码的脉冲,机载雷达平台的速度为Vp,目标相对于平台的径向速度为Vt,目标相对于平台飞行方向的空间锥角为φt;
其中,(·)T表示矩阵转置操作,ξt表示目标的散射幅度, fd,t=2[Vpcos(φt)+Vt]/(λfr)表示目标的归一化多普勒频率,λ表示波长,fr表示脉冲重复频率,和分别表示发射空域导向矢量和接收空域导向矢量,它们的具体表示形式为和其中, fs,t=dTcos(φt)/λ表示目标的空域频率;
S13:对Xt,k进行向量化操作,可以得到
S15:机载MIMO雷达系统中,杂波回波可以看做是信号相关的干扰,而且单个距离单元的杂波回波通常由许多相互独立且是信号相关的杂波散射块构成;令接收到的待检测单元的杂波回波为
其中,ξc,i表示第i个杂波块的复幅度,表示相应的归一化多普勒频率,vc,i表示目标的空时导向矢量,具体表达形式为其中,b(φc,i)和a(φc,i)分别表示第i个杂波块对应的时域导向矢量,接收导向矢量和发射导向矢量,它们具有与目标相同的表达形式;
S16:将uc,i,bc,i和ac,i分别表示b(φc,i)和a(φc,i),则机载MIMO雷达系统总的接收信号为x=xt+xc+xn,其中,xn表示KNRLc×1维的噪声向量,服从零均值且协方差矩阵为的复高斯分布。
进一步的,步骤S2构建具有恒模约束的稳健发射接收联合设计模型的具体操作包括以下步骤,
S24:当目标参数不确定时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型可以表示为
进一步的,步骤S23中用平均输出SCNR来解决目标参数信息的不确定性,可表示为
和
其中,(n1,n2)∈{1,2,…,NR}2,(q1,q2)∈{1,2,…,NT}2。
进一步的,步骤S3的具体操作为:固定波形矩阵S,获取最优接收滤波器w。
进一步的,步骤S4的具体操作为:固定接收滤波器w,利用分式规划和Powermethod-like迭代获取最优波形矩阵S。
进一步的,固定接收滤波器w,利用分式规划和Power method-like迭代获取最优波形矩阵S的具体操作包括,
S41:固定接收滤波器w时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型的优化问题可以表示为
S42:令s(n,r)为步骤S41中模型的第(n,r)次迭代时的解,对于固定的 w(n),f(n,r)表示为由s(n,r)计算得到的步骤S41中模型的目标函数值;其中, n表示全局外部迭代次数,r表示内部波形优化迭代次数;
其中,Q(n,r)=Ξt(W(n))-f(n,r)Rdl(W(n));
S45:将步骤S44中的模型优化问题,等价转化为
进一步的,步骤S4的具体操作为:固定接收滤波器w,基于半正定规划和随机化的方法来获取最优波形矩阵S。
进一步的,固定接收滤波器w,基于半正定规划和随机化的方法来获取最优波形矩阵S的具体操作包括,
S41:固定接收滤波器w时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型的优化问题可以表示为
S42:将步骤S41中的模型优化问题转换成等价形式
S43:将步骤S41中对波形s的优化问题转化为对波形协方差矩阵Rs的优化,从而将步骤S41表示成基于SDP的优化问题,即
s.t. Diag(Rs)=ps·1,
rank(Rs)=1,Rs≥0;
其中,rank(·)表示矩阵的秩;
S44:去掉秩一约束,步骤S43中的优化问题可以转化为
s.t. Diag(Rs)=ps·1,Rs≥0;
S45:采用Charnes-Cooper变换求解步骤S44中的优化问题,则
s.t. tr[RDL(W)Rs]=1
本发明的有益效果是:
本发明公开了一种目标参数不确定情况下恒模约束时机载多输入多输出(Multiple-input-multiple-output,MIMO)雷达稳健的发射接收联合优化方法,解决了目标参数先验信息不准确时发射波形和接收滤波器的联合设计问题。考虑目标多普勒频率和空间锥角存在着不确定性,以平均输出 SCNR作为优化目标,即最大化平均输出SCNR,同时保持发射波形的恒模约束,提出了基于分式规划和Power method-like迭代,或者基于半正定规划和随机化的稳健恒模发射波形优化方法,该发明方法能够获得迭代递增的输出SCNR性能,而且对目标参数的不确定性具有稳健性,可用于雷达和通信等领域。
附图说明
图1为本发明机载MIMO-STAP雷达几何结构示意图;
图2a为本发明目标参数确定性和随机性情形下PML迭代方法的输出 SCNR随着迭代次数的变化情况曲线图;
图2b为本发明目标参数确定性和随机性情形下SDP随机化方法的输出SCNR随着迭代次数的变化情况曲线图;
图3a为本发明目标参数确定情形下的SDP随机化方法的空时响应图;
图3b为本发明目标参数不确定情形下的SDP随机化方法的空时响应图;
图4a为本发明目标参数确定情形下的PML迭代方法的空时响应图;
图4b为本发明目标参数不确定情形下的PML迭代方法的空时响应图;
图5a为本发明SDP随机化方法和PML迭代方法输出SCNR随着目标多普勒不确定性的变化关系示意图;
图5b为本发明SDP随机化方法和PML迭代方法输出SCNR随着目标空域频率不确定性的变化关系示意图。
具体实施方式
为了使本领域的普通技术人员能更好的理解本发明的技术方案,下面结合附图和实施例对本发明的技术方案做进一步的描述。
实施例一:
恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,包括以下步骤,
S1:构建机载MIMO雷达波形信号模型;
具体的,考虑正侧视机载MIMO雷达系统,其几何结构示意图如附图 1所示,对快时间变化的机载MIMO雷达系统,假设机载雷达的发射阵元和接收阵元个数分别为NT和NR,dT和dR分别表示发射阵元之间和接收阵元之间的间距,它们之间满足如下关系dT=γdR,每个发射阵元在一个相干处理间隔内发射K个慢时间编码的脉冲,机载雷达平台的速度为Vp,目标相对于平台的径向速度为Vt,目标相对于平台飞行方向的空间锥角为φt;
其中,(·)T表示矩阵转置操作,ξt表示目标的散射幅度, fd,t=2[Vpcos(φt)+Vt]/(λfr)表示目标的归一化多普勒频率,λ表示波长,fr表示脉冲重复频率,和分别表示发射空域导向矢量和接收空域导向矢量,它们的具体表示形式为和其中, fs,t=dTcos(φt)/λ表示目标的空域频率;
对Xt,k进行向量化操作,可以得到
对于波形相关的杂波回波,不是一般性,而且单个距离单元的杂波回波通常由许多相互独立且是信号相关的杂波散射块构成;考虑待检测单元的杂波信号,则接收到的待检测单元的杂波回波为
其中,ξc,i表示第i个杂波块的复幅度,表示相应的归一化多普勒频率,vc,i表示目标的空时导向矢量,具体表达形式为其中,b(φc,i)和a(φc,i)分别表示第i个杂波块对应的时域导向矢量,接收导向矢量和发射导向矢量,它们具有与目标相同的表达形式;
进一步的,步骤S2:基于步骤S1中的波形信号模型,构建具有恒模约束的稳健发射接收联合设计模型;
具体的,对步骤S1中接收到的目标信号进行滤波处理得到对S1中接收到的杂波和噪声信号进行处理得到最大化输出信杂噪比,得到STAP处理后的SCNR,令表示接收机载雷达 STAP系统的权矢量,则经过处理后输出SCNR为
其中,Rc表示波形相关的杂波协方差矩阵,具体形式为
则式(6)中输出SCNR可表示为下面的等价形式
时,采用平均输出SCNR来解决目标参数信息的不确定性,也即
和
其中,Γ(vt)可通过下式来求得
其中,Ψ可以写成块矩阵的形式
当目标参数不确定时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型可以表示为
进一步,步骤S3:固定波形矩阵S,对稳健发射接收联合设计模型中的接收滤波器进行优化,获取最优接收滤波器w。
具体的,固定波形矩阵S时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型可以表示为
问题(24)属于经典的广义瑞利商问题,其最优解wopt可以表示为
进一步的,步骤S4:固定接收滤波器w,利用分式规划和Power method-like迭代对稳健发射接收联合设计模型中的发射波形进行优化,获取最优波形矩阵S。
具体的,固定接收滤波器w时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型的优化问题可以表示为
令s(n,r)表示优化问题(26)中第(n,r)次迭代时的解,对于固定的 w(n),f(n,r)表示为由s(n,r)计算得到的步骤S41中模型的目标函数值;其中, n表示全局外部迭代次数,r表示内部波形优化迭代次数;
因此,第(n,r+1)次迭代时模型的优化问题转化为
其中,Q(n,r)=Ξt(W(n))-f(n,r)Rdl(W(n));
基于Power method-like迭代方法,可通过求解下面的问题来迭代更新波形向量S
则问题(29)的解为
该解(30)也即稳健发射接收联合设计模型中最优发射波形S。
实施例二:
实施例二与实施例一的区别之处在于,步骤S4中固定接收滤波器w,基于半正定规划和随机化的方法来获取最优波形矩阵S,其他操作均相同。
固定接收滤波器w,基于半正定规划和随机化的方法来获取最优波形矩阵S的具体操作为:
固定接收滤波器w时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型的优化问题可以表示为
将(31)中的模型优化问题转换成等价形式
将波形s的优化问题转化为对波形协方差矩阵Rs的优化,从而将步骤模型(32)表示成基于SDP的优化问题,即
s.t. Diag(Rs)=ps·1,
rank(Rs)=1,Rs≥0 (33)
其中,rank(·)表示矩阵的秩;
去掉秩一约束,基于SDP的优化问题可以转化为
s.t. Diag(Rs)=ps·1,Rs≥0 (34)
采用Charnes-Cooper变换求解(34)中的优化问题,则
s.t. tr[RDL(W)Rs]=1
仿真试验:
通过仿真试验来证明本发明的稳健联合优化方法的有效性,机载MIMO雷达的主要参数如下:发射阵元和接收阵元个数分别为NT=6和 NR=6,脉冲个数为K=12,Nc=360,采用正交线性调频波形作为初始化波形和参考波形,则波形矩阵的第(nt,k)个输入为
其中,nt=1,…,NT,k=1,…,K,参考波形向量为s0=vec(S0)。
进一步的,对输出SCNR的对比进行仿真,假设目标的归一化多普勒频率和空域频率的不确定性大小分别为κt=0.03和θt=0.05,目标在角度-多普勒平面的位置为(-0.2,0.1)。将实施例一中基于Power method-like迭代方法记为PML迭代方法,将实施例二中基于半正定规划和随机化的方法记为SDP随机化方法。图2a和图2b分别给出了目标参数确定性和随机性情形下PML迭代方法和SDP随机化方法的输出SCNR随着迭代次数的变化情况。
从附图2a和2b中可以看出,本发明中的输出SCNR在目标参数确定性情形优于目标参数随机性情形,两种算法在两种情况下最优输出SCNR 的差约为2.95dB。
进一步的,对空时响应图对比进行仿真,附图3a为目标参数确定情形下的SDP随机化方法的空时响应图,附图3b为目标参数不确定情形下的SDP随机化方法的空时响应图。附图4a为目标参数确定情形下的PML 迭代方法的空时响应图,附图4b为目标参数不确定情形下的PML迭代方法的空时响应图。
对比附图3a、3b和附图4a、4b可以看出,两种情况下两种方法的联合优化方法均能够有效地抑制杂波,而且空时响应图的主瓣均位于(-0.2, 0.1)处。本发明中的联合优化方法空时响应图的主瓣在目标参数不确定情形下要比目标参数确定性情形下要略大,这意味着所提方法在目标参数不确定时具有更好的稳健性。
进一步的,对输出SCNR随目标参数不确定的变化进行仿真,附图 5a为SDP随机化方法和PML迭代方法输出SCNR随着目标多普勒不确定性的变化关系,其中θt=0;附图5b为SDP随机化方法和PML迭代方法输出SCNR随着目标空域频率不确定性的变化关系,其中κt=0。
对比附图5a和附图5b可以看出,整体上不同算法的输出SCNR随着目标多普勒频率和空域频率不确定性的增加而下降,尽管性能有所下降,但是这也反映了本发明中的联合优化方法对目标参数不确定性所表现出的稳健性。
进一步的,下表1给出了确定性场景下SDP随机化方法和PML迭代方法的最优输出SCNR值,达到收敛时需要的迭代次数和计算时间。
表1不同算法的迭代次数和计算时间
从表1可以看出,两种方法的最优SCNR值相同,尽管PML迭代方法比SDP随机化方法需要更多的迭代次数,但是PML迭代方法总的计算时间比SDP随机化方法显著低很多,因此,在具有相同输出SCNR性能的情况下,PML迭代方法在计算效率方面具有显著优势。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (8)
1.恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,其特征在于,包括以下步骤,
S1:构建机载MIMO雷达快时间波形相关信号模型;
S2:基于步骤S1中的快时间波形相关信号模型,构建具有恒模约束的稳健发射接收联合设计模型;
S3:对稳健发射接收联合设计模型中的接收滤波器进行优化;
S4:对稳健发射接收联合设计模型中的发射波形进行优化;
步骤S1中构建机载MIMO雷达快时间波形相关信号模型的方法包括以下步骤,
S11:对快时间变化的机载MIMO雷达系统,假设机载雷达的发射阵元和接收阵元个数分别为NT和NR,dT和dR分别表示发射阵元之间和接收阵元之间的间距,它们之间满足如下关系dT=γdR,每个发射阵元在一个相干处理间隔内发射K个慢时间编码的脉冲,机载雷达平台的速度为Vp,目标相对于平台的径向速度为Vt,目标相对于平台飞行方向的空间锥角为φt;
其中,(·)T表示矩阵转置操作,ξt表示目标的散射幅度,fd,t=2[Vpcos(φt)+Vt]/(λfr)表示目标的归一化多普勒频率,λ表示波长,fr表示脉冲重复频率,和分别表示发射空域导向矢量和接收空域导向矢量,它们的具体表示形式为和其中,fs,t=dTcos(φt)/λ表示目标的空域频率;
S13:对Xt,k进行向量化操作,得到
S15:机载MIMO雷达系统中,杂波回波看做是信号相关的干扰,而且单个距离单元的杂波回波通常由许多相互独立且是信号相关的杂波散射块构成;令接收到的待检测单元的杂波回波为
其中,ξc,i表示第i个杂波块的复幅度,表示相应的归一化多普勒频率,vc,i表示目标的空时导向矢量,具体表达形式为其中,b(φc,i)和a(φc,i)分别表示第i个杂波块对应的时域导向矢量,接收导向矢量和发射导向矢量,它们具有与目标相同的表达形式;
S16:将uc,i,bc,i和ac,i分别表示b(φc,i)和a(φc,i),则机载MIMO雷达系统总的接收信号为x=xt+xc+xn,其中,xn表示KNRLc×1维的噪声向量,服从零均值且协方差矩阵为的复高斯分布;
步骤S2构建具有恒模约束的稳健发射接收联合设计模型的具体操作包括以下步骤,
S21:对步骤S1中接收到的目标信号进行滤波处理得到E{|wHxt|2},对S1中接收到的杂波和噪声信号进行处理得到最大化输出信杂噪比,得到STAP处理后的SCNR,其中,Rc表示波形相关的杂波协方差矩阵,具体为
S24:当目标参数不确定时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型表示为
3.根据权利要求2所述的恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,其特征在于,步骤S3的具体操作为:固定波形矩阵S,获取最优接收滤波器w。
5.根据权利要求3所述的恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,其特征在于,步骤S4的具体操作为:固定接收滤波器w,利用分式规划和Power method-like迭代获取最优波形矩阵S。
6.根据权利要求5所述的恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,其特征在于,固定接收滤波器w,利用分式规划和Power method-like迭代获取最优波形矩阵S的具体操作包括:
S41:固定接收滤波器w时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型的优化问题表示为
S42:令s(n,r)为步骤S41中模型的第(n,r)次迭代时的解,对于固定的w(n),f(n,r)表示为由s(n,r)计算得到的步骤S41中模型的目标函数值;其中,n表示全局外部迭代次数,r表示内部波形优化迭代次数;
其中,Q(n,r)=Ξt(W(n))-f(n,r)Rdl(W(n));
S45:将步骤S44中的模型优化问题,等价转化为
7.根据权利要求3所述的恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,其特征在于,步骤S4的具体操作为:固定接收滤波器w,基于半正定规划和随机化的方法来获取最优波形矩阵S。
8.根据权利要求7所述的恒模约束下机载MIMO雷达的稳健发射接收联合优化方法,其特征在于,固定接收滤波器w,基于半正定规划和随机化的方法来获取最优波形矩阵S的具体操作包括:
S41:固定接收滤波器w时,恒模约束下MIMO-STAP雷达快时间下的发射波形和接收滤波器稳健发射接收联合设计模型的优化问题表示为
S42:将步骤S41中的模型优化问题转换成等价形式
S43:将步骤S41中对波形s的优化问题转化为对波形协方差矩阵Rs的优化,从而将步骤S41表示成基于SDP的优化问题,即
s.t.Diag(Rs)=ps·1,
rank(Rs)=1,Rs≥0;
其中,rank(·)表示矩阵的秩;
S44:去掉秩一约束,步骤S43中的优化问题转化为
s.t.Diag(Rs)=ps·1,Rs≥0;
S45:采用Charnes-Cooper变换求解步骤S44中的优化问题,则
s.t.tr[RDL(W)Rs]=1
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010705704.4A CN111856410B (zh) | 2020-07-21 | 2020-07-21 | 恒模约束下机载mimo雷达的稳健发射接收联合优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010705704.4A CN111856410B (zh) | 2020-07-21 | 2020-07-21 | 恒模约束下机载mimo雷达的稳健发射接收联合优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111856410A CN111856410A (zh) | 2020-10-30 |
CN111856410B true CN111856410B (zh) | 2022-08-05 |
Family
ID=73000746
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010705704.4A Expired - Fee Related CN111856410B (zh) | 2020-07-21 | 2020-07-21 | 恒模约束下机载mimo雷达的稳健发射接收联合优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111856410B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112558021B (zh) * | 2020-11-18 | 2022-03-04 | 中国人民解放军国防科技大学 | 用于脉冲压缩雷达系统的波形及滤波器设计方法 |
CN113960583B (zh) * | 2021-09-24 | 2024-08-13 | 中国人民解放军国防科技大学 | 基于发射波束域的机载mimo雷达发射接收稳健联合优化方法 |
CN114675238B (zh) * | 2022-02-24 | 2023-11-03 | 中国人民解放军国防科技大学 | 一种雷达通信一体化波形直接优化方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105807275A (zh) * | 2016-04-28 | 2016-07-27 | 大连大学 | 基于部分杂波先验知识的mimo-ofdm-stap稳健波形设计 |
CN107831480A (zh) * | 2017-10-13 | 2018-03-23 | 西安电子科技大学 | 弹载雷达和差通道稳健自适应杂波抑制方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008157646A2 (en) * | 2007-06-18 | 2008-12-24 | University Of Florida Research Foundation , Inc. | Mimo transmit beamforming under uniform elemental peak power constant |
CN101907702A (zh) * | 2010-06-24 | 2010-12-08 | 西安电子科技大学 | 用于mimo雷达的两维多脉冲对消器 |
CN105974391B (zh) * | 2016-04-28 | 2018-09-25 | 大连大学 | 非确知目标先验知识条件下mimo雷达稳健波形设计方法 |
CN107656257B (zh) * | 2017-09-28 | 2021-03-23 | 西安电子科技大学 | 一种弹载mimo雷达波形协方差矩阵的优化设计方法 |
CN108469601B (zh) * | 2018-01-17 | 2021-11-30 | 大连大学 | 改善mimo-stap最差检测性能的收发联合稳健优化方法 |
CN110456334B (zh) * | 2019-07-27 | 2023-04-28 | 南京理工大学 | 基于优化稀疏阵的tdm-mimo雷达系统及其信号处理方法 |
CN111025238B (zh) * | 2019-12-20 | 2022-07-22 | 西北工业大学 | 空间合成信号特性可控的mimo雷达波形设计方法 |
-
2020
- 2020-07-21 CN CN202010705704.4A patent/CN111856410B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105807275A (zh) * | 2016-04-28 | 2016-07-27 | 大连大学 | 基于部分杂波先验知识的mimo-ofdm-stap稳健波形设计 |
CN107831480A (zh) * | 2017-10-13 | 2018-03-23 | 西安电子科技大学 | 弹载雷达和差通道稳健自适应杂波抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111856410A (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111856410B (zh) | 恒模约束下机载mimo雷达的稳健发射接收联合优化方法 | |
CN105699945B (zh) | 频控阵mimo雷达系统中的波形优化设计方法 | |
CN101770022B (zh) | 基于遗传算法的mimo雷达阵列位置误差自校正方法 | |
CN104251989B (zh) | 基于压缩空间谱的单基地mimo雷达目标波达方向估计方法 | |
Boyer | Performance bounds and angular resolution limit for the moving colocated MIMO radar | |
CN105717496B (zh) | 一种基于矩阵填充的频控阵mimo雷达系统的实现方法 | |
Yan et al. | Two-stage ESPRIT for unambiguous angle and range estimation in FDA-MIMO radar | |
CN108303683B (zh) | 单基地mimo雷达实值esprit非圆信号角度估计方法 | |
CN103901395B (zh) | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 | |
CN105319545B (zh) | 提高stap检测性能的mimo-ofdm雷达波形设计方法 | |
CN103983958A (zh) | 基于多测量矢量稀疏表示的mimo雷达连续目标角度估计方法 | |
CN109738856B (zh) | 一种基于时间反演的噪声抑制Capon有源目标DOA估算方法 | |
Yu et al. | Compressive sensing for MIMO radar | |
CN105487054B (zh) | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 | |
CN105182313A (zh) | 一种基于不完全杂波先验知识的mimo-stap稳健波形设计方法 | |
CN103513250A (zh) | 一种基于鲁棒自适应波束形成原理的模基定位方法及系统 | |
CN108594194A (zh) | 基于四线性分解的双基地mimo雷达角度估算方法 | |
CN104730517A (zh) | 双基地多输入多输出雷达多目标跟踪方法 | |
CN113376607A (zh) | 机载分布式雷达小样本空时自适应处理方法 | |
Zhai et al. | Weak target detection in massive mimo radar via an improved reinforcement learning approach | |
Liu et al. | $\ell _ {0} $-Norm Minimization-Based Robust Matrix Completion Approach for MIMO Radar Target Localization | |
Wu et al. | Waveform design for higher-resolution localization with MIMO radar | |
CN109901103A (zh) | 基于非正交波形的mimo雷达doa估算方法及设备 | |
Hassanien et al. | Direction finding for MIMO radar with colocated antennas using transmit beamspace preprocessing | |
Pandey et al. | Space Time Adaptive Processing for High Signal to Clutter Ratio Target's Detection Radar System |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220805 |
|
CF01 | Termination of patent right due to non-payment of annual fee |