【故障诊断】基于最小熵反卷积、最大相关峰度反卷积和最大二阶环平稳盲反卷积等盲反卷积方法在机械故障诊断中的应用研究（Matlab代码实现）

目录

💥1 概述

📚2 运行结果

2.1 稀疏最大谐波噪声比反卷积

2.2 最大相关峰度反卷积

2.3 最大二阶环平稳盲反卷积

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

在机械故障诊断中，盲反卷积方法可以用于信号处理和故障特征提取。最小熵反卷积、最大相关峰度反卷积和最大二阶环平稳盲反卷积是几种常用的盲反卷积技术。它们通过对故障信号进行盲反卷积，可以恢复出信号中的故障特征，帮助识别和定位故障。

最小熵反卷积是一种非线性逆滤波方法，通过最小化重构信号的熵来恢复原始信号。它在处理非线性失真和信号叠加时表现出较好的性能。

最大相关峰度反卷积是一种基于峭度（kurtosis）准则的盲反卷积方法，将观测信号的峭度最大化，来寻找潜在信号的时延和幅度。最大相关峰度反卷积在处理非平稳信号和部分遮挡信号时具有较好的鲁棒性。

最大二阶环平稳盲反卷积是一种在二阶统计特性约束下的盲反卷积方法，通过最大化信号的二阶平稳度来恢复原始信号。它能有效处理多路径传输和多时变信号的特点。

这些盲反卷积方法在故障诊断中的应用可以帮助提取故障信号中的关键特征，例如冲击、共振等，进而对故障进行定位和诊断。然而，需要注意的是，在实际应用中，选择合适的盲反卷积方法需要根据具体的故障类型、噪声情况和信号特点进行权衡和判断，以达到更好的诊断效果。

最小熵反卷积（MED）、最大相关峰度反卷积（MCKD）和最大二阶环平稳盲反卷积（CYCBD）等盲反卷积方法在机械故障诊断中得到了广泛的应用。同时考虑周期性和脉冲性，MCKD，CYCBD可以解决MED偏向于随机脉冲而不是周期性故障脉冲的问题。然而，MCKD和CYCBD的优越性在很大程度上取决于先前的故障周期。在工业应用中，由于转速波动和测量问题，很难准确获得故障周期。

因此，我们首先提出利用迭代算法估计迭代周期，以解决盲反卷积方法中前期的问题。根据自相关原理，即当时延满足周期或其倍数时，将显示较高的值，在MCKD中选择具有局部最大值的位置作为迭代周期。包络谐波积谱（EHPS）最初是为估计CYCBD中的特征频率而定制的。基于BDM迭代算法的周期估计可以帮助MCKD和CYCBD在无需先验知识的情况下应用于机械故障诊断。

📚2 运行结果

2.1 稀疏最大谐波噪声比反卷积

 

 

2.2 最大相关峰度反卷积

2.3 最大二阶环平稳盲反卷积

 

 

 

 

 部分代码：

