轴承振动信号处理

% fs=12000;
% load X156_DE_time;
% xdata=X156_DE_time;
% xdata=xdata（1:10240）;
% N=10240;
% 
% 
% xdata=（xdata-mean（xdata））/std（xdata1）;%mean平均值 std标准偏差   对原始信号进行归一处理
% 
% %时域波形
% figure（1）;
% subplot（211）;
% plot（1:Nxdata）;
% xlabel（‘时间 t/n‘）;
% ylabel（‘电压 V/v‘）;

%采样频率
fs=12000;
%轴承外环故障信号
% fid=fopen（‘b.mat‘‘r‘）;%故障
N=10240;
% xdata=fread（fidN‘int16‘）;%fread：以二进制形式读取文件中数据，int16：16字节数，整数型
% fclose（fid）;
load X156_DE_time;
xdata=X156_DE_time;
xdata=xdata（1:10240）;
xdata=（xdata-mean（xdata））/std（xdata1）;%mean平均值 std标准偏差   对原始信号进行归一处理

%时域波形
figure（1）;
subplot（211）;
plot（1:Nxdata）;
xlabel（‘时间 t/n‘）;
ylabel（‘电压 V/v‘）;
%原始信号的功率谱
nfft=10240;
Y1=abs（fft（xdatanfft））;
subplot（212）;
plot（（0:nfft/2-1）/nfft*fsY1（1:nfft/2））;
xlabel（‘频率 f/Hz‘）;
ylab