CN108387884A - 基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法 - Google Patents

Abstract

本发明公开了一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，其主要思路为：确定机载雷达，并且分别获取N_max个距离门雷达回波数据和N_max个待检测单元数据；然后确定杂波脊；令l∈{1,2,…,N_max}，计算第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵利用第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵计算第l个待检测单元数据x_l的重构空时二维协方差矩阵令l的值加1，直到得到第N_max个待检测单元数据x_l的重构空时二维协方差矩阵并利用第N_max个待检测单元数据的重构空时二维协方差矩阵计算用于处理第N_max个待检测单元数据的权进而得到空时自适应处理结果，所述空时自适应处理结果为基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制结果。

Description

基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法

技术领域

本发明属于雷达技术领域，特别涉及一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，适用于非均匀环境下机载预警雷达抑制强地杂波以及检测地面慢速运动目标。

背景技术

机载预警雷达在接收目标回波信号的同时，不可避免地会受到复杂强地杂波的影响。由于复杂强地杂波严重影响了机载预警雷达对地面运动目标的检测，如何减轻或消除地杂波对检测性能的影响，是机载预警雷达所面临的难题；空时自适应处理技术能够从空域和时域两维上同时区分运动目标和杂波，能够有效改善机载预警雷达的杂波抑制和运动目标检测性能，尤其是对慢速目标的检测更为有效，因而得到了广泛的关注。

由于待检测单元杂波信号的统计特性通常事先是不知道的，传统空时自适应处理STAP需要选取待检测单元附近的回波数据作为训练样本对其进行估计；为获得较优的STAP性能(相对最优处理的性能损失不超过3dB),所需要的训练样本数不应少于系统自由度的两倍，此处假设所选用的训练样本与待检测单元数据满足独立同分布条件，即这些训练样本是均匀的且和待检测单元数据的统计特性相同；然而，由于实际场景中存在复杂的地表类型、目标污染等情况，机载雷达通常都工作在非均匀杂波环境中；另外，当雷达天线旋转时，这也会导致机载雷达回波数据非均匀，使得上述独立同分布假设在实际环境中很难满足，进而影响待检测单元杂波协方差矩阵估计精度，导致STAP的杂波抑制性能显著下降。

为解决以上技术问题，可以采用基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制的技术。基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法实际上就是利用待检测单元的数据重构待检测单元的杂波协方差矩阵，该方法可以有效克服训练样本的非均匀性，并且不需要人为设置参数，在工程上比较容易实现。

发明内容

针对上述现有技术存在的问题，本发明的目的在于提出一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，该种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法能够有效抑制非均匀环境下机载预警雷达的强地杂波，提高待检测单元杂波协方差矩阵的估计精度，工程上易于实现，同时又能取得良好的杂波抑制和目标检测性能；本发明方法能有效抑制地杂波，并能明显改善对微弱慢速运动的目标检测性能。

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

一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，包括以下步骤：

步骤1，确定机载雷达，并且分别获取N_max个距离门雷达回波数据和N_max个待检测单元数据；然后确定杂波脊；N_max为大于1的正整数；

步骤2，分别确定第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵初始值和第l个待检测单元数据x_l的噪声功率初始值

初始化：令i表示第i次迭代，i的初始值为1，l∈{1,2,...,N_max}，l的初始值为1；N_max表示雷达回波数据包含的距离门雷达回波数据总个数，与待检测单元数据总个数取值相等；

步骤3，计算第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵和第i次迭代后第l个待检测单元数据x_l的噪声功率

步骤4，若i<γ，则令i的值加1，返回步骤3；其中，γ为预先设定的最大迭代次数；

若i≥γ，则迭代停止，并将迭代停止时对应的第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵记为第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵

步骤5，利用第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵计算第l个待检测单元数据x_l的重构空时二维协方差矩阵

步骤6，将i的值初始化为1；然后令l的值加1，返回步骤3，直到得到第N_max个待检测单元数据x_l的重构空时二维协方差矩阵并利用第N_max个待检测单元数据的重构空时二维协方差矩阵计算用于处理第N_max个待检测单元数据的权进而得到空时自适应处理结果，所述空时自适应处理结果为基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制结果。

本发明的有益效果：

第一，本发明方法利用基于知识辅助稀疏渐进最小方差方法的优势，充分考虑了训练样本与待检测单元数据之间的差异，通过将待检测单元数据分解成多个空域快拍数据，再利用空域快拍数据重构待检测单元的空域协方差矩阵，最后利用杂波脊先验信息得到待检测单元的空时协方差矩阵，使得本发明方法能够有效抑制复杂强地杂波，从而显著改善对运动目标的检测性能。

第二，本发明方法通过利用知识辅助稀疏渐进最小方差谱重构技术重构待检测单元的杂波协方差矩阵，该方法不需要人为设置参数，运算量小，相比传统的杂波谱重构方法谱估计精度更高，因此该方法的实时性更好，实用价值更高。

附图说明

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

图1为本发明的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法流程图；

图2为机载雷达几何构型图；

图3(a)为正侧阵时最优处理空时二维谱示意图；

图3(b)为正侧阵时本发明方法空时二维谱示意图

图4(a)为75度偏航角时最优处理空时二维谱示意图；

图4(b)为75度偏航角时本发明方法空时二维谱示意图

图5(a)为正侧阵时本发明方法与最优处理方法改善因子对比示意图；

图5(b)为75度偏航角时本发明方法与最优处理方法改善因子对比示意图。

具体实施方式

参照图1，为本发明的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法流程图；所述基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，包括以下步骤：

步骤1，确定机载雷达，机载雷达的检测范围内存在杂波散射体S和目标，所述杂波散射体S为机载雷达波束照射到地面上时能够反射机载雷达不需要的回波的物体。

机载雷达发射信号并接收雷达回波数据，所述雷达回波数据包含N_max个距离门雷达回波数据，依次记为第1个距离门的雷达回波数据、第2个距离门的雷达回波数据、...、第N_max个距离门的雷达回波数据，每个距离门的雷达回波数据是机载雷达发射M个脉冲后接收到的雷达回波数据。

参照图2，为机载雷达几何构型图；以机载雷达载机平台在地面上的投影点为原点O建立OXYZ坐标系，其中X轴为机载雷达载机平台飞行方向，Y轴为地面上与机载雷达载机平台飞行方向垂直的方向，将XOY平面作为地平面，Z轴为垂直于地平面的方向，所述地平面未考虑地表的曲率问题。

机载雷达载机平台以高度H、速度v沿着X轴方向飞行；杂波散射体S为相对于机载雷达载机平台斜距为R_l的杂波散射体；将机载雷达载机平台飞行速度v的方向与机载雷达天线阵列轴向的夹角，记为偏航角α；将杂波散射体S和机载雷达天线阵列中心之间的连线与地平面水平方向的夹角，记为俯仰角将杂波散射体S和机载雷达天线阵列中心之间的连线在地平面的投影与机载雷达天线阵列轴向的夹角，记为方位角θ。

其中，将机载雷达载机平台到杂波散射体S的直线距离记为斜距，第l个距离门的雷达回波数据对应斜距为R_l，R_l＝H+(l-1)△R，△R表示机载雷达的距离分辨率，B表示机载雷达的接收机带宽，c表示光传播速度，l＝1,2,…,N_max。

当l＝1时，R_l＝H，l＝N_max，对应第N_max个距离门的雷达回波数据；当l＝N_max时，R_l＝H+(N_max-1)△R，对应第N_max个距离门的雷达回波数据，且取值与机载雷达的最大不模糊距离R_u取值相等，T_r表示机载雷达的脉冲重复周期。

将每个距离门的雷达回波数据分别记为一个待检测单元数据，进而得到N_max个待检测单元数据。

杂波散射体S的空域频率表示为

杂波散射体S的归一化空域频率表示为f_s：

其中，d′表示机载雷达天线阵列的阵元间距，ψ为机载雷达载机平台飞行速度方向与杂波散射体S之间的夹角，λ为机载雷达发射信号的波长，表示杂波散射体S的空域频率，f_s表示杂波散射体S的归一化空域频率，cos表示余弦函数；将机载雷达载机平台飞行速度v的方向与机载雷达天线阵列轴向的夹角，记为偏航角α；将杂波散射体S和机载雷达天线阵列中心之间的连线与地平面的夹角，记为俯仰角将杂波散射体S和机载雷达天线阵列中心之间的连线在地平面的投影与机载雷达天线阵列轴向的夹角，记为方位角θ。

