最大相关峭度解卷积结合稀疏编码收缩的齿轮微弱.DOC

1、最大相关峭度解卷积结合稀疏编码收缩的齿轮微弱故障特征提取唐贵基，王晓龙（华北电力大学能源动力与机械工程学院，河北 保定 071003）摘要：针对强背景噪声环境下齿轮早期故障诊断问题，提出了最大相关峭度解卷积结合稀疏编码收缩的微弱故障特征提取方法。由于最大相关峭度解卷积算法的处理结果同时受滤波器长度参数及解卷积周期参数的影响，为自适应地实现最佳的解卷积效果，利用粒子群算法优良的寻优特性，对最大相关峭度解卷积算法的最佳影响参数组合进行搜索。原故障信号经影响参数优化的最大相关峭度解卷积算法处理后，冲击特征会明显增强，为剔除剩余噪声，对所获解卷积信号做进一步稀疏编码收缩降噪处理，并通过分析降噪信号的包

2、络谱来识别故障特征频率成分。实例分析结果验证了该方法的有效性和可靠性。关键词：齿轮；微弱特征；粒子群优化；最大相关峭度解卷积；稀疏编码收缩中图分类号：TH132.417；TP206 文献标识码：AWeak feature extraction of gear fault based on maximum correlated kurtosisdeconvolution and sparse code shrinkageTANG Gui-ji， WANG Xiao-long(School of Energy，Power and Mechanical Engineering，North China

3、 Electric Power University，Baoding 071003，China)Abstract： Aiming at incipient fault diagnosis problem of gear under strong background noise, an feature extraction method for weak fault based on maximum correlated kurtosis deconvolution and sparse code shrinkage was proposed. As the processing result

4、 of maximum correlated kurtosis deconvolution algorithm was affected by filter length parameter and deconvolution period parameter, in order to achieve the best deconvolution result adaptively, particle swarm optimization algorithm with excellent optimization characteristic was used to search for th

5、e optimal combination of influencing parameters of maximum correlated kurtosis deconvolution algorithm. The impact characteristic of original fault signal could be enhanced after processed by maximum correlated kurtosis deconvolution algorithm with optimized parameters, in order to eliminate the res

6、idual noise, the deconvolution signal was further processed by sparse code shrinkage de-noising algorithm, then fault characteristic frequency components could be identified by analyzing the envelope spectrum of de-noising signal. The analysis results verified the effectiveness and reliability of th

7、is method.Key words： gear； weak feature； particle swarm optimization； maximum correlated kurtosis deconvolution； sparse code shrinkage引 言齿轮是机械设备中广泛应用的传动部件,但故障发生概率高，容易损坏，准确有效的提取齿轮早期故障微弱特征对于避免发生严重故障、 减少经济损失意义重大。在实际工程应用中，由于振动传输路径复杂多变、工作环境噪声干扰严重、多振动源的激励和响应互相耦合等诸多因素的影响，使得齿轮早期微弱故障特征提取相对困难, 这也一直是故障诊断领域的一个研

8、究热点 1-2。针对该问题，不少学者进行了深入研究，文献3利用经验模态分解方法对精轧机齿轮故障信号进行自适应分解，并对得到的固有模态分量做Hilbert 包络解调来获取每个分量的幅频信息，最终实现强噪声背景下齿轮早期裂纹故障的准确诊断；文献4利用小波包的高分辨率分解及重构能力，将齿轮微弱故障信号分解到不同频段，然后选择有效频段进行信号重构，顺利提取出微弱故障特征信息；文献5将共振稀疏分解与包络谱相结合，采用共振稀疏分解将微弱冲击从齿轮振动信号中分离出来，然后对冲击序列做包络分析，可对齿轮早期故障做出准确判断；文献6提出离散小波变换与神经网络相结合的锥齿轮早期故障智能识别方法，首先通过离散小波变

