CN106772253A - 一种非均匀杂波环境下的雷达杂波抑制方法 - Google Patents

Abstract

本发明公开了一种非均匀杂波环境下的雷达杂波抑制方法，其主要思路为：确定机载雷达，并获取机载雷达的原始雷达回波数据，然后确定机载雷达的时域滑窗孔径；计算机载雷达第k个多普勒通道的最终空时导向矢量，并计算第l个距离门的扩展雷达回波数据，然后依次计算机载雷达第k个多普勒通道、第l个距离门对应的自适应权值和空时处理过程后第k个多普勒通道、第l个距离门的对应输出矢量和机载雷达第k个多普勒通道、第l个距离门对应的滤波输出，并计算机载雷达第l个距离门对应的滤波输出；依次令k∈{1,2,…,K}，令l∈{1,2,…,L}，进而得到杂波抑制后机载雷达的距离—多普勒雷达回波数据。

Description

一种非均匀杂波环境下的雷达杂波抑制方法

技术领域

本发明属于雷达杂波抑制技术领域，特别涉及一种非均匀杂波环境下的雷达杂波抑制方法，是一种基于非均匀杂波环境下的级联二维空时处理方法，适用于斜侧阵机载雷达杂波抑制或非正侧阵机载雷达杂波抑制。

背景技术

机载预警雷达已经成为现代战争的预警情报、战场监视和作战指挥系统的重要组成部分，且能够显著增加低空突防目标的探测距离，大大提高武器系统的效能，使得在现在战争中扮演着越来越重要的角色，并受到各军事大国的高度重视。机载预警雷达的波束指向一般是水平或稍下方向，使得机载预警雷达不可避免地会接收到地杂波。由于技术水平限制，使得阵列机载雷达天线的旁瓣电平较高，且阵列机载雷达回波中地杂波较强，可达90dB，远强于目标回波；此外，由于飞机平台的运动，地杂波的多普勒范围大大展宽，使得可供检测的多普勒范围大大降低，有很大一部分目标都“湮没”在地杂波中。因此，为了可靠检测远距离目标，首先要解决机载雷达的杂波抑制问题。

空时自适应处理(STAP)技术的诞生追溯到20世纪70年代，在杂波协方差矩阵和目标信号分别确知的条件下，Brennan和Reed等在期刊Aerospace and Electronic Systems(AES)上，提出了全空时自适应处理(STAP)的概念和理论，其思想是将阵列信号处理的基本原理推广到由脉冲和阵元采样的两维场中；全空时自适应处理(STAP)能够取得比较理想的机载雷达杂波抑制效果，但机载雷达杂波协方差矩阵精确已知的条件在工程实现中难以满足，通常情况下是由参考单元估计得到；为了使空时自适应处理(STAP)器的性能损失不超过3dB，空时自适应处理(STAP)器训练样本个数至少应为机载雷达系统自由度的两倍，而且，机载雷达杂波协方差矩阵须从独立同分布的训练样本估计而来，所以全空时自适应处理(STAP)的运算量和设备复杂度令人难以接受。

为了降低运算量，德国的R.Klemm博士于1987年在期刊Signal Processing上提出了辅助通道法，将空时自适应处理(STAP)器维数从NM降至N+M-1，此处N表示机载雷达天线阵元个数，M表示一个相干处理间隔内机载雷达接收的脉冲个数；在没有误差的理想情况下，该方法能够取得比较理想的杂波抑制效果，性能接近最优空时自适应处理(STAP)处理，并且降低了估计杂波协方差矩阵所需要的独立同分布参考单元数目，但在有误差情况下，该方法的杂波抑制性能并不理想。

国内的西安电子科技大学在降维空时自适应处理(STAP)研究方面也展开了大量工作，并提出了局域化处理法，即因子化方法(1DT)和扩展因子化方法(mDT)；该局域化处理法首先利用低旁瓣的多普勒滤波器对机载雷达回波进行局域化处理，然后在空域自适应处理或空时域联合自适应处理将雷达杂波抑制掉；但这种局域化的处理方法依然需要大量满足独立同分布的训练样本来构造协方差矩阵，在非均匀杂波环境下，几乎很难得到大量满足独立同分布要求的样本，所以在非均匀杂波环境下的杂波抑制性能并没有得到大的改善。

