基于FFT实现轴承故障诊断含Matlab代码

% ********** 内圈故障300r/min ********** clear; clc; t=0:1:16383; %采样步长 fs=10000; %采样频率 y=importdata('neiquan_10k_32768_5.dat'); %读取数据 N=length(t); z=y-mean(y); %零均值处理 figure(1); subplot(211); plot(t,z); %原始数据图形 title('内圈300r/min');xlabel('采样数');ylabel('幅值'); subplot(212); f=fft(z,N); plot((0:N/2-1)/N*fs,abs(f(1:N/2))); title('FFT变换');xlabel('频率/Hz');ylabel('幅值'); % ********** 内圈故障450r/min ********** clear; clc; t=0:1:16383; %采样步长 fs=10000; %采样频率 y=importdata('neiquan_10k_32768_7.5.dat'); %读取数据 N=length(t); z=y-mean(y); %零均值处理 figure(2); subplot(211); plot(t,z); %原始数据图形 title('内圈450r/min');xlabel('采样数');ylabel('幅值'); subplot(212); f=fft(z,N); plot((0:N/2-1)/N*fs,abs(f(1:N/2))); title('FFT变换');xlabel('频率/Hz');ylabel('幅值'); % ********** 外圈故障300r/min ********** clear; clc; t=0:1:16383; %采样步长 fs=10000; %采样频率 y=importdata('waiquan_10k_32768_5.dat'); %读取数据 N=length(t); z=y-mean(y); %零均值处理 figure(3); subplot(211); plot(t,z); %原始数据图形 title('外圈300r/min');xlabel('采样数');ylabel('幅值'); subplot(212); f=fft(z,N); plot((0:N/2-1)/N*fs,abs(f(1:N/2))); %包络谱 title('FFT变换');xlabel('频率/Hz');ylabel('幅值'); % ********** 外圈故障450r/min ********** clear; clc; t=0:1:16383; %采样步长 fs=10000; %采样频率 y=importdata('waiquan_10k_32768_7.5.dat'); %读取数据 N=length(t); z=y-mean(y); %零均值处理 figure(4); subplot(211); plot(t,z); %原始数据图形 title('外圈450r/min');xlabel('采样数');ylabel('幅值'); subplot(212); f=fft(z,N); plot((0:N/2-1)/N*fs,abs(f(1:N/2))); title('FFT变换');xlabel('频率/Hz');ylabel('幅值');