一种无人机实时轨迹预测方法与流程

文档序号：21643665发布日期：2020-07-29 02:56阅读：918来源：国知局

本发明属于空中交通流量管理技术领域，具体涉及一种无人机实时轨迹预测方法。

背景技术：

无人驾驶飞机简称“无人机”，英文缩写为“uav”，是利用无线电遥控设备和自备的程序控制装置操纵的不载人飞机，或者由车载计算机完全地或间歇地自主地操作。

无人机结构简单、机动性强、飞行状态多变，因此，其飞行状态及轨迹难以预测。随着无人机需求与数量的不断增加，为避免无人机之间的冲突与碰撞，保证众多的无人机在低空空域安全、有序地飞行，需要对无人机轨迹预测方法进一步研究，从而为无人机冲突检测、异常行为预警和空域状态科学评估提供决策支持。

目前，传统的无人机轨迹预测方法是基于无人机动力学模型的，通过建立动力学方程来解释无人机飞行时的力学原理，在一定的假设条件下，对无人机轨迹进行预测，这些预测方法需要复杂的动力学建模以及较多的假设与约束，而且不能充分利用无人机的历史信息。在无人机的轨迹预测的过程中，无人机的状态是多变的。因此，为了提高预测精度，首先需要对无人机的状态进行判断，再充分利用历史轨迹点信息对无人机的位置进行预测。

技术实现要素：

针对于上述现有技术的不足，本发明的目的在于提供一种无人机实时轨迹预测方法，以解决现有技术中无人机飞行状态难辨识、轨迹预测精度低、轨迹预测速度慢的问题；本发明的方法能捕捉无人机的飞行状态，充分利用历史轨迹点信息，避免无人机的冲突与碰撞，保证众多的无人机在低空空域安全、有序地飞行。

为达到上述目的，本发明采用的技术方案如下：

本发明的一种无人机实时轨迹预测方法，步骤如下：

(1)获取数据：采集无人机历史广播式自动相关监视系统(ads-b)数据，得到无人机时间戳、无人机id、纬度、经度、高度、水平速度、飞行时间的数据信息；

(2)数据预处理：将上述获取到的数据信息按飞行时间进行排序，将纬度和经度转换成直角坐标；并删除获取到的数据中非固定时间间隔和存在不明显位置信息的数据；

(3)生成无人机轨迹各个变量数据集：对上述步骤(2)预处理后的数据进行筛选，得到纬度、经度、高度和水平速度数据，并计算无人机的垂直速度，生成相应的时间序列数据；

(4)利用马尔科夫模型确定无人机状态；

(5)根据无人机的状态，增添bn层，建立基于lstm网络的无人机轨迹预测模型；

(6)预测无人机的经度、纬度及高度。

进一步地，所述步骤(3)中纬度、经度、高度、水平速度的时间序列由ads-b数据预处理及筛选后得到，垂直速度的时间序列通过计算高度变化值除以飞行时间间隔得到，最终的无人机轨迹数据集包括纬度、经度、高度、水平速度和垂直速度五个变量。

进一步地，所述步骤(4)的具体过程如下：

(41)对数据集进行分组：

将步骤(3)中的数据集按不同的时间间隔分组，生成相应时间间隔的无人机轨迹时间序列数据；

(42)划分无人机的飞行状态：

根据无人机的水平速度vx和垂直速度vy，对无人机的飞行状态进行划分；当vy＞0.1米/秒，无人机处于上升状态s1；当米/秒且vy∈[-0.1,0.1]米/秒，无人机处于平飞状态s2；当vy＜-0.1米/秒，无人机处于下降状态s3；当vx∈[-0.1,0.1]米/秒且vy∈[-0.1,0.1]米/秒，无人机处于悬停状态s4；

(43)计算初始状态概率：

假设无人机速度的时间序列的取值空间为(vxi,vyi),i＝1,2…n，一共有n个观测值，根据无人机速度的时间序列，无人机被划分为4个状态s＝{s1,s2,s3,s4}；(vxn,vyn)作为最后一个速度的观测值，不考虑(vxn,vyn)在n时刻对应的转移状态；假设根据无人机速度的时间序列的前n-1个已知值，有ci个无人机处于状态si(i＝1,2,3,4)，有那么无人机状态si的初始状态概率为

