基于ADS-B系统航迹点的修正方法和系统技术方案

一种基于ADS-B系统航迹点的修正方法和系统，其方法包括步骤：对航迹点按接收时间进行排列；根据航迹点的状态信息预测当前时刻航迹点的状态预测值以及计算当前时刻航迹点的状态预测协方差矩阵；根据航迹点的状态信息预测当前时刻航迹点的测量预测值以及计算当前时刻航迹点的残差向量协方差矩阵；识别当前时刻航迹点是否为野值，并根据当前时刻航迹点是否为野值确定当前时刻航迹点的衰减因子；根据当前时刻航迹点的衰减因子计算当前时刻滤波器增益；修正当前时刻航迹点的状态预测值以及当前航迹点的状态预测协方差矩阵。根据本发明专利技术方案，能够有效识别航迹点中的野值，实现对航迹点中的野值进行实时动态的修正，提高航迹点数据的准确度。

【技术实现步骤摘要】
基于ADS-B系统航迹点的修正方法和系统
本专利技术涉及基于ADS-B的民用航空空中监管领域，特别是涉及一种基于ADS-B系统航迹点的修正方法和系统。
技术介绍
ADS-B(AutomaticDependentSurveillanceBroadcast，广播式自动相关监视)监控设备在民用航空空中监管的使用过程中，由于噪声干扰、窜扰、混扰、多径干扰等因素的影响，得到的航迹点数据往往包含了较大的随机误差，造成航迹点数据严重偏离目标真实值，即成为野值。目前，传统的航迹点野值处理方法是通过模糊控制函数判断新的航迹值是否为野值，如果不是野值，则不对航迹估计值进行调整；如果野值较小，则调整增益实现对野值的修正；如果野值较大，则对野值进行剔除。然而，这种野值的处理方法，存在野值判断不灵敏，无法实现对野值实时动态的修正；对于连续成片的野值，还有可能存在野值的判定和处理误差较大的情况，且对于较大的野值，直接将该野值剔除，无法实现航迹点野值的修正，而且降低了航迹点数据的准确度。
技术实现思路
基于此，有必要针对
技术介绍
中，现有的航迹点野值处理方法出现野值判断不灵敏，无法实现对野值实时动态的修正，降低了航迹点数据的完整性的问题，提供一种基于ADS-B系统航迹点的修正方法，能够有效识别航迹点中的野值，实现对航迹点中的野值进行实时动态的修正，提高航迹点数据的准确度。为实现上述目的，本专利技术实施例采用的技术方案如下：一种基于ADS-B系统航迹点的修正方法，包括步骤：对航迹点按接收时间进行排列；根据航迹点的状态信息预测当前时刻航迹点的状态预测值以及计算当前时刻航迹点的状态预测协方差矩阵；根据航迹点的状态信息预测当前时刻航迹点的测量预测值以及计算当前时刻航迹点的残差向量协方差矩阵；识别当前时刻航迹点是否为野值，并根据当前时刻航迹点是否为野值确定当前时刻航迹点的衰减因子；根据当前时刻航迹点的衰减因子计算当前时刻滤波器增益；根据当前时刻航迹点的测量预测值、当前航迹点的残差向量协方差矩阵以及当前时刻滤波器增益修正当前时刻航迹点的状态预测值以及当前航迹点的状态预测协方差矩阵。根据上述的一种基于ADS-B系统航迹点的修正方法，本专利技术提供一种基于ADS-B系统航迹点的修正系统，包括预处理单元、计算修正单元；所述预处理单元对航迹点按接收时间进行排列；所述计算修正单元根据航迹点的状态信息预测当前时刻航迹点的状态预测值以及计算当前时刻航迹点的状态预测协方差矩阵；根据航迹点的状态信息预测当前时刻航迹点的测量预测值以及计算当前时刻航迹点的残差向量协方差矩阵；识别当前时刻航迹点是否为野值，并根据当前时刻航迹点是否为野值确定当前时刻航迹点的衰减因子；根据当前时刻航迹点的衰减因子计算当前时刻滤波器增益；根据当前时刻航迹点的测量预测值、当前航迹点的残差向量协方差矩阵以及当前时刻滤波器增益修正当前时刻航迹点的状态预测值以及当前航迹点的状态预测协方差矩阵。根据本专利技术方案，先对航迹点按接收时间进行排列；然后根据航迹点的状态信息预测当前时刻航迹点的状态预测值以及计算当前时刻航迹点的状态预测协方差矩阵；并预测当前时刻航迹点的测量预测值以及计算当前时刻航迹点的残差向量协方差矩阵；进而识别当前时刻航迹点是否为野值，并根据当前时刻航迹点是否为野值确定当前时刻航迹点的衰减因子；再进一步根据当前时刻航迹点的衰减因子计算当前时刻滤波器增益；最后修正当前时刻航迹点的状态预测值以及当前航迹点的状态预测协方差矩阵。进而达到有效识别航迹点中的野值，实现对航迹点中的野值进行实时动态的修正，提高航迹点数据的准确度。附图说明图1为本专利技术基于ADS-B系统航迹点的修正方法的实施例流程图；图2为本专利技术基于ADS-B系统航迹点的修正方法的一维仿真模型图；图3为本专利技术基于ADS-B系统航迹点的修正方法的一维仿真修正结果图；图4为本专利技术基于ADS-B系统航迹点的修正系统结构图。具体实施方式为使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步的详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本专利技术，并不限定本专利技术的保护范围。请参阅图1，为本专利技术基于ADS-B系统航迹点的修正方法的实施例流程图：步骤S101：对航迹点按接收时间进行排列；具体地，所述按接收时间排列航迹点是指将接收到的航迹点由大地坐标系转化为直角坐标系，再将航迹点的数据按照接收的时间是先后顺序进行排列。步骤S102：根据航迹点的状态信息预测当前时刻航迹点的状态预测值以及计算当前时刻航迹点的状态预测协方差矩阵；所述当前时刻航迹点的状态预测值是指当前时刻飞行目标的状态信息，包括飞行目标的位置、速度、加速度等；所述当前时刻航迹点的状态预测协方差矩阵主要是体现当前时刻航迹点的状态预测值的误差，通过调整当前时刻航迹点的状态预测协方差矩阵能够减少当前时刻航迹点的状态预测值的误差。具体地，在一个实施例中，所述根据航迹点的状态信息预测当前时刻航迹点的状态预测值以及计算当前航迹点的状态预测协方差矩阵通过下述方程式计算得到：P(k/k-1)＝Φ(k-1)P(k-1/k-1)ΦT(k-1)+Q(k-1)其中，为当前时刻航迹点的状态预测值，P(k/k-1)为当前时刻航迹点的状态预测协方差矩阵，为前一时刻航迹点的状态预测值，P(k-1/k-1)为前一时刻航迹点的状态预测协方差矩阵；所述航迹点状态信息包括前一时刻状态转移矩阵Φ(k-1)、前一时刻控制矩阵U(k-1)、前一时刻机动加速度均值前一时刻航迹点的状态预测协方差矩阵P(k/k-1)、前一时刻状态转移矩阵的转置矩阵ΦT(k-1)、前一时刻噪声方差矩阵Q(k-1)。进一步地：前一时刻状态转移矩阵α为机动频率，转弯机动频率α＝1/60，逃避机动频率α＝1/20，大气扰动机动频率α＝1；T为采样时间；e为欧拉数；前一时刻控制矩阵当前时刻噪声方差矩阵式中为k-2时刻的加速度预测值。步骤S103：根据航迹点的状态信息预测当前时刻航迹点的测量预测值以及计算当前时刻航迹点的残差向量协方差矩阵；所述当前时刻航迹点的测量预测值指的是当前时刻航迹点的位置信息，所述当前时刻航迹点的残差向量协方差矩阵主要是体现当前时刻航迹点的测量预测值的误差，通过调整当前时刻航迹点的残差向量协方差矩阵能够减少当前时刻航迹点的测量预测值的误差。具体地，在一个实施例中，所述根据航迹点的状态信息预测当前时刻航迹点的测量预测值以及计算当前时刻航迹点的残差向量协方差矩阵通过下述方程式计算得到：S(k)＝H(k)P(k/k-1)HT(k)+R(k)其中，为当前时刻航迹点的测量预测值，S(k)为当前时刻航迹点的残差向量协方差矩阵；所述航迹点状态信息还包括当前时刻测量矩阵H(k)、当前时刻测量矩阵的转置矩阵HT(k)、当前时刻观测噪声R(k)。较优地，在一个实施例中，当前时刻观测噪声R(k)的值可以为30，当前时刻测量矩阵H(k)＝[100]。步骤S104：识别当前时刻航迹点是否为野值，并根据当前时刻航迹点是否为野值确定当前时刻航迹点的衰减因子；所述野值是指该航迹点数据严重偏离目标真实值，即为野值；由于要对野值进行修正，所以先判定该航迹点是否为野值，再确定该航迹点的衰减因子。具体地，在一个实施例中，所述识别当前时刻航迹点是否为野值本文档来自技高网...

