基于观测器的纯电动智能汽车纵向车速控制方法与流程

本发明属于汽车控制技术领域。

背景技术：

为了减少交通事故的发生，降低内燃机汽车对能源消耗及环境污染的影响，伴随着互联网、信息、电子和智能技术的发展，汽车的智能化及电动化技术已经成为解决上述问题的有效途径。近年来，大众、宝马、奥迪等知名汽车制造企业以及百度、谷歌等著名互联网企业，都在不断的增加对智能驾驶汽车领域的人力、财力的投入，以抢占智能驾驶的前沿技术。智能驾驶技术的发展必然会引领汽车行业新一轮的重大变革。纵向车速控制作为纯电动智能汽车的底层控制算法，其控制效果直接影响着智能汽车的安全性和乘车舒适性等性能。对于纯电动智能汽车，由于电机的响应速度快，电机的扭矩和转速容易获取，这为纯电动智能汽车的纵向车速控制提供了良好的基础条件。针对于集中式纯电动智能汽车纵向车速控制，主要有以下问题：

1.智能汽车纵向车速的控制，即通过设计控制器合理的产生驾驶需求力矩（包括驱动力矩和制动力矩），从而实现纵向车速的跟踪控制。

2.纯电动智能汽车纵向车速控制系统具有非线性。同时，控制器的输出要满足执行器电机和制动器的硬性约束，即驱动和制动力矩信号不能超过电机的实际最大输出力矩和制动器的最大制动力矩。

3.电动汽车需要动力电源，现在常用的是锂电池组，来给电机供电，电机的供电电压也影响着电机的最大输出力矩，因此在考虑电机的最大输出力矩时也必须考虑电池组输出电压的影响，即电机实际最大输出力矩是一个变约束。

4.系统参数的不可测，整车质量作为影响系统模型的重要参数，其随着乘车人的多少及体重而改变，而其并不可测。

技术实现要素：

本发明的目的是利用基于观测器的滚动时域优化控制算法来设计控制器，通过控制器优化出驾驶员需求的力矩，然后进行驱动和制动力矩分配，从而实现纵向车速有效控制的基于观测器的纯电动智能汽车纵向车速控制方法。

本发明实现matlab/simulink和amesim的联合仿真，

①必须设置pc电脑的环境变量，让两者相互关联；

②在amesim界面中添加与simulink通信的接口模块，将matlab/simulink和amesim间需要通信的变量连接到这个模块；

③经过系统编译之后，amesim中的模型信息以s-function的形式保留在simulink中，从而实现两者的联合仿真与通信。

本发明的步骤是：

一、集中式电动汽车仿真模型搭建：

电动汽车仿真模型包括电驱动模块、传动模块、轮胎模块以及车辆纵向动力学，整车模型参数如表一

表一电动汽车参数表

二、基于观测器的滚动时域优化控制器：

2.1面向控制器设计模型搭建

2.1.1车辆纵向动力学模型

在不考虑横向力的情况下，车辆坡路上行驶纵向受力根据牛顿第二定律有：

（1）

其中：为行车质量，为驱动力、为行驶阻力；包括空气阻力、路面摩擦阻力、坡度阻力及机械制动力；

车重与行车质量关系用式下式表示：

（2）

其中为一个车轮的惯量，为车轮半径；

在坡路上行驶的汽车受到坡度阻力为：

（3）

其中为重力加速度；

路面上行驶的汽车受到的空气阻力为：

（4）

其中为空气黏性密度，为风阻系数，为车辆的迎风面积，为风速，为车速；

忽略汽车风速影响，故空气阻力表示为：

（5）

摩擦阻力是道路与轮胎间的摩擦力，通过下式：

（6）

其中为路面摩擦系数，为粘滞摩擦系数；

机械制动力，为制动力矩；

得到车辆受到的行驶阻力为：

（7）；

2.1.2传动系建模

2.1.2.1离合器

由刚性假设，其传递的转矩为：

（8）

其中为电机输出转矩，为离合器输出转矩，为电机输出转速，为离合器输出转速；

2.1.2.2变速器

变速器输出转矩建模如下式：

（9）

