CN107776559A

CN107776559A - 基于趋势补偿的汽车防抱死系统轮速信号处理的方法

Info

Publication number: CN107776559A
Application number: CN201711044178.6A
Authority: CN
Inventors: 邹浙湘; 温开元; 黄宝山; 李兵; 吴光德; 郭佳韵
Original assignee: Beijing Institute of Technology Zhuhai
Current assignee: Beijing Institute of Technology Zhuhai
Priority date: 2017-10-31
Filing date: 2017-10-31
Publication date: 2018-03-09
Anticipated expiration: 2037-10-31
Also published as: CN107776559B

Abstract

本发明公开并提供了一种联合中值均值加权和经验模函数分解来滤除温漂干扰并无损轮速信息的基于趋势补偿的汽车防抱死系统轮速信号处理的方法。本发明中所提供的算法避免了传统滤波器中由于固定截止频率而带来的误差。MEM滤波器以一种非线性的方式将原始信号的低频部分滤出，EMD方法将MEM的输出分解成一系列IMF。在这些IMF中，t检验会选出温度漂移干扰的成分。通过上述方法，该算法能够有效滤除温度漂移干扰，并能保存好轮速信号的有用信息。从而实现高精度的汽车制动控制，减少发动机运行故障，提高车辆运载效能。本发明适用于汽车控制领域。

Description

基于趋势补偿的汽车防抱死系统轮速信号处理的方法

技术领域

本发明涉及一种基于趋势补偿的汽车防抱死系统轮速信号处理的方法。

背景技术

混合动力汽车是目前汽车发展的新趋势，其具有绿色环保、节能高效等诸多优点。然而，多个能源提供动力也带来了汽车电气控制复杂，发动机协调控制难度以及故障增加的问题。尤其是混合动力公交等公共交通工具，因启动和停止的次数频繁，对于汽车防抱死系统（Anti-lock Braking System，ABS）的制动精度及准确度均有更高的要求。因此，有效提高ABS系统的信号检测水平，减小或者抑制ABS的控制误差，成为了汽车控制领域的新研究热点^[6]。

现有的ABS系统一般均是通过对轮速信号的获取和处理，来得到进行反馈控制的信号。常见的用于拾取轮速信号的传感器主要包括电涡流感测型、磁电型、光电型以及霍尔型等。其中，电涡流及磁电型的传感器在检测精度上相对粗糙，检测信号存在迟滞或者延后的特点；光电型传感器在检测灵敏度方面虽然最高，却对于使用环境的要求比较严格，污垢以及外界光强均会成为检测精度的影响因素，其在车辆运行的复杂环境下干扰因素过多；霍尔传感器^[8]作为一种高灵敏、高精度的轮速传感器，其相对于磁电类的传感器具有更突出的检测性能，同时对于外界环境影响又相对鲁棒，是一种更加适于混合动力汽车应用的轮速检测传感器。

然而，霍尔传感器在使用过程中，也有其自身的特点。环境温度变化会引起霍尔传感器检测信号出现温度漂移，具体体现在轮速检测信号上为一个低频的趋势波动信号。本文针对混合动力汽车的ABS霍尔传感器进行温度趋势估计和补偿，通过对霍尔传感器实时检测得到的轮速信号进行滤波，提取出联合中值均值加权(MEM)和经验模函数分解(EMD)表征温度漂移部分的趋势信号。并从轮速信号^[14]中进行消除，从而实现高精度的汽车制动控制，减少发动机运行故障，提高车辆运载效能。

发明内容

本发明所要解决的技术问题是克服现有技术的不足，提供了一种联合中值均值加权和经验模函数分解来滤除温漂干扰并无损轮速信息的基于趋势补偿的汽车防抱死系统轮速信号处理的方法。

本发明所采用的技术方案是：本发明包括以下步骤：

A.以原信号的首尾对原信号进行延拓处理，可以根据如下方程进行延拓：

其中X(t)为滤波的轮速信号，X(0)为滤波的轮速信号的首值，X(L-1)为滤波的轮速信号的尾值，L为长度为，W为该滑窗的长度，为长度为W的窗滑过延拓后的信号；

B.对每个窗内的向量进行中值和均值加权估计，该估计值为中值跟均值的一个凸组合

，

其中，和表示的是当前窗里向量的均值与中值，是权重系数，则具体中值均值加权的处理是得到下列方程的数值解：

，

其中

，

是示性函数，在＞0时取1，在＜0时取-1；

C.将步骤b中所得的所有估计值组成一个新的向量，作为温漂干扰的初步估计，对于该初步估计的初步修正的方法如下：

，

其中

，

表示的是迭代n次之后得到的数值解，表示的是精确解；

D.对步骤C中经初步修正后的温漂干扰进行经验模函数分解，具体步骤如下：

（1）确定的所有局部极大值点和局部极小值点；