发明内容

针对上述现有技术存在的不足，本发明的目的在于提出一种非均匀杂波环境下的雷达杂波抑制方法，该种非均匀杂波环境下的雷达杂波抑制方法能够充分利用待检样本附近的少量样本作为局部训练样本，采用相干积累的形式保证信号增益，同时又保证了足够多的自由度进行雷达杂波的充分抑制。

为达到上述技术目的，本发明采用如下技术方案予以实现。

一种非均匀杂波环境下的雷达杂波抑制方法，包括以下步骤：

步骤1，确定机载雷达，并获取机载雷达的原始雷达回波数据X，然后确定机载雷达的时域滑窗孔径K_t；

步骤2，初始化：令k∈{1,2,…,K}，k表示第k个多普勒通道，K表示机载雷达的原始雷达回波数据X在频率域包含的多普勒总个数，k的初始值为1；令l∈{1,2,…,L}，l表示第l个距离门，L表示机载雷达的原始雷达回波数据X包含的距离门总个数，l的初始值为1；

步骤3，依次计算机载雷达第k个多普勒通道进行空时处理的时域导向矢量S_kt，以及机载雷达第k个多普勒通道的最终空时导向矢量S_k；

步骤4，将第l个距离门的原始雷达回波数据记为x_l，并计算第l个距离门的扩展雷达回波数据Q_l，然后依次计算第l个距离门对应的训练样本数据Train_l、时域滑窗后第l个距离门对应的扩展训练样本数据Train_l′，以及第l个距离门对应的协方差矩阵R_l；

步骤5，计算机载雷达第k个多普勒通道、第l个距离门对应的自适应权值w_kl，进而依次计算空时处理过程后第k个多普勒通道、第l个距离门的对应输出矢量y_kl和机载雷达第k个多普勒通道、第l个距离门对应的滤波输出z_kl；

步骤6，令k加1，依次重复步骤3至步骤5，直到得到机载雷达第K个多普勒通道、第l个距离门对应的滤波输出z_Kl，并将此时得到的机载雷达第1个多普勒通道、第l个距离门对应的滤波输出z_1l至机载雷达第K个多普勒通道、第l个距离门对应的滤波输出z_Kl，作为机载雷达第l个距离门对应的滤波输出z_l，z_l＝[z_1l z_2l … z_kl … z_Kl]^H，上标H表示共轭转置；

步骤7，令l加1,依次重复步骤3至步骤6，直到得到机载雷达第L个距离门对应的滤波输出z_L，并将此时得到的机载雷达第l个距离门对应的滤波输出z_l至机载雷达第L个距离门对应的滤波输出z_L，作为杂波抑制后机载雷达的距离—多普勒雷达回波数据Z，Z＝[z₁ z₂… z_l … z_L]。

本发明与现有技术相比具有以下优点：

第一，本发明中提出的方法在训练样本较少的情况下，能够通过时域滑窗的方式提供充足的自由度来抑制雷达杂波。

第二，使用常规的扩展因子化方法(mDT)估计机载雷达杂波协方差矩阵时，选取的训练样本距离维跨度较大，导致在非均匀杂波环境下，不能很好地估计当前距离门对应的协方差矩阵，而本发明方法使得训练样本选择局域化，修正了扩展因子化方法协方差矩阵估计不准的缺陷，m＝3。

附图说明

下面结合附图和具体实施方式对本发明作进一步详细说明。

图1是本发明的一种非均匀杂波环境下的雷达杂波抑制方法流程图；

图2(a)是使用扩展因子化方法(mDT)进行杂波抑制后得到的距离—多普勒图，其中横坐标为多普勒通道，纵坐标为距离门，m＝3；

