CN111538328B - 一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法 - Google Patents

一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法 Download PDF

Info

Publication number
CN111538328B
CN111538328B CN202010258390.8A CN202010258390A CN111538328B CN 111538328 B CN111538328 B CN 111538328B CN 202010258390 A CN202010258390 A CN 202010258390A CN 111538328 B CN111538328 B CN 111538328B
Authority
CN
China
Prior art keywords
vehicle
control
autonomous driving
longitudinal
obstacle avoidance
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
CN202010258390.8A
Other languages
English (en)
Other versions
CN111538328A (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN202010258390.8A priority Critical patent/CN111538328B/zh
Publication of CN111538328A publication Critical patent/CN111538328A/zh
Application granted granted Critical
Publication of CN111538328B publication Critical patent/CN111538328B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0212Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
    • G05D1/0214Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory in accordance with safety or protection criteria, e.g. avoiding hazardous areas
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0212Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
    • G05D1/0223Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory involving speed control of the vehicle

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法,本方法基于离散时间车辆动力学模型,首先定义出描述行车环境的势场函数,再将离散时间动力学模型作为预测模型,然后建立有限时域优先级分层优化控制问题,并结合滚动优化控制原理,在每个采样时刻计算得到最优控制解,实现基于目标优先级驱动的自主车辆避障轨迹规划与跟踪控制的预测控制。本发明方法的优点是理解简单,通用性强,控制量满足自主驾驶车辆安全避障规划要求和车辆操纵稳定性。

Description

一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分 层预测控制方法
技术领域
本发明属于自主驾驶车辆轨迹规划与跟踪控制领域,涉及一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法。
背景技术
自主无人驾驶是智能车辆发展的重要方向,自主车辆的轨迹规划在环境感知的基础上,为车辆规划一条安全无碰撞且可行的行驶路径,是实现车辆自主无人驾驶的关键技术。自主驾驶车辆轨迹规划与跟踪控制应满足车辆动力学/运动学特性以及各种约束,预测控制算法能够很好地解决此类问题。通过对现有基于预测控制的自主驾驶车辆轨迹规划与跟踪控制的文献检索,发现现有方法将避障目标函数和跟踪目标函数直接相加,转换为单个综合目标函数,并用加权系数表示避障目标和跟踪目标的相对重要性,但避障目标函数和跟踪目标函数不同,度量的单位不同,而且避障目标和跟踪目标存在相互冲突特性,因此,将避障目标函数和跟踪目标函数直接相加的方法不能真实反映自主驾驶车辆的避障高优先级要求,尤其是对于车辆这种复杂非线性动力学系统而言,固定两者的加权系数适应性差,但行驶过程在线调整两者的加权系数将增加驾驶危险性而且难以实现。本方法把自主驾驶车辆避障轨迹规划与跟踪控制过程建立为一个目标优先级明确的多目标控制问题,将车辆的避障轨迹规划目标作为最高优先级,跟踪控制目标作为低优先级,然后从高优先级避障轨迹规划目标要求开始计算,逐层实现自主车辆避障轨迹与跟踪控制目标,保证对更高优先级避障目标的控制效果。本方法与现有自主驾驶车辆避障轨迹规划与跟踪控制方法相比,不需要将避障目标函数和跟踪目标函数相加,从而不存在两者的加权系数,因而本发明方法使用上更具灵活性和实用性。
发明内容
相对现有自主驾驶车辆避障轨迹规划与跟踪控制的加权预测控制方法,为了克服现有技术的不足,本发明设定自主驾驶车辆的避障轨迹规划目标为高优先级控制目标,自主驾驶车辆的跟踪为次优先级控制目标,根据排序结果,依次优化自主驾驶车辆的避障轨迹,同时确保自主驾驶车辆对优化避障轨迹中的安全跟踪控制,避免了避障轨迹规划目标和跟踪控制目标的加权系数复杂的在线调整和重新配置过程。
本发明解决其技术问题所采用的技术方案是:
一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法,所述方法包括如下步骤:
1)、考虑自主驾驶车辆离散时间动力学模型,参见式(1):
Figure BDA0002438340670000021
其中,k为采样时间,Ts为采样周期,X、Y、
Figure BDA0002438340670000022
分别为车辆在全局坐标系下的纵向位置、横向位置、横摆角,vx、vy、r分别为车辆在质心处的纵向速度、横向速度、横摆角速度,m为车辆质量,Iz为车辆绕其垂直轴的转动惯量,df和dr分别是车辆质心到前、后轴的距离,δf和ax分别为前轮转向角和纵向加速度,Ccf和Ccr分别表示前、后轮的侧偏刚度,定义车辆行驶环境势场函数,参见式(2):
Figure BDA0002438340670000023
其中,UA为自主驾驶车辆总势场,常数e为自然指数,Al为分道线势场强度,σl为分道线势场收敛幅度,yl,i为第i条道路分道线的横向位置,Ar为道路边界势场强度,yr,i为第i条道路边界的横向位置,Xgoal为目标点纵向位置,组合参数ε和β为正整数,D是自主驾驶车辆与障碍车辆之间的横向距离,σo为横向势场收敛幅度,Uo为障碍车辆纵向势场函数,参见式(3):
Figure BDA0002438340670000031
其中,Ac,1和Ac,2分别为自主驾驶车辆在区域S1和S2内的纵向势场强度,K为自主驾驶车辆与障碍车辆之间的纵向距离,ρ为自主驾驶车辆在区域S1内的纵向势场收敛幅度,符号“∈”表示“属于”,区域S1定义为自主驾驶车辆距离障碍车辆在纵向安全距离以外区域,区域S2定义为自主驾驶车辆距离障碍车辆在纵向安全距离以内区域和等于纵向安全距离区域;
2)、定义状态列向量
Figure BDA0002438340670000033
和控制列向量u=[δf,ax]T,其中,符号“T”表示“向量转置”,将状态列向量和控制列向量代入式(1),得到状态方程,参见式(4):
x(k+1)=x(k)+Tsf(x(k),u(k)) (4)
其中,f(x(k),u(k))为函数列向量,参见式(5):
Figure BDA0002438340670000032
3)、考虑模型式(4),建立自主驾驶车辆的动态预测模型,参见式(6):
x(i+1|k)=x(i|k)+Tsf(x(i|k),u(i|k)),i=0,1,...,N-1 (6)
其中,x(i|k)为自主驾驶车辆在时刻k对未来时刻k+i的预测状态列向量,u(i|k)为自主驾驶车辆在时刻k对未来时刻k+i的预测控制列向量,正整数N为预测时间窗口;
4)、考虑自主驾驶车辆的动态预测模型式(6),分别定义表征车辆操纵稳定性的车辆前轮侧偏角αf和后轮侧偏角αr,分别参见式(7)和(8):
αf(i|k)=δf(i|k)-[vy(i|k)+ax(i|k)r(i|k)]/vx(i|k) (7)
αr(i|k)=[drr(i|k)-vy(i|k)]/vx(i|k) (8)
再分别定义控制列向量约束、控制增量约束、前轮侧偏角αf和后轮侧偏角αr的边界约束,分别参见式(9)至(12):
umin≤u(i|k)≤umax,i=0,1,...,N-1 (9)
Δumin≤Δu(i|k)≤Δumax,i=0,1,...,N-1 (10)
αfmin≤αf(i|k)≤αfmax,i=0,1,...,N-1 (11)
αrmin≤αr(i|k)≤αrmax,i=0,1,...,N-1 (12)
其中,Δu(i|k)=u(i|k)-u(i-1|k)为相邻预测时刻k+i和k+i-1的预测控制增量,umin、Δumin、αfmin、αrmin分别为控制列向量、控制增量、前轮侧偏角、后轮侧偏角的约束下界,umax、Δumax、αfmax、αrmax分别为控制列向量、控制增量、前轮侧偏角、后轮侧偏角的约束上界;
5)、考虑自主驾驶车辆的动态预测模型式(6),分别定义避障目标函数J1和跟踪目标函数J2,参见式(13)和(14):
Figure BDA0002438340670000041
Figure BDA0002438340670000042
其中,xgoal为车辆轨迹规划目标点,vx,goal为车辆轨迹规划的纵向目标速度,矩阵S和W都为正定的加权矩阵,系数Q和R都为正的加权系数,||·||为2范数;
6)、自主驾驶车辆避障轨迹规划与跟踪控制首先要保证车辆行车安全性,再尽可能恒速跟踪行驶,因此,检测当前时刻k的状态x(k),定义避障轨迹规划与跟踪控制优先级分层预测控制问题,参见式(15)和(16):
Figure BDA0002438340670000043
Figure BDA0002438340670000044
其中,符号“|”表示约束,up(k)为当前时刻k的决策控制变量序列up(k)={u(0|k),u(1|k),…,u(N-1|k)},等式x(0|k)=x(k)为式(15)和(16)的初始条件,J1 *为式(15)的最优值,J2 *为式(16)的最优值;
7)、求解式(16)得对应J2 *的最优解up *(k)={u*(0|k),u*(1|k),…,u*(N-1|k)},参见式(17):
Figure BDA0002438340670000051
取第一个分量u*(0|k)作用于自主驾驶车辆;在下一个采样时刻k+1到达后,检测车辆的运动状态x(k+1),并以该状态更新优化控制问题式(15)和(16)的初始条件,然后优化计算当前时刻的最优解up *(k+1),周而复始,直到自主驾驶车辆运动到目标位置为止。
进一步,所述步骤1)中,纵向安全距离为3m。当然,也可以设置为其他数值。
本发明的技术构思为:用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法基于离散时间车辆动力学模型,定义出描述行车环境的势场函数,将离散时间动力学模型作为预测模型,然后建立有限时域优先级分层优化控制问题,结合滚动优化控制原理,在每个采样时刻计算得到最优控制解,实现基于目标优先级驱动的自主车辆避障轨迹规划与跟踪控制的预测控制。本发明方法的优点是理解简单,通用性强,控制量满足自主驾驶车辆安全避障规划要求和车辆操纵稳定性。
本发明主要执行部分在自主驾驶车辆运动控制计算机上运行实施。本方法应用过程可以大致分为3个阶段:
1、参数设置:分为车辆物理参数和控制器参数,在车辆物理参数导入界面中,输入采样周期Ts,整车质量m,转动惯量Iz,车辆质心到前、后轴的距离df和dr,前轮和后轮的侧偏刚度Ccf和Ccr,分道线势场强度Al,分道线势场值收敛幅度σl,道路边界势场强度Ar,目标点纵向位置Xgoal,正整数ε和β,车辆横向势场收敛幅度σo,纵向势场收敛幅度ρ,区域S1和S2纵向势场强度Ac,1和Ac,2,控制量、控制增量、前轮侧偏角、后轮侧偏角约束下界umin、Δumin、αfmin、αrmin,控制量、控制增量、前轮侧偏角、后轮侧偏角约束上界umax、Δumax、αfmax、αrmax,目标点xgoal和纵向目标速度vx,goal;在控制器参数导入界面中,输入预测时间窗口N,矩阵W、S和系数R、Q;输入参数确认后,由控制计算机将设置数据送入计算机存储单元RAM中保存;
2、离线调试:点击组态界面中的“调试”按钮,进入控制器离线仿真调试阶段,调整组态界面中的预测时间窗口N、矩阵W、S和系数R、Q,观察自主驾驶车辆位置、横摆角、横摆角速度、纵向速度和横向速度的变化,选取自主驾驶车辆路径光滑安全的一组预测时间窗口值和加权矩阵值;预测时间窗口N越大,避障规划路径越平滑;Q越大,越容易避开障碍物;R越大,跟踪误差越小;S越大规划的路径越平稳;W越大规划的路径与道路中心线的位置偏差越小;因此,实际调试R、W、S、Q时,应权衡自主驾驶车辆状态避障效果、跟踪效果、路径平滑、路径靠近中心线之间的综合控制效果;
3、在线运行:点击组态界面“运行”按钮,启动自主驾驶车辆运动控制计算机的CPU读取存储单元RAM中车辆物理参数和控制器参数,并执行“车辆避障轨迹规划与跟踪控制程序”,通过在线测量自主驾驶车辆的纵横向位置、横摆角、横摆角速度和纵横速度,控制输入车辆的纵向加速度和前轮转角,实现自主驾驶车辆安全平稳行驶;在下一个采样时刻到达时,在线测量自主驾驶车辆的实际纵横向位置、横摆角、横摆角速度和纵横速度,重复整个执行过程;周而复始,实现自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制。
本发明明确自主驾驶车辆避障规矩规划与跟踪控制的目标优先级,采用优先级分层排序,优先规划车辆的避障目标轨迹,同时跟踪控制规划的轨迹,以确保自主驾驶车辆在避障过程中的最重要安全性,避免避障目标和跟踪目标复杂的加权系数调整和重新配置的需要。
用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法可以在自主驾驶车辆控制系统组态界面上完成,此过程可以参考本说明书下文中提供的实例应用。与现有自主驾驶车辆避障轨迹规划与跟踪控制方法相比,本发明明确自主驾驶车辆避障目标是最高优先级,跟踪目标是次优先级,优先规划车辆的避障目标轨迹,同时跟踪控制规划的轨迹,以确保自主驾驶车辆在避障过程中的最重要安全性,避免避障目标和跟踪目标复杂的加权系数调整和重新配置的需要。下文具体实施方法以用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法为例说明本发明的实际效果,但本发明的应用范围并不以本实施例中的自主驾驶车辆避障轨迹规划与跟踪控制为限。
本发明的有益效果主要表现在:1、设计简单、容易理解、通用性强;2、在避障轨迹规划与跟踪控制过程中,确保车辆在避障过程中的最重要安全性,避免避障目标和跟踪目标复杂的加权系数调整和重新配置的需要。
附图说明
图1为自主驾驶车辆在X-Y平面运动轨迹曲线的示意图,其中,虚线表示道路中心线,实直线表示车道线;长方体表示障碍车辆,弯曲线表示规划的自主驾驶车辆避障轨迹;
图2为自主驾驶车辆避障轨迹跟踪控制的纵向加速度曲线的示意图;
图3为自主驾驶车辆避障轨迹跟踪控制的前轮转向角曲线的示意图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图3,一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法,所述方法包括如下步骤:
1)、考虑自主驾驶车辆离散时间动力学模型,参见式(1):
Figure BDA0002438340670000081
其中,k为采样时间,Ts为采样周期,X、Y、
Figure BDA0002438340670000082
分别为车辆在全局坐标系下的纵向位置、横向位置、横摆角,vx、vy、r分别为车辆在质心处的纵向速度、横向速度、横摆角速度,m为车辆质量,Iz为车辆绕其垂直轴的转动惯量,df和dr分别是车辆质心到前、后轴的距离,δf和ax分别为前轮转向角和纵向加速度,Ccf和Ccr分别表示前、后轮的侧偏刚度,定义车辆行驶环境势场函数,参见式(2):
Figure BDA0002438340670000083
其中,UA为自主驾驶车辆总势场,常数e为自然指数,Al为分道线势场强度,σl为分道线势场收敛幅度,yl,i为第i条道路分道线的横向位置,Ar为道路边界势场强度,yr,i为第i条道路边界的横向位置,Xgoal为目标点纵向位置,组合参数ε和β为正整数,D是自主驾驶车辆与障碍车辆之间的横向距离,σo为横向势场收敛幅度,Uo为障碍车辆纵向势场函数,参见式(3):
Figure BDA0002438340670000084
其中,Ac,1和Ac,2分别为自主驾驶车辆在区域S1和S2内的纵向势场强度,K为自主驾驶车辆与障碍车辆之间的纵向距离,ρ为自主驾驶车辆在区域S1内的纵向势场收敛幅度,符号“∈”表示“属于”,区域S1定义为自主驾驶车辆距离障碍车辆在纵向安全距离(例如3米)以外区域,区域S2定义为自主驾驶车辆距离障碍车辆在纵向安全距离(例如3米)以内区域和等于纵向安全距离(例如3米)区域,即区域S2为自主驾驶车辆距离障碍车辆小于等于纵向安全距离(例如3米)的区域;
2)、定义状态列向量
Figure BDA0002438340670000091
和控制列向量u=[δf,ax]T,其中,符号“T”表示“向量转置”,将状态列向量和控制列向量代入式(1),得到状态方程,参见式(4):
x(k+1)=x(k)+Tsf(x(k),u(k)) (4)
其中,f(x(k),u(k))为函数列向量,参见式(5):
Figure BDA0002438340670000092
3)、考虑模型式(4),建立自主驾驶车辆的动态预测模型,参见式(6):
x(i+1|k)=x(i|k)+Tsf(x(i|k),u(i|k)),i=0,1,...,N-1 (6)
其中,x(i|k)为自主驾驶车辆在时刻k对未来时刻k+i的预测状态列向量,u(i|k)为自主驾驶车辆在时刻k对未来时刻k+i的预测控制列向量,正整数N为预测时间窗口;
4)、考虑自主驾驶车辆的动态预测模型式(6),分别定义表征车辆操纵稳定性的车辆前轮侧偏角αf和后轮侧偏角αr,分别参见式(7)和(8):
αf(i|k)=δf(i|k)-[vy(i|k)+ax(i|k)r(i|k)]/vx(i|k) (7)
αr(i|k)=[drr(i|k)-vy(i|k)]/vx(i|k) (8)
再分别定义控制列向量约束、控制增量约束、前轮侧偏角αf和后轮侧偏角αr的边界约束,分别参见式(9)至(12):
umin≤u(i|k)≤umax,i=0,1,...,N-1 (9)
Δumin≤Δu(i|k)≤Δumax,i=0,1,...,N-1 (10)
αfmin≤αf(i|k)≤αfmax,i=0,1,...,N-1 (11)
αrmin≤αr(i|k)≤αrmax,i=0,1,...,N-1 (12)
其中,Δu(i|k)=u(i|k)-u(i-1|k)为相邻预测时刻k+i和k+i-1的预测控制增量,umin、Δumin、αfmin、αrmin分别为控制列向量、控制增量、前轮侧偏角、后轮侧偏角的约束下界,umax、Δumax、αfmax、αrmax分别为控制列向量、控制增量、前轮侧偏角、后轮侧偏角的约束上界;
5)、考虑自主驾驶车辆的动态预测模型式(6),分别定义避障目标函数J1和跟踪目标函数J2,参见式(13)和(14):
Figure BDA0002438340670000101
Figure BDA0002438340670000102
其中,xgoal为车辆轨迹规划目标点,vx,goal为车辆轨迹规划的纵向目标速度,矩阵S和W都为正定的加权矩阵,系数Q和R都为正的加权系数,||·||为2范数;
6)、自主驾驶车辆避障轨迹规划与跟踪控制首先要保证车辆行车安全性,再尽可能恒速跟踪行驶,因此,检测当前时刻k的状态x(k),定义避障轨迹规划与跟踪控制优先级分层预测控制问题,参见式(15)和(16):
Figure BDA0002438340670000103
Figure BDA0002438340670000104
其中,符号“|”表示约束,up(k)为当前时刻k的决策控制变量序列up(k)={u(0|k),u(1|k),…,u(N-1|k)},等式x(0|k)=x(k)为式(15)和(16)的初始条件,J1 *为式(15)的最优值,J2 *为式(16)的最优值;
7)、求解式(16)得对应J2 *的最优解up *(k)={u*(0|k),u*(1|k),…,u*(N-1|k)},参见式(17):
Figure BDA0002438340670000105
取第一个分量u*(0|k)作用于自主驾驶车辆;在下一个采样时刻k+1到达后,检测车辆的运动状态x(k+1),并以该状态更新优化控制问题式(15)和(16)的初始条件,然后优化计算当前时刻的最优解up *(k+1),周而复始,直到自主驾驶车辆运动到目标位置为止。
本实例为用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法的实施过程,具体操作过程如下:
1、在参数设置界面中,输入采样周期Ts=50.0ms,整车质量m=1027.0Kg,转动惯量Iz=1536.7Kg·m2,车辆质心到前轴的距离df=1.015m,车辆质心到后轴的距离dr=1.895m,前轮侧偏刚度Ccf=49072.0N/rad,后轮侧偏刚度Ccr=49072.0N/rad,分道线势场强度Al=2.0,分道线势场值收敛幅度σl=1.2,道路边界势场强度Ar=0.6,目标点纵向位置Xgoal=150.0m,正整数ε=1.0,正整数β=0.1,车辆横向势场收敛幅度σo=0.6,纵向势场收敛幅度ρ=0.5,区域S1纵向势场强度Ac,1=15,区域S2纵向势场强度Ac,2=10,控制量、控制增量、前轮侧偏角、后轮侧偏角约束下界umin=[-4,-2]T、Δumin=[-2,-9.4]T、αfmin=αrmin=-2,控制量、控制增量、前轮侧偏角、后轮侧偏角约束上界umax=[2,25]T、Δumax=[2,9.4]T、αfmax=αrmax=2、目标点xgoal=[150,-2,0,25,0,0]T和纵向目标速度vx,goal=25m/s;在控制器参数导入界面中,输入预测时间窗口初始值N=5,矩阵初始值W=5I6和S=5I2,其中I6为矩阵维数是6的单位矩阵,I2为矩阵维数是2的单位矩阵,系数初始值R=5和Q=5;
2、在组态界面上点击“调试”按钮进入调试界面,启动自主车辆运动控制计算机的CPU调用事先编制好的“车辆避障轨迹规划与跟踪控制程序”求解控制器,具体计算过程如下:根据N和R、W、S、Q值,求解优化控制问题式(17),得到最优解up *(k)={u*(0|k),u*(1|k),…,u*(N-1|k)},取第一个分量u*(0|k)作用于自主驾驶车辆;根据N和R、W、S、Q的取值和调整规则,比较自主驾驶车辆行驶路径结果和控制量计算结果,调试预测时间窗口N=15,调试R、W、S、Q得到矩阵W=50I6和S=2000I2和系数R=50和Q=5,将调试结果保存到计算机存储单元RAM中;
3、在线运行:点击组态界面“运行”按钮,启动自主驾驶车辆运动控制计算机的CPU读取存储单元RAM中车辆物理参数和控制器参数,并执行“车辆避障轨迹规划与跟踪控制程序”,通过在线测量自主驾驶车辆的纵横向位置、横摆角、横摆角速度和纵横速度,控制输入车辆的纵向加速度和前轮转角,实现自主驾驶车辆安全平稳行驶;在下一个采样时刻到达时,在线测量自主驾驶车辆的实际纵横向位置、横摆角、横摆角速度和纵横速度,重复整个执行过程;周而复始,实现自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制。
以上阐述的是本发明给出的一个实施例所表现出优良的自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制效果。需要指出的是,上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明做出的任何修改,都落入本发明的保护范围。

