文献阅读-Ensemble Empirical Mode Decomposition-Based Teager Energy Spectrum for Bearing Fault Diagnostic

本文介绍了北科大冯志鹏老师的一项研究,提出了一种结合EEMD和Teager能量算子的轴承故障诊断方法。通过EEMD分解振动信号,选取相关性强且峰度高的本征模态函数(IMF),再用Teager能量算子检测故障脉冲,最后通过傅立叶变换识别故障重复频率,有效提取弱故障和复合故障特征。实验表明,这种方法优于传统谱分析,在轴承故障诊断中表现出色。
摘要由CSDN通过智能技术生成

引言:
这篇文章是北科大-冯志鹏老师的早期研究工作,工作主要是面向轴承的故障诊断问题,创新性的提出了一种基于EEMD的方法用于微弱故障信号特征提取。


Ensemble Empirical Mode Decomposition-Based Teager Energy Spectrum for Bearing Fault Diagnostics

1. 摘要

振动信号中的周期性脉冲及其重复频率是诊断滚动轴承局部损伤的关键指标。提出了一种基于集成经验模态分解(EEMD)和Teager能量算子的轴承故障特征频率提取方法。首先利用Teager能量算子对信号进行EEMD分解,以满足单分量要求。然后,根据本征模态函数与原始信号的相关性及其峰度来选择感兴趣的本征模态函数(IMF)。接下来,将Teager能量算子应用于选定的IMF,以检测故障诱发的脉冲。最后,将傅立叶变换应用于获得的Teager能量序列,以识别故障引起的周期性脉冲的重复频率,从而诊断轴承故障。通过对模拟轴承振动信号的分析,说明了该方法的原理。通过种子故障试验和故障运行试验的实验信号分析,验证了该方法在提取轴承故障特征频率方面的有效性,尤其是在**识别弱故障和复合故障症状方面的性能。**比较研究表明,该方法比传统的谱分析方法和平方包络谱分析方法具有更好的性能,或者说是对传统谱分析方法的补充。

2. 方法建模分析部分

2.1 动力学方程部分

这部分算是一点题外话

以下是计算原始信号减去EEMD、SSA和LS季节模型所产生的残差的频谱指数的matlab代码: ```matlab % 假设原始信号为signal,采样率为fs % EEMD分解 imf = ceemdan(signal); % 进行CEEMDAN分解,得到各个本征模态函数imf residual = signal - sum(imf, 2); % 计算残差 % 计算频谱指数 [pxx, f] = pwelch(residual, [], [], [], fs); % 计算残差的功率谱密度pxx和对应频率f fi = find(f > 0 & f < 0.1); % 选择感兴趣的频率区间 SI_EEMD = trapz(pxx(fi)) / trapz(pxx); % 计算频谱指数SI_EEMD % SSA分解 L = 100; % 选择时间窗口长度 K = floor(length(signal) / L); % 计算分解后的子序列数 X = zeros(L, K); for i = 1:K X(:, i) = signal((i-1)*L+1:i*L); end imf = ssa(X); % 进行SSA分解,得到各个本征模态函数imf residual = signal - sum(imf, 2); % 计算残差 % 计算频谱指数 [pxx, f] = pwelch(residual, [], [], [], fs); % 计算残差的功率谱密度pxx和对应频率f fi = find(f > 0 & f < 0.1); % 选择感兴趣的频率区间 SI_SSA = trapz(pxx(fi)) / trapz(pxx); % 计算频谱指数SI_SSA % LS季节模型 t = (1:length(signal))'; % 构造时间序列t X = [ones(length(signal), 1), sin(2*pi*t/365), cos(2*pi*t/365)]; % 构造基函数矩阵X beta = X \ signal; % 进行最小二乘拟合,得到拟合系数beta residual = signal - X*beta; % 计算残差 % 计算频谱指数 [pxx, f] = pwelch(residual, [], [], [], fs); % 计算残差的功率谱密度pxx和对应频率f fi = find(f > 0 & f < 0.1); % 选择感兴趣的频率区间 SI_LS = trapz(pxx(fi)) / trapz(pxx); % 计算频谱指数SI_LS ``` 其中,ceemdan和ssa是EEMD和SSA的matlab函数,可在MATLAB File Exchange中下载。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Robo-网络矿产提炼工

你的鼓励将是我最大的动力!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值