CN111258323B - 一种智能车辆轨迹规划与跟踪的联合控制方法 - Google Patents

一种智能车辆轨迹规划与跟踪的联合控制方法 Download PDF

Info

Publication number
CN111258323B
CN111258323B CN202010239134.4A CN202010239134A CN111258323B CN 111258323 B CN111258323 B CN 111258323B CN 202010239134 A CN202010239134 A CN 202010239134A CN 111258323 B CN111258323 B CN 111258323B
Authority
CN
China
Prior art keywords
vehicle
state
model
control
constraint
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
CN202010239134.4A
Other languages
English (en)
Other versions
CN111258323A (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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202010239134.4A priority Critical patent/CN111258323B/zh
Publication of CN111258323A publication Critical patent/CN111258323A/zh
Application granted granted Critical
Publication of CN111258323B publication Critical patent/CN111258323B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/0088—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot characterized by the autonomous decision making process, e.g. artificial intelligence, predefined behaviours
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/02—Control of position or course in two dimensions
    • G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/12—Target-seeking control

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Artificial Intelligence (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
  • Steering Control In Accordance With Driving Conditions (AREA)

Abstract

本发明提供了一种智能车辆轨迹规划与跟踪的联合控制方法，包括：对道路进行参数化；建立横纵向车辆动力学模型；扩充车辆模型，并经过线性化、离散化、增量式转换，形成预测模型；进而根据模型预测控制方法，设计轨迹规划与跟踪的联合控制问题的目标函数；获取当前车辆状态，并更新预测模型中的矩阵参数；根据车辆状态与上一时刻的预测输出状态，建立道路边界、行驶稳定性、状态量约束；将问题转化为二次规划问题形式，并求解最优控制增量组合，并将最优的前轮转角δ*和纵向控制值D*作用于自动驾驶车辆。本发明实现了最优行车线规划与控制的联合处理，降低自动驾驶模块的复杂程度，减少整体内存占用与处理时间，有助于提高自动驾驶赛车的操纵稳定性。

Description

一种智能车辆轨迹规划与跟踪的联合控制方法
技术领域
本发明属于智能汽车控制技术领域，涉及一种智能汽车运动规划与控制方法，特别是涉及一种智能车辆轨迹规划与跟踪的联合控制方法。
背景技术
伴随着计算机信息处理技术和传感器技术的进步，自动驾驶车辆技术迅速发展，并在军事及商业应用领域取得了广泛应用与长足进步。作为对自动驾驶系统的能力检验，自动驾驶和方程式赛车融合的发展方向进一步促进了大学生方程式竞赛的发展，推动了智能驾驶行业对汽车底盘的控制能力与算法鲁棒性的进步。与自动驾驶乘用车的平稳性、安全性、舒适性的目标不同，无人驾驶方程式赛车具有独特的动力学特性和面向竞速的比赛要求。
当赛车在已知赛道上行驶时，自动驾驶系统需执行两个重要过程——轨迹规划和轨迹追踪：轨迹规划模块将可行驶区域(赛道边界信息)规划出最优行驶线(racing line)，并转化为目标路径点；轨迹跟踪模块根据目标路径与自身状态计算控制量，实现赛道行驶控制。
轨迹规划是一类路径搜索问题，常用算法有图搜索、树搜索、智能优化算法等。Juraj Kabzan等根据当前车辆状态和运动模型，生成一系列车辆运行轨迹，并设计代价函数，通过优化方法筛选出最优的运行轨迹，作为赛车行驶线。而Paul A.Theodosis等根据赛车手驾驶经验，将赛车过弯分成三个阶段：制动入弯、最大转向、加速出弯，并表示成三段clothoid螺旋线，进一步结合附着椭圆求解螺旋线系数，得到最优行驶线，提升路径的可执行性与赛车的操控稳定性。
在轨迹跟踪控制方法中，滑模控制、PurePursuit、LQR等方法都能在低速工况下达到较好的效果；而系统地考虑车辆的非线性动力学特性以及复杂道路条件的影响，研究在中高速、大曲率工况下车辆的操控稳定性，成为轨迹跟踪方法发展的关键技术。KrisadaKritayakirana等结合根据赛车手驾驶经验的规划方法，根据道路曲率与摩擦圆，提出自动驾驶汽车在极限工况下的横纵向控制方法，保证了在追踪最佳行车线的同时，轮胎不超过附着极限。在行驶稳定性方面，Matthew Brown等从横向控制的角度提出了车辆状态的滑移稳定性区域，并通过在模型预测控制器中引入对质心侧偏角及航向角变化率的约束，将车辆保持在稳定的状态空间内。但是上述方法假设纵向速度为恒定值，仅限于横向控制。在此基础上，本文推导了考虑纵向控制的稳定性约束，将其作为非等式约束加入控制器中，使其为控制方法提供车辆稳定性保障。
如上所述，大多数轨迹规划方法无法考虑车辆动力学，导致车辆在追踪路径时难以保持行驶稳定。另一方面，基于车辆模型，轨迹规划和轨迹追踪可以使用模型预测控制的方法处理，但尚未有方法考虑了转向稳定性与轮胎附着情况。
发明内容
针对现有技术的不足，本发明提供一种智能车辆轨迹规划与跟踪的联合控制方法，将规划与跟踪控制结合到一个问题中，使得在已知道路边界的条件下，智能汽车可以通过单个模块，实现最优路径规划与轨迹跟踪控制，并当处于轮胎-路面附着极限时，可以稳定安全行驶。
本申请将轨迹规划和轨迹追踪两个模块的功能结合到一个优化问题中，将预测状态作为规划状态，并同时考虑稳定性约束与道路约束，求解最优控制量。如此，可以避免在规划与追踪模块分离时，所规划路径不具备动力学可执行性的问题，保证了所规划轨迹的实时更新性与可执行性；同时降低了自动驾驶模块的复杂程度，减少整体内存占用与处理时间；有助于提高智能车辆的操纵稳定性与驾驶表现。
一种智能车辆轨迹规划与跟踪的联合控制方法，包括以下步骤：
S1，获取道路信息，并对其进行参数化；
S2，建立世界坐标系与车辆坐标系，建立横纵向车辆动力学模型；
S3，车辆模型预处理，得到车辆预测模型；
S4，设计模型预测目标函数；
S5，获取当前时刻T下的车辆状态
Figure BDA0002431969310000031
以
Figure BDA0002431969310000032
更新S3中车辆预测模型中的系数矩阵；
S6，根据
Figure BDA0002431969310000033
和上一时刻T-1的预测输出状态
Figure BDA0002431969310000034
建立道路边界约束、行驶稳定性约束和状态量约束；
S7，对目标函数和各约束进行预处理，将问题转化为二次规划问题并求解，并最终得到T+1时刻的车辆控制量；
S8，输出车辆控制向量对智能车辆进行联合控制。
进一步地，步骤S1具体包括：
①通过智能车辆的感知、建图模块，获取目标道路的两侧边界信息，以离散的坐标点表示；
②将道路两侧边界点按照距离两两配对，获取离散的道路中心线点集P；定义路程sc(下标c表示路径中心线)为路径上任一点到起点的路径曲线总长度，起点处sc＝0；通过sc对P进行三次样条插值计算，分别得出世界坐标Xc(sc),Yc(sc)和航向角
Figure BDA0002431969310000041
即获得道路中心点连续轨迹fc(sc)。
③计算中心线各路程sc处对应的道路边界点的坐标Xl(s),Yl(s),Xr(s),Yr(s)，完成对赛道中线与边界的参数化，其中，坐标下标表示路径点位置，c,l,r分别为中线、道路左边界、道路右边界。
进一步地，步骤S2中所述横纵向车辆动力学模型如下：
Figure BDA0002431969310000042
Figure BDA0002431969310000043
Figure BDA0002431969310000044
Figure BDA0002431969310000045
Figure BDA0002431969310000046
Figure BDA0002431969310000047
其中，Fx＝CmD-Cf-Cw(vx)2
式中，X,Y分别为车辆质心在世界坐标系中的坐标，
Figure BDA0002431969310000048
为车身航向角；m为整车质量、Iz为转动惯量；vx,vy分别为车辆纵向、横向速度；r为横摆角速度；Ff,y,Fr,y分别为前轮、后轮侧偏力，通过等效轮胎模型转换成与侧偏角α的线性关系；δ为前轮转角；lf,lr分别为车辆质心到前轴、后轴的距离；Fx为纵向合外力，Cm,Cf,Cw为车辆参数，通过标定获得。
进一步地，步骤S3具体包括：
①扩充车辆动力学模型，建立车辆状态与参考路径的关系；扩充后的车辆动力学模型可表示为：
Figure BDA0002431969310000051
Figure BDA0002431969310000052
Figure BDA0002431969310000053
Figure BDA0002431969310000054
Figure BDA0002431969310000055
Figure BDA0002431969310000056
Figure BDA0002431969310000057
其中状态变量为
Figure BDA0002431969310000058
输入为u＝[D,δ,vs]T，该动力学模型可以简化表示为
Figure BDA0002431969310000059
②模型线性化：选择线性化参考点，将模型在所述线性化参考点处泰勒公式展开；
③模型离散化；
④将控制量转换成控制增量的形式；
经上述预处理后，得到的车辆预测模型如下：
Figure BDA00024319693100000510
Figure BDA00024319693100000511
其中，
Figure BDA00024319693100000512
为与线性化参考点有关的系数矩阵，
Figure BDA00024319693100000513
为与参考点有关的常数矩阵；系统状态量为
Figure BDA00024319693100000514
表示车辆状态与控制量状态；控制量为Δu＝[ΔD,Δδ,Δvs]T，表示控制增量，其中，状态变量、控制量后的括号“()”表示“时刻”，如x(k)表示k时刻下车辆模型的状态(包括
Figure BDA0002431969310000061
)，u(k-1)表示k-1时刻下车辆模型的控制量状态(包括[D,δ,vs]T)，则ξ(k)表示由k时刻下预测模型的状态(由x(k),u(k-1)组成)。
进一步地，所述扩充车辆动力学模型，具体包括：
①定义当前距离车辆最近的路径点为最近点(Nearest Point)，进一步地，定义该最近点的路程sc为最近点路程sn，在状态量中引入当前路程s，表示当前车辆在参考路径上行进的距离(路程)；
②在控制量中引入路程速度vs，定义为当前路程s对时间的导数，即
Figure BDA0002431969310000062
表示车辆按照参考路径的路程方向上的运行速度；
③定义等效纵向误差
Figure BDA0002431969310000063
(为简便记为
Figure BDA0002431969310000064
)：
Figure BDA0002431969310000065
其中，X,Y为车辆当前位置坐标，Xref(s),Yref(s)为参考路径在路程s处的路径点坐标，
Figure BDA0002431969310000066
为参考路径在路程s处的航向角；各变量的右上角标表示时刻，以k表示当前。等效纵向误差
Figure BDA0002431969310000067
表示当前状态xk中路程sk与最近点路程
Figure BDA0002431969310000068
的相近程度。
进一步地，将上一时刻所解预测状态
Figure BDA0002431969310000069
的时间对应点选择为所述线性化参考点。
进一步地，步骤S4中所述目标函数如下：
Figure BDA00024319693100000610
其中，等号右边第一项为在预测周期末(第Np步，即k＝Np)时车辆运行的路程增量Δs，第二项表示对等效纵向误差el的惩罚，第三项表示对控制增量Δu的惩罚；Np为模型预测控制计算步长，p,Q,R分别表示各个指标的权重。
进一步地，步骤S6中：
所述道路边界约束，将道路边界转化为可线性表示的可行驶区域，具体为：对于每一个预测状态，获取车辆位置与路径中心线上的最近点，将该点处的切线(即航向角方向)分别平移至该点对应的两侧道路边界点上，两相互平行的切线之间形成一个平面可行驶区域，形成对该预测状态ξ的线性约束
Figure BDA0002431969310000071
即为道路边界约束；其中
Figure BDA0002431969310000072
为对状态ξ的坐标加权矩阵，
Figure BDA0002431969310000073
为上下边界矩阵。
所述行驶稳定性约束包括对滑移稳定性约束和轮胎极限约束，并转化成线性不等式约束
Figure BDA0002431969310000074
其中，
Figure BDA0002431969310000075
为稳定状态加权矩阵，
Figure BDA0002431969310000076
为稳定状态边界矩阵。
所述状态量约束包括对状态量ξ(k)和控制增量Δu(k)的约束：|ξ(k)|≤ξmax,|Δu(k)|≤Δumax。ξmax为各状态量的极限值向量，Δumax为各控制增量的极限值向量。所述行驶稳定性约束包括对滑移稳定性约束和轮胎极限约束，并转化成线性不等式约束
Figure BDA0002431969310000077
所述状态量约束包括对状态量ξ(k)和控制增量Δu(k)的约束：|ξ(k)|≤ξmax,|Δu(k)|≤Δumax。
进一步地，所述滑移稳定性约束，对车辆模型中的横向动力学部分进行分析，绘制vy-r的相轨迹图，并根据相轨迹图中的发散区域，得到车辆的滑移稳定性区域，具体为：
Figure BDA0002431969310000081
其中，第一式表示在稳态转向时，轮胎力饱和情况下的极限航向角变化率，即极限横摆角速度
Figure BDA0002431969310000082
第二式是对极限后轮侧偏角αpeak的约束；第三式为对前轮侧偏角的约束；
根据轮胎摩擦圆，所述轮胎极限约束为：
(ax)2+(ay)2≤(gμ)2
其中ax,ay分别为车辆质心处的纵向、横向加速度，μ为路面摩擦系数。
进一步地，步骤S7具体包括：
①将目标函数在“线性化参考点”处线性化，与S2中线性化方法一致；并转化为二次规划问题目标函数的形式；
②建立二次规划问题的约束，具体为：
ξ(0)＝ξmeasure (a)
Figure BDA0002431969310000092
Figure BDA0002431969310000093
Figure BDA0002431969310000094
Figure BDA0002431969310000095
Figure BDA0002431969310000096
(k＝1,...,Np)
其中，ξ(0)为当前车辆状态，(a)式表示当前车辆状态的测量值ξmeasure作为首状态点ξ(0)；(b)式为车辆预测模型；(c)式为道路边界约束；(d)式为行驶稳定性约束；(e),(f)式为对车辆极限状态和底层执行机构的约束；
Figure BDA0002431969310000097
ξ分别代表各状态量的上下界，
Figure BDA0002431969310000098
Δu分别代表各控制增量的上下界；
③求解上述二次规划问题，解得目标控制增量Δu*(包括预测时域内的多组前轮转角增量Δδ和纵向控制值增量ΔD，即Δu*＝[Δδ*,ΔD*])；将Δu*的第一组Δu*(1)＝[Δδ*(1),ΔD*(1)]叠加到当前时刻T的车辆控制状态量
Figure BDA0002431969310000099
上，得到当前最优控制量u*。
即：
Figure BDA00024319693100000910
与现有技术相比，本发明的有益效果为：
1.本发明提出的方法基于模型预测控制理论，将智能车辆的运动规划、轨迹跟踪问题转化为统一的优化问题。相比于规划、跟踪分离的方法，该方法可以避免所规划路径不具备动力学可执行性的问题，保证了所规划轨迹的实时更新性与可执行性；降低自动驾驶模块的复杂程度，减少整体内存占用与处理时间；有助于提高智能车辆的操纵稳定性与驾驶表现。
2.本发明在预测模型中引入了“等效纵向误差”，替代了通过“投影”(即对离散路径点进行遍历，求解距离车辆最近的路径点的过程)求解路径参考点，避免因“投影”操作产生的多重迭代问题，优化了模型预测控制问题的求解过程，较大程度地提高了求解效率。
3.本发明分别引入了横向滑移稳定性约束与纵向轮胎极限约束，避免了横纵向联合控制时因附着极限导致的车辆失稳问题，提高了自动驾驶汽车自主规划与跟踪联合控制的稳定性。
4.本发明建立参考路径与预测模型的关系，通过车辆预测状态，保证所规划轨迹在车辆动力学上的可执行性。
附图说明
图1为本发明实施例的无人驾驶系统构架图。
图2为本发明实施例的等效纵向误差示意图。
图3为本发明实施例的线性化参考点选择方法流程示意图。
图4为本发明实施例的第k步道路边界约束示意图。
图5为本发明实施例的滑移稳定性约束示意图。
图6为本发明实施例中轨迹规划与跟踪联合控制方法的流程图。
实施方式
为了使本发明的目的，技术方案及优点更加清楚明白，以下结合附图及实施实例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
为方便表述，设当前为T时刻，时刻信息在状态变量ξ、解U的右下标表示。模型预测控制方法的控制周期为Ts，预测时域(步长)为Np，即预测控制步数k∈[1，Np]。
一种智能车辆轨迹规划与跟踪的联合控制方法可以在极限情况下使用。该方法可以应用于前轮转向的普通四轮车辆上；如应用于道路交通环境下的乘用车时，仅需将不含障碍物的可行驶区域作为道路信息，即可实现联合控制。为了方便理解，本实施例提供的联合控制方法应用于自动驾驶方程式赛车上，该赛车的无人驾驶系统参见图1。
所述方法包括以下步骤：
S1：获取道路的边界信息，并对其进行参数化：通过道路中心线与两侧边界的离散坐标点Xc(s)，Yc(s)，Xl(s)，Yl(s)，Xr(s)，Yr(s)表示，具体为：
①通过智能车辆的感知、建图模块，获取目标道路的两侧边界信息，以离散的坐标点表示；
②将道路两侧边界点按照距离两两配对，获取离散的道路中心线点集P；定义路程sc(下标c表示路径中心线)为路径上任一点到起点的路径曲线总长度，起点处sc＝0；通过sc对P进行三次样条插值计算，分别得出世界坐标Xc(sc)，Yc(sc)和航向角
Figure BDA0002431969310000111
即获得道路中心点连续轨迹fc(sc)。
③计算中心线各路程sc处对应的道路边界点的坐标(Xc(s)，Yc(s)，Xl(s)，Yl(s)，Xr(s)，Yr(s))。其中，坐标下标表示路径点位置，c，l，r分别为中线、道路左边界、道路右边界。
S2：建立横纵向车辆动力学模型，并完成模型参数的测量与标定，所述模型如下：
Figure BDA0002431969310000121
Figure BDA0002431969310000122
Figure BDA0002431969310000123
Figure BDA0002431969310000124
Figure BDA0002431969310000125
Figure BDA0002431969310000126
其中，X,Y分别为车辆质心在世界坐标系中的坐标，
Figure BDA0002431969310000127
为车身航向角；m为整车质量、Iz为转动惯量；vx,vy分别为车辆纵向、横向速度；r为横摆角速度；Ff,y,Fr,y分别为前轮、后轮侧偏力；δ为前轮转角；lf,lr分别为车辆质心到前轴、后轴的距离；纵向合外力Fx，Fx＝CmD-Cf-Cw(vx)2，Cm,Cf,Cw为车辆参数，通过标定获得。(世界坐标的方向遵循世界坐标系ENU规定，vx,vy,r遵循车身坐标系FLU规定，轮胎侧向力、纵向力的方向遵循轮胎坐标系规定。)
采用Pacejka的Brush轮胎模型，并作小侧偏角、小侧向加速度假设，将其简化为等效侧偏刚度
Figure BDA0002431969310000128
与侧偏角α#的线性模型：
Figure BDA0002431969310000129
式中，
Figure BDA00024319693100001210
为前、后轮的等效侧偏刚度。αf,αr分别代表前、后轮侧偏角。
对m,lf,lr等参数进行测量；根据四个车轮的载荷，估算Iz等参数的值；根据轮胎性能曲线参数，计算
Figure BDA0002431969310000131
设计纵向控制标定实验，完成Cm,Cf,Cw的测算。
S3：车辆模型预处理，得到车辆预测模型，所述车辆模型预处理包括：模型扩充、线性化、离散化和增量式转换，具体为：
①扩充模型，建立车辆状态与参考路径的关系，同时避免计算最近点时的多重迭代问题，具体为：
a)定义当前距离车辆最近的路径点为最近点(Nearest Point)，进一步地，定义该最近点的路程sc为最近点路程sn；在状态量中引入当前路程s，表示当前车辆在参考路径上行进的距离(路程)；
b)在控制量中引入路程速度vs，定义为当前路程s对时间的导数，即
Figure BDA0002431969310000132
表示车辆按照参考路径的路程方向上的运行速度；
c)定义等效纵向误差
Figure BDA0002431969310000133
(为简便记为
Figure BDA0002431969310000134
)：
Figure BDA0002431969310000135
其中，X,Y为车辆当前位置坐标，Xref(s),Yref(s)为参考路径在路程s处的路径点坐标，
Figure BDA0002431969310000136
为参考路径在路程s处的航向角；各变量的右上角标表示时刻，以k表示当前。等效纵向误差
Figure BDA0002431969310000137
表示当前xk中路程s与路径最近点sn的相近程度，如图2中的粗实线箭头。在目标函数中设置对
Figure BDA0002431969310000138
的惩罚，保证MPC解中的s(k)与
Figure BDA0002431969310000139
的误差较小。将s引入状态量xk，建立了车辆状态xk与参考路径的关系，便于计算各项误差；通过计算vs对sk进行更新，避免多重迭代问题。
②线性化：获取当前状态的线性化参考点，并将模型线性化——在该点处泰勒公式展开，忽略高阶项。获取线性化参考点的方法为“将上一时刻所解预测状态
Figure BDA0002431969310000141
的时间对应点作为线性化参考点”，以下以“T时刻”为例，具体说明线性化参考点的获取方法：
如图3所示：从上一时刻(T-1时刻)的解
Figure BDA0002431969310000142
中，获得车辆预测状态
Figure BDA0002431969310000143
(三角形点)；在T时刻时，将从T-1时刻解中的第k+1步预测状态
Figure BDA0002431969310000144
作为T时刻第k步的线性化参考点(空心圆点)，并在预测时域Np-1内一一对应；T时刻第Np步的参考点可与第Np-1步一致。
③使用零阶保持进行离散化；
④将控制量转换成控制增量的形式；
如上处理后，得到的车辆预测模型如下：
Figure BDA0002431969310000145
Figure BDA0002431969310000146
其中，
Figure BDA0002431969310000147
为与线性化参考点有关的系数矩阵，
Figure BDA0002431969310000148
为与参考点有关的常数矩阵；系统状态量为
Figure BDA0002431969310000149
表示车辆状态与控制量状态；控制量为Δu＝[ΔD,Δδ,Δvs]T，表示控制增量。其中，系统状态量、控制量后的括号“()”表示“时刻”，如x(k)表示k时刻下车辆模型的状态(包括
Figure BDA00024319693100001410
)，u(k-1)表示k-1时刻下车辆模型的控制量状态(包括[D,δ,vs]T)，则ξ(k)表示由k时刻下预测模型的状态(由x(k),u(k-1)组成)。
S4：根据模型预测控制方法，建立轨迹规划与跟踪的联合控控制问题，包括：
为了“根据已知赛道边界，实现轨迹规划与跟踪的联合控制”，设计模型预测目标函数为：
Figure BDA0002431969310000151
其中，等号右边第一项为在预测周期末(第Np步，即k＝Np)时车辆运行的路程增量Δs，第二项表示对等效纵向误差el的惩罚，第三项表示对控制增量Δu的惩罚；Np为模型预测控制计算步长，p,Q,R分别表示各个指标的权重，每个指标的权重可根据实际情况赋值。
目标函数对应如下控制目标：在预测周期内，使车辆在参考路径方向的运行路程s最大化，进而达到全局行驶路径最短、圈速最快的目的；其次，要使等效纵向误差
Figure BDA0002431969310000157
的值尽量小，保证路程sk求解的准确性；另外，为减少频繁的操作、确保车辆执行机构的稳定执行，对控制增量Δu进行惩罚。
S5：经过传感器测量与滤波处理，获取当前时刻(相关变量以下标“T”表示)下的车辆状态
Figure BDA0002431969310000152
以
Figure BDA0002431969310000153
更新S3中预测模型的系数矩阵。
S6：根据
Figure BDA0002431969310000154
和上一时刻(T-1时刻)的预测输出状态
Figure BDA0002431969310000155
建立道路边界约束、行驶稳定性约束和状态量约束：
①道路边界约束：将道路边界转化为可线性表示的可行驶区域，并转化成线性不等式约束
Figure BDA0002431969310000156
建立道路边界约束的具体方法是：对于每一个预测状态，获取车辆位置在中心线上的最近点，将该点处的切线分别平移至该点对应的两侧道路边界点上，两相互平行的切线之间形成一个平面可行驶区域形成对该预测状态的线性约束，即为道路边界约束。具体为：
对照附图4，过程如下，：
a)与S3②中获取线性化参考点方法类似，在形成T时刻第k步的道路边界约束时，使用T-1时刻解出的对应预测状态点
Figure BDA0002431969310000161
作为状态参考点
Figure BDA0002431969310000162
根据其路程
Figure BDA0002431969310000163
获取路径中线上的路径参考点，获取对应坐标
Figure BDA0002431969310000164
与航向角
Figure BDA0002431969310000165
b)将路径参考点处的切线平移至该点对应的(两侧)道路边界点上，形成第k步的上下边界约束(粗实线)。
c)获取一系列预测时域内的状态参考点与路径参考点，生成一系列的两端边界约束，组合成为T时刻下的环境约束。
因此，将道路边界约束表示不等式约束：
Figure BDA0002431969310000166
Figure BDA0002431969310000167
Figure BDA0002431969310000168
Figure BDA0002431969310000169
Figure BDA00024319693100001610
Figure BDA00024319693100001611
其中，(Xl,Yl)、(Xr,Yr)为T-1时刻的解中
Figure BDA00024319693100001612
里
Figure BDA00024319693100001613
所对应的一组(两侧)道路边界点坐标。
Figure BDA00024319693100001614
表示可行区域的切线斜率，bu,bl分别表示可行区域上下界。
②行驶稳定性约束：包括对滑移稳定性、轮胎附着极限的约束，并转化成线性不等式约束
Figure BDA0002431969310000171
具体如下：
a)滑移稳定性约束：对车辆模型中的横向动力学部分进行分析，绘制vy-r的相轨迹图，并根据相轨迹图中的发散区域，得到车辆的滑移稳定性区域，如图5。具体为：
Figure BDA0002431969310000172
其中，第一式表示在稳态转向时，轮胎力饱和情况下的极限航向角变化率，即极限横摆角速度
Figure BDA0002431969310000173
对应图5的①、③边界；第二式是对极限后轮侧偏角αpeak的约束，目的为在r较大时限制vy，避免进入发散区域，对应图5的②、④边界；第三式为对前轮侧偏角的约束，将其限制在轮胎模型的线性区域，确保等效轮胎模型的准确性。
b)轮胎极限约束：
如果将车辆视为一个质点，根据轮胎摩擦圆，车辆的极限加速度将取决于轮胎与路面的摩擦系数μ，即轮胎摩擦极限约束为：
(ax)2+(ay)2≤(gμ)2
其中ax,ay分别为车辆质心处的纵向、横向加速度。
为便于将模型预测控制问题转化为二次规划问题求解、将其表示为线性约束的形式，在求解控制量Dk时，将ξk-1的车辆状态量vx k-1,vy k-1,rk-1视为常量，将轮胎附着极限约束转化为仅对纵向控制量D的线性约束，如下：
Dmin (k)≤CmDk≤Dmax (k)
其中：
Dmax (k)＝Dmax(ξ(k-1))
＝[maf(ξ(k-1))-Flon(ξ(k-1))]
Dmin (k)＝Dmin(ξ(k-1))
＝[-maf(ξ(k-1))-Flon(ξ(k-1))]
Figure BDA0002431969310000181
Flon＝-Cf-Cw(vx)2-Ff,ysin(δ)+mvyr
其中，ξ(k-1)表示上一时刻解得的预测状态，Flon为除纵向控制量外的纵向合力。
该约束表示通过牺牲“纵向控制”来配合“横向控制”，从而保证轮胎力不超过附着极限。
而横向控制的稳定性由滑移稳定性约束保证，由此形成横纵向的行驶稳定性约束，使用线性不等式约束表示：
Figure BDA0002431969310000182
③状态量约束：包括对状态量ξ(k)、控制增量Δu(k)的约束：
|ξ(k)|≤ξmax,|Δu(k)|≤Δumax。
S7：对目标函数和约束进行预处理，将问题转化为二次规划问题形式；求解QP问题，求解最优控制增量组合Δu*。
①将目标函数在S3的参考点处线性化，与S2中线性化方法一致；并转化为QP问题目标函数的形式，具体为：
Figure BDA0002431969310000191
式中待优化的变量是控制量增量Δu＝[ΔD,Δδ,Δvs]T、控制量所对应的预测状态ξ、滑移约束的松弛因子ssh。松弛因子的存在是为了通过允许高速车辆出现暂时性滑移现象以保证此优化问题始终存在可行解。对车辆滑移约束施加一定的松弛因子，可保证当道路的无碰撞轨迹与稳定性约束发生冲突时,优先保证车辆在可行区域内行驶,即允许车辆暂时发生滑移。
Figure BDA0002431969310000192
矩阵为S4中三个控制目标在线性化后的一次项与二次项组合。
②将S3所得车辆预测模型表示为等式约束，加上S6所得非等式线性约束，联立成为QP问题的约束，具体为：
ξ(0)＝ξmeasure (a)
Figure BDA0002431969310000193
Figure BDA0002431969310000194
Figure BDA0002431969310000195
Figure BDA0002431969310000196
Figure BDA0002431969310000197
(k＝1,...,Np)
其中，ξ(0)为当前车辆状态，(a)式表示使用S3所得当前车辆状态的测量值ξmeasure作为首状态点ξ(0)；(b)式为S4中所得预处理后的车辆动力学模型；(c)式为道路边界约束；(d)式为行驶稳定性约束；(e),(f)式为对车辆极限状态和底层执行机构的约束。
③求解上述QP问题，解得目标控制增量Δu*，包括预测时域内的多组前轮转角增量Δδ和纵向控制值增量ΔD，即Δu*＝[Δδ*,ΔD*]；将Δu*的第一组Δu*(1)＝[Δδ*(1),ΔD*(1)]叠加到当前时刻T的车辆控制状态量
Figure BDA0002431969310000201
上，得到当前最优控制量u*(u*＝[δ*,D*]，包括最优的前轮转角δ*和纵向控制值D*)。即：
Figure BDA0002431969310000202
S8：输出车辆控制向量u*，将最优的前轮转角δ*和纵向控制值D*通过底层控制执行器作用于自动驾驶车辆，并重复S5～S8，实现赛车的自主轨迹规划与跟踪的联合控制。
本实施例实现了最优行车线规划与控制的联合处理，降低自动驾驶模块的复杂程度，减少整体内存占用与处理时间，有助于提高自动驾驶赛车的操纵稳定性与驾驶表现。
本领域普通技术人员可以理解：以上实施例仅用于说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求限定的范围。

