Matlab汽车运动控制系统设计

1绪论

1.1选题背景与意义

汽车已经成为人们日常生活不可缺少的代步交通工具，在汽车发达国家，旅客运输的60％以上，货物运输的50％以上由汽车来完成，汽车工业水平和家庭平均拥有汽车数量已经成为衡量一个国家工业发达程度的标志。进行汽车运动性能研究时．一般从操纵性、稳定性和乘坐舒适性等待性着手。但近年来．随着交通系统的日趋复杂，考虑了道路环境在内的汽车运动性能开始受到关注。因此，汽车运动控制系统的研究也显得尤为重要，在文中，首先对汽车的运动原理进行分析，建立控制系统简化模型，确定期望的静态指针(稳态误差)和动态指针(超调量和上升时间)。然后对汽车运动控制系统进行设计分析。从而确定系统的最佳静态和动态指针。

2 论文基本原理分析

2.1.1汽车运动横向控制

（1）绝对位置的获得方法

汽车横向方向的控制使用GPS(全球定位系统)的绝对位置信息。GPS信息的精度与采样周期、时间滞后等有关。为提高GPS的数据精度和平滑数据．采用卡尔曼滤波对采样数据进行修正。GPS的采样周期为200ms相对应控制的周期采用50ms。另外考虑通信等的滞后、也需要进行补偿，采用航位推测法(dead reckoning)解决此问题。通过卡尔曼滤波和航位推测法推算出的值作为汽车的绝对位置使用来控制车速、横摆角速度等车辆的状态量。GPS 的数据通过卡尔曼滤波减少偏差、通过航位推测法进行误差和迟滞补偿．提高了位置数据推算的精度。

（2）前轮转角变化量的算出方法

这里对前轮目标转角变化量(∆δ)的算出方法作简要说明，横方向控制采用预见控制，可以从现在汽车的状态预测经过时间t p秒后的汽车位置，由t p秒后的预测位置和目标路径

的位置可以算出t p秒后为沿着目标路径行驶所需要的汽车横摆角速度(ωr)。这个数值前回馈或者从与现在值的目标路径的误差的反馈来推算前轮目标转角变化量(式(1))．

∆δ[k]=(k1∙ωr[k]+k2∙εr[k])∙T c式(1)

式中T c为控制周期，k1, k2根据与现在目标路径的误差(ε)最小的原则来求解。

2.1.2汽车运动纵方向的控制

建立一个合理的传动系统模型是设计高性能汽车纵向运动控制系统的基础。目前纵向运动控制器设计过程中采用的传动系统简化模型主要有两类：一类是忽略传动系统的部分动态特性得到简化模型：另一类是通过对输入输出特性辨识得到简化模型。本文借鉴文献，忽略传动系统的部分动态特性，将车辆简化为两轮模型，对于前轮驱动车辆，整车受力如图1所示。前后车轮运动方程分别为