（2）通过至少为三次样条插值分别拟合局部极大值点和局部极小值点，得到上下两个包络线；

（3）从中减去（2）里上下包络线的均值，得；

（4）将作为一个新的重复上述2和3得到，直到SD值介于0.2~0.3之间。其中，SD值通过下式计算：

；

（5）将（4）中得到的定义为第一个固态模函数，从中减去该固态模函数得到第一个残基：

；

（6）将残基看成一个新，重复上述所有步骤，得到一系列和 … ，直到最后的残基是一个常量或单调直线或单极值点的函数，通常将最后的残基看成最后一阶固态模函数，故整个经验模函数的分解相当于完成了如下方程的替换：

；

E.根据中值均值加权的固有缺陷，采用采取从高阶固态模函数到低阶固态模函数进行部分求和值的t检验对步骤C中的固态模函数进行筛选，t检验的两个假设如下：

vs ；

F.将筛选出的固态模函数重构后，从中减去即可得到较准确的温漂干扰；再从原轮速信号中减去得到的准确温漂干扰，即可得到准确的轮速信号，具体如下式：

，

其中为中值均值加权引入的误差，为最终修正后的轮速信号。

在步骤B中，k的取值范围为[1.14, 1.95]，每次数值求解的迭代次数至少为30次。

在步骤C以及D中，当滑窗滑过整个数据信号后得到对应数量的估计值，组合在一起就是对原信号中温漂干扰的估计。

本发明提出了一种利用联合中值均值加权和经验模函数分解估计温度漂移趋势信号的算法。其首先通过联合中值均值加权估计出温度漂移趋势成分后，这一估计过程会引入一个零均值的近似正态分布噪声，表现在频域为高频噪声，因此再对估计温度漂移趋势进行自适应固态模函数分解。使用t检验的方法，对从低向高阶固态模函数进行累加判断检验，判断出各阶固态模函数中不属于温度漂移趋势的成分，再从联合中值均值加权估计的基线信息中剔除，继而得到温度漂移趋势的精确估计。最后从轮速信号中进行消除，从而实现高精度的汽车制动控制，减少发动机运行故障，提高车辆运载效能。

附图说明

图1是基于实验采集的ABS制动时标准数据库的本文算法跟形态学滤波方法的对比图；

图2是温度漂移修正轮速信号对比图；

图3是温度漂移趋势信号图。

具体实施方式

本发明包括以下步骤：

A. 以原信号的首尾对原信号进行延拓处理。假设需要滤波的轮速信号X(t)的长度为L，利用其首值X(0)以及尾值X(L-1)进行延拓处理。因此，后续步骤2中的滑窗长度设定为该采样率值在0.3~0.6之间。设定该滑窗的长度为W，则可以根据如下方程进行延拓

；

B. 以长度为W的窗滑过延拓后的信号，对每个窗内的向量进行中值和均值加权估计。该估计值可近似的理解为中值跟均值的一个凸组合

，

其中，和表示的是当前窗里向量的均值与中值，是权重系数。则具体MEM的处理是得到下列方程的数值解

，

其中

，

是示性函数，在＞0时取1，在＜0时取-1。通常，k的取值范围为[1.14,1.95]。本例算法中取k值为1.5。每次数值求解的迭代次数至少为30次，这样才能保证最终的解收敛到可以承受的求解精度内。本例算法的迭代次数n为50次。

C. 所有估计值组成一个新的向量，作为温漂干扰的一个初步估计。并对其进行EMD分解，得到一系列固态模函数（IMF）。当滑窗滑过整个数据信号后得到对应数量的估计值，组合在一起就是对原信号中温漂干扰的估计。需要补充说明的是，每个估计值对精确值是有误差的，这个误差反映了MEM方法的固有缺陷，且有一定的统计规律，即服从一个0均值的渐近正态分布，从而方便了算法的进一步修正。具体表述如下：

,

其中

,

表示的是迭代n次之后得到的数值解，表示的是精确解。

D. 为了进一步的修正这部分的误差，对滤出的进行EMD分解。具体步骤如下：

（1）确定的所有局部极大值点和局部极小值点；

（2）通过3次样条插值分别拟合局部极大值点和局部极小值点，得到上下两个包络线；

（3）从中减去（2）里上下包络线的均值，得；

（4）将作为一个新的重复上述2和3得到，直到SD值介于0.2~0.3之间。其中，SD值通过下式计算：；

（5）将（4）中得到的定义为第一个IMF，从中减去该IMF得到第一个残基

；

（6）将残基看成一个新，重复上述所有步骤，得到一系列和 … ，直到最后的残基是一个常量或单调直线或单极值点的函数。通常将最后的残基看成最后一阶IMF，故整个EMD的分解相当于完成了如下方程的替换

。