Claims (7)

1.一种智能车辆轨迹规划与跟踪的联合控制方法，其特征在于，包括以下步骤：
S1，获取道路信息，并对其进行参数化；
S2，建立世界坐标系与车辆坐标系，建立横纵向车辆动力学模型；
S3，车辆模型预处理，得到车辆预测模型；
S4，设计模型预测目标函数；
S5，获取当前时刻T下的车辆状态
Figure FDA0003201822010000011
以
Figure FDA0003201822010000012
更新S3中车辆预测模型中的系数矩阵；
S6，根据
Figure FDA0003201822010000013
和上一时刻T-1的预测输出状态
Figure FDA0003201822010000014
建立道路边界约束、行驶稳定性约束和状态量约束；
S7，对目标函数和各约束进行预处理，将问题转化为二次规划问题并求解，并最终得到T+1时刻的车辆控制量；
S8，输出车辆控制向量对智能车辆进行联合控制；
其中，步骤S3具体包括：
①扩充车辆动力学模型，建立车辆状态与参考路径的关系；扩充后的车辆动力学模型表示为：
Figure FDA0003201822010000015
Figure FDA0003201822010000016
Figure FDA0003201822010000017
Figure FDA0003201822010000018
Figure FDA0003201822010000019
Figure FDA00032018220100000110
Figure FDA00032018220100000111
其中状态变量为
Figure FDA0003201822010000021
输入为u＝[D,δ,vs]T，该动力学模型简化表示为
Figure FDA0003201822010000022
②模型线性化：选择线性化参考点，将模型在所述线性化参考点处泰勒公式展开；
③模型离散化；
④将控制量转换成控制增量的形式；
经①～④预处理后，得到的车辆预测模型如下：
Figure FDA0003201822010000023
Figure FDA0003201822010000024
其中，X,Y分别为车辆质心在世界坐标系中的坐标，
Figure FDA0003201822010000025
为车身航向角；m为整车质量、Iz为转动惯量；vx,vy分别为车辆纵向、横向速度；r为横摆角速度；Ff,y,Fr,y分别为前轮、后轮侧偏力，通过等效轮胎模型转换成与侧偏角α的线性关系；δ为前轮转角；lf,lr分别为车辆质心到前轴、后轴的距离；Fx为纵向合外力，
Figure FDA0003201822010000026
为与线性化参考点有关的系数矩阵，
Figure FDA0003201822010000027
为与参考点有关的常数矩阵；系统状态量为
Figure FDA0003201822010000028
表示车辆状态与控制量状态，控制量为Δu＝[ΔD,Δδ,Δvs]T，表示控制增量，其中，系统状态量、控制量后的括号“()”表示“时刻”，x(k)表示k时刻下车辆模型的状态，u(k-1)表示k-1时刻下车辆模型的控制量状态，ξ(k)表示由k时刻下预测模型的状态，D为纵向控制量，s为当前路程，
Figure FDA0003201822010000029
为s的导数，vs为路程速度；
步骤S4中所述目标函数如下：
Figure FDA0003201822010000031
式中，等号右边第一项为在预测周期末第Np步，即k＝Np时车辆运行的路程增量Δs，第二项表示对等效纵向误差el的惩罚，第三项表示对控制增量Δu的惩罚；Np为模型预测控制计算步长，p,Q,R分别表示各个指标的权重，s为当前路程；
其中，所述扩充车辆动力学模型，具体包括：
①定义当前距离车辆最近的路径点为最近点，定义该最近点的路程sc为最近点路程sn；在状态量中引入当前路程s，表示当前车辆在参考路径上行进的距离；
②在控制量中引入路程速度vs，定义为当前路程s对时间的导数，即
Figure FDA0003201822010000032
表示车辆按照参考路径的路程方向上的运行速度；
③定义等效纵向误差
Figure FDA0003201822010000033
为简便记为
Figure FDA0003201822010000034
Figure FDA0003201822010000035
其中，X,Y为车辆当前位置坐标，Xref(s),Yref(s)为参考路径在路程s处的路径点坐标，
Figure FDA0003201822010000036
为参考路径在路程s处的航向角；各变量的右上角标表示时刻，以k表示当前，等效纵向误差
Figure FDA0003201822010000037
表示当前状态xk中路程sk与最近点路程
Figure FDA0003201822010000038
的相近程度。
2.根据权利要求1所述的一种智能车辆轨迹规划与跟踪的联合控制方法，其特征在于，步骤S1具体包括：
①通过智能车辆的感知、建图模块，获取目标道路的两侧边界信息，以离散的坐标点表示；
②将道路两侧边界点按照距离两两配对，获取离散的道路中心线点集P；定义路程sc为路径上任一点到起点的路径曲线总长度，起点处sc＝0，其中，下标c表示路径中心线；通过sc对P进行三次样条插值计算，分别得出世界坐标Xc(sc)，Yc(sc)和航向角
Figure FDA0003201822010000041
即获得道路中心点连续轨迹fc(sc)，
③计算中心线各路程sc处对应的道路边界点的坐标Xl(s)，Yl(s)，Xr(s)，Yr(s)，完成对赛道中线与边界的参数化，其中，坐标下标表示路径点位置，c，l，r分别为中线、道路左边界、道路右边界。
3.根据权利要求1所述的一种智能车辆轨迹规划与跟踪的联合控制方法，其特征在于，步骤S2中所述横纵向车辆动力学模型如下：
Figure FDA0003201822010000042
Figure FDA0003201822010000043
Figure FDA0003201822010000044
Figure FDA0003201822010000045
Figure FDA0003201822010000046
Figure FDA0003201822010000047
其中，Fx＝CmD-Cf-Cw(vx)2
式中，X，Y分别为车辆质心在世界坐标系中的坐标，
Figure FDA0003201822010000048
为车身航向角；m为整车质量、Iz为转动惯量；vx，vy分别为车辆纵向、横向速度；r为横摆角速度；Ff，y，Fr，y分别为前轮、后轮侧偏力，通过等效轮胎模型转换成与侧偏角α的线性关系；δ为前轮转角；lf，lr分别为车辆质心到前轴、后轴的距离；Fx为纵向合外力，Cm,Cf,Cw为车辆参数，通过标定获得，D为纵向控制量。
4.根据权利要求1所述的一种智能车辆轨迹规划与跟踪的联合控制方法，其特征在于：将上一时刻所解预测状态
Figure FDA0003201822010000051
的时间对应点选择为所述线性化参考点。
5.根据权利要求1所述的一种智能车辆轨迹规划与跟踪的联合控制方法，其特征在于，步骤S6中，
所述道路边界约束，将道路边界转化为可线性表示的可行驶区域，具体为：对于每一个预测状态，获取车辆位置与路径中心线上的最近点，将该点处的切线即航向角方向分别平移至该点对应的两侧道路边界点上，两相互平行的切线之间形成一个平面可行驶区域，形成对该预测状态ξ的线性约束
Figure FDA0003201822010000052
即为道路边界约束；其中
Figure FDA0003201822010000053
为对状态ξ的坐标加权矩阵，
Figure FDA0003201822010000054
为上下边界矩阵；
所述行驶稳定性约束包括对滑移稳定性约束和轮胎极限约束，并转化成线性不等式约束
Figure FDA0003201822010000055
其中，
Figure FDA0003201822010000056
为稳定状态加权矩阵，
Figure FDA0003201822010000057
为稳定状态边界矩阵；
所述状态量约束包括对状态量ξ(k)和控制增量Δu(k)的约束：|ξ(k)|≤ξmax,|Δu(k)|≤Δumax，ξmax为各状态量的极限值向量，Δumax为各控制增量的极限值向量。
6.根据权利要求5所述的一种智能车辆轨迹规划与跟踪的联合控制方法，其特征在于，
所述滑移稳定性约束，对车辆模型中的横向动力学部分进行分析，绘制vy-r的相轨迹图，并根据相轨迹图中的发散区域，得到车辆的滑移稳定性区域，具体为：
Figure FDA0003201822010000061
其中，第一式表示在稳态转向时，轮胎力饱和情况下的极限航向角变化率，即极限横摆角速度
Figure FDA0003201822010000062
第二式是对极限后轮侧偏角αpeak的约束；第三式为对前轮侧偏角的约束，g为重力加速度；
根据轮胎摩擦圆，所述轮胎极限约束为：
(ax)2+(ay)2≤(gμ)2
其中ax,ay分别为车辆质心处的纵向、横向加速度，μ为路面摩擦系数，vx,vy分别为车辆纵向、横向速度，r为横摆角速度，lf,lr分别为车辆质心到前轴、后轴的距离，δ为前轮转角。
7.根据权利要求6所述的一种智能车辆轨迹规划与跟踪的联合控制方法，其特征在于，步骤S7具体包括：
①将目标函数在“线性化参考点”处线性化，与S2中线性化方法一致；并转化为二次规划问题目标函数的形式；
②建立二次规划问题的约束，具体为：
ξ(0)＝ξmeasure (a)
Figure FDA0003201822010000071
Figure FDA0003201822010000072
Figure FDA0003201822010000073
Figure FDA0003201822010000074
Figure FDA0003201822010000075
其中，k＝1,...,Np
其中，ξ(0)为当前车辆状态；(a)式表示当前车辆状态的测量值ξmeasure作为首状态点ξ(0)；(b)式为车辆预测模型；(c)式为道路边界约束；(d)式为行驶稳定性约束；(e),(f)式为对车辆极限状态和底层执行机构的约束；
Figure FDA0003201822010000076
ξ分别代表各状态量的上下界，
Figure FDA0003201822010000077
Δu分别代表各控制增量的上下界，
Figure FDA0003201822010000078
为与线性化参考点有关的系数矩阵，
Figure FDA00032018220100000712
为与参考点有关的常数矩阵；
③求解上述二次规划问题，解得目标控制增量Δu*，包括预测时域内的多组前轮转角增量Δδ和纵向控制值增量ΔD，即Δu*＝[Δδ*,ΔD*]；将Δu*的第一组Δu*(1)＝[Δδ*(1),ΔD*(1)]叠加到当前时刻T的车辆控制状态量
Figure FDA00032018220100000710
上，得到当前最优控制量u*，即：
Figure FDA00032018220100000711
CN202010239134.4A 2020-03-30 2020-03-30 一种智能车辆轨迹规划与跟踪的联合控制方法 Active CN111258323B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010239134.4A CN111258323B (zh) 2020-03-30 2020-03-30 一种智能车辆轨迹规划与跟踪的联合控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010239134.4A CN111258323B (zh) 2020-03-30 2020-03-30 一种智能车辆轨迹规划与跟踪的联合控制方法