图2(b)是使用本发明方法进行杂波抑制后得到的距离—多普勒图，其中横坐标为多普勒通道，纵坐标为距离门；

图3(a)是分别使用扩展因子化方法(mDT)和本发明方法进行杂波抑制处理后1-150号距离门对应的杂波剩余图，其中横坐标表示多普勒通道，纵坐标表示杂波剩余，单位为dB，m＝3；

图3(b)是分别使用扩展因子化方法(mDT)和本发明方法进行杂波抑制处理后150-349号距离门对应的杂波剩余图；其中横坐标表示多普勒通道，纵坐标表示杂波剩余，单位为dB，m＝3。

具体实施方式

参考图1，为本发明的一种非均匀杂波环境下的雷达杂波抑制方法流程图；本发明的一种机载雷达杂波抑制方法，包括以下步骤：

步骤1，确定机载雷达，并获取机载雷达的原始雷达回波数据X，然后确定机载雷达的时域滑窗孔径K_t。

具体地，确定机载雷达，并获取机载雷达的原始雷达回波数据X，所述原始雷达回波数据X为N×M×L维矩阵，N表示机载雷达的天线阵面方位向均匀包含的阵元个数，M表示机载雷达在每个相干处理周期内发射的脉冲个数；L表示机载雷达的原始雷达回波数据X包含的距离门总个数；定义机载雷达的原始雷达回波数据X在频率域包含的多普勒总个数为K，且K得取值为大于M且为2的幂次方整数；然后确定机载雷达的滑窗时域孔径K_t，使其满足下式：

2×(M-K_t+1)≥2×K_t×N

或者K_t为整数。

步骤2，初始化：令k∈{1,2,…,K}，k表示第k个多普勒通道，k的初始值为1；K表示机载雷达的原始雷达回波数据X在频率域包含的多普勒总个数；令l∈{1,2,…,L}，l表示第l个距离门，L表示机载雷达的原始雷达回波数据X包含的距离门总个数，l的初始值为1。

步骤3，依次计算机载雷达第k个多普勒通道进行空时处理的时域导向矢量S_kt，以及机载雷达第k个多普勒通道的最终空时导向矢量S_k。

具体地，将机载雷达空时处理过程中的空域导向矢量记为S_s，其表达式为：

其中，n∈{0,1,2,…,N-1}，表示第n个阵元接收机载雷达的脉冲回波数据时相对第1个阵元的相位偏移，f_s表示机载雷达的空域频率，N表示机载雷达天线阵面方位向均匀包含的阵元个数，上标H表示共轭转置。

将机载雷达第k个多普勒通道进行空时处理的时域导向矢量记为S_kt，其表达式为：

其中，上标H表示共轭转置，m∈{0,1,2,…,M-1}，M表示机载雷达在每个相干处理周期内发射的脉冲个数，表示机载雷达第m个脉冲相对第1个脉冲的相位偏移，f_dk表示机载雷达第k个多普勒通道归一化的时域多普勒频率，且f_dk＝k-1/K；机载雷达的滑窗时域孔径K_t远远小于M，则将机载雷达第k个多普勒通道进行空时处理的时域导向矢量S_kt中第1到第K_t个元素，作为机载雷达第k个多普勒通道进行空时处理的优化时域导向矢量S′_kt，其表达式为：

S′_kt＝S_kt[1:K_t]

其中，S_kt[1:K_t]表示取机载雷达第k个多普勒通道进行空时处理的时域导向矢量S_kt中第1到第K_t个元素构成的矢量；进而计算得到机载雷达第k个多普勒通道的最终空时导向矢量S_k，S_s表示机载雷达空时处理过程中的空域导向矢量，表示Kronecker乘积。

步骤4，将第l个距离门的原始雷达回波数据记为x_l，并计算第l个距离门的扩展雷达回波数据Q_l，然后依次计算第l个距离门对应的训练样本数据Train_l、时域滑窗后第l个距离门对应的扩展训练样本数据Train_l′，以及第l个距离门对应的协方差矩阵R_l。

