CN115042770A - 一种基于分布式鲁棒模型预测的车辆队列横向控制方法 - Google Patents

一种基于分布式鲁棒模型预测的车辆队列横向控制方法 Download PDF

Info

Publication number
CN115042770A
CN115042770A CN202210681057.7A CN202210681057A CN115042770A CN 115042770 A CN115042770 A CN 115042770A CN 202210681057 A CN202210681057 A CN 202210681057A CN 115042770 A CN115042770 A CN 115042770A
Authority
CN
China
Prior art keywords
vehicle
lateral
angle
following
queue
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.)
Granted
Application number
CN202210681057.7A
Other languages
English (en)
Other versions
CN115042770B (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.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN202210681057.7A priority Critical patent/CN115042770B/zh
Publication of CN115042770A publication Critical patent/CN115042770A/zh
Application granted granted Critical
Publication of CN115042770B publication Critical patent/CN115042770B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W30/00Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units, or advanced driver assistance systems for ensuring comfort, stability and safety or drive control systems for propelling or retarding the vehicle
    • B60W30/02Control of vehicle driving stability
    • B60W30/04Control of vehicle driving stability related to roll-over prevention
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W30/00Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units, or advanced driver assistance systems for ensuring comfort, stability and safety or drive control systems for propelling or retarding the vehicle
    • B60W30/02Control of vehicle driving stability
    • B60W30/04Control of vehicle driving stability related to roll-over prevention
    • B60W2030/043Control of vehicle driving stability related to roll-over prevention about the roll axis
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W50/00Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
    • B60W2050/0001Details of the control system
    • B60W2050/0019Control system elements or transfer functions
    • B60W2050/0028Mathematical models, e.g. for simulation
    • B60W2050/0031Mathematical model of the vehicle
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W2520/00Input parameters relating to overall vehicle dynamics
    • B60W2520/12Lateral speed
    • 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

Abstract

本发明公开了一种基于分布式鲁棒模型预测的车辆队列横向控制方法,其包括以下步骤:S1、确定车辆的参数;S2、基于确定的车辆参数,建立队列车辆的横向动力学模型;S3、基于步骤S2建立的队列车辆的横向动力学模型,设计模型预测控制器用于车辆系统的控制。本发明考虑了外部道路环境对队列车辆的干扰,通过鲁棒模型预测的方式保证系统实际状态的稳定性。通过基于路径的跟随方式,减小一定范围内的通信时延给车辆状态的影响。

Description

一种基于分布式鲁棒模型预测的车辆队列横向控制方法
技术领域
本发明涉及智慧交通的技术领域,尤其涉及到一种基于分布式鲁棒模型预测的车辆队列横向控制方法。
背景技术
随着汽车保有量的不断增加,交通网络已经接近其所能承载的能力极限。通过增加基础设施建设来增加交通网络容量的方法,无法作为一种能够长期使用的手段。因此,能够提高道路利用率的自动车辆队列驾驶技术,在过去的几十年中得到了发展。车辆队列由一系列能够近距离驾驶的车辆构成,通过车载测量元件或者是车辆间的通信来协调车辆的运动状态。允许近距离驾驶的车辆队列能够增大交通网络的密度,同时也能够减少燃料的消耗。
目前车辆队列研究大多关注于车辆的纵向运动控制,然而在日常行驶场景中,除了直线行驶,变道、拐弯等转向运动也极为常见。在这类驾驶场景中,车辆的转向操作通常交由驾驶员进行,但是在高速低间距的条件下,人为的转向操作将可能导致安全问题,因此对车辆队列进行横向控制显得尤为重要。其中,车辆的横向控制过程涉及车辆的质心侧偏角,如果该角过大,车辆将有侧倾的风险,因此控制策略需保证车辆的横向稳定性。同时,由于队列车辆之间具有一定的耦合,前方车辆的状态将影响到后方的车辆,为避免车辆与其它车道的车辆发生碰撞的风险,控制策略需保证车辆队列的横向队列稳定性。实际的道路环境复杂多样,因此控制策略必须具备一定的鲁棒性。
发明内容
本发明专利的目的是在扰动存在的情况下,提供一种基于分布式鲁棒模型预测的车辆队列横向控制方法,以满足车辆的横向稳定性及队列的稳定性要求。
为实现上述目的,本发明所提供的技术方案为:
一种基于分布式鲁棒模型预测的车辆队列横向控制方法,包括以下步骤:
S1、确定车辆的参数;
S2、基于确定的车辆参数,建立队列车辆的横向动力学模型;
S3、基于步骤S2建立的队列车辆的横向动力学模型,设计模型预测控制器用于车辆系统的控制。
进一步地,所述步骤S1中,车辆的参数包括固定参数和实时参数;
其中,
固有参数包括:车辆前后轴到车辆质心的距离Lf和Lr、车辆的质量m、车辆绕z轴的惯性矩Iz、车辆前轮和后轮的侧偏刚度系数Cαf和Cαr、系统采样时间τ、车辆受到的外界扰动wi(k)的范围
Figure BDA0003698403920000021
车辆转向角的上界ui,max、车辆转向角的下界ui,min、车辆横摆角速度的上界
Figure BDA0003698403920000022
车辆横摆角速度的下界
Figure BDA0003698403920000023
可接受的最大横向速度vy,i,max、可接受的最小横向速度vy,i,min、最大偏航角ψe,i,max、最小偏航角ψe,i,min、最大横向位移误差ye,i,max、最小横向位移误差ye,i,min、系统状态的权重矩阵Pi、控制输入的权重矩阵Qi以及反馈控制率Ki
实时参数包括:通过车载速度传感器和陀螺仪获得的车辆横向速度vy以及横摆角速度
Figure BDA0003698403920000024
通过车载处理器获得的车辆与预期路径的横向位移误差ye和偏航角ψe
进一步地,所述步骤S2包括:
以车辆队列的第一辆车作为领航车辆,编号0,队列中的其他车辆作为跟随车辆;面对复杂的路况,领航车交由人类驾驶,而跟随车辆则通过控制器实现轨迹跟踪;假设车辆在平滑的路面行驶,忽略车辆的俯仰运动,以二自由度的车辆动力学模型作为控制器设计的依据;车辆以速度v行驶,相较于自身纵轴和横轴的速度分别为vx和vy;车辆的前轮转向角为δ,横摆角为ψ,质心侧偏角为β;Fxf和Fxr分别表示轮胎前轮和后轮上的纵向力,Fyf和Fyr分别表示轮胎前轮和后轮上的横向力;则单个车辆的横向动力学模型如下:
Figure BDA0003698403920000031
Figure BDA0003698403920000032
Figure BDA0003698403920000033
其中,
Figure BDA0003698403920000034
为纵轴加速度,
Figure BDA0003698403920000035
为横轴加速度,
Figure BDA0003698403920000036
为横摆角加速度;
假设车辆在正常的道路条件下行驶,轮胎工作在侧偏特性的线性区域内,此时轮胎上的横向受力与其侧滑角度α的关系为:
Fyf=Cαfαf
Fyr=Cαrαr
车辆前后轮的侧滑角α和车辆的质心侧偏角β之间满足:
Figure BDA0003698403920000037
Figure BDA0003698403920000038
αf为车辆前轮的侧滑角,αr为车辆后轮的侧滑角;
车辆的质心侧偏角β表示为:
Figure BDA0003698403920000039
车辆的质心侧偏角小,进行近似,有:
Figure BDA00036984039200000310
Figure BDA00036984039200000311
假设车辆的纵向速度由其纵向控制器保持恒定值,
Figure BDA00036984039200000312
车辆的横向动力学模型表示为:
Figure BDA0003698403920000041
Figure BDA0003698403920000042
上述模型描述单个车辆的横向运动状态,而在车辆行驶过程中,还要求车辆能够跟踪预期轨迹;路径K为车辆的预期路径,车辆的实际位置与预期路径K之间存在偏差;
定义车辆的偏航角ψe为车辆方向与最近路径点S的切线方向之间的夹角:
ψe=ψ-θs
则偏航角速度
Figure BDA0003698403920000043
为:
Figure BDA0003698403920000044
定义车辆的横向位移误差ye为车辆的实际位置与期望位置之间的误差,通过矢量分解可知,横向位移误差的变化率为:
Figure BDA0003698403920000045
车辆的偏航角ψe小,进行小角度近似有:
Figure BDA0003698403920000046
从而得到单车横向跟随模型为:
Figure BDA0003698403920000047
Figure BDA0003698403920000048
Figure BDA0003698403920000049
Figure BDA00036984039200000410
在基于路径的车辆跟随方式下,车辆通过传感器监测前方车辆的位置,并通过车辆间的通信获得在前车在相应位置的横摆角速度信号,以此生成前方车辆的行驶轨迹;车辆通过对路径的跟踪实现对前方车辆的跟随;第i辆车的跟随目标是其当前位置与第i-1车的行驶轨迹最近的点
Figure BDA00036984039200000411
此时车辆将前车在
Figure BDA00036984039200000412
处的横摆角ψi-1作为预期路径的切线方向与水平线之间的夹角,并将前车在s处的横摆角速度信号
Figure BDA0003698403920000051
作为
Figure BDA0003698403920000052
信号,输入到模型中;
以s0作为行驶的起点,第i辆车自s0点起行驶的路程si满足:
Figure BDA0003698403920000053
其中,
Figure BDA0003698403920000054
在行驶过程中vx,i>>vy,i,有
Figure BDA0003698403920000055
由vx,i(t)=c,c为任意常数,得:
si(t)=ct
由于车辆的偏航角小,因此车辆在当前时刻的参考信号近似为前方车辆自相同起点下,相同行驶路程的横摆角速度信号;
假定该信号与当前时刻存在Δt的时间间隔,则:
Figure BDA0003698403920000056
车辆队列的领航车在给定转向角输入后生成行驶轨迹,控制器控制跟随车辆实现队列行驶;取
Figure BDA0003698403920000057
ui(t)=δi(t),
Figure BDA00036984039200000511
假设对任意t<0,都有xi(t)=0,则第i辆跟随车辆的系统状态表述为:
Figure BDA0003698403920000059
其中,
Figure BDA00036984039200000510
由于系统的采样时间为τ,因此第i辆车的系统离散后的结果为:
xi(k+1)=Gixi(k)+Hiui(k)+Iidi(k)
其中,
Figure BDA0003698403920000061
在实际应用场景中,车辆将不可避免的受到外界扰动的影响;假设外界扰动w(k)有界,则加入扰动后第i辆车的模型为:
xi(k+1)=Gixi(k)+Hiui(k)+Iidi(k)+wi(k)
受车辆自身的限制,车辆的转向角有界,定义ui,max为车辆转向角的上界,ui,min为车辆转向角的下界,有:
ui,min≤ui≤ui,max
上述约束构成控制辆的约束集合
Figure BDA0003698403920000062
在车辆队列行驶过程中,为确保车辆行驶不发生侧倾,质心侧偏角接近于零;因此,给定横向速度vy的上下界,以确保行驶的安全性;同时,为确保车辆的跟随性能,给定横向误差与偏航角的上下界;受路面条件限制,车辆的横摆角速度存在上下界;定义
Figure BDA0003698403920000063
为系统状态的上界,
Figure BDA0003698403920000064
为系统状态的下界,有:
xi,min≤xi≤xi,max
上述约束构成状态量的约束集合
Figure BDA0003698403920000065
进一步地,通过模型预测控制器对车辆系统进行控制时,还包括用到标称系统、误差系统以及反馈控制器;其中,标称系统和误差系统皆源于车辆系统;模型预测控制器对去扰动作用的标称系统进行控制;反馈控制器对误差系统进行控制。
进一步地,对车辆系统进行控制的过程如下:
跟随车辆的标称系统不受外界扰动的影响,其表达式如下:
Figure BDA0003698403920000066
其中,
Figure BDA0003698403920000071
为的标称状态,
Figure BDA0003698403920000072
为标称系统下的控制输入,取误差
Figure BDA0003698403920000073
则误差系统表示为:
Figure BDA0003698403920000074
误差系统输出
Figure BDA0003698403920000075
满足:
Figure BDA0003698403920000076
为确保系统状态在扰动作用下的收敛性,误差系统采用反馈控制:
Figure BDA0003698403920000077
此时,误差系统转化为:
Figure BDA0003698403920000078
Figure BDA0003698403920000079
则标称系统的输入和状态约束为:
Figure BDA00036984039200000710
其中
Figure BDA00036984039200000711
为标称系统的最小鲁棒正不变集:
Figure BDA00036984039200000712
其中α满足α∈[0,1);α与p间满足
Figure BDA00036984039200000713
通过模型预测控制器控制标称系统,令状态序列
Figure BDA00036984039200000714
控制序列
Figure BDA00036984039200000715
其中,
Figure BDA00036984039200000716
表示k时刻预测的k+j时刻的系统状态,
Figure BDA00036984039200000717
表示k时刻预测的k+j时刻的系统输入;
模型预测控制器的优化目标如下:
Figure BDA00036984039200000718
满足如下动力学方程及约束:
Figure BDA00036984039200000719
Figure BDA0003698403920000081
Figure BDA0003698403920000082
Figure BDA0003698403920000083
是使得车辆队列横向弦稳定性得到满足时第i辆车的标称系统状态约束集合:
Figure BDA0003698403920000084
其中对于任意i∈[2,N],定义
Figure BDA0003698403920000085
属于集合
Figure BDA0003698403920000086
Figure BDA0003698403920000087
满足不等式:
Figure BDA0003698403920000088
k*与时间间隔Δt间满足
Figure BDA0003698403920000089
β为
Figure BDA00036984039200000810
类函数;
对于第1辆跟随车辆,其状态约束为:
Figure BDA00036984039200000811
集合
Figure BDA00036984039200000812
Figure BDA00036984039200000813
满足:
Figure BDA00036984039200000814
式中σ为
Figure BDA00036984039200000815
类函数;优化目标中每一步的代价函数l以及终端约束Vf满足:
Figure BDA00036984039200000816
Figure BDA00036984039200000817
上式中,Pi和Qi分别为优化目标中预测状态和输入相对应的加权矩阵,
Figure BDA00036984039200000818
为保证控制器稳定性的终端约束;模型预测控制器在每一步时求解优化问题,获得相应的控制序列,其控制输出
Figure BDA00036984039200000819
为:
Figure BDA00036984039200000820
通过控制输出
Figure BDA00036984039200000821
进行车辆系统控制。
与现有技术相比,本方案原理及优点如下:
1)本方案考虑了外部道路环境对队列车辆的干扰,通过鲁棒模型预测的方式保证系统实际状态的稳定性。
2)本方案通过基于路径的跟随方式,减小一定范围内的通信时延给车辆状态的影响。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的服务作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一种基于分布式鲁棒模型预测的车辆队列横向控制方法的原理流程图;
图2为二自由度车辆动力学模型;
图3为车辆跟随预期路径示意图;
图4为车辆跟随示意图;
图5为包括有模型预测控制器的控制框架图;
图6为领航车行驶轨迹示意图;
图7为领航车横摆角速度示意图;
图8为跟随车辆状态及行驶轨迹示意图。
具体实施方式
下面结合具体实施例对本发明作进一步说明:
如图1所示,本实施例所述的一种基于分布式鲁棒模型预测的车辆队列横向控制方法,包括以下步骤:
S1、确定车辆的参数;
本步骤中,车辆参数包括固定参数和实时参数;
其中,
固有参数包括:车辆前后轴到车辆质心的距离Lf和Lr、车辆的质量m、车辆绕z轴的惯性矩Iz、车辆前轮和后轮的侧偏刚度系数Cαf和Cαr、系统采样时间τ、车辆受到的外界扰动wi(k)的范围
Figure BDA0003698403920000101
车辆转向角的上界ui,max、车辆转向角的下界ui,min、车辆横摆角速度的上界
Figure BDA0003698403920000102
车辆横摆角速度的下界
Figure BDA0003698403920000103
可接受的最大横向速度vy,i,max、可接受的最小横向速度vy,i,min、最大偏航角ψe,i,max、最小偏航角ψe,i,min、最大横向位移误差ye,i,max、最小横向位移误差ye,i,min、系统状态的权重矩阵Pi、控制输入的权重矩阵Qi以及反馈控制率Ki
实时参数包括:通过车载速度传感器和陀螺仪获得的车辆横向速度vy以及横摆角速度
Figure BDA0003698403920000104
通过车载处理器获得的车辆与预期路径的横向位移误差ye和偏航角ψe
S2、基于确定的车辆参数,建立队列车辆的横向动力学模型;具体包括:
以车辆队列的第一辆车作为领航车辆,编号0,队列中的其他车辆作为跟随车辆;面对复杂的路况,领航车交由人类驾驶,而跟随车辆则通过控制器实现轨迹跟踪;假设车辆在平滑的路面行驶,忽略车辆的俯仰运动,以二自由度的车辆动力学模型作为控制器设计的依据,车辆的前后轮分别由单个轮胎代替,如图2所示;车辆以速度v行驶,相较于自身纵轴和横轴的速度分别为vx和vy;车辆的前轮转向角为δ,横摆角为ψ,质心侧偏角为β;Fxf和Fxr分别表示轮胎前轮和后轮上的纵向力,Fyf和Fyr分别表示轮胎前轮和后轮上的横向力;则单个车辆的横向动力学模型如下:
Figure BDA0003698403920000105
Figure BDA0003698403920000106
Figure BDA0003698403920000107
其中,
Figure BDA0003698403920000111
为纵轴加速度,
Figure BDA0003698403920000112
为横轴加速度,
Figure BDA0003698403920000113
为横摆角加速度;
假设车辆在正常的道路条件下行驶,轮胎工作在侧偏特性的线性区域内,此时轮胎上的横向受力与其侧滑角度α的关系为:
Fyf=Cαfαf
Fyr=Cαrαr
车辆前后轮的侧滑角α和车辆的质心侧偏角β之间满足:
Figure BDA0003698403920000114
Figure BDA0003698403920000115
αf为车辆前轮的侧滑角,αr为车辆后轮的侧滑角;
车辆的质心侧偏角β表示为:
Figure BDA0003698403920000116
车辆的质心侧偏角小,进行近似,有:
Figure BDA0003698403920000117
Figure BDA0003698403920000118
假设车辆的纵向速度由其纵向控制器保持恒定值,
Figure BDA0003698403920000119
车辆的横向动力学模型表示为:
Figure BDA00036984039200001110
Figure BDA00036984039200001111
上述模型描述单个车辆的横向运动状态,而在车辆行驶过程中,还要求车辆能够跟踪预期轨迹,如图3所示;路径K为车辆的预期路径,车辆的实际位置与预期路径K之间存在偏差;
定义车辆的偏航角ψe为车辆方向与最近路径点S的切线方向之间的夹角:
ψe=ψ-θs
则偏航角速度
Figure BDA0003698403920000121
为:
Figure BDA0003698403920000122
定义车辆的横向位移误差ye为车辆的实际位置与期望位置之间的误差,通过矢量分解可知,横向位移误差的变化率为:
Figure BDA0003698403920000123
车辆的偏航角ψe小,进行小角度近似有:
Figure BDA0003698403920000124
从而得到单车横向跟随模型为:
Figure BDA0003698403920000125
Figure BDA0003698403920000126
Figure BDA0003698403920000127
Figure BDA0003698403920000128
在基于路径的车辆跟随方式下,车辆通过传感器监测前方车辆的位置,并通过车辆间的通信获得在前车在相应位置的横摆角速度信号,以此生成前方车辆的行驶轨迹;车辆通过对路径的跟踪实现对前方车辆的跟随,如图4所示;第i辆车的跟随目标是其当前位置与第i-1车的行驶轨迹最近的点
Figure BDA0003698403920000129
此时车辆将前车在
Figure BDA00036984039200001210
处的横摆角ψi-1作为预期路径的切线方向与水平线之间的夹角,并将前车在
Figure BDA00036984039200001211
处的横摆角速度信号
Figure BDA00036984039200001212
作为
Figure BDA00036984039200001213
信号,输入到模型中;
以s0作为行驶的起点,第i辆车自s0点起行驶的路程si满足:
Figure BDA00036984039200001214
其中,
Figure BDA0003698403920000131
在行驶过程中vx,i>>vy,i,有
Figure BDA0003698403920000132
由vx,i(t)=c,c为任意常数,得:
si(t)=ct
由于车辆的偏航角小,因此车辆在当前时刻的参考信号近似为前方车辆自相同起点下,相同行驶路程的横摆角速度信号;
假定该信号与当前时刻存在Δt的时间间隔,则:
Figure BDA0003698403920000133
车辆队列的领航车在给定转向角输入后生成行驶轨迹,控制器控制跟随车辆实现队列行驶;取
Figure BDA0003698403920000134
ui(t)=δi(t),
Figure BDA0003698403920000139
假设对任意t<0,都有xi(t)=0,则第i辆跟随车辆的系统状态表述为:
Figure BDA0003698403920000136
其中,
Figure BDA0003698403920000137
由于系统的采样时间为τ,因此第i辆车的系统离散后的结果为:
xi(k+1)=Gixi(k)+Hiui(k)+Iidi(k)
其中,
Figure BDA0003698403920000138
在实际应用场景中,车辆将不可避免的受到外界扰动的影响;假设外界扰动w(k)有界,则加入扰动后第i辆车的模型为:
xi(k+1)=Gixi(k)+Hiui(k)+Iidi(k)+wi(k)
受车辆自身的限制,车辆的转向角有界,定义ui,max为车辆转向角的上界,ui,min为车辆转向角的下界,有:
ui,min≤ui≤ui,max
上述约束构成控制辆的约束集合
Figure BDA0003698403920000141
在车辆队列行驶过程中,为确保车辆行驶不发生侧倾,质心侧偏角接近于零;因此,给定横向速度vy的上下界,以确保行驶的安全性;同时,为确保车辆的跟随性能,给定横向误差与偏航角的上下界;受路面条件限制,车辆的横摆角速度存在上下界;定义
Figure BDA0003698403920000142
为系统状态的上界,
Figure BDA0003698403920000143
为系统状态的下界,有:
xi,min≤xi≤xi,max
上述约束构成状态量的约束集合
Figure BDA0003698403920000144
S3、基于步骤S2建立的队列车辆的横向动力学模型,设计模型预测控制器用于车辆系统的控制。
本步骤中,通过模型预测控制器对车辆系统进行控制时,还包括用到标称系统、误差系统以及反馈控制器;其中,标称系统和误差系统皆源于车辆系统;模型预测控制器对去扰动作用的标称系统进行控制;反馈控制器对误差系统进行控制,如图5所示。
具体控制过程如下:
跟随车辆的标称系统不受外界扰动的影响,其表达式如下:
Figure BDA0003698403920000145
其中,
Figure BDA0003698403920000146
为的标称状态,
Figure BDA0003698403920000147
为标称系统下的控制输入,取误差
Figure BDA0003698403920000148
则误差系统表示为:
Figure BDA0003698403920000149
误差系统输出
Figure BDA0003698403920000151
满足:
Figure BDA0003698403920000152
为确保系统状态在扰动作用下的收敛性,误差系统采用反馈控制:
Figure BDA0003698403920000153
此时,误差系统转化为:
Figure BDA0003698403920000154
Figure BDA0003698403920000155
则标称系统的输入和状态约束为:
Figure BDA0003698403920000156
其中
Figure BDA0003698403920000157
为标称系统的最小鲁棒正不变集:
Figure BDA0003698403920000158
其中α满足α∈[0,1);α与p间满足
Figure BDA0003698403920000159
通过模型预测控制器控制标称系统,令状态序列
Figure BDA00036984039200001510
控制序列
Figure BDA00036984039200001511
其中,
Figure BDA00036984039200001512
表示k时刻预测的k+j时刻的系统状态,
Figure BDA00036984039200001513
表示k时刻预测的k+j时刻的系统输入;
模型预测控制器的优化目标如下:
Figure BDA00036984039200001514
满足如下动力学方程及约束:
Figure BDA00036984039200001515
Figure BDA00036984039200001516
Figure BDA00036984039200001517
Figure BDA00036984039200001518
是使得车辆队列横向弦稳定性得到满足时第i辆车的标称系统状态约束集合:
Figure BDA0003698403920000161
其中对于任意i∈[2,N],定义
Figure BDA0003698403920000162
属于集合
Figure BDA0003698403920000163
Figure BDA0003698403920000164
满足不等式:
Figure BDA0003698403920000165
k*与时间间隔Δt间满足
Figure BDA0003698403920000166
β为
Figure BDA0003698403920000167
类函数;
对于第1辆跟随车辆,其状态约束为:
Figure BDA0003698403920000168
集合
Figure BDA0003698403920000169
Figure BDA00036984039200001610
满足:
Figure BDA00036984039200001611
式中σ为
Figure BDA00036984039200001612
类函数;优化目标中每一步的代价函数l以及终端约束Vf满足:
Figure BDA00036984039200001613
Figure BDA00036984039200001614
上式中,Pi和Qi分别为优化目标中预测状态和输入相对应的加权矩阵,
Figure BDA00036984039200001615
为保证控制器稳定性的终端约束;模型预测控制器在每一步时求解优化问题,获得相应的控制序列,其控制输出
Figure BDA00036984039200001616
为:
Figure BDA00036984039200001617
通过控制输出
Figure BDA00036984039200001618
进行车辆系统控制。
为证实本发明的有效性,作如下仿真:
假设车辆队列由6辆车构成,车辆参数如下:Lr=1.6m,Lf=1.1m,Cαr=143000N/rad,Cαf=117000N/rad,m=1650kg,Iz=2900kgm2,vx=72km/h。领航车的初始位置(X0,Y0)=(20,0),基于路径的跟随方式下的时间间隔Δt=1s,连续系统转化为离散系统的采样间隔τ=0.02s,车辆间通信时延τ*=0.1s。预测步长Np=40,扰动wi为区间[-0.00002,0.00002]内均匀分布的随机噪声。
领航车行驶轨迹如图6所示,领航车横摆角速度如图7所示,模型预测控制器作用下的跟随车辆状态及行驶轨迹如图8所示。
以上所述之实施例子只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。