9、换提取故障特征，继而利用决策树来筛选故障特征，最后通过神经网络实现故障分类；针对齿轮微弱故障特征提取困难的问题,文献7提出基于稀疏信号分解的多阶分数阶傅里叶变换自适应滤波方法，该方法能很好地从频率呈曲线变化的多分量信号中剥离出啮合频率包络调制信号，对滤波后的信号做进一步解调分析，可顺利提取出微弱故障特征，上述各方法在齿轮早期故障诊断应用中均取得了一定的效果。最大相关峭度解卷积(Maximum correlated kurtosis deconvolution, MCKD)8可通过循环迭代2实现解卷积运算、突出信号中包含的连续脉冲，但 该算法的处理结果受滤波器长度及解卷积周期两个参数的影响，如果

10、任何一个参数设置不合适，都难以达到理想的分析效果。为了实现参数的自适应筛选，本文提出利用粒子群算法对 MCKD 算法的最佳影响参数组合进行搜寻，并将影响参数优化的 MCKD 算法与稀疏编码收缩算法相结合，用于处理齿轮早期故障信号，提出了最大相关峭度解卷积结合稀疏编码收缩的微弱故障特征提取方法，实例分析结果表明该方法可实现齿轮早期故障的精确诊断。1 基于粒子群优化的 MCKD 算法1.1 MCKD 算法零均值信号 ( )的相关峭度表ny1,2N达式为：(1)213()nTnNyCKT（ ）相关峭度是在峭度基础上提出的概念，当解卷积周期参数 为 0 时，相关峭度即退化为峭度，然而峭度对少数异常突出

11、脉冲过于敏感，相关峭度则充分考虑了冲击成分的连续性，与峭度相比能更准确的衡量信号中周期性脉冲序列所占的比重。假设 为一个连续冲击信号， 为冲击信号yx经周围环境及路径传输后实际采集到的响应信号，上述过程可通过如下公式表示为： (2)xhye为便于分析，在此先不考虑噪声 的影响，则MCKD 算法本质是寻找一个 FIR 滤波器，通过输出信号 恢复输入信号 ，即：xy(3)11Lknfxf其中， 是长度为 的12=TL L滤波器系数。为突出解卷积运算结果中的连续尖脉冲，该算法以信号的相关峭度为评定标准，并将相关峭度最大化作为目标函数，即：(4)21()maxaxNnTyCKTff（ ）上述优化求解问

12、题等价于求解方程：（ ） (5)0kdCKTf（ ） 1,2kL由公式(3)- (5)可归纳得到如下表达式：(6)621142211()+2(+)NNnTnkkNnkkTyyxdfx（ ）以矩阵的形式将上式重新表述为：(7)246010()()4TdCKTyXyXf 其中，； 22201 TTNyy；2；13211210rrrNrrrNLrxxXx 。12 TTTyy整理式(7)可得：(8)2001()TXX由于已知：(9)0Tyf则滤波器系数可通过下式得到：(10)21001()()TTfXX将滤波器系数带入式(3)，即可得到实际采集信号 的解卷积信号 （原冲击信号） 。xy1.2 最佳影响

13、参数组合的自适应选取MCKD 算法可 有 效 增 强 信 号 中 隐 含 的 冲 击 脉冲 ， 适 合 处 理 信噪比低并具有周期性冲击特性的信号， 然 而 滤波器长度参数 和解卷积周期参数L的选择对解卷积结果有着至关重要的影响。如T果保持一个影响参数不变，仅以另一个影响参数为优化对象，讨论该参数对运算结果的影响，这3种局部寻优方式由于忽略了两个参数间的交互作用，因此得到的也只是相对最优影响参数。粒子群算法9作为一种群体智能优化算法，具有良好的全局寻优能力，本文采用粒子群算法对 MCKD算 法 的 两 个 影 响 参 数 进 行 并 行 搜 索 优 化 ， 可 实 现滤 波 器 长 度 参 数

