matlab 轴承故障诊断代码

以下是一份简单的 MATLAB 代码，用于诊断轴承故障：

clear all
close all

% 读取数据
load('bearing_data.mat');
t = bearing_data(:,1);
x = bearing_data(:,2);

% 绘制原始信号
figure;
plot(t,x)
xlabel('时间 (s)')
ylabel('加速度 (g)')
title('原始信号')

% 频谱分析
L = length(x);
Fs = L/(t(end)-t(1));
f = Fs*(0:(L/2))/L;
X = fft(x-mean(x));
P = abs(X/L).^2;
P2 = P(1:L/2+1);
P2(2:end-1) = 2*P2(2:end-1);

% 绘制频谱图
figure;
plot(f,P2)
xlabel('频率 (Hz)')
ylabel('功率谱密度')
title('频谱图')

% 提取特征
[~,I] = max(P2);
f_peak = f(I);
[~,locs] = findpeaks(x,'MinPeakHeight',0.1,'MinPeakDistance',1000);
n_peaks = length(locs);

% 判断故障类型
if f_peak > 200 && n_peaks > 5
    fprintf('滚动体故障');
elseif f_peak > 100 && f_peak < 200 && n_peaks > 10
    fprintf('内圈故障');
elseif f_peak < 100 && n_peaks > 10
    fprintf('外圈故障');
else
    fprintf('正常');
end

注意，这只是一个简单的示例代码，需要根据具体情况进行修改和优化。同时，还需要结合实际情况选择合适的特征提取方法和故障诊断算法。