Walsh变换与数学形态学的特征提取对比研究

径 哲郭利

（军械工程学院导弹工程系，河北石家庄050003）

摘要：为了从液压泵振动信号中准确提取冲击成分的频率特征，研究了Walsh变换和数学形态学在振动信号特征提取中的应用。Walsh变换具有能够分析微弱脉冲信号的特点，且在特征提取中具有较好的敏感性和抗噪能力。数学形态学利用其独特的结构元素探针，在信号冲击特征的提取方面效果较好。仿真和实验结果表明，数学形态学方法比Walsh变换更能够清晰地提取振动信号的频率特征。

关键词：Walsh变换傅里叶变换数学形态学功率谱液压泵特征提取

中图分类号：TH16;TP15 文献标志码：A DOI：10.  16086/j．cnki.  issnl000  - 0380. 201512006

Abstract: In order to clear out of the hydraulic pump vibration signals to extract the pulse characteristic frequency, the method of hydraulic pump vibration signals feature extraction is presented based on the Walsh transform and mathematics morphology. Walsh transform has advantages in analyzing the weak pulse signals, and shows the sensitivity and anti-noise ability to extract the feature. Mathematics morphology has advantage in pulse signals feature extraction with using its unique structural element probe. Through digital simulation test and hydraulic pump simulation experiment, the results show that mathematics morphology can more clearly than the Walsh transform to extract the frequency characteristics of vibration signals.

Keywords:Walsh transformFourier transform Mathematics morphologyPower spectrumHydraulic pumpFeature extraction

0引言

 液压泵是液压系统的“心脏”，其性能好坏直接决定着整个液压系统能否正常运行。液压泵的振动信号因包含丰富的状态信息，且易于获取，是分析液压泵故障诊断重要的特征信号之一。对于轴向柱塞式液压泵而言，其振动信号常含有较多的冲击成分，因此，对液压泵振动信号中冲击成分的频率特征进行提取，是实现液压泵故障诊断的基础。

 Walsh变换与数学形态学近年来常被用于信号的特征提取。Walsh变换与Fourier变换相似，也是信号的处理方法之一，与建立在复数乘法之上的Fourier变换相比，Walsh变换中只有简单的加减法，所以Walsh变换的计算速度极快，其相应的功率谱较Fourier功率谱而言，敏感性和抗噪能力更强。数学形态学通过结构元素探针在信号中的移动可以提取信号的冲击特征，和时域、频域的其他数学方法不同，其抑制脉冲干扰能力较强。

 针对液压泵振动信号常含有冲击成分的问题，本文将Walsh变换和数学形态学引入到信号特征提取中来，通过仿真和实验分析，数学形态学在提取信号冲击特征频率方面优于Walsh变换。

1 Walsh变换

 Walsh函数系是一组完备的正交函数系，由+1和-1两个元素构成。在实际工程应用中，常利用离散Walsh变换处理有限长序列离散信号。离散Hadamard变换和离散Walsh变换本质上相同，只是排列次序略有差异。Hadamard变换具有简单的递推关系，对于高阶矩阵的计算其更加简单直接，因此，常利用Hadamard变换进行高阶矩阵的计算。

离散Hadamard变换的代数式可定义为：

式中：{B(K)}为Hadamard变换后的系数序列；K=O，1,2，…，N-1;n=0,1,2，…，N-1.

 根据矩阵形式的定义，可得Hadamard变换的矩阵定义式：

式中：日。为Ⅳ阶Hadamard矩阵。

 常见的Walsh变换有三种排序方式，分别为Walsh序、Paley序和Hadamard序。三种序列之间可以互相转换，关系如下：

  在计算功率谱时，常选用Walsh序的变换系数：C(O)，C(I)，…，C(N-1)，则Walsh变换的功率谱为：

2数学形态学

2.1数学形态学基本运算

 数学形态学的基本运算包含腐蚀和膨胀两种算子。设一维多值信号x(n)的定义域为x={0，1，2，…，N一1}，一维结构元素序列y(n)的定义域为Y={O，1，2，…，M-1}，其中Ⅳ和M都是整数，且Ⅳ>M，则x(n)关于y(n)的膨胀和腐蚀分别定义为：

 数学形态学中的腐蚀运算有抑制正向冲击和平滑负向冲击的作用，而膨胀运算则相反，在平滑正向冲击的同时抑制了负向冲击。

 利用腐蚀和膨胀的概念，可将x(n)关于y(n)的开、闭运算分别定义为：

 数学形态学中的开运算常用来去除信号边缘的毛刺，同时可以滤除信号上方的峰值噪声；闭运算常用于平滑或抑制信号下方的波谷噪声，起到填补信号漏洞和裂纹的作用。