14、 和 解 卷 积 周 期 参 数 的 自 适 应LT筛 选 。假 设 在 一 个 维 空 间 中 ， 由 个 粒 子 组 成 种DM群 ， 其 中 第 个 粒 子 表 示 一 个),(21MXXi维 向 量 ， 代 表 第 个 粒 子 在D,(iiix i维 搜 索 空 间 中 的 位 置 （ 即 优 化 问 题 的 潜 在 解 ） 。第 个 粒 子 的 速 度 为 ， 其 个 体 局i ),v,(ViDii21部 极 值 为 ， 种 群 全 局 极 值 为),(21iipP， 每 个 粒 子 通 过 个 体 局 部 极 值 和),(21DgG种 群 全 局 极 值 迭 代 更 新 自 身 的

15、 速 度 和 位 置 ， 公 式表 述 如 下 10：(11)112)()(kidikid kidkidivxxgcxpcv其 中 ， 为 惯 性 权 重 ； ；D,1； 为 当 前 迭 代 次 数 ； 和 为 加 速Mi,21 kc2度 因 子 ； 为 介 于 间 的 随 机 数 。 1,0利 用 粒 子 群 算 法 搜 寻 影响参数时，需确定一个适应度函数，粒子每次更 新 位 置 时 计 算 一 次 该位 置 对 应 的 适 应 度 值 ， 通 过 适 应 度 值 对 比 进 行 更新 。 由 于 包 络 谱 熵 能 有 效 衡 量 信 号 中 的 周 期 性 冲击 成 分 11， 本 文

16、 便 将 解 卷 积 信 号 的 包 络 谱 熵 作 为寻 优 过 程 中 的 适 应 度 函 数 ， 用 于 评定 MCKD 算法 解 卷 积 运 算 结 果 的 优 劣 。 零 均 值 信 号的 包 络 谱 熵 可由 公 式 表 示 成 如)21)(,NjyEp下 形 式 ：(12)Njjjja()()p110log其 中 ， 是 的 归 一 化 形 式 ， 是 信 号jja)(j包 络 谱 的 离 散 化 形 式 。)(jy以 任 意 粒 子 的 位 置 为 MCKD 算法的滤 波iX器 长 度 参 数 和 解 卷 积 周 期 参 数 ， 设 定 好 这 两 个 影响 参 数 后 ， 利

17、 用 MCKD 算 法 对 原 信 号 进 行 处 理并 计 算 所 获 解 卷 积 信 号 的 包 络 谱 熵 ， 熵 值即为该粒子位置的适应度值。当 解 卷 积 结)Ep(Xi果 中 周 期 性 脉 冲 成 分 突 出 时 ， 包 络 谱 熵 值 较 小 ，而 解 卷 积 效 果 欠 佳 、 规 律 性 冲 击 成 分 不 明 显 时 ，包 络 谱 熵 值 则 相 对 较 大 ， 因 此 包 络 谱 熵 值 最 小 化即 为 最 终 搜 索 目 标 。 利 用 粒 子 群 算 法 搜 寻 优 化MCKD 算法 最 佳 影 响 参 数 组 合 的 具 体 步 骤 如 下 ：(1)初 始 化

18、粒 子 群 算 法 的 各 项 参 数 。(2)初 始 化 粒 子 种 群 ， 以 影 响 参 数 组 合 , LT作 为 粒 子 的 位 置 ， 每 个 位 置 对 应 一 对 影 响 参 数 ，随 机 产 生 一 定 数 量 的 参 数 组 合 作 为 粒 子 的 初 始 位置 ， 随 机 初 始 化 每 个 粒 子 的 移 动 速 度 。(3)以 粒 子 的 位 置 为 MCKD 算 法 的 两 个 影 响参 数 ， 并 对 原 信 号 做 解 卷 积 运 算 ， 计 算 每 个 粒 子的 位 置 对 应 的 适 应 度 值 。iX)Ep(Xi(4)对 比 适 应 度 值 大 小 并 更

19、 新 个 体 局 部 极 值 和种 群 全 局 极 值 。(5)利 用 公 式 (11)更 新 粒 子 的 速 度 和 位 置 。 (6)循 环 迭 代 , 转 至 步 骤 (3)， 直 至 迭 代 次 数 达到 最 大 设 定 值 ， 输 出 最 佳 适 应 度 值 及 粒 子 的 位 置 。2 稀疏编码收缩稀 疏 编 码 收 缩 算 法 可 从 高 斯 信 号 中 推 测 得 到隐 含 的 非 高 斯 成 分 ， Hyvarinen 提出如下模 型 来表 示 非 高 斯 信 号 的 概 率 密 度 函 数 12：(13)/2)1(2)( )3(12dxdx）其 中 ， 是 信 号 的 标

20、准 差 ， 是 控 制 概 率 密度 函 数 稀 疏 性 的 参 数 ， 对 于 具 有 冲 击 特 性 的 非 高斯 信 号 来 说 ， 参数 可 设 定 为 0.113。假 设 为 原 始 信 号 ， 其 统 计 特 征 呈 现 出 非 高x斯 性 质 ， 是 均 值 为 零 的 高 斯 白 噪 声 信 号 ， 则 观v测 信 号 可 表 示 为 ：(14)vxy根 据 公 式 (13)中 的 先 验 分 布 模 型 ， 利 用 最 大似 然 估 计 原 理 可 从 观 测 信 号 中 推 测 得 到 原 始 信号 的 估 计 信 号 ：xx(15) )3(4)(21,0ma)( 2dyd

21、ysign其 中 ， ； 是 白 噪 声 信 号 的/)v标 准 差 ， 可 由 公 式 估 计 得 到 ，）（ yMAD代 表 观 测 信 号 的 绝 对 偏 差 的 中 值 ； 原）（ yMAD始 信 号 的 标 准 差 可 通 过 公 式 得 到 ，xd2yd为 观 测 信 号 的 标 准 差 。y43 诊断流程分 析 实 际 信 号 时 ， 选 择 合适的影响参数、有效突出信号中隐含的规律性故障冲击，是利用MCKD 算法处理齿轮早期故障信号的关键所在。本文利用粒子群算法自适应地筛选最佳的影响参数组合，可避免参数设定时人为主观盲目选择带来的弊端。齿轮故障信号经影响参数优化的 MCKD算法