Publications (2)

Publication Number Publication Date
CN111258323A CN111258323A (zh) 2020-06-09
CN111258323B true CN111258323B (zh) 2021-10-26

Family

ID=70949933

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010239134.4A Active CN111258323B (zh) 2020-03-30 2020-03-30 一种智能车辆轨迹规划与跟踪的联合控制方法

Country Status (1)

Country Link
CN (1) CN111258323B (zh)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111665853B (zh) * 2020-07-07 2023-03-31 中国人民解放军国防科技大学 一种面向规划控制联合优化的无人车辆运动规划方法
CN111845774B (zh) * 2020-07-20 2021-12-03 上海大学 一种基于横纵向协调的自动驾驶汽车动态轨迹规划及跟踪方法
CN111857152A (zh) * 2020-07-31 2020-10-30 北京京东乾石科技有限公司 用于生成车辆控制信息的方法和装置
CN114185337B (zh) * 2020-09-15 2023-08-04 宇通客车股份有限公司 一种车辆、车辆预碰撞检测方法及装置
CN112193318A (zh) * 2020-10-15 2021-01-08 北京航天发射技术研究所 车辆路径控制方法、装置、设备及计算机可读存储介质
CN114442601A (zh) * 2020-11-06 2022-05-06 郑州宇通客车股份有限公司 一种无人驾驶车辆跟踪控制方法及装置
CN112578672B (zh) * 2020-12-16 2022-12-09 吉林大学青岛汽车研究院 基于底盘非线性的无人驾驶汽车轨迹控制系统及其轨迹控制方法
CN112731925B (zh) * 2020-12-21 2024-03-15 浙江科技学院 用于无人驾驶方程式赛车锥桶识别和路径规划及控制方法
CN112572443B (zh) * 2020-12-22 2021-12-07 北京理工大学 一种高速公路车辆换道实时避撞轨迹规划方法及系统
CN112706785B (zh) * 2021-01-29 2023-03-28 重庆长安汽车股份有限公司 一种自动驾驶车辆行驶环境认知目标选择方法、装置及存储介质
CN112965498A (zh) * 2021-02-23 2021-06-15 上海电机学院 一种智能扫路车运行轨迹跟踪控制方法
CN113029151B (zh) * 2021-03-15 2023-04-14 齐鲁工业大学 一种智能车辆路径规划方法
CN113076596B (zh) * 2021-03-31 2023-03-28 同济大学 融合车辆模型和轮胎力曲线的轨迹跟随并行预测方法
CN113211456B (zh) * 2021-04-09 2023-06-06 江苏大学 一种喷砂除锈并联机器人移动平台轨迹跟踪控制方法
CN113211438B (zh) * 2021-05-08 2023-06-16 东方红卫星移动通信有限公司 基于预瞄距离自适应的轮式机器人控制方法及控制系统
CN113867330B (zh) * 2021-05-11 2023-05-02 吉林大学 一种基于多自由度预测模型在任意路径下实现车辆漂移的控制方法
CN113320542B (zh) * 2021-06-24 2022-05-17 厦门大学 一种自动驾驶车辆的跟踪控制方法
CN113788008B (zh) * 2021-07-20 2023-06-27 浙江万安科技股份有限公司 一种基于纳什-mpc的集成底盘轨迹跟踪控制方法
CN113467480B (zh) * 2021-08-09 2024-02-13 广东工业大学 一种用于无人驾驶方程式的全局路径规划算法
CN113561994B (zh) * 2021-08-13 2022-04-05 北京三快在线科技有限公司 轨迹规划方法、装置、存储介质及电子设备
CN113985889B (zh) * 2021-11-11 2024-04-30 北京三快在线科技有限公司 一种轨迹规划方法、装置、存储介质及电子设备
CN113985895B (zh) * 2021-11-29 2024-05-03 佛山市毕佳索智能科技有限公司 一种基于最优化的agv路径跟踪方法
CN114200838B (zh) * 2021-12-08 2022-07-12 青岛中鸿重型机械有限公司 一种智能电动铲运机控制方法
CN114435399B (zh) * 2022-01-27 2023-09-12 上海工程技术大学 基于预测模型的自动驾驶汽车稳定性路径跟踪方法
DE102022102118A1 (de) 2022-01-31 2023-08-03 Dr. Ing. H.C. F. Porsche Aktiengesellschaft Verfahren und Vorrichtung zum Erzeugen von Trajektorien auf einer gegebenen Strecke für mehrere Fahrzeuge mit gegebenen Leistungsspezifikationen
CN115465332B (zh) * 2022-10-27 2023-12-29 北京交通大学 一种动车组列车控制方法、系统、电子设备及存储介质
CN116382260B (zh) * 2023-02-13 2023-11-28 大连理工大学 一种基于安全航行通道的水面舰艇停泊轨迹规划方法
CN116300478B (zh) * 2023-05-22 2023-08-18 上海友道智途科技有限公司 基于变量分离微分动态规划的轨迹生成方法、介质及设备
CN116279596B (zh) * 2023-05-26 2023-08-04 禾多科技（北京）有限公司 车辆控制方法、装置、电子设备和计算机可读介质
CN117131612B (zh) * 2023-10-26 2024-01-16 浙江大学 车辆倾覆预测方法、系统和电子设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9669828B2 (en) * 2012-06-01 2017-06-06 Toyota Motor Engineering & Manufacturing North America, Inc. Cooperative driving and collision avoidance by distributed receding horizon control
EP2848487B1 (en) * 2013-09-12 2020-03-18 Volvo Car Corporation Manoeuvre generation for automated drive
CN105584479B (zh) * 2016-01-18 2018-10-19 北京理工大学 一种面向脑控车辆的模型预测控制方法和利用该方法的脑控车辆
US10591926B2 (en) * 2017-09-18 2020-03-17 Baidu Usa Llc Smooth road reference for autonomous driving vehicles based on 2D constrained smoothing spline
CN110329255B (zh) * 2019-07-19 2020-11-13 中汽研（天津）汽车工程研究院有限公司 一种基于人机协同策略的车道偏离辅助控制方法

