Page 131 - 《振动工程学报》2025年第8期

第 8 期 姚容华，等： 奇异谱分解和最大相关峭度解卷积在轴承故障声学诊断中的应用 1771

图 13 信号分解及包络谱
图 11 分量包络峭度值 Fig. 13 Signal decomposition and envelope spectrum
Fig. 11 Component envelope kurtosis values
如图 14 所示，图 14（a）为多点峭度谱计算过程，图 14
（b）为峰度计算过程，计算得到 MCKD 的最佳周期，
内圈为 T o = 87，外圈为 T i = 127。
其次计算滤波器长度，本文选择包络排列熵为
目标对长度进行计算。定义迭代区间为［0，1000］，
包 络 排 列 熵 最 值 对 应 的 长 度 为 滤 波 器 最 优 长 度 。
图 12 故障特征能量幅值比 计算过程如图 15 所示，图 15（a）为内圈包络排列熵
Fig. 12 Fault characteristic energy amplitude ratio
计 算 过 程 ，图 15（b）为 外 圈 包 络 排 列 熵 计 算 过 程 ，
计 算 得 到 内 圈 的 长 度 为 L o = 536，外 圈 的 长 度 为
晰看出故障特征频率及其倍频。为了突出 SSD 的
分解能力，与 EMD 进行对比。图 13（b）为 EMD 分
解包络谱，通过与 OSSD 分解结果进行对比，EMD
的分解结果中内、外圈故障信息被大量噪声包围，很
难 用 于 后 续 处 理 ，结 果 证 明 SSD 的 分 解 结 果 优 于
EMD。
对敏感分量进行滤波和增强冲击，选择 MCKD
对信号进行滤波。利用本文所提方法对 MCKD 进
行优化。首先计算滤波器的周期，利用多点峭度谱
确定 MCKD 的滤波周期 T 的估计值，确定估计值后
缩小区间，以峰度为目标计算出精确值。计算过程

图 14 周期 T 的计算过程
Fig. 14 Calculation process of cycle T