具体的，定义一种滑窗处理过程(SlideWinProc)如下：

首先，将第l个距离门的原始雷达回波数据记为x_l，其维度为N×M，表达式为：

x_l＝[x₁′_l x′_2l … x′_ml … x′_Ml]

其中，x′_ml表示第m个脉冲、第l个距离门的原始雷达回波数据，且是N×1维列矢量；m＝1,2,…,M，l＝1,2,…,L；由于机载雷达的时域滑窗孔径K_t满足且K_t为整数，则任意取一个满足条件的整数r，作为机载雷达的时域滑窗孔径值，进而对第l个距离门的原始雷达回波数据x_l进行孔径值为r的时域滑窗，得到第l个距离门的扩展雷达回波数据Q_l，维度为N_r×V，其表达式为：

Q_l＝[q_1l q_2l … q_vl … q_Vl]

其中，N_r＝N×r，V表示第l个距离门的扩展雷达回波数据Q_l包含的样本个数，V＝M-r+1；q_vl表示第v个样本、第l个距离门的扩展雷达回波数据，q_vl＝[x′_vl x′_(v+1)l … x′_(v+r-1)l]^H，x′_vl表示第v个样本、第l个距离门的原始雷达回波数据，上标H表示共轭转置，v＝1,2,…,V；所述第v个样本与对机载雷达的原始雷达回波数据X进行滑窗的起始脉冲编号对应。

根据第l个距离门的原始雷达回波数据x_l，分别将第l-1个距离门的原始雷达回波数据x_l-1、将第l+1个距离门的原始雷达回波数据x_l+1，进而得到第l个距离门对应的训练样本数据Train_l，Train_l＝[x_l-1 x_l+1]，其维度为N×M×2。

然后，对第l个距离门对应的训练样本数据Train_l中包含的所有元素分别进行孔径值为r的时域滑窗，得到时域滑窗后第l个距离门对应的扩展训练样本数据Train_l′，维度为N_r×2V，N_r＝N×r；Train_l′＝[Q_l-1 Q_l+1]，Q_l-1表示第l-1个距离门的扩展雷达回波数据，Q_l+1表示第l+1个距离门的扩展雷达回波数据。

最后，根据时域滑窗后第l个距离门对应的扩展训练样本数据Train_l′，计算第l个距离门对应的协方差矩阵R_l，其表达式为：

其中，上标H表示共轭转置。

步骤5，计算机载雷达第k个多普勒通道、第l个距离门对应的自适应权值w_kl，进而依次计算空时处理过程后第k个多普勒通道、第l个距离门的对应输出矢量y_kl和机载雷达第k个多普勒通道、第l个距离门对应的滤波输出z_kl。

具体地，所述机载雷达第k个多普勒通道、第l个距离门对应的自适应权值w_kl，其计算式为：

其中，R_l表示第l个距离门对应的协方差矩阵，S_k表示机载雷达第k个多普勒通道的最终空时导向矢量，上标H表示共轭转置，上标-1表示求逆操作。

将第l个距离门的扩展雷达回波数据Q_l，作为空时处理过程中第l个距离门的输入数据，然后将第k个多普勒通道、第l个距离门对应的自适应权值记为w_kl，进而计算空时处理过程后第k个多普勒通道、第l个距离门的对应输出矢量y_kl，其表达式为：

其中，上标H表示共轭转置。

根据机载雷达第k个多普勒通道空时处理的时域导向矢量S_kt和空时处理过程后第k个多普勒通道、第l个距离门的对应输出矢量y_kl，计算得到机载雷达第k个多普勒通道、第l个距离门对应的滤波输出z_kl，其表达式为：

z_kl＝y_kl×conj(S_kt[1:V])

其中，conj[·]表示取共轭操作，S_kt[1:V]表示取机载雷达第k个多普勒通道空时处理的时域导向矢量S_kt中第1到第V个元素构成的矢量；V表示第l个距离门的扩展雷达回波数据Q_l包含的样本个数。