Claims (5)

1.一种基于分布式鲁棒模型预测的车辆队列横向控制方法,其特征在于,包括以下步骤:
S1、确定车辆的参数;
S2、基于确定的车辆参数,建立队列车辆的横向动力学模型;
S3、基于步骤S2建立的队列车辆的横向动力学模型,设计模型预测控制器用于车辆系统的控制。
2.根据权利要求1所述的一种基于分布式鲁棒模型预测的车辆队列横向控制方法,其特征在于,所述步骤S1中,车辆的参数包括固定参数和实时参数;
其中,
固有参数包括:车辆前后轴到车辆质心的距离Lf和Lr、车辆的质量m、车辆绕z轴的惯性矩Iz、车辆前轮和后轮的侧偏刚度系数Cαf和Cαr、系统采样时间τ、车辆受到的外界扰动wi(k)的范围
Figure FDA0003698403910000011
车辆转向角的上界ui,max、车辆转向角的下界ui,min、车辆横摆角速度的上界
Figure FDA0003698403910000012
车辆横摆角速度的下界
Figure FDA0003698403910000013
可接受的最大横向速度vy,i,max、可接受的最小横向速度vy,i,min、最大偏航角ψe,i,max、最小偏航角ψe,i,min、最大横向位移误差ye,i,max、最小横向位移误差ye,i,min、系统状态的权重矩阵Pi、控制输入的权重矩阵Qi以及反馈控制率Ki
实时参数包括:通过车载速度传感器和陀螺仪获得的车辆横向速度vy以及横摆角速度
Figure FDA0003698403910000014
通过车载处理器获得的车辆与预期路径的横向位移误差ye和偏航角ψe
3.根据权利要求2所述的一种基于分布式鲁棒模型预测的车辆队列横向控制方法,其特征在于,所述步骤S2包括:
以车辆队列的第一辆车作为领航车辆,编号0,队列中的其他车辆作为跟随车辆;面对复杂的路况,领航车交由人类驾驶,而跟随车辆则通过控制器实现轨迹跟踪;假设车辆在平滑的路面行驶,忽略车辆的俯仰运动,以二自由度的车辆动力学模型作为控制器设计的依据;车辆以速度v行驶,相较于自身纵轴和横轴的速度分别为vx和vy;车辆的前轮转向角为δ,横摆角为ψ,质心侧偏角为β;Fxf和Fxr分别表示轮胎前轮和后轮上的纵向力,Fyf和Fyr分别表示轮胎前轮和后轮上的横向力;则单个车辆的横向动力学模型如下:
Figure FDA0003698403910000021
Figure FDA0003698403910000022
Figure FDA0003698403910000023
其中,
Figure FDA0003698403910000024
为纵轴加速度,
Figure FDA0003698403910000025
为横轴加速度,
Figure FDA0003698403910000026
为横摆角加速度;
假设车辆在正常的道路条件下行驶,轮胎工作在侧偏特性的线性区域内,此时轮胎上的横向受力与其侧滑角度α的关系为:
Fyf=Cαfαf
Fyr=Cαrαr
车辆前后轮的侧滑角α和车辆的质心侧偏角β之间满足:
Figure FDA0003698403910000027
Figure FDA0003698403910000028
αf为车辆前轮的侧滑角,αr为车辆后轮的侧滑角;
车辆的质心侧偏角β表示为:
Figure FDA0003698403910000029
车辆的质心侧偏角小,进行近似,有:
Figure FDA0003698403910000031
Figure FDA0003698403910000032
假设车辆的纵向速度由其纵向控制器保持恒定值,
Figure FDA0003698403910000033
车辆的横向动力学模型表示为:
Figure FDA0003698403910000034
Figure FDA0003698403910000035
所述模型描述单个车辆的横向运动状态,而在车辆行驶过程中,还要求车辆能够跟踪预期轨迹;路径K为车辆的预期路径,车辆的实际位置与预期路径K之间存在偏差;
定义车辆的偏航角ψe为车辆方向与最近路径点S的切线方向之间的夹角:
ψe=ψ-θs
则偏航角速度
Figure FDA0003698403910000036
为:
Figure FDA0003698403910000037
定义车辆的横向位移误差ye为车辆的实际位置与期望位置之间的误差,通过矢量分解可知,横向位移误差的变化率为:
Figure FDA0003698403910000038
车辆的偏航角ψe小,进行小角度近似有:
Figure FDA0003698403910000039
从而得到单车横向跟随模型为:
Figure FDA0003698403910000041
Figure FDA0003698403910000042
Figure FDA0003698403910000043
Figure FDA0003698403910000044
在基于路径的车辆跟随方式下,车辆通过传感器监测前方车辆的位置,并通过车辆间的通信获得在前车在相应位置的横摆角速度信号,以此生成前方车辆的行驶轨迹;车辆通过对路径的跟踪实现对前方车辆的跟随;第i辆车的跟随目标是其当前位置与第i-1车的行驶轨迹最近的点
Figure FDA0003698403910000045
此时车辆将前车在
Figure FDA0003698403910000046
处的横摆角ψi-1作为预期路径的切线方向与水平线之间的夹角,并将前车在
Figure FDA0003698403910000047
处的横摆角速度信号
Figure FDA0003698403910000048
作为
Figure FDA0003698403910000049
信号,输入到模型中;
以s0作为行驶的起点,第i辆车自s0点起行驶的路程si满足:
Figure FDA00036984039100000410
其中,
Figure FDA00036984039100000411
在行驶过程中vx,i>>vy,i,有
Figure FDA00036984039100000412
由vx,i(t)=c,c为任意常数,得:
si(t)=ct
由于车辆的偏航角小,因此车辆在当前时刻的参考信号近似为前方车辆自相同起点下,相同行驶路程的横摆角速度信号;
假定该信号与当前时刻存在Δt的时间间隔,则:
Figure FDA00036984039100000413
车辆队列的领航车在给定转向角输入后生成行驶轨迹,控制器控制跟随车辆实现队列行驶;取
Figure FDA0003698403910000051
ui(t)=δi(t),
Figure FDA0003698403910000052
假设对任意t<0,都有xi(t)=0,则第i辆跟随车辆的系统状态表述为:
Figure FDA0003698403910000053
其中,
Figure FDA0003698403910000054
由于系统的采样时间为τ,因此第i辆车的系统离散后的结果为:
xi(k+1)=Gixi(k)+Hiui(k)+Iidi(k)
其中,
Figure FDA0003698403910000055
在实际应用场景中,车辆将不可避免的受到外界扰动的影响;假设外界扰动w(k)有界,则加入扰动后第i辆车的模型为:
xi(k+1)=Gixi(k)+Hiui(k)+Iidi(k)+wi(k)
受车辆自身的限制,车辆的转向角有界,定义ui,max为车辆转向角的上界,ui,min为车辆转向角的下界,有:
ui,min≤ui≤ui,max
所述约束构成控制辆的约束集合
Figure FDA0003698403910000056
在车辆队列行驶过程中,为确保车辆行驶不发生侧倾,质心侧偏角接近于零;因此,给定横向速度vy的上下界,以确保行驶的安全性;同时,为确保车辆的跟随性能,给定横向误差与偏航角的上下界;受路面条件限制,车辆的横摆角速度存在上下界;定义
Figure FDA0003698403910000057
为系统状态的上界,
Figure FDA0003698403910000058
为系统状态的下界,有:
xi,min≤xi≤xi,max
上述约束构成状态量的约束集合
Figure FDA0003698403910000061
4.根据权利要求3所述的一种基于分布式鲁棒模型预测的车辆队列横向控制方法,其特征在于,通过模型预测控制器对车辆系统进行控制时,还包括用到标称系统、误差系统以及反馈控制器;其中,标称系统和误差系统皆源于车辆系统;模型预测控制器对去扰动作用的标称系统进行控制;反馈控制器对误差系统进行控制。
5.根据权利要求4所述的一种基于分布式鲁棒模型预测的车辆队列横向控制方法,其特征在于,对车辆系统进行控制的过程如下:
跟随车辆的标称系统不受外界扰动的影响,其表达式如下:
Figure FDA0003698403910000062
其中,
Figure FDA0003698403910000063
为的标称状态,
Figure FDA0003698403910000064
为标称系统下的控制输入,取误差
Figure FDA0003698403910000065
则误差系统表示为:
Figure FDA0003698403910000066
误差系统输出
Figure FDA0003698403910000067
满足:
Figure FDA0003698403910000068
为确保系统状态在扰动作用下的收敛性,误差系统采用反馈控制:
Figure FDA0003698403910000069
此时,误差系统转化为:
Figure FDA00036984039100000610
Figure FDA00036984039100000611
则标称系统的输入和状态约束为:
Figure FDA00036984039100000612
其中
Figure FDA00036984039100000613
为标称系统的最小鲁棒正不变集:
Figure FDA0003698403910000071
其中α满足α∈[0,1);α与p间满足
Figure FDA0003698403910000072
通过模型预测控制器控制标称系统,令状态序列
Figure FDA0003698403910000073
控制序列
Figure FDA0003698403910000074
其中,
Figure FDA0003698403910000075
表示k时刻预测的k+j时刻的系统状态,
Figure FDA0003698403910000076
表示k时刻预测的k+j时刻的系统输入;
模型预测控制器的优化目标如下:
Figure FDA0003698403910000077
满足如下动力学方程及约束:
Figure FDA0003698403910000078
Figure FDA0003698403910000079
Figure FDA00036984039100000710
Figure FDA00036984039100000711
是使得车辆队列横向弦稳定性得到满足时第i辆车的标称系统状态约束集合:
Figure FDA00036984039100000712
其中对于任意i∈[2,N],定义
Figure FDA00036984039100000713
属于集合
Figure FDA00036984039100000714
Figure FDA00036984039100000715
满足不等式:
Figure FDA00036984039100000716
k*与时间间隔Δt间满足
Figure FDA00036984039100000717
β为
Figure FDA00036984039100000718
类函数;
对于第1辆跟随车辆,其状态约束为:
Figure FDA00036984039100000719
集合
Figure FDA00036984039100000720
满足:
Figure FDA0003698403910000081
式中σ为
Figure FDA0003698403910000082
类函数;优化目标中每一步的代价函数l以及终端约束Vf满足:
Figure FDA0003698403910000083
Figure FDA0003698403910000084
上式中,Pi和Qi分别为优化目标中预测状态和输入相对应的加权矩阵,
Figure FDA0003698403910000085
为保证控制器稳定性的终端约束;模型预测控制器在每一步时求解优化问题,获得相应的控制序列,其控制输出
Figure FDA0003698403910000086
为:
Figure FDA0003698403910000087
通过控制输出
Figure FDA0003698403910000088
进行车辆系统控制。
CN202210681057.7A 2022-06-16 2022-06-16 一种基于分布式鲁棒模型预测的车辆队列横向控制方法 Active CN115042770B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210681057.7A CN115042770B (zh) 2022-06-16 2022-06-16 一种基于分布式鲁棒模型预测的车辆队列横向控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210681057.7A CN115042770B (zh) 2022-06-16 2022-06-16 一种基于分布式鲁棒模型预测的车辆队列横向控制方法

