CN105487054B - 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 - Google Patents
提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 Download PDFInfo
- Publication number
- CN105487054B CN105487054B CN201510756071.9A CN201510756071A CN105487054B CN 105487054 B CN105487054 B CN 105487054B CN 201510756071 A CN201510756071 A CN 201510756071A CN 105487054 B CN105487054 B CN 105487054B
- Authority
- CN
- China
- Prior art keywords
- clutter
- matrix
- stap
- mimo
- waveform
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000013461 design Methods 0.000 title abstract description 11
- 238000005457 optimization Methods 0.000 claims abstract description 45
- 238000011068 loading method Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000009795 derivation Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 39
- 239000013598 vector Substances 0.000 claims description 24
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000011002 quantification Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000002592 echocardiography Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- 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
-
- 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
- G01S2013/0236—Special technical features
- G01S2013/0281—LPI, Low Probability of Intercept radar
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于信号处理领域,涉及一种提高基于MIMO‑OFDM雷达STAP最差检测性能的稳健波形设计方法,并提出通过将参数不确定性显式地包含进波形优化问题以此来改善基于MIMO‑OFDM雷达的STAP最差检测性能的稳健波形优化问题。其实现步骤包括:(1)建立MIMO‑OFDM‑STAP模型以获得接收单元的数据表达式;(2)通过对目标函数的推导,得出最优输出SINR的表达式;(3)基于稳健波形优化模型,得到最大化最差情况下的输出SINR;(4)提出基于对角加载(DL)方法求解稳健波形优化问题。
Description
技术领域
本发明属于信号处理领域,涉及一种提高基于MIMO-OFDM雷达STAP最差检测性能的稳健波形设计方法。可改善MIMO-OFDM雷达***性能的波形设计对初始参数估计误差较为敏感的问题,可显著提高最差初始参数估计下的***稳健检测性能。
背景技术
随着MIMO通信蓬勃的发展,以及雷达为突破自身限制对新理论以及新技术的需求,MIMO雷达概念应运而生。与只能发送相干波形的相控阵雷达相比,MIMO雷达可以利用多个发射单元发射几乎任意波形。基于阵列天线间距,MIMO雷达***可分为以下两类:(1)分置雷达,(2)共置雷达。前者采用分置较远的收发单元发射所需信号,同时从不同角度观察目标,从而可利用空间分集以克服由于目标闪烁造成的性能下降。相反,后者使用距离很近的发射单元以增加接收阵列的虚拟孔径,从而使得其性能优于相控阵雷达。
正交频分复用(Orthogonal frequency division multiplexing,OFDM)信号作为一种宽带低截获雷达波形受到越来越多的关注。OFDM雷达利用多个正交的子载波并行进行探测,从而能够有效对抗多径传播引起的频率选择性衰落,提高***的抗干扰特性。将OFDM与MIMO技术结合起来,可以充分发挥MIMO和OFDM的优势,从而能够显著提高对目标的检测性能。
空时自适应处理(STAP)是从上个世纪九十年代初发展起来的,用于对机载雷达(airborne radar)数据进行处理的技术。STAP技术在军事和民用中都有着广泛的应用,比如,地质监测,预警,地面动目标检测(GMTI),动目标检测(MTI),区域侦查等。对于传统的相控阵雷达,STAP基础理论研究已相当成熟。许多用于改善STAP复杂性以及收敛性的算法业已被提出。这些算法稍微经过修改就可以应用于MIMO雷达。
MIMO雷达发射波形设计通常基于目标和环境的先验知识进行,而此先验知识通过估计得到,因而不可避免的存在估计误差。此时,波形设计难以实现最优匹配,造成***检测性能下降。对角加载采样矩阵求逆算法(loaded sample matrix inversion,LSMI)是常见的自适应稳健方法之一,该方法通过对采样协方差矩阵进行对角加载,可提高自适应算法的收敛速度和稳健性,但其加载量由经验参数给出,没有解析解。J.Li等基于CRB准则研究了改善参数估计性能的波形设计问题。但是该波形设计问题的求解需要某些参数的确切值,比如目标位置,反射系数等。因此,优化波形的确定将依赖于这些值。工程应用中,这些参数值由于通过估计得到,因而存在不确定性。由于参数估计的最终精度对这些不确定性比较敏感,所以基于某个参数估计值得到的优化波形可能导致较差的参数估计精度。
发明内容
针对复杂环境下地面慢速目标的空时联合处理问题,本发明将参数的不确定性融入优化模型,研究了改善MIMO-OFDM雷达参数估计性能的稳健波形优化问题。本发明在恒模约束下,构建稳健波形优化模型,基于对角加载(DL)技术,将此非线性优化问题转化为可以获得高效求解的半定规划问题,用比较成熟的优化工具进行求解,仿真验证了能够有效改善输出SINR,进而最大化***检测性能。实现本发明的基本思路是,首先建立MIMO-OFDM-STAP***模型,对目标函数进行推导,然后建立稳健波形优化模型,最后完成稳健波形优化问题的求解。
本发明的技术方案是:提高基于MIMO-OFDM雷达STAP最差检测性能的稳健波形设计的方法,包括如下步骤:
一、建立MIMO-OFDM-STAP***模型
(1)MIMO-OFDM-STAP接收信号描述
考虑如说明书附图之图1所示MIMO-OFDM-STAP场景。此场景中,第l个脉冲重复间隔PRI内的接收数据可表示为:
式中,为第m个发射阵元在每个PRI内发射的复基带信号的离散形式,K为波形采样数,且am为相应的信号幅度,fm=f0+mΔf,f0为信号载频,Δf为频率间隔,满足TΔf=1;ρt和ρ(θ)分别为所考虑的距离环内目标的复幅度以及位于θ的杂波反射系数;v、vt分别代表雷达平台和目标的移动速度,λ为波形中心波长;此外,表示第n个接收阵元在第l个PRI内接收的干扰以及噪声。
如果把目标距离单元中的杂波回波建模为若干独立杂波块的叠加,第l个PRI内的接收数据在接收端进行下变频处理,第l个脉冲重复间隔PRI内的接收数据可改写为:
其中,和分别表示目标及位于θi杂波的接收导向矢量,NC(NC>>NML)为杂波环采样数目,S=[s1,s2,…,sM]T表示每个PRI中的信号矩阵。 以及 和分别为目标及位于θi杂波的发射导向矢量。假设Zl的列是独立同分布的圆对称复高斯随机向量,其均值为0,协方差矩阵为未知矩阵
(2)感兴趣距离环内空时快拍表述
利用SH(SSH)-1/2作为匹配滤波器,且则相应的矢量化匹配滤波输出可表示为:
其中,IN是N×N的单位矩阵,Φ=SSH(SSH)-1/2=diag{|a1| |a2| … |aM|},diag{·}表示对角矩阵。
由上式我们可得所感兴趣距离环内总的空时快拍为:
其中,和分别表示目标及位于θi杂波的多普勒导向矢量。
二、目标函数推导
(1)最优MIMO-OFDM-STAP处理器条件下输出SINR表述
基于最小方差无畸变准则(MVDR),可得最优输出SINR可表示为:
式中,
(2)杂波高斯分布,且与干扰不相关条件下输出SINR表述简化
假设杂波独立同分布,且服从均值为0,方差为的高斯分布,则在杂波与干扰加噪声项不相关的假设下,输出SINR可简化为如下表达式:
其中,
三、稳健波形优化模型
高斯噪声环境下,可以证明最大化检测概率等价于最大化输出信干噪比。由此,基于以上分析可得,在恒模和发射总功率约束下,通过构造一个凸集来优化波形协方差矩阵(WCM)来最大化检测概率的波形优化问题可表述为
|am|=Cm
||am||2≥0
式中,P代表发射总功率。
四、稳健波形优化求解
(1)基于DL方法的正定化Φ
此优化问题包含恒模约束,显然是一个非线性优化(NP)问题,求解全局最优解时容易陷入局部最优解。同时,由于我们无法确定的性质,因此,我们不能够利用凸优化方法来解。针对此问题,我们采用对角加载方法对Φ进行对角加载,使得
式中,ρ<<λmax(Φ)即所谓的加载因子(loading factor),λmax(·)表示矩阵的最大特征值。
(2)基于正定化Φ简化输出SINR
将代入输出SINR表达式,替换为并利用矩阵求逆定理,目标函数可重新表示为:
从上式可以知道,当真实空时导向矢量位于对应于矩阵(IMNL+RTSRC)-1RTS最小特征值的特征矢量所指方向时,就会出现MIMO-OFDM-STAP检测性能最差的情况。因此,上式可重新写为
式中,λmin(·)表示矩阵的最小特征值。
由于并利用矩阵特征值性质,上式目标函数可以重新转化为
(3)基于凸优化求解稳健波形优化问题
基于上述讨论,则稳健波形优化问题可以转化为如下的SDP问题:
||am||2≥0
本发明与现有技术相比具有以下优点:
第一,针对输出SINR对参数估计误差敏感的问题,本发明考虑通过显式地将参数不确定凸集包含进波形优化模型以缓解输出SINR对参数估计误差的敏感性,从而提高MIMO-OFDM-STAP***的检测性能。
第二,提出一种新的对角加载技术求解复杂的非线性波形优化问题,将非线性波形优化问题转化为半正定规划问题,从而可以利用比较成熟的优化工具箱获得高效求解。
附图说明
图1MIMO-OFDM-STAP模型;
图2为本发明实现的流程图;
图3.在ASNR=30dB,CNR=30dB条件下得到的最优稳健发射方向图;
图4.在ASNR=30dB条件下得到的SINR随着CNR的变化曲线;
本发明的效果可通过以下仿真进一步说明:
具体实施方式
下面结合附图2对本发明的实现步骤做进一步详细描述:
一、建立MIMO-OFDM-STAP***模型
(1)MIMO-OFDM-STAP接收信号描述
考虑如说明书附图之图1所示MIMO-OFDM-STAP场景。此场景中,第l个脉冲重复间隔PRI内的接收数据可表示为:
式中,为第m个发射阵元在每个PRI内发射的复基带信号的离散形式,K为波形采样数,且am为相应的信号幅度,fm=f0+mΔf,f0为信号载频,Δf为频率间隔,满足TΔf=1;ρt和ρ(θ)分别为所考虑的距离环内目标的复幅度以及位于θ的杂波反射系数;v、vt分别代表雷达平台和目标的移动速度,λ为波形中心波长;此外,表示第n个接收阵元在第l个PRI内接收的干扰以及噪声。
如果把目标距离单元中的杂波回波建模为若干独立杂波块的叠加,第l个PRI内的接收数据在接收端进行下变频处理,第l个脉冲重复间隔PRI内的接收数据可改写为:
其中,和分别表示目标及位于θi杂波的接收导向矢量,NC(NC>>NML)为杂波环采样数目,S=[s1,s2,…,sM]T表示每个PRI中的信号矩阵。 以及 和分别为目标及位于θi杂波的发射导向矢量。假设Zl的列是独立同分布的圆对称复高斯随机向量,其均值为0,协方差矩阵为未知矩阵
(2)感兴趣距离环内空时快拍表述
利用SH(SSH)-1/2作为匹配滤波器,且则相应的矢量化匹配滤波输出可表示为:
其中,IN是N×N的单位矩阵,Φ=SSH(SSH)-12=diag{|a1| |a2| … |aM|},diag{·}表示对角矩阵。
由上式我们可得所感兴趣距离环内总的空时快拍为:
其中,和分别表示目标及位于θi杂波的多普勒导向矢量。
二、目标函数推导
(1)最优MIMO-STAP处理器条件下输出SINR表述
基于最小方差无畸变准则(MVDR),可得最优输出SINR可表示为:
式中,
(2)杂波高斯分布,且与干扰不相关条件下输出SINR表述简化
假设杂波独立同分布,且服从均值为0,方差为的高斯分布,则在杂波与干扰加噪声项不相关的假设下,输出SINR可简化为如下表达式:
其中,
三、稳健波形优化模型
高斯噪声环境下,可以证明最大化检测概率等价于最大化输出信干噪比。由此,基于以上分析可得,在恒模和发射总功率约束下,通过构造一个凸集来优化波形协方差矩阵(WCM)来最大化检测概率的波形优化问题可表述为
|am|=Cm
||am||2≥0
式中,P代表发射总功率。
四、稳健波形优化求解
(1)基于DL方法的正定化Φ
此优化问题包含恒模约束,显然是一个非线性优化(NP)问题,求解全局最优解时容易陷入局部最优解。同时,由于我们无法确定的性质,因此,我们不能够利用凸优化方法来解。针对此问题,我们采用对角加载方法对Φ进行对角加载,使得
式中,ρ<<λmax(Φ)即所谓的加载因子(loading factor),λmax(·)表示矩阵的最大特征值。
(2)基于正定化RC简化输出SINR
将代入输出SINR表达式,替换为并利用矩阵求逆定理,目标函数可重新表示为:
从上式可以知道,当真实空时导向矢量位于对应于矩阵(IMNL+RTSRC)-1RTS最小特征值的特征矢量所指方向时,就会出现MIMO-STAP检测性能最差的情况。因此,上式可重新写为
式中,λmin(·)表示矩阵的最小特征值。
由于并利用矩阵特征值性质,上式目标函数可以重新转化为
(3)基于凸优化求解稳健波形优化问题
基于上述讨论,则稳健波形优化问题可以转化为如下的SDP问题:
||am||2≥0
本发明的效果可通过以下仿真进一步说明:
仿真条件:MIMO雷达是4发2收,接收阵元间距为半波长,发射阵元间距为2倍波长,脉冲数为3,采用两个MIMO雷达对目标进行检测,分别为A(0.5,0.5),B(1.5,0.5),阵列信噪比的定义为其中,P指总发射功率,指所附加的白色热噪声的方差,阵列的信噪比从10到50分贝变化,杂波块数为10000,杂噪比从10到50分贝变化,干扰噪声比值为60分贝,采样点数为256。本文假设在4°的方向有目标,杂波的建模使用离散点,其RCS建模为独立同分布的高斯随机变向量,均值为零,方差为并假设固定在相干处理间隔。仿真中,把本文提出的算法和不相关波形进行对比,可以看到信噪比的改善情况。
仿真内容:
仿真1:用本发明方法勾画出最优稳健发射方向图,如图3。从图3中看出,本发明提出的方法在目标附近产生一个尖峰。换言之,MIMO-OFDM-STAP在参数误差凸集上的最差检测性能可以通过所提方法得到改善。此外,还可以看到在MIMO雷达(1.5,0.5)中出现了栅瓣,这是由于此雷达稀疏的发射阵元布置。
仿真2:用本发明方法勾画出SINR随着CNR的变化曲线,如图4。从图4可以看出,与不相关波形相比,本发明所提方法能得到更大的最坏情况下的输出SINR,即所提方法可以显著提高不相关波形条件下的最坏情况下的MIMO-OFDM-STAP检测性能。
综上,本发明提出了一种稳健波形设计方法,将参数不确定凸集显式地包含进波形优化模型以最大化最坏情况下的输出SINR。为求解复杂的非线性优化问题,本发明基于对角加载(DL)方法,将此非线性优化问题转化为可以获得高效求解的半定规划问题,从而最大化输出SINR,进而最大化***检测性能。仿真表明,与非相关发射波形相比,本发明所提方法得到的发射波形可显著改善***检测性能。基于以上讨论可知,本发明所提方法可为工程应用中通过设计发射波形提高雷达检测***的稳健性能提供坚实的理论与实现依据。
Claims (1)
1.提高基于MIMO-OFDM雷达STAP最差检测性能的稳健波形设计方法,其特征在于,包括如下步骤:
一、建立MIMO-OFDM-STAP***模型
(1)MIMO-OFDM-STAP接收信号描述
在MIMO-OFDM-STAP场景中,第l个脉冲重复间隔PRI内的接收数据可表示为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>t</mi>
</msub>
<msub>
<mi>s</mi>
<mi>m</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mo>/</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>sin&theta;</mi>
<mi>t</mi>
</msub>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>r</mi>
</msub>
<mi>n</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>t</mi>
</msub>
<mi>m</mi>
<mo>+</mo>
<mn>2</mn>
<mi>v</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>m</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</msup>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msubsup>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</msub>
<msub>
<mi>s</mi>
<mi>m</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mo>/</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>r</mi>
</msub>
<mi>n</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>t</mi>
</msub>
<mi>m</mi>
<mo>+</mo>
<mn>2</mn>
<mi>v</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>m</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</msup>
<mo>+</mo>
<mi>d</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,n表示第n个接收阵元,M为发射阵元个数,θt为目标所处位置,dt,dr分别表示发射接收阵列中相邻阵元间距,t为阵列接收时间,为第m个发射阵元在每个PRI内发射的复基带信号的离散形式,K为波形采样数,且am为相应的信号幅度,fm=f0+mΔf,f0为信号载频,Δf为频率间隔,满足TΔf=1,T为信号重复周期;ρt和ρ(θ)分别为所考虑的距离环内目标的复幅度以及位于θ的杂波反射系数;v、vt分别代表雷达平台和目标的移动速度,λ为波形中心波长;此外,表示第n个接收阵元在第l个PRI内接收的干扰以及噪声;
如果把目标距离单元中的杂波回波建模为若干独立杂波块的叠加,第l个PRI内的接收数据在接收端进行下变频处理,第l个脉冲重复间隔PRI内的接收数据可改写为:
<mrow>
<msub>
<mi>X</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>t</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>D</mi>
</msub>
<mi>l</mi>
</mrow>
</msup>
<msup>
<mi>ab</mi>
<mi>T</mi>
</msup>
<mi>S</mi>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;&beta;f</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mi>l</mi>
</mrow>
</msup>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<msup>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mi>T</mi>
</msup>
<mi>S</mi>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mi>l</mi>
</msub>
</mrow>
其中,NC为杂波环采样数目,ρi为位于θi杂波的复幅度, 和分别表示目标及位于θi杂波的接收导向矢量,和分别为目标及位于θi杂波的发射导向矢量, 以及S=[s1,s2,…,sM]T表示每个PRI中的信号矩阵,假设Zl的列是独立同分布的圆对称复高斯随机向量,其均值为0,协方差矩阵为未知矩阵
(2)感兴趣距离环内空时快拍表述
利用SH(SSH)-1/2作为匹配滤波器,且δmq为狄拉克算子,则相应的矢量化匹配滤波输出可表示为:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mi>l</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>t</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>D</mi>
</msub>
<mi>l</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>&CircleTimes;</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;&beta;f</mi>
<mrow>
<mi>s</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mi>l</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>v</mi>
<mi>e</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>Z</mi>
<mo>~</mo>
</mover>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中,Φ=SSH(SSH)-1/2=diag{|a1| |a2|…|aM|},diag{·}表示对角矩阵,IN表示N×N的单位矩阵,N为接收阵元个数;
由上式可得所感兴趣距离环内总的空时快拍为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mi>C</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>t</mi>
</msub>
<msub>
<mi>U</mi>
<mi>D</mi>
</msub>
<mo>&CircleTimes;</mo>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>(</mo>
<mrow>
<mi>b</mi>
<mo>&CircleTimes;</mo>
<mi>a</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&CircleTimes;</mo>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>v</mi>
<mi>e</mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>Z</mi>
<mo>~</mo>
</mover>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mi>D</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>b</mi>
<mo>&CircleTimes;</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>v</mi>
<mi>e</mi>
<mi>c</mi>
<mo>(</mo>
<msub>
<mover>
<mi>Z</mi>
<mo>~</mo>
</mover>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,和分别表示目标及位于θi杂波的多普勒导向矢量,L为发射脉冲个数,fD,i为位于θi杂波的多普勒频率;
二、目标函数推导
(1)最优输出SINR表述
基于最小方差无畸变准则,可得最优输出SINR可表示为:
<mrow>
<mi>S</mi>
<mi>I</mi>
<mi>N</mi>
<mi>R</mi>
<mo>=</mo>
<msup>
<msub>
<mi>&rho;</mi>
<mi>t</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mi>D</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>b</mi>
<mo>&CircleTimes;</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>H</mi>
</msup>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<msup>
<mi>&alpha;</mi>
<mi>T</mi>
</msup>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mi>D</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>b</mi>
<mo>&CircleTimes;</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
式中,
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>v</mi>
<mi>e</mi>
<mi>c</mi>
<mo>(</mo>
<msub>
<mover>
<mi>Z</mi>
<mo>~</mo>
</mover>
<mi>l</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
</mrow>
<mo>)</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>v</mi>
<mi>e</mi>
<mi>c</mi>
<mo>(</mo>
<msub>
<mover>
<mi>Z</mi>
<mo>~</mo>
</mover>
<mi>l</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,E(·)为期望算子;
(2)杂波高斯分布,且与干扰不相关条件下输出SINR表述简化
假设杂波独立同分布,且服从均值为0,方差为的高斯分布,则在杂波与干扰加噪声项不相关的假设下,输出SINR可简化为如下表达式:
<mrow>
<mi>S</mi>
<mi>I</mi>
<mi>N</mi>
<mi>R</mi>
<mo>=</mo>
<mo>|</mo>
<msub>
<mi>&rho;</mi>
<mi>t</mi>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mi>H</mi>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>+</mo>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<msub>
<mi>R</mi>
<mi>C</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
</mrow>
其中,
<mrow>
<msub>
<mi>R</mi>
<mi>C</mi>
</msub>
<mo>=</mo>
<msup>
<mi>V&Xi;V</mi>
<mi>H</mi>
</msup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
<mo>,</mo>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mo>=</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>M</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>Q</mi>
<mo>></mo>
<mn>0</mn>
<mo>,</mo>
<mi>A</mi>
<mo>=</mo>
<msub>
<mi>I</mi>
<mi>L</mi>
</msub>
<mo>&CircleTimes;</mo>
<mi>&Phi;</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>I</mi>
<mi>N</mi>
</msub>
<mo>;</mo>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<msup>
<mi>A</mi>
<mi>H</mi>
</msup>
<msup>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>A</mi>
<mo>;</mo>
</mrow>
三、稳健波形优化模型
高斯噪声环境下,可以证明最大化检测概率等价于最大化输出信干噪比,由此,基于以上分析可得,在恒模和发射总功率约束下,通过构造一个凸集来优化波形协方差矩阵来最大化检测概率的波形优化问题可表述为
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
</munder>
<munder>
<mi>min</mi>
<msub>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>t</mi>
</msub>
</munder>
<msubsup>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>t</mi>
<mi>H</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>M</mi>
<mi>N</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>+</mo>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<msub>
<mi>R</mi>
<mi>C</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<msub>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>t</mi>
</msub>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>t</mi>
</msub>
<mo>&Element;</mo>
<mi>V</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
|am|=Cm
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msubsup>
<mi>a</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mi>P</mi>
</mrow>
||am||2≥0
式中,为真实目标空时导向矢量;Cm为第m个发射信号的幅度常数;IMNL为维数为MNL的单位矩阵;P代表发射总功率;
四、稳健波形优化求解
(1)基于DL方法的正定化Φ
此优化问题包含恒模约束,显然是一个非线性优化问题,求解全局最优解时容易陷入局部最优解,同时,由于无法确定的性质,因此,不能够利用凸优化方法来解,针对此问题,采用对角加载方法对Φ进行对角加载,使得
<mrow>
<mover>
<mi>&Phi;</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mi>&Phi;</mi>
<mo>+</mo>
<mi>&rho;</mi>
<mi>I</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
式中,ρ<<λmax(Φ)即所谓的加载因子—loading factor,λmax(·)表示矩阵的最大特征值;
(2)基于正定化RC简化输出SINR
将代入输出SINR表达式,替换为并利用矩阵求逆定理,目标函数可重新表示为:
<mrow>
<msubsup>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>t</mi>
<mi>H</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>M</mi>
<mi>N</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>+</mo>
<mover>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mo>~</mo>
</mover>
<msub>
<mi>R</mi>
<mi>C</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mover>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mo>~</mo>
</mover>
<msub>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mi>t</mi>
</msub>
</mrow>
从上式可以知道,当真实空时导向矢量位于对应于矩阵最小特征值的特征矢量所指方向时,就会出现MIMO-STAP检测性能最差的情况,因此,上式可重新写为
<mrow>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>M</mi>
<mi>N</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>+</mo>
<mover>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mo>~</mo>
</mover>
<msub>
<mi>R</mi>
<mi>C</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mover>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
式中,λmin(·)表示矩阵的最小特征值,η为真实目标空时导向矢量,的模的上界;
由于并利用矩阵特征值性质,上式目标函数可以重新转化为
<mrow>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mo>~</mo>
</mover>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>C</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
(3)基于凸优化求解稳健波形优化问题
基于上述讨论,则稳健波形优化问题可以转化为如下的SDP问题:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>&Psi;</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</munder>
</mtd>
<mtd>
<mrow>
<msup>
<mi>&eta;</mi>
<mn>2</mn>
</msup>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>C</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>I</mi>
</mtd>
<mtd>
<mover>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mo>~</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msubsup>
<mi>a</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<msub>
<mi>D</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msubsup>
<mi>a</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mi>P</mi>
</mrow>
||am||2≥0
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510756071.9A CN105487054B (zh) | 2015-11-09 | 2015-11-09 | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510756071.9A CN105487054B (zh) | 2015-11-09 | 2015-11-09 | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105487054A CN105487054A (zh) | 2016-04-13 |
CN105487054B true CN105487054B (zh) | 2018-05-04 |
Family
ID=55674163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510756071.9A Active CN105487054B (zh) | 2015-11-09 | 2015-11-09 | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105487054B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105807275A (zh) * | 2016-04-28 | 2016-07-27 | 大连大学 | 基于部分杂波先验知识的mimo-ofdm-stap稳健波形设计 |
CN107167782B (zh) * | 2017-06-27 | 2020-04-10 | 西安电子科技大学 | 基于信杂噪比最大的雷达三维异构阵稀疏重构方法 |
CN108008364B (zh) * | 2017-11-22 | 2021-03-19 | 大连大学 | 改善mimo-stap检测性能的发射波形及接收权联合优化方法 |
CN108469601B (zh) * | 2018-01-17 | 2021-11-30 | 大连大学 | 改善mimo-stap最差检测性能的收发联合稳健优化方法 |
WO2019159112A1 (en) * | 2018-02-14 | 2019-08-22 | Tiejun Shan | Method for location approximation |
CN110031821B (zh) * | 2019-03-25 | 2020-11-17 | 白犀牛智达(北京)科技有限公司 | 一种车载避障激光雷达波形提取方法、激光雷达及介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8488724B2 (en) * | 2009-05-14 | 2013-07-16 | Silvus Technologies, Inc. | Wideband interference mitigation for devices with multiple receivers |
CN103852749A (zh) * | 2014-01-28 | 2014-06-11 | 大连大学 | 提高mimo-stap检测性能的稳健波形优化方法 |
-
2015
- 2015-11-09 CN CN201510756071.9A patent/CN105487054B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105487054A (zh) | 2016-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105487054B (zh) | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 | |
CN105319545B (zh) | 提高stap检测性能的mimo-ofdm雷达波形设计方法 | |
Xu et al. | Joint range and angle estimation using MIMO radar with frequency diverse array | |
CN105699945B (zh) | 频控阵mimo雷达***中的波形优化设计方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN107656257B (zh) | 一种弹载mimo雷达波形协方差矩阵的优化设计方法 | |
Tohidi et al. | Sparse antenna and pulse placement for colocated MIMO radar | |
CN109782243B (zh) | 基于块Hankel矩阵填充的阵元故障MIMO雷达角度估计方法 | |
CN105223560B (zh) | 基于杂波俯仰方位谱稀疏恢复的机载雷达目标检测方法 | |
CN107167783B (zh) | 一种共形阵列杂波协方差矩阵的稀疏重构方法 | |
CN105137410B (zh) | 基于ofdm的高分辨率雷达通信一体化的波形优化方法 | |
CN107015213B (zh) | 基于music算法的双基地mimo雷达角度估算方法 | |
CN107390197B (zh) | 基于特征空间的雷达自适应和差波束测角方法 | |
CN105807275A (zh) | 基于部分杂波先验知识的mimo-ofdm-stap稳健波形设计 | |
CN102520400B (zh) | 非高斯杂波环境下的mimo雷达目标检测的仿真方法 | |
CN107703489B (zh) | 一种mimo雷达恒模波形和接收机的联合设计方法 | |
CN104833959A (zh) | 基于目标先验信息的mimo雷达波形优化方法 | |
CN105301580A (zh) | 一种基于***阵互谱相位差方差加权的被动探测方法 | |
CN112612006B (zh) | 基于深度学习的机载雷达非均匀杂波抑制方法 | |
CN105182313A (zh) | 一种基于不完全杂波先验知识的mimo-stap稳健波形设计方法 | |
Chen et al. | Sparse reconstruction based target localization with frequency diverse array MIMO radar | |
CN108828504B (zh) | 基于部分相关波形的mimo雷达目标方向快速估计方法 | |
CN110082744A (zh) | 多普勒分级复用的mimo机载双基雷达杂波抑制方法 | |
CN108896974B (zh) | 一种改进的mimo阵列高分辨空间谱估计方法 | |
CN111856410B (zh) | 恒模约束下机载mimo雷达的稳健发射接收联合优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |