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

权利要求书

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的噪声功率

在所述步骤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的空域协方差矩阵；

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

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

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

在所述步骤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表示对杂波脊进行划分后包括的子杂波散射体总个数；

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