基于空时采样矩阵的机载MIMO雷达自适应杂波抑制方法

基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法
技术领域
1.本发明属于信号处理技术领域，更具体地，涉及一种基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法。


背景技术：

2.mimo雷达是将无线通信系统中的多输入多输出技术引入到雷达领域，并和数字阵列技术相结合的一种新体制雷达。mimo雷达的概念于2003年由美国林肯实验室的bliss和forsythe首次提出。mimo 雷达中多输入通常指多个阵元同时发射非相干的波形，多输出指多个阵元同时接收，用匹配滤波器组来分离回波信号中的各发射分量。相反，传统相控阵雷达多个阵元同时发射相同的波形，在空间合成一个发射信号，因此可看作为单输入多输出雷达，即simo。根据发射阵元和接收阵元的间距大小，可将mimo雷达分为分布式mimo雷达（或统计mimo雷达）和集中式mimo雷达（或相干mimo雷达）两类。分布式mimo雷达中收发天线各阵元相距较远，使得各阵元可分别从不同的视角观察目标，可获得空间分集得益，克服目标rcs闪烁效应，提高雷达目标探测性能。集中式mimo雷达中收发天线各阵元相距较近，对目标的视角近似相同，每个阵元发射不同的信号波形，从而获得波形分集得益。集中式mimo雷达可以产生大的虚拟孔径，提高雷达的空间分辨率。本发明的研究对象是集中式mimo体制机载雷达，简称为机载mimo雷达。
3.传统的降维stap方法和降秩stap方法同样适用于机载mimo雷达，但是存在两方面问题：一是采用mimo体制后由于波形分集导致杂波自由度显著增大，杂波抑制更具挑战性；二是系统自由度的增加导致运算量和样本需求量增大。因此，如何在降低运算量和样本需求的同时确保杂波抑制性能是机载mimo雷达需解决的关键技术之一。ka-stap方法是一种利用数字地形图、地表覆盖数据和人工道路与建筑物位置信息等先验知识增强传统stap收敛性的方法。ka-stap方法分为两类：一是间接利用先验知识，通过先验知识对训练样本进行选取，以保证选取的训练样本满足独立同分布条件；二是直接利用先验知识，将利用回波数据估计出的杂波协方差矩阵和通过先验知识构造的杂波协方差矩阵相结合，通过先验知识对估计的杂波协方差矩阵进行适当修正。当先验知识存在误差时，ka-stap方法的杂波抑制性能改善有限。
4.因此，急需发明一种有效的基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，以提高机载mimo雷达对运动目标的检测性能。


技术实现要素：