所述机载雷达第k个多普勒通道、第l个距离门单元对应的滤波输出z_kl完成了第k个多普勒通道、第l个距离门的杂波抑制。

步骤6，令k加1，依次重复步骤3至步骤5，直到得到机载雷达第K个多普勒通道、第l个距离门对应的滤波输出z_Kl，并将此时得到的机载雷达第1个多普勒通道、第l个距离门对应的滤波输出z_1l至机载雷达第K个多普勒通道、第l个距离门对应的滤波输出z_Kl，作为机载雷达第l个距离门对应的滤波输出z_l，z_l＝[z_1l z_2l … z_kl … z_Kl]^H，上标H表示共轭转置。

步骤7，令l加1,依次重复步骤3至步骤6，直到得到机载雷达第L个距离门对应的滤波输出z_L，并将此时得到的机载雷达第l个距离门对应的滤波输出z_l至机载雷达第L个距离门对应的滤波输出z_L，作为杂波抑制后机载雷达的距离—多普勒雷达回波数据Z，Z＝[z₁ z₂… z_l … z_L]。

至此，本发明的一种非均匀杂波环境下的雷达杂波抑制方法结束。

通过以下仿真实验对本发明效果作进一步验证说明。

(一)实验参数及实验条件

机载雷达的天线阵面结构采用椭圆斜侧阵面，并且各个阵元均匀排列在椭圆斜侧阵面上，机载雷达的天线阵面的方位向和俯仰向分别均匀排列56个阵元和10个阵元，阵元间距为d＝λ/2，λ为机载雷达的载波波长，机载雷达在每个相干处理周期内发射的脉冲个数M为220，机载雷达的原始雷达回波数据X包含的距离门个数为349。

(二)实验内容及结果分析

分别采用扩展因子化方法(mDT)和本发明方法对原始雷达回波数据进行杂波抑制，得到的两幅距离—多普勒图分别如图2(a)、图2(b)所示；图2(a)是使用扩展因子化方法(mDT)进行杂波抑制后得到的距离—多普勒图，其中横坐标为多普勒通道，纵坐标为距离门，m＝3；图2(b)是使用本发明方法进行杂波抑制后得到的距离—多普勒图，其中横坐标为多普勒通道，纵坐标为距离门；图2(a)、图2(b)中的白色区域分别表示杂波抑制后的杂波剩余功率分布，黑色区域分别表示杂波抑制后的噪声剩余功率分布。

从图2(a)中可以看出，采用扩展因子化方法(mDT)对原始雷达回波数据进行杂波抑制后，杂波剩余功率较强，尤其在旁瓣杂波区，m＝3；从图2(b)中可以看出，本发明方法能够有效降低杂波剩余，对旁瓣杂波的抑制性能较好，有利于后续对目标进行检测。

分别使用扩展因子化方法(mDT)和本发明方法对对原始雷达回波数据进行杂波抑制，分别得到两种方法的杂波剩余比较图，如图3(a)和图3(b)所示；图3(a)是分别使用扩展因子化方法(mDT)和本发明方法进行杂波抑制处理后1-150号距离门对应的杂波剩余图，其中横坐标表示多普勒通道，纵坐标表示杂波剩余，单位为dB，m＝3；图3(b)是分别使用扩展因子化方法(mDT)和本发明方法进行杂波抑制处理后150-349号距离门对应的杂波剩余图；其中横坐标表示多普勒通道，纵坐标表示杂波剩余，单位为dB，m＝3。

从图3(a)和图3(b)可以看出，与扩展因子化方法(mDT)相比，本发明方法在保证信号增益不变的前提下有效抑制旁瓣杂波，m＝3。

综上所述，仿真实验验证了本发明的正确性，有效性和可靠性。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围；这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (5)

1.一种非均匀杂波环境下的雷达杂波抑制方法，其特征在于，包括以下步骤：

步骤1，确定机载雷达，并获取机载雷达的原始雷达回波数据X，然后确定机载雷达的时域滑窗孔径K_t；