【技术保护点】
一种基于ADS‑B系统航迹点的修正方法，其特征在于，包括步骤:对航迹点按接收时间进行排列；根据航迹点的状态信息预测当前时刻航迹点的状态预测值以及计算当前时刻航迹点的状态预测协方差矩阵；根据航迹点的状态信息预测当前时刻航迹点的测量预测值以及计算当前时刻航迹点的残差向量协方差矩阵；识别当前时刻航迹点是否为野值，并根据当前时刻航迹点是否为野值确定当前时刻航迹点的衰减因子；根据当前时刻航迹点的衰减因子计算当前时刻滤波器增益；根据当前时刻航迹点的测量预测值、当前航迹点的残差向量协方差矩阵以及当前时刻滤波器增益修正当前时刻航迹点的状态预测值以及当前航迹点的状态预测协方差矩阵。

【技术特征摘要】
1.一种基于ADS-B系统航迹点的修正方法，其特征在于，包括步骤:对航迹点按接收时间进行排列；根据航迹点的状态信息预测当前时刻航迹点的状态预测值以及计算当前时刻航迹点的状态预测协方差矩阵，通过下述方程式计算得到：S(k)＝H(k)P(k/k-1)HT(k)+R(k)其中，为当前时刻航迹点的测量预测值，S(k)为当前时刻航迹点的残差向量协方差矩阵；所述航迹点状态信息还包括当前时刻测量矩阵H(k)、当前时刻测量矩阵的转置矩阵HT(k)、当前时刻观测噪声R(k)；根据航迹点的状态信息预测当前时刻航迹点的测量预测值以及计算当前时刻航迹点的残差向量协方差矩阵，通过下述方程式计算得到：P(k/k-1)＝Φ(k-1)P(k-1/k-1)ΦT(k-1)+Q(k-1)其中，为当前时刻航迹点的状态预测值，P(k/k-1)为当前时刻航迹点的状态预测协方差矩阵，为前一时刻航迹点的状态预测值，P(k-1/k-1)为前一时刻航迹点的状态预测协方差矩阵；所述航迹点状态信息包括前一时刻状态转移矩阵Φ(k-1)、前一时刻控制矩阵U(k-1)、前一时刻机动加速度均值当前时刻航迹点的状态预测协方差矩阵P(k/k-1)、前一时刻状态转移矩阵的转置矩阵ΦT(k-1)、前一时刻噪声方差矩阵Q(k-1)；识别当前时刻航迹点是否为野值，并根据当前时刻航迹点是否为野值确定当前时刻航迹点的衰减因子，通过下述方程式计算得到：其中,为当前时刻航迹点的新息，当时，当前时刻航迹点不是野值，当时，当前时刻航迹点为野值；Y(k)为当前时刻航迹点的测量值，为当前时刻航迹点的测量预测值，d设置为3，r(k)为当前时刻航迹点的衰减因子；根据当前时刻航迹点的衰减因子计算当前时刻滤波器增益，通过下述方程式计算得到：K(k)＝r(k)P(k)HT(k)S-1(k)其中，K(k)为当前时刻滤波器增益，S-1(k)为当前时刻航迹点的残差向量协方差矩阵的逆矩阵；根据当前时刻航迹点的测量预测值、当前航迹点的残差向量协方差矩阵以及当前时刻滤波器增益修正当前时刻航迹点的状态预测值以及当前时刻航迹点的状态预测协方差矩阵，通过下述方程式计算得到：P(k/k)＝P(k/k-1)-K(k)S(k)KT(k)其中，为修正后的当前时刻航迹点的状态预测值，P(k/k)为修正后的当前航迹点的状态预测协方差矩阵，KT(k)为当前时刻滤波器增益转置矩阵。2.根据权利要求1所述的基于ADS-B系统航迹点的修正方法，其特征在于，在对航迹点按接收时间进行排列之后，预测当前时刻航迹点的状态预测值以及计算当前时刻航迹点的状态预测协方差矩阵之前还包括步骤：将每两个航迹点的时间间隔与滤...

【专利技术属性】
技术研发人员：赵峙岳，刘晖，窦军华，陈德亚，蒲强，
申请(专利权)人：广州海格通信集团股份有限公司，
类型：发明
国别省市：广东;44