杂波散射体S的多普勒频率表示为

其中，v表示机载雷达载机平台飞行速度。

杂波散射体S的归一化多普勒频率表示为f_d：

其中，v为机载雷达载机平台飞行速度，ψ为机载雷达载机平台飞行速度方向与杂波散射体S之间的夹角，λ为机载雷达发射信号的波长，为杂波散射体S的多普勒频率，f_d为杂波散射体S的归一化多普勒频率。

确定杂波脊，所述杂波脊是由机载雷达载机平台运动引起的反映多普勒频率和空域频率之间耦合关系的轨迹，其轨迹确定方程为：

其中，f_r表示机载雷达的脉冲重复频率，f_s表示杂波散射体S的归一化空域频率，d′表示机载雷达天线阵列的阵元间距，表示俯仰角，cos表示余弦函数，sin表示正弦函数，α表示偏航角。

所述杂波脊由机载雷达系统参数和杂波块相对于机载雷达载机平台的角度唯一确定，其与机载雷达天线阵列布置形式、地表覆盖内型、目标污染等因素无关，因此可以将杂波脊当作一种先验信息来运用。

将所述杂波脊划分为Nc个杂波块，每个杂波块分别记为一个子杂波散射体，进而得到Nc个子杂波散射体。

定义杂波脊上第m个子杂波散射体的归一化多普勒频率为f_d,m，然后得到杂波脊上第m个子杂波散射体的归一化多普勒频率f_d,m的时域导向矢量为b(f_d,m)，其表达式为：

其中，m∈{1,2,…,Nc}，f_d,m表示杂波脊上第m个子杂波散射体的归一化多普勒频率，M表示一个相干处理间隔内机载雷达发射的脉冲个数，上标T表示矩阵转置运算，e表示自然对数的底数；j表示虚数单位，

定义杂波脊上第m个子杂波散射体的归一化空域频率为f_s,m，然后得到杂波脊上第m个子杂波散射体归一化空域频率f_s,m的空域导向矢量为a(f_s,m)，其表达式为：

其中，N表示机载雷达天线阵列包含的阵元总个数，上标T表示矩阵转置运算，m∈{1,2,…,Nc}，f_s,m表示杂波脊上第m个子杂波散射体的归一化空域频率。

定义杂波脊上第m个子杂波散射体的归一化多普勒频率为f_d,m，则杂波脊上第m个子杂波散射体的空时二维导向矢量为v(f_d,m,f_s,m)，其表达式为：

其中，f_d,m表示杂波脊上第m个子杂波散射体的归一化多普勒频率，f_s,m表示杂波脊上第m个子杂波散射体的归一化空域频率，m∈{1,2,…,Nc}，b(f_d,m)表示杂波脊上第m个子杂波散射体的归一化多普勒频率f_d,m的时域导向矢量，a(f_s,m)表示杂波脊上第m个子杂波散射体归一化空域频率f_s,m的空域导向矢量，表示Kronecker积，Nc表示对杂波脊进行划分后包括的子杂波散射体总个数；进而得到N×Nc维空域导向矢量矩阵A和MN×Nc维空时导向矢量矩阵V。

所述N×Nc维空域导向矢量矩阵A，其表达式为：

所述MN×Nc维空时导向矢量矩阵V，其表达式为：

其中，a(f_s,m)表示杂波脊上第m个子杂波散射体归一化空域频率f_s,m的空域导向矢量，v(f_d,m,f_s,m)表示杂波脊上第m个子杂波散射体的空时二维导向矢量，m∈{1,2,…,Nc}，Nc表示对杂波脊进行划分后包括的子杂波散射体总个数。

步骤2，分别确定第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵初始值和噪声功率初始值其中第l个待检测单元数据x_l(MN×1)可以表示如下：

其中，y_lk表示第l个待检测单元数据x_l中第k个脉冲的雷达回波数据，第l个待检测单元数据x_l中每个脉冲的雷达回波数据分别为一个空域采样数据；上标T表示矩阵转置运算，l＝1,2,…,N_max，k＝1,2,…,M，M表示一个相干处理间隔内机载雷达发射的脉冲个数。