%%
load sig3
x = x - mean(x);
addpath('..\00 subfunction\')

%%
fs = 20000;
N = length(x);
t = (0:N - 1) / fs;
t = t(:);
BPFI = 38;

%% Raw data
figure;
plot(t, x, 'b');
xlabel('Time [s]')
ylabel('Amplitude')
title('Raw data')
legend(['Kurtosis=', num2str(kurtosis(x))])
setfontsize(20);
set(gcf, 'position', [100, 100, 800, 400])
axis tight
ylim([-2 2.5])

envelope_x = abs(hilbert(x)) - mean(abs(hilbert(x)));
ff = 0:fs / N:fs - fs / N;
amp_envelope_x = abs(fft(envelope_x, N)) * 2 / fs;
figure;
plot(ff, amp_envelope_x, 'b')
xlabel('Frequency [Hz]')
ylabel('Amplitude')
setfontsize(20);
set(gcf, 'position', [100, 100, 800, 400])
axis tight
xlim([0, 200]);
ylim([0 0.025])

%% SMHD

[y_final, f_final, kurtIter] = smhd(fs, x, 100, 30, 1.5 * rms(x), [], 0);

%% Filtered signal
figure;
plot(t, y_final, 'b');
xlabel('Time [s]')
ylabel('Amplitude')
title('Filtered signal by SMHD')
legend(['Kurtosis=', num2str(kurtosis(y_final))])
setfontsize(20);
set(gcf, 'position', [100, 100, 800, 400])
axis tight
ylim([-3.5 4.5])

envelope_y = abs(hilbert(y_final)) - mean(abs(hilbert(y_final)));
amp_envelope_y = abs(fft(envelope_y, N)) * 2 / fs;
figure;
plot(ff, amp_envelope_y, 'b')
xlabel('Frequency [Hz]')
ylabel('Amplitude')
setfontsize(20);
set(gcf, 'position', [100, 100, 800, 400])
axis tight
xlim([0, 200]);
ylim([0 0.3])
 

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1] 苗彦， 赵明， 林建军， 雷尹， 改进的最大相关峰度反卷积方法在滚动体轴承故障诊断中的应用， 机械系统与信号处理， 92 （2017） 173-195.

[2] 苗彦， 赵明， 梁国强， 林建军， 改进MCKDA在基于编码器信号的风力发电机齿轮故障检测中的应用， 可再生能源， 151 （2020） 192-203.

[3] 张斌， 苗彦， 林建军， 易玲玲， 自适应最大二阶环平稳性盲反卷积及其在机车轴承故障诊断中的应用， 机械系统与信号处理， 158 （2021） 107736.

[4] 苗彦， 张斌， 林建军， 赵明， 刘海， 刘忠， 李海， 盲反卷积在机械故障诊断中的应用研究进展， 机械系统与信号处理， 163 （2022） 108202.

[5] 苗彦， 赵明， 林建， 徐晓， 用于轴承弱故障特征检测的稀疏最大谐波噪声比反卷积， 测量科学与技术， 27 （2016） 105004.

[6] 苗彦， 张斌， 林建军， 赵明， 刘海， 刘忠， 李海， 盲反卷积在机械故障诊断中的应用研究进展， 机械系统与信号处理， 163 （2022） 108202.

[1] Y. Miao, M. Zhao, J. Lin, Y. Lei, Application of an improved maximum correlated kurtosis deconvolution method for fault diagnosis of rolling element bearings, Mechanical Systems and Signal Processing, 92 (2017) 173-195.

[2] Y. Miao, M. Zhao, K. Liang, J. Lin, Application of an improved MCKDA for fault detection of wind turbine gear based on encoder signal, Renewable Energy, 151 (2020) 192-203.

The matlab codes of period estimation using EHPS permit to reproduce some results in the papers:

[1] B. Zhang, Y. Miao, J. Lin, Y. Yi, Adaptive maximum second-order cyclostationarity blind deconvolution and its application for locomotive bearing fault diagnosis, Mechanical Systems and Signal Processing, 158 (2021) 107736.

[2] Y. Miao, B. Zhang, J. Lin, M. Zhao, H. Liu, Z. Liu, H. Li, A review on the application of blind deconvolution in machinery fault diagnosis, Mechanical Systems and Signal Processing, 163 (2022) 108202.

In addition, the matlab codes of the deconvolution method, Sparse maximum harmonics-to-noise-ratio deconvolution (SMHD), permit to reproduce some results in the papers:

[1] Y. Miao, M. Zhao, J. Lin, X. Xu, Sparse maximum harmonics-to-noise-ratio deconvolution for weak fault signature detection in bearings, Measurement Science and Technology, 27 (2016) 105004.

[2] Y. Miao, B. Zhang, J. Lin, M. Zhao, H. Liu, Z. Liu, H. Li, A review on the application of blind deconvolution in machinery fault diagnosis, Mechanical Systems and Signal Processing, 163 (2022) 108202.

🌈4 Matlab代码实现