其中为扭转阻尼系数，为输出转速，为传动比，为档位传动比，为主减速比；

2.1.2.3驱动轴

（10）

其中为驱动轴输出，为驱动轴输出转速；

将式8和式9代入式10整理得：

（11）

用驱动力，用表示车轮半径，则由力和力矩之间的关系，同时车速，故结合式11得：

（12）

其中车轮半径由下式求得，式中为轮毂半径，为轮胎扁平比，为轮胎宽度；

（13）；

2.2联合观测器：

2.2.1递归最小二乘法质量辨识

结合式1和式7，得到如下等式：

（14）

同时结合式（2）、（3）、（4）、（5）、（6）整理成最小二乘格式，驱动轴力矩估计值，转化为驱动力，得

（15）

其中表示车辆纵向加速度，等效旋转质量，其值，其中为一个车轮的惯量，为车轮半径；表示包含驱动轴估计量系统输出量，表示可获取的数据向量，为待辨识量，为系统的过程白噪声；

根据最小二乘法原理，分别定义k-1、k时刻系统辨识得到的整车质量为、，则得到质量辨识模型：

（16）

式中，为第k时刻的遗忘因子；

遗忘因子规律为：

（17）；

2.2.2驱动轴力矩观测器

设驱动轴两端转速可测，则搭建驱动轴力矩生成模型：

（18）

其中为开环计算驱动轴力矩，、为驱动轴两端转速，为变速箱输出转速，为车轮转速，为驱动轴等效刚度系数，为驱动轴等效阻尼系数；

得到等效的车轮动力学模型：

（19）

其中为车轮转速估计值，为驱动力矩，为阻力矩，为驱动半轴末端的转动惯量，

（20）

驱动力矩通过阻力数学模型求得，即，结合式7得

（21）

其中表示坡度阻力矩，表示空气阻力矩，表示摩擦阻力矩，为机械制动力矩，表示行车阻力矩；驱动轴力矩估计器时用代替；

定义偏差，即车轮转速估计值减去实际值，对偏差求导，同时结合式19和式21得：

（22）

取控制输入为如下形式：

（23）

其中，为虚拟控制输入，利用反馈线性化方法，式22线性化形式

（24）

将虚拟控制输入设计成pi形式，即：

（25）

其中，为比例系数，为积分系数；

联合式23和25，得到控制器：

（26）

定义lyapunov函数如式（27）：

（27）

对两侧求导得：

（28）

将式（25）带入上式并整理得：

（29）

因此有时，故而；

因此驱动轴力矩为：

（30）。

本发明基于观测器的纵向车速控制器：

电机力矩和机械制动力矩之间在数值上通过传动比进行转化，即：

（31）

即制动力矩和驱动力矩统一表示为，

通过式(1)、(7)和(12)整理得到如下关系式：

（32）

其中选取纵向车速为状态量，即，选取驾驶需求为控制变量，即，同样选取纵向车速作为输出量，即；通过欧拉方法对状态空间方程进行离散化，用表示采样步长，则在时刻，离散化得到系统离散模型如下：

（33）

为输出系数矩阵，定义系统的预测时域为，控制时域为，需要满足，那么在时刻的预测输出序列表示为：

（34）

同时，k时刻的优化控制输入序列表示为：

（35）

在采样时刻，状态量的值为，推导出状态量和输出量的预测过程如公式(19)和(20)所示：

（36）

（37）

同时，参考输入为期望的车速，故得到参考输入序列：

（38）

在控制器设计时需要考虑如下约束：

（39）

同时为了保障纵向车速的跟踪，同时提高乘坐的舒适性，因此控制器选取性能指标目标为：

（40）

其中；

进而描述为(41)式的优化问题，即使目标函数值最小：

（41）

式(41)中，反映了实际输出车速与期望车速的偏差，反映了驾驶需求的强弱，和分别为输出信号序列和控制信号序列的权重因子；使用nag求解41式的优化问题，优化出系统的控制输入序列，然后将序列中第一个元素作用于系统；下一时刻，重复上述的优化求解过程，即对自主驾驶纵向车速的闭环优化控制，

（42）；

