一种基于振动加速度信号辨识最大压力升高率的方法

摘要

本发明公开了一种基于振动加速度信号辨识最大压力升高率的方法，包括以下步骤：步骤一，利用切比雪夫低通滤波器滤除振动加速度信号中高频干扰；步骤二，利用经验模态分解算法去除振动加速度信号中的低频干扰；步骤三，利用处理曲线上从燃烧始点到最大压力升高率出现时刻内的面积表征最大压力升高率。通过低通滤波器去除实测加速度信号的高频干扰，利用EMD方法去除其中的低频干扰，有效的提高了实测信号的信噪比。本发明利用处理曲线上从燃烧始点到最大压力升高率出现时刻间的面积表征最大压力升高率。本发明提出了建立的最大压力升高率与IMF面积间的线性关系，基于该线性关系对其辨识结果进行修正。

说明书

技术领域

本发明涉及信号处理技术领域，尤其涉及一种基于振动加速度信号辨识最大压力升高率的方法。

背景技术

最大压力升高率作为燃烧过程控制中的重要特征参数，一般通过对实测的振动加速度信号进行积分来获取。常用的时域积分方法有梯形公式积分法和Simpson积分法，梯形公式积分法算法简单，所需计算时间较短，但精度较差，且收敛速度缓慢；Simpson积分法相对于梯形公式来说，计算精确高，然而预处理后的振动加速度信号中仍存在干扰，尤其是燃烧始点附近，导致积分的初始值存在偏差，影响后续的积分值，且Simpson积分法的计算效率较低。

发明内容

本发明的目的就是为了解决上述问题，提供一种基于振动加速度信号辨识最大压力升高率的方法，利用缸盖表面振动加速度速度信号辨识最大压力升高率，通过低通滤波器去除实测加速度信号的高频干扰，利用EMD方法去除其中的低频干扰，有效的提高了实测信号的信噪比。

为了实现上述目的，本发明采用如下技术方案：

一种基于振动加速度信号辨识最大压力升高率的方法，包括以下步骤：

步骤一，利用切比雪夫低通滤波器滤除振动加速度信号中高频干扰；

步骤二，利用经验模态分解(EMD)算法去除振动加速度信号中的低频干扰；

步骤三，利用处理曲线上从燃烧始点到最大压力升高率出现时刻内的面积表征最大压力升高率。

所述步骤一中切比雪夫低通滤波器纹波系数为0.001，阶次为10。

所述步骤一中滤除振动加速度信号中150倍基频以上的与燃烧激励响应信号无关的高频干扰。

所述步骤三中的处理曲线是计算各阶本征模式函数与缸内压力信号二次导数的相关性，将前两阶本征模式函数求和后组成处理曲线。

本发明的有益效果：

通过低通滤波器去除实测加速度信号的高频干扰，利用EMD方法去除其中的低频干扰，有效的提高了实测信号的信噪比。

本发明利用处理曲线上从燃烧始点到最大压力升高率出现时刻间的面积表征最大压力升高率。

本发明提出了建立的最大压力升高率与IMF面积间的线性关系，基于该线性关系对其辨识结果进行修正。

本发明对利用振动信号对缸内燃烧过程特征参数的研究，评价缸内燃烧状态，对内燃机燃烧过程在线监测、闭环控制及故障诊断具有重要的理论意义及实用价值。

附图说明

图1为缸压、压力升高率和缸压二次导数曲线；

图2为处理曲线与缸压二次导数曲线；

图3(a)为800r/min识别结果；图3(b)为1000r/min识别结果；

图3(c)为1200r/min识别结果；图3(d)为1400r/min识别结果。

具体实施方式

下面结合附图与实施例对本发明作进一步说明。

一种基于振动加速度信号辨识最大压力升高率的方法，利用切比雪夫低通滤波器滤除高频干扰,利用EMD经验模态分解去除振动加速度信号中的低频干扰，利用处理曲线上从燃烧始点到最大压力升高率出现时刻内的面积表征最大压力升高率。

上述切比雪夫低通滤波器纹波系数为0.001，阶次为10，滤除振动加速度信号中150倍基频以上的与燃烧激励响应信号无关的高频干扰。

上述EMD是经验模式分解(Empirical Mode Decomposition，EMD)，EMD是根据非线性、非平稳信号本身的特征时间尺度将信号分解成一组稳态和线性的数据序列集，即本征模式函数(IMF)。IMF反映信号的内部特征，残余分量表示信号的趋势，需满足以下两个条件：整个信号中零点数与极点数相等或者至多相差1；对于信号上任意一点，由局部极大值确定的包络线和局部极小值确定的包络线的均值均为0。其具有识别局部特征的优势。

上述处理曲线是计算各阶IMF与缸内压力信号二次导数的相关性，前两阶IMF与平均缸内压力信号二次导数相关性较高，将前两阶IMF求和后组成处理曲线。

图1显示了缸压、压力升高率曲线和缸内压力信号二次导数曲线，内燃机中常将压力升高率在压缩冲程中的折点作为燃烧始点，在缸压二次导数曲线中，燃烧始点对应峰值前的第一个过零点，即A点，B对应的是最大压力升高率出现时刻。图2显示了一个循环的缸内压力信号二次导数和经过EMD分解后得到的IMF1与IMF2之和的振动加速度曲线(定义为处理曲线)，前两个的圆圈即是图1中的A、B两点。与之相对应的，处理曲线中的前两个三角形分别用于表征燃烧始点出现时刻、最大压力升高率出现时刻。

故本发明利用振动加速度曲线上从燃烧始点到最大压力升高率出现时刻内的面积，即图2中所示阴影面积(定义为IMF面积)来表征最大压力升高率。图3(a)-图3(d)给出了不同工况下最大压力升高率与IMF面积，为了便于分析，图中同时给出了利用线性回归方法离线建立的最大压力升高率与IMF面积间的线性关系。由图可见，不同工况下，拟合度分别为0.97、0.97、0.91、0.91，二者之间有良好的线性关系。故实际应用过程中，可基于处理后的振动加速度信号计算相应的IMF面积，结合离线标定的不同工况下的两者之间的关系，进而估算实际工况下的最大压力升高率。其中，800r/min、1000r/min、1200r/min,10N·m及1400r/min,10N·m工况的最大偏差分别为2.2％、0.72％、3.41％及3.97％，表明IMF面积可以用于表征最大压力升高率。

上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。