基于GS-SVM的发动机异响在线检测平台.pdf

第2 0 期AUTOMOBILE APPLIED TECHNOLOGY设i计研究10.16638/ki.1671-7988.2023.020.009基于GS-SVM的发动机异响在线检测平台杨兴国，余瑶（重庆工商职业学院智能制造与汽车学院，重庆40 152 0）摘要：为了解决人工听诊法进行发动机异响识别时产生的劳动强度大、工作效率低与准确率波动等问题，文章提出一种基于网格搜索-支持向量机（GS-SVM）的发动机异响在线检测技术。该平台主要包括转速监测、信号采集、信号去噪、特征提取和模式识别等功能，LabVIEW软件负责发动机转速监测和信号采集，并将信号传输至MATLAB接口。在MATLAB软件中，首先利用小波相关系数滤波法去除背景噪声；然后分别利用小波包变换和双谱估计提取信号特征，经归一化处理的信号特征作为支持向量机进行模式识别的输入向量；接着选择分类器-支持向量机（C-SVC）和径向基核函数（RBF），并采用改进的网格搜索法优化参数和g，建立分类模型；最后利用训练成熟的模型预测发动机异响类型。经过测试，该方法的准确率在90%以上，具有一定的工程意义。关键词：发动机异响；小波包；双谱估计；GS-SVM；在线检测平台中图分类号：U469.7The Online Detection Platform of Engine Abnormal Noise Based on GS-SVMChongqing Technology and Business Institute,Chongqing 401520,China)Abstract:In order to solve the problems of high labor intensity,low work efficiency andfluctuating accuracy rate in engine abnormal noise recognition using manual auscultation,thispaper proposes an online detection technology for engine abnormal noise based on grid search-support vector machine(GS-SVM).The platform mainly includes functions such as speed monito-ring,signal acquisition,signal denoising,feature extraction,and pattern recognition,and theLabVIEW software is responsible for monitoring engine speed and collecting signals,andtransmitting the signals to the MATLAB interface.Firstly,the wavelet correlation filtering methodis used to remove background noise in MATLAB software;Then,wavelet packet transform andbispectral estimation are used to extract signal features,and the normalized signal features are usedas input vectors for support vector machine pattern recognition;Next,select the classifier-supportvector machine(C-SVC)and radial basis function(RBF),and use an improved grid search method作者简介：杨兴国（198 6 一），男，博士研究生，讲师，研究方向为汽车噪声与振动控制，E-mail:。基金项目：重庆市教育委员会科学技术研究项目（KJQN202004005）；重庆工商职业学院重点科研项目（NDZD2020-02）。文献标识码：A文章编号：16 7 1-7 98 8(2 0 2 3)2 0-39-0 7YANGXingguo,YUYao(College of Intelligent Manufacturing and Automobile,40to optimize the parameters c and g to establish a classification model;Finally,use the trained andmature model to predict the type of abnormal noise.After testing,the accuracy of this method isover 90%,which has a certain engineering significance.Keywords:Engine abnormal noise;Wavelet packet;Bispectrum estimation;GS-SVM;Online dete-ction platform发动机是汽车的核心部件，其声品质很大程度上决定着消费者的购买欲望和驾驶体验。汽车企业在发动机装配工序完成后常采用人工听诊法对其进行异响检测，避免异响发动机流入汽车市场1-2。人工听诊法操作简便，灵活性强，但识别准确率受环境噪声、技术水平和生理状况等因素影响而不稳定，且劳动强度大。国内外学者针对发动机故障诊断开展了大量研究3-41，但研究对象主要为汽车后市场的发动机，很少以车企刚下线的新发动机作为研究对象。杨诚等5利用LabVIEW测试软件和特征频率识别技术搭建了一个摩托车发动机异响检测平台，取得了较好的效果。然而，仅以故障特征频率作为发动机异响识别依据会受噪声或偶然因素的影响，难以保证准确率。因此，寻找能够多方位表征发动机异响特征的新方法势在必行。杨诚等2 提出运用最小均方算法（Least MeanSquare，LM S）和对称点图形（Symmetrized DotPattern，SD P）相结合的方法提取发动机声音信号特征，识别异响类型。贾继德6 提出利用小波对称极坐标法研究发动机异响信号特征。李皓等7 提出运用小波时频偏向干分析识别车内异响声源并做出相应改善。安周鹏等8 提取水电机组信号的小波能量谱及功率谱特征，提高了故障诊断的准确率。王昱翔等9提取发动机声信号的短时能量和功率谱特征，并将其通过梅尔滤波器得到梅尔频率倒谱系数，该方法得到的时频特征能够较好地反应发动机的运转状态。周全10 提出一种基于二分法的变分模态分解方法提取发动机缸盖表面振动信号，弄清了“皖皖”声异响的产生机理。文献11通过提取油气管道振动信号的小波包能量谱和高阶谱特征，建立了一套油气管网安全监测系统；程静等12 利用双谱估计识别风机轴承振动信号中的非线性相位耦合信息，便于诊断轴承故障。此外，包络分析、希尔伯特变换、边频谱、神经网络、自由回归（Auto Regressive,AR）模型等方法也常应用于轴承、齿轮、机械设备故障诊断。为了更全面地描述发动机信号特征，本文采用小波包变换和双谱估计联合提取发动机信号特征，将汽车实用技术其结果作为发动机异响识别的输入向量。常用的异响模式识别方法有统计模式识别、句法模式识别、模糊模式识别、人工神经网络和支持向量机（Support Vector Machine,SVM)，其中后两者在发动机异响模式识别中应用最为广泛13。相对于人工神经网络，SVM具有决策参数少、收敛速度快和鲁棒性好等优点，是一种有坚实理论基础的新颖的小样本学习方法。因此，本文选择SVM作为发动机异响模式识别方法。1在线检测系统的框架1.1在线检测平台设备为了实现发动机异响的在线检测功能，需要设计专用的自动化检测平台，仅对发动机声音信号的采集、特征提取及模式识别展开研究，并未探讨发动机传送结构及控制装置等内容。发动机异响在线检测平台（见图1）主要由检测台架、传声器、信号采集前端、计算机及相关软件和附属部件组成。其中，本文分别采用北京声望声电技术有限公司的自由声场传感器MPA201和美国国家仪器公司的套件NI-USB9234作为信号采集传感器和数据采集前端，它们的性能参数喷油器发动机喷油器图1发发动机异响在线检测平台表1传声器MPA201的性能参数参数数值灵敏度/(mV/Pa)45.7频率/Hz2020.000动态范围/dB(A)16134本底噪声/dB(A)16输出阻抗/250湿度/%095输出接口BNC注：BNC：卡扣配合型连接器（BritishNavalConnector）。2023年信号采计算机集前端第2 0 期如表1、表2 所示。表2 信号采集前端UI-USB9234的性能参数参数数值通道数/个4位分辨率24采样速率/(ks/s)51.2电压/V-55使用温度/-40701.2在线检测平台工艺1.2.1 发动机工况对于大多数发动机异响而言，异响信号强弱与温度关系不大，为了尽量提高发动机异响在线检测效率，本文选择冷态下采集发动机声音信号。通常，发动机负载越大，燃烧噪声能量越高，异响信号容易被淹没，且新发动机普遍有较高的声品质，异响信号常在低负荷下就突显出来，故本文选择在空载下采集发动机声信号。理论研究和实际经验表明，过高的转速会使燃烧噪声和进、排气噪声成为发动机声信号的主要成分，对异响信号产生掩蔽效应，不利于提取发动机异响信号特征14-1。同时，为了与人工听诊的工艺保持一致，本文选择在转速为150 0 17 0 0 r/min下采集发动机信号。根据四冲程发动机的工作特点，除了利用发动机转速传感器和凸轮轴位置传感器提取发动机转速外，喷油器和点火线圈的控制信号也可以提取发动机转速信息。由脉冲信号频率就能计算出发动机转速。点火线圈传感器引线常被包裹，不便于接线，因此，本文采用喷油器信号监测发动机转速。信号采集前端监测发动机喷油器的喷油信号，并传输给计算机，同时通过LabVIEW软件计算发动机转速（见图1）。如果转速偏离设定值，软件自动调整发动机转速，直至转速处于150 0 1700 r/min。同时，LabVIEW软件也承担着发动机信号采集系统的开始、暂停和结束，以及信号分析与图形显示等任务。1.2.2测点选择发动机异响源于内部零件的不正常受迫振动，振动能量经过不同路径向外辐射噪声。为了获取较强的发动机声信号，测点应尽可能靠近辐射噪声能量最大的零件，且安装方便、不干扰发动机正常工作。由于发动机的薄壁零件具有很强的噪声辐射能力，本文选择在发动机汽缸盖（汽车）或者曲轴箱盖（摩托车）附近布置传声器，杨兴国，等：基于GS-SVM的发动机异响在线检测平台1.2.3信号采集与分析流程为了实现发动机异响在线检测功能，专门设计参数数值激励电流/mA2动态/dB102输入阻抗/k2305连接接口BNC规格/mm130903341且传声器断面距离发动机表面约2 0 cm。了检测台架和控制装置，通过软件可以控制检测台架的运行和信号传输，操作简便。利用LabVIEW软件监测发动机转速检测和采集信号，MATLAB软件负责信号处理和图形化显示，LabVIEW软件有MATLAB接口，可以实现信号传输功能。发动机的信号采集、特征提取和模式识别流程如图2所示。启动或设置采暂停检样参数测平台下一循环正常模式识别异常记录异常信息图2 信号采集与分析流程图2发发动机异响信号分析2.1发动机信号的去噪相比于独立分量分析法、经验模式分解法和盲源分离法等传统信号消噪法，小波滤波具有良好的时频局部化性能，即突出信号、削弱噪声的特点，使其成为在机械故障诊断领域应用较为广泛的分析方法之一。小波分析经过40 余年的发展，产生了不同的降噪理论和降噪方法，截至目前，有三种小波降噪法应用广泛，即模极大值法、小波阈值法和小波系数相关法。本文采用之前的小波相关系数滤波法16 去除发动机声信号的背景噪声，相比于模极大值法和阈值法，它的滤波性能稳定，滤波后信噪比高、均方误差小、平滑度高，能够较为准确地还原真实信号。图3为箱体异响机消噪前后时域信号，通过小波相关系数法消噪后，噪声信号基本被去除，真实信号得以重构。是发动机进启动入检测窗发动机否否转速位于等待15001700 r/min是信号结束信分析号采集开始采集信号正常机421.00.50-0.5-1.001.00.50-0.5-1.002.2发动机信号的特征提取2.2.1小波包能量谱相对于连续小波变换，小波包变换不仅细分低频部分，而且对高频部分进行追加分解，得到更加精确的信号估计，具有更高的频率分辨率。小波包系数的大小表示小波包基函数与信号的相似程度，具有能量的量纲8。因此，可以利用小波包变换来表征信号能量在频域的分布情况。由小波包理论可知，不同的小波基函数对信号作小波包变换的结果影响较大。为了更精确地表征发动机信号的能量，结合理论分析和实践经验，选择“db8”小波基函数，并进行3层小波包分解。经过3层小波包变换后，全分析频段（0 10240Hz）被平均分成8 个区域，每个区域占据1280Hz，信号的小波包能量谱定义为(1)n=1式中，j为小波包分解层数，取3；k为频段数；n为频段内的空间位置。由于不同频段的小波包能量谱差异较大，方便后续分析数据，需要对信号的小波包能量谱进行归一化处理，方法如式（2）所示。（2）2Tk-1汽车实用技术图4为经消噪和归一化处理后的发动机声信号的小波包能量谱，4种类型的发动机小波包能量主要集中在低频段。但是，异响发动机在中频段或高频段也有明显的能量分布，比如箱体异响机在第2 一4频段的小波包能量谱较高；右盖异响机在第2 一4频段和第7 频段的小波能量谱较高；而尖叫机在第2 一4频段和第7 一8 频段都有较明显0.10.2t/s(a)消噪前0.10.2t/s(b)消噪后图3箱体异响机的时域信号Tik212023年0.30.40.30.40.50.5的小波包能量分布。不同类型发动机的小波包能量谱分布存在差异，因此，发动机声信号的小波包能量谱可以作为信号特征16。1.00.80.60.40.202.2.2双谱功率谱描述信号能量在频域内的分布情况，它是时域信号自相关函数的傅里叶变换15。与功率谱求解过程类似，而高阶谱是信号高阶累积量的多维傅里叶变换16。双谱具有高阶谱所有特征，且求解过程简单，工程应用最为广泛。设x是均值为0 的3阶平稳随机序列（过程），它的3阶累积量为C3.x=Ex(n)x(n+t)x(x+t2)满足绝对可和条件，即ZZC.,(ti,t2)+0T2对式（3）的3阶累积量作二维傅里叶变换可得B,(w,wa)=ZZC.(1,t,)e-0m+m.)由双谱理论的对称性可知，同时为便于计算，可令式（3）一式（5）的Ti=t2=t，得到随机序列（信号）x(n)三阶累积量的对角变换（简称对角切片），再进行一维傅里叶变换，得到该随机序列的箱体异响机右盖异响机尖叫机12图4发动机信号的小波包能量谱TT234频段5678(3)(4)(5)频率，双谱能量在整个频域内均匀分布16；箱体第2 0 期一4维谱。设时间序列为2x(t)=Zcos(2 f,+p.)+w(n)6i=1式中，w(n)为模拟信号采集现场的背景噪声而加载的高斯白噪声；f包含了6 个频率信息，fi-fs分别为0.52 8 1、0.8 456、1.932 3、1.2 90 9、1.7 6 7 1,且满足f6=f4+fs；P1,2*95均为0,2 元上独立均匀分布的相位参数，且它们也满足关系式96=4+5。由此可见，式（6）包含了一组二次相位耦合频率信息。仿真信号的采样频率为16；样本数为2 56；采样点为512；傅里叶变换的数据长度为10 2 4。图5(a)的时间序列x(t)的功率谱中包含信号所有频率成分，并在其他频率处有较多的噪声能量。图5(b)所为x(t)双谱的对角切片只在f4、f s、f c 出现能量峰值，而在f、J2、处没有图5(a)的能量峰值。这说明信号做双谱估计可以检测到相位耦合频0.50.40.30.20.1002015杨兴国，等：基于GS-SVM的发动机异响在线检测平台0.51.01.5Frequency/Hz(a)功率谱43率。同时，信号作双谱估计后，噪声得到有效消除，这是功率谱所不具备的特点。机械设备发生故障后，会出现偏离高斯过程的非线性信号，这(6)种非线性关系一般为两个频率的和频率或者差频率，有时也表现为基频的各阶倍频率。采用双谱估计可以精确地提取到信号的二次相位耦合频率信息，便于查找故障点和故障模式识别。常用的双谱估计有参数化法和非参数化法两种，其中参数化估计利用AR、移动平均（MovingAverage,MA）、自回归移动平均（Auto Regressiveand MovingAverage,ARMA）等模型进行估计。参数化估计计算精度高，但运算量大而复杂，且模型参数确定困难；而非参数化估计运算量小，计算方法相对简单，计算精度适中，工程应用性更好。目前，非参数化估计主要有两种方法，即间接法和直接法。间接法双谱分析要求首先对随机序列进行均值处理，然后计算均值序列的三阶累积量，最后再对三阶累积量作二维傅里叶变换得到随机序列的双谱估计；而直接法先对随机序列进行均值分段处理，接着进行一次傅里叶变换，最后再将变换结果作双重相关运算，从而得到随机序列的双谱估计结果。文献11,15-16 对发动机声音信号作双谱估计。由双谱估计理论可知，如果信号存在二次相位耦合频率，在耦合频率处出现能量峰值，对应一组和频率或者差频率。图6 为2.2.1节的消噪后的发动机的双谱。正常机没有出现二次相位耦合2.02.53.03.54.0异响机在频率为8 6 0、17 2 0、2 58 0 Hz时，均出现双谱极值，且三个频率构成二次相位和频率；右盖异响机在92 0、132 0、2 2 40、2 56 0、2 7 8 0、5410325000.5图5信号的傅里叶变换1.01.5Frequency/Hz(b)双谱2.02.53.0 3.54.0020004000Frequency/Hz(a)正常机6 00080001000044654338012103.02.52.01.51.00.506543205340Hz均构成二次相位和频率；而尖叫异响机的2 6 0 0、4140、6 7 40 Hz均出现双谱极值，并构成一组二次相位耦合和频率16。因此，对发动机声信号作双谱变换可以提取信号特征。为了提取方便，选择上述12 个二次相位耦合频率为中心频汽车实用技术率，计算其10 Hz区间内的双谱累加值作为该中860心频率的特征值。2.3发动机异响的模式识别2.3.1样本选择经过小波包变换和双谱估计均可以提取发动机信号特征，为了体现信号特征的多样性，本文1720采用两种方法分别提取发动机信号特征，小波系数相关法去除发动机声信号的背景噪声16，小波2.580包变换和双谱估计提取信号特征，得到8 个小波包能量谱特征值和12 个双谱特征值，构成一个有20个元素的特征向量。2.0004000Frequency/Hz(b)箱体异响机22409202560132027802.0004000Frequency/Hz(c)右盖异响机27203.7801351014.1402.6002000图6 发动机声信号的双谱估计2023年6 0008 0001000053406.000674040006 000Frequency/Hz(d)尖叫机4种类型发动机各取30 个样本，总共构成12 0个样本，每个样本包含2 0 个特征值和1个属性值，其中属性值0、1、2 和3分别代表正常机、箱体异响机、右盖异响机和尖叫机。从12 0 个样本中随机抽取8 0 个样本构成训练集，剩余40 个样本构成测试集。2.3.2支持向量机的参数选择SVM由统计学习理论演变而来，它建立了一套基于统计理论，完善且标准的机器学习方法，改善了分类识别方案设计的随意性，利用结构风险最小化原则保证决策函数的精度，是一种广泛应用的模式识别方法15。理论研究表明，SVM进行建模分类运算时如果引入核函数可以提高向量800010000800010000机的计算速度和分类精度。由经验可知，径向基核函数（Radial Basis Function Kernel,RBF）分类精度高、泛化一般，为了提高RBF的泛化能力，需要采用合适的方法确定两个重要参数c和g，同时需要将样本数据进行归一化处理进一步提高分类精度。目前，比较流行的SVM核函数参数寻优方法有网格搜索法（Grid Search,GS）、遗传算法寻优和粒子群算法寻优17。其中，遗传算法和粒子群算法属于启发式算法，算法比较复杂17。遗传算法具有较强的全局搜索能力，但是局部搜索能力弱；粒子群算法收敛速度快，但是容易陷入局部最优解17。网格搜索法原理简单，首先确定c和g的搜索范围和搜索步长，假设参数有m个取值，g有n个取值，一共构成mxn 个(c,g)组合；然后运用交叉验证原理计算每一个(c，g)参数组合的分类精度；最后选择分类精度最高的一组参数作为最优参数组合17。为了搜索全局最优参数组合，需要将搜索范围设置的足够大，且尽可能减少搜索步长，但是遍历全部参数组合非常耗时，可行性差。实践经验表明，（c，g)参数组合只在较小区间内有较高的分类精度，而在其他区域分类精度很低。因此，本文采用一种改进的网格参数寻优第2 0 期法，大幅度缩短参数寻优时间，且不降低分类精度。首先采用较大的搜索步长确定最优参数组合的小区间，然后在该区间内采用较小的搜索步长进一步确定全局最优参数组合。本文采用台湾大学的林智仁教授团队开发的LIBSVM工具箱进行发动机异响的模式识别17。首先选择C-SVC分类器，采用3折交又验证训练模型。然后确定参数c和g的初始搜索范围均为2 ,21；搜索步长为8；训练精度为8 3.7 5%。接着，设置参数c和g的细搜范围均为2 3,2 3,步长设置为0.12 5，训练精度约为95%。采用改进的网格寻优法确定最优参数c为0.12 5，g为0.2 5。采用这组参数组合对测试集进行测试，测试结果如图7所示，测试精度为97.5%。13.02.52.0型米1.51.00.502.3.3发动机异响的模式识别为了实现发动机异响快速检测功能，需要将信号处理过程进行集成化处理。因此，将信号去噪、特征提取、模式识别和图形显示集成为图形用户界面（GraphicalUser Interface,GUI）模块，方便检测人员查看异响类型和特征图像。值得注意的是，本文的样本数据来源于人工检测，样本数量有限，模型的测试精度可能存在偏差。通过不断地测试，本文提出的基于GS-SVM发动机异响检测方法的准确率在90%以上，具有一定的工程实用价值。3结论特征提取是发动机异响识别的关键环节，本文采用的小波包变换和双谱估计均能提取到发动机信号特征，但两种方法相结合，更能体现信号特征的多样性，提高了模式识别的准确率。SVM是解决小样本模式识别问题的有效方法，通过归一化样本数据，选择合适的分类器和核函数，并采用改进的网格搜索法确优化参数和g，提高了分类器的训练精度。杨兴国，等：基于GS-SVM的发动机异响在线检测平台声学技术,2 0 11(6):513-517.2杨诚,李爽,冯焘,等.基于LMS 和 SDP 的发动机异响诊断方法研究J.汽车工程,2 0 14,35(11):1410-1414.3郑海波，李志远,陈心昭，等.基于时频分布的发动机异响特征分析及故障诊断研究J.内燃机学报,2 0 0 2,20(3):267-272.O实际类型米预测类型美美51015202530354045样本序列图7测试集的精度45将发动机异响识别由人工听诊改为软件自动识别，既解放了劳动力，也提高了异响识别准确率的稳定性，且准确率与人工听诊相当，达到90%以上。因此，本文提出的方法进行发动机异响识别，可行性好，准确率较高，具有一定的工程实践意义。参考文献1杨诚,魏国军,王中方.机油泵齿轮异响检测方法.4臧玉萍,张德江,王维正.基于小波变换技术的发动机异响故障诊断J.机械工程学报,2 0 0 9,45(6):2 39-2 45.5杨诚,黄金旺,杨振冬，等.摩托车发动机异响在线检测系统平台开发J.汽车工程学报,2 0 13,3(5):37 4-37 9.6贾继德.基于小波对称极坐标方法的发动机异响故障特征提取研究J.汽车工程,2 0 13,35(11):10 35-1037,1029.7李皓,陈长征,王洁.基于小波时频偏向干分析的车内异响声源识别机改善J.噪声与振动控制,2 0 2 2,42(1):173-177,186.8安周鹏,肖志怀,陈宇凡,等.小波包能量谱和功率谱分析在水电机组故障诊断中应用.水力发电学报,2 0 15,34(6):18 2-190.9王昱翔,殷希梅,陈凯，等.车辆发动机异响信号的故障特征参数提取方法J.兵器装备工程学报,2 0 2 0,41(2):44-48.10周全.汽车发动机异响诊断算法及产生机理研究D.杭州：浙江大学,2 0 2 1.11孙洁娣，靳世友.基于小波包能量谱及高阶谱的特征提取方法J.天津大学学报,2 0 10,43(6):56 2-56 6.12 程静,王维庆,何山.双谱分析法在风机轴承故障诊断中的应用J.自动化仪表,2 0 16,37(7):2 7-2 9.13范会敏,王浩.模式识别方法概述J.电子设计工程，2012,20(19):48-51.14李爽.基于小波滤波和极坐标镜像法的发动机异响在线检测研究与应用D.重庆：重庆大学,2 0 15.15杨兴国.发动机异响的特征提取与模式识别D.重庆：重庆大学,2 0 17.16 杨兴国.基于小波包和双谱的发动机异响模式识别J.噪声与振动控制,2 0 18,38(6):146-149,17 7.17王健峰,张磊,陈国兴,等,基于改进的网格搜索法SVM参数优化J.应用科技,2 0 12,39(3):2 8-31.