Publications (2)

Publication Number Publication Date
CN115042770A true CN115042770A (zh) 2022-09-13
CN115042770B CN115042770B (zh) 2024-04-16

Family

ID=83161128

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210681057.7A Active CN115042770B (zh) 2022-06-16 2022-06-16 一种基于分布式鲁棒模型预测的车辆队列横向控制方法

Country Status (1)

Country Link
CN (1) CN115042770B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116729407A (zh) * 2023-08-15 2023-09-12 蘑菇车联信息科技有限公司 车辆横向速度检测方法、电子设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040059480A1 (en) * 2002-08-01 2004-03-25 Hrovat Davorin David System and method for detecting roll rate sensor fault
US20060085112A1 (en) * 2004-10-15 2006-04-20 Ford Global Technologies, Llc System and method for dynamically determining vehicle loading and vertical loading distance for use in a vehicle dynamic control system
US20170329348A1 (en) * 2016-05-13 2017-11-16 Macau University Of Science And Technology System and method for operating a follower vehicle in a vehicle platoon
CN113848942A (zh) * 2021-10-18 2021-12-28 青岛德智汽车科技有限公司 一种约束导向的智能网联汽车鲁棒换道合流控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040059480A1 (en) * 2002-08-01 2004-03-25 Hrovat Davorin David System and method for detecting roll rate sensor fault
US20060085112A1 (en) * 2004-10-15 2006-04-20 Ford Global Technologies, Llc System and method for dynamically determining vehicle loading and vertical loading distance for use in a vehicle dynamic control system
US20170329348A1 (en) * 2016-05-13 2017-11-16 Macau University Of Science And Technology System and method for operating a follower vehicle in a vehicle platoon
CN113848942A (zh) * 2021-10-18 2021-12-28 青岛德智汽车科技有限公司 一种约束导向的智能网联汽车鲁棒换道合流控制方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116729407A (zh) * 2023-08-15 2023-09-12 蘑菇车联信息科技有限公司 车辆横向速度检测方法、电子设备及存储介质
CN116729407B (zh) * 2023-08-15 2023-11-28 蘑菇车联信息科技有限公司 车辆横向速度检测方法、电子设备及存储介质