Claims (2)

1.一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法,其特征在于,所述方法包括如下步骤:
1)、考虑自主驾驶车辆离散时间动力学模型,参见式(1):
Figure FDA0003645789150000011
其中,k为采样时间,Ts为采样周期,X、Y、
Figure FDA0003645789150000012
分别为车辆在全局坐标系下的纵向位置、横向位置、横摆角,vx、vy、r分别为车辆在质心处的纵向速度、横向速度、横摆角速度,m为车辆质量,Iz为车辆绕其垂直轴的转动惯量,df和dr分别是车辆质心到前、后轴的距离,δf和ax分别为前轮转向角和纵向加速度,Ccf和Ccr分别表示前、后轮的侧偏刚度,定义车辆行驶环境势场函数,参见式(2):
Figure FDA0003645789150000013
其中,UA为自主驾驶车辆总势场,常数e为自然指数,Al为分道线势场强度,σl为分道线势场收敛幅度,yl,i为第i条道路分道线的横向位置,Ar为道路边界势场强度,yr,i为第i条道路边界的横向位置,Xgoal为目标点纵向位置,组合参数ε和β为正整数,D是自主驾驶车辆与障碍车辆之间的横向距离,σo为横向势场收敛幅度,Uo为障碍车辆纵向势场函数,参见式(3):
Figure FDA0003645789150000014
其中,Ac,1和Ac,2分别为自主驾驶车辆在区域S1和S2内的纵向势场强度,K为自主驾驶车辆与障碍车辆之间的纵向距离,ρ为自主驾驶车辆在区域S1内的纵向势场收敛幅度,符号“∈”表示“属于”,区域S1定义为自主驾驶车辆距离障碍车辆在纵向安全距离以外区域,区域S2定义为自主驾驶车辆距离障碍车辆在纵向安全距离以内区域和等于纵向安全距离区域;
2)、定义状态列向量
Figure FDA0003645789150000021
和控制列向量u=[δf,ax]T,其中,符号“T”表示“向量转置”,将状态列向量和控制列向量代入式(1),得到状态方程,参见式(4):
x(k+1)=x(k)+Tsf(x(k),u(k)) (4)
其中,f(x(k),u(k))为函数列向量,参见式(5):
Figure FDA0003645789150000022
3)、考虑模型式(4),建立自主驾驶车辆的动态预测模型,参见式(6):
x(i+1|k)=x(i|k)+Tsf(x(i|k),u(i|k)),i=0,1,...,N-1 (6)
其中,x(i|k)为自主驾驶车辆在时刻k对未来时刻k+i的预测状态列向量,u(i|k)为自主驾驶车辆在时刻k对未来时刻k+i的预测控制列向量,正整数N为预测时间窗口;
4)、考虑自主驾驶车辆的动态预测模型式(6),分别定义表征车辆操纵稳定性的车辆前轮侧偏角αf和后轮侧偏角αr,分别参见式(7)和(8):
αf(i|k)=δf(i|k)-[vy(i|k)+ax(i|k)r(i|k)]/vx(i|k) (7)
αr(i|k)=[drr(i|k)-vy(i|k)]/vx(i|k) (8)
再分别定义控制列向量约束、控制增量约束、前轮侧偏角αf和后轮侧偏角αr的边界约束,分别参见式(9)至(12):
umin≤u(i|k)≤umax,i=0,1,...,N-1 (9)
Δumin≤Δu(i|k)≤Δumax,i=0,1,...,N-1 (10)
αfmin≤αf(i|k)≤αfmax,i=0,1,...,N-1 (11)
αrmin≤αr(i|k)≤αrmax,i=0,1,...,N-1 (12)
其中,Δu(i|k)=u(i|k)-u(i-1|k)为相邻预测时刻k+i和k+i-1的预测控制增量,umin、Δumin、αfmin、αrmin分别为控制列向量、控制增量、前轮侧偏角、后轮侧偏角的约束下界,umax、Δumax、αfmax、αrmax分别为控制列向量、控制增量、前轮侧偏角、后轮侧偏角的约束上界;
5)、考虑自主驾驶车辆的动态预测模型式(6),分别定义避障目标函数J1和跟踪目标函数J2,参见式(13)和(14):
Figure FDA0003645789150000023
Figure FDA0003645789150000024
其中,xgoal为车辆轨迹规划目标点,vx,goal为车辆轨迹规划的纵向目标速度,矩阵S和W都为正定的加权矩阵,系数Q和R都为正的加权系数,||·||为2范数;
6)、自主驾驶车辆避障轨迹规划与跟踪控制首先要保证车辆行车安全性,再尽可能恒速跟踪行驶,因此,检测当前时刻k的状态x(k),定义避障轨迹规划与跟踪控制优先级分层预测控制问题,参见式(15)和(16):
Figure FDA0003645789150000031
Figure FDA0003645789150000032
其中,符号“|”表示约束,up(k)为当前时刻k的决策控制变量序列up(k)={u(0|k),u(1|k),…,u(N-1|k)},等式x(0|k)=x(k)为式(15)和(16)的初始条件,J1 *为式(15)的最优值,J2 *为式(16)的最优值;
7)、求解式(16)得对应J2 *的最优解up *(k)={u*(0|k),u*(1|k),…,u*(N-1|k)},参见式(17):
Figure FDA0003645789150000033
取第一个分量u*(0|k)作用于自主驾驶车辆;在下一个采样时刻k+1到达后,检测车辆的运动状态x(k+1),并以该状态更新优化控制问题式(15)和(16)的初始条件,然后优化计算当前时刻的最优解up *(k+1),周而复始,直到自主驾驶车辆运动到目标位置为止。
2.如权利要求1所述的用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法,其特征在于,所述步骤1)中,纵向安全距离为3m。
CN202010258390.8A 2020-04-03 2020-04-03 一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法 Active CN111538328B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010258390.8A CN111538328B (zh) 2020-04-03 2020-04-03 一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010258390.8A CN111538328B (zh) 2020-04-03 2020-04-03 一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法

