小波分析对轴承进行故障诊断

%采样频率
fs=10000;
%轴承滚动体故障信号
load bearingout.mat
fid=fopen（‘bearingout.mat‘‘r‘）;%故障
xdata=fread（fid1024‘int16‘）;
fclose（fid）;
xdata=（xdata-mean（xdata））/std（xdata1）;
 
%时域波形
figure（1）;
plot（1:1024xdata）;
xlabel（‘时间 t/n‘）;
ylabel（‘电压 V/v‘）;
 
%db10小波进行4层分解
%一维小波分解
[cl] = wavedec（xdata4‘db10‘）;
%重构第1～4层细节系数
d4 = wrcoef（‘d‘cl‘db10‘4）;
d3 = wrcoef（‘d‘cl‘db10‘3）;
d2 = wrcoef（‘d‘cl‘db10‘2）;
d1 = wrcoef（‘d‘cl‘db10‘1）;