根据第l个待检测单元数据x_l的第k个脉冲空域采样数据可以计算第l个待检测单元数据x_l的空域协方差矩阵其表达式为：

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

确定第l个待检测单元数据x_l在杂波脊上第m个子杂波散射体的杂波功率初始值为其表达式为：

其中，上标H表示共轭转置操作，|·|²表示取绝对值后再平方操作。

初始化：令i表示第i次迭代，i的初始值为1，l∈{1,2,…,N_max}，l的初始值为1；N_max表示雷达回波数据包含的距离门雷达回波数据总个数，与待检测单元数据总个数取值相等。

进而得到第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵初始值为 其表达式为：

其中，R^Nc×Nc表示Nc×Nc维实数矩阵，m∈{1,2,…,Nc}，Nc表示对杂波脊进行划分后包括的子杂波散射体总个数，表示第l个待检测单元数据x_l在杂波脊上第1个子杂波散射体的杂波功率初始值，表示第l个待检测单元数据x_l在杂波脊上第2个子杂波散射体的杂波功率初始值，表示第l个待检测单元数据x_l在杂波脊上第Nc-1个子杂波散射体的杂波功率初始值，表示第l个待检测单元数据x_l在杂波脊上第Nc个子杂波散射体的杂波功率初始值，Nc表示对杂波脊进行划分后包括的子杂波散射体总个数。

确定第l个待检测单元数据x_l的空域重构协方差矩阵初始值

A表示N×Nc维空域导向矢量矩阵，上标H表示共轭转置操作，I_N表示N×N维单位矩阵，N表示机载雷达天线阵列包含的阵元总个数，表示第l个待检测单元数据x_l的噪声功率初始值。

所述第l个待检测单元数据x_l的噪声功率初始值为其表达式为：

其中，下标n指代噪声，||·||²表示矩阵的2范数的平方。

步骤3，计算第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵

计算第i次迭代后第l个待检测单元数据x_l在杂波脊上第m个子杂波散射体的杂波功率其表达式为：

其中，m＝1,…,Nc，表示第i-1次迭代后第l个待检测单元数据x_l在杂波脊上第m个子杂波散射体的杂波功率，表示求逆操作，表示第i-1次迭代后第l个待检测单元数据x_l的空域重构协方差矩阵，上标-1表示矩阵求逆操作。

进而计算得到第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵其表达式为：

其中，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第1个子杂波散射体的杂波功率，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第2个子杂波散射体的杂波功率，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第Nc-1个子杂波散射体的杂波功率，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第Nc个子杂波散射体的杂波功率。

然后，根据第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵和第i次迭代后第l个待检测单元数据x_l的噪声功率计算得到第i次迭代后第l个待检测单元数据x_l的空域重构协方差矩阵其表达式为：

其中，A表示MN×Nc维空域导向矢量矩阵，A∈C^N×Nc，C^N×Nc表示N×Nc维复数矩阵；∈表示属于，C^N×N表示N×N维复数矩阵；上标H表示共轭转置操作，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵，I_N表示N×N维单位矩阵，N表示机载雷达天线阵列包含的阵元总个数，表示第i次迭代后第l个待检测单元数据x_l的噪声功率。

所述第i次迭代后第l个待检测单元数据x_l的噪声功率为其表达式为：

其中，a(f_s,m)表示杂波脊上第m个子杂波散射体归一化空域频率f_s,m的空域导向矢量，表示与相乘操作，表示求逆操作，表示第i-1次迭代后第l个待检测单元数据x_l的空域重构协方差矩阵，上标H表示共轭转置操作，tr表示求迹运算，表示第l个待检测单元数据x_l的空域协方差矩阵。

步骤4，若i<γ，则令i的值加1，返回步骤3；其中，γ为预先设定的最大迭代次数，通常取10。

若i≥γ，则迭代停止，并将迭代停止时对应的第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵记为第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵所述第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵是基于知识辅助稀疏渐进最小方差方法计算得到的杂波功率矩阵。

步骤5，利用杂波脊先验信息和第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵计算第l个待检测单元数据x_l的重构空时二维协方差矩阵其表达式为：

