基于数学形态学的汽车轮毂单元振动信号消噪算法研究

汽车轮毂单元是将与轮毂轴承相配套的零件即轮毂或转向节和轴承圈套制成-体的结构。汽车轮毂轴承单元的主要作用是承受汽车的重量，并传输动力至轮胎。其受力包括径向载荷与轴向载荷，是非常重要的安全件。因此，实时监测轮毂单元的运转状况 ，对于提高行车安全、防止安全事故具有重要意义。

收稿日期：2012-t0-26基金项目：国家自然科学基金资助项 目(51105124)；浙江省自然科学基金资助项目(Y1090199)作者简介 ：孟庆华 (1977-)，男 ，山东泰安人 ，副教授 ，博士 ，硕士生导师，主要从事 电动汽车设计、机械故障诊断等方面的研究.E-mailmengqinghua###gmail.corn机 电 工 程 第30卷在运转过程中，包含局部缺陷的汽车轮毂单元会产生间断性的脉冲冲击激振，使所测振动信号包含快速衰减的脉冲响应分量。这种瞬时频率突变的持续时间极短，同时所测振动信号包含测量噪声、表面波纹度引起的振动等干扰，脉冲冲击特征往往被淹没。

因此，在强噪背景下有效提取信号中的脉冲冲击成分是轮毂单元故障诊断的关键。随着小波分析与奇异值分解在信号降噪中的应用n ，信号降噪技术取得了- 定的进展。但是小波降噪与奇异值分解分别由于阈值选韧奇异值选择的不确定性，导致降噪效果受到影响。

数学形态学 (Mathematical morphology)是基于积分几何和随机集论建立起来的非线性图像(信号)处理、分析工具，已在图像处理 、模式识别、形状分析等方面得到广泛应用。数学形态学还具有十分有效的非线性滤波功能，其信号处理的效果只撒于待处理信号的局部形状特征，而不需要信号的频域信息，并且其算法只有加、减法和全值计算，不涉及乘除法，运算速度快↑年来逐渐被引用至-维信号分析领域 ，在电力系统 。 、振动 、心电 以及语音信号n 处理中得到了逐步应用。

为了对汽车轮毂单元振动信号进行有效消噪，本研究提出基于数学形态学运算和软阈值的振动信号消噪方法。该算法只涉及加减和全值计算，运算简单且执行高效，适用于汽车轮毂单元故障的在线监测与诊断。

1 数学形态学基本原理数学形态滤波是基于信号的几何特征，通过结构元素探测目标信号的位置，以实现对信号的匹配，从而达到信号提娶细节保持及噪声抑制的目的。数学形态学基本运算包括腐蚀、膨胀运算，以及由此构造出的开、闭运算等。

设原始信号f(n)为定义在F(0，1，，Ⅳ-1)上的离散函数，定义结构元素b(n)为B(0，1，， -1)上的离散函数，且N>M，则f(n)关于b(n)的腐蚀0”和膨胀0”分别定义为：6)(n)minf(n -b(m) m0，l，2，M-1) L.。6) )maxf(n-m)十6( l ，m0，1，2，M-1) zf(n)关于6(n)的开运算。”和闭运算·”分别定义为：Or 。6)( orob0 6) ) (3)Or6)(几)OrbOb)(n) (4)开、闭运算以不同方式平滑信号，在信号滤波中应用普遍。其中，开运算抑制信号的正脉冲(波峰)，闭运算抑制信号的负脉冲(波谷)。

2 基于形态学的消噪算法研究形态学滤波的性能主要撒于结构元素b(n)的宽度 ，当信号中噪声脉冲的宽度不超过结构元素宽度时，噪声脉冲可被开、闭运算去除。本研究为有效抑制噪声正、负脉冲，可采用开- 闭(F。 )、闭-开(F 。)运算组合的形态学滤波运算(Morphology Filter，MF)，如下式所示：MF Or) [F ㈤)F 。 )] (5)其中：F。 (-厂㈣) (凡)。 )6 (6)F 。 ) ㈣。6(n)。6 (7)由于轮毂轴承单元振动信号与噪声的时频特性不同，若结构元素宽度合适，其形态学滤波运算结果中将保留信号大部分的特征波形，而脉冲冲击的峰谷信息和高频噪声则被过滤。由原始信号减去形态学滤波运算结果可得到过滤信号，即形态学峰谷提柔果，通过对其进-步处理就可以得到原始信号中被除去的峰谷信号。定义信号 关于6( )的形态学峰谷提取运算(peak valey extract，PVE)为：P -MF Or) -寺[F。 (n)Fco (凡))]， 。