5.为此，本发明提供一种基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，用以克服现有技术中存在的问题。
6.为实现上述目的，本发明提供一种基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，包括如下步骤，步骤1：针对机载mimo雷达回波数据，利用空时滑窗处理估计杂波功率；
步骤2：结合空时等效阵列形式和步骤1估计的杂波功率计算矩阵g；步骤3：根据空时等效阵列形式计算空时采样矩阵b；步骤4：利用步骤2和3的结果构造杂波噪声协方差矩阵r1；步骤5：利用回波数据估计杂波噪声协方差矩阵r2；步骤6：根据步骤4和5的结果形成最终的杂波噪声协方差矩阵r；步骤7：根据步骤6的结果形成空时自适应权值，并进行杂波抑制处理；设定所述机载mimo雷达发射天线包含m个阵元，阵元间距为d
t
，接收天线包含n个阵元，阵元间距为dr，杂波自由度为rc，距离模糊次数为nr，发射阵元间距与接收阵元间距比为γ，归一化多普勒频率与空间频率比为β，相参处理脉冲数为k。
7.进一步地，所述步骤1中利用空时滑窗处理得到待检测距离单元回波对应的空时杂波谱为（1）其中表示空时滑窗后子孔径级空时导向矢量，表示通过空时滑窗后的数据估计得到的杂波协方差矩阵，m1、n1和k1分别表示发射空域、接收空域和时域子孔径长度，表示取共轭转置，对上述杂波谱在空时平面上求平均后得到杂波功率估计值，数学表达式为（2）其中表示第(n,m)个单元对应的杂波功率值，ks、kd分别表示空域和时域划分的单元数目。
8.进一步地，在所述步骤2中，当空时等效阵列为均匀线阵时，各独立采样点的位置位于半波长的整数倍处，此时（3）其中sinc(
·
)表示辛格函数，矩阵g的维数为rc×
rc，u=1,2,
…
,rc，v=1,2,
…
,rc，φ
l
表示第l个距离单元对应的俯仰角，λ表示雷达工作波长；当空时等效阵列为密集非均匀阵时，对空时等效阵列中等效阵元间距小于半波长的阵元进行凝聚处理，此时（4）其中集合表示凝聚处理后的独立采样点位置；
当空时等效阵列为稀疏非均匀阵时，将各独立采样点的位置选取为等效阵元位置，矩阵g的求取方式同式（4）。
9.进一步地，在所述步骤3中，当空时等效阵列为均匀线阵和密集非均匀阵时，空时采样矩阵b的各元素为（5）其中m=1,2,
…
,m，n=1,2,
…
,n，k=1,2,
…
,k；当空时等效阵列为稀疏非均匀阵时，矩阵b为由维数为mn
×
mn的矩阵a组成的块对角矩阵，其形式为（6）其中矩阵a的各元素为（7）。
10.进一步地，在所述步骤4中，利用步骤2和3的结果构造杂波噪声协方差矩阵为（8）其中σ2表示噪声功率，i表示维数为mnk
×
mnk的单位矩阵。
11.进一步地，在所述步骤5中，利用回波数据估计杂波噪声协方差矩阵为（9）其中x
l
表示第l个距离单元的机载mimo雷达回波数据，l表示训练样本数目。
12.进一步地，在所述步骤6中，根据步骤4和5的结果形成最终的杂波噪声协方差矩阵为（10）。
13.进一步地，在所述步骤7中，形成的空时自适应权值为（11）其中，经过杂波抑制处理后的输出信号为（12）其中x表示机载mimo雷达回波数据。
14.总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有如下有益效果：（1）本发明通过提供一种基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，首先利用杂波自由度信息和系统参数分别形成空时采样矩阵b和杂波子空间矩阵g，基于待检测单元数据通过空时滑窗处理估计杂波功率，从而构造得到杂波噪声协方差矩阵；其次利用少量训练样本数据估计杂波噪声协方差矩阵；再次将构造得到的杂波噪声协方差矩阵与利用少量训练样本数据估计的杂波噪声协方差矩阵相结合，增强方法的稳健性；最后形成空时自适应权值，完成杂波抑制处理，提高了机载mimo雷达的运动目标检测性能。
15.（2）本发明当β为整数时不需要训练样本，仅通过待检测单元数据即可完成对杂波噪声协方差矩阵的构造，因此适用于极端非均匀杂波环境；（3）本发明当β为非整数且在均匀杂波环境下时，本发明仅需要少量训练样本即可实现快速收敛；（4）本发明同时利用了系统参数信息和回波数据，具有强的稳健性；（5）本发明的空时采样矩阵和独立采样点位置可离线计算，因此运算量显著降低，尤其适用于高系统自由度的mimo体制机载雷达。
附图说明
16.图1为本发明所述基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法的结构框图；图2为本发明所述基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法的流程示意图。
具体实施方式
17.为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
18.请参阅图1-2所示，本发明通过提供一种基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，包括如下步骤：步骤1：针对机载mimo雷达回波数据，利用空时滑窗处理估计杂波功率；杂波功率估计单元利用空时滑窗处理得到待检测距离单元回波对应的空时杂波谱为（1）其中表示空时滑窗后子孔径级空时导向矢量，表示通过空时滑窗后的数据估计得到的杂波协方差矩阵，m1、n1和k1分别表示发射空域、接收空域和时域子孔径长度，表示取共轭转置，对上述杂波谱在空时平面上求平均后得到杂波功率估
计值，数学表达式为（2）其中表示第(n,m)个单元对应的杂波功率值，ks、kd分别表示空域和时域划分的单元数目。
19.步骤2：结合空时等效阵列形式和步骤1估计的杂波功率计算矩阵g；矩阵g计算单元中当空时等效阵列为均匀线阵时，各独立采样点的位置位于半波长的整数倍处，此时（3）其中sinc(
·
)表示辛格函数，矩阵g的维数为rc×
rc，u=1,2,
…
,rc，v=1,2,
…
,rc，φ
l
表示第l个距离单元对应的俯仰角，λ表示雷达工作波长；当空时等效阵列为密集非均匀阵时，对空时等效阵列中等效阵元间距小于半波长的阵元进行凝聚处理，此时（4）其中集合表示凝聚处理后的独立采样点位置；当空时等效阵列为稀疏非均匀阵时，将各独立采样点的位置选取为等效阵元位置，矩阵g的求取方式同式（4）。
20.步骤3：根据空时等效阵列形式计算空时采样矩阵b；空时采样矩阵b单元中当空时等效阵列为均匀线阵和密集非均匀阵时，空时采样矩阵b的各元素为（5）其中m=1,2,
…
,m，n=1,2,
…
,n，k=1,2,
…
,k；当空时等效阵列为稀疏非均匀阵时，矩阵b为由维数为mn
×
mn的矩阵a组成的块对角矩阵，其形式为（6）其中矩阵a的各元素为
（7）。
21.步骤4：利用步骤2和3的结果构造杂波噪声协方差矩阵；杂波噪声协方差矩阵构造单元利用步骤2和3的结果构造杂波噪声协方差矩阵为（8）其中σ2表示噪声功率，i表示维数为mnk
×
mnk的单位矩阵。
22.步骤5：利用回波数据估计杂波噪声协方差矩阵；杂波噪声协方差矩阵估计单元利用回波数据估计杂波噪声协方差矩阵为（9）其中x
l
表示第l个距离单元的机载mimo雷达回波数据，l表示训练样本数目。
23.步骤6：根据步骤4和5的结果形成最终的杂波噪声协方差矩阵；最终杂波噪声协方差矩阵计算单元根据步骤4和5的结果形成最终的杂波噪声协方差矩阵为（10）步骤7：根据步骤6的结果形成空时自适应权值，并进行杂波抑制处理；杂波抑制处理单元形成的空时自适应权值为（11）其中，经过杂波抑制处理后的输出信号为（12）其中x表示机载mimo雷达回波数据。
24.本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，其特征在于，包括如下步骤，步骤1：针对机载mimo雷达回波数据，利用空时滑窗处理估计杂波功率；步骤2：结合空时等效阵列形式和步骤1估计的杂波功率计算矩阵g；步骤3：根据空时等效阵列形式计算空时采样矩阵b；步骤4：利用步骤2和3的结果构造杂波噪声协方差矩阵r1；步骤5：利用回波数据估计杂波噪声协方差矩阵r2；步骤6：根据步骤4和5的结果形成最终的杂波噪声协方差矩阵r；步骤7：根据步骤6的结果形成空时自适应权值，并进行杂波抑制处理；设定所述机载mimo雷达发射天线包含m个阵元，阵元间距为d
t
，接收天线包含n个阵元，阵元间距为d
r
，杂波自由度为r
c
，距离模糊次数为n
r
，发射阵元间距与接收阵元间距比为γ，归一化多普勒频率与空间频率比为β，相参处理脉冲数为k。2.根据权利要求1所述的基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，其特征在于，所述步骤1中利用空时滑窗处理得到待检测距离单元回波对应的空时杂波谱为（1）其中表示空时滑窗后子孔径级空时导向矢量，表示通过空时滑窗后的数据估计得到的杂波协方差矩阵，m1、n1和k1分别表示发射空域、接收空域和时域子孔径长度，表示取共轭转置，对上述杂波谱在空时平面上求平均后得到杂波功率估计值，数学表达式为（2）其中表示第(n,m)个单元对应的杂波功率值，k
s
、k
d
分别表示空域和时域划分的单元数目。3.根据权利要求1所述的基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，其特征在于，在所述步骤2中，当空时等效阵列为均匀线阵时，各独立采样点的位置位于半波长的整数倍处，此时（3）其中sinc(
·
)表示辛格函数，矩阵g的维数为r
c
×
r
c
，u=1,2,
…
,r
c
，v=1,2,
…
,r
c
，φ
l
表示第l个距离单元对应的俯仰角，λ表示雷达工作波长；当空时等效阵列为密集非均匀阵时，对空时等效阵列中等效阵元间距小于半波长的阵
元进行凝聚处理，此时（4）其中集合表示凝聚处理后的独立采样点位置；当空时等效阵列为稀疏非均匀阵时，将各独立采样点的位置选取为等效阵元位置，矩阵g的求取方式同式（4）。4.根据权利要求1所述的基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，其特征在于，在所述步骤3中，当空时等效阵列为均匀线阵和密集非均匀阵时，空时采样矩阵b的各元素为（5）其中m=1,2,
…
,m，n=1,2,
…
,n，k=1,2,
…
,k；当空时等效阵列为稀疏非均匀阵时，矩阵b为由维数为mn
×
mn的矩阵a组成的块对角矩阵，其形式为（6）其中矩阵a的各元素为（7）。5.根据权利要求1所述的基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，其特征在于，在所述步骤4中，利用步骤2和3的结果构造杂波噪声协方差矩阵为（8）其中σ2表示噪声功率，i表示维数为mnk
×
mnk的单位矩阵。6.根据权利要求1所述的基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，其特征在于，在所述步骤5中，利用回波数据估计杂波噪声协方差矩阵为（9）其中x
l
表示第l个距离单元的机载mimo雷达回波数据，l表示训练样本数目。7.根据权利要求1所述的基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，其特征在于，在所述步骤6中，根据步骤4和5的结果形成最终的杂波噪声协方差矩阵为
（10）。8.根据权利要求1所述的基于空时采样矩阵的机载mimo雷达自适应杂波抑制方法，其特征在于，在所述步骤7中，形成的空时自适应权值为（11）其中，经过杂波抑制处理后的输出信号为（12）其中x表示机载mimo雷达回波数据。

技术总结
本发明公开了一种基于空时采样矩阵的机载MIMO雷达自适应杂波抑制方法，包括，步骤1：针对机载MIMO雷达回波数据，利用空时滑窗处理估计杂波功率；步骤2：结合空时等效阵列形式和步骤1估计的杂波功率计算矩阵G；步骤3：根据空时等效阵列形式计算空时采样矩阵B；步骤4：利用步骤2和3的结果构造杂波噪声协方差矩阵；步骤5：利用回波数据估计杂波噪声协方差矩阵；步骤6：根据步骤4和5的结果形成最终的杂波噪声协方差矩阵；步骤7：根据步骤6的结果形成空时自适应权值，并进行杂波抑制处理。并进行杂波抑制处理。并进行杂波抑制处理。


技术研发人员：谢文冲 熊元燚 王永良 柳成荫 陈威 田步秋
受保护的技术使用者：中国人民解放军空军预警学院
技术研发日：2022.04.24
技术公布日：2022/5/25