其中，上标H表示共轭转置操作，V表示MN×Nc维空时导向矢量矩阵，表示第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵，表示第l个待检测单元数据x_l在杂波脊上第1个子杂波散射体的最终杂波功率，表示第l个待检测单元数据x_l在杂波脊上第2个子杂波散射体的最终杂波功率，表示第l个待检测单元数据x_l在杂波脊上第Nc-1个子杂波散射体的最终杂波功率，表示第l个待检测单元数据x_l在杂波脊上第Nc个子杂波散射体的最终杂波功率，I_MN表示MN×MN维单位矩阵；β表示事先设定的对角加载系数，1<β<10，本实施例中β取值为2。

步骤6，将i的值初始化为1；然后令l的值加1，返回步骤3，直到得到第N_max个待检测单元数据x_l的重构空时二维协方差矩阵并利用第N_max个待检测单元数据的重构空时二维协方差矩阵计算用于处理第N_max个待检测单元数据的权其表达式为：

其中，μ表示系数标量，上标H表示共轭转置操作，()^-1表示矩阵求逆操作，v_t(f_d,t,f_s,t)表示目标的空时二维导向矢量，f_d,t表示目标的归一化多普勒频率，f_s,t表示目标的归一化空域频率，下标t指代目标。

所述f_s,t表示目标的归一化空域频率，其表达式为：

其中，d′表示机载雷达天线阵列的阵元间距，ψ₀为机载雷达载机平台飞行速度方向与目标之间的夹角，λ为机载雷达发射信号的波长，cos表示余弦函数。

所述f_d,t表示目标的归一化多普勒频率，其表达式为：

其中，v表示机载雷达载机平台飞行速度，f_r表示机载雷达的脉冲重复频率，cos表示余弦函数，λ为机载雷达发射信号的波长，ψ₀为机载雷达载机平台飞行速度方向与目标之间的夹角。

所述v_t(f_d,t,f_s,t)表示目标的空时二维导向矢量，其表达式为：

其中，表示Kronecker积；a(f_s,t)表示目标的归一化空域频率f_s,t的空域导向矢量，b(f_d,t)表示目标的归一化多普勒频率f_d,t的时域导向矢量，其表达式分别为：

其中，N表示机载雷达天线阵列包含的阵元总个数，上标T表示矩阵转置运算，M表示一个相干处理间隔内机载雷达发射的脉冲个数，e表示自然对数的底数；j表示虚数单位，

然后对所述用于处理第N_max个待检测单元数据的权与第N_max个待检测单元数据进行空时自适应处理，即将所述用于处理第N_max个待检测单元数据的权与第N_max个待检测单元数据进行内积，进而得到空时自适应处理结果，所述空时自适应处理结果为基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制结果；经过本发明方法处理后杂波抑制效果明显，对微弱慢速运动的目标检测性能明显提高。

至此，本发明的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制的方法结束。

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

(一)仿真实验数据说明。

为了验证本发明方法的准确性，通过仿真实验予以证明；仿真实验参数如下：

(二)仿真结果及分析

本发明的仿真结果分别如图3(a)、图3(b)、图4(a)、图4(b)、图5(a)和图5(b)所示，图3(a)为正侧阵时最优处理方法空时二维谱示意图，图3(b)为正侧阵时本发明方法空时二维谱示意图，图4(a)为非正侧阵时最优处理方法空时二维谱示意图，图4(b)为非正侧阵时本发明方法空时二维谱示意图，其中横坐标都为归一化多普勒频率，纵坐标都为归一化空域频率；图5(a)为正侧阵时本发明方法与最优处理方法改善因子对比示意图，图5(b)为非正侧阵时本发明方法与最优处理方法改善因子对比示意图；其中横坐标都为归一化多普勒频率，纵坐标都为改善因子(dB)。

从图3(b)和图4(b)可以看出，无论机载雷达天线以正侧阵方式还是非正侧阵方式安置时，运用本发明方法估计出来的待检测单元杂波谱都与最优杂波谱很相似，这说明本发明方法能够有效估计待检测单元的杂波协方差矩阵，从而保证了机载雷达对慢速运动目标的检测性能。

从图5(a)和图5(b)可以看出，无论机载雷达天线以正侧阵方式还是非正侧阵方式安置时，本发明方法的改善因子和最优处理方法的改善因子基本重合，这说明采用本发明方法的机载雷达能获得优异的杂波抑制性能。

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

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