2.2形态学滤波器

 利用形态学的开、闭运算可以构建两种常用的滤波器：平均滤波器( AVG)和差值滤波器(DIF)。二者分别表示如下：

 AVG可以在消除信号正、负冲击的同时平滑信号，而DIF可以提取信号中的正、负冲击。

2.3形态学结构元素

 结构元素形状和大小的选择对滤波效果有重要的影响，是形态学的基本算子之一。常用的形状有三角、半圆、菱形等。结构元素越复杂，计算量越大，提取信号的效果也越好。考虑到提取效果和计算量等综合因素，处理信号时常选用扁平的圆盘形作为结构元素。

3仿真信号分析

 在对比研究中，采用如下形式的仿真信号模拟液压泵出现故障时的振动信号：

 在信号s(t)中，冲击分量包含两部分：xi(t)模拟液压泵振动信号中的固有冲击成分，X2(f)模拟液压泵出现故障时故障引起的冲击成分。n(t)代表白噪声成分，加噪后信号的信噪比为-5dB。仿真信号的采样频率为2 048Hz，采样点数4 096。图1为s(t)的时域图和频谱图。

 由仿真信号易得，信号冲击成分的能量主要集中在270 Hz。现分别利用Fourier功率谱、Walsh功率谱和形态学差值滤波器对信号进行特征提取，结果如图2所示。由仿真图可以看出，Fourier功率谱在提取10 Hz和30 Hz的冲击频率特征的同时，能够提取10 Hz的2倍频，但是15 Hz和26 Hz混淆了频率特征的提取。Walsh功率谱提取的9.75 Hz与29.75 Hz与冲击成分频率十分接近，但是其提取出46. 75 Hz与87 Hz亦对频率特征的提取造成一定的影响。由形态学差值滤波器对频率特征进行提取，能准确地提取10 Hz的l、2、3、4倍频和30 Hz的1、2倍频，10 Hz的3倍频和30 Hz相重合，所以30 Hz处的频谱幅值较大。综上所述，形态学差值滤波器的特征提取效果最好。

4液压泵实测信号分析

 在液压泵的常见故障模式中，斜盘磨损和滑靴松动是最容易发生的两种故障。滑靴与斜盘之间的摩擦副是斜盘式轴向柱塞泵多对摩擦副中最复杂的一对，对液压泵的性能影响较大。在液压泵工作过程中，出现斜盘磨损故障或出现松靴故障时，都会影响其正常工作，降低液压泵寿命，甚至发生重大事故。

本文采用的试验器件名称及型号如表1所示。

 其中，驱动电机的额定转速为1 480 r/min；液压泵的柱塞数为7，理论排量为10 ml/r，额定转速为1 500 r/min，液压泵主溢流阀压力为10 MPa，采样频率为20 kHz，采样点数为10 000。由于泵轴的转速为1 500 r/min，单个柱塞附加冲击的基频为f=n/60（n为电机的转速），则本试验的液压泵的冲击振动基频为I 480/60 x7 =172. 667 Hz。

  采集到的松靴和斜盘磨损的信号的时频图分别如图3、图4所示。

 从采集信号的频谱图中只能看出液压泵振动信号的能量主要集中在几个频率上，并不能提取信号基频等频率特征。下面分别对振动信号求取Fourier功率谱、Walsh功率谱以及通过形态学差值滤波器提取信号的特征频率，结果如图5所示。

 图5中，A代表振幅。由图5对比分析可得，Fourier功率谱提取出的170.9 Hz、354 Hz和517.6 Hz与理论分析的基频172. 667 Hz、2倍频345. 334 Hz和3倍频518. 001 Hz较为接近，但频谱幅值较低，效果一般。Walsh功率谱提取出的接近基频的170.9 Hz幅值较高，频率特征较明显，但358.9 Hz与859.4 Hz与2倍频345. 334 Hz、5倍频863. 335 Hz相差较多，且提取出的频率有干扰成分。由形态学差值滤波器提取的172 Hz、344 Hz、516 Hz与688 Hz与理论分析的基频及其2、3、4倍频较为接近，效果最好。松靴信号特征提取曲线如图6所示。

 由图6对比分析可得，Fourier功率谱提取的344.2 Hz与振动信号基频的2倍频较为接近；而Walsh功率谱提取的168.5 Hz与178.2 Hz在基频附近，但是基频已被湮没，48.83 Hz等频率对特征频率的分析噪声一定的影响。形态学差值滤波器提取的172 Hz、344 Hz与振动信号的基频及其2倍频较为接近，虽然50 Hz等频率也有一定的影响，但基频幅值较高，可以提取信号的特征频率。

5结束语

 传统的Fourier功率谱与Walsh功率谱、形态学差值滤波器均能提取液压泵振动信号的特征频率，为故障预测和诊断打下良好基础。通过仿真试验和液压泵模拟故障实验，结果表明，形态学差值滤波器在提取频率特征方面优于Fourier功率谱和Walsh功率谱。