厶. (8) 1 、u，，-吉l厂(n)。6 (n)b2(n)厂(n)6 (n)。b2(凡)l经处理后的形态学运算提取的峰谷信号与形态学滤波的结果相加，可以作为汽车轮毂单元振动信号消噪处理的最终结果。该结果在最大限度去除噪声的前提下，较完整地保留了有用信号▲-步分析该结果可实现轮毂轴承单元的故障检测与诊断。

综上分析 ，本研究提出基于形态学运算和软阈值的轮毂单元振动信号消噪算法，并在这-基础上对信号进行特征提取，基本流程如图1所示。

2.1 形态学滤波及峰谷提取结构元素在数学形态运算中的作用类似于-般信号处理中的过滤窗口，其形状和大小的选取对于信号分析至关重要♂构元素的结构主要有扁平型、直线型、三角形、圆形及其他多边形等，本研究针对信号受随机噪声和脉冲冲击影响的特点，采用菱形和椭圆形结构元素，菱形结构元素对脉冲噪声有较好的滤波效果，椭圆形结构元素用于对随机噪声信号的过滤，宽度取信号最大值的1/50～1/10 。

第4期 孟庆华，等：基于数学形态学的汽车轮毂单元振动信号消噪算法研究图1 基于数学形态学和软阈值处理的轮毂单元振动信号消噪算法流程图本研究将菱形结构元素bi(n)和椭圆形结构元素6：( )代人式(6，7)，可得到：y(n)MFbif)：去[ )F 。 )]， 二 (9) ' ，吉I (n)。6 (n)b2(n) (n)bl(n)。6 (n)l则峰谷提取信号结果为：∽ FVE if). ， 1n 、f-专 (n)。6。(n)·b2(n) (n)·bl(rt)。6 (凡)]u在形态学运算时，结构元素的宽度越大 ，所滤除的噪声干扰越多，但同时也损失了大量有用信号 ；反之，结构元素的宽度越小，保留的有用信号越多，但滤除的噪声干扰也越少。通过多次试验 ，本研究选取 b。

的最大宽度为5，b 的最大宽度为 12时，取得了较好的试验效果。

2.2 基于3r准则的软阈值处理本研究采用适当宽度的结构元素对汽车轮毂轴承单元振动信号进行形态学滤波运算后，脉冲冲击在形态学峰谷提柔果中表现为少量尖峰，而噪声的峰谷信息则为主导因素，且呈现高斯白噪声的分布规律。由统计学3仃准则n 可知，如果峰谷信号中噪声信号的均值为 ，方差为 ，那么其取值几乎全部(约占99.7%)集中在 -3er， 3tr内。另外 ，由于峰谷信号中脉冲冲击成分的瞬时幅值通常明显高于噪声幅值，可看作噪声信号中的突变点。本研究将3叮准则应用于峰谷信号阈值处理中，可在保留有用信号的同时，将故障特征(脉冲冲击)和噪声区分开。

该算法的具体步骤如下：(1)计算峰谷提取信号的均值和方差。首先，剔除有用信号干扰，对峰谷提取信号分别求均值 和方差or ，将信号中数值在 3or ， 3o'q外的点剔除，并用 -3o- 或 3o" 代替，得到新的峰谷提取信号；然后对新的峰谷提取信号再次求均值 和方差 or 。

(2)对峰谷信号z(n)进行阈值处理。本研究将峰谷提取信号中幅值在 -3o"，1.3卅内的点认为是噪声而舍弃，其他点认为是轮毂单元故障信号予以保留。具体计算如下式所示 ：，， 、 Jz ， 3tr t，q10， 其他 (11)1，2， ，式中： (n)-经阈值处理后的峰谷信号。

2-3 形态学滤波结果与阈值化峰谷信号求和本研究把形态学滤波结果y(n)和经阈值处理后的峰谷信号 ，(n)相加，得到最终消噪处理结果， (n)，如下式所示：厂 (n))， )z (凡) ( 1，2，，Ⅳ) (12)3 仿真与试验信号分析为检验所提出算法的有效性，笔者分别对仿真信号和轮毂单元振动试验信号进行消噪。本研究采用Matlab对所提出的算法进行了仿真实验验证。

