%?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