{J fωf=T s−rW f f−rF f J rωr=rF r−rW r f

上式中J f和J r，分别为前后轮转动惯量(左右轮之和)，ωr为后轮转速，W r和W f分别为前后车轮的垂直载荷(左右轮之和)，F r和F f分别为前后轮切向力(左右轮之和)，r为车轮半径，f为滚动阻力系数。对于汽车纵向运动控制系统，不会出现非常大的加减速度，采用线性化轮胎模型，得到切向力与滑移率关系为：F=K w s。式中K w为轮胎纵向刚度，s 为滑移率。驱动时s=l-v／(r)，制动时s=l-(r)／v。认为风阻F w作用于汽车质心，则前后轮垂直载荷分别为

{W f=M g

b

L

−Mv

h g

L

W r=M g a

+Mv

h g

式中M为整车品质；a和b分别为前后轴到质心的距离，L=a+b；h g为质心至地面的高

度。整车运动方程为Mv=F f−F r−C

αAv2 ,式中C

α

风阻系数，A为等效迎风面积。

发动机转矩、发动机转速、涡轮转速、半轴转矩、前后轮转速、车速7状态的非线性传动系统模型，在低频带内，发动机动态对传动系统特性基本无影响，如果控制系统只涉及较低频段．可以忽略发动机动态。忽略了半轴、轮胎滑移以及载荷转移和发动机转矩，只包括发动机转速、车速2个状态。

飞轮运动方程为J eωe=T edes−T p

整车运动方程为Mv=T i R g R d

r

−CαAv2−Mgf

2.2汽车运动控制系统的模型简化分析

考虑图2所示的汽车运行控制系统。如果忽略车轮的转动惯量．并且假定汽车受到的摩擦阻力大小与运动速度成正比，方向与汽车运动方向相反，则该系统可以简化成简单的质量阻尼系统。根据牛顿运动定律，该系统的模型(亦即系统的运动力方程)表示为：

{m dv

dt

+bv=u

y=v

（3-1）

其中，u为汽车的驱动力。

为了得到控制系统的传递函数，对式(3-1)进行Laplace变换。假定系统的初始条件为零，则动态系统的Laplace变换为：

{msV (s )+bV (s )=U(s)Y (s )=V(s)

由于系统输出是汽车的运动速度，用Y(S)替代V(s)，得到：

msY (s )+bY (s )=U(s)

因此．汽车运动控制系统模型的传递函数为：

Y(s)U(s)=1

ms +b

2.3汽车控制系统PID 控制器的校正

根据阶跃响应曲线．利用串联校正的原理．以及参数变化对系统响应的影响，对静态和动态性能指针进行具体的分析，最终设计出满足我们需要的控制系统。系统在未加入任何校正环节时的开环传递函数，在MATLAB 环境下对系统未加校正时开环阶跃响应曲线进行仿真．绘制如图3阶跃响应曲线，图中系统的开环响应曲线未产生振荡，其上升时间约100秒，稳态误差达到98％，远不能满足跟随设定值的要求。

图 3

图 4

（1）首先选择P校正，也就是在系统中加入一个比例放大器，为了大幅度降低系统的稳态误差，同时减小上升时间。

P校正后系统的闭环传递函数为：

Y(s) U(s)=

K p

ms+(b+K p)

此时控制系统的稳态值为Kp/(b+Kp)=Kp／(50+Kp)。本系统的比例增益Kp=800。即稳态值为800／(50+800)=O.941，这样可以把系统的稳态误差降低到0．06左右。加入P 校正后控制系统的死循环阶跃响应曲线如图3所示。图中，系统的稳态值约为0.941．稳态误差约为5.9％，这和最初的设计要求仍有差距，并且上升时间在7秒左右，不能达到设计的需要。因此我们选择PI校正。

（2）加入PI校正器后系统的闭环单位反馈传递函数为：

Y(s) U(s)=

K p s+K i

ms2+(b+K p)s+K i

考虑到Ki的作用．我们可以大幅度降低kp，取kp=200。Ki=70，在MATLAB环境下仿真得出的系统响应曲线如图4(中)所示。从图4(中)中可以得知，加入PI校正后系统的上升时间有所下降，但仍大于5秒。同时又产生了另一个问题，系统的超调量达到了26．43％．这是使用积分器带来的副作用。因此适当地加入微分量。

（3）可以选择PD校正，此时系统的闭环单位反馈传递函数为：

Y(s) U(s)=

K D s+K p

(m+K D)s+(b+K p)

鉴于K D对上升时间和稳态误差影响不大．我们在P校正的基础上．将K D降低少许，给出K D=10。系统响应曲线如图4(中)所示。

（4）加入PID校正，此时系统的闭环单位反馈传递函数为：

Y(s) U(s)=

K D s2+K p s+K i

(m+K D)s2+(b+K p)s+K i

K p，K i和K D的选择一般先根据经验确定一个大致的范围，然后通过MATLAB绘制的图形逐步校正。这里我们取K p=700，K i=100，K D=100。得到加入PLD校正后系统的死循环阶跃