第一作者:刘伯鸿(1968—),男,甘肃临洮人,副教授. 研究方向为交通信息工程及控制.email:liubohong@mail.lzjtu.cn.
针对高速列车在复杂环境运行时,传统预测控制器出现的动力学模型失配、在线辨识实时性不强等问题,提出一种基于RBF-ARX模型的列车预测控制方法.首先,构建RBF-ARX模型作为列车控制器的预测模型,使用SNPOM算法处理列车运行中的大量历史数据,得到预测模型的模型阶次及优化参数,并验证了模型的准确性.然后使用该模型构造预测控制器对列车目标曲线进行在线跟踪仿真.结果表明:基于RBF-ARX模型的高速列车预测控制器不仅解决了动力学模型失配的问题,并且在控制精度及舒适度等方面都优于传统预测控制器.
When high-speed train are operating in the complex area, its tradition predictive controller will mismatch its dynamic model and not identify its model on-line.In order to solve this, a train predictive control method based on RBF-ARX model is proposed. Firstly, RBF-ARX model is constructed as the predictive model of train controllers, and SNPOM algorithm is used to process a large number of historical data in train operation to obtain the model order and parameters are optimized, and the accuracy of the model is verified. Then the model is used to construct a predictive controller to track the train target curve online. The results show that the predictive controller based on RBF-ARX model not only solves the problem of dynamic model mismatch, but also is superior to the traditional predictive controller in control accuracy and comfort index.
随着高速列车的快速发展, 由司机人工操作加列车自动防护系统(ATP)保障的控车方式开始暴露出其行车效率低、控车精度低等问题, 控车方式正在逐步向列车自动驾驶系统(ATO)代替人工操作的方向发展, 要实现这一转化, 核心是对列车控制算法和建模方式进行改进, 使列车能在保证行车安全的前提下, 实时追踪列车的目标-距离曲线, 并满足各项性能指标的要求[1].目前, 针对列车控制器的改进方面, 文献[2]设计了一种基于AM-VFF-RLS辨识算法的列车预测控制器, 该控制器具有在线辨识准确的优点, 但辨识算法过于繁琐, 导致在线模型辨识实时性不强.文献[3]提出了基于滑模预测控制的高速列车控制器, 该控制算法在满足一定条件下可以实现精确控车, 但由于控制过程中线路参数固定不变, 导致该控制器鲁棒性不强, 难以适应高速列车复杂多变的运行环境.文献[4]采用单质点、多质点模型作为预测控制器的预测模型, 分别分析了高速列车在SISO和MIMO情况下改进广义预测控制器的控制曲线, 在运行初期有良好的控制效果, 但存在列车模型难建立、计算量大、后期控制器易出现模型失配的问题, 实用性不佳.文献[5]研究了多混合模型下的预测控制器, 建立了基于列车动力学原理的MLD模型作为预测模型实时控制, 但该控制器的多种约束条件使得计算量庞大, 因此不具有实用性.文献[6]提出了基于模糊预测控制的ATO控制器, 使用模糊理论对列车系统建模, 取得了良好的控制效果, 但是建立模糊隶属度函数过于依赖专家经验, 难以应对列车在多变环境下的运行情况并且实时更新隶属函数过于繁琐.目前列车控制器的预测模型多采用单质点、多质点动力学模型以及受控自回归积分滑动平均模型(CARIMA)[7, 8, 9], 但单质点、多质点模型阶数过大, CARIMA模型辨识需要计算Dioaphantine公式, 存在计算数据量大、实时性不强的问题.
针对上述问题, 本文作者设计了一种基于RBF-ARX模型的预测控制器, 利用列车运行过程中ATP系统记录的大量控制器输入、输出数据进行预测模型设计, 摆脱了传统预测控制器依赖动力学模型的现状, 通过与传统预测控制器在控车精度和牵引/制动力曲线平滑度等方面进行对比, 验证了该控制器的可行性.
RBF-ARX模型采用RBF三级神经网络对有外部输入的线性自回归模型(Auto-Regressive eXogenous)进行参数辨识, 结合了RBF神经网络无限逼近和ARX模型全局描述的特点, 在工业领域取得了广泛的应用[10, 11, 12].该模型采用了 “ 黑箱” 原理[13], 因此无需对列车的内部机理进行精确描述, 只需研究列车的输入输出数据, 分析得出和实际列车系统等效的模型, 代替物理模型进行研究, 可以解决列车等非线性系统在复杂环境运行时的建模难题.
RBF神经网络是前馈人工神经网络的一种, 因其设计方法简单且原理上可以近似逼近任意函数的特性, 在数据处理领域得到了广泛的应用[14, 15].本文使用了三级RBF神经网络, 三级神经网络层分别是输入、映射、输出层, 其中输入层到映射层的变换函数是非线性函数, 映射层到输出层的转换函数采用径向基函数, 一般情况下径向基函数选用高斯函数.
ARX模型是时间序列模型中的一种, 融合了自回归模型(AR)和移动平均模型(MA)的结构[16], 能够准确反映列车的输入输出特性, 具有结构清晰、全局描述的特点, 如下
式中:
鉴于直接求解参数的过程比较繁琐且耗时, 利用神经网络计算速度快, 自学习能力强的优点, 使用RBF神经网络的中心点表示式(1)中的参数, 得到的结构即为RBF-ARX模型, 如下
$\begin{array}{* {35}{l}} y\left( t \right)={{\varphi }_{0}}\left[\widetilde{\omega} \left( t-1 \right) \right]+\overset{p}{\mathop{\underset{i=1}{\mathop{\sum }}\, }}\, {{\varphi }_{y, i}}\left[\widetilde{\omega} \left( t-1 \right) \right] \\ y\left( t-1 \right)+\overset{q}{\mathop{\underset{j=1}{\mathop{\sum }}\, }}\, {{\varphi }_{u, j}}\left[\widetilde{\omega} \left( t-1 \right) \right]u\left( t-1 \right)+\xi \left( t \right) \\ \end{array}$(2)
其中相关参数如下
${{\varphi }_{0}}\left[\widetilde{\omega} \left( t-1 \right) \right]=c_{0}^{0}+\overset{m}{\mathop{\underset{k=1}{\mathop{\sum }}\, }}\, c_{k}^{0}\text{exp}\left[ -\lambda _{k}^{y}\|\widetilde{\omega} \left( t-1 \right)-\text{z}_{\text{k}}^{\text{y}}{{\|}^{2}} \right], \\ {{\varphi }_{y, i}}\left[\widetilde{\omega} \left( t-1 \right) \right]=c_{i, 0}^{y}+\overset{m}{\mathop{\underset{k=1}{\mathop{\sum }}\, }}\, c_{i, k}^{y}\text{exp}\left[ -\lambda _{k}^{y}\|\widetilde{\omega} \left( t-1 \right)-\text{z}_{\text{k}}^{\text{y}}{{\|}^{2}} \right], \\ {{\varphi }_{u, i}}\left[\widetilde{\omega} \left( t-1 \right) \right]=c_{j, 0}^{u}+\overset{m}{\mathop{\underset{k=1}{\mathop{\sum }}\, }}\, c_{i, k}^{u}\text{exp}\left[ -\lambda _{k}^{u}\|\widetilde{\omega} \left( t-1 \right)-\text{z}_{\text{k}}^{\text{u}}{{\|}^{2}} \right], \\ \widetilde{\omega} \left( t-1 \right)={{\left[\widetilde{\omega} \left( t-1 \right), \widetilde{\omega} \left( t-2 \right), \ldots , \widetilde{\omega} \left( t-d \right) \right]}^{\text{T}}}, \\ z_{k}^{h}={{\left( z_{k, 1}^{h}, z_{k, 2}^{h}, \ldots , z_{k, d}^{h} \right)}^{\text{T}}}; h=y, u$(3)
式中: $\widetilde{\omega}(t-1)$表示列车控制器的状态向量, 选择系列车控制器的输出值作为特征向量;
建立列车RBF-ARX模型后, 需要使用列车历史正常运行数据离线辨识模型阶数及模型参数, 为了保障模型尽可能地逼近列车实际运行的状态, 本文使用SNPOM算法求解并优化模型阶数及模型参数[17], 该算法全称为结构化的非线性参数辨识, 是基于线性最小二乘法(LSM)和列维布格-马奎尔特(LMM)两种算法提出的, 具有收敛速度快、预测误差小的特点, 具体求解流程如图1所示.
1)首先对式(2)中的参数进行分类, 分为线性
使用式(4)对式(2)进行简化, 可以方便后续计算, 结果如下
$y\left( t \right)=\varphi {{[{{\text{ }\!\!\theta\!\!\text{ }}_{N}}, \widetilde{\omega}\left( t-1 \right)]}^{\text{T}}}{{\theta }_{\text{L}}}+\xi \left( t \right)$ (5)
2)判定模型阶数, 即判定
式中: V表示辨识数据集的均方误差; N为辨识数据集的总数, 并且满足N> > max(p, q); k是参数总数, 求解(p, q)时采用不同的(p, q)组合依次按照AIC准则, 循环计算出不同阶次组合对应的AIC值, 并据最小AID值确定
3)设定目标函数
式中:y(t+1|t)表示列车预测控制向前一步预测的输出值; F($θ _L$, $θ _N$)表示列车RBF-ARX模型的预测值与列车真实数据之间的差.
4)参数优化.
优化参数的实质为求解使得目标函数
式中:
5)判断条件
预测控制通过利用其预测输出的特性来设计控制器对输出进行优化, 因此相比传统的控制方法, 具有控制精度高、计算量小的优点[18, 19], 目前已成功应用到列车领域.
预测控制通常具有预测模型、滚动优化、反馈矫正3个环节, 本文设计RBF-ARX模型作为预测模型, 具体结构如图2所示.
图2中输入变量是列车的牵引/制动力
高速列车一般采用“ 牵引-匀速-惰性-制动” 的多工况运行机制, 本文假定高速列车目标速度曲线如式(10), 作为预测控制器的参考轨迹.
式(2)中RBF-ARX模型得到的是差分形式, 而控制器一般处理状态空间形式, 因此需要将差分形式转换成状态空间形式.定义状态变量
式中:A、B、C、D的具体求法可以参照文献[21].
本文使用了从济南到泰安段高速列车正常运行过程中, ATP设备记录的列车1 000组输入输出数据, 在Matlab/Simulink平台上进行仿真分析, 选择列车型号为CRH380A, 其主要参数设置为:动力配置6M2T, 空车质量为388.64 t, 最高运营速度为400 km/h, 受流电压为AC25 kV, 50 Hz,
本文选择能够表征模型状态的输出量, 即列车速度作为模型的状态变量.高速列车是单输入单输出系统, 根据经验选择
从图3(a)中可以明显看出RBF-ARX模型的输出和实际输出基本重叠, 图3(b)表明模型的误差在[-0.2, 0.2]之间波动, 可以认为是由噪声干扰造成的, 整体来说,
结果表明该模型的误差范围位于[-0.23, 0.25]之间, 满足系统辨识的目标要求, 极大程度地逼近了复杂环境下高速列车运行的非线性模型, 因此选择该模型作为预测控制器的预测模型, 能够提高控制器的控制精度以及多变环境下的鲁棒性.
本文在Simulink平台下搭建基于RBF-ARX模型预测控制器和传统控制器的仿真模型, 为了反映相同条件下两种控制器的性能差别, 设置两控制器的参数都为:最大预测步长P=18, 控制步长M=9, 控制加权系数λ =0.05, 柔化系数α =0.4, 阶梯因子β =0.05, 控制器的控制对象是列车, 为单输入单输出系统, 其输入量为列车牵引/制动力, 输出控制量为列车实时速度, 预测控制器的参考轨迹如式(10), 对得到的RBF-ARX模型方程进行转换, 得到的3阶状态空间方程参数如下:
得到控制器的状态空间方程后, 按照图2所示结构建立控制器, 传统控制器的预测模型选择单质点动力学模型, 分别对两种控制器进行目标-速度曲线的追踪仿真, 设置列车运行线路为平直区段, 线路长度为250 km, 总运行时长为1 000 s, 模拟高速列车在复杂环境运行的速度控制, 比较了基于RBF-ARX模型的预测控制和传统的预测控制器的追踪效果, 结果如图5所示.
图5中红线表示列车的目标-距离追踪曲线, 而黑线表示基于RBF-ARX模型的预测控制器的速度追踪曲线, 蓝线表示的是传统预测控制器的追踪曲线, 从图5中可以看出, 列车高速运行在复杂环境下, 传统预测控制器的追踪效果不佳, 在873 ms处甚至因为模型失配会导致列车紧急制动, 造成极大的损失, 危及行车安全, 相比之下, 基于RBF-ARX模型的预测控制器, 不仅能够实时追踪列车的运行, 且节能效果良好, 不存在紧急制动的隐患, 为了定量描述, 分别计算两种控制器仿真的速度跟踪误差的统计量, 如表1所示.
![]() | 表1 速度追踪误差统计 Tab.1 Speed tracking error statistics |
如表1所示, 传统控制器的均值为-2.365 1、方差为15.642 10, 远大于基于RBF-ARX模型的预测控制器误差的均值和方差, 均值和方差越靠近0值表示控制器的追踪效果越好, 和目标-距离曲线之间的误差越小, 可以得出基于RBF-ARX模型的预测控制器在曲线追踪方面比传统控制器有所改善.
列车在运行过程中乘客所感受到的舒适度也是评价列车控制效果的重要指标, 本文使用控制力的变化率为
式中:γ 为变化率, m/s3;
分别求解两种控制器控制力的
1)改变了传统预测控制器以物理模型作为控制器预测模型的做法, 从列车运行过程中产生的输入、输出数据入手, 设计了一种基于数据驱动的高速列车预测控制器.
2)利用从现场采集的真实数据, 进行RBF-ARX模型的构建和模型参数、阶次的优化求解, 并且对得到模型的可行性进行了验证, 结果证明了该模型具有可行、误差小的特点.
3)使用得到的RBF-ARX模型构建了新型预测控制器, 在相同控制参数和目标-距离曲线下, 与传统的预测控制器进行比较, 分析了两类控制器的相关性能差异, 仿真结果表明基于RBF-ARX模型的预测控制器具有追踪误差小、舒适度高的特点.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|