3.1 仿真实验设采样频率 为10 kHz，采样长度 为2 000点，则采样时间为：t(L-1)lZ(s)正常振动信号为：0.7 sin(2'r.f，t)sin(2 )式中： ， -50 Hz、120 Hz；故障信号p(t)-幅值为3 mV的脉冲信号；n( )-白噪声信号。

则染噪故障信号为：y ) ( )p( ) ( )以上仿真信号时域波形如图2所示。

用传统形态学滤波方法、小波降噪方法(db3小波基、软阈值)及本研究方法的信号消噪结果以及频域· 414 · 机 电 工 程 第30卷-45童。

- 5t/s(a)正常信号 t/s(b)故障信号p420-2- 4tIs(c)白噪声n0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2t/s(d)染噪故障信号y(t)图2 仿真试验信号建模分析结果如图3、图4所示。其中，用传统形态滤波器滤波后的信号的时域波形如图3(a)所示；用bd3小波基分解重构的信号时域波形如图3(b)所示；用本研究方法处理得到的信号的时频波形如图3(c)所示，箭头标注位置为识别的故障信号位置，大多数可被明显观察到。含噪故障信号频域图如图4(a)所示；用db3小波基分解重构的信号频谱图如图4(b)所示；用本研究方法消噪处理后所得到信号的频域分析图如4(c)所示。无噪故障信号的时域波形及频谱分析结果如图5所示。

从仿真结果中可以看出：传统形态学滤波器在过滤大部分噪声的同时，也滤掉了脉冲信号，即故障信号无法识别；小波消噪方法在阈值选择上存在-定难度 ，本研究采用软阈值消噪，得到的结果中仅保留了5耄。

- 5。5- 5t/s(a)用传统形态滤波器滤波后的信号的时域波形0 0.02 0.04 0.06 0.08 O.1 O.12 O.14 O.16 O.18 02tIs(b)用bd3小波基3层分解重构的信号时域波形t/s(e)本研究方法消噪后信号的时域波形图3 受随机噪声干扰的故障信号消噪仿真试验f/Hz(a)含噪故障信号的单边频谱图f/I-Iz(b)用db3小波基分解重构信号的单边频谱图f/Hz(c)本研究方法消噪后信号单边频谱图图4 受随机噪声干扰的故障信号及消噪后信号的频谱分析结果tls(a)无噪故障信号时域波形j/Hz(b)无噪故障信号频谱分析结果图5 无噪故障信号的时域仿真和频谱分析结果部分故障信号，但是由于能量泄漏较大，故障点位置的幅值较小，不易识别；本研究提出的方法在时域可以明显观察到故障信号，但是由于结构元素的宽度没达到最优，在滤波过程中有部分噪声信号未能消除，反映在频域上的频谱分析结果不如小波消噪的好；观察图3(c)中频率800 Hz～1 000 Hz的范围，具有与图4(b)相同的脉冲信号的频谱特征，由于脉冲信号能量较小，特征不够明显。从总体分析，该算法只有加减法和全值计算，不涉及乘除法，运算速度快，且对故障的识别效果很明显，因此具有较大的优势。

3.2 试验信号分析试验对象为浙江某企业制造的轮毂单元，本研究对其在不同工况下的运行状况，分别采集了磨合期、磨合后及运转95 h后在不同工况下的数据信号。具]-九 - - 二V-O 订- 打～- V-。 0第4期 孟庆华，等：基于数学形态学的汽车轮毂单元振动信号消噪算法研究 .415 。

体工况为300 r/min、500 r/min、700 r/min等不同转速和0.2 g、0.4 g、0.6 g等不同负载下的9种组合。采样频率为10 240 Hz，采样时间各 10 S。轮毂单元试验台如图6所示。

> 50500图6 轮毂单元试验台0 0.1 0.2 0-3 0.4 0.5 0.6 0.7 0.8 0.9 1(a)原始信号童 8- 50- - - - L----J-----J----------1----- -----L-----L-----L0 0.1 0.2 0-3 0.4 0.5 0.6 0.7 0.8 0.9 1(b)用传统形态学滤波器消噪后信号O O.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1(c)用小波方法消噪后信号喜- 580 I- . 。 ：△ .!曼 .! 。 。