Publications (2)

Publication Number Publication Date
CN111538328A CN111538328A (zh) 2020-08-14
CN111538328B true CN111538328B (zh) 2022-07-26

Family

ID=71974916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010258390.8A Active CN111538328B (zh) 2020-04-03 2020-04-03 一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法

Country Status (1)

Country Link
CN (1) CN111538328B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112230659B (zh) * 2020-10-16 2024-09-20 深圳安途智行科技有限公司 精准规划运动轨迹的方法、智能控制设备及自动驾驶车辆
CN113625702B (zh) * 2020-10-28 2023-09-22 北京科技大学 基于二次规划的无人车同时路径跟踪与避障方法
CN112965471B (zh) * 2021-02-10 2022-02-18 大连理工大学 一种考虑角速度约束和改进斥力场的人工势场路径规划方法
EP4296133A4 (en) * 2021-03-01 2024-04-24 Huawei Technologies Co., Ltd. INTELLIGENT DRIVE METHOD AND DEVICE AS WELL AS STORAGE MEDIUM AND COMPUTER PROGRAM
CN113306573B (zh) * 2021-06-23 2022-08-16 湖南大学 一种自动驾驶车辆学习型路径跟踪预测控制方法
CN113268817B (zh) * 2021-06-25 2022-11-01 上海交通大学 基于动态优先级的路口自主车辆调度与控制方法
CN113778082B (zh) * 2021-08-23 2023-08-22 北京科技大学 一种基于自触发机制的无人车轨迹跟踪控制方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107577231A (zh) * 2017-08-28 2018-01-12 驭势科技(北京)有限公司 车辆的控制决策的制定方法、装置和无人驾驶车辆
CN107885932A (zh) * 2017-11-07 2018-04-06 长春工业大学 一种考虑人机和谐的汽车紧急避撞分层式控制方法
CN108099908A (zh) * 2017-12-07 2018-06-01 浙江工业大学 一种车辆自适应巡航优化控制计算方法
CN108909703A (zh) * 2018-06-27 2018-11-30 聊城大学 一种自动驾驶紧急避让的失稳可控域的判定方法
US10168674B1 (en) * 2013-04-22 2019-01-01 National Technology & Engineering Solutions Of Sandia, Llc System and method for operator control of heterogeneous unmanned system teams
CN110126837A (zh) * 2018-02-09 2019-08-16 通用汽车环球科技运作有限责任公司 用于自主车辆运动规划的系统和方法
CN110427827A (zh) * 2019-07-08 2019-11-08 辽宁工程技术大学 一种多尺度感知及全局规划下的自主驾驶网络
CN110780674A (zh) * 2019-12-04 2020-02-11 哈尔滨理工大学 一种提高自动驾驶轨迹跟踪控制的方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10168674B1 (en) * 2013-04-22 2019-01-01 National Technology & Engineering Solutions Of Sandia, Llc System and method for operator control of heterogeneous unmanned system teams
CN107577231A (zh) * 2017-08-28 2018-01-12 驭势科技(北京)有限公司 车辆的控制决策的制定方法、装置和无人驾驶车辆
CN107885932A (zh) * 2017-11-07 2018-04-06 长春工业大学 一种考虑人机和谐的汽车紧急避撞分层式控制方法
CN108099908A (zh) * 2017-12-07 2018-06-01 浙江工业大学 一种车辆自适应巡航优化控制计算方法
CN110126837A (zh) * 2018-02-09 2019-08-16 通用汽车环球科技运作有限责任公司 用于自主车辆运动规划的系统和方法
CN108909703A (zh) * 2018-06-27 2018-11-30 聊城大学 一种自动驾驶紧急避让的失稳可控域的判定方法
CN110427827A (zh) * 2019-07-08 2019-11-08 辽宁工程技术大学 一种多尺度感知及全局规划下的自主驾驶网络
CN110780674A (zh) * 2019-12-04 2020-02-11 哈尔滨理工大学 一种提高自动驾驶轨迹跟踪控制的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Defeng He.Dual-mode Nonlinear MPC via Terminal Control Laws With Free-parameters.《IEEE/CAA JOURNAL OF AUTOMATICA SINICA》.2017,第4卷(第3期), *
何德峰 等.四驱电动汽车驱动力分配阶梯式模型预测控制.《浙江工业大学学报》.2020,第48卷(第1期), *
王兴华 等.一种基于行为的自主/遥控水下机器人共享控制方法.《舰船科学技术》.2020,第42卷(第1期), *
谢朔 等.基于混合式体系结构的无人船控制系统设计.《船海工程》.2016,(第03期), *

