基于EEMD（集合经验模态分解）对故障信号进行分解

一、EEMD分解核心流程

graph TD A[原始信号] --> B[添加高斯白噪声] B --> C{多次EMD分解} C --> D[提取IMF分量] D --> E[噪声平均] E --> F[最终IMF]

1. 算法步骤详解

1. 噪声注入

   向原始信号x(t)添加高斯白噪声wi(t)，生成扰动信号xi(t)=x(t)+wi(t)

   noise_level = 0.2; % 噪声幅值系数（通常为信号标准差的0.1-0.3倍）
   noisy_signal = x + noise_level*randn(size(x));
   

2. 多次分解

   对加噪信号进行N次EMD分解（推荐N=50−100次），得到多组IMF分量

   num_trials = 100; % 分解次数
   all_IMFs = cell(num_trials,1);
   for i=1:num_trials[all_IMFs{i}, ~] = emd(noisy_signal);
   end
   

3. 分量平均

   对相同阶次的IMF分量取均值，消除噪声影响

   final_IMFs = zeros(size(all_IMFs{1}));
   for i=1:size(all_IMFs{1},1)final_IMFs(i,:) = mean(cell2mat(arrayfun(@(x) x(i,:), all_IMFs, 'UniformOutput', false)));
   end
   

---

二、关键参数优化

参数	推荐范围	影响分析
噪声幅值	0.1-0.3σ	过小无法抑制模态混叠，过大引入噪声
分解次数	50-200	增加次数提高稳定性，但计算量上升
停止阈值	1e-5-1e-6	过高导致IMF不完整，过低残留噪声
包络线插值	三次样条	线性插值易产生过冲，样条更平滑

---

三、工程应用案例

1. 轴承故障诊断

• 信号特征：内圈故障信号包含周期性冲击（频率ffault=40Hz）

• 分解结果：

  % 加载轴承振动信号（采样率4kHz）
  [x,fs] = load_bearing_signal('inner_race_fault.mat');% EEMD分解
  [IMFs, res] = eemd(x, 0.2, 100);% 提取故障特征IMF（第3阶）
  fault_IMF = IMF(:,3);
  

• 时频分析：通过Hilbert变换提取瞬时频率，定位故障频率成分

2. 齿轮箱状态监测

• 多分量分离：分解出包含边频带的IMF分量（图3）

  % 绘制IMF时频图
  figure;
  subplot(3,1,1);
  spectrogram(IMF(:,1),256,250,256,fs);
  title('IMF1 - 低频噪声');subplot(3,1,2);
  spectrogram(IMF(:,2),256,250,256,fs);
  title('IMF2 - 齿轮啮合频率');subplot(3,1,3);
  spectrogram(IMF(:,3),256,250,256,fs);
  title('IMF3 - 故障边频带');
  

---

四、MATLAB实现代码

function [IMFs,res] = eemd(x, noise_level, num_trials)% 输入参数：% x: 原始信号% noise_level: 噪声幅值系数（默认0.2）% num_trials: 分解次数（默认100）N = length(x);IMFs = cell(num_trials,1);% 多次分解for i=1:num_trials% 添加白噪声noise = noise_level*std(x)*randn(N,1);noisy_x = x + noise;% 执行EMD分解imf = emd(noisy_x);IMFs{i} = imf;end% 平均处理IMFs_avg = zeros(size(IMFs{1}));for i=1:size(IMFs{1},1)IMFs_avg(i,:) = mean(cell2mat(arrayfun(@(x) x(i,:), IMFs, 'UniformOutput', false)));end% 残差计算res = x - sum(IMFs_avg,1);
end% 示例调用
[x,fs] = load_signal('fault_signal.mat');
[IMFs,res] = eemd(x,0.2,100);

参考代码 使用EEMD对故障信号进行分解 www.youwenfan.com/contentcnj/77954.html

五、常见问题解决方案

1. 端点效应

   • 采用镜像延拓法：将信号两端对称延拓

   x_extended = [flipud(x(1:100)), x, flipud(x(end-99:end))];
   

2. 计算效率优化 并行计算：使用parfor替代循环 GPU加速：将信号转换为gpuArray

3. 分解质量验证 检查IMF的平稳性：计算各阶IMF的Hurst指数（应接近0.5） 能量分布验证：故障能量应集中在特定IMF分量

---

六、参考文献

1. 王锡凡. 电力系统可靠性分析基础[M]. 清华大学出版社, 2018.
2. 张琛等. 基于EEMD奇异值熵的滚动轴承故障诊断[J]. 振动．测试与诊断, 2019.
3. MathWorks官方文档: Empirical Mode Decomposition Toolbox

通过合理设置参数并配合特征提取算法，EEMD可有效提取故障信号中的微弱特征，为复杂工业设备状态监测提供可靠支持。实际应用中需根据信号特性调整噪声水平和分解次数，并通过交叉验证确保分解质量。