CN114819068A - 一种混合型目标航迹预测方法及系统 - Google Patents
一种混合型目标航迹预测方法及系统 Download PDFInfo
- Publication number
- CN114819068A CN114819068A CN202210372743.6A CN202210372743A CN114819068A CN 114819068 A CN114819068 A CN 114819068A CN 202210372743 A CN202210372743 A CN 202210372743A CN 114819068 A CN114819068 A CN 114819068A
- Authority
- CN
- China
- Prior art keywords
- target
- target track
- track prediction
- prediction
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000013598 vector Substances 0.000 claims description 44
- 238000004422 calculation algorithm Methods 0.000 claims description 42
- 238000001914 filtration Methods 0.000 claims description 36
- 238000012549 training Methods 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 27
- 230000006870 function Effects 0.000 claims description 21
- 230000015654 memory Effects 0.000 claims description 11
- 210000002569 neuron Anatomy 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 238000003062 neural network model Methods 0.000 claims description 6
- 230000007787 long-term memory Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000001537 neural effect Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 230000006872 improvement Effects 0.000 abstract description 2
- 238000013528 artificial neural network Methods 0.000 description 18
- 238000004088 simulation Methods 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 7
- 238000013135 deep learning Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 210000004027 cell Anatomy 0.000 description 3
- 238000010801 machine learning Methods 0.000 description 3
- 230000006403 short-term memory Effects 0.000 description 3
- 230000004913 activation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000009194 climbing Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000008034 disappearance Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000306 recurrent effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 208000025174 PANDAS Diseases 0.000 description 1
- 208000021155 Paediatric autoimmune neuropsychiatric disorders associated with streptococcal infection Diseases 0.000 description 1
- 240000000220 Panda oleosa Species 0.000 description 1
- 235000016496 Panda oleosa Nutrition 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000004793 poor memory Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/047—Probabilistic or stochastic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/049—Temporal neural networks, e.g. delay elements, oscillating neurons or pulsed inputs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及目标轨迹预测技术领域，公开了一种混合型目标航迹预测方法及系统，一种混合型目标航迹预测方法，将卡尔曼滤波算法和长短时记忆神经网络模型相结合，对目标轨迹进行预测。本发明解决了现有技术存在的飞行目标轨迹预测的有效性和准确性还有较大提升空间等问题。
Description
技术领域
本发明涉及目标轨迹预测技术领域，具体是一种混合型目标航迹预测方法及系统。
背景技术
目标轨迹预测是保证目标航行安全、规划飞行航迹和搜寻空中目标等任务的关键技术，作为航空、交通管制、军事等领域的重点研究问题，一直受到学界和工业界广泛关注。随着近现代信息科学技术的不断发展，目标轨迹预测的方法和模型也在不断地更新和改进。目标轨迹预测即利用当前时刻目标的状态信息预测目标下一时刻的状态，该任务可等效地看作为针对时间序列的预测问题。飞行目标轨迹预测模型分为动力学模型、卡尔曼滤波算法模型和可学习模型。
动力学模型根据目标运动的空气动力学方程，对目标的轨迹进行预测，但在现实中，飞行目标往往受到天气、环境等诸多因素的影响，而且目标飞行过程所受力复杂多变，根据空气动力学建立的目标运动模型对实际情况进行了较大程度的简化，由此导致这类模型对目标轨迹预测的精度较低。
卡尔曼滤波器(Kalman filter,KF)是一种应用非常广泛且有效的目标轨迹预测方法。SJ Julier等在1997年提出扩展卡尔曼滤波(extend Kalman filter,EKF)算法，不仅继承了卡尔曼滤波算法的优势，而且能够处理非线性系统模型。吕波等利用扩展卡尔曼滤波嵌入非线性控制系统，针对航空器的航迹预测任务建立了多信息扩展卡尔曼滤波预测模型。Chen等提出了基于无迹卡尔曼滤波的航迹预测模型并完成了对飞行目标的航迹预测。然而，在现实的许多目标轨迹预测任务中，运动模型一般不能预先指定，在这种情况下，卡尔曼滤波器的应用受到了严重的限制。此外，卡尔曼滤波器通常是在恒定加速度的假设下应用的，这也是对现实目标运动情况的粗略逼近，这样建立的算法模型对目标轨迹预测的精度有限。为了克服传统卡尔曼滤波算法的局限性，许多学者尝试直接从训练数据中学习运动模型，利用可学习模型来学习目标运动方程可以避免手工设计卡尔曼滤波转移和测量矩阵，从而提高目标轨迹预测精度。
可学习模型(即机器学习模型)包括以神经网络为代表的机器学习方法，神经网络可通过大量数据学习数据的本质特征，其深层非线性的网络结构能够实现对复杂函数的逼近，在目标轨迹预测方面得到了广泛应用。如钱夔等用BP神经网络对飞机航迹数据进行训练学习，实现对飞行目标轨迹的预测，但BP神经网络比较简单，预测精度有限，Wu等提出了基于卷积神经网络的飞行航迹预测模型，根据飞行目标的高度和经纬度信息预测飞行目标的位置信息。
以上所提的模型对目标轨迹预测各有优劣，但都没有利用目标轨迹数据在时序上具有的高度相关性。循环神经网络(recurrent neural network,RNN)是一种能够提取序列数据特征的经典神经网络模型，然而由于RNN自身简单的网络结构对于长时间序列的记忆能力较差，于是Hochreiter等提出了由RNN改进的长短时记忆(Long Short Memory,LSTM)神经网络，LSTM模型由于其自身复杂的结构使其能够记忆更长的序列数据中蕴含的信息，并在后续处理中加以运用，如今LSTM网络已广泛应用在各类目标的轨迹预测任务中并取得了良好的效果。
以上现有技术飞行目标轨迹预测的有效性和准确性还有较大提升空间，亟待改善。
发明内容
为克服现有技术的不足，本发明提供了一种混合型目标航迹预测方法及系统，解决现有技术存在的飞行目标轨迹预测的有效性和准确性还有较大提升空间等问题。
本发明解决上述问题所采用的技术方案是：
一种混合型目标航迹预测方法，将卡尔曼滤波算法和长短时记忆神经网络模型相结合，对目标轨迹进行预测。
作为一种优选的技术方案，其特征在于，包括以下步骤：
S1，构建目标轨迹预测网络模型：
首先，根据卡尔曼滤波算法，建立目标轨迹预测模型，如下：
yt＝f(yt-1)+wt (18)，
zt＝yt+vt (19)，
其中，yt为目标状态向量，zt为观测向量，wt为过程噪声，vt为观测噪声，f(·)为状态转移函数；
然后，目标轨迹预测模型根据目标前一时刻的状态估计值预测出下一时刻目标的状态，同时更新当前预测状态的协方差矩阵，预测公式为：
其中，
为目标前一时刻的状态向量估计值，
为目标当前时刻状态向量估计值，
为当前时刻状态向量估计值的协方差矩阵，F为f(·)关于
的雅克比矩阵，
为可学习参数，T表示矩阵的转置；
再然后，更新计算得出目标在t时刻的最终状态向量估计值
更新公式为：
其中，Kt为卡尔曼增益，
为可学习参数，
是目标在t时刻的观测值，
为目标在t时刻的最终状态向量估计值，
为目标在t时刻的协方差矩阵；
S2：训练目标轨迹预测网络模型：对目标轨迹预测网络模型进行训练，直至训练次数≥设定的次数阈值；
S3：输出目标轨迹预测结果：将训练后的目标轨迹预测网络模型中的目标轨迹输出。
作为一种优选的技术方案，步骤S1中，模型中的长短时记忆神经模块包括双隐含层LSTM网络层，隐层神经元个数为64。
作为一种优选的技术方案，步骤S2中，目标轨迹预测网络模型的损失函数为：在均方误差函数的基础上增加目标状态估计的惩罚项。
作为一种优选的技术方案，损失函数公式为：
其中，θ代表网络中的所有参数，λ为超参数，T表示训练所用序列总点数。
作为一种优选的技术方案，步骤S2中，训练策略为：网络中所有的权重矩阵用Xavier随机初始化，偏差均初始化为0，学习率设定为0.0001，采用Adam优化器，以及，设定训练批尺寸为16。
作为一种优选的技术方案，步骤S2中，取80％的数据点作为目标轨迹预测网络模型的训练数据，剩下20％的数据点作为目标轨迹预测网络模型的测试数据。
作为一种优选的技术方案，还包括设于步骤S1之前的以下步骤：
S0，数据预处理：对输入的样本数据进行标准化处理，公式如下：
其中，x为原样本向量，
为样本均值，σ为样本标准差，y为经标准化处理后的样本向量。
作为一种优选的技术方案，还包括设于步骤S3之后的以下步骤：
S4，反标准化处理：在对目标轨迹预测结束后，将输出的预测结果
进行反标准化处理，公式如下：
一种混合型目标航迹预测系统，基于所述的一种混合型目标航迹预测方法，包括依次电相连的目标轨迹预测网络模型构建模块、目标轨迹预测网络模型训练模块、目标轨迹预测结果输出模块；
目标轨迹预测网络模型构建模块：用以构建目标轨迹预测网络模型，构建步骤如下：
首先，根据卡尔曼滤波算法，建立目标轨迹预测模型，如下：
yt＝f(yt-1)+wt (28)，
zt＝yt+vt (29)，
其中，yt为目标状态向量，zt为观测向量，wt为过程噪声，vt为观测噪声，f(·)为状态转移函数；
然后，目标轨迹预测模型根据目标前一时刻的状态估计向量预测出下一时刻目标的状态向量，同时更新当前预测状态的协方差矩阵，预测公式为：
其中，
为目标前一时刻的状态向量估计值，
为目标当前时刻的状态向量估计值，
为当前预测状态的协方差矩阵，为f(·)关于
的雅克比矩阵，
为可学习参数，T表示矩阵的转置；
再然后，更新计算得出目标在t时刻的最终状态向量估计值
更新公式为：
其中，Kt为卡尔曼增益，
为可学习参数，
是目标在t时刻的观测值，
为目标在t时刻的最终状态向量估计值，
为目标在t时刻的协方差矩阵；
目标轨迹预测网络模型训练模块：用以对目标轨迹预测网络模型进行训练，直至训练次数≥设定的次数阈值；
目标轨迹预测结果输出模块：用以将训练后的目标轨迹预测网络模型中的目标轨迹输出。
本发明相比于现有技术，具有以下有益效果：
(1)本发明构建了传统卡尔曼滤波算法和深度学习之间联系的桥梁，将卡尔曼滤波算法和深度网络有机地结合起来，避免了显式建模目标运动轨迹方程，还能利用深度神经网络强大的学习表征能力从数据中学习模型中的超参数；
(2)本发明克服了传统模型建模不准确和深度学习方法对数据的严重依赖性；
(3)本发明提高了模型与算法的可解释性，加快了网络的收敛速度；
(4)本发明实现了对目标复杂的非线性运动轨迹的良好估计和预测。
附图说明
图1为现有技术的LSTM网络结构图；
图2为本发明KFU-LSTM-Net结构图；
图3为本发明KFU-LSTM模块内部结构图；
图4为卡尔曼滤波算法展开图；
图5为LSTM模块结构图；
图6为利用本发明对目标轨迹进行预测的一个实施例的流程图；
图7为本发明所述的一种混合型目标航迹预测方法的步骤示意图。
具体实施方式
下面结合实施例及附图，对本发明作进一步的详细说明，但本发明的实施方式不限于此。
实施例1
如图1至图7所示，将卡尔曼滤波算法和长短时记忆神经网络模型相结合，对目标轨迹进行预测。
卡尔曼滤波算法描述如下：
1.1卡尔曼滤波算法
卡尔曼滤波算法是一种基于时域离散的自回归优化算法，其被广泛应用于动态目标的位置、速度等状态信息的预测任务中。设目标状态向量为yt，观测变量为zt，根据卡尔曼滤波算法，有如下关系式：
yt＝Ayt-1+wt (35)，
zt＝Hyt+vt (36)，
式中：A是目标状态转移矩阵，代表该系统从t-1时刻到t时刻状态之间的转移关系；H是目标观测矩阵，wt和vt是过程噪声和观测噪声，其分别满足正态分布wt～N(0，Qt)，vt～N(0，Rt)。
卡尔曼滤波算法通过反馈环路迭代更新求解目标状态向量，该反馈环路包含两个更新步骤：预测步骤和更新步骤。在预测步骤中，卡尔曼滤波算法用上一时刻目标的状态向量来预测目标下一时刻的状态：
在更新步骤中，需要首先计算卡尔曼增益Kt，然后通过观测变量zt来修正目标在t时刻的状态估计值
和协方差矩阵
1.2长短时记忆神经网络
不同于前馈神经网络，递归神经网络(RNN)是针对于时序数据进行处理的网络模型。然而，由于梯度消失问题，早期的RNN模型都很难被训练。为了解决梯度消失或梯度爆炸问题并使模型可以长期存储时间序列中的有用信息，Hochreiter等人提出了长短时记忆神经网络(LSTM)，LSTM是RNN神经网络的一种特殊的变体，其更改了RNN的神经元结构，在神经元内部增加了遗忘门、输入门、输出门等单元，由此可以通过LSTM神经元获得长时间序列的依赖信息。LSTM网络结构如图1所示。
长短期记忆(LSTM)神经网络是专门为解决长时间序列难以记忆和梯度消失的问题而引入的，并已成为最广泛使用的网络架构之一。LSTM网络的隐藏层中包含了遗忘门、输入门、输出门，它们不仅控制着信息的流动模式，而且对网络预测的性能方面至关重要。LSTM神经元内部计算公式如下：
ft＝σ(Wfhht-1+Wfxxt+bf) (42)
it＝σ(Wihht-1+Wixxt+bi) (43)
ot＝σ(Wohht-1+Woxxt+bo) (44)
ht＝ot⊙tanh(ct) (47)
上式中σ(·)表示sigmoid激活函数，⊙代表相同位置元素对应相乘。
(1)遗忘门：将xt和ht-1作为输入，经过sigmoid激活函数σ将信息输入到细胞c中，由此来控制细胞对短期记忆的遗忘程度，遗忘门的计算如公式(8)所示。
(2)输入门：输入门包含两个部分，第一部分负责产生新记忆信息
第二部分负责确定新记忆信息的流入比例it，分别如公式(9)和公式(11-12)所示。
(3)输出门：输出门通过将细胞状态ct以及输出比例ot相乘，输出得到该LSTM神经元最终的输出值ht，如公式(10)和公式(13)所示。
2基于卡尔曼滤波算法展开的LSTM网络
2.1算法模型
在本节中，本发明将卡尔曼滤波算法和长短时记忆神经网络模型相结合，提出了基于卡尔曼滤波算法展开的长短时记忆神经网络(KFU-LSTM-Net)模型来对目标轨迹进行预测。
根据卡尔曼滤波算法，本发明建立的目标轨迹预测模型为如下形式：
yt＝f(yt-1)+wt (48)
zt＝yt+vt (49)
模型根据目标前一时刻的状态估计值
预测出下一时刻目标的状态
同时更新当前预测状态的协方差矩阵
模型的预测步骤公式为：
公式(16)中的映射函数f(·)通过LSTM网络学习拟合，利用LSTM强大的非线性拟合能力，使得模型可以学习到目标复杂的运动轨迹；公式(17)中F是f关于
的雅克比矩阵，
为模型中的可学习参数，T表示矩阵的转置。
模型的更新步骤公式为：
公式(18)中
是模型中的可学习参数，式(19)中
是目标在t时刻的观测值。
本发明所建网络模型整体架构如图2所示，模型中的KFU-LSTM模块内部结构如图3所示。在KFU-LSTM模块中本发明设计了LSTM网络来学习状态转移函数f(·)，即将上一时刻目标的状态向量
作为LSTM网络的输入得到目标在t时刻的初步状态估计值
再将
和测量状态向量
作为输入，通过卡尔曼滤波算法更新计算得出目标在t时刻的最终状态预测向量
从算法展开的角度来看，在每一个时间步t，LSTM将上一时刻(即t-1时刻)的目标状态估计值
作为输入计算得到
由协方差矩阵
和
根据公式(17-18)计算出卡尔曼增益Kt；接着，结合LSTM输出的目标状态初步估计值
和观测值
根据公式(19)最终计算得出目标在t时刻的状态估计值
卡尔曼滤波算法展开示意图如图4所示。
2.2网络模型结构及训练参数设定
(1)本发明所提网络结构及参数：本发明所建网络模型整体架构如图2所示，模型中的KFU-LSTM模块内部结构如图3所示。对于KFU-LSTM中的LSTM模块，本发明建立双隐含层的LSTM网络，其中隐层神经元个数为64，其后添加两个全连接层。具体LSTM的结构如图5所示。
(2)网络的损失函数：不同于标准的均方误差损失函数，本发明在均方误差函数的基础上增加了目标状态估计的惩罚项，以此来加强网络对于目标状态的估计能力，本发明所提出的损失函数定义如下：
式(21)中θ代表网络中的所有参数，本发明实验中超参数λ设定值为0.5。
(3)训练策略：网络中所有的权重矩阵用Xavier随机初始化，偏差均初始化为0。学习率设定为0.0001，采用Adam优化器，设定训练批尺寸为16。本发明所有实验均在Windows10(64bit)操作系统的Pytorch平台中运行，电脑设备硬件配置为Intel Core i7-9750H CPU和NVIDIA GTX 1660Ti。实验环境如表1所示。
表1实验环境列表
本发明所提KFU-LSTM-Net模型针对目标轨迹预测整体流程如图6所示。
3实验及结果分析
为了验证本发明所提出的基于卡尔曼滤波算法展开的深度LSTM网络模型对于目标轨迹预测的准确性和有效性，本发明模型对目标轨迹仿真数据进行预测，根据目标轨迹预测结果图和预测误差进行了详细的分析，并在实验部分对比了基于RNN和LSTM网络的目标轨迹预测方法。
3.1数据描述
(1)目标轨迹仿真数据为三种带有随机高斯噪声的三维曲线，其中80％的数据点作为网络模型的训练数据，剩下20％的数据点作为模型的测试数据。参数方程式及取值范围如下表2所示：
表2仿真数据的参数方程及变量取值
(2)实验使用的真实数据是在FlightAware上的航线日志，在FlightAware平台上导出csv格式的飞行数据再利用Python的Pandas库加载实验数据。飞机航迹分为爬升、转弯和平飞三个阶段，数据包含200个采集点，每个采集点包含时间、纬度、经度、高度、航向、节、航速、爬升角共八维数据，本发明从中选择经度、纬度和高度这三个航迹指标作为网络模型的输入进行目标的航迹预测。取80％的数据点作为网络模型的训练数据，剩下20％的数据点作为模型的测试数据。
3.2数据预处理
为了消除输入数据量纲的影响和提高网络的收敛速率，本发明对输入的样本数据进行标准化处理，以提高网络模型预测的准确性和鲁棒性，公式如下：
式(22)中x为原样本向量，
为样本均值，σ为样本标准差，y为经标准化处理后的样本向量。
在对目标轨迹预测结束后，可根据如下公式(23)将网络输出的预测结果
进行反标准化处理，以便同目标轨迹真实值相比较：
3.3对比实验
为验证本发明所提KFU-LSTM-Net模型对于目标轨迹预测的有效性和优越性，本发明针对三种不同的飞行仿真数据和真实飞行航迹数据进行实验，分别用RNN、LSTM和本发明所提KFU-LSTM-Net模型进行预测，RNN和LSTM模型均采用2层隐含层且隐层神经元个数为64的网络架构，训练采用Adam训练策略，学习率设定为0.001。
为了测试本发明所提模型对目标状态预测的精度和不同场景下的表现情况，本发明考察了模型对目标轨迹的预测值与真实值的均方误差(root mean square error,RMSE)指标，各模型在仿真数据和真实数据集上预测的RMSE对比表格如表3所示。从表3中结果可以直观地得出，本发明所提KFU-LSTM模型较RNN和LSTM模型的RMSE有明显减小，KFU-LSTM模型在三个仿真数据和真实数据集上对飞行目标轨迹预测的精度都超过了RNN和LSTM模型。由此表明了本发明所提的模型在目标轨迹预测精度上的显著提升，验证了本发明KFU-LSTM模型相较RNN和LSTM网络在目标轨迹时序预测上的优越性和有效性。
表3模型目标轨迹预测结果的均方误差
本发明提出并建立了卡尔曼滤波算法展开的LSTM神经网络模型，并将该模型应用到飞行目标轨迹预测的任务中。本发明所建立的模型能够学习飞行目标复杂的运动模型。相较于传统的卡尔曼滤波算法只能用于估计线性运动模型，本发明所提模型将目标状态转移函数建模为可学习的网络，再利用卡尔曼滤波对目标预测状态进行修正更新，因此能够对目标复杂的非线性运动轨迹进行良好地估计和预测。另一方面，通过以卡尔曼滤波算法更新的方式将LSTM预测和观测向量结合起来，即隐式学习了目标所有可能的运动路径。不仅避免了显式建模目标运动轨迹方程，还能利用深度神经网络强大的学习表征能力从数据中学习模型中的超参数。相较于RNN、LSTM等为代表的深度学习方法，本发明建立的模型在目标轨迹预测的精度更高。
本发明构建了传统卡尔曼滤波算法和深度学习之间联系的桥梁，将卡尔曼滤波算法和深度网络有机地结合起来，不仅克服了传统模型建模不准确和深度学习方法对数据的严重依赖性，而且提高了模型与算法的可解释性，加快了网络的收敛速度，在仿真数据和真实数据上的实验结果验证了本发明所提的基于卡尔曼滤波算法展开网络模型对飞行目标轨迹预测的有效性和优越性。
本发明提出了基于卡尔曼滤波算法展开的长短时记忆神经网络(Kalman filterunfolding LSTM network,KFU-LSTM-Net)。所提的模型将目标状态转移函数建模为可学习的网络，再利用卡尔曼滤波对目标状态向量进行修正更新得出目标状态预测值。本发明所建模型从目标轨迹数据中学习目标运动轨迹和卡尔曼滤波器中的超参数，让模型能够更好地利用深度学习的强大表征能力，在仿真和真实目标轨迹数据实验上验证了本发明所提模型对于飞行目标轨迹预测的有效性和准确性。
针对传统飞行目标轨迹预测方法模型较为简化且预测精度较低的问题，提出了一种数据和模型混合驱动的目标航迹预测方法，将基于卡尔曼滤波算法展开的深度神经网络(KFU-LSTM-Net)模型用于飞行目标的轨迹预测任务。该模型通过长短时记忆(long short-term memory,LSTM)网络从目标的航迹数据中学习目标的运动状态，再利用卡尔曼滤波算法对LSTM预测的目标状态估计值进行动态修正，其有效结合了卡尔曼滤波算法和深度神经网络各自的优势。在仿真数据和真实数据上的实验验证了所提KFU-LSTM-Net模型较其他网络模型对飞行目标轨迹预测的准确性和有效性。
如上所述，可较好地实现本发明。
本说明书中所有实施例公开的所有特征，或隐含公开的所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以以任何方式组合和/或扩展、替换。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，依据本发明的技术实质，在本发明的精神和原则之内，对以上实施例所作的任何简单的修改、等同替换与改进等，均仍属于本发明技术方案的保护范围之内。
Claims (10)
1.一种混合型目标航迹预测方法，其特征在于，将卡尔曼滤波算法和长短时记忆神经网络模型相结合，对目标轨迹进行预测。
2.根据权利要求1所述的一种混合型目标航迹预测方法，其特征在于，包括以下步骤：
S1，构建目标轨迹预测网络模型：
首先，根据卡尔曼滤波算法，建立目标轨迹预测模型，如下：
yt＝f(yt-1)+wt (1)，
zt＝yt+vt (2)，
其中，yt为目标状态向量，zt为观测向量，wt为过程噪声，vt为观测噪声，f(·)为状态转移函数；
然后，目标轨迹预测模型根据目标前一时刻的状态估计值预测出下一时刻目标的状态，同时更新当前预测状态的协方差矩阵，预测公式为：
其中，
为目标前一时刻的状态向量估计值，
为目标当前时刻状态向量估计值，
为当前时刻状态向量估计值的协方差矩阵，F为f(·)关于
的雅克比矩阵，
为可学习参数，T表示矩阵的转置；
再然后，更新计算得出目标在t时刻的最终状态向量估计值
更新公式为：
其中，Kt为卡尔曼增益，
为可学习参数，
是目标在t时刻的观测值，
为目标在t时刻的最终状态向量估计值，
为目标在t时刻的协方差矩阵；
S2：训练目标轨迹预测网络模型：对目标轨迹预测网络模型进行训练，直至训练次数≥设定的次数阈值；
S3：输出目标轨迹预测结果：将训练后的目标轨迹预测网络模型中的目标轨迹输出。
3.根据权利要求2所述的一种混合型目标航迹预测方法，其特征在于，步骤S1中，模型中的长短时记忆神经模块包括双隐含层LSTM网络层，隐层神经元个数为64。
4.根据权利要求3所述的一种混合型目标航迹预测方法，其特征在于，步骤S2中，目标轨迹预测网络模型的损失函数为：在均方误差函数的基础上增加目标状态估计的惩罚项。
5.根据权利要求4所述的一种混合型目标航迹预测方法，其特征在于，损失函数公式为：
其中，θ代表网络中的所有参数，λ为超参数，T表示训练所用序列总点数。
6.根据权利要求5所述的一种混合型目标航迹预测方法，其特征在于，步骤S2中，训练策略为：网络中所有的权重矩阵用Xavier随机初始化，偏差均初始化为0，学习率设定为0.0001，采用Adam优化器，以及，设定训练批尺寸为16。
7.根据权利要求6所述的一种混合型目标航迹预测方法，其特征在于，步骤S2中，取80％的数据点作为目标轨迹预测网络模型的训练数据，剩下20％的数据点作为目标轨迹预测网络模型的测试数据。
8.根据权利要求2至7任一项所述的一种混合型目标航迹预测方法，其特征在于，还包括设于步骤S1之前的以下步骤：
S0，数据预处理：对输入的样本数据进行标准化处理，公式如下：
其中，x为原样本向量，
为样本均值，σ为样本标准差，y为经标准化处理后的样本向量。
9.根据权利要求8所述的一种混合型目标航迹预测方法，其特征在于，还包括设于步骤S3之后的以下步骤：
S4，反标准化处理：在对目标轨迹预测结束后，将输出的预测结果
进行反标准化处理，公式如下：
10.一种混合型目标航迹预测系统，其特征在于，基于权利要求2至9任一项所述的一种混合型目标航迹预测方法，包括依次电相连的目标轨迹预测网络模型构建模块、目标轨迹预测网络模型训练模块、目标轨迹预测结果输出模块；
目标轨迹预测网络模型构建模块：用以构建目标轨迹预测网络模型，构建步骤如下：
首先，根据卡尔曼滤波算法，建立目标轨迹预测模型，如下：
yt＝f(yt-1)+wt (11)，
zt＝yt+vt (12)，
其中，yt为目标状态向量，zt为观测向量，wt为过程噪声，vt为观测噪声，f(·)为状态转移函数；
然后，目标轨迹预测模型根据目标前一时刻的状态估计向量预测出下一时刻目标的状态向量，同时更新当前预测状态的协方差矩阵，预测公式为：
其中，
为目标前一时刻的状态向量估计值，
为目标当前时刻的状态向量估计值，
为当前预测状态的协方差矩阵，F为f(·)关于
的雅克比矩阵，
为可学习参数，T表示矩阵的转置；
再然后，更新计算得出目标在t时刻的最终状态向量估计值
更新公式为：
其中，Kt为卡尔曼增益，
为可学习参数，
是目标在t时刻的观测值，
为目标在t时刻的最终状态向量估计值，
为目标在t时刻的协方差矩阵；
目标轨迹预测网络模型训练模块：用以对目标轨迹预测网络模型进行训练，直至训练次数≥设定的次数阈值；
目标轨迹预测结果输出模块：用以将训练后的目标轨迹预测网络模型中的目标轨迹输出。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210372743.6A CN114819068A (zh) | 2022-04-11 | 2022-04-11 | 一种混合型目标航迹预测方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210372743.6A CN114819068A (zh) | 2022-04-11 | 2022-04-11 | 一种混合型目标航迹预测方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114819068A true CN114819068A (zh) | 2022-07-29 |
Family
ID=82534446
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210372743.6A Pending CN114819068A (zh) | 2022-04-11 | 2022-04-11 | 一种混合型目标航迹预测方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114819068A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115512221A (zh) * | 2022-09-22 | 2022-12-23 | 中国人民解放军海军航空大学 | 基于gnn的同步航迹鲁棒性关联方法 |
| CN116050278A (zh) * | 2023-02-17 | 2023-05-02 | 中国兵器装备集团自动化研究所有限公司 | 一种目标轨迹预测方法、装置、设备及存储介质 |
| CN117192063A (zh) * | 2023-11-06 | 2023-12-08 | 山东大学 | 基于耦合卡尔曼滤波数据同化的水质预测方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110232169A (zh) * | 2019-05-09 | 2019-09-13 | 北京航空航天大学 | 基于双向长短时记忆模型和卡尔曼滤波的轨迹去噪方法 |
| CN111829505A (zh) * | 2020-01-08 | 2020-10-27 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 多传感器航迹质量外推航迹融合方法 |
| JP6873519B1 (ja) * | 2020-04-24 | 2021-05-19 | 中国人民解放軍国防科技大学 | 軌跡予測方法及びシステム |
| CN113408392A (zh) * | 2021-06-15 | 2021-09-17 | 西安电子科技大学 | 基于卡尔曼滤波和神经网络的航迹补全方法 |
-
2022
- 2022-04-11 CN CN202210372743.6A patent/CN114819068A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110232169A (zh) * | 2019-05-09 | 2019-09-13 | 北京航空航天大学 | 基于双向长短时记忆模型和卡尔曼滤波的轨迹去噪方法 |
| CN111829505A (zh) * | 2020-01-08 | 2020-10-27 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 多传感器航迹质量外推航迹融合方法 |
| JP6873519B1 (ja) * | 2020-04-24 | 2021-05-19 | 中国人民解放軍国防科技大学 | 軌跡予測方法及びシステム |
| CN113408392A (zh) * | 2021-06-15 | 2021-09-17 | 西安电子科技大学 | 基于卡尔曼滤波和神经网络的航迹补全方法 |
Non-Patent Citations (1)
| Title |
|---|
| HUSEYIN COSKUN等: "Long Short-Term Memory Kalman Filters:Recurrent Neural Estimators for Pose Regularization" * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115512221A (zh) * | 2022-09-22 | 2022-12-23 | 中国人民解放军海军航空大学 | 基于gnn的同步航迹鲁棒性关联方法 |
| CN115512221B (zh) * | 2022-09-22 | 2024-02-27 | 中国人民解放军海军航空大学 | 基于gnn的同步航迹鲁棒性关联方法 |
| CN116050278A (zh) * | 2023-02-17 | 2023-05-02 | 中国兵器装备集团自动化研究所有限公司 | 一种目标轨迹预测方法、装置、设备及存储介质 |
| CN117192063A (zh) * | 2023-11-06 | 2023-12-08 | 山东大学 | 基于耦合卡尔曼滤波数据同化的水质预测方法及系统 |
| CN117192063B (zh) * | 2023-11-06 | 2024-03-15 | 山东大学 | 基于耦合卡尔曼滤波数据同化的水质预测方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111351488B (zh) | 飞行器智能轨迹重构再入制导方法 | |
| CN108803321B (zh) | 基于深度强化学习的自主水下航行器轨迹跟踪控制方法 | |
| CN106600059B (zh) | 基于改进rbf神经网络的智能电网短期负荷预测方法 | |
| CN114819068A (zh) | 一种混合型目标航迹预测方法及系统 | |
| CN114048889B (zh) | 基于长短期记忆网络的飞行器轨迹预测的方法 | |
| CN111695299B (zh) | 一种中尺度涡轨迹预测方法 | |
| CN108563119B (zh) | 一种基于模糊支持向量机算法的无人艇运动控制方法 | |
| CN112859898B (zh) | 一种基于双通道双向神经网络的飞行器轨迹预测方法 | |
| CN111913175A (zh) | 一种传感器短暂失效下带补偿机制的水面目标跟踪方法 | |
| CN114199248B (zh) | 一种基于混合元启发算法优化anfis的auv协同定位方法 | |
| CN111027627A (zh) | 一种基于多层感知机的振动信息地形分类识别方法 | |
| CN114357872A (zh) | 一种基于stacking模型融合的船舶运动黑箱辨识建模与运动预测方法 | |
| Fan et al. | Towards verification-aware knowledge distillation for neural-network controlled systems | |
| CN115033022A (zh) | 面向移动平台基于专家经验的ddpg无人机降落方法 | |
| CN116700327A (zh) | 一种基于连续动作优势函数学习的无人机轨迹规划方法 | |
| CN114167295B (zh) | 基于多算法融合的锂离子电池soc估算方法与系统 | |
| CN115265532A (zh) | 一种用于船用组合导航中的辅助滤波方法 | |
| CN116520281B (zh) | 一种基于ddpg的扩展目标跟踪优化方法和装置 | |
| CN113341696A (zh) | 一种运载火箭姿态控制参数智能整定方法 | |
| CN116432539A (zh) | 一种时间一致性协同制导方法、系统、设备及介质 | |
| CN115222023A (zh) | 基于神经网络的非合作目标控制策略识别方法及系统 | |
| Zhu et al. | An adaptive interactive multiple-model algorithm based on end-to-end learning | |
| Poddar et al. | Tuning of GPS aided attitude estimation using evolutionary algorithms | |
| CN115562345A (zh) | 一种基于深度强化学习的无人机侦测轨迹规划方法 | |
| CN116449360A (zh) | 一种基于长短时间记忆网络的机动目标跟踪方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220729 |
|
| RJ01 | Rejection of invention patent application after publication |