需要对优化出的力矩进行力矩分配将驾驶员力矩需求中大于0的部分化为驱动力矩，小于等于0的部分化为制动力矩，即

（43）

其中为优化的驾驶力矩需求，为电机期望力矩，为制动力矩，将制动力矩通过下式转化为机械制动信号：

（44）。

本发明主要针对纯电动智能汽车纵向车速控制问题，针对系统的重要参数设计观测器，滚动时域优化控制算法能够很好的完成在线优化求解，同时能够显性的处理约束。具体的说通过建立纵向车速控制系统的机理模型，得到驾驶力矩需求的预测方程，然后构造代价函数，同时充分考虑约束条件，优化求解得到最优的驾驶员力矩需求。

本发明的有益效果是：

1.在纵向车速控制方面，传统的控制算法多数都是不基于模型的，而在实际道路上行驶的车辆工况复杂多变，很难找到一组控制器参数满足所有工况。同时将估计器加入到控制器设计过程中，抑制了车辆质量这一参数对汽车纵向车速控制的影响，基于观测器的滚动时域优化控制算法在控制器设计时基于系统的机理模型，将观测的车辆质量等一些车辆运行工况信息直接引入机理模型，从而使车辆纵向车速控制更为精确。

2.本发明中设计的纵向车速控制系统是一个非线性系统，并且考虑到了电机、电池组和制动器的执行器硬性约束，传统的控制算法并不能有效的处理系统的约束，而滚动时域优化控制算法能有效地处理带约束的控制问题，在求解时直接将约束条件编译到simulink中的s_function中在线求解。

附图说明

图1是实施本发明所述的基于滚动时域优化纵向车速控制框图；

图2是实施本发明所述的集中式电动汽车amesim整车模型；

图3是本发明所述的车辆坡路上行驶纵向受力示意图；

图4是本发明所述驱动轴力矩估计方案；

图5是本发明所述电机最大输出转矩map图；

图6是本发明中设计扭基于滚动优化算法的纵向车速控制器设计流程图；

图7是本发明进行控制器效果验证时采用的复杂城市工况的期望车速，横坐标为时间，单位s，纵坐标是纵向车速，单位为m/s;

图8是本发明所述的在平路工况下的仿真结果，从上到下依次为质量辨识对比曲线、驱动轴力矩估计对比曲线和车速跟踪曲线。其中在质量辨识对比曲线中实线是实际车辆质量，虚线是辨识质量。力矩对比曲线中虚线为估计力矩，实线为实际力矩。实际车速和期望车速对比图中，虚线为实际车速，实线为期望车速；

图9是本发明所述的在恒定大坡度工况下的仿真结果，从上到下依次为质量辨识对比曲线、驱动轴力矩估计对比曲线和车速跟踪曲线。其中在质量辨识对比曲线中实线是实际车辆质量，虚线是辨识质量。力矩对比曲线中虚线为估计力矩，实线为实际力矩。实际车速和期望车速对比图中，虚线为实际车速，实线为期望车速；

图10是本发明所述的变坡度仿真时的道路坡度变化情况；

图11是本发明所述的在接近真实路面变坡度工况下的仿真结果，从上到下依次为质量辨识对比曲线、驱动轴力矩估计对比曲线和车速跟踪曲线。其中在质量辨识对比曲线中实线是实际车辆质量，虚线是辨识质量。力矩对比曲线中虚线为估计力矩，实线为实际力矩。实际车速和期望车速对比图中，虚线为实际车速，实线为期望车速。

具体实施方式

本发明中基于数据驱动预测控制的电动汽车扭矩优化方法实施的控制框图如图1所示，图中车速优化控制器是在simulink中搭建的，控制器的输入为期望车速，实际车速作为可测信号，质量和坡度作为可测干扰实时反馈给控制器，tmax是电机最大驱动力矩，它是由电机的机械特性和电池的输出电压共同决定的，既体现了电机本身的执行器硬性约束条件，又体现出了随着电池放电时间增加电压减小对整车性能的影响。整车控制器得到的驱动力矩必须小于等于tmax，因此tmax是作为约束给到控制器的。图2中集中式纯电动汽车模型是在amesim中搭建的，用来模拟实车的运行。控制器优化驾驶需求力矩通过力矩分配为驱动和制动力矩信号分别给到电机和制动模块，控制车辆的运行，而车辆的实际车速作为反馈信号反馈给控制器。