Claims (7)

1.一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，其特征在于，包括以下步骤：

步骤1，确定机载雷达，并且分别获取N_max个距离门雷达回波数据和N_max个待检测单元数据；然后确定杂波脊；N_max为大于1的正整数；

步骤2，分别确定第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵初始值和第l个待检测单元数据x_l的噪声功率初始值

初始化：令i表示第i次迭代，i的初始值为1，l∈{1,2,…,N_max}，l的初始值为1；N_max表示雷达回波数据包含的距离门雷达回波数据总个数，与待检测单元数据总个数取值相等；

步骤3，计算第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵和第i次迭代后第l个待检测单元数据x_l的噪声功率

步骤4，若i<γ，则令i的值加1，返回步骤3；其中，γ为预先设定的最大迭代次数；

若i≥γ，则迭代停止，并将迭代停止时对应的第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵记为第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵

步骤5，利用第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵计算第l个待检测单元数据x_l的重构空时二维协方差矩阵

步骤6，将i的值初始化为1；然后令l的值加1，返回步骤3，直到得到第N_max个待检测单元数据x_l的重构空时二维协方差矩阵并利用第N_max个待检测单元数据的重构空时二维协方差矩阵计算用于处理第N_max个待检测单元数据的权进而得到空时自适应处理结果，所述空时自适应处理结果为基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制结果。

2.如权利要求1所述的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，其特征在于，在步骤1中，所述N_max个距离门雷达回波数据和N_max个待检测单元数据，其获取过程为：

确定机载雷达，机载雷达的检测范围内存在杂波散射体S和目标；机载雷达发射信号并接收雷达回波数据，所述雷达回波数据包含N_max个距离门雷达回波数据，依次记为第1个距离门的雷达回波数据、第2个距离门的雷达回波数据、…、第N_max个距离门的雷达回波数据，每个距离门的雷达回波数据是机载雷达发射M个脉冲后接收到的雷达回波数据；

将每个距离门的雷达回波数据分别记为一个待检测单元数据，进而得到N_max个待检测单元数据；

所述杂波脊是由机载雷达载机平台运动引起的反映多普勒频率和空域频率之间耦合关系的轨迹，其轨迹确定方程为：

其中，将机载雷达载机平台飞行速度v的方向与机载雷达天线阵列轴向的夹角，记为偏航角α；将杂波散射体S和机载雷达天线阵列中心之间的连线与地平面的夹角，记为俯仰角λ为机载雷达发射信号的波长，f_r表示机载雷达的脉冲重复频率，f_d表示杂波散射体S的归一化多普勒频率，v为机载雷达载机平台飞行速度，f_s表示杂波散射体S的归一化空域频率，d′表示机载雷达天线阵列的阵元间距，表示俯仰角，cos表示余弦函数，sin表示正弦函数。

3.如权利要求2所述的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，其特征在于，在步骤2中，所述第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵初始值和第l个待检测单元数据x_l的噪声功率初始值其表达式分别为：

其中，Nc表示对杂波脊进行划分后包括的子杂波散射体总个数，表示第l个待检测单元数据x_l在杂波脊上第1个子杂波散射体的杂波功率初始值，表示第l个待检测单元数据x_l在杂波脊上第2个子杂波散射体的杂波功率初始值，表示第l个待检测单元数据x_l在杂波脊上第Nc-1个子杂波散射体的杂波功率初始值，表示第l个待检测单元数据x_l在杂波脊上第Nc个子杂波散射体的杂波功率初始值，Nc表示对杂波脊进行划分后包括的子杂波散射体总个数；

令表示第l个待检测单元数据x_l在杂波脊上第m个子杂波散射体的杂波功率初始值，其表达式为：

其中，m＝1,…,Nc，a(f_s,m)表示杂波脊上第m个子杂波散射体归一化空域频率f_s,m的空域导向矢量，上标H表示共轭转置操作，表示第l个待检测单元数据x_l的空域协方差矩阵，N表示机载雷达天线阵列包含的阵元总个数，上标T表示矩阵转置运算，M表示一个相干处理间隔内机载雷达发射的脉冲个数，k＝1,2,…,M，y_lk表示第l个待检测单元数据x_l中第k个脉冲的雷达回波数据，f_s,m表示杂波脊上第m个子杂波散射体归一化空域频率，下标n指代噪声，||·||²表示矩阵的2范数的平方，|·|²表示取绝对值后再平方操作。