22、处理后，信噪比会显著提高，原信号中被噪声所淹没的周期性脉冲能够有效突显出来，但仍然会存在较多噪声干扰。由于齿轮故障冲击成分通常具有非高斯性质，而噪声成分则呈现出高斯分布特性，稀疏编码收缩算法可实现高斯信号和非高斯信号的有效分离，因此利用该算法对所获的解卷积信号做进一步处理，保留冲击成分的同时可有效剔除剩余噪声。综合上述分析，本文提出最大相关峭度解卷积结合稀疏编码收缩的齿轮微弱故障特征提取方法，有望实现齿轮早期故障的准确识别，诊断流程如图 1 所示，具体实现过程如下：载入原始信号初始化各项参数粒子群算法并行寻优MCKD 算法处理信号稀疏编码收缩降噪Hilbert 包络解调运算包络信号频谱分析齿轮

23、故障识别图 1 诊 断 流 程 图Fig.1 Diagnosis flow chart(1)设定粒子群算法各项参数并搜寻 MCKD 算法的两个影响参数。在粒子群参数设置过程中，如果最大进化代数及种群规模这两个参数值设置的太大，则计算负担过重，如果设置的过小，则无法保证搜寻到最佳的影响参数组合，参考文献10中的参数设定值以及多次试验对比分析，本文设定的粒子群算法各项参数值如表 1 所示。表 1 粒子群算法各项参数Tab.1 Each parameter of particle swarm algorithm 最大进化代数 xG种群规模 M加速度因子 1加速度因子 2惯性权重 20 30 2 2