本发明的控制目标是，纵向车速控制器根据实时反馈回来的实际车速与期望车速信号对比，在满足约束条件前提下，优化得到驾驶需求力矩，然后通过力矩分配，从而得到驱动力矩和制动力矩信号，并给到整车模型中的电机和制动模块，控制车辆的运行，最终让实际车速跟踪上期望车速。

本发明提供了一套基于以上运行原理和运行过程的装置。即基于pc机的离线电动汽车扭矩优化设计试验平台。搭建以及运行过程如下：

软件选择

该控制系统的被控对象和控制器的仿真模型分别通过软件matlab/simulink和amesim进行搭建，软件版本分别为matlabr2009a和amesimr10，求解器分别选择为ode3和euler。仿真步长为定步长，步长选择为0.01s。

本发明要实现matlab/simulink和amesim的联合仿真。

①首先必须按照要求设置pc电脑的环境变量，让两者相互关联。

②然后在amesim界面中添加与simulink通信的接口模块，将matlab/simulink和amesim间需要通信的变量连接到这个模块；

③最后经过系统编译之后，amesim中的模型信息以s-function的形式保留在simulink中，从而实现两者的联合仿真与通信。

在运行simulink仿真模型时，amesim模型也在同时进行计算和求解。仿真过程中两者之间不断进行数据的交换。如果对amesim中的模型结构或者参数设置进行了修改，则需要重新编译。值得注意的是，两者的仿真步长必须一致。

本发明步骤是：

一、集中式电动汽车仿真模型搭建：

如图2所示，整个电动汽车仿真模型包括电驱动模块、传动模块、轮胎模块以及车辆纵向动力学等几个部分，整车模型参数如表一所示。

表一电动汽车参数表

。

电驱动系统包括电池部分和电机部分，纯电动汽车的电池组是锂电池组，由多个单体电池串并联构成。电池组对外输出的终端电压为单个电池输出电压和，电池系统输出终端电压即电池组提供给电机的电压；本发明中采用的是永磁同步电机。

传动系统包括变速器，差速器和驱动轴三个部分。电机输出的动力经由变速器通过不同的齿轮半径产生不同速比从而进行减速增矩，本发明忽略了车辆的侧向动态，那么差速器两侧的输出转速相同，即差速器不起作用，差速器输出转速即为驱动轴的输入转速，而驱动轴的输出转速与车轮转速相等，轴上传递的力矩通过驱动轴两端转速差计算。电机输出的力矩经由变速器通过不同的齿轮半径产生不同速比从而进行减速增矩，本模型主减速比为2.2786，档位减速比为3.9431，即传动比为8.6847。

车辆纵向动力学部分，其中考虑了车辆在行驶过程中受到驱动力，制动力和行驶阻力的作用，其中行驶阻力包括空气阻力，滚动阻力和摩擦阻力。在这个模块中可以设置车辆的整体质量、坡度、风速等参数。

二、基于观测器的滚动时域优化控制器：

2.1面向控制器设计模型搭建

2.1.1车辆纵向动力学模型

为了实现纵向车速控制和车辆质量参数估计的研究，需要建立车辆纵向动力学模型。在不考虑横向力的情况下，车辆坡路上行驶纵向受力示意图如图3所示。图3中为道路坡度，为行车质量，则重力为，为驱动力、为空气阻力以，为路面摩擦阻力。根据牛顿第二定律有：

（1）

其中：为行车质量，为驱动力、为行驶阻力；包括空气阻力、路面摩擦阻力、坡度阻力及机械制动力。

需要指出的是车重与行车质量并不相等，计算行车质量时受到行驶方向上惯性效果的影响，其关系可以近似用式下式表示：

（2）

其中为一个车轮的惯量，为车轮半径；驱动力/制动力在下一小节做详细介绍，下面对车辆行驶过程中受到的行驶阻力进行分析介绍。