Also Published As

Publication number Publication date
CN115042770B (zh) 2024-04-16

Similar Documents

Publication Publication Date Title
CN110187639B (zh) 一种基于参数决策框架的轨迹规划控制方法
CN109855639B (zh) 基于障碍物预测与mpc算法的无人驾驶轨迹规划方法
CN109669461B (zh) 一种复杂工况下自动驾驶车辆决策系统及其轨迹规划方法
CN111717192B (zh) 一种自动驾驶车辆的控制方法及系统
CN109144076B (zh) 一种多车辆横纵向耦合协同控制系统及控制方法
CN109131312B (zh) 一种智能电动汽车acc/esc集成控制系统及其方法
CN108454623A (zh) 一种四轮独立驱动无人驾驶电动车辆轨迹跟踪控制方法
CN112445229B (zh) 一种领航车队协同的单车道多队列分层控制方法
Hima et al. Trajectory tracking for highly automated passenger vehicles
CN110962849B (zh) 一种弯道自适应巡航方法
CN109969183A (zh) 基于安全可控域的弯道跟车控制方法
CN111409641A (zh) 一种基于纵横向协调的轨迹跟踪控制方法及系统
CN103085816A (zh) 一种用于无人驾驶车辆的轨迹跟踪控制方法及控制装置
CN110703754A (zh) 一种自动驾驶车辆路径与速度高度耦合的轨迹规划方法
CN110827535A (zh) 非线性车辆队列协同自适应抗扰纵向控制方法
CN112068445B (zh) 自动驾驶车辆路径规划与路径跟踪集成控制方法及系统
CN104881030A (zh) 基于快速终端滑模原理的无人车侧纵向耦合跟踪控制方法
CN115042770B (zh) 一种基于分布式鲁棒模型预测的车辆队列横向控制方法
CN110703775A (zh) 一种基于改进人工势场的商用车车道保持路径规划方法
CN112937551B (zh) 一种考虑驾驶员输入特征的车辆控制方法及系统
CN114312848A (zh) 基于双层mpc的智能驾驶汽车轨迹规划与跟踪控制方法
CN108594830A (zh) 一种基于空间域的网联智能车编队行驶控制方法
CN111661048A (zh) 多铰接式车辆及其轨迹跟随控制方法与系统
CN114357631B (zh) 基于协同分布式优化的集成式车辆防侧翻主动控制方法
CN115447615A (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