24、1(2)粒子群寻优结束后会得到一组最佳影响参数组合 , ，设定 MCKD 算法的滤 波 器 长 度oLT参 数 为 ，解 卷 积 周 期 参 数 为 ，并利用影响参oT数优化的 MCKD 算法对故障信号做预处理。(3)利用稀疏编码收缩算法对所得的解卷积信号做进一步降噪处理。(4)对降噪信号做 Hilbert 包络 解调运算，并计算包络谱。(5)将齿轮故障特征频率的理论计算值与包络谱中的谱线进行对比，实现故障最终判别。4 应用实例4.1 点蚀故障在齿轮箱故障试验台从动齿轮的齿面上加工出一个微小凹痕来模拟齿轮轻微点蚀故障，图 2为试验平台的结构简图，在齿轮箱端盖靠近轴承处布置四只加速度传感器，设置

25、的采样频率为5120Hz，采样点数为 8192 点。其中，电机转速为850r/min,输入轴主动轮齿数为 55，输出轴从动轮齿数为 75，通过计算可知输入、输出轴转频分别为 fi=14.17Hz、fo=10.39Hz 。图 2 试验平台Fig.2 Experimental platform实测信号波形如图 3 所示，时域波形中隐约可见少量冲击,但毫无规律可寻，对实测信号直接做包络谱分析，结果如图 4 所示，没有发现幅值突出的频率成分，由此可见传统包络谱分析手段对于该故障信号来说显得无能为力。00.20.40.60.811.21.41.6-5-2502550时 间 t/s幅值A/(m/s2)5图

26、 3 实测信号波形Fig.3 Waveform of measured signal010203040506070809010123频 率 f/Hz幅值A/(m/s2)图 4 实测信号包络谱Fig.4 Envelope spectrum of measured signal下面利用本文提出方法对实测信号进行分析，结果如图 5 所示。图 5(a)为包络谱熵随进化代数变化的关系曲线，包络谱熵最小值 7.1495 出现在第9 代进化种群中，粒子群算法寻优所得最佳影响参数组合为56,486，由此设定 MCKD 算法的滤波器长度参数为 56，解卷积周期参数为 486，利用影响参数优化的 MCKD 算法对

27、实测信号做预处理，结果如图 5(b)所示。通过分析解卷积信号的波形可发现，原本被强烈噪声所掩盖的冲击成分被有效突显出来。图 5(c)是采用稀疏编码收缩算法对解卷积信号做进一步处理所获的降噪信号，解卷积信号经过降噪后，大量冗余干扰噪声被剔除，而与齿轮故障相关的冲击成分则被保留下来，图 5(d)是降噪信号的包络谱分析结果，图中主要成分为输出轴转频及其倍频 ，由此判断输出轴上的off9从动齿轮存在局部损伤，分析结果与实际情况相符，实例分析表明本文提出方法可顺利识别出齿轮早期故障的微弱特征频率信息。12345678910112131415161718192077.27.47.67.88进 化 代 数包

28、络谱熵值(a) 包络谱熵随进化代数变化的关系曲线(a) Relation curve of envelope spectrum entropy changing with evolutional generation00.20.40.60.811.21.41.6-6-40-200204060时 间 t/s幅值A/(m/s2)(b) 解卷积信号波形(b) Waveform of deconvolution signal00.20.40.60.811.21.41.6-6-40-200204060时 间 t/s幅值A/(m/s2)(c) 降噪信号波形(c) Waveform of denoised

29、signal010203040506070809010123频 率 f/Hz幅值A/(ms2) 9fo8fo7fo6fo5fo3fo2fofo 4fo(d) 降噪信号包络谱(d) Envelope spectrum of denoised signal图 5 本 文 提 出 方 法 的 实 测 信 号 分 析 结 果Fig.5 Analysis results of measured signal by proposed method下面分别单独应用参数优化 MCKD 算法和稀疏编码收缩算法对实测信号进行处理来验证二者结合的必要性。图 6 为原信号经 MCKD 算法处理所得解卷积信号的包络谱（