E.根据MEM的固有缺陷，采用t检验对步骤3中的IMF进行筛选。由于MEM的固有误差混在中，所以需要对步骤3中最后等式右边的各IMF进行筛选。EMD分解后得到的各个IMF从低阶到高阶有如下频率分布规律：低阶IMF分量含有较多的高频分量和较少的低频分量。在频域中，对于MEM引入的误差分布在比温漂干扰频率高的频带。另外，这部分误差在时域服从一个0均值的渐近正态分布。因此，本算法采取从高阶IMF到低阶IMF进行部分求和值的t检验。该检验是为了筛选出哪些阶的IMF组合为MEM引入的误差。t检验的两个假设如下

vs 。

F. 将筛选出的IMF重构后，从中减去即可得到较准确的温漂干扰。再从原轮速信号中减去得到的准确温漂干扰，即可得到准确的轮速信号。假设在t检验过程中，在第P阶接受了假设，则部分和即为MEM引入的误差。这部分的和需要从中减去，这样便得到了最终修正后的：

。

需要补充说明的是，温漂干扰本身也有可能是0均值。故在迭代t检验时，本算法会在当前阶等于IMF总个数的一半时强行终止。在这种特殊情况下会引入部分误差，但误差不会过大，最后从原信号中减去得到修正后的轮速信号。

基于本发明所描述的算法所做出的实验及其实验结果：

轮速信号取自实验采集的ABS制动时标准数据库，其采样率为360 Hz。将10 mV的振幅范围进行11位的离散化处理，振幅在0~2 047间，其中1 024对应0值。实验样本截取3 000点数据，相应地MEM滑动窗口长度L在120~200点之间。另外，为了进一步测试算法的有效性，通过低通一个随机序列来产生一系列的人工噪声。随机序列的振幅范围区间为[0, c]，服从均匀分布。图2为受人工噪声（c=500）污染的信号和修正后的信号，可看出所提算法即使在存在明显温度漂移的情况下，仍能有效恢复原始轮速信号。

图2中虚线为受温度漂移干扰的轮速信号，实线为修正后的信号。图3进一步给出了，估计得到的温度漂移趋势。本文通过一个标准来评价算法有效性，并计算信噪比SNR。数值的计算方法为。给定一个向量X，表示该向量的标准差。横向对比了本文算法和基于形态学滤波的方法，表1列出了在不同振幅c值下的平均SNR提高值，本文所提出算法的提高值最高可达21.84 dB。

实验结果表明，本文所提算法能够有效估计温度漂移干扰，从而消除ABS制动系统的轮速监测误差。因为所提出的算法包含了一个自修正的过程，从而能更加合理、精确地选择出温度漂移干扰的分量，由此便可提高整体的滤出效果。

作为一个完全非线性的滤除温度漂移干扰方式，本发明中所提供的算法避免了传统滤波器中由于固定截止频率而带来的误差。MEM滤波器以一种非线性的方式将原始信号的低频部分滤出，EMD方法将MEM的输出分解成一系列IMF。在这些IMF中，t检验会选出温度漂移干扰的成分。通过上述方法，该算法能够有效滤除温度漂移干扰，并能保存好轮速信号的有用信息。从而实现高精度的汽车制动控制，减少发动机运行故障，提高车辆运载效能。

本发明适用于汽车控制领域。

Claims

1.一种基于趋势补偿的汽车防抱死系统轮速信号处理的方法，其特征在于：它包括以下步骤：

其中X(t)为滤波的轮速信号，X(0)为滤波的轮速信号的首值，X(L-1)为滤波的轮速信号的尾值，L为长度为，W为该滑窗的长度，为长度为W的窗滑过延拓后的信号；

,

，

其中

，

是示性函数，在＞0时取1，在＜0时取-1；

，

其中

，

表示的是迭代n次之后得到的数值解，表示的是精确解；

（1）确定的所有局部极大值点和局部极小值点；

（3）从中减去（2）里上下包络线的均值，得

；

（4）将作为一个新的重复上述2和3得到，直到SD值介于0.2~0.3之间，其中，SD值通过下式计算：

；

（6）将残基看成一个新，重复上述所有步骤，得到一系列和…，直到最后的残基是一个常量或单调直线或单极值点的函数，通常将最后的残基看成最后一阶固态模函数，故整个经验模函数的分解相当于完成了如下方程的替换：

；

vs ；

，

2.根据权利要求1所述的基于趋势补偿的汽车防抱死系统轮速信号处理的方法，其特征在于：在步骤B中，k的取值范围为[1.14, 1.95]，每次数值求解的迭代次数至少为30次。

3.根据权利要求1所述的基于趋势补偿的汽车防抱死系统轮速信号处理的方法，其特征在于：在步骤C以及D中，当滑窗滑过整个数据信号后得到对应数量的估计值，组合在一起就是对原信号中温漂干扰的估计。