CN107167782A - 基于信杂噪比最大的雷达三维异构阵稀疏重构方法 - Google Patents
基于信杂噪比最大的雷达三维异构阵稀疏重构方法 Download PDFInfo
- Publication number
- CN107167782A CN107167782A CN201710499046.6A CN201710499046A CN107167782A CN 107167782 A CN107167782 A CN 107167782A CN 201710499046 A CN201710499046 A CN 201710499046A CN 107167782 A CN107167782 A CN 107167782A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- battle array
- msubsup
- dimensional
- 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.)
- Granted
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/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
-
- 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/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
Abstract
本发明公开了一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，思路为：确定三维圆柱阵雷达，该三维圆柱阵雷达的检测范围内存在目标，获取三维异构阵机载雷达杂波空时导向矢量矩阵；将三维异构阵机载雷达杂波空时导向矢量矩阵转换为二维导向矢量矩阵，并计算目标加杂波的协方差矩阵的行列式；计算三维异构阵机载雷达的输出信杂噪比和三维异构阵机载雷达的输出信杂噪比的行列式表示形式；定义二进制挑选向量，计算三维异构阵机载雷达的输出信杂噪比的行列式表示形式的重写式，并计算三维异构阵机载雷达的稀疏降维矩阵以及三维异构阵机载雷达的滤波结果，所述三维异构阵机载雷达的滤波结果为基于信杂噪比最大的雷达三维异构阵稀疏重构结果。
Description
技术领域
本发明属于雷达技术领域，涉及一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，适用于机载三维异构阵雷达的杂波抑制。
背景技术
与机身共形的天线具有诸多优点，比如更好的气动特性、更好的隐身性能、更大的天线面积等，因此其应用前景广阔；但机身的多曲面特性使得天线也具有多曲面结构，进而产生三维异构阵问题；相对于二维平面阵来说，三维异构阵的杂波特性将会变得十分复杂，这种复杂性体现在单元间的互耦特性更加复杂，单元的极化改变更加复杂；三维异构阵给杂波带来的复杂性也给杂波抑制带来了很大的困难包括阵列流型更加复杂，目标和杂波的导向矢量计算困难，导向矢量间的相关性更难计算，杂波呈现距离相关性非均匀。
在杂波协方差矩阵和目标信号均确知的条件下，Brennan等人在1973年，根据线性约束最小方差准则，推导出输出信杂噪比最大的全空时二维自适应处理(Space TimeAdaptive Processing，STAP)的概念和理论，全STAP能够取得比较理想的杂波抑制效果，但是最优的STAP处理器需要用大量的训练样本来估计杂波加噪声协方差矩阵。
在三维异构阵的情况下，由于雷达接收的空时数据维数往往很大，全维STAP处理器在实际应用中面临两个问题：1)三维阵的阵面变大，阵元数变多使空时两维自由度庞大导致杂波协方差矩阵求逆的运算量十分巨大，无法满足实时处理的要求；2)缺乏用于估计杂波协方差矩阵的独立同分布(Independent Identically Distributed,IID)样本单元。根据Reed,Mallett,and Brennan三个人提出的准则，即RMB准则可知，为使因估计不准确而带来的输出信杂噪比损失在3dB以内，要求样本数目至少为协方差矩阵阶数的两倍，但实际工程中会面临强烈的非均匀杂波环境，不同距离的杂波样本分布特性各异，这将造成独立同分布IID样本支撑困难。因此，这两个问题使得在三维阵列的情况下使用全维STAP性能降低，杂波抑制效果变差。
为了摆脱在三维阵情况下雷达接接收数据很大时杂波协方差矩阵得估计对于样本的过度依赖，导致的杂波抑制效果降低的问题，有必要对三维阵列的空时自适应处理降维处理进行研究。
发明内容
针对上述现有技术存在的问题，本发明的目的在于提出一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，该种基于信杂噪比最大的雷达三维异构阵稀疏重构方法能够解决三维异构阵情况下全空时自适应处理的杂波协方差矩阵估计需要大量训练样本的问题，能够降低杂波协方差矩阵估计所需的训练样本数，提高杂波抑制性能。
本发明的主要思路：利用三维异构阵杂波回波数据得到杂波的协方差矩阵和杂波加目标的协方差矩阵，推导出输出信杂噪比的行列式表示形式，用使输出信杂噪比最大的阵元选择方法得到最优的阵元选择矢量进而得到降维矩阵，得到降维后的改善因子和杂波抑制结果。
为达到上述技术目的，本发明采用如下技术方案予以实现。
一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，包括以下步骤：
步骤1，确定三维圆柱阵雷达，所述三维圆柱阵雷达的检测范围内存在目标，并获取三维异构阵机载雷达杂波空时导向矢量矩阵AN×K×L；
其中，N为三维圆柱阵雷达包含的阵元个数，K为每个相干处理间隔内三维圆柱阵雷达发射的脉冲个数，L表示表示三维圆柱阵雷达接收到的杂波包含的距离门总个数，N、K、L分别为大于0的正整数；
步骤2，将三维异构阵机载雷达杂波空时导向矢量矩阵AN×K×L转换为NK×L的二维导向矢量矩阵BNK×L，并计算得到目标加杂波的协方差矩阵Rs的行列式；
步骤3，计算三维异构阵机载雷达的输出信杂噪比SCNRout；
步骤4，计算得到三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式；
步骤5，定义二进制挑选向量y，所述二进制挑选向量y为NK×1维列矢量，并根据三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式，计算得到三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式；
步骤6，根据三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式,得到三维异构阵机载雷达的稀疏降维矩阵；
步骤7，根据三维异构阵机载雷达的稀疏降维矩阵,计算得到三维异构阵机载雷达的滤波结果；所述三维异构阵机载雷达的滤波结果为基于信杂噪比最大的雷达三维异构阵稀疏重构结果。
本发明的有益效果：在三维异构阵情况下，由于阵元数增多使得全空时自适应处理的自由度明显增大，从而需要更多的训练样本来估计杂波协方差矩阵，本发明方法通过在信杂噪比最大的情况下进行阵元挑选，挑选后自由度降低为全空时的一半，能够有效抑制杂波，并能够很好解决三维异构阵情况下杂波协方差矩阵估计的训练样本数有限的问题。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1是本发明的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法流程图；
图2(a)为在少样本数情况下使用本发明方法与全空时处理方法分别得到的改善因子曲线对比图；
图2(b)为在多样本数情况下使用本发明方法与全空时处理方法分别得到的改善因子曲线对比图；；
图3(a)为少样本数情况下使用全空时处理方法进行杂波抑制得到的结果示意图；
图3(b)为少样本数情况下使用本发明方法进行杂波抑制得到的结果示意图。
具体实施方式
参照图1，为本发明的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法流程图；其中所述基于信杂噪比最大的雷达三维异构阵稀疏重构方法，包括以下步骤：
步骤1，确定三维圆柱阵雷达，所述三维圆柱阵雷达的检测范围内存在目标，并获取三维异构阵机载雷达杂波空时导向矢量矩阵AN×K×L；本实施例是对三维圆柱阵雷达进行仿真得到三维异构阵机载雷达杂波空时导向矢量矩阵。
具体地：AN×K×L为N×K×L维三维数据矩阵，N为三维圆柱阵雷达包含的阵元个数，三维圆柱阵雷达共五层，每层六个阵元，K为每个相干处理间隔CPI内三维圆柱阵雷达发射的脉冲个数，CPI表示相干处理间隔，L表示表示三维圆柱阵雷达接收到的杂波包含的距离门总个数，N、K、L分别为大于0的正整数。
步骤2，将三维异构阵机载雷达杂波空时导向矢量矩阵AN×K×L转换为NK×L的二维导向矢量矩阵BNK×L，并计算得到目标加杂波的协方差矩阵Rs的行列式。
具体地，将三维异构阵机载雷达杂波空时导向矢量矩阵AN×K×L转换为NK×L的二维导向矢量矩阵BNK×L，计算NK×L的二维导向矢量矩阵BNK×L的秩Na，Na＝rank(BNK×L)，0≤Na≤min{NK,L}，rank表示求秩运算，min表示求最小值操作。
然后对NK×L的二维导向矢量矩阵BNK×L进行奇异值分解，得到NK×NK维酉矩阵UNK×NK，取NK×NK维酉矩阵UNK×NK的前Na列，记为NK×Na维中间矩阵VC，
vi表示NK×NK维酉矩阵UNK×NK中的第i列矢量，且vi为NK×1维列矢量。
根据NK×Na维中间矩阵VC，计算得到三维异构阵机载雷达的杂波协方差矩阵Rc，其定义式为：
分别设定目标的方位角为φ，目标的俯仰角为θ，目标的多普勒频率为fd，则将目标的空间频率记为fs，将目标的归一化多普勒频率记为 N为三维圆柱阵雷达包含的阵元个数,λ表示三维圆柱阵雷达发射信号波长，Fr为脉冲重复频率，cos表示求余弦操作，sin表示求正弦操作，那么目标的空域导向矢量为a，a为N×1维列矢量，n∈{0,1,…,N-1}，an表示目标的空域导向矢量a中第n个元素，N为三维圆柱阵雷达包含的阵元个数，与目标的空域导向矢量a中包含的元素个数取值相等；上标T表示转置操作，目标的时域导向矢量为b，b为K×1维列矢量，bm表示目标的时域导向矢量b中第m个元素，K为每个相干处理间隔CPI内三维圆柱阵雷达发射的脉冲个数，与目标的时域导向矢量b中包含的元素个数取值相等；CPI表示相干处理间隔；进而计算得到目标的空时导向矢量s，其计算表达式为：
其中，表示kronecker积，s为NK×1维列矢量，目标的空时导向矢量s的定义表达式为：s＝[s1,s2,...,si',...,sNK]T，上标T为转置操作，i'∈{1,2,…,NK}，si'表示目标的空时导向矢量s中第i'个元素，si'＝anbm。
定义目标加杂波的信号矩阵为Vs，Vs＝[s,Vc]＝[s,v1,v2,...,vi,...,vNa]，进而计算得到目标加杂波的协方差矩阵为Rs，sHs＝NK，上标H表示共轭转置操作；根据块矩阵求行列式的公式，得到目标加杂波的协方差矩阵Rs的行列式|Rs|，
步骤3，计算三维异构阵机载雷达的输出信杂噪比SCNRout。
具体地：(1)将三维异构阵机载雷达的杂波加噪声协方差矩阵表示为R，且R为NK×NK的二维矩阵，表示高斯白噪声的功率，Rc为三维异构阵机载雷达的杂波协方差矩阵，INK表示N×K维全1矩阵。
由步骤1和步骤2可知三维异构阵机载雷达的杂波协方差矩阵Rc的计算表达式为：vi表示NK×NK维酉矩阵UNK×NK中的第i列矢量，且vi为NK×1维列矢量；qi表示NK×NK维酉矩阵UNK×NK中第i列矢量的杂波功率，上标H表示共轭转置操作。
定义Na×1维矢量并对Na×1维矢量进行对角化，即是将Na×1维矢量对角化，进而得到Na×Na维二维矩阵Q，diag表示对角化操作，得到三维异构阵机载雷达的杂波加噪声协方差矩阵的简化形式
根据矩阵求逆公式计算得到三维异构阵机载雷达的杂波加噪声协方差矩阵的简化形式的逆
假定NK×NK维酉矩阵UNK×NK中每一列矢量的杂波功率都远大于噪声功率，即qi>>则得到三维异构阵机载雷达的杂波加噪声协方差矩阵的简化形式的逆的简化式 进而计算得到三维异构阵机载雷达的空时自适应处理的最优权Wopt，
(2)将目标的空时导向矢量s分解为sc和s⊥两个正交的子空间，s＝sc+s⊥，sc和s⊥正交，并且互为补空间，即sc∪s⊥＝s， 表示空集，sc表示杂波子空间，s⊥表示杂波子空间sc的正交补空间，VC表示NK×Na维中间矩阵，INK表示N×K维全1矩阵，上标H表示共轭转置，K为每个相干处理间隔CPI内三维圆柱阵雷达发射的脉冲个数，N为三维圆柱阵雷达包含的阵元个数，Na表示NK×L的二维导向矢量矩阵BNK×L的秩，L表示三维圆柱阵雷达接收到的杂波包含的距离门总个数。
由步骤1可知，NK×Na维中间矩阵vi表示NK×NK维酉矩阵UNK×NK中的第i列矢量，且vi为NK×1维列矢量，Na个NK×1维列矢量张成了杂波子空间sc，P表示正交投影矩阵，根据正交投影的概念，上标H表示共轭转置操作，上标-1表示求逆操作；由此计算得到目标的空时导向矢量s与杂波子空间sc之间夹角α的余弦值cos(α)，其中，|·|表示求绝对值操作，||s||2表示目标的空时导向矢量s的二范数，将sc的表达式代入|cos(α)|的表达式后得到|cos(α)|的展开式，再对|cos(α)|的展开式进行求平方操作，得到目标的空时导向矢量s与杂波子空间sc之间夹角α的余弦平方值|cos(α)|2，
(3)根据目标的空时导向矢量s与杂波子空间sc之间夹角α的余弦平方值|cos(α)|2，计算得到三维异构阵机载雷达的输出信杂噪比SCNRout，
s表示目标的空时导向矢量，表示三维异构阵机载雷达的杂波加噪声协方差矩阵的简化形式的逆的简化式，SNR表示目标的信噪比， 表示目标的信号功率，代表高斯白噪声的功率，N为三维圆柱阵雷达包含的阵元个数，三维圆柱阵雷达共五层，每层六个阵元，K表示每个相干处理间隔CPI内三维圆柱阵雷达发射的脉冲个数。
步骤4：计算得到三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式 Rc为三维异构阵机载雷达的杂波协方差矩阵，Rs表示目标加杂波的协方差矩阵，|Rs|表示目标加杂波的协方差矩阵Rs的行列式，|Rc|表示三维异构阵机载雷达的杂波协方差矩阵Rc的行列式。
具体地，对步骤2得到的目标加杂波的协方差矩阵Rs的行列式
进行变形，得到然后将代入目标的空时导向矢量s与杂波子空间sc之间夹角α的余弦平方值|cos(α)|2的表达式中，得到最后将代入到SCNRout的表达式中，可将三维异构阵机载雷达的输出信杂噪比SCNRout改写为用Rs和Rc的行列式表示的形式，进而得到三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式
步骤5：定义一个二进制挑选向量y，所述二进制挑选向量y为NK×1维列矢量，N为三维圆柱阵雷达包含的阵元个数，K为每个相干处理间隔CPI内三维圆柱阵雷达发射的脉冲个数，NK表示一共有NK个阵元脉冲对；二进制挑选向量y中包含NK个元素，每个元素的取值分别为0或1，1代表选择了对应的阵元脉冲，0代表没有选择任何阵元脉冲，那么将三维异构阵机载雷达的杂波协方差矩阵Rc和目标加杂波的协方差矩阵Rs分别表示为Rc(y)和Rs(y)，diag表示对角化操作，Rc(y)表示被挑选后的杂波协方差矩阵，Rs(y)表示被挑选后的杂波加目标协方差矩阵；进而得到三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式
因此，选择二进制挑选向量y的取值问题就转化为使三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式最大的问题，也就是让最小的问题。
步骤6，在三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式取值最大的情况下进行阵元挑选。
具体地：假设在总体的N个阵元与K个脉冲的乘积NK个阵元脉冲对中挑选M个阵元脉冲对，则二进制挑选向量y里面有M个1，NK-M个0；所以，在三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式取值最大的情况下进行阵元挑选的目的是使得二进制挑选向量y里面有M个1，NK-M个0。
6.1令二进制挑选向量为y，并令将二进制挑选向量y的初始值全部为1，记为二进制挑选向量初始值y(0)，即y(0)＝[1,1,...,1]NK，二进制挑选向量初始值y(0)的下标集为η(0)，η(0)＝[1,2,...,NK]。
6.2定义外层循环变量k，k为第k次迭代，k的初始值为1，k∈{1,2,…,NK-M}，M<NK；定义内层循环变量g，g为第g次迭代，g随着k的变化而变化,g∈{1,2,…,NK-k+1},g的初始值为1；当k取值为1时，g取值从1循环到NK；当k取值为2时，g取值从1循环到NK-1；当k取值为M时，g取值从1循环到NK-M+1；外层循环变量每循环一次，就会将第k次迭代后的二进制挑选向量y(k)中的一个1置为0，当外层循环变量循环M次时，就会将第M次迭代后的二进制挑选向量y(M)中的M个1都置为0。
6.3根据第k次迭代后的二进制挑选向量y(k),得到第k次迭代后的二进制挑选向量y(k)的下标集η(k),所述第k次迭代后的二进制挑选向量y(k)的下标集η(k)中与第k次迭代后的二进制挑选向量y(k)中元素为0的对应位置处元素为0,y(k)的下标集η(k)中其余位置处的元素按照自然数顺序从1排列至NK-k+1；所述第k次迭代后的二进制挑选向量y(k)为第k-1次迭代后的二进制挑选向量y(k-1)。
6.4令第g次迭代后下标集η(k)中与g取值相等的元素在η(k)中的位置为h，即η(k)(h)＝g，h∈{1,2,...,NK-k+1}，然后将第k次迭代后的二进制挑选向量y(k)中位置h处的元素置为0，即y(k)(h)＝0，并将第k次迭代后的二进制挑选向量y(k)中位置h处的元素置为0后得到的向量,记为第g次迭代后在位置h处置零的向量且满足即向量中位置h处的元素为0。
然后计算第g次迭代后在位置h处置零的向量的杂波协方差矩阵和第g次迭代后在位置h处置零的向量的杂波加目标协方差矩阵以及分别计算杂波协方差矩阵的行列式和杂波加目标协方差矩阵的行列式进而计算杂波协方差矩阵的行列式和杂波加目标协方差矩阵的行列式的比值μk(g)，其表达式分别为:
其中，杂波协方差矩阵的行列式和杂波加目标协方差矩阵的行列式的维数都是NK×NK,diag表示对角化操作；然后,将杂波协方差矩阵的行列式和杂波加目标协方差矩阵的行列式的比值μk(g),记为第k次迭代后经过第g次迭代得到的第g个比值。
6.5令g分别取1至NK-k+1，重复执行6.3和6.4，进而分别得到第k次迭代后经过第1次迭代得到的第1个比值至第k次迭代后经过第NK-k+1次迭代得到的第NK-k+1个比值，记为第k次迭代后得到的NK-k+1个比值,比较第k次迭代后得到的NK-k+1个比值，得到比值最小值，并得到比值最小值时g的对应取值在第k次迭代后的二进制挑选向量y(k)的下标集η(k)中的位置，记为hmin，hmin∈{1,2,…,NK-k+1}，然后将第k次迭代后的二进制挑选向量y(k)中位置hmin处的1置为0，即y(k)(hmin)＝0，进而得到第k次迭代后的二进制挑选向量y(k)，所述第k次迭代后的二进制挑选向量y(k)中有k个0，NK-k个1。
6.6令k加1，返回6.3，直到得到第NK-M次迭代后的二进制挑选向量y(NK-M)，迭代停止，此时得到的第NK-M次迭代后的二进制挑选向量y(NK-M)中NK-M个位置处的值分别已置为0，其余M个位置处的值分别为1，即迭代的最终结果使得第NK-M次迭代后的二进制挑选向量y(NK-M)中有M个1，NK-M个0。
然后计算得到对角矩阵Y，Y＝diag(y(NK-M))，diag表示对角化操作，对角矩阵Y为NK行NK列的二维矩阵，由于第NK-M次迭代后的二进制挑选向量y(NK-M)中有NK-M个元素为0，因此对角矩阵Y中有NK-M列全部为零，剔除对角矩阵Y中全部为0的列，并将对角矩阵Y中剔除全部为0的列后的矩阵，记为三维异构阵机载雷达的稀疏降维矩阵Z，Z是NK行M列的二维矩阵。
步骤7：由步骤3得到的空时自适应处理的最优权为使用三维异构阵机载雷达的稀疏降维矩阵Z对三维异构阵机载雷达的杂波加噪声协方差矩阵R和目标的空时导向矢量s分别进行降维处理，分别计算得到稀疏后的自适应权值和稀疏后的改善因子IF，进而计算得到三维异构阵机载雷达的滤波结果out，其表达式分别为：
其中，E表示三维异构阵雷达的杂波功率,上标H表示共轭转置操作，上标-1表示求逆操作；所述三维异构阵机载雷达的滤波结果out即为基于信杂噪比最大的雷达三维异构阵稀疏重构结果。
下面结合仿真实验对本发明效果作进一步验证说明：
(一)仿真参数：
三维异构阵机载雷达天线采用圆柱阵进行仿真，圆柱阵为5层，每层6个阵元，共30个阵元，波长为0.2m，阵元间距为半波长0.1m，相干积累脉冲数20，X轴向阵元间隔为0.1m，三维异构阵机载雷达的载机高度为6000m，三维异构阵机载雷达的载机速度为200m/s，载机速度方向沿X轴，三维异构阵机载雷达发射的信号带宽为2.5MHz，采样频率2.5MHz，脉冲重复频率为5000Hz，信号带宽为2.5MHz，主波束指向为，方位角90°，俯仰角0°，噪声系数3dB，杂噪比40dB。
(二)仿真结果及分析：
为了说明本发明性能的优越性，给出了传统方法以及本发明方法的处理结果对比图。
参照图2(a)，为在少样本数情况下使用本发明方法与全空时处理方法分别得到的改善因子曲线对比图；其中，纵坐标表示改善因子的大小，单位为分贝(dB)，横坐标表示多普勒通道，共20个，实线为全空时处理方法，虚线为本发明方法。
从图2(a)可以看出，由于采用了阵元挑选，选择了一半的阵元脉冲数，也就是NK/2＝300个阵元脉冲数，所以本发明方法的自由度较全空时处理的自由度NK＝600降低了一半，因此所需的训练样本数也是全空时的一半，因此在样本数为700时，对于本发明方法已经满足了两倍自由度的要求，而对于全空时处理方法并没有满足，因此在少样本情况下，本发明方法的改善因子要优于全空时处理的改善因子。
图2(b)为在多样本数情况下使用本发明方法与全空时处理方法分别得到的改善因子曲线对比图，选择的样本数为全空时情况下的五倍自由度，在多样本情况下，由于本发明方法直接去掉了一半的阵元脉冲，因此结果没有全空时的改善因子好，这也说明了，本发明方法更加适用于样本数缺乏的情况下进行杂波抑制。
图3(a)为少样本数情况下使用全空时处理方法进行杂波抑制得到的结果示意图，图3(b)为少样本数情况下使用本发明方法进行杂波抑制得到的结果示意图，图3(a)和图3(b)的横坐标均为多普勒通道数，共20个；纵坐标均为距离门数，从0到3000，共3000个距离门。
对比图3(a)和图3(b)可以看出在少样本情况下全空时处理方法的杂波剩余比较多，本发明方法在少样本情况下的滤波结果比全空时处理结果好，杂波剩余更少，因此杂波抑制性能更好。
综上所述，仿真实验验证了本发明的正确性，有效性和可靠性。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围；这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (8)
1.一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，其特征在于，包括以下步骤：
步骤1，确定三维圆柱阵雷达，所述三维圆柱阵雷达的检测范围内存在目标，并获取三维异构阵机载雷达杂波空时导向矢量矩阵AN×K×L；
其中，N为三维圆柱阵雷达包含的阵元个数，K为每个相干处理间隔内三维圆柱阵雷达发射的脉冲个数，L表示表示三维圆柱阵雷达接收到的杂波包含的距离门总个数，N、K、L分别为大于0的正整数；
步骤2，将三维异构阵机载雷达杂波空时导向矢量矩阵AN×K×L转换为NK×L的二维导向矢量矩阵BNK×L，并计算得到目标加杂波的协方差矩阵Rs的行列式；
步骤3，计算三维异构阵机载雷达的输出信杂噪比SCNRout；
步骤4，计算得到三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式；
步骤5，定义二进制挑选向量y，所述二进制挑选向量y为NK×1维列矢量，并根据三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式，计算得到三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式；
步骤6，根据三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式,得到三维异构阵机载雷达的稀疏降维矩阵；
步骤7，根据三维异构阵机载雷达的稀疏降维矩阵,计算得到三维异构阵机载雷达的滤波结果；所述三维异构阵机载雷达的滤波结果为基于信杂噪比最大的雷达三维异构阵稀疏重构结果。
2.如权利要求1所述的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，其特征在于，在步骤2中，所述目标加杂波的协方差矩阵Rs的行列式，其得到过程为：
将三维异构阵机载雷达杂波空时导向矢量矩阵AN×K×L转换为NK×L的二维导向矢量矩阵BNK×L，计算NK×L的二维导向矢量矩阵BNK×L的秩Na，Na＝rank(BNK×L)，0≤Na≤min{NK,L}，rank表示求秩运算，min表示求最小值操作；
然后对NK×L的二维导向矢量矩阵BNK×L进行奇异值分解，得到NK×NK维酉矩阵UNK×NK，取NK×NK维酉矩阵UNK×NK的前Na列，记为NK×Na维中间矩阵VC， vi表示NK×NK维酉矩阵UNK×NK中的第i列矢量，且vi为NK×1维列矢量；
根据NK×Na维中间矩阵VC，计算得到三维异构阵机载雷达的杂波协方差矩阵Rc，其定义式为：上标H表示共轭转置操作；
定义目标加杂波的信号矩阵为Vs，s表示目标的空时导向矢量，并计算得到目标加杂波的协方差矩阵为Rs， 进而计算得到目标加杂波的协方差矩阵Rs的行列式|Rs|，
3.如权利要求2所述的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，其特征在于，所述s表示目标的空时导向矢量，其定义表达式为：s＝[s1,s2,...,si',...,sNK]T，上标T为转置操作，i'∈{1,2,…,NK}，si'表示目标的空时导向矢量s中第i'个元素，si'＝anbm；
所述s表示目标的空时导向矢量，其计算表达式为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mi>a</mi>
<mo>&CircleTimes;</mo>
<mi>b</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mi>b</mi>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mi>b</mi>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mi>b</mi>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<mi>b</mi>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mi>b</mi>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msub>
<mi>b</mi>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>b</mi>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mi>n</mi>
</msub>
<msub>
<mi>b</mi>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mo>&rsqb;</mo>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，表示kronecker积，s为NK×1维列矢量；a表示目标的空域导向矢量，a为N×1维列矢量， 上标T表示转置操作，an表示目标的空域导向矢量a中第n个元素，N为三维圆柱阵雷达包含的阵元个数，与目标的空域导向矢量a中包含的元素个数取值相等；b表示目标的时域导向矢量，b为K×1维列矢量，bm表示目标的时域导向矢量b中第m个元素，K为每个相干处理间隔内三维圆柱阵雷达发射的脉冲个数，与目标的时域导向矢量b中包含的元素个数取值相等；fs表示目标的空间频率，φ表示设定的目标的方位角，θ表示设定的目标的俯仰角，cos表示求余弦操作，sin表示求正弦操作，λ表示三维圆柱阵雷达发射信号波长。
4.如权利要求3所述的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，其特征在于，步骤3的子步骤为：
(1)将三维异构阵机载雷达的杂波加噪声协方差矩阵表示为R，且R为NK×NK的二维矩阵，表示高斯白噪声的功率，Rc为三维异构阵机载雷达的杂波协方差矩阵，INK表示N×K维全1矩阵；
三维异构阵机载雷达得杂波协方差矩阵Rc的计算表达式为： vi表示NK×NK维酉矩阵UNK×NK中的第i列矢量，且vi为NK×1维列矢量；qi表示NK×NK维酉矩阵UNK×NK中第i列矢量的杂波功率，上标H表示共轭转置操作；
定义Na×1维矢量并对Na×1维矢量进行对角化，即进而得到Na×Na维二维矩阵Q，diag表示对角化操作；计算三维异构阵机载雷达的杂波加噪声协方差矩阵的简化形式
<mrow>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<msubsup>
<mi>QV</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mi>Q</mi>
</mrow>
<mo>)</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
然后计算三维异构阵机载雷达的杂波加噪声协方差矩阵的简化形式的逆
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mi>Q</mi>
</mrow>
<mo>)</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>&lsqb;</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mi>K</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mi>K</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mi>Q</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<msubsup>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mi>K</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mi>K</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>Q</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
假定NK×NK维酉矩阵UNK×NK中每一列矢量的杂波功率都远大于噪声功率，即 则得到三维异构阵机载雷达的杂波加噪声协方差矩阵的简化形式的逆的简化式 进而计算三维异构阵机载雷达的空时自适应处理的最优权Wopt，
(2)将目标的空时导向矢量s分解为sc和s⊥两个正交的子空间，s＝sc+s⊥，sc和s⊥正交，并且互为补空间，即sc∪s⊥＝s， 表示空集，sc表示杂波子空间，s⊥表示杂波子空间sc的正交补空间，VC表示NK×Na维中间矩阵，INK表示N×K维全1矩阵，上标H表示共轭转置，K为每个相干处理间隔内三维圆柱阵雷达发射的脉冲个数，N为三维圆柱阵雷达包含的阵元个数，Na表示NK×L的二维导向矢量矩阵BNK×L的秩，L表示三维圆柱阵雷达接收到的杂波包含的距离门总个数；
所述sc表示杂波子空间，且P表示正交投影矩阵，上标H表示共轭转置操作，上标-1表示求逆操作；由此计算得到目标的空时导向矢量s与杂波子空间sc之间夹角α的余弦值cos(α)，其中，|·|表示求绝对值操作，||s||2表示目标的空时导向矢量s的二范数，进而计算得到目标的空时导向矢量s与杂波子空间sc之间夹角α的余弦平方值|cos(α)|2，
<mrow>
<msup>
<mrow>
<mo>|</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msup>
<mi>s</mi>
<mi>H</mi>
</msup>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mi>s</mi>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mi>N</mi>
<mi>K</mi>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mi>s</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mi>K</mi>
</mrow>
</mfrac>
<msup>
<mi>s</mi>
<mi>H</mi>
</msup>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mi>s</mi>
<mo>;</mo>
</mrow>
(3)根据目标的空时导向矢量s与杂波子空间sc之间夹角α的余弦平方值|cos(α)|2，计算得到三维异构阵机载雷达的输出信杂噪比SCNRout，
s表示目标的空时导向矢量，表示三维异构阵机载雷达的杂波加噪声协方差矩阵的简化形式的逆的简化式，SNR表示目标的信噪比， 表示目标的信号功率，代表高斯白噪声的功率，N为三维圆柱阵雷达包含的阵元个数，K表示每个相干处理间隔内三维圆柱阵雷达发射的脉冲个数。
5.如权利要求4所述的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，其特征在于，在步骤4中，将所述三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式记为其表达式为：
<mrow>
<mi>S</mi>
<mi>C</mi>
<mi>N</mi>
<msub>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>R</mi>
<mi>c</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
其中，Rc为三维异构阵机载雷达的杂波协方差矩阵，Rs表示目标加杂波的协方差矩阵，|Rs|表示目标加杂波的协方差矩阵Rs的行列式，|Rc|表示三维异构阵机载雷达的杂波协方差矩阵Rc的行列式，上标H表示共轭转置操作，上标-1表示求逆操作，s表示目标的空时导向矢量，K为每个相干处理间隔内三维圆柱阵雷达发射的脉冲个数，N为三维圆柱阵雷达包含的阵元个数，SNR表示目标的信噪比， 表示目标的信号功率，代表高斯白噪声的功率。
6.如权利要求5所述的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，其特征在于，在步骤5中，所述二进制挑选向量y为NK×1维列矢量，N为三维圆柱阵雷达包含的阵元个数，K为每个相干处理间隔CPI内三维圆柱阵雷达发射的脉冲个数，NK表示一共有NK个阵元脉冲对；二进制挑选向量y中包含NK个元素，每个元素的取值分别为0或1，1代表选择了对应的阵元脉冲，0代表没有选择对应的阵元脉冲；
将所述三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式记为
<mrow>
<mi>S</mi>
<mi>C</mi>
<mi>N</mi>
<msub>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>R</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>=</mo>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msubsup>
<mi>V</mi>
<mi>s</mi>
<mi>H</mi>
</msubsup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中，Rc(y)表示被挑选后的杂波协方差矩阵，Rs(y)表示被挑选后的杂波加目标协方差矩阵，上标H表示共轭转置操作，VC表示NK×Na维中间矩阵，Vs表示目标加杂波的信号矩阵。
7.如权利要求6所述的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，其特征在于，步骤6的子步骤为：
在三维异构阵机载雷达的输出信杂噪比SCNRout的行列式表示形式的重写式取值最大的情况下进行阵元挑选，使得二进制挑选向量y里面有M个1，NK-M个0；
6.1令二进制挑选向量为y，并令将二进制挑选向量y的初始值全部为1，记为二进制挑选向量初始值y(0)，即y(0)＝[1,1,...,1]NK，二进制挑选向量初始值y(0)的下标集为η(0)，η(0)＝[1,2,...,NK]；
6.2定义外层循环变量k，k为第k次迭代，k的初始值为1，k∈{1,2,…,NK-M}，M<NK；定义内层循环变量g，g为第g次迭代，g随着k的变化而变化,g∈{1,2,…,NK-k+1},g的初始值为1；当k取值为1时，g取值从1循环到NK；当k取值为2时，g取值从1循环到NK-1；当k取值为M时，g取值从1循环到NK-M+1；外层循环变量每循环一次，就会将第k次迭代后的二进制挑选向量y(k)中的一个1置为0，当外层循环变量循环M次时，就会将第M次迭代后的二进制挑选向量y(M)中的M个1都置为0；
6.3根据第k次迭代后的二进制挑选向量y(k),得到第k次迭代后的二进制挑选向量y(k)的下标集η(k),所述第k次迭代后的二进制挑选向量y(k)的下标集η(k)中与第k次迭代后的二进制挑选向量y(k)中元素为0的对应位置处元素为0,y(k)的下标集η(k)中其余位置处的元素按照自然数顺序从1排列至NK-k+1；所述第k次迭代后的二进制挑选向量y(k)为第k-1次迭代后的二进制挑选向量y(k-1)；
6.4令第g次迭代后下标集η(k)中与g取值相等的元素在η(k)中的位置为h，即η(k)(h)＝g，h∈{1,2,...,NK-k+1}，然后将第k次迭代后的二进制挑选向量y(k)中位置h处的元素置为0，即y(k)(h)＝0，并将第k次迭代后的二进制挑选向量y(k)中位置h处的元素置为0后得到的向量,记为第g次迭代后在位置h处置零的向量且满足即向量中位置h处的元素为0；
然后计算第g次迭代后在位置h处置零的向量的杂波协方差矩阵和第g次迭代后在位置h处置零的向量的杂波加目标协方差矩阵以及分别计算杂波协方差矩阵的行列式和杂波加目标协方差矩阵的行列式进而计算杂波协方差矩阵的行列式和杂波加目标协方差矩阵的行列式的比值μk(g)，其表达式分别为:
<mrow>
<msub>
<mi>R</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>y</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>V</mi>
<mi>c</mi>
<mi>H</mi>
</msubsup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>y</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>V</mi>
<mi>c</mi>
</msub>
<mo>,</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>y</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>V</mi>
<mi>s</mi>
<mi>H</mi>
</msubsup>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>y</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<msup>
<mi>&mu;</mi>
<mi>k</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>g</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>R</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>y</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>y</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
其中，杂波协方差矩阵的行列式和杂波加目标协方差矩阵的行列式的维数都是NK×NK,diag表示对角化操作；然后,将杂波协方差矩阵的行列式和杂波加目标协方差矩阵的行列式的比值μk(g),记为第k次迭代后经过第g次迭代得到的第g个比值；
6.5令g分别取1至NK-k+1，重复执行6.3和6.4，进而分别得到第k次迭代后经过第1次迭代得到的第1个比值至第k次迭代后经过第NK-k+1次迭代得到的第NK-k+1个比值，记为第k次迭代后得到的NK-k+1个比值,比较第k次迭代后得到的NK-k+1个比值，得到比值最小值，并得到比值最小值时g的对应取值在第k次迭代后的二进制挑选向量y(k)的下标集η(k)中的位置，记为hmin，hmin∈{1,2,…,NK-k+1}，然后将第k次迭代后的二进制挑选向量y(k)中位置hmin处的1置为0，即y(k)(hmin)＝0，进而得到第k次迭代后的二进制挑选向量y(k)，所述第k次迭代后的二进制挑选向量y(k)中有k个0，NK-k个1；
6.6令k加1，返回6.3，直到得到第NK-M次迭代后的二进制挑选向量y(NK-M)，迭代停止，此时得到的第NK-M次迭代后的二进制挑选向量y(NK-M)中NK-M个位置处的值分别已置为0，其余M个位置处的值分别为1，即迭代的最终结果使得第NK-M次迭代后的二进制挑选向量y(NK-M)中有M个1，NK-M个0；
然后计算得到对角矩阵Y，Y＝diag(y(NK-M))，对角矩阵Y为NK行NK列的二维矩阵，剔除对角矩阵Y中全部为0的列，并将对角矩阵Y中剔除全部为0的列后的矩阵，记为三维异构阵机载雷达的稀疏降维矩阵Z，Z是NK行M列的二维矩阵。
8.如权利要求7所述的一种基于信杂噪比最大的雷达三维异构阵稀疏重构方法，其特征在于，在步骤7中，将三维异构阵机载雷达的滤波结果记为out，其表达式为：
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mover>
<mi>W</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>o</mi>
<mi>p</mi>
<mi>t</mi>
</mrow>
<mi>H</mi>
</msubsup>
<mo>&times;</mo>
<mo>(</mo>
<mrow>
<msup>
<mi>Z</mi>
<mi>H</mi>
</msup>
<mi>R</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
</mrow>
<mrow>
<msub>
<mover>
<mi>W</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>o</mi>
<mi>p</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<msup>
<mi>Z</mi>
<mi>H</mi>
</msup>
<mi>s</mi>
</mrow>
<mo>)</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>Z</mi>
<mi>H</mi>
</msup>
<mi>R</mi>
<mi>Z</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>(</mo>
<mrow>
<msup>
<mi>Z</mi>
<mi>H</mi>
</msup>
<mi>s</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>Z</mi>
<mi>H</mi>
</msup>
<mi>R</mi>
<mi>Z</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>(</mo>
<mrow>
<msup>
<mi>Z</mi>
<mi>H</mi>
</msup>
<mi>s</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
其中，R表示三维异构阵机载雷达的杂波加噪声协方差矩阵，s表示目标的空时导向矢量，表示稀疏后的自适应权值，E表示三维异构阵雷达的杂波功率,上标H表示共轭转置操作，上标-1表示求逆操作。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710499046.6A CN107167782B (zh) | 2017-06-27 | 2017-06-27 | 基于信杂噪比最大的雷达三维异构阵稀疏重构方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710499046.6A CN107167782B (zh) | 2017-06-27 | 2017-06-27 | 基于信杂噪比最大的雷达三维异构阵稀疏重构方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107167782A true CN107167782A (zh) | 2017-09-15 |
| CN107167782B CN107167782B (zh) | 2020-04-10 |
Family
ID=59826166
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710499046.6A Active CN107167782B (zh) | 2017-06-27 | 2017-06-27 | 基于信杂噪比最大的雷达三维异构阵稀疏重构方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107167782B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112800497A (zh) * | 2020-12-28 | 2021-05-14 | 西安电子科技大学 | 基于稀疏谱恢复的机载三维异构阵杂波抑制方法 |
| CN114726385A (zh) * | 2022-04-21 | 2022-07-08 | 电子科技大学 | 基于功率估计的卫星导航接收机空域抗干扰方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6292592B1 (en) * | 1998-10-19 | 2001-09-18 | Raytheon Company | Efficient multi-resolution space-time adaptive processor |
| CN103018727A (zh) * | 2011-09-27 | 2013-04-03 | 中国科学院电子学研究所 | 一种基于样本训练的机载雷达非平稳杂波抑制方法 |
| CN103399303A (zh) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | 机载雷达抗密集转发式欺骗干扰方法与系统 |
| CN104237883A (zh) * | 2014-09-15 | 2014-12-24 | 西安电子科技大学 | 一种采用稀疏表示的机载雷达空时自适应处理方法 |
| CN105487054A (zh) * | 2015-11-09 | 2016-04-13 | 大连大学 | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 |
-
2017
- 2017-06-27 CN CN201710499046.6A patent/CN107167782B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6292592B1 (en) * | 1998-10-19 | 2001-09-18 | Raytheon Company | Efficient multi-resolution space-time adaptive processor |
| CN103018727A (zh) * | 2011-09-27 | 2013-04-03 | 中国科学院电子学研究所 | 一种基于样本训练的机载雷达非平稳杂波抑制方法 |
| CN103399303A (zh) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | 机载雷达抗密集转发式欺骗干扰方法与系统 |
| CN104237883A (zh) * | 2014-09-15 | 2014-12-24 | 西安电子科技大学 | 一种采用稀疏表示的机载雷达空时自适应处理方法 |
| CN105487054A (zh) * | 2015-11-09 | 2016-04-13 | 大连大学 | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112800497A (zh) * | 2020-12-28 | 2021-05-14 | 西安电子科技大学 | 基于稀疏谱恢复的机载三维异构阵杂波抑制方法 |
| CN112800497B (zh) * | 2020-12-28 | 2023-08-11 | 西安电子科技大学 | 基于稀疏谱恢复的机载三维异构阵杂波抑制方法 |
| CN114726385A (zh) * | 2022-04-21 | 2022-07-08 | 电子科技大学 | 基于功率估计的卫星导航接收机空域抗干扰方法 |
| CN114726385B (zh) * | 2022-04-21 | 2023-02-24 | 电子科技大学 | 基于功率估计的卫星导航接收机空域抗干扰方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107167782B (zh) | 2020-04-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110275166B (zh) | 基于admm的快速稀疏孔径isar自聚焦与成像方法 | |
| CN104237883B (zh) | 一种采用稀疏表示的机载雷达空时自适应处理方法 | |
| CN103713288B (zh) | 基于迭代最小化稀疏贝叶斯重构线阵sar成像方法 | |
| CN108008385B (zh) | 基于稀疏贝叶斯学习的干扰环境isar高分辨成像方法 | |
| CN105403856A (zh) | 基于嵌套式最小冗余阵列的波达方向估计方法 | |
| CN103744076B (zh) | 基于非凸优化的mimo雷达动目标检测方法 | |
| CN106405548A (zh) | 基于多任务贝叶斯压缩感知的逆合成孔径雷达成像方法 | |
| CN101887117B (zh) | 基于三迭代的机载mimo雷达空时降维自适应处理方法 | |
| CN108828547A (zh) | 基于深度神经网络的米波雷达低仰角测高方法 | |
| CN106468770A (zh) | K分布杂波加噪声下的近最优雷达目标检测方法 | |
| CN110109050A (zh) | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 | |
| CN106680815A (zh) | 基于张量稀疏表示的mimo雷达成像方法 | |
| CN107167783A (zh) | 一种共形阵列杂波协方差矩阵的稀疏重构方法 | |
| CN106772304A (zh) | 基于空域多级分解的机载mimo雷达后多普勒自适应处理方法 | |
| CN104345299A (zh) | 基于简化ec的机载mimo雷达空时自适应处理方法 | |
| CN107576947A (zh) | 基于时间平滑的l型阵对相干信源二维波达方向估计方法 | |
| CN103399308B (zh) | 主瓣和旁瓣干扰背景下雷达目标角度快速估计方法 | |
| CN103792523B (zh) | 基于张量积的uhf波段多通道雷达径向速度检测方法 | |
| CN109298382A (zh) | 一种基于期望极大算法的非均匀直线阵波达方向角估计方法 | |
| Kang et al. | SAR imaging based on deep unfolded network with approximated observation | |
| CN107167782A (zh) | 基于信杂噪比最大的雷达三维异构阵稀疏重构方法 | |
| CN107064896A (zh) | 基于截断修正sl0算法的mimo雷达参数估计方法 | |
| Shixin et al. | Joint 2D DOA and Doppler frequency estimation for L-shaped array using compressive sensing | |
| CN104035078A (zh) | 一种基于阵元阶数递推的降维空时自适应权值计算方法 | |
| CN104459660A (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 |