步骤2，初始化：令k∈{1,2,…,K}，k表示第k个多普勒通道，K表示机载雷达的原始雷达回波数据X在频率域包含的多普勒总个数，k的初始值为1；令l∈{1,2,…,L}，l表示第l个距离门，L表示机载雷达的原始雷达回波数据X包含的距离门总个数，l的初始值为1；

步骤3，依次计算机载雷达第k个多普勒通道进行空时处理的时域导向矢量S_kt，以及机载雷达第k个多普勒通道的最终空时导向矢量S_k；

步骤4，将第l个距离门的原始雷达回波数据记为x_l，并计算第l个距离门的扩展雷达回波数据Q_l，然后依次计算第l个距离门对应的训练样本数据Train_l、时域滑窗后第l个距离门对应的扩展训练样本数据Train_l′，以及第l个距离门对应的协方差矩阵R_l；

步骤5，计算机载雷达第k个多普勒通道、第l个距离门对应的自适应权值w_kl，进而依次计算空时处理过程后第k个多普勒通道、第l个距离门的对应输出矢量y_kl和机载雷达第k个多普勒通道、第l个距离门对应的滤波输出z_kl；

步骤6，令k加1，依次重复步骤3至步骤5，直到得到机载雷达第K个多普勒通道、第l个距离门对应的滤波输出z_Kl，并将此时得到的机载雷达第1个多普勒通道、第l个距离门对应的滤波输出z_1l至机载雷达第K个多普勒通道、第l个距离门对应的滤波输出z_Kl，作为机载雷达第l个距离门对应的滤波输出z_l，z_l＝[z_1l z_2l … z_kl … z_Kl]^H，上标H表示共轭转置；

步骤7，令l加1,依次重复步骤3至步骤6，直到得到机载雷达第L个距离门对应的滤波输出z_L，并将此时得到的机载雷达第l个距离门对应的滤波输出z_l至机载雷达第L个距离门对应的滤波输出z_L，作为杂波抑制后机载雷达的距离—多普勒雷达回波数据Z，Z＝[z₁ z₂ …z_l … z_L]。

2.如权利要求1所述的一种非均匀杂波环境下的雷达杂波抑制方法，其特征在于，在步骤1中，所述原始雷达回波数据X为N×M×L维矩阵，N表示机载雷达的天线阵面方位向均匀包含的阵元个数，L表示机载雷达的原始雷达回波数据X包含的距离门个数，M表示机载雷达在每个相干处理周期内发射的脉冲个数；

所述机载雷达的时域滑窗孔径K_t，其满足条件为：

2×(M-K_t+1)≥2×K_t×N

或者K_t为整数。

3.如权利要求1所述的一种非均匀杂波环境下的雷达杂波抑制方法，其特征在于，在步骤3中，所述机载雷达第k个多普勒通道进行空时处理的时域导向矢量S_kt，以及所述机载雷达第k个多普勒通道的最终空时导向矢量S_k，其表达式分别为：

$S_{kt}={\left[\begin{array}{cccccc}1& e^{j2{\mathrm{\πf}}_{dk}}& ...& e^{j2{\mathrm{\πf}}_{dk}m}& ...& e^{j2{\mathrm{\πf}}_{dk}(M-1)}\end{array}\right]}^H$

$S_k={S_{kt}}^{\′}\⊗S_s$

其中，S_s表示机载雷达空时处理过程中的空域导向矢量，表示Kronecker乘积，n∈{0,1,2,…,N-1}，表示第n个阵元接收机载雷达的脉冲回波数据时相对第1个阵元的相位偏移，f_s表示机载雷达的空域频率，N表示机载雷达天线阵面方位向均匀包含的阵元个数，上标H表示共轭转置，

m∈{0,1,2,…,M-1}，M表示机载雷达在每个相干处理周期内发射的脉冲个数，表示机载雷达第m个脉冲相对第1个脉冲的相位偏移，f_dk表示机载雷达第k个多普勒通道归一化的时域多普勒频率，且f_dk＝k-1/K；机载雷达的滑窗时域孔径K_t远远小于M，则将机载雷达第k个多普勒通道进行空时处理的时域导向矢量S_kt中第1到第K_t个元素，作为机载雷达第k个多普勒通道进行空时处理的优化时域导向矢量S′_kt，S′_kt＝S_kt[1:K_t]。

