[发明专利]基于知识辅助最大似然的机载雷达杂波抑制方法

权利要求书

1.一种基于知识辅助最大似然的机载雷达杂波抑制方法，其特征在于，包括以下步骤：

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

步骤2，分别确定第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵初始值和第l个待检测单元数据x_l的空时重构协方差矩阵初始值

初始化：令i表示内层循环第i次迭代，i的初始值为1；令j′表示外层循环第j′次迭代，j′的初始值为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，|·|表示取绝对值运算，表示取平方根运算，表示第j′-1次迭代后第l个待检测单元数据x_l在杂波脊上第m个子杂波散射体的中间变量矩阵，表示第i-1次迭代后第l个待检测单元数据x_l在杂波脊上的过渡矩阵，表示第i-1次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵，V_l表示第l个待检测单元数据x_l对应的空时导向矢量矩阵，表示第i-1次迭代后第l个待检测单元数据x_l的空域重构协方差矩阵，()^-1表示矩阵求逆运算；

所述步骤3还包括：

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

其中，上标H表示共轭转置操作，V_l表示第l个待检测单元数据x_l对应的空时导向矢量矩阵，表示第i次迭代后第l个待检测单元数据x_l在杂波脊上的杂波功率矩阵，Nc表示对杂波脊进行划分后包括的子杂波散射体总个数；

步骤4，若i＜γ，则令i的值加1，返回步骤3；其中，γ为预先设定的内层循环最大迭代次数；若i≥γ并且j′＜ε，则j′的值加1，并将i的值初始化为1，返回步骤3；其中，ε为预先设定的外层循环最大迭代次数；

若j′≥ε，则迭代停止，并将内层循环迭代停止时对应的第γ次迭代后第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个子杂波散射体的最终杂波功率，V_l表示第l个待检测单元数据x_l对应的空时导向矢量矩阵；

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