(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211140732.1 (22)申请日 2022.09.20 (71)申请人 中国电子科技 集团公司第十四研究 所 地址 210039 江苏省南京市雨 花台区国睿 路8号 (72)发明人 郝艳军　 (74)专利代理 机构 南京知识律师事务所 32 207 专利代理师 康翔　高娇阳 (51)Int.Cl. G01S 13/50(2006.01) G06F 30/20(2020.01) (54)发明名称 一种STAP样本 选择方法 (57)摘要 本发明公开了一种STAP样 本选择方法， 将每 个参考距离单元的多通道信号建模为前向和后 向多通道AR模型， 用Burg算法求解AR模型系数， 根据AR模型系数， 求解每个参考距离单元的多 通 道信号协方差矩阵， 用梯度下降法计算多通道协 方差矩阵的Kahler中值， 从矩阵流形计算每个参 考距离单元的协方差矩阵与其Kahler中值的黎 曼距离， 根据黎曼距离判断样本的奇异性， 剔除 奇异样本， 选择均匀样本， 只需要对每帧数据分 别处理， 计算效率高， 在非均匀杂波的环境下， 当 可用的参考距离单元数量较少， 或非均匀样本的 奇异性不强时， 也可 以有效的挑选出奇异样本， 获得可以用于估计杂波协方差矩 阵的均匀杂波 单元， 计算结果 准确。 权利要求书4页 说明书9页 CN 115453513 A 2022.12.09 CN 115453513 A 1.一种STAP样本 选择方法， 其特 征在于， 包括： 多通道信号AR模型建模： 将每个参考距离单元的多通道信号建模为前向和后向多通道 AR模型， 用Burg算法求 解AR模型系数； 估值协方差矩阵： 根据AR模型系数， 求 解每个参考距离单 元的多通道信号协方差矩阵； 计算Kahler中值： 用梯度下降法计算多通道 协方差矩阵的Kahler中值； 计算检测统计量： 从矩阵流形计算每个参考距离单元的协方差矩阵与其Kahler中值的 黎曼距离； 筛选样本： 根据黎曼距离判断样本的奇异性， 剔除奇异样本， 选择均匀样本 。 2.根据权利要求1所述的STAP样本选择方法， 其特征在于， 所述多通道信号AR模型建 模， 包括： 将每个参考距离单元的二维空时信号视为p个通道的复向量时间序列的采样， 用 表示， 其中 为时间序列， △t为时间间隔， yi(t) 为第i个通道的时间序列采样， y(t)为广义平稳过程， 多通道AR模型的N阶前向滤波器用 表示， 其中eN(t)为p阶向量， 表示前向滤波器的输 出， FkN为 p×p复矩阵， 表示前向滤波器的系数， 为FkN的共轭转置， N阶后向滤波器用 表示， 其中bN(t)表示后向滤波器的输出， BkN表示后向滤波 器的系数。 3.根据权利要求2所述的STAP样本选择方法， 其特征在于， 所述用Burg算法求解AR模型 系数， 包括： 设Ri为yi(t)的协方差矩阵， R‑i为Ri的共轭转置， PN ＝E[eN(t)eN(t)H]为前向滤波器的功率矩阵， P ′N＝E[bN(t)bN(t)H]为后向滤波器的功率矩 阵， 用 表示前向滤波器功率矩阵和协方差权　利　要　求　书 1/4 页 2 CN 115453513 A 2矩阵、 滤波器系数之间的关系， 用 表示后 向滤波器功率矩阵和协方差矩阵、 滤波器系数之间的关系， 用 表示前向滤波器的AR模型， 用 表示后向滤波器的AR模型， 用 表示N阶滤波器的输出， 其中eN(m)、 bN(m)和YN(m) 分 别 表 示 第 m 个 前 向 滤 波 器 输 出 、后 向 滤 波 器 输 出 和 信 号 输 出 ，用 和 表示滤波器输出的递推关系， 其中CN和C′N分别表示前向和后向滤波器的反射系数矩阵， 用 表示极小化滤波器输出的均方误差时的反射系数矩阵， 其中 分别表示N ‑1阶滤波器的输出统计量， M为信号总权　利　要　求　书 2/4 页 3 CN 115453513 A 3