30、MCKD 算法的滤波器长度及解卷积周期参数与上述分析相同，分别为56、486） ，虽然谱图中存在多个与故障特征相关的频率成分，但由于存在较多干扰谱线，特征频率基频及 5 倍频成分处谱峰并不突出。图 7 为原信号经稀疏编码收缩算法处理所得降噪信号的包络谱，图中仅特征频率基频处谱线幅值偏高，但峰值仍不够明显。通过对比可发现，将 MCKD 算法和稀疏编码收缩算法结合起来分析故障信号，可有效剔除原信号中的强烈噪声干扰成分，最大限度突出微弱故障信息。0102030405060708090100.511.522.53频 率 f/Hz幅值A/(m/s2) 2fo3fo4fo 6fo7fo8fo9fo图 6

31、解 卷 积 信 号 包 络 谱Fig.6 Envelope spectrum of deconvolution signal0102030405060708090100.511.522.53频 率 f/Hz幅值A/(m/s2) fo图 7 降 噪 信 号 包 络 谱6Fig.7 Envelope spectrum of denoising signal为验证所述方法优势，利用文献3所用的经验模态分解(Empirical mode decomposition, EMD)方法和文献4 所用的小波包分解方法分别对实测信号进行处理，并做包络解调分析，各取其最佳包络谱（所有分量均做包络谱，取效果最佳的分

32、量）与本文方法进行对比。图 8 为实测信号 EMD处理结果，自适应分解后共获得 7 个 IMF 分量，其中 IMF1 分量的包络解调效果最佳，图 9 为该分量的时域波形及包络谱。对实测信号做 2 层小波包分解，小波函数为文献4采用的 db5 小波，结果如图 10 所示，d20d23 是小波包分解所获 4 个节点系数的重构信号，其中 d21 节点重构信号的包络解调效果最佳，重构信号波形及包络谱如图11 所示。对比分析后发现，图 9(b)、图 11(b)所示包络谱中，齿轮故障相关频率成分的幅值水平以及倍频成分的连续性与图 5(d)相比，均存在较大差距，由此表明本文方法在齿轮微弱故障特征提取上具有一

33、定优势。 IMF1IMF2IMF3IMF4IMF5IMF600.20.40.60.811.21.41.6时 间 t/sIMF7图 8 EMD 处理的结果Fig.8 Processing results of EMD00.20.40.60.811.21.41.6-4-30-20-10010203040时 间 t/s幅值A/(m/s2)(a) IMF1 分量波形(a) Waveform of IMF1010203040506070809010123频 率 f/Hz幅值A/(m/s2) fo 3fo4fo5fo 8fo7fo(b) IMF1 分量包络谱(b) Envelope spectrum of

34、 IMF1图 9 IMF1 分量的波形及包络谱Fig.9 Waveform and envelope spectrum of IMF1d20d21d200.20.40.60.811.21.41.6时 间 t/sd23图 10 小波包处理的结果Fig.10 Processing results of wavelet packet00.20.40.60.811.21.41.6-3-20-100102030时 间 t/s幅值A/(ms2)(a) d21 节点重构信号波形(a) Waveform of reconstructed signal of d21 node010203040506070809

35、010123频 率 f/Hz幅值A/(m/s2) fo 3fo 7fo(b) 重构信号包络谱(b) Envelope spectrum of reconstructed signal图 11 重构信号的波形及包络谱Fig.11 Waveform and envelope spectrum of reconstructed signal4.2 磨损故障试验数据来自美国学者 J.Rafiee 的个人网站 14，试验装置包括齿轮箱、电机、加载机构、脉动分析仪、三轴加速计、转速计及减震器。电机实测转速为 1397r/min，输入轴主动轮齿数为 29，输出轴从动轮齿数为 24，输入、输出轴转频分别为fi

