CN114148318B - 冰雪环境下基于反馈线性化与lqr的车辆路径跟踪方法 - Google Patents
冰雪环境下基于反馈线性化与lqr的车辆路径跟踪方法 Download PDFInfo
- Publication number
- CN114148318B CN114148318B CN202111576809.5A CN202111576809A CN114148318B CN 114148318 B CN114148318 B CN 114148318B CN 202111576809 A CN202111576809 A CN 202111576809A CN 114148318 B CN114148318 B CN 114148318B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- formula
- unit
- model
- nonlinear
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims abstract description 17
- 238000013461 design Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000006073 displacement reaction Methods 0.000 claims description 9
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical group C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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/00—Purposes 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
- B60W30/02—Control of vehicle driving stability
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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/00—Details 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
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Human Computer Interaction (AREA)
- Steering Control In Accordance With Driving Conditions (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
本发明提供了一种冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法,根据冰雪环境下车辆轮胎力不满足线性关系,建立冰雪环境下仿射形式的非线性车辆系统模型,基于此模型,采用反馈线性化的方法对复杂的非线性车辆系统模型进行线性化处理,得到较简单的线性系统模型和虚拟控制输入;根据得到的简单线性车辆系统模型,使用线性二次型调节器的设计方法设计路径跟踪控制器,保证系统稳定性和目标的最优性,实现无人驾驶车辆路径跟踪。
Description
技术领域
本发明涉及一种冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法,属于汽车驾驶控制领域。
背景技术
随着经济的迅速发展,汽车已经普及到每家每户,成为我们日常出行方式中一个重要的组成部分。虽然汽车给我们的生活带来了极大的便利,但是汽车的使用也引发了一系列的问题。我国地域辽阔,环境复杂,南北方的道路环境情况差别巨大。在我国北方很多地区的公路长期受到冰雪天气的危害,一年中有近半年的道路存在积雪,对道路的通行能力以及安全性均会造成严重的影响。冰雪天气会造成路面附着系数下降,道路环境复杂,造成严重的安全隐患。与此同时冰雪天气也会造成车辆轮胎侧向稳定性降低,使车辆在行驶时更容易发生侧滑的现象。所以本发明将关注点放在冰雪环境下无人驾驶技术的研究上。
发明内容
本发明所述的冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法,考虑到冰雪环境下汽车的动力学具有非线性特性,采用反馈线性化的方法将非线性的车辆模型转化为线性的车辆模型,再根据线性二次调节器的方法设计控制器,研究冰雪环境下非线性车辆系统路径跟踪的方法。
本发明是采用如下技术方案实现的:
冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法,根据冰雪环境下车辆轮胎力不满足线性关系的问题,建立冰雪环境下仿射形式的非线性车辆系统模型,基于此模型,采用反馈线性化的方法对复杂的非线性车辆系统模型进行线性化处理,得到较简单的线性系统模型和虚拟控制输入;根据得到的简单线性车辆系统模型,使用线性二次型调节器的设计方法设计路径跟踪控制器,保证系统稳定性和目标的最优性,实现无人驾驶车辆路径跟踪;本方法的具体步骤如下:
步骤一、建立仿射形式的非线性车辆系统模型
假定车辆是一个刚性体,其中车辆装置着四个不会发生形变的车轮,并以车辆的前轮作为转向轮,根据车辆运动的几何关系以及车辆的运动学方程得到车辆运动学模型,如式(1)所示:
其中ψ为车辆横摆角,单位:rad;vo为车辆质心o处的速度,单位:m/s;β为车辆质心侧偏角,单位:rad;xO为车辆质心o的纵向位移,单位:m;yO为车辆质心o的横向位移,单位:m;
考虑到车辆在运行过程中主要是纵向行驶,因此vx>>vy,所以车辆行驶的速度v近似用纵向速度vx表示,同时,车辆的横摆角ψ和质心侧偏角β都较小,因此采用如下近似关系:
将式(2)代入到式(1)中进行简化处理,得到简化的车辆运动学模型:
对于车辆的动力学模型,以车辆质心o为坐标原点,以沿着车身向前的方向为x轴的正方向,垂直于横轴向上的方向为y轴的正方向,通过质心指向上方的方向为z轴的正方向,忽略车辆的纵向动力学,只考虑车辆的横摆方向动力学和侧向动力学,根据牛顿第二定律与力矩平衡关系推出如式(4)所示的单轨模型数学表达式:
其中,Fyf为车辆的前轮侧向力,单位:N;Fyr为车辆的后轮侧向力,单位:N;Fxf为车辆的前轮纵向力,单位:N;Fxr为车辆的后轮纵向力,单位:N;vy为车辆质心处的侧向速度,单位:m/s;vx为车辆质心处的纵向速度,单位:m/s;r为车辆的横摆角速度,单位:rad/s;δ为车辆前轮转角,单位:rad;m为车辆的质量,单位:kg;Iz为车辆绕z轴的转动惯量,单位:kg·m2;a为车辆质心o到车辆前轴的距离,单位:m;b为车辆质心o到车辆后轴的距离,单位:m;
本方法的控制变量为车辆的前轮转角,考虑车辆的侧向动力学特性,所以假设车辆的纵向速度不变,则式(4)可以写成:
行驶稳定的车辆,其前轮转角δ很小,所以采用如下近似关系:
将式(6)代入到式(5)中得:
考虑到冰雪环境下,路面的附着系数低,当车辆行驶在冰雪路面时,轮胎力呈非线性的关系,所以选用如下的轮胎模型来表示轮胎侧向力与轮胎侧偏角的关系:
其中,Cf表示前轮的侧偏刚度,单位:N·rad;Cr表示后轮的侧偏刚度,单位:N·rad;Ka表示前轮的拟合参数;Kb表示后轮的拟合参数;αf表示前轮的侧偏角,单位:rad;αr表示后轮的侧偏角,单位:rad;Cf的取值在28000~90000的范围内,Cr的取值在28000~90000的范围内,Ka的取值大于0,Kb的取值大于0;
为了便于计算,将式(8)写成:
根据坐标系的规定,分别描述车辆前轮的轮胎侧偏角与后轮的轮胎侧偏角为:
将式(10)代入到式(9)中得:
将式(11)代入到式(7)中可得:
其中a1=(aAf-bAr)/(Izvx),a3=(a2Af+b2Ar)/(Izvx),b1=(Af+Ar)/(mvx),b3=(aAf-bAr)/(mvx)-vx,p0=(-aAf)/Iz,q0=(-Af)/m,Δ1≈Δ2=0;
将式(12)和式(13)中的所有vy转换成β得到车辆的动力学模型:
其中c1=(aAf-bAr)/(Iz),c2=(aBf-bBr)/(Iz),c3=(a2Af+b2Ar)/(Izvx),c6=3(a2Bf+b2Br)/(Izvx),d1=(Af+Ar)/(mvx), l0=(-aAf)/Iz,l1=(-3aBf)/Iz,l2=(-6a2Bf)/(Izvx),w0=(-Af)/(mvx),w1=(-3Bf)/(mvx),Δ3≈Δ4=0;
结合式(3)、式(14)和式(15),整理得到符合以下仿射形式的非线性车辆系统模型:
其中f1(x)=vx(x2+x4),f2(x)=x3, h(x)=yo,δ为前轮转角,
步骤二、采用反馈线性化的方法对车辆系统模型线性化
根据单输入单输出反馈线性化的基本原理,为了建立y与δ之间的关系,对式(16)的输出关于时间求导,利用李导数整理得:
Lfh(x)=vx(x2+x4),Lgh(x)=0 (18)
在x=0附近Lgh(x)=0,则再次关于时间求导,整理后得到:
在x=0附近LgLfh(x)≠0,所以车辆系统的相对阶为2;
对式(16)中的仿射形式的非线性车辆系统模型进行状态变换:
定义新的状态量z1=h(x),所以新的状态方程可以表示为:
非线性车辆系统转变为:
所以式(16)可以转化为以下形式:
其中C=[10];
由于(A,B)在x的取值范围内是可控的,γ(x)为非奇异矩阵,因此可对步骤一中所提出的仿射形式的非线性车辆系统模型进行反馈线性化;
设虚拟控制输入为v,并且求出δ与v的关系:
所以,非线性系统可以转换为以下线性系统:
其中C=[1 0];
步骤三、基于反馈线性化与LQR设计路径跟踪控制器
对于线性化后的车辆模型,利用线性二次型调节器的方法设计跟踪控制器,设理想的侧向位移为yd,则理想值与实际值之间的误差为yd-y=e,求最优控制输入v,使下列的二次型性能指标最小:
其中正定对称矩阵R为控制输入变量v的加权矩阵,半正定矩阵Q为状态变量z的加权矩阵,t0为控制时间的起点,tf为控制时间的终点;
当通过反馈线性化后的车辆系统最优控制输入为v=R-1BTPe=-Ke,二次型性能指标最小;
其中e可以表示为:
所以v可以进一步地表示为:
最优状态反馈控制增益矩阵为:
K=[K1 K2]=R-1BTP (30)
其中P是一个正定矩阵满足代数Riccati方程:
PA+ATP+Q-PBR-1BTP=0 (31)
此时δ与v的关系可以进一步地表示为:
将得到的控制量δ关系式输入到车辆系统中,实现无人驾驶车辆的路径跟踪。
与现有技术相比,本发明的有益效果为:
本发明提供了一种冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法,因为对非线性动力学模型进行了化简,所以车辆非线性模型可以整理为仿射非线性的形式。因为本方法建立了仿射形式的非线性车辆系统模型,所以能够采用反馈线性化的方法对复杂的非线性车辆系统模型进行线性化处理,得到较简单的线性系统模型和虚拟控制输入;因为本方法使用反馈线性化的方法得到了较简单的线性车辆系统模型,所以控制器的设计过程可以使用线性二次型调节器的设计方法。因为本方法对控制器的设计过程可以使用线性二次型调节器的设计方法,控制器可以保证系统稳定性和目标的最优性,实现无人驾驶车辆路径跟踪。
附图说明
图1为本发明提供的冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法的流程图;
图2为本发明提供的车辆运动学模型示意图;
图3为本发明提供的单轨模型示意图;
图4为本发明提供的冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法车辆路径跟踪控制框图;
具体实施方式
下面对本发明做详细的描述:
本发明提供了一种冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法,具体步骤过程如下:
步骤一、建立仿射形式的非线性车辆系统模型
描述车辆运行状态的模型按其具体功能可分为运动学模型和动力学模型两类。通过研究表明,在高速行驶时,车辆的动态特性对其自身的运行状态影响更大,在低速行驶时,车辆的运动学特性对其自身的运行状态影响更大;所以,本文同时考虑了车辆的运动学和动力学特性,将无人驾驶汽车的运动视为平面运动。
对于车辆的运动学模性,假定车辆是一个刚性体,其中车辆装置着四个不会发生形变的车轮,并以车辆的前轮作为转向轮,根据车辆运动的几何关系以及车辆的运动学方程得到车辆运动学模型,如式(1)所示:
其中ψ为车辆横摆角,单位:rad;vo为车辆质心o处的速度,单位:m/s;β为车辆质心侧偏角,单位:rad;xO为车辆质心o的纵向位移,单位:m;yO为车辆质心o的横向位移,单位:m;
考虑到车辆在运行过程中主要是纵向行驶,因此vx>>vy,所以车辆行驶的速度v近似用纵向速度vx表示,同时,车辆的横摆角ψ和质心侧偏角β都较小,因此采用如下近似关系:
将式(2)代入到式(1)中进行简化处理,得到简化的车辆运动学模型:
对于车辆的动力学模型,由于汽车的动力学特性非常复杂,如果精确描述车辆的运行状态,则需要建立高自由度的模型。高自由度的模型虽然可以精确描述车辆的运行状态,但是高自由度模型复杂,使用他们不利于进行控制器设计的研究。本发明选用简单的单轨模型。单轨模型忽略了转向系统的影响,直接使用前轮转角作为输入;它也忽略了悬架的作用,空气动力学的作用等。单轨模型可以用一个等效的前轮与一个等效的后轮来代替前轮与后轮。具体步骤为以车辆质心o为坐标原点,以沿着车身向前的方向为x轴的正方向,垂直于横轴向上的方向为y轴的正方向,通过质心指向上方的方向为z轴的正方向,忽略车辆的纵向动力学,只考虑车辆的横摆方向动力学和侧向动力学,根据牛顿第二定律与力矩平衡关系推出如式(4)所示的单轨模型数学表达式:
其中,Fyf为车辆的前轮侧向力,单位:N;Fyr为车辆的后轮侧向力,单位:N;Fxf为车辆的前轮纵向力,单位:N;Fxr为车辆的后轮纵向力,单位:N;vy为车辆质心处的侧向速度,单位:m/s;vx为车辆质心处的纵向速度,单位:m/s;r为车辆的横摆角速度,单位:rad/s;δ为车辆前轮转角,单位:rad;m为车辆的质量,单位:kg;Iz为车辆绕z轴的转动惯量,单位:kg·m2;a为车辆质心o到车辆前轴的距离,单位:m;b为车辆质心o到车辆后轴的距离,单位:m;
本发明的控制变量为车辆的前轮转角,主要考虑车辆的侧向动力学特性,所以假设车辆的纵向速度不变,则式(4)可以写成:
行驶稳定的车辆,其前轮转角δ很小,所以采用如下近似关系:
将式(6)代入到式(5)中得:
考虑到冰雪环境下,路面的附着系数低,当车辆行驶在冰雪路面时,轮胎力呈非线性的关系,所以选用如下的轮胎模型来表示轮胎侧向力与轮胎侧偏角的关系:
其中,Cf表示前轮的侧偏刚度,单位:N·rad;Cr表示后轮的侧偏刚度,单位:N·rad;Ka表示前轮的拟合参数;Kb表示后轮的拟合参数;αf表示前轮的侧偏角,单位:rad;αr表示后轮的侧偏角,单位:rad;Cf的取值在28000~90000的范围内,Cr的取值在28000~90000的范围内,Ka的取值大于0,Kb的取值大于0;
为了便于计算,将式(8)写成:
根据坐标系的规定,分别描述车辆前轮的轮胎侧偏角与后轮的轮胎侧偏角为:
将式(10)代入到式(9)中得:
将式(11)代入到式(7)中可得:
其中a1=(aAf-bAr)/(Izvx),a3=(a2Af+b2Ar)/(Izvx),b1=(Af+Ar)/(mvx),b3=(aAf-bAr)/(mvx)-vx,p0=(-aAf)/Iz,q0=(-Af)/m,Δ1≈Δ2=0;
将式(12)和式(13)中的所有vy转换成β得到车辆的动力学模型:
其中c1=(aAf-bAr)/(Iz),c2=(aBf-bBr)/(Iz),c3=(a2Af+b2Ar)/(Izvx),c6=3(a2Bf+b2Br)/(Izvx),d1=(Af+Ar)/(mvx), l0=(-aAf)/Iz,l1=(-3aBf)/Iz,l2=(-6a2Bf)/(Izvx),w0=(-Af)/(mvx),w1=(-3Bf)/(mvx),Δ3≈Δ4=0;
结合式(3)、式(14)和式(15),整理得到符合以下仿射形式的非线性车辆系统模型:
其中f1(x)=vx(x2+x4),f2(x)=x3, h(x)=yo,δ为前轮转角,
步骤二、采用反馈线性化的方法对车辆系统模型线性化
根据单输入单输出反馈线性化的基本原理,为了建立y与δ之间的关系,对式(16)的输出关于时间求导,利用李导数整理得:
Lfh(x)=vx(x2+x4),Lgh(x)=0 (18)
在x=0附近Lgh(x)=0,则再次关于时间求导,整理后得到:
在x=0附近LgLfh(x)≠0,所以车辆系统的相对阶为2。
对式(16)中的仿射形式的非线性车辆系统模型进行状态变换:
定义新的状态量z1=h(x),所以新的状态方程可以表示为:
非线性车辆系统转变为:
所以式(16)可以转化为以下形式:
其中C=[10];
由于(A,B)在x的取值范围内是可控的,γ(x)为非奇异矩阵,因此可对步骤一中所提出的仿射形式的非线性车辆系统模型进行反馈线性化;
设虚拟控制输入为v,并且求出δ与v的关系:
所以,非线性系统可以转换为以下线性系统:
其中C=[10];
步骤三、基于反馈线性化与LQR设计路径跟踪控制器
对于线性化后的车辆模型,利用线性二次型调节器的方法设计跟踪控制器,设理想的侧向位移为yd,则理想值与实际值之间的误差为yd-y=e,求最优控制输入v,使下列的二次型性能指标最小:
其中正定对称矩阵R为控制输入变量v的加权矩阵,半正定矩阵Q为状态变量z的加权矩阵,t0为控制时间的起点,tf为控制时间的终点;
当通过反馈线性化后的车辆系统最优控制输入为v=R-1BTPe=-Ke,二次型性能指标最小;
其中e可以表示为:
所以v可以进一步地表示为:
最优状态反馈控制增益矩阵为:
K=[K1 K2]=R-1BTP (30)
其中P是一个正定矩阵满足代数Riccati方程:
PA+ATP+Q-PBR-1BTP=0 (31)
此时δ与v的关系可以进一步地表示为:
将得到的控制量δ关系式输入到车辆系统中,实现无人驾驶车辆的路径跟踪。
Claims (1)
1.冰雪环境下基于反馈线性化与LQR的车辆路径跟踪方法,其特征在于,根据冰雪环境下车辆轮胎力不满足线性关系的问题,建立冰雪环境下仿射形式的非线性车辆系统模型,基于此模型,采用反馈线性化的方法对复杂的非线性车辆系统模型进行线性化处理,得到较简单的线性系统模型和虚拟控制输入;根据得到的简单线性车辆系统模型,使用线性二次型调节器的设计方法设计路径跟踪控制器,保证系统稳定性和目标的最优性,实现无人驾驶车辆路径跟踪;本方法的具体步骤如下:
步骤一、建立仿射形式的非线性车辆系统模型
假定车辆是一个刚性体,其中车辆装置着四个不会发生形变的车轮,并以车辆的前轮作为转向轮,根据车辆运动的几何关系以及车辆的运动学方程得到车辆运动学模型,如式(1)所示:
其中ψ为车辆横摆角,单位:rad;vo为车辆质心o处的速度,单位:m/s;β为车辆质心侧偏角,单位:rad;xO为车辆质心o的纵向位移,单位:m;yO为车辆质心o的横向位移,单位:m;
考虑到车辆在运行过程中主要是纵向行驶,因此vx>>vy,所以车辆行驶的速度v近似用纵向速度vx表示,同时,车辆的横摆角ψ和质心侧偏角β都较小,因此采用如下近似关系:
将式(2)代入到式(1)中进行简化处理,得到简化的车辆运动学模型:
对于车辆的动力学模型,以车辆质心o为坐标原点,以沿着车身向前的方向为x轴的正方向,垂直于横轴向上的方向为y轴的正方向,通过质心指向上方的方向为z轴的正方向,忽略车辆的纵向动力学,只考虑车辆的横摆方向动力学和侧向动力学,根据牛顿第二定律与力矩平衡关系推出如式(4)所示的单轨模型数学表达式:
其中,Fyf为车辆的前轮侧向力,单位:N;Fyr为车辆的后轮侧向力,单位:N;Fxf为车辆的前轮纵向力,单位:N;Fxr为车辆的后轮纵向力,单位:N;vy为车辆质心处的侧向速度,单位:m/s;vx为车辆质心处的纵向速度,单位:m/s;r为车辆的横摆角速度,单位:rad/s;δ为车辆前轮转角,单位:rad;m为车辆的质量,单位:kg;Iz为车辆绕z轴的转动惯量,单位:kg·m2;a为车辆质心o到车辆前轴的距离,单位:m;b为车辆质心o到车辆后轴的距离,单位:m;
本方法的控制变量为车辆的前轮转角,考虑车辆的侧向动力学特性,所以假设车辆的纵向速度不变,则式(4)可以写成:
行驶稳定的车辆,其前轮转角δ很小,所以采用如下近似关系:
将式(6)代入到式(5)中得:
考虑到冰雪环境下,路面的附着系数低,当车辆行驶在冰雪路面时,轮胎力呈非线性的关系,所以选用如下的轮胎模型来表示轮胎侧向力与轮胎侧偏角的关系:
其中,Cf表示前轮的侧偏刚度,单位:N·rad;Cr表示后轮的侧偏刚度,单位:N·rad;Ka表示前轮的拟合参数;Kb表示后轮的拟合参数;αf表示前轮的侧偏角,单位:rad;αr表示后轮的侧偏角,单位:rad;Cf的取值在28000~90000的范围内,Cr的取值在28000~90000的范围内,Ka的取值大于0,Kb的取值大于0;
为了便于计算,将式(8)写成:
根据坐标系的规定,分别描述车辆前轮的轮胎侧偏角与后轮的轮胎侧偏角为:
将式(10)代入到式(9)中得:
将式(11)代入到式(7)中可得:
其中a1=(aAf-bAr)/(Izvx),a3=(a2Af+b2Ar)/(Izvx),b1=(Af+Ar)/(mvx),b3=(aAf-bAr)/(mvx)-vx,p0=(-aAf)/Iz,q0=(-Af)/m,Δ1≈Δ2=0;
将式(12)和式(13)中的所有vy转换成β得到车辆的动力学模型:
其中c1=(aAf-bAr)/(Iz),c2=(aBf-bBr)/(Iz),c3=(a2Af+b2Ar)/(Izvx),c6=3(a2Bf+b2Br)/(Izvx),d1=(Af+Ar)/(mvx), l0=(-aAf)/Iz,l1=(-3aBf)/Iz,l2=(-6a2Bf)/(Izvx),w0=(-Af)/(mvx),w1=(-3Bf)/(mvx),Δ3≈Δ4=0;
结合式(3)、式(14)和式(15),整理得到符合以下仿射形式的非线性车辆系统模型:
其中f1(x)=vx(x2+x4),f2(x)=x3, h(x)=yo,δ为前轮转角,
步骤二、采用反馈线性化的方法对车辆系统模型线性化
根据单输入单输出反馈线性化的基本原理,为了建立y与δ之间的关系,对式(16)的输出关于时间求导,利用李导数整理得:
Lfh(x)=vx(x2+x4),Lgh(x)=0 (18)
在x=0附近Lgh(x)=0,则再次关于时间求导,整理后得到:
在x=0附近LgLfh(x)≠0,所以车辆系统的相对阶为2;
对式(16)中的仿射形式的非线性车辆系统模型进行状态变换:
定义新的状态量z1=h(x),所以新的状态方程可以表示为:
非线性车辆系统转变为:
所以式(16)可以转化为以下形式:
其中C=[1 0];
由于(A,B)在x的取值范围内是可控的,γ(x)为非奇异矩阵,因此可对步骤一中所提出的仿射形式的非线性车辆系统模型进行反馈线性化;
设虚拟控制输入为v,并且求出δ与v的关系:
所以,非线性系统可以转换为以下线性系统:
其中C=[1 0];
步骤三、基于反馈线性化与LQR设计路径跟踪控制器
对于线性化后的车辆模型,利用线性二次型调节器的方法设计跟踪控制器,设理想的侧向位移为yd,则理想值与实际值之间的误差为yd-y=e,求最优控制输入v,使下列的二次型性能指标最小:
其中正定对称矩阵R为控制输入变量v的加权矩阵,半正定矩阵Q为状态变量z的加权矩阵,t0为控制时间的起点,tf为控制时间的终点;
当通过反馈线性化后的车辆系统最优控制输入为v=R-1BTPe=-Ke,二次型性能指标最小;
其中e可以表示为:
所以v可以进一步地表示为:
最优状态反馈控制增益矩阵为:
K=[K1 K2]=R-1BTP (30)
其中P是一个正定矩阵满足代数Riccati方程:
PA+ATP+Q-PBR-1BTP=0 (31)
此时δ与v的关系可以进一步地表示为:
将得到的控制量δ关系式输入到车辆系统中,实现无人驾驶车辆的路径跟踪。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111576809.5A CN114148318B (zh) | 2021-12-22 | 2021-12-22 | 冰雪环境下基于反馈线性化与lqr的车辆路径跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111576809.5A CN114148318B (zh) | 2021-12-22 | 2021-12-22 | 冰雪环境下基于反馈线性化与lqr的车辆路径跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114148318A CN114148318A (zh) | 2022-03-08 |
CN114148318B true CN114148318B (zh) | 2023-10-27 |
Family
ID=80451681
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111576809.5A Active CN114148318B (zh) | 2021-12-22 | 2021-12-22 | 冰雪环境下基于反馈线性化与lqr的车辆路径跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114148318B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115447603B (zh) * | 2022-04-26 | 2024-05-24 | 重庆大学 | 一种应用于自动驾驶车辆的有限时间路径跟踪控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110539752A (zh) * | 2019-06-26 | 2019-12-06 | 江苏大学 | 一种智能汽车多预测范围模型预测轨迹跟踪控制方法及系统 |
CN112051851A (zh) * | 2020-09-09 | 2020-12-08 | 清华大学 | 极限工况下电动四驱车辆的自主漂移控制方法及系统 |
WO2021168669A1 (zh) * | 2020-02-25 | 2021-09-02 | 华为技术有限公司 | 车辆控制方法及装置 |
CN113386781A (zh) * | 2021-05-24 | 2021-09-14 | 江苏大学 | 一种基于数据驱动车辆动力学模型的智能车辆轨迹跟踪控制方法 |
CN113553726A (zh) * | 2021-08-06 | 2021-10-26 | 吉林大学 | 一种冰雪环境下的主从博弈型人机协同转向控制方法 |
-
2021
- 2021-12-22 CN CN202111576809.5A patent/CN114148318B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110539752A (zh) * | 2019-06-26 | 2019-12-06 | 江苏大学 | 一种智能汽车多预测范围模型预测轨迹跟踪控制方法及系统 |
WO2021168669A1 (zh) * | 2020-02-25 | 2021-09-02 | 华为技术有限公司 | 车辆控制方法及装置 |
CN112051851A (zh) * | 2020-09-09 | 2020-12-08 | 清华大学 | 极限工况下电动四驱车辆的自主漂移控制方法及系统 |
CN113386781A (zh) * | 2021-05-24 | 2021-09-14 | 江苏大学 | 一种基于数据驱动车辆动力学模型的智能车辆轨迹跟踪控制方法 |
CN113553726A (zh) * | 2021-08-06 | 2021-10-26 | 吉林大学 | 一种冰雪环境下的主从博弈型人机协同转向控制方法 |
Non-Patent Citations (2)
Title |
---|
基于模型预测控制的智能车辆路径跟踪控制器设计;王艺;蔡英凤;陈龙;王海;李健;储小军;;汽车技术(第10期);全文 * |
基于非线性模型预测控制的智能车路径跟踪算法;吴飞龙;郭世永;;汽车技术(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114148318A (zh) | 2022-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111890951B (zh) | 智能电动汽车轨迹跟踪与运动控制方法 | |
CN108674414B (zh) | 一种极限工况的智能汽车轨迹跟踪控制方法 | |
CN104977933A (zh) | 一种自主驾驶车辆的区域型路径跟踪控制方法 | |
CN107272692A (zh) | 基于微分平坦和自抗扰的无人车路径规划与跟踪控制方法 | |
CN113050651B (zh) | 智能车自主行驶路径跟踪的时滞控制方法及系统 | |
CN104881030A (zh) | 基于快速终端滑模原理的无人车侧纵向耦合跟踪控制方法 | |
CN106515716A (zh) | 轮式驱动电动汽车底盘集成控制系统协调控制装置及方法 | |
CN111796522B (zh) | 一种车辆状态估计方法 | |
CN113682282A (zh) | 一种车辆稳定性控制方法、系统、车辆和存储介质 | |
CN103600744A (zh) | 四轮转向/驱动车辆的路径保持和车轮侧滑的控制方法 | |
CN110239519A (zh) | 一种车辆爆胎过程的协调控制方法 | |
CN114148318B (zh) | 冰雪环境下基于反馈线性化与lqr的车辆路径跟踪方法 | |
Li et al. | Adaptive sliding mode control of lateral stability of four wheel hub electric vehicles | |
CN114851857B (zh) | 一种分布式驱动电动客车转矩控制方法 | |
CN109292018A (zh) | 基于同轴式轮腿结构的四轮转向轨迹跟踪控制方法 | |
CN109849898A (zh) | 基于遗传算法混合优化gpc的车辆横摆稳定性控制方法 | |
CN110968911B (zh) | 一种基于新型趋近律的汽车abs滑模控制器设计方法 | |
CN112346337A (zh) | 极限工况下基于后轮主动转向的车辆稳定性控制方法 | |
CN111994085A (zh) | 复杂路况下车辆行驶稳定区域的估计方法 | |
CN113602278B (zh) | 四轮独立驱动电动汽车分布式模型预测路径跟踪控制方法 | |
CN113044047B (zh) | 一种基于类pid-stsm的afs/dyc集成控制方法 | |
CN112835292B (zh) | 一种自动驾驶电动车在崎岖路面的横向控制方法 | |
Sun et al. | Intelligent vehicle path tracking based on feedback linearization and LOR under extreme conditions | |
CN110723200B (zh) | 一种转向回正及中间位置控制系统及其控制方法 | |
Wang et al. | MHE-based friction coefficient estimation and compensation method for path tracking control |
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 |