利用可控悬架对车辆行驶性能进行优化控制时, 悬架系统状态的精确程度直接影响可控悬架的控制效果, 但是由于技术与成本的限制, 有些控制参数无法准确测量或传感器价格昂贵, 并且可靠性得不到保证, 难以获取所需的状态参数.基于上述情况, 需要搭建合适的观测器对悬架系统进行状态估计.考虑到典型路面干扰产生的随机激励所引起的垂向车辆动力学的随机性, 多采用卡尔曼滤波来进行车辆的状态估计[1-2].但是传统的卡尔曼滤波器有严格的条件假设:模型应该是精确的, 统计噪声应该是高斯白噪声.因此, 仅采用在不同工作条件下的传统卡尔曼滤波器难以得到更高精度的状态估计.
目前普遍使用的多模型估计方法是基于多模型的卡尔曼滤波算法[3],例如:伪贝叶斯算法[4]、交互式多模型(IMM)算法[5]和交互式多模型自适应估计(MMAE)[6].交互式多模型算法已经被证实可以以较低计算成本实现良好的计算效果[7].交互式多模型算法通常将模式转移建模为与连续状态无关的恒定模式概率转移的马尔科夫过程; 在很多情况下, 当连续状态位于连续状态空间的子集内时就会发生模式转移[8].
在基于标称系统模型进行状态观测估计时, 不确定系统参数可能造成估计质量的大幅降低.状态估计的误差就会变大.因此, 本文使用递归最小二乘法(RLS)滤波来进行参数识别, 并且这种方法已被证实是一种有效的参数识别算法[9].文献[10]分析了一种用于参数识别的时变系统的并行自动谐调递归最小二乘算法; 文献[11]研究了一种引入遗忘因子的不稳定递归最小二乘算法与卡尔曼滤波器相结合来估计低频机电模型的观测器.
根据上述分析, 本文通过考虑不同模型之间的模式转移, 提出了一种新的基于交互式多模型卡尔曼滤波与递归最小二乘算法耦合的状态估计方法.
1 道路激励模型本研究假定路面不平度是一种均匀的各向同性的高斯随机过程, 它的统计特征可以用功率谱密度(PSD)来描述.根据ISO 8601[12], 路面功率谱密度定义为
(1) |
其中:n(单位为m-1)为空间频率, n0=0.1 m-1是参考空间频率; Gq(n0)(单位为m-3)为参考频率n0下的路面功率谱密度值, 称为路面不平度系数; w为频率指数, 它决定路面功率谱密度的频率结构, 通常取w=2[13].
考虑到角频率ω, 式(1)可写成
(2) |
基于频谱分解定理, G(ω)在功率谱密度函数稳定随机控制过程中起作用.式(2)因此可以表示成式(3)和式(4).
(3) |
(4) |
根据路面生成的时间频率公式[14]:
(5) |
其中σ是白噪声信号的系数.
式(4)和式(5)可以表示成式(6)和式(7).
(6) |
(7) |
(8) |
四分之一线性车辆模型如图 1所示.
模型动力学方程可表示为
(9) |
其中:Fd表示可控阻尼力; ks, kt表示悬架弹簧刚度和轮胎弹簧刚度; xb,
系统的状态向量和输出向量为
其中, 4种状态分别是频响空间、轮胎侧偏、簧上质量速度和簧下质量速度.
因此, 状态空间变量的状态空间方程可写为式(10):
(10) |
其中:
其中:c表示悬架阻尼; w和v表示假定独立的过程高斯噪声和观测高斯噪声;
状态空间离散表达式为
(11) |
式中:k表示离散时间的采样时刻; T表示采样时间.
2.2 广义车辆悬架模型分析悬架的非线性特性主要由弹簧特性、递减非对称的悬架阻尼力:速度的特性、库仑摩擦和悬架支点与车轮平面的运动传递系数产生.这些非线性特性可由相关的非线性阻尼c(k)和变化的弹簧质量mb(k)导出.基于图 1采用一系列不同的附加载荷可以推导出车辆悬架的簧上质量运动方程mb(k)=mb(k-1)+Δm(k).
本研究使用马尔科夫线性四分之一车辆模型.因此, 广义的车辆四分之一状态矢量空间模型可由式(12)表示.
(12) |
其中:
其中, f1i和f2i分别表示第i个不变的最小区域空间和最大区域空间的组合.对于路面速度来说, 未知的干扰可以用高斯白噪声过程来近似.模型的不确定性部分也可以假设成高斯白噪声.每一个线性子模型i都可依据表 1所示的数值结合状态空间表达式联立表示.
递归最小二乘算法是根据已知数据来估计系统参数.假设系统为
(13) |
其中:Θ表示输入矩阵; Y表示输出矩阵; θ表示待识别参数矩阵.
本文提出了可以变换为递归最小二乘法的应用程序.k时刻的簧上质量可以写为
(14) |
其中, Θ(k)定义为在第k步中可测量的数据因子.y(k), mb(k)和e(k)分别表示第k步中输出的测量值、被识别的参数和高斯白噪声过程.
基于上述分析, k时刻簧上质量的递归格式可写为
(15) |
联立更新的观测变量和误差协方差矩阵R, 最小二乘公式可写为
(16) |
(17) |
其中, Kk(k)代表增益矩阵.
本文提出了一个时间离散随机线性混合系统, 如式(18)所示:
(18) |
其中:k=1, 2, …; x(k)∈Rn和z(k)∈Rp分别是离散状态向量和观测向量; m(k)∈{1, 2, …, r}是离散状态(或模式); Am(k)和Cm(k)是与k时刻的m(k)一致的状态转移矩阵[13].过程噪声wm(k)(k)和观测噪声vm(k)(k)是互不相关的零均值高斯白噪声.高斯序列:Q(k)= E(wm(k)wm(k)T)和R(k)= E(vm(k)vm(k)T)分别是它们相应的过程噪声方差矩阵和观测噪声协方差矩阵.模式转移即m(1), m(2), …, m(k)是一个马尔科夫链, 可以写成一个连续态非独立的模态转移矩阵:
(19) |
其中,πij(x(k))是一个在连续态x(k)的条件下从i模态到j模态的条件概率.
(20) |
本研究用P[·|·]来表示条件概率密度分布函数(PDF).
所以, 交互式多模型卡尔曼滤波方法步骤如下:
第一步, 为j模态相同等级的模态定义交互式输入.
由前一步的目标状态估计
预测模态j的概率:
(21) |
从i模态到j模态的混合概率:
(22) |
j模态的混合状态估计:
(23) |
j模态的混合协方差估计:
(24) |
其中:Pij是从i模态到j模态的转移概率;μj(k-1)是j模态在t=k-1步的概率.
第二步, 计算j模态的卡尔曼滤波.
为了更新悬架的预测状态估计(
预测从k-1时刻到k时刻的状态改变:
(25) |
其中Aj为系统离散矩阵.
预测从k-1时刻到k时刻的协方差改变:
(26) |
计算k时刻卡尔曼滤波的增益方程:
(27) |
计算卡尔曼滤波方程:
(28) |
计算卡尔曼滤波器的协方差:
(29) |
第三步, 更新j模态的概率.
采用似然函数来更新模型的概率μj(k).j模态的表达式可改写为
(30) |
在此基础上, j模态的概率:
(31) |
第四步, 定义j模态相同等级的交互式输出.
基于j模态的概率, 整体的模态估计可通过加权每一个滤波器的状态估计结果得到.
计算整体的模态估计:
(32) |
计算整体的协方差估计:
(33) |
交互式多模型卡尔曼滤波器的流程图如图 3所示.
在本次仿真中, 用第3节推导出的交互式多模型卡尔曼滤波算法和递归最小二乘算法耦合的状态估计方案来估计悬架系统垂向动力学的广义状态.本研究假设条件为:轮胎不会与地面脱离接触; 并且假定马尔科夫转移矩阵在仿真之前已知; 道路激励是30 km/h的ISO-C级谱.
本研究在仿真中比较上文提到的状态参数交互式多模型卡尔曼滤波联合算法和单一状态的卡尔曼滤波器状态估计性能.计算出状态估计误差的均方根值, 并在表 2中总结了两种状态估计方法的仿真结果.通过比较仿真的状态估计精度, 所有相关信号的性能的变化即均方根误差波动始终保持在19%内.
研究了增加Δmb=200 kg的额外车体负载的情况.这一载荷近似增加了约45%的车体质量, 可以看作是在汽车上坐了8名乘客.图 4阐明了卡尔曼滤波和交互式多模型卡尔曼滤波的相应结果.如表 2和图 4所示, 交互式多模型卡尔曼滤波可以提高多模型悬架的状态估计精度.
卡尔曼滤波和交互式多模型卡尔曼滤波算法性能估计在图 5所示的四分之一车辆悬架的试验台架上进行.在状态估计的试验过程中, 道路激励由一个液压油缸施加在车轮上, 簧上质量和簧下质量的加速度数据来自道路等级ISO-C级激励谱.此试验中, 数量估计的道路激励参考信号是使用在试验设备中的第二部分通过离线计算得到的.
道路模拟机的参数如表 3所示.
在状态估计中, 仿真结果均方根误差和试验数据均方根误差如表 4所示.在表 4中, 仿真结果与试验结果吻合良好, 可见仿真和试验结果之间的最大均方根误差小于21%, 位于合理范围内.
比较试验结果和仿真结果可见, 使用本文提出的算法可以提高悬架系统的状态估计精度.
5 结论1) 利用交互式多模型卡尔曼滤波算法与递归最小二乘算法耦合得到的自适应卡尔曼滤波方法, 可以在有效识别悬架参数的基础上, 精确得到悬架系统状态.仿真结果说明, 通过与传统卡尔曼滤波方法的比较, 所提出的新方法能够获得更高精度的悬架系统状态估计.
2) 仿真与试验结果表明, 在车辆四分之一悬架台架试验中, 当悬架簧载质量发生变化时, 传统卡尔曼滤波方法估计性能会大幅恶化, 所以此工况下可以利用本文所提出的新方法进行悬架实际状态的实时精确估计.通过仿真与试验可知, 此方法的仿真和试验结果之间的状态估计均方根误差小于21%.
[1] |
Yu F, Crolla D.
An optimal self-tuning controller for an active suspension[J]. Vehicle System Dynamics, 1998, 29: 51–65.
DOI:10.1080/00423119808969366 |
[2] |
Rigatos G, Siano P, Pessolano S.
Design of active suspension control system with the use of Kalman filter-based disturbances estimator[J]. Cybernetics and Physics, 2012, 1: 259–275.
|
[3] |
Li X R, Jilkov V.
Survey of maneuvering target tracking.part Ⅴ:multiple-model methods[J]. IEEE Transactions on Aerospace and Electronic Systems, 2005, 41: 1255–1321.
DOI:10.1109/TAES.2005.1561886 |
[4] |
Maybeck P.
Stochastic models, estimation, and control[M]. New York: Academic Press, 1982.
|
[5] |
Blom H, Bar-Shalom Y.
The interacting multiple model algorithm for systems with Markovian switching coefficients[J]. IEEE Transactions on Automatic Control, 1988, 33: 780–783.
DOI:10.1109/9.1299 |
[6] |
Hanlon P, Maybeck P.
Multiple-model adaptive estimation using a residual correlation Kalman filter bank[J]. IEEE Transactions on Aerospace and Electronic Systems, 2000, 36: 393–406.
DOI:10.1109/7.845216 |
[7] |
Bar-Shalom Y, Kirubarajan X R, Li T.
Estimation with applications to tracking and navigation[M]. New York: Wiley, 2001.
|
[8] |
Campo L, Mookerjee P, Bar-Shalom Y.
State estimation for systems with sojourn-time-dependent Markov model switching[J]. IEEE Transactions on Automatic Control, 1991, 36: 238–243.
DOI:10.1109/9.67304 |
[9] |
Zhu Y.
Efficient recursive state estimator for dynamic systems without knowledge of noise covariance[J]. Mechanic Engineering, 2003, 8: 64–66.
|
[10] |
Li J S, Zheng Y J, Lin Z P.
Recursive identification of time-varying systems:self-tuning and matrix RLS algorithms[J]. Systems & Control Letters, 2014, 66: 104–110.
|
[11] |
Moreno I, Messina A R.
Adaptive tracking of system oscillatory modes using an extended RLS algorithm[J]. Electric Power Systems Research, 2014, 114: 28–38.
DOI:10.1016/j.epsr.2014.03.034 |
[12] |
ISO.Mechanical vibration-road surface profiles-reporting of measured data: ISO8601-1995[S].Melbourne: ARRB Group Ltd., 1995.
|
[13] |
余志生.
汽车理论[M]. 北京: 机械工业出版社, 2009: 207-212.
( Yu Zhi-sheng. Auto theory[M]. Beijing: China Machine Press, 2009: 207-212. ) |
[14] |
Zhang G, Fang Z, Chen S, et al.
White noise simulation for road roughness based on power function[J]. Automotive Engineering, 2008, 30: 44–47.
|