Also Published As

Publication number Publication date
CN111538328A (zh) 2020-08-14

Similar Documents

Publication Publication Date Title
CN111538328B (zh) 一种用于自主驾驶车辆避障轨迹规划与跟踪控制的优先级分层预测控制方法
CN111413966B (zh) 一种递进式模型预测无人驾驶规划跟踪协同控制方法
CN111845774B (zh) 一种基于横纵向协调的自动驾驶汽车动态轨迹规划及跟踪方法
CN113320542B (zh) 一种自动驾驶车辆的跟踪控制方法
CN107561942B (zh) 基于模型补偿的智能车辆轨迹跟踪模型预测控制方法
CN111103798B (zh) 一种基于反演滑模控制的agv路径跟踪方法
Awad et al. Model predictive control with fuzzy logic switching for path tracking of autonomous vehicles
Boukens et al. Design of an intelligent optimal neural network-based tracking controller for nonholonomic mobile robot systems
CN110162046B (zh) 基于事件触发型模型预测控制的无人车路径跟随方法
CN110267866B (zh) 用于控制车辆轨迹的装置
Kebbati et al. Lateral control for autonomous wheeled vehicles: A technical review
CN113126623B (zh) 一种考虑输入饱和的自适应动态滑模自动驾驶车辆路径跟踪控制方法
CN113581201A (zh) 一种基于势场模型的无人汽车避撞控制方法及系统
CN114326709B (zh) 一种基于模型预测控制的具有预览特性的智能车辆路径跟踪方法
CN114942642A (zh) 一种无人驾驶汽车轨迹规划方法
CN114779641A (zh) 基于新航向误差定义的环境自适应mpc路径追踪控制方法
CN115343950A (zh) 一种适应复杂路面的车辆路径跟踪控制方法及控制系统
CN117002520A (zh) 一种智能车辆滚动优化非合作博弈决策方法
CN115167135A (zh) 反馈与模型前馈串级的无人驾驶车辆自趋优位姿控制系统
CN115447615A (zh) 基于车辆运动学模型预测控制的轨迹优化方法
Liu et al. Path tracking control of automated vehicles based on adaptive MPC in variable scenarios
CN114415522B (zh) 一种履带车辆平台路径跟踪控制方法及系统
Gao et al. Research on Self-driving Vehicle Path Tracking Adaptive Method Based on Predictive Control
CN115129046B (zh) 一种基于滑模神经网络控制的自动驾驶车辆路径跟踪方法
Ma et al. Optimal Design of MPC Autonomous Vehicle Trajectory Tracking Controller Considering Variable Time Domain

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