(44)计算一步转移概率矩阵：

根据已知的无人机速度的时间序列{(vxt,vyt),t∈t}，有ci个无人机处于状态si，假设其中下一时刻转移到状态sj的数据有cij个；则无人机从状态si转移到状态sj的一步转移概率为由pij组成的无人机状态转移矩阵p为一步状态转移概率矩阵；

(45)无后效性检验：

采用χ²统计量检验无人机速度的时间序列的无后效性；无人机有4种飞行状态，将边际概率定义为：

统计量服从自由度为(m-1)²的χ²分布；选定置信度α，查χ²分布临界值表得令：

若则无人机速度的时间序列具有无后效性；

(46)利用一步转移概率矩阵预测无人机飞行状态：

若在t时刻无人机处于状态si，且一步状态转移概率矩阵p的第i行无人机状态转移概率的最大值为pij，则在t+1时刻，无人机处于状态sj，故预测在t+1时刻，无人机将转移到状态sj。

进一步地，所述步骤(5)的具体过程如下：

(51)初始化lstm(长短期记忆网络)网络各层参数；

(511)初始化lstm输入层参数；

当搭建无人机轨迹预测lstm网络的输入层时，每次批量训练40个样本，与步骤(41)对应，时间步长初始值设为10、20、30、40、50；

(512)初始化lstm隐含层参数；

当搭建lstm网络的隐含层时，隐含层层数设置为2，每层50个神经元；

(513)初始化lstm网络层参数；

激活函数选择sigmoid函数和tanh函数，初始化网络层的权重和偏置；

(514)损失函数设置；

在无人机轨迹预测lstm网络中采用平方差损失函数；

(52)确定lstm的输入样本特征与输出样本特征；

(521)构建输入样本特征；

将步骤(3)得到的五个变量的数据集：纬度、经度、高度、水平速度和垂直速度，按照步骤(511)设置的时间步长划分为时间序列；

(522)构建输出样本特征；

神经网络输出样本特征为时间序列最后一个时刻的纬度、经度、高度；

(523)数据标准化；

输入样本特征通过增添的bn层，进行批量离差标准化。

进一步地，所述步骤(6)的具体过程如下：

(61)模型训练及测试；

随机抽取上述步骤(3)得到的无人机轨迹数据集的75％作为训练集，剩余的25％作为测试集；

(62)优化器选择；

采用adam优化器，并使学习率自动调整；

(63)数据去标准化；

对lstm模型的输出值，通过bn层，进行逆离差标准化，得到的纬度、经度、高度即为预测的无人机轨迹点；

(64)预测精度评价；

采用平均绝对误差(mae)作为评价指标。

本发明的有益效果：

本发明的方法能够将马尔科夫模型与长短期记忆神经网络模型相结合，在判断无人机飞行状态的前提下，预测无人机的轨迹，提高无人机轨迹的预测精度。较传统的预测方法，既充分利用历史信息，避免复杂的动力学建模，又能捕捉无人机状态。本方法能为无人机冲突检测、异常行为预警和空域状态科学评估提供决策支持，进一步完善无人机管理体系。

附图说明

图1为本发明方法的流程图。

图2为无人机轨迹数据示例图。

图3为本发明方法的具体实现过程示意图。

具体实施方式

为了便于本领域技术人员的理解，下面结合实施例与附图对本发明作进一步的说明，实施方式提及的内容并非对本发明的限定。

参照图1、图3所示，本发明的一种无人机实时轨迹预测方法，步骤如下：

(1)获取数据：采集无人机历史ads-b数据，参见图2所示，得到无人机时间戳、无人机id、纬度、经度、高度、水平速度、飞行时间的数据信息。

(2)数据预处理：将上述获取到的数据信息按飞行时间进行排序，将纬度和经度转换成直角坐标；并删除获取到的数据中非固定时间间隔和存在不明显位置信息的数据。

所述步骤(3)中纬度、经度、高度、水平速度的时间序列由ads-b数据预处理及筛选后得到，垂直速度的时间序列通过计算高度变化值除以飞行时间间隔得到，最终的无人机轨迹数据集包括纬度、经度、高度、水平速度和垂直速度五个变量。

(4)利用马尔科夫模型确定无人机状态，具体过程如下：