一种用于指控系统的协方差矩阵的分解方法技术方案

本发明专利技术公开一种用于指控系统的协方差矩阵的分解方法，包括：S10、获取采集到的目标航迹信息；S20、将目标航迹信息进行融合，得到目标航迹信息的协方差矩阵matrixP；S30、对协方差矩阵matrixP进行预处理后，判断协方差矩阵matrixP中是否存在异常行；S40、对存在异常行的协方差矩阵matrixP进行修正并将修正后的协方差矩阵matrixP进行矩阵分解。该分解方法在矩阵的非正定条件下，仍然能够进行分解，满足UKF滤波器计算Sigma点集的需求，达到提升指控系统航迹滤波器鲁棒性的目的。系统航迹滤波器鲁棒性的目的。系统航迹滤波器鲁棒性的目的。

【技术实现步骤摘要】
一种用于指控系统的协方差矩阵的分解方法


[0001]本专利技术涉及数据处理
更具体地，涉及一种用于指控系统的协方差矩阵的分解方法。

技术介绍

[0002]在指控系统中，需要对于多个信息源观测生成的航迹信息进行准确的关联，融合生成综合航迹，以产生统一的空情情报态势，其中，Kalman滤波方法是一种常用的融合方法。KF(卡尔曼滤波)和EKF(扩展卡尔曼滤波)、UKF(无迹卡尔曼滤波)等线性和非线性滤波方法，作为实时状态估算的有力工具，已经广泛应用于工程实践的各个领域。
[0003]UKF的核心思想，是通过采样策略逼近非线性函数分布，利用计算的Sigma点集逼近真实概率密度函数，并直接利用非线性函数传递Sigma点集，从而避免了因线性化非线性函数引入的误差。在计算Sigma采样点集的过程中，需要对协方差矩阵进行分解，常用的分解方法为Cholesky(楚列斯基分解)。
[0004]协方差矩阵能够进行Cholesky分解的前提是矩阵正定，但在信息源的实际测量过程中，由于测量环境、目标状态、设备能力等多方面条件的制约，可能出现滤波过本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种用于指控系统的协方差矩阵的分解方法，其特征在于，包括：S10、获取采集到的目标航迹信息；S20、将所述目标航迹信息进行融合，得到所述目标航迹信息的协方差矩阵matrixP；S30、对所述协方差矩阵matrixP进行预处理后，判断所述协方差矩阵matrixP中是否存在异常行；S40、对存在异常行的协方差矩阵matrixP进行修正并将修正后的协方差矩阵matrixP进行矩阵分解。2.根据权利要求1所述的分解方法，其特征在于，对所述协方差矩阵matrixP进行预处理包括：遍历所述协方差矩阵matrixP中所有元素p
i,j
如果，p
i,j
＜0，则令p
i,j
＝|p
i,j
|；如果，则令p
i,j
＝Δicpt息；其中，i代表所述协方差矩阵matrixP的行数，j代表所述协方差矩阵matrixP的行数的列数。3.根据权利要求1所述的分解方法，其特征在于，对于不存在异常行的协方差矩阵matrixP采用如下公式进行分解：对于首行元素：对于第二行及第二行之后的元素：得到分解后的矩阵其中，i代表所述协方差矩阵matrixP的行数，j代表所述协方差矩阵matrixP的行数的列数。4.根据权利要求1所述的分解方法，其特征在于，所述步骤S20包括：采用凸组合的方法将目标航迹信息进行融合。5.根据权利要求3所述的分解方法，其特征在于，在所述步骤S30之前，还包括：设置计算截断误差参数Δicpt＝10
‑
10
，矩阵最大修正次数ncorrect
max
＝n
...

【专利技术属性】
技术研发人员：刘鑫，张耀东，王中文，刘明阳，
申请(专利权)人：北京电子工程总体研究所，
类型：发明
国别省市：