CN115762182A - 一种基于运动学约束的车辆智能轨迹预测方法 - Google Patents

一种基于运动学约束的车辆智能轨迹预测方法 Download PDF

Info

Publication number
CN115762182A
CN115762182A CN202211368939.4A CN202211368939A CN115762182A CN 115762182 A CN115762182 A CN 115762182A CN 202211368939 A CN202211368939 A CN 202211368939A CN 115762182 A CN115762182 A CN 115762182A
Authority
CN
China
Prior art keywords
vehicle
model
acceleration
future
term memory
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
Application number
CN202211368939.4A
Other languages
English (en)
Inventor
邬海杰
张旭
蒋锟
周泰
马莹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Tage Idriver Technology Co Ltd
Original Assignee
Beijing Tage Idriver Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Tage Idriver Technology Co Ltd filed Critical Beijing Tage Idriver Technology Co Ltd
Priority to CN202211368939.4A priority Critical patent/CN115762182A/zh
Publication of CN115762182A publication Critical patent/CN115762182A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于运动学约束的车辆智能轨迹预测方法,包括步骤:构建长短期记忆网络模型,根据车辆过去t1内的横向误差、航向误差、横摆角速度,预测未来t2内车辆的横摆角速度;基于当前实时速度与实时加速度,构建衰减加速度模型,在t2内加速度衰减到0,迭代求出未来t2内纵向速度,作为未来t2内车辆的纵向速度测量值;建立车辆运动学模型,将长短期记忆网络预测出的横摆角速度和衰减加速度模型输出的纵向速度作为中间变量,预测车辆的未来轨迹。该方法考虑预测算法的泛化性及容错性,在综合考虑车辆的可执行性的情形下,采用结合长短期记忆网络与车辆运动学模型的车辆轨迹预测方法,使预测出的轨迹在满足车辆可执行的条件下,更加准确。

Description