在坡路上行驶的汽车受到坡度阻力为：

（3）

其中为重力加速度；为车重。考虑到amesim模型中坡度为按百分比计算的道路坡度，为了后续仿真实验时坡度形式的统一，这里我们统一按百分制计算坡度，那么需要转化。

路面上行驶的汽车受到的空气阻力为：

（4）

其中为空气黏性密度，为风阻系数，为车辆的迎风面积，为风速，为车速；这里由于空气阻力相对于坡度阻力和摩擦阻力数值较小，然后在城市中行驶的汽车风速相对于车速又较小，在控制器设计时忽略了风速的影响，故空气阻力表示为：

（5）

摩擦阻力是道路与轮胎间的摩擦力，通过下式：

（6）

其中为路面摩擦系数，为粘滞摩擦系数；

机械制动力，为制动力矩；

得到车辆受到的行驶阻力为：

（7）；

2.1.2传动系建模

建模时对离合器、传动轴、驱动轴进行刚性假设，同时忽略了主减速器和档位变速器之间传递的力矩损失。

2.1.2.1离合器

由刚性假设，其传递的转矩为：

（8）

其中为电机输出转矩，为离合器输出转矩，为电机输出转速，为离合器输出转速。

2.1.2.2变速器

这里由于我们忽略了主减速器和档位变速器之间传递的力矩损失，故在建模时我们对其进行了统一，变速器输出转矩建模如下式：

（9）

其中为扭转阻尼系数，为输出转速，扭转阻尼系数和输出转速的乘积用于近似摩擦力矩损失，为传动比，为档位传动比，为主减速比。

2.1.2.3驱动轴

（10）

其中为驱动轴输出，为驱动轴输出转速；即为车轮转速。

将式8和式9代入式10整理得：

（11）

用驱动力，用表示车轮半径，则由力和力矩之间的关系，同时车速，故结合式11得：

（12）

其中车轮半径由下式求得，式中为轮毂半径，为轮胎扁平比，为轮胎宽度；

（13）；

2.2联合观测器：

为了准确估计车辆质量，我们设计质量和驱动轴力矩联合观测器。我们考虑到质量和驱动轴力矩之间存在耦合关系，即估计质量是需要用到驱动轴力矩信息，估计力矩是需要使用车辆质量信息。因而在确定估计方案时，分析两个待观测量各自的特性，汽车的质量的改变主要是由乘车人数、油箱中油的多少以及装载货物的多少引起，质量在车辆行驶过程中相对稳定，因而车辆质量是慢变量，可以用辨识的方法进行辨识，同时考虑这一特性用当前时刻的质量辨识结果对下一时刻的力矩进行估计对力矩估计结果产生影响很小，这样可以提高算法的效率，同时也解决了二者之间的耦合关系。

2.2.1递归最小二乘法质量辨识

分析车辆在坡路上行驶的动力学模型，结合式1和式7，得到如下等式：

（14）

同时结合式（2）、（3）、（4）、（5）、（6）整理成最小二乘格式，同时结合分析，估计质量时用到驱动轴力矩估计值，转化为驱动力，得

（15）

其中表示车辆纵向加速度，等效旋转质量，其值，其中为一个车轮的惯量，为车轮半径；表示包含驱动轴估计量系统输出量，表示可获取的数据向量，为待辨识量，为系统的过程白噪声；需要指出的是在辨识车辆质量过程中我们将驱动轴力矩当成可测量，那么根据力和力矩之间的关系，驱动力也是可测量，同时认为纵向加速度和机械制动力矩为实时可测参数。

根据前面介绍的最小二乘法原理，分别定义k-1、k时刻系统辨识得到的整车质量为、，则得到质量辨识模型：

（16）

式中，为第k时刻的遗忘因子。

车辆质量是慢变量，由于设定的质量初始值与实际质量可能存在较大的偏差，因此在辨识刚开始时，需要设定一个较大的信任度衰减，即此时选取的遗忘因子数值较小，随着辨识的不断进行，辨识结果会收敛到实际车辆质量附近，那么就需要一个较大的遗忘因子，以获取较小的信任度衰减，故本文中选取的遗忘因子规律为：