4.如权利要求3所述的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，其特征在于，步骤2还包括：

计算第l个待检测单元数据x_l的空域重构协方差矩阵初始值

上标H表示共轭转置操作，I_N表示N×N维单位矩阵，N表示机载雷达天线阵列包含的阵元总个数，表示第l个待检测单元数据x_l的噪声功率初始值；A表示N×Nc维空域导向矢量矩阵，其表达式为：

其中，a(f_s,m)表示杂波脊上第m个子杂波散射体归一化空域频率f_s,m的空域导向矢量，N表示机载雷达天线阵列包含的阵元总个数。

5.如权利要求4所述的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，其特征在于，在步骤3中，所述第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵其表达式为：

其中，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第1个子杂波散射体的杂波功率，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第2个子杂波散射体的杂波功率，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第Nc-1个子杂波散射体的杂波功率，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第Nc个子杂波散射体的杂波功率；

令表示第i次迭代后第l个待检测单元数据x_l在杂波脊上第m个子杂波散射体的杂波功率，其表达式为：

其中，m＝1,…,Nc，表示第i-1次迭代后第l个待检测单元数据x_l在杂波脊上第m个子杂波散射体的杂波功率，表示求逆操作，表示第i-1次迭代后第l个待检测单元数据x_l的空域重构协方差矩阵，上标-1表示矩阵求逆操作；

所述第i次迭代后第l个待检测单元数据x_l的噪声功率其表达式为：

其中，a(f_s,m)表示杂波脊上第m个子杂波散射体归一化空域频率f_s,m的空域导向矢量，表示与相乘操作，表示求逆操作，表示第i-1次迭代后第l个待检测单元数据x_l的空域重构协方差矩阵，上标H表示共轭转置操作，tr表示求迹运算，表示第l个待检测单元数据x_l的空域协方差矩阵。

6.如权利要求5所述的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，其特征在于，步骤3还包括：

根据第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵和第i次迭代后第l个待检测单元数据x_l的噪声功率计算得到第i次迭代后第l个待检测单元数据x_l的空域重构协方差矩阵其表达式为：

其中，上标H表示共轭转置操作，A表示MN×Nc维空域导向矢量矩阵，上标H表示共轭转置操作，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵，I_N表示N×N维单位矩阵，N表示机载雷达天线阵列包含的阵元总个数。

7.如权利要求6所述的一种基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法，其特征在于，在步骤5中，所述第l个待检测单元数据x_l的重构空时二维协方差矩阵其表达式为：

其中，上标H表示共轭转置操作，表示第l个待检测单元数据x_l在杂波脊上的最终杂波功率矩阵，表示第l个待检测单元数据x_l在杂波脊上第1个子杂波散射体的最终杂波功率，表示第l个待检测单元数据x_l在杂波脊上第2个子杂波散射体的最终杂波功率，表示第l个待检测单元数据x_l在杂波脊上第Nc-1个子杂波散射体的最终杂波功率，表示第l个待检测单元数据x_l在杂波脊上第Nc个子杂波散射体的最终杂波功率，I_MN表示MN×MN维单位矩阵；β表示事先设定的对角加载系数，1<β<10；V表示MN×Nc维空时导向矢量矩阵，其表达式为：

其中，v(f_d,m,f_s,m)表示杂波脊上第m个子杂波散射体的空时二维导向矢量，其表达式为：

其中，f_d,m表示杂波脊上第m个子杂波散射体的归一化多普勒频率，f_s,m表示杂波脊上第m个子杂波散射体的归一化空域频率，m∈{1,2,…,Nc}，b(f_d,m)表示杂波脊上第m个子杂波散射体的归一化多普勒频率f_d,m的时域导向矢量，M表示一个相干处理间隔内机载雷达发射的脉冲个数，上标T表示矩阵转置运算，a(f_s,m)表示杂波脊上第m个子杂波散射体归一化空域频率f_s,m的空域导向矢量，表示Kronecker积，Nc表示对杂波脊进行划分后包括的子杂波散射体总个数。