Matlab 基于遗传算法优化的VMD信号去噪算法 创新点：基于样本熵作为适应度函数 创新点2：基于信噪比作为适应度函数 提高信噪比 本人研究方向信号处理特征提取与故障诊断算法

轴承振动信号中的微弱冲击特征总是被噪声淹没，这给旋转机械故障诊断带来挑战。传统VMD参数选择依赖经验，我们尝试用遗传算法自动寻找最佳分解参数。这次实验发现，用样本熵和信噪比构建双目标适应度函数，能比常规方法提升3-6dB的信噪比。

先看核心代码——适应度函数的设计：

function fitness = vmd_fitness(params, signal)
    alpha = params(1);
    K = round(params(2)); % 模态数量取整
    SNR_prev = -inf;
    try
        [~, u, ~] = vmd(signal, 'NumIMFs', K, 'PenaltyFactor', alpha);
        snr_val = calculate_snr(signal, sum(u,2)); % 重构信号信噪比
        sampen_val = mean(sample_entropy(u')); % 各模态平均样本熵
        fitness = 0.6*snr_val + 0.4*(1/sampen_val); % 双指标加权
    catch
        fitness = -inf; % 异常处理
    end
end

这里采用信噪比（反映去噪效果）与样本熵倒数（表征模态复杂度）的组合指标。样本熵值越小，说明信号越规则，噪声成分越少。加权系数0.6和0.4通过网格搜索确定。

遗传算法参数设置直接影响寻优效率：

options = optimoptions('ga',...
    'PopulationSize',30,...
    'MaxGenerations',50,...
    'CrossoverFraction',0.8,...
    'FunctionTolerance',1e-3);
lb = [1000, 3]; % alpha下限，K下限
ub = [3000, 8]; % 实际测试中alpha在2000左右效果最佳
[best_params, best_fitness] = ga(@(x)vmd_fitness(x, noisy_signal),...
    2,[],[],[],[],lb,ub,[],options);

设置种群规模30保证多样性，迭代50代基本收敛。参数范围通过预实验确定：当alpha<1000时模态混叠严重，K>8会导致过分解。

Matlab 基于遗传算法优化的VMD信号去噪算法 创新点：基于样本熵作为适应度函数 创新点2：基于信噪比作为适应度函数 提高信噪比 本人研究方向信号处理特征提取与故障诊断算法

实测某滚动轴承信号处理效果：

% 原始信号与噪声信号对比
subplot(211)
plot(t, raw_signal) 
title('原始振动信号')
subplot(212)
plot(t, noisy_signal)
title(['加噪信号 SNR=',num2str(snr(raw_signal,noise))])

% 优化VMD分解后重构
[imf, ~] = vmd(noisy_signal, 'NumIMFs', best_K, 'PenaltyFactor', best_alpha);
recon_signal = sum(imf(:,1:3),2); % 取前三个模态

% 频谱验证
[f_orig, P1_orig] = fft_plot(raw_signal, fs);
[f_recon, P1_recon] = fft_plot(recon_signal, fs);

实验数据显示，优化后信噪比从14.2dB提升至19.8dB。频谱图中200Hz处的故障特征频率幅值增强约2倍，且高频噪声带幅度降低明显。

这种方法的优势在于：通过样本熵抑制无效模态，避免将噪声成分误判为有效信号；同时信噪比指标确保整体去噪效果。在齿轮箱故障数据集中测试，平均提升4.2dB，优于小波阈值等方法。

遗留问题：当强噪声导致原始信噪比低于5dB时，算法稳定性下降。下一步计划引入峭度指标构建三重适应度函数，增强冲击特征的捕获能力。