（17）；

2.2.2驱动轴力矩观测器

本技术设计了闭环驱动轴力矩观测器，通过力矩生成模型开环计算结果，然后通过车速偏差对开环计算结果进行矫正，从而将驱动轴力矩观测问题转化为转速跟踪问题。估计方案如图4所示。考虑到负载力矩的非线性的特点，结合反馈线性化方法的优势，故采用反馈线性化方法进行力矩观测器设计。

设驱动轴两端转速可测，则搭建驱动轴力矩生成模型：

（18）

其中为开环计算驱动轴力矩，、为驱动轴两端转速，为变速箱输出转速，为车轮转速，为驱动轴等效刚度系数，为驱动轴等效阻尼系数。

根据拉格朗日动力学方程，并且假设汽车行驶中车轮做纯滚动无滑动运动，故得到等效的车轮动力学模型：

（19）

其中为车轮转速估计值，为驱动力矩，为阻力矩，为驱动半轴末端的转动惯量，近似用式20计算得到

（20）

同时考虑力和力矩之间的关系，驱动力矩通过阻力数学模型求得，即，结合式7得

（21）

其中表示坡度阻力矩，表示空气阻力矩，表示摩擦阻力矩，为机械制动力矩，表示行车阻力矩；同时考虑到在进行驱动轴力矩估计时采用的是上一时刻质量的表示结果，即在计算时采用的是，故在设计驱动轴力矩估计器时用代替。

定义偏差，即车轮转速估计值减去实际值，对偏差求导，同时结合式19和式21得：

（22）

取控制输入为如下形式：

（23）

其中，为虚拟控制输入，利用反馈线性化方法，式22线性化形式

（24）

将虚拟控制输入设计成pi形式，即：

（25）

其中，为比例系数，为积分系数；

联合式23和25，得到控制器：

（26）

定义lyapunov函数如式（27）：

（27）

对两侧求导得：

（28）

将式（25）带入上式并整理得：

（29）

因此有时，故而；

综上所述，估计得到的驱动轴力矩为：

（30）。

本发明基于观测器的纵向车速控制器：

本文的控制目标是实现自主驾驶过程中纵向车速的跟踪控制，在不同工况下，通过控制器优化出期望的驾驶需求，然后通过电机和机械制动等执行机构实现实际车速跟踪上期望车速。

为了方便控制器设计优化，我们将机械制动力矩和电机制动力矩统一进行优化，即在考虑电机为理想电机的情况下，优化出达到期望车速时的驾驶员需求，然后在通过一定的控制策略将其分为电机力矩需求和机械力矩需求。

为了实现统一优化驾驶员力矩需求，需要建立机械制动力矩和电机力矩之间关系，这里我们忽略传动系的惯性损失，那么电机力矩和机械制动力矩之间在数值上通过传动比进行转化，即：

（31）

即制动力矩和驱动力矩统一表示为，综上分析集合前面建立的系统模型，同时为了更准确的实现车速控制，我们在设计控制器时整车质量采用联合估计器估计的质量。

通过式(1)、(7)和(12)整理得到如下关系式：

（32）

其中选取纵向车速为状态量，即，选取驾驶需求为控制变量，即，同样选取纵向车速作为输出量，即；通过欧拉方法对状态空间方程进行离散化，用表示采样步长，则在时刻，离散化得到系统离散模型如下：

（33）

为输出系数矩阵，根据模型预测控制理论，定义系统的预测时域为，控制时域为，需要满足，那么在时刻的预测输出序列表示为：

（34）

同时，k时刻的优化控制输入序列表示为：

（35）。

在采样时刻，状态量的值为，根据模型预测控制的基本原则和相关理论，推导出状态量和输出量的预测过程如公式(19)和(20)所示：

（36）

（37）。

通过分析计算当前时刻的状态变量值和上一时刻系统输入值，优化求解出控制量序列，而只将优化求解出的控制量序列的第一个量作用于系统。而在下一采样时刻，电动汽车模型会反馈回新的输入变量和状态量，控制器重新优化求解控制问题。

同时，参考输入为期望的车速，故得到参考输入序列：

（38）

