CN109855639B - 基于障碍物预测与mpc算法的无人驾驶轨迹规划方法 - Google Patents

基于障碍物预测与mpc算法的无人驾驶轨迹规划方法 Download PDF

Info

Publication number
CN109855639B
CN109855639B CN201910034446.9A CN201910034446A CN109855639B CN 109855639 B CN109855639 B CN 109855639B CN 201910034446 A CN201910034446 A CN 201910034446A CN 109855639 B CN109855639 B CN 109855639B
Authority
CN
China
Prior art keywords
lane
longitudinal
vehicle
time
speed
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.)
Active
Application number
CN201910034446.9A
Other languages
English (en)
Other versions
CN109855639A (zh
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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201910034446.9A priority Critical patent/CN109855639B/zh
Publication of CN109855639A publication Critical patent/CN109855639A/zh
Application granted granted Critical
Publication of CN109855639B publication Critical patent/CN109855639B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种基于障碍物预测与MPC算法的无人驾驶轨迹规划方法:周围车辆运动轨迹预测;碰撞预测分析与最小安全距离的计算;考轨迹与约束条件的生成;利用模型预测控制(MPC)进行轨迹优化,得到最终轨迹输出。本发明能够利用障碍物运动预测与模型预测控制(MPC)原理的轨迹最优化进行无人驾驶车辆的运动决策与轨迹规划。

Description

基于障碍物预测与MPC算法的无人驾驶轨迹规划方法
技术领域
本发明属于无人驾驶决策规划领域,更具体的说,是涉及一种基于障碍物预测与MPC算法的无人驾驶轨迹规划方法。
背景技术
随着社会的发展,人们对生活品质的要求越来越高,汽车成为人类生活不可或缺的交通工具,然而汽车数量的增加以及驾驶员自身安全意识的薄弱导致全球交通事故发生率急剧增长。面对日益严峻的交通安全与交通拥堵问题,构建智能交通***的任务变得更加急切。无人驾驶车辆作为构建智能交通***的关键部分,近年来受到各方的关注。无人驾驶车辆融合了环境感知与定位、决策规划与运动控制等多项功能,从而代替驾驶员的眼、脑和手,具有反应迅速、行驶安全可靠等优点。目前美国、英国和德国等一些国家的无人驾驶技术已发展较为成熟,我国的无人驾驶技术起步较晚,在一些关键技术的发展方面距离世界先进水平还有一定的差距。
模型预测控制(MPC)是20世纪70年代产生于工业过程控制领域的一类计算机控制算法,其基本要素为:预测模型、滚动优化、反馈校正和参考轨迹。MPC通过不断滚动的局部优化,使模型失配、畸变、干扰等引起的不确定性及时得到弥补,从而取得较好的动态控制性能。MPC最明显的优点是对模型的精确度要求不高并能在控制过程中增加多种约束,因此非常适合用在无人驾驶的规划和控制中。在近几年无人驾驶领域的发展中,MPC得到了很好的应用,例如加拿大滑铁卢大学采用多约束模型预测控制进行车辆避障与轨迹跟踪,瑞典沃尔沃汽车集团采用模型预测控制进行换道策略规划与车辆控制。在无人驾驶领域,常采用的预测模型有物理运动学规律、车辆运动学模型和车辆动力学模型等,由于运动规划对于车辆本身的特性要求不高,在本发明中采用物理运动学规律作为预测模型。
之前出现的一些基于MPC的规划算法只考虑了应用于直路的避障换道策略,或仅为给定路径的跟踪算法,本发明在根据周围环境做出轨迹预测后,除了考虑避障换道部分的轨迹规划外,还考虑路口转弯部分的轨迹规划。
发明内容
本发明的目的是为了克服现有技术中的不足,提供一种基于障碍物预测与MPC算法的无人驾驶轨迹规划方法,能够利用障碍物运动预测与模型预测控制(MPC)原理的轨迹最优化进行无人驾驶车辆的运动决策与轨迹规划。
本发明的目的是通过以下技术方案实现的。
本发明的基于障碍物预测与MPC算法的无人驾驶轨迹规划方法,包括以下步骤:
第一步:周围车辆运动轨迹预测
周围车辆包括同车道前方车辆、相邻车道前方和后方的车辆,通过车载传感器设备及相应的图像识别算法得到周围车辆的运动速度与方向;假设在预测时域内,周围车辆沿当前车道按照当前速度行驶,得到周围车辆运动轨迹;
第二步:碰撞预测分析与最小安全距离的计算
无人驾驶车辆VA,当前运动速度为vA,由决策层输出的参考速度vdes和运动指令可以估计车辆的大概行驶轨迹;为避免无人驾驶车辆VA与周围车辆发生碰撞,进行碰撞预测分析与最小安全距离(MSS)的计算;换道场景坐标系为建立在无人驾驶车辆后轴中心的车载坐标系,换道场景中包括同车道前方车辆V1,速度为v1,与VA的距离为L1;目标车道前方同方向行驶的车辆V2,速度为v2,与VA的距离为L2;目标车道后方同方向行驶的车辆V3,速度为v3,与VA的距离为L3。VA换道过程从时刻t0=0开始准备换道,到时刻tadj为施加横向加速度前的调整时间,tlat为施加横向加速度的时间,故时刻tadj+tlat施加横向加速度结束,VA完全到达目标车道,tfinish为换道策略预测时间;
①VA与V1碰撞分析:
对VA施加横向加速度,其右前角通过V1左侧纵向切线之前易于V1发生碰撞,两车之间的最小安全距离为:
Figure GDA0003573619380000031
其中,aA和a1分别为VA与V1的加速度,
Figure GDA0003573619380000032
为VA施加横向加速度开始到其右前角通过V1左侧纵向切线的时间;
②VA与V2碰撞分析:
对VA施加横向加速度,其左前角通过V2右侧纵向切线之后易于V2发生碰撞,两车之间的最小安全距离为:
Figure GDA0003573619380000033
其中,a2为V2的加速度,
Figure GDA0003573619380000034
为VA施加横向加速度开始到其左前角通过V2右侧纵向切线的时间;
③VA与V3碰撞分析:
对VA施加横向加速度,其左后角通过V3右侧纵向切线之后易于V3发生碰撞,两车之间的最小安全距离为:
Figure GDA0003573619380000035
其中,a3为V3的加速度,
Figure GDA0003573619380000036
为VA施加横向加速度开始到左后角通过V3右侧纵向切线的时间;
第三步:参考轨迹与约束条件的生成
采用换道场景坐标系,分为纵向和横向两部分,约束条件分别为位置、速度、加速度、加加速度的限制,其中加速度与加加速度的约束为确定范围,根据车辆的一般特性确定;考虑直路避障与十字路口转弯两种情况:
①直路避障:
在本发明中均设定为纵向匀速换道,一般换道过程航向角φ为3~5°,故得到施加横向加速度的时间为:
Figure GDA0003573619380000037
其中,l为车道宽度;
在时刻tadj前参考轨迹为VA沿当前车道的匀速和匀减速过程,对应横向位置约束为当前车道左右两条车道线;时间段
Figure GDA0003573619380000041
部分的换道轨迹用一段与纵向夹角为3~5°的直线代替,对应横向位置约束为当前车道和目标车道两侧的车道线;在时刻tadj+tlat之后为沿目标车道的匀加速和匀速过程,对应横向位置约束为目标车道的左右两条车道线;向左换道的参考轨迹与横向位置约束中包括横向位置最小值ymin和横向位置最大值ymax;由决策层输出参考速度vdes,产生换道行为一般是在当前车道前方车辆速度较小时,故在时刻tadj前VA的纵向速度约束为[v1,vdes];到达目标车道后VA为保持与V2和V3的安全距离,故纵向速度最大值为vxmax=min{vdes,v2},最小值为vxmin=max{v1,v3};VA的横向速度约束由tlat决定;
②十字路口转弯:
车辆需要在低速情况下进行转弯较安全,故在转弯前需要进行一段减速过程,减速至安全转弯速度vturn,故参考轨迹即为沿当前车道中心线的匀速和匀减速过程,对应的横向位置约束为当前车道的左右两侧车道线;转弯时采用一段圆弧,可根据转弯所到达的车道初始点确定转弯半径和圆弧圆心;转弯结束后,车辆先进行匀加速过程,直至速度达到vdes,对应位置约束仍为车道两侧车道线。由于VA整体的速度大小范围为[vturn,vdes],根据各时刻车辆运动方向确定其纵横向速度方向;十字路口右转的参考轨迹与位置约束中包括转弯前横向位置最小值ymin和横向位置最大值ymax,由于转弯后车辆沿横向行驶,还包括转弯后纵向位置最小值xmin、纵向位置最大值xmax以及转弯部分位置约束;
第四步:利用模型预测控制(MPC)进行轨迹优化,得到最终轨迹输出
本发明采用的模型为物理运动学规律为:
Figure GDA0003573619380000042
其中,x0和v0分别表示初始时刻车辆的位置和速度,x和v分别表示预测时刻t的车辆位置和速度,ax为相应加速度;利用MPC算法进行优化时,纵向和横向分别进行计算,采用相同的模型,下面以纵向为例介绍计算过程:
①选取状态量为
Figure GDA0003573619380000051
xx、vx、ax分别为纵向位置、纵向速度和纵向加速度,控制量为u(k)=Δax(k),也即Δax为纵向加加速度,故得到新的状态空间表达式为:
Figure GDA0003573619380000052
相应的
Figure GDA0003573619380000053
②由状态方程得到未来时刻的输出为:
Figure GDA0003573619380000054
③选取代价函数为:
Figure GDA0003573619380000055
其中,
Figure GDA0003573619380000056
I为单位矩阵,Np为预测时域,Nc为控制时域,xref为参考轨迹转换到当前车载坐标系下的纵向坐标,a和b均为正的权重系数。约束条件为:
Figure GDA0003573619380000061
其中,xmax和xmin分别为纵向位置的约束,vxmax和vxmin分别为纵向速度的约束,axmax和axmin分别为纵向加速度的约束,Δaxmax和Δaxmin分别为纵向加加速度的约束;
④通过二次规划求解得到最优控制量为U*=Δa(t|t);
⑤将得到的最优控制量作为下一步预测的输入,重复①-④的过程,直到得到t0~tfinish所有时刻的轨迹输出。
设定VA进行纵向匀速换道,故第二步中得到的三个最小安全距离为
sMSS(VA,V1)=(vA-v1)(tadj+tc1)
Figure GDA0003573619380000062
Figure GDA0003573619380000063
先计算vA=vdes是否满足换道条件,如不满足则VA需要提前进行减速至vA=v1后按t=t0重新进行计算。
与现有技术相比,本发明的技术方案所带来的有益效果是:
本发明采用MPC算法进行轨迹规划,对模型的精确度要求不高,对于车辆运动约束可以较明确地体现在MPC求解的约束条件中。本发明在进行车辆的轨迹规划时考虑了周围障碍物的运动轨迹影响,保证了轨迹的安全性。由MPC算法得到的轨迹输出满足了平滑性的要求。本发明将车辆运动分为纵向和横向分别进行规划,采用不同的预测时域,由于纵向运动主要为速度规划,故采用较大的预测时域使得速度变化较平缓,而横向运动主要为位置规划,主要体现在换道部分,故采用较小的预测时域,使得换道动作反应灵敏。在每一步的MPC预测求解时均利用上一步预测得到的位置点与速度建立新的车载坐标系,由此方便约束条件的处理,并且得到的纵横向速度和加速度更为直观。除换道策略外,本发明考虑了十字路口转弯部分的轨迹规划,得到较好的规划结果。本发明将轨迹规划与跟踪合为一体,简化了无人驾驶车辆规划层的工作。
附图说明
图1为换道场景示意图。
图2为车辆换道过程示意图。
图3为向左换道的参考轨迹与位置约束示意图。
图4为十字路口右转的参考轨迹与位置约束示意图。
图5为MPC求解过程中坐标系转换示意图。
图6为向左换道仿真结果图。
图7为十字路口右转仿真结果图。
具体实施方式
下面结合附图对本发明作进一步的描述。
无人驾驶车辆的决策规划层的三个基本要求为:安全性、实时性、舒适性。面对实时交通环境,为满足以上三点要求,本发明做出以下几个方面的研究:(1)周围运动车辆的轨迹预测;(2)换道策略:碰撞预测分析与最小安全距离的计算;(3)十字路口车辆转弯的轨迹规划;(4)模型预测控制算法参数设置。
本发明的基于障碍物预测与MPC算法的无人驾驶轨迹规划方法,包括以下步骤:
第一步:周围车辆运动轨迹预测
周围车辆包括同车道前方车辆、相邻车道前方和后方的车辆,通过车载传感器设备及相应的图像识别算法可以得到周围车辆的运动速度与方向。假设在预测时域内,周围车辆沿当前车道按照当前速度行驶,得到周围车辆运动轨迹。
第二步:碰撞预测分析与最小安全距离的计算
无人驾驶车辆VA,当前运动速度为vA,由决策层输出的参考速度vdes和运动指令可以估计车辆的大概行驶轨迹。为避免无人驾驶车辆VA与周围车辆发生碰撞,需要进行碰撞预测分析与最小安全距离(MSS)的计算。换道场景如图1所示,坐标系为建立在无人驾驶车辆后轴中心的车载坐标系。换道场景中包括同车道前方车辆V1,速度为v1,与VA的距离为L1;目标车道前方同方向行驶的车辆V2,速度为v2,与VA的距离为L2;目标车道后方同方向行驶的车辆V3,速度为v3,与VA的距离为L3。VA换道过程如图2所示,从时刻t0=0开始准备换道,到时刻tadj为施加横向加速度前的调整时间,tlat为施加横向加速度的时间,故时刻tadj+tlat施加横向加速度结束,VA完全到达目标车道,tfinish为换道策略预测时间。
①VA与V1碰撞分析:
对VA施加横向加速度,其右前角通过V1左侧纵向切线之前易于V1发生碰撞,两车之间的最小安全距离为:
Figure GDA0003573619380000081
其中,aA和a1分别为VA与V1的加速度,
Figure GDA0003573619380000082
为VA施加横向加速度开始到其右前角通过V1左侧纵向切线的时间。
②VA与V2碰撞分析:
对VA施加横向加速度,其左前角通过V2右侧纵向切线之后易于V2发生碰撞,两车之间的最小安全距离为:
Figure GDA0003573619380000083
其中,a2为V2的加速度,
Figure GDA0003573619380000084
为VA施加横向加速度开始到其左前角通过V2右侧纵向切线的时间。
③VA与V3碰撞分析:
对VA施加横向加速度,其左后角通过V3右侧纵向切线之后易于V3发生碰撞,两车之间的最小安全距离为:
Figure GDA0003573619380000091
其中,a3为V3的加速度,
Figure GDA0003573619380000095
为VA施加横向加速度开始到左后角通过V3右侧纵向切线的时间。
本发明中设定VA进行纵向匀速换道,故可以得到的三个最小安全距离为:
Figure GDA0003573619380000092
先计算vA=vdes是否满足换道条件,如不满足则VA需要提前进行减速至vA=v1后按t=t0重新进行计算。
第三步:参考轨迹与约束条件的生成
此部分坐标系为图1所示的采用换道场景坐标系,分为纵向和横向两部分,约束条件分别为位置、速度、加速度、加加速度的限制,其中加速度与加加速度的约束为确定范围,根据车辆的一般特性确定。在本发明中考虑直路避障与十字路口转弯两种情况,以下将分别进行说明。
①直路避障:
在本发明中均设定为纵向匀速换道,一般换道过程航向角φ为3~5°,故可以得到施加横向加速度的时间为:
Figure GDA0003573619380000093
其中,l为车道宽度。
在时刻tadj前参考轨迹为VA沿当前车道的匀速和匀减速过程,对应横向位置约束为当前车道左右两条车道线。时间段
Figure GDA0003573619380000094
部分的换道轨迹用一段与纵向夹角为3~5°的直线代替,对应横向位置约束为当前车道和目标车道两侧的车道线。在时刻tadj+tlat之后为沿目标车道的匀加速和匀速过程,对应横向位置约束为目标车道的左右两条车道线。向左换道的参考轨迹与横向位置约束如图3所示,其中虚线为横向位置最小值ymin,点划线为横向位置最大值ymax。由决策层输出参考速度vdes,产生换道行为一般是在当前车道前方车辆速度较小时,故在时刻tadj前VA的纵向速度约束为[v1,vdes]。到达目标车道后VA为保持与V2和V3的安全距离,故纵向速度最大值为vxmax=min{vdes,v2},最小值为vxmin=max{v1,v3}。VA的横向速度约束由tlat决定。
②十字路口转弯:
车辆需要在低速情况下进行转弯较安全,故在转弯前需要进行一段减速过程,减速至安全转弯速度vturn,故参考轨迹即为沿当前车道中心线的匀速和匀减速过程,对应的横向位置约束为当前车道的左右两侧车道线。转弯时采用一段圆弧,可根据转弯所到达的车道初始点确定转弯半径和圆弧圆心。转弯结束后,车辆先进行匀加速过程,直至速度达到vdes,对应位置约束仍为车道两侧车道线。由于VA整体的速度大小范围为[vturn,vdes],根据各时刻车辆运动方向确定其纵横向速度方向。如图4所示为十字路口右转的参考轨迹与位置约束,其中水平虚线为转弯前横向位置最小值ymin,水平点划线为横向位置最大值ymax,由于转弯后车辆沿横向行驶,竖直虚线为转弯后纵向位置最小值xmin,竖直点划线为纵向位置最大值xmax,转弯部分位置约束如虚线所示。
第四步:利用模型预测控制(MPC)进行轨迹优化,得到最终轨迹输出
本发明采用的模型为物理运动学规律为:
Figure GDA0003573619380000101
其中,x0和v0分别表示初始时刻车辆的位置和速度,x和v分别表示预测时刻t的车辆位置和速度,ax为相应加速度。由于换道过程中航向角较小,横向与纵向间速度相互影响不计。利用MPC算法进行优化时,纵向和横向分别进行计算,采用相同的模型,下面以纵向为例介绍计算过程。
①选取状态量为
Figure GDA0003573619380000111
xx、vx、ax分别为纵向位置、纵向速度和纵向加速度,控制量为u(k)=Δax(k),也即Δax为纵向加加速度,故得到新的状态空间表达式为:
Figure GDA0003573619380000112
相应的
Figure GDA0003573619380000113
②由状态方程可以得到未来时刻的输出为:
Figure GDA0003573619380000114
③选取代价函数为:
Figure GDA0003573619380000115
其中,
Figure GDA0003573619380000116
I为单位矩阵,Np为预测时域,Nc为控制时域,xref为参考轨迹转换到当前车载坐标系下的纵向坐标,a和b均为正的权重系数。约束条件由第三步得到为:
Figure GDA0003573619380000121
其中,xmax和xmin分别为纵向位置的约束,vxmax和vxmin分别为纵向速度的约束,axmax和axmin分别为纵向加速度的约束,Δaxmax和Δaxmin分别为纵向加加速度的约束。
④通过二次规划求解得到最优控制量为U*=Δa(t|t)。
求解时令
Figure GDA0003573619380000122
故代价函数可以整理为:
Figure GDA0003573619380000123
接下来根据二次规划的思想进行求解即可。
⑤将得到的最优控制量作为下一步预测的输入,重复①-④的过程,直到得到t0~tfinish所有时刻的轨迹输出。在每一步MPC求解时均需转换为当前位置的车载坐标系,如图5所示。原坐标系中任意点(xbefore,ybefore)与转换到新坐标系中的坐标值(xafter,yafter)的关系为:
Figure GDA0003573619380000131
其中(x1,y1)为上一时刻计算得到的位置点,也即转换坐标系后的原点,θ为其航向角。每一时刻得到的位置点都需转换到原始坐标系中进行输出,对于速度、加速度和加加速度无需进行转换。
在纵向和横向分别进行规划时采用不同的预测时域,由于纵向运动主要为速度规划,故采用较大的预测时域使得速度变化较平缓,而横向运动主要为位置规划,主要体现在换道部分,故采用较小的预测时域,使得换道动作反应灵敏。
多车道向左换道场景和十字路口转弯场景仿真结果如图6和图7所示。其中点状线为规定好的车道线,虚线为参考轨迹,实线为MPC优化后的轨迹。
尽管上面结合附图对本发明的功能及工作过程进行了描述,但本发明并不局限于上述的具体功能和工作过程,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以做出很多形式,这些均属于本发明的保护之内。

Claims (2)

1.一种基于障碍物预测与MPC算法的无人驾驶轨迹规划方法,其特征在于,包括以下步骤:
第一步:周围车辆运动轨迹预测
周围车辆包括同车道前方车辆、相邻车道前方和后方的车辆,通过车载传感器设备及相应的图像识别算法得到周围车辆的运动速度与方向;假设在预测时域内,周围车辆沿当前车道按照当前速度行驶,得到周围车辆运动轨迹;
第二步:碰撞预测分析与最小安全距离的计算
无人驾驶车辆VA,当前运动速度为vA,由决策层输出的参考速度vdes和运动指令估计车辆的大概行驶轨迹;为避免无人驾驶车辆VA与周围车辆发生碰撞,进行碰撞预测分析与最小安全距离(MSS)的计算;换道场景坐标系为建立在无人驾驶车辆后轴中心的车载坐标系,换道场景中包括同车道前方车辆V1,速度为v1,与VA的距离为L1;目标车道前方同方向行驶的车辆V2,速度为v2,与VA的距离为L2;目标车道后方同方向行驶的车辆V3,速度为v3,与VA的距离为L3;VA换道过程从时刻t0=0开始准备换道,到时刻tadj为施加横向加速度前的调整时间,tlat为施加横向加速度的时间,故时刻tadj+tlat施加横向加速度结束,VA完全到达目标车道,tfinish为换道策略预测时间;
①VA与V1碰撞分析:
对VA施加横向加速度,其右前角通过V1左侧纵向切线之前易与V1发生碰撞,两车之间的最小安全距离为:
Figure FDA0003583862420000011
其中,aA和a1分别为VA与V1的加速度,
Figure FDA0003583862420000012
为VA施加横向加速度开始到其右前角通过V1左侧纵向切线的时间;
②VA与V2碰撞分析:
对VA施加横向加速度,其左前角通过V2右侧纵向切线之后易与V2发生碰撞,两车之间的最小安全距离为:
Figure FDA0003583862420000021
其中,a2为V2的加速度,
Figure FDA0003583862420000022
为VA施加横向加速度开始到其左前角通过V2右侧纵向切线的时间;
③VA与V3碰撞分析:
对VA施加横向加速度,其左后角通过V3右侧纵向切线之后易与V3发生碰撞,两车之间的最小安全距离为:
Figure FDA0003583862420000023
其中,a3为V3的加速度,
Figure FDA0003583862420000024
为VA施加横向加速度开始到左后角通过V3右侧纵向切线的时间;
第三步:参考轨迹与约束条件的生成
采用换道场景坐标系,分为纵向和横向两部分,约束条件分别为位置、速度、加速度、加加速度的限制,其中加速度与加加速度的约束为确定范围,根据车辆的特性确定;考虑直路避障与十字路口转弯两种情况:
①直路避障:
设定VA进行纵向匀速换道,换道过程航向角φ为3~5°,得到施加横向加速度的时间为:
Figure DEST_PATH_FDA0003573619370000025
其中,l为车道宽度;
在时刻tadj前参考轨迹为VA沿当前车道的匀速和匀减速过程,对应横向位置约束为当前车道左右两条车道线;时间段[tadj,tadj+tlat]部分的换道轨迹用一段与纵向夹角为3~5°的直线代替,对应横向位置约束为当前车道和目标车道两侧的车道线;在时刻tadj+tlat之后为沿目标车道的匀加速和匀速过程,对应横向位置约束为目标车道的左右两条车道线;向左换道的参考轨迹对应的横向位置约束中包括横向位置最小值ymin和横向位置最大值ymax;由决策层输出参考速度vdes,换道行为是在当前车道前方车辆速度较小时产生,故在时刻tadj前VA的纵向速度约束为[v1,vdes];到达目标车道后,为保持VA与V2、V3的安全距离,纵向速度最大值为vxmax=min{vdes,v2},最小值为vxmin=max{v1,v3};VA的横向速度约束由tlat决定;
②十字路口转弯:
车辆在低速情况下进行转弯较安全,故在转弯前进行一段减速过程,减速至安全转弯速度vturn,参考轨迹即为沿当前车道中心线的匀速和匀减速过程,对应的横向位置约束为当前车道的左右两侧车道线;转弯时采用一段圆弧,根据转弯所到达的车道初始点确定转弯半径和圆弧圆心;转弯结束后,车辆先进行匀加速过程,直至速度达到vdes,对应位置约束仍为车道两侧车道线;由于VA整体的速度大小范围为[vturn,vdes],根据各时刻车辆运动方向确定其纵横向速度方向;十字路口右转的参考轨迹对应的位置约束中包括转弯前横向位置最小值ymin和横向位置最大值ymax,由于转弯后车辆沿横向行驶,还包括转弯后纵向位置最小值xmin、纵向位置最大值xmax以及转弯部分位置约束;
第四步:利用模型预测控制(MPC)进行轨迹优化,得到最终轨迹输出;
其中,采用的模型为物理运动学规律:
Figure FDA0003583862420000031
其中,x0和v0分别表示初始时刻车辆的位置和速度,x和v分别表示预测时刻t的车辆位置和速度,ax为相应加速度;利用MPC算法进行优化时,纵向和横向分别进行计算,采用相同的模型,下面以纵向为例介绍计算过程:
①选取状态量为
Figure FDA0003583862420000041
xx、vx、ax分别为纵向位置、纵向速度和纵向加速度,控制量为u(k)=Δax(k),也即Δax为纵向加加速度,故得到新的状态空间表达式为:
Figure FDA0003583862420000042
相应的
Figure FDA0003583862420000043
②由状态方程得到未来时刻的输出为:
Figure FDA0003583862420000044
③选取代价函数为:
Figure DEST_PATH_FDA0003573619370000044
其中,
Figure FDA0003583862420000046
I为单位矩阵,Np为预测时域,Nc为控制时域,xref为参考轨迹转换到当前车载坐标系下的纵向坐标,a和b均为正的权重系数,约束条件为:
Figure DEST_PATH_FDA0003573619370000051
其中,xmax和xmin分别为纵向位置的约束,vxmax和vxmin分别为纵向速度的约束,axmax和axmin分别为纵向加速度的约束,Δaxmax和Δaxmin分别为纵向加加速度的约束;
④通过二次规划求解得到最优控制量为U*=Δa(t|t);
⑤将得到的最优控制量作为下一步预测的输入,重复①-④的过程,直到得到t0~tfinish所有时刻的轨迹输出。
2.根据权利要求1所述的基于障碍物预测与MPC算法的无人驾驶轨迹规划方法,其特征在于,设定VA进行纵向匀速换道,故第二步中得到的三个最小安全距离为
Figure FDA0003583862420000052
Figure FDA0003583862420000053
Figure FDA0003583862420000054
设定vA=vdes,判断vA是否满足换道条件,如不满足则VA需要提前进行减速至vA=v1后按t=t0重新进行计算。
CN201910034446.9A 2019-01-15 2019-01-15 基于障碍物预测与mpc算法的无人驾驶轨迹规划方法 Active CN109855639B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910034446.9A CN109855639B (zh) 2019-01-15 2019-01-15 基于障碍物预测与mpc算法的无人驾驶轨迹规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910034446.9A CN109855639B (zh) 2019-01-15 2019-01-15 基于障碍物预测与mpc算法的无人驾驶轨迹规划方法

Publications (2)

Publication Number Publication Date
CN109855639A CN109855639A (zh) 2019-06-07
CN109855639B true CN109855639B (zh) 2022-05-27

Family

ID=66894626

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910034446.9A Active CN109855639B (zh) 2019-01-15 2019-01-15 基于障碍物预测与mpc算法的无人驾驶轨迹规划方法

Country Status (1)

Country Link
CN (1) CN109855639B (zh)

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110539752B (zh) * 2019-06-26 2020-12-18 江苏大学 一种智能汽车多预测范围模型预测轨迹跟踪控制方法及***
CN110329263B (zh) * 2019-07-05 2020-08-21 西南交通大学 自动驾驶车辆自适应换道轨迹规划方法
CN110298131B (zh) * 2019-07-05 2021-07-13 西南交通大学 一种混合驾驶环境下自动驾驶换道决策模型建立方法
CN110244746B (zh) * 2019-07-11 2020-02-18 肇庆学院 一种基于视觉注意的机器人动态障碍物避开方法及***
CN110488816B (zh) * 2019-08-06 2021-12-31 华为技术有限公司 自动驾驶纵向规划方法及相关设备
CN110647801A (zh) * 2019-08-06 2020-01-03 北京汽车集团有限公司 设置感兴趣区域的方法、装置、存储介质及电子设备
CN110362096B (zh) * 2019-08-13 2021-05-18 东北大学 一种基于局部最优性的无人驾驶车辆动态轨迹规划方法
CN110466514B (zh) * 2019-08-30 2020-10-27 北京小马慧行科技有限公司 车辆控制方法及装置
CN110466513A (zh) * 2019-08-30 2019-11-19 北京小马慧行科技有限公司 车辆控制方法及装置
CN110703754B (zh) * 2019-10-17 2021-07-09 南京航空航天大学 一种自动驾驶车辆路径与速度高度耦合的轨迹规划方法
CN111785062B (zh) * 2020-04-01 2021-09-14 北京京东乾石科技有限公司 在无信号灯路口实现车路协同的方法和装置
WO2021212379A1 (zh) * 2020-04-22 2021-10-28 华为技术有限公司 车道线检测方法及装置
CN111710159B (zh) * 2020-05-29 2021-09-03 同济大学 一种基于虚拟车道线的交叉路口车辆路径规划方法和装置
CN113805572B (zh) * 2020-05-29 2023-12-15 华为技术有限公司 运动规划的方法与装置
CN111932881A (zh) * 2020-08-12 2020-11-13 广西科技大学 交通路口管理方法、装置、终端设备及存储介质
CN114148344B (zh) * 2020-09-08 2023-06-02 华为技术有限公司 一种车辆行为预测方法、装置及车辆
CN112193244B (zh) * 2020-09-30 2021-07-20 浙江大学 基于线性约束的自动驾驶车辆运动规划方法
CN112269384B (zh) * 2020-10-23 2021-09-14 电子科技大学 一种结合障碍物行为意图的车辆动态轨迹规划方法
CN112249019A (zh) * 2020-10-30 2021-01-22 上海智驾汽车科技有限公司 一种基于纵向安全距离的车辆换道决策方法
CN112631276A (zh) * 2020-12-03 2021-04-09 长沙行深智能科技有限公司 一种无人车动态障碍决策方法、***、介质及设备
CN112622930A (zh) * 2020-12-22 2021-04-09 北京百度网讯科技有限公司 无人车的行驶控制方法、装置、设备以及自动驾驶车辆
CN114802245B (zh) * 2021-01-11 2024-04-05 广东科学技术职业学院 根据交通灯进行并道的方法
CN112965489B (zh) * 2021-02-05 2022-04-08 北京理工大学 一种基于碰撞检测的智能车高速换道规划方法
CN112874509B (zh) * 2021-03-03 2022-04-29 知行汽车科技(苏州)有限公司 基于智能驾驶员模型idm的轨迹规划方法、装置及存储介质
CN113655789A (zh) * 2021-08-04 2021-11-16 东风柳州汽车有限公司 路径跟踪方法、装置、车辆及存储介质
CN113867365B (zh) * 2021-10-28 2024-05-14 广州文远知行科技有限公司 无人驾驶车辆变加速度的确定方法、装置和相关设备
CN114084155B (zh) * 2021-11-15 2023-10-20 清华大学 预测型智能汽车决策控制方法、装置、车辆及存储介质
CN115123310B (zh) * 2022-08-31 2022-12-13 中汽数据(天津)有限公司 无人车避障局部路径规划方法、设备和存储介质

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8190330B2 (en) * 2009-03-06 2012-05-29 GM Global Technology Operations LLC Model based predictive control for automated lane centering/changing control systems
CN105857294B (zh) * 2016-05-04 2018-12-04 常州工学院 一种汽车换道避撞控制方法
EP3291202B1 (en) * 2016-08-29 2019-04-17 Volvo Car Corporation Method of road vehicle trajectory planning
CN106926844B (zh) * 2017-03-27 2018-10-19 西南交通大学 一种基于实时环境信息的动态自动驾驶换道轨迹规划方法
CN107168319B (zh) * 2017-06-01 2020-07-28 同济大学 一种基于模型预测控制的无人车避障方法
CN108387242B (zh) * 2018-02-07 2021-04-09 西南交通大学 自动驾驶换道准备和执行一体化轨迹规划方法

Also Published As

Publication number Publication date
CN109855639A (zh) 2019-06-07

Similar Documents

Publication Publication Date Title
CN109855639B (zh) 基于障碍物预测与mpc算法的无人驾驶轨迹规划方法
CN110597245B (zh) 基于二次型规划和神经网络的自动驾驶换道轨迹规划方法
CN110329263B (zh) 自动驾驶车辆自适应换道轨迹规划方法
CN110187639B (zh) 一种基于参数决策框架的轨迹规划控制方法
CN109035862B (zh) 一种基于车车通信的多车协同换道控制方法
CN111338340B (zh) 基于模型预测的无人驾驶汽车局部路径规划方法
JP6651642B2 (ja) 車両制御装置
WO2020143288A1 (zh) 一种复杂工况下自动驾驶车辆决策***及其轨迹规划方法
CN103754221B (zh) 一种车辆自适应巡航控制***
CN106218638B (zh) 一种智能网联汽车协同换道控制方法
CN110377039A (zh) 一种车辆避障轨迹规划与跟踪控制方法
CN110362096A (zh) 一种基于局部最优性的无人驾驶车辆动态轨迹规划方法
CN112965476B (zh) 一种基于多窗口模型的高速无人车轨迹规划***及方法
CN112233413B (zh) 一种面向智能网联车辆的多车道时空轨迹优化方法
CN103085816A (zh) 一种用于无人驾驶车辆的轨迹跟踪控制方法及控制装置
JP6630443B2 (ja) 車両制御装置
WO2019003302A1 (ja) 車両制御装置
CN111899509B (zh) 一种基于车路信息耦合的智能网联汽车状态向量计算方法
Guo et al. Anticipative and predictive control of automated vehicles in communication-constrained connected mixed traffic
CN114194215A (zh) 一种智能车辆避障换道轨迹规划方法及***
CN117141489B (zh) 一种基于最小作用量原理的智能车辆轨迹分层规划方法
Schitz et al. Model-predictive cruise control for direct teleoperated driving tasks
CN114442630A (zh) 一种基于强化学习和模型预测的智能车规划控制方法
CN115042770B (zh) 一种基于分布式鲁棒模型预测的车辆队列横向控制方法
CN116494976A (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
GR01 Patent grant
GR01 Patent grant