36、=23.29Hz、fo=28.14Hz。加速度计位于齿轮箱外部靠近输入轴位置处，采样频率为 16384Hz。试验过程采集了主动齿轮发生轻微磨损、中度磨损和断齿三种不同程度故障的振动信号，本文选择对轻微磨损阶段的实测信号进行分析，分析点数为16384 点。图 12 所示实测信号波形中，基本没有出现明显冲击脉冲,对该信号做进一步包络谱分析，结果如图 13 所示。图中存在多个谱峰，除了在输入轴二倍转频 处找到一根幅值偏高谱线外，再没有if2发现其它故障相关频率成分，仅通过这一有限信息很难准确判断齿轮箱是否出现故障。700.10.20.30.40.50.60.70.80.91-15-10-500501

37、0150时 间 t/s幅值A/(ms2)图 12 实测信号波形Fig.12 Waveform of measured signal0102030405060708090101012013014015012345频 率 f/Hz幅值A/(m/s2) 2fi图 13 实测信号包络谱Fig.13 Envelope spectrum of measured signal图14是利用本文方法对实测信号进行分析所得结果。图14(a)中，粒子种群仅通过3代进化就搜寻到了滤波器长度参数和解卷积周期参数的最佳组合212,91。通过观察图 14(b)解卷积信号波形可发现，经影响参数优化的MCKD 算法处理后，实测

38、信号中隐含的规律性冲击被顺利剥离出来，对解卷积信号做进一步稀疏编码收缩降噪处理后，噪声成分被剔除殆尽，脉冲序列的周期特征变的更加明显。图14(d)所示降噪信号包络谱中，输入轴转频及其谐波 处谱线峰值十分突出，基ifi6本不存在其它干扰成分，谱图相对干净。通过综合分析降噪信号的波形及包络谱可知，冲击脉冲出现的频率与输入轴转频相同，由此断定输入轴上的主动齿轮发生了故障。1234567891011213141516171819207.27.47.67.88进 化 代 数包络谱熵值(a) 包络谱熵随进化代数变化的关系曲线(a) Relation curve of envelope spectrum e

39、ntropy changing with evolutional generation00.10.20.30.40.50.60.70.80.91-15-10-5005010150时 间 t/s幅值A/(ms2)(b) 解卷积信号波形(b) Waveform of deconvolution signal00.10.20.30.40.50.60.70.80.91-15-10-5005010150时 间 t/s幅值A/(m/s2)(c) 降噪信号波形(c) Waveform of denoised signal0102030405060708090101012013014015051015频 率

40、f/Hz幅值A/(ms2) fi 2fi 3fi 4fi 6fi5fi(d) 降噪信号包络谱(d) Envelope spectrum of denoised signal图 14 本 文 提 出 方 法 的 实 测 信 号 分 析 结 果Fig.14 Analysis results of measured signal by proposed method同样，分别利用经验模态分解和 db5 小波包分解对实测信号进行处理（篇幅所限，所获各分量时域波形省略） ，并取最佳分量的包络谱与本文方法进行对比。经验模态分解处理后共获得 6 个IMF 分量，其中 IMF3 分量的包络谱效果最佳，图15

41、是该分量的波形及包络谱。2 层小波包分解重构后，发现 d23 节点重构信号的包络分析结果最佳，图 16 是重构信号的波形及包络谱。通过对比可知，图 15(b)和图 16(b)均没有图 14(d)效果明显，齿轮轻微磨损故障信号对比分析结果再次验证了本文方法的优势。00.10.20.30.40.50.60.70.80.91-2-1001020时 间 t/s幅值A/(m/s2)(a) IMF3 分量波形(a) Waveform of IMF3010203040506070809010101201301401500.511.522.53频 率 f/Hz幅值A/(m/s2) 2fi(b) IMF3 分量