在自主驾驶过程中，为了保证驾驶安全，需要对状态量进行约束，同时考虑到电机本身特性，其输出力矩具有约束，同时考虑到机械结构的限制，又有最大制动力矩的约束，综上所述，在控制器设计时需要考虑如下约束：

（39）

其中电机提供的最大力矩由图5中的map查表得到。

同时为了保障纵向车速的跟踪，同时提高乘坐的舒适性（即加速和制动过程中控制动作尽可能小），，因此控制器选取性能指标目标为：

（40）

其中；

进而描述为(41)式的优化问题，即使目标函数值最小：

（41）

式(41)中，反映了实际输出车速与期望车速的偏差，反映了驾驶需求的强弱，和分别为输出信号序列和控制信号序列的权重因子。的大小反映速度跟踪精度的要求，越大，速度跟踪的偏差越接近零。则反映对控制动作的要求，越大，控制动作越小，乘车舒适性更好。使用nag(一种滚动时域优化算法matlab求解工具箱)求解41式的优化问题，优化出系统的控制输入序列，然后将序列中第一个元素作用于系统；下一时刻，重复上述的优化求解过程，即对自主驾驶纵向车速的闭环优化控制，

（42）；

滚动时域优化控制器设计的流程如图6所示：滚动时域优化控制器优化出驾驶员力矩需求，但在控制中需要的控制信号是电机力矩需求和机械制动信号，故需要对优化出的力矩进行力矩分配将驾驶员力矩需求中大于0的部分化为驱动力矩，小于等于0的部分化为制动力矩，即

（43）

其中为优化的驾驶力矩需求，为电机期望力矩，为制动力矩，将制动力矩通过下式转化为机械制动信号：

（44）。

实验验证

反复调节控制参数，分别选取输出信号序列和控制信号序列的权重因子γy=100，γu=2，采样时间为0.01s，我们选择了加减速频繁的城市工况，期望车速如图7。设定整车质量为1500kg，分别在平路工况、恒定大坡度和贴近真实路面的变坡度工况下对控制器进行验证。

1）平路工况仿真验证

首先我们选择水平路面验证，设定道路坡度为0，在仿真时设定风速为0，车辆质量为1500kg，即车辆空载运行，仿真结果如图8所示，图中给出从上到下依次为质量辨识对比曲线、驱动轴力矩估计对比曲线和车速跟踪曲线，从图中可以看出估计器和控制器都具有良好的效果。

2）恒定大坡度工况仿真验证

在仿真环境下，我们设定道路坡度为10%，即恒定较大坡路上，验证在长时间持续较大坡路行驶控制器的控制效果是否稳定，仿真结果如图9所示。图中给出从上到下依次为质量辨识对比曲线、驱动轴力矩估计对比曲线和车速跟踪曲线。仿真结果表明，在较大坡路上行驶，估计器估计效果很好，实际车速在大部分时间可以跟踪上期望车速，但在200~300s左右，实际车速并没有跟踪上期望车速，而是车速维持在20m/s，这是由于电机期望力矩受到了控制器设计过程中我们电机最大电机力矩约束的作用，通过图5电机最大输出力矩map可知，此时的电机期望力矩刚好等于电机最大输出力矩。说明控制器效果良好，同时系统约束起到了很好的作用。

3）变坡度工况仿真验证

在实际车辆运行环境中，道路坡度并不会保持不变，因而我们设定更接近实际运行工况的变坡度（道路坡度如图10）进行验证。仿真结果如图10一所示，从仿真结果中可以看出，联合估计器估计效果良好，同时在变坡度工况下车辆纵向车速很好的跟踪期望车速。

本发明针对纯电动智能汽车设计了基于滚动时域优化方法的纵向车速控制器，这种方法很好的实现了在线优化，同时显性的处理约束。为了验证纵向车速优化控制器的有效性，在amesim高级仿真软件中搭建了集中式电动汽车模型，并且在复杂城市道路，分别在平路工况、恒定大坡度工况和贴近实际路面的变坡度工况下验证了控制器的性能。仿真结果表明，滚动时域优化纵向车速控制器在不同行驶工况下，具有良好的控制性能。