Also Published As

Publication number Publication date
CN111258323A (zh) 2020-06-09

Similar Documents

Publication Publication Date Title
CN111258323B (zh) 一种智能车辆轨迹规划与跟踪的联合控制方法
CN111845774B (zh) 一种基于横纵向协调的自动驾驶汽车动态轨迹规划及跟踪方法
hwan Jeon et al. Optimal motion planning with the half-car dynamical model for autonomous high-speed driving
Zhang et al. Evaluating model predictive path following and yaw stability controllers for over-actuated autonomous electric vehicles
Altché et al. High-speed trajectory planning for autonomous vehicles using a simple dynamic model
CN112622903B (zh) 一种车辆跟随驾驶环境下自主车辆的纵向和横向控制方法
CN108107732A (zh) 主动前轮转向和直接横摆力矩联合的汽车稳定性控制方法
Chen et al. Design and experimental evaluation of an efficient MPC-based lateral motion controller considering path preview for autonomous vehicles
CN110217227A (zh) 一种适用于冰雪道路工况的转向制动联合避撞控制方法
CN116627044B (zh) 行驶轨迹预测控制方法
Kebbati et al. Lateral control for autonomous wheeled vehicles: A technical review
Kapania Trajectory planning and control for an autonomous race vehicle
CN112578672A (zh) 基于底盘非线性的无人驾驶汽车轨迹控制系统及其轨迹控制方法
Hegedüs et al. Motion planning for highly automated road vehicles with a hybrid approach using nonlinear optimization and artificial neural networks
CN116552550A (zh) 基于参数不确定性和横摆稳定性的车辆轨迹跟踪控制系统
Feraco et al. Combined lane keeping and longitudinal speed control for autonomous driving
Kone Lateral and longitudinal control of an autonomous racing vehicle.
CN116834754A (zh) 一种自动驾驶车辆自适应调速的横纵协同控制方法
Kovacs et al. Integrated path planning and lateral-longitudinal control for autonomous electric vehicles
Li et al. Path tracking control based on the prediction of tire state stiffness using the optimized steering sequence
CN112298193B (zh) 快速实时的后轮主动转向预测控制方法
CN114879698A (zh) 基于改进人工势场和mpc的机器人驾驶车辆避障方法
Goel In Complete Control Simultaneous Path Speed and Sideslip Angle Control of a Drifting Automobile
CN114995426A (zh) 一种基于神经网络动力学模型的无人驾驶车辆轨迹跟踪控制方法、系统及车载控制设备
CN111857112B (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