一种基于运动学约束的车辆智能轨迹预测方法
技术领域
本发明涉及自动驾驶车辆技术领域,尤其涉及自动驾驶车辆轨迹预测领域。具体讲,涉及一种基于运动学约束的车辆智能轨迹预测方法。
背景技术
现有的关于车辆轨迹预测的研究可以大致分为两类。分别是基于浅层学习的方法和基于深度学习的轨迹预测算法。基于浅层学习的方法主要是利用基于动力学(运动学)特征与贝叶斯滤波器相结合的方法,该模型能够预测出符合车辆实际的预测轨迹,但是该模型以恒定速度为假设,不符合实际情况。基于深度学习的方法主要利用大规模的数据集学习出合理的映射关系,其比较典型的深度学习网络有RNN、长短期记忆网络(LSTM)等。深度学习的方法不需要对模型进行复杂的建模,但是该方法会预测出一些车辆无法执行的轨迹。
基于上述两种方法的局限性,本发明提出了一种结合长短期记忆网络与车辆运动学模型的车辆轨迹预测方法,该方法首先利用长短期记忆网络对大量的行车数据进行训练,得到一个泛化性较高的网络模型,利用该模型预测车辆的横摆角速度。其次,构建衰减加速度模型预测未来的纵向车速,由于车辆在行车过程中是一个变加速运动,因此采用衰减加速度模型能够更符合车辆实际行驶过程。最后,利用预测出的横摆角速度与纵向速度,结合从车辆动力学模型生成符合车辆运动学的轨迹。
发明内容
为克服现有技术的不足,本发明旨在提出一种能够应用于车辆实际行驶过程中的轨迹预测算法。具体而言,考虑预测算法的泛化性及容错性,在综合考虑车辆的可执行性的情形下,采用结合长短期记忆网络与车辆运动学模型的车辆轨迹预测方法,利用长短期记忆网络的能够使预测的车辆横摆角速度更具有泛化性,同时利用车辆运动学模型对规划轨迹进行约束,从而使预测出的轨迹在满足车辆可执行的条件下,更加准确。
本发明采用的技术方案是:一种基于运动学约束的车辆智能轨迹预测方法,包括以下步骤,
第一步,构建长短期记忆网络模型,根据车辆过去t1内的横向误差、航向误差、横摆角速度,预测未来t2内车辆的横摆角速度;
第二步,基于当前实时速度与实时加速度,构建衰减加速度模型,在t2内加速度衰减到0,迭代求出未来t2内纵向速度,作为未来t2内车辆的纵向速度测量值;
第三步,建立车辆运动学模型,将长短期记忆网络预测出的横摆角速度和衰减加速度模型输出的纵向速度作为中间变量,预测车辆的未来轨迹。
进一步,所述第一步,所述长短期记忆网络模型的网络架构包括输入层、LSTM网络层、全链接层以及回归层;其中输入层为一层,输入维度为特征数目;LSTM网络层设置有250层;全链接层设置为1层,输出维度为1,表示输出特征数目。
进一步,所述第三步包括以下子步骤:
(1)建立车辆的状态矢量:
Figure BDA0003924215250000021
其中,xk,yk,θk,vk,γk,ak
Figure BDA0003924215250000022
分别为全局坐标x,全局坐标y,车辆航向,纵向速度,横摆角速度,纵向加速度,横摆角加速度;
(2)建立过程模型:基于车辆的运动学方程,建立(1)中状态矢量下一时刻与当前时刻关系的过程模型,
xk+1=f(xk)+qk=[f1,f2,f3,f4,f5,f6,f7]+qk
其中,f(xk)是状态矩阵;qk是过程噪声,qk~N(0,Qk),
Figure BDA0003924215250000023
非线性状态转移模型:
f1
Figure BDA0003924215250000031
f2
Figure BDA0003924215250000032
f3
Figure BDA0003924215250000033
f4:vk+1=vk+akT
f5
Figure BDA0003924215250000034
f6:ak+1=ak
f7
Figure BDA0003924215250000035
其中,采样周期T=0.05s;
(3)获取模拟测量值:将第一步与第二步中预测出的速度与横摆角速度作为模拟测量值,用于后续的EKF计算,
yk=[vk,γk]T
其中γk、vk分别为第一步中利用长短期记忆网络神经网络获得的横摆角速度和第二步中利用衰减加速度模型获得的纵向速度;
(4)利用非线性卡尔曼滤波对(1)中的车辆状态矢量xk进行预测、更新:
预测当前时刻先验状态矢量及协方差矩阵:
Figure BDA0003924215250000036
Figure BDA0003924215250000037
其中,
Figure BDA0003924215250000038
是当前时刻先验估计状态,xk-1是上一时刻后验估计状态,f(·)是状态转移方程,
Figure BDA0003924215250000039
是当前时刻的先验协方差矩阵,F是状态转移矩阵,Pk-1是上一时刻的后验协方差矩阵,Qk-1为上一时刻的过程噪声;
更新当前时刻状态矢量及协方差矩阵:
A.卡尔曼增益更新:
Figure BDA00039242152500000310
Figure BDA00039242152500000311
其中,Kk是当前时刻卡尔曼增益,
Figure BDA0003924215250000041
为观测矩阵,Rk是当前时刻观测噪声协方差矩阵,Rk=diag{0.01,0.001};
B.状态矢量更新:
Figure BDA0003924215250000042
Figure BDA0003924215250000043
C.协方差矩阵更新:
Figure BDA0003924215250000044
其中,Pk是当前时刻的后验协方差矩阵,I是单位矩阵。
本发明与现有技术相比所具有的有益效果:
(1)将长短期记忆网络作为车辆状态观测器,基于车辆历史轨迹预测未来横摆角速度。由于基于神经网络可以充分逼近任意复杂的非线性关系的特性,可以在无法建立准确的车辆预测模型的情况下准确预测出车辆未来的横摆角速度。
(2)建立衰减加速度模型,将预测车辆在未来5s内为匀速运动转为变加速运动,能够更加符合真实的纵向车速变化;
(3)建立车辆运动学模型,将第一步预测出的横摆角速度以及第二步经过衰减加速度模型输出的速度作为输入,通过对长短期记忆网络加入运动学约束来预测未来的车辆轨迹,解决了长短期记忆网络预测结果不符合车辆运动学轨迹的问题,提高了预测算法的适应性及鲁棒性;
附图说明
图1是结合长短期记忆网络与车辆运动学模型的车辆轨迹预测方法流程图。
图2是长短期记忆网络预测横摆角速度流程。
具体实施方式
本发明采用的技术方案是,利用长短期记忆网络与衰减加速度模型预测中间变量横摆角速度和纵向速度,将中间变量输入到车辆运动学模型中进行车辆轨迹的预测。一种基于运动学约束的车辆智能轨迹预测方法,如图1所示,具体步骤如下:
第一步,车辆横向参数预测——长短期记忆网络预测车辆横摆角速度:为了能够得到准确的车辆的横摆角速度,采用长短期记忆网络模型进行预测;长短期记忆网络通过大量的数据学习,再结合车辆的历史轨迹就能够预测出比较准确的横摆角速度。
建立长短期记忆网络模型,设置网络架构为输入层,LSTM网络层,全链接层以及回归层。其中输入层为一层,输入维度为特征数目;LSTM网络层设置有250层;全链接层设置为1层,输出维度为1,表示输出特征数目。
模型的输入为:横向误差、航向误差、横摆角速度,输出为:横摆角速度。即根据车辆过去3s内的横向误差、航向误差、横摆角速度,预测未来5s内的车辆横摆角速度。其中,横向误差为参考路径关键点相对于车辆实际位置的横向距离,取左正右负。航向误差为参考路径关键点相对于车辆实际位置的航向角偏差。长短期记忆网络预测横摆角速度流程如图2所示。
第二步,车辆纵向参数预测——衰减加速度模型预测车辆纵向车速:车辆在运动过程中车速并不是恒定不变的,而未来的加速度难以预测,因此采用衰减加速度模型,模拟车辆的变加速度运动。
车辆在运动过程中车速并不是恒定不变的,而未来的加速度难以预测,因此采用衰减加速度模型,模拟车辆的变加速度运动。输入参数为当前速度与当前加速度,输出参数为未来5s内的速度。
具体过程为:基于当前实时速度与实时加速度,构建衰减加速度模型,在5s内加速度衰减到0,迭代求出未来5s内纵向速度,作为未来有限时域(5s)内的纵向速度测量值。
第三步,车辆轨迹预测——结合车辆运动学模型预测车辆未来轨迹:单纯的基于数据模型的方法预测出的轨迹有时候不符合车辆的运动学特性。因此建立车辆运动学模型,将长短期记忆网络预测出的结果和衰减加速度模型输出的结果作为中间变量,在保证预测轨迹的可行性的同时,也能够防止网络预测失败导致的轨迹预测结果偏大。
基于EKF(扩展卡尔曼滤波)对预测时域内的车辆轨迹进行预测,由此将轨迹预测问题转化为状态估计问题。
(1)建立车辆的状态矢量:
Figure BDA0003924215250000061
其中,xk,yk,θk,vk,γk,ak
Figure BDA0003924215250000062
分别为全局坐标x,全局坐标y,车辆航向,纵向速度,横摆角速度,纵向加速度,横摆角加速度。
(2)建立过程模型:基于车辆的运动学方程,建立(1)中状态矢量下一时刻与当前时刻关系的过程模型,
xk+1=f(xk)+qk=[f1,f2,f3,f4,f5,f6,f7]+qk
其中,f(xk)是状态矩阵,采样周期T=0.05s;qk是过程噪声,qk~N(0,Qk),
Figure BDA0003924215250000063
非线性状态转移模型:
f1
Figure BDA0003924215250000064
f2
Figure BDA0003924215250000065
f3
Figure BDA0003924215250000066
f4:vk+1=vk+akT
f5
Figure BDA0003924215250000067
f6:ak+1=ak
f7
Figure BDA0003924215250000068
其中,采样周期T=0.05s。
(3)获取模拟测量值:将第一步与第二步中预测出的速度与横摆角速度作为模拟测量值,用于后续的EKF计算,
yk=[vk,γk]T
其中γk、vk分别为第一步中利用长短期记忆网络神经网络获得的横摆角速度和第二步中利用衰减加速度模型获得的纵向速度。
(4)利用非线性卡尔曼滤波对(1)中的车辆状态矢量xk进行预测、更新:
预测当前时刻先验状态矢量及协方差矩阵:
Figure BDA0003924215250000071
Figure BDA0003924215250000072
其中,
Figure BDA0003924215250000073
是当前时刻先验估计状态,xk-1是上一时刻后验估计状态,f(·)是状态转移方程,
Figure BDA0003924215250000074
是当前时刻的先验协方差矩阵,F是状态转移矩阵,Pk-1是上一时刻的后验协方差矩阵,Qk-1为上一时刻的过程噪声。
更新当前时刻状态矢量及协方差矩阵:
A.卡尔曼增益更新:
Figure BDA0003924215250000075
Figure BDA0003924215250000076
其中,Kk是当前时刻卡尔曼增益,
Figure BDA0003924215250000077
为观测矩阵,Rk是当前时刻观测噪声协方差矩阵,Rk=diag{0.01,0.001}。
B.状态矢量更新:
Figure BDA0003924215250000078
Figure BDA0003924215250000079
C.协方差矩阵更新:
Figure BDA00039242152500000710
其中,Pk是当前时刻的后验协方差矩阵,I是单位矩阵。
以上所述仅为本发明的具体实施方式,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于运动学约束的车辆智能轨迹预测方法,其特征在于,包括以下步骤:
第一步,构建长短期记忆网络模型,根据车辆过去t1内的横向误差、航向误差、横摆角速度,预测未来t2内车辆的横摆角速度;
第二步,基于当前实时速度与实时加速度,构建衰减加速度模型,在t2内加速度衰减到0,迭代求出未来t2内纵向速度,作为未来t2内车辆的纵向速度测量值;
第三步,建立车辆运动学模型,将长短期记忆网络预测出的横摆角速度和衰减加速度模型输出的纵向速度作为中间变量,预测车辆的未来轨迹。
2.根据权利要求1所述的一种基于运动学约束的车辆智能轨迹预测方法,其特征在于,所述第一步,所述长短期记忆网络模型的网络架构包括输入层、LSTM网络层、全链接层以及回归层;其中输入层为一层,输入维度为特征数目;LSTM网络层设置有250层;全链接层设置为1层,输出维度为1,表示输出特征数目。
3.根据权利要求1或2所述的一种基于运动学约束的车辆智能轨迹预测方法,其特征在于,所述第三步包括以下子步骤:
(1)建立车辆的状态矢量:
Figure FDA0003924215240000011
其中,xk,yk,θk,vk,γk,ak
Figure FDA0003924215240000012
分别为全局坐标x,全局坐标y,车辆航向,纵向速度,横摆角速度,纵向加速度,横摆角加速度;
(2)建立过程模型:基于车辆的运动学方程,建立(1)中状态矢量下一时刻与当前时刻关系的过程模型,
xk+1=f(xk)+qk=[f1,f2,f3,f4,f5,f6,f7]+qk
其中,f(xk)是状态矩阵;qk是过程噪声,qk~N(0,Qk),
Figure FDA0003924215240000013
非线性状态转移模型:
Figure FDA0003924215240000021
Figure FDA0003924215240000022
Figure FDA0003924215240000023
f4:vk+1=vk+akT
Figure FDA0003924215240000024
f6:ak+1=ak
Figure FDA0003924215240000025
其中,采样周期T=0.05s;
(3)获取模拟测量值:将第一步与第二步中预测出的速度与横摆角速度作为模拟测量值,用于后续的EKF计算,
yk=[vk,γk]T
其中γk、vk分别为第一步中利用长短期记忆网络神经网络获得的横摆角速度和第二步中利用衰减加速度模型获得的纵向速度;
(4)利用非线性卡尔曼滤波对(1)中的车辆状态矢量xk进行预测、更新:
预测当前时刻先验状态矢量及协方差矩阵:
Figure FDA0003924215240000026
Figure FDA0003924215240000027
其中,
Figure FDA0003924215240000028
是当前时刻先验估计状态,xk-1是上一时刻后验估计状态,f(·)是状态转移方程,
Figure FDA0003924215240000029
是当前时刻的先验协方差矩阵,F是状态转移矩阵,Pk-1是上一时刻的后验协方差矩阵,Qk-1为上一时刻的过程噪声;
更新当前时刻状态矢量及协方差矩阵:
A.卡尔曼增益更新:
Figure FDA00039242152400000210
Figure FDA00039242152400000211
其中,Kk是当前时刻卡尔曼增益,
Figure FDA0003924215240000031
为观测矩阵,Rk是当前时刻观测噪声协方差矩阵,Rk=diag{0.01,0.001};
B.状态矢量更新:
Figure FDA0003924215240000032
Figure FDA0003924215240000033
C.协方差矩阵更新:
Figure FDA0003924215240000034
其中,Pk是当前时刻的后验协方差矩阵,I是单位矩阵。
CN202211368939.4A 2022-11-03 2022-11-03 一种基于运动学约束的车辆智能轨迹预测方法 Pending CN115762182A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211368939.4A CN115762182A (zh) 2022-11-03 2022-11-03 一种基于运动学约束的车辆智能轨迹预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211368939.4A CN115762182A (zh) 2022-11-03 2022-11-03 一种基于运动学约束的车辆智能轨迹预测方法