4.如权利要求1所述的一种非均匀杂波环境下的雷达杂波抑制方法，其特征在于，步骤4的具体过程为：

首先，将第l个距离门的原始雷达回波数据x_l，其维度为N×M，表达式为：

x_l＝[x′_1l x′_2l … x′_ml … x′_Ml]

其中，x′_ml表示第m个脉冲、第l个距离门的原始雷达回波数据，且是N×1维列矢量；在满足且K_t为整数的条件下任意取一个整数r，作为机载雷达的时域滑窗孔径值，进而对第l个距离门的原始雷达回波数据x_l进行孔径值为r的时域滑窗，得到第l个距离门的扩展雷达回波数据Q_l，维度为N_r×V，其表达式为：

Q_l＝[q_1l q_2l … q_vl … q_Vl]

其中，N_r＝N×r，V表示第l个距离门的扩展雷达回波数据Q_l包含的样本个数，V＝M-r+1；q_vl表示第v个样本、第l个距离门的扩展雷达回波数据，

q_vl＝[x′_vl x′_(v+1)l … x′_(v+r-1)l]^H，x′_vl表示第v个样本、第l个距离门的原始雷达回波数据，上标H表示共轭转置，v＝1,2,…,V；所述第v个样本与对机载雷达的原始雷达回波数据X进行滑窗的起始脉冲编号对应；

根据第l个距离门的原始雷达回波数据x_l，分别将第l-1个距离门的原始雷达回波数据x_l-1、将第l+1个距离门的原始雷达回波数据x_l+1，进而得到第l个距离门对应的训练样本数据Train_l，Train_l＝[x_l-1 x_l+1]，其维度为N×M×2；

然后，对第l个距离门对应的训练样本数据Train_l中包含的所有元素分别进行孔径值为r的时域滑窗，得到时域滑窗后第l个距离门对应的扩展训练样本数据Train_l′，维度为N_r×2V，N_r＝N×r；Train_l′＝[Q_l-1 Q_l+1]，Q_l-1表示第l-1个距离门的扩展雷达回波数据，Q_l+1表示第l+1个距离门的扩展雷达回波数据；

最后，根据时域滑窗后第l个距离门对应的扩展训练样本数据Train_l′，计算第l个距离门对应的协方差矩阵R_l，其表达式为：

5.如权利要求1所述的一种非均匀杂波环境下的雷达杂波抑制方法，其特征在于，在步骤5中，所述机载雷达第k个多普勒通道、第l个距离门对应的自适应权值w_kl，其计算式为：

$w_{kl}=\frac{{R_l}^{-1}{S}_k}{S_k^H{R_l}^{-1}{S}_k}$

其中，R_l表示第l个距离门对应的协方差矩阵，S_k表示机载雷达第k个多普勒通道的最终空时导向矢量，上标H表示共轭转置，上标-1表示求逆操作；

根据机载雷达第k个多普勒通道空时处理的时域导向矢量S_kt和空时处理过程后第k个多普勒通道、第l个距离门的对应输出矢量y_kl，计算得到机载雷达第k个多普勒通道、第l个距离门对应的滤波输出z_kl，其表达式为：

z_kl＝y_kl×conj(S_kt[1:V])

其中，conj[·]表示取共轭操作，S_kt[1:V]表示取机载雷达第k个多普勒通道空时处理的时域导向矢量S_kt中第1到第V个元素构成的矢量；V表示第l个距离门的扩展雷达回波数据Q_l包含的样本个数；y_kl表示空时处理过程后第k个多普勒通道、第l个距离门的对应输出矢量，Q_l表示第l个距离门的扩展雷达回波数据，同时将第l个距离门的扩展雷达回波数据作为空时处理过程中第l个距离门的输入数据。