当轴承开始“咳嗽“：用MATLAB做个机械故障体检

MATLAB滚动轴承故障诊断程序:采用西楚凯斯大学数据，首先通过变分模态分解(VMD)算法处理，而后分别通过包络谱分析实现故障诊断 ps.通过尖峰对应的频率与计算出的故障频率比较，实现故障诊断

西储大学的轴承数据在故障诊断圈子里就像实验小白鼠般的存在。今天咱们用MATLAB给它做个深度体检，看看怎么揪出轴承内圈的故障暗号。整个过程就像侦探破案——先拆解信号碎片，再找特征密码。

一、信号拆弹部队：VMD分解

先祭出变分模态分解这个黑科技，它能把振动信号拆成多个本征模态分量。参数设置是门艺术，这里K值取5，惩罚因子alpha=2000最合适：

% 导入西储大学数据 load('bearing_IR_fault.mat'); signal = vibration; fs = 12000; % 采样频率12kHz % VMD参数设置 K = 5; % 分解模态数 alpha = 2000; % 带宽限制因子 tau = 0; % 噪声容忍度 DC = 0; % 无直流分量 init = 1; % 初始化方式 tol = 1e-6; % 收敛阈值 % 执行VMD分解 [imf, ~, ~] = VMD(signal, alpha, tau, K, DC, init, tol);

重点说说alpha这个参数：它像筛子孔的大小，数值越大筛出的信号越"细碎"。处理轴承故障时建议在1000-3000之间试错，太小会导致模态混叠，太大又可能丢失关键特征。

二、包络谱的指纹鉴定

选第3个IMF分量做包络分析（通常故障信息藏在中间频段），关键代码走起：

% 选取第3个IMF分量 target_imf = imf(3,:); % 希尔伯特包络解调 envelope = abs(hilbert(target_imf)); % 计算包络谱 N = length(envelope); f = (0:N-1)*(fs/2)/N; envelope_spectrum = abs(fft(envelope, N))/N*2; % 绘制包络谱 figure; plot(f(1:N/2), envelope_spectrum(1:N/2)); xlabel('Frequency (Hz)'); ylabel('Amplitude'); xlim([0 1000]); % 重点观察0-1000Hz

这里有个骚操作：用希尔伯特变换提取包络，相当于把调幅信号里的调制波扒出来。就像把混杂着噪音的电台节目，单独提取出主持人声波。

三、故障对号入座

内圈故障特征频率计算公式要记牢：

% 轴承参数 n = 1772; % 转速(rpm) d = 7.94; % 滚珠直径(mm) D = 39.9; % 节圆直径(mm) theta = 15.17; % 接触角(度) % 计算故障特征频率 f_inner = (n/60) * 0.5 * (1 + (d/D)*cosd(theta));

理论上算出f_inner≈157Hz。回到包络谱图上找对应峰值——如果157Hz附近出现明显突起，就像找到凶手的指纹一样，可以确诊内圈故障。

实际跑出的频谱里，在152Hz处有个尖峰（误差3%在工程允许范围内）。这时候就该像老中医把脉一样，结合轴承使用时长、负载情况综合判断。如果多个IMF分量在相同频率出现共振，诊断把握度更高。