0 0.1 O.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 l(d)用本研究提出方法消噪后信号f/Hz(e)原始信号的单边频谱№ (f)用本研究提出方法消噪后信号的单边频谱图7 转速500 r/min、负载0.4 g时信号的时域波形和频谱分析结果通过对其中运转95 h后的两种工况：500 r/min、0.4 g负载及500 r/min、v.0 6 g负载下的采集信号，分别采用传统形态学滤波器消噪、小波消噪及本研究提出方法进行信号处理，得到结果如图7、图8所示。

其中，转速500 r/rain、负载0.4 g时信号的时域波形及频谱分析结果如图7所示。其中，原始信号如图7(a)所示，用传统形态学滤波器消噪后信号如图7(b)所示，用小波方法消噪后信号如图7(c)所示 ，用本研究提出方法消噪后信号如图7(d)所示；原始信号的单边频谱如图7(e)所示，用本研究提出方法消噪后信号的单边频谱如图7(f)所示 ；转速 500 r/min、负载0.6 g时信号的时域波形及频谱分析结果如图8所示，其中原始信号如图8Ca)所示，用传统形态学滤波器消噪后信号如图8(b)所示，用小波方法消噪后信号如图8(c)所示，用本研究提出方法消噪后信号如图8(d)所示；> 100县 0- 1oo0 0.1 0.2 0.3 0.4 0.5 0.6 O.7 0.8 0.9 1(a)原始信号圣- 100L----- --.-------J------L--- - ----J-----J------L-----L-0 0.1 0.2 0.3 0.4 0.5 O.6 0.7 0.8 O.9 lf，s(b)用传统形态学滤波器消噪后信号100量 0- 1000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1t/s(c)用小波方法消噪后信号巫0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1(d)用本研究提出方法消噪后信号f/rlz(e)原始信号的单边频谱f/Hz(f)用本研究提出方法消噪后信号的单边频谱图8 转速500 r/rain、负载0.6 g时信号的时域波形和频谱分析结果机 电 工 程 第30卷原始信号的单边频谱如图8(e)所示，用本研究提出方法消噪后信号的单边频谱如图8(f)所示。从频谱图中，得到如下信息：在转速500 r/min、负载0.4 g这-工况下 ，在频率范围50 Hz～8O Hz之间有-段能量集中的频率段 ，尤其在频率为60 Hz附近有-个明显的能量集中点，此外，在频率范围350 Hz～400 Hz之间也有- 段能量集中的频率段；在转速500 r/min、负载0.6 g的工况下，本研究同样观察到在频率范围50 Hz～80 Hz之间有-段能量集中的频率段 ，并且频率为60 Hz附近有-个明显的能量集中点，而其他频率上没有能量集中的能量带。在时域上，本研究提出的消噪算法在信号消噪结果中可以观察到-定时间间隔的脉冲信号，如图7(d)中A、B、c 3点位置，脉冲信号较其他滤波方法明显。

通过对时域波形和频域波形观察分析可知，本研究提出的消噪算法在保留大部分有用信号的同时，有效地抑制了噪声，且故障信号识别率提高了20%左右。该方法优于形态学滤波和小波消噪方法。

4 结束语针对汽车轮毂单元振动信号的有用信号往往被背景噪声淹没及故障信号不易识别的特点，本研究提出基于形态学运算和自适应阈值的汽车轮毂单元振动信号消噪算法。通过仿真和试验信号分析，验证了该方法的有效性。该方法与形态学滤波相比，消噪效果更好 ，相比小波消噪方法，本研究方法运算更简单、高效。

基于形态学运算和 自适应阈值的消噪算法有其突出的优点，但还存在如下不足：①该方法在应用过程中，结构元素的选择是个难点，尤其是结构元素的结构和宽度选择，本研究方法所采用的菱形和椭圆形结构元素的宽度值选取参考经验值，-般取信号最大幅值的1/50～1/10，关于结构元素宽度的确定方法还需进-步研究；②基于3 准则的自适应阈值提取过程中可能将实验数据中的坏点作为故障信号，因此 ，研究者应将待分析的信号预处理或作进-步分析，需要提出适当的算法对其处理。