Publications (1)

Publication Number Publication Date
CN115762182A true CN115762182A (zh) 2023-03-07

Family

ID=85357560

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211368939.4A Pending CN115762182A (zh) 2022-11-03 2022-11-03 一种基于运动学约束的车辆智能轨迹预测方法

Country Status (1)

Country Link
CN (1) CN115762182A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117932234A (zh) * 2024-03-25 2024-04-26 苏州观瑞汽车技术有限公司 一种用于制作刹车标定表的数据处理方法及***

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117932234A (zh) * 2024-03-25 2024-04-26 苏州观瑞汽车技术有限公司 一种用于制作刹车标定表的数据处理方法及***
CN117932234B (zh) * 2024-03-25 2024-06-07 苏州观瑞汽车技术有限公司 一种用于制作刹车标定表的数据处理方法及***

Similar Documents

Publication Publication Date Title
WO2021103834A1 (zh) 换道决策模型生成方法和无人车换道决策方法及装置
CN111624992B (zh) 一种基于神经网络的搬运机器人的路径跟踪控制方法
CN111752274B (zh) 一种基于强化学习的激光agv的路径跟踪控制方法
CN111428317B (zh) 一种基于5g和循环神经网络的关节摩擦力矩补偿方法
CN113671962B (zh) 一种无人驾驶铰接式清扫车的横向控制方法
Yang et al. Longitudinal tracking control of vehicle platooning using DDPG-based PID
CN114199248B (zh) 一种基于混合元启发算法优化anfis的auv协同定位方法
CN111159642B (zh) 一种基于粒子滤波的在线轨迹预测方法
CN115762182A (zh) 一种基于运动学约束的车辆智能轨迹预测方法
CN110716575A (zh) 基于深度双q网络强化学习的uuv实时避碰规划方法
CN114169230A (zh) 一种机器人动力学参数的辨识方法
CN113093548A (zh) 基于事件触发机制的移动机器人轨迹跟踪最优控制方法
WO2024087766A1 (zh) 一种考虑时变时延的无人车横纵向协同控制方法及装置
CN103399488B (zh) 基于自学习的多模型控制方法
CN112416021A (zh) 一种基于学习的旋翼无人机路径跟踪预测控制方法
CN110039537B (zh) 一种基于神经网络的在线自学习多关节运动规划方法
CN114217603B (zh) 一种多无人船编队安全优化控制***
CN113391553B (zh) 具有执行器饱和的异构cacc***的自适应最优控制方法
CN116819973B (zh) 一种轨迹跟踪控制方法
Liu et al. Robust Adaptive Self‐Structuring Neural Network Bounded Target Tracking Control of Underactuated Surface Vessels
Li et al. Performance analysis of deep learning supported Kalman filter
CN114012733B (zh) 一种用于pc构件模具划线的机械臂控制方法
CN114967472A (zh) 一种无人机轨迹跟踪状态补偿深度确定性策略梯度控制方法
Liu et al. Safe model-based control from signal temporal logic specifications using recurrent neural networks
Shi et al. AH‐Infinity Control for Path Tracking with Fuzzy Hyperbolic Tangent Model

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