%采樣頻率
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）;