42、包络谱(b) Envelope spectrum of IMF3图 15 IMF3 分量的波形及包络谱Fig.15 Waveform and envelope spectrum of IMF3800.10.20.30.40.50.60.70.80.91-8-60-40-20020406080时 间 t/s幅值A/(m/s2)(a) d23 节点重构信号波形(a) Waveform of reconstructed signal of d23 node0102030405060708090101012013014015012345频 率 f/Hz幅值A/(m/s2) fi 2fi 3fi(b)

43、重构信号包络谱(b) Envelope spectrum of reconstructed signal图 16 重构信号的波形及包络谱Fig.16 Waveform and envelope spectrum of reconstructed signal5 结 论(1)齿轮早期故障特征信号微弱，并受环境噪声及传递路径衰减影响，特征提取相对困难，传统包络谱方法可能诊断失效，利用 MCKD 算法分析早期故障信号，可有效增强微弱冲击特征。(2)滤波器长度参数和解卷积周期参数的选择对 MCKD 算法处理结果起着决定性作用，利用粒子群算法自适应地搜索、优化 MCKD 算法的影响参数是行之有效的。 (

44、3)利用稀疏编码收缩算法对解卷积信号进行处理，可实现非高斯故障冲击信号与高斯噪声干扰信号的有效隔离，达到降噪目的。(4)实例分析及对比结果表明，最大相关峭度解卷积结合稀疏编码收缩的微弱故障特征提取方法可实现齿轮早期故障诊断，并具有一定可靠性和优势。参考文献：1 梅检民，肖云魁，贾继德，等.基于改进阶比的变速器微弱故障特征提取J.振动工程学报，2012，25(3)：317-322.MEI Jianmin，XIAO Yunkui，JIA Jide，et al. Weak fault characteristics extraction for automobile transmission bas

45、ed on improved order analysisJ.Journal of Vibration Engineering，2012，25(3) ：317-322.2 雷亚国,韩冬,林京，等.自适应随机共振新方法及其在故障诊断过中的应用J.机械工程学报，2012，48(7)：62-67.LEI Yaguo，HAN Dong，LIN Jing，et al. New adaptive stochastic resonance method and its application to fault diagnosisJ.Journal of Mechanical Engineering，2012

46、，48(7) ：62-67.3 高立新, 吴丽娟, 张建宇.基于EMD解调方法的齿轮早期故障诊断J.北京工业大学学报，2009，35(7)：876881.GAO Lixin, WU Lijuan, ZHANG Jianyu. An application in gear incipient failure diagnosis based on EMD demodulation methodJ.Journal of Beijing University of Technology，2009,35(7):876 881.4 杨洁明，熊诗波.小波包分析方法在齿轮早期故障特征提取中的应用J.振动、测试

47、与诊断，2000，20(4)：269-272.YANG Jieming，XIONG Shibo. Extraction of early fault information of gears using wavelet packetJ.Journal of Vibration, Measurement&Diagnosis，2000，20(4)：269-272.5 张文义，于德介，陈向民.基于信号共振稀疏分解与包络谱的齿轮故障诊断J.中国机械工程，2013，24(24)：3349-3354.ZHANG Wenyi，YU Dejie，CHEN Xiangmin. Fault diagnosis o

48、f gears based on resonance-based sparse signal decomposition and envelope spectrumJ. China Mechanical Engineering，2013，24(24)：3349-3354.6 SARAVANAN N, RAMACHANDRAN K I. Incipient gear box fault diagnosis using discrete wavelet transform(DWT) for feature extraction and classification using artificial

49、 neural network(ANN)J.Expert Systems with Applications，2010，37(6)：4168-4181. 7 梅检民，肖云魁，沈虹，等. 多阶 FRFT 自适应滤波及齿轮微弱故障特征提取J.振动工程学报，2013，26(5)：771-778.MEI Jianmin，XIAO Yunkui，SHEN Hong，et al. Multi-order FRFT adaptive filter and its application to extract weak fault feature based on sparse signal decomposition J.Jour