采样频率

fs=10000

轴承外环故障信号

fid=fopen( bearingout.dat , r ) 故障

N=1024

xdata=fread(fid,N, int16 )

fclose(fid)

xdata=(xdata-mean(xdata))/std(xdata,1)



时域波形

figure(1)

plot(1:N,xdata)

xlabel( 时间 t/n )

ylabel( 电压 V/v )



db10小波进行4层分解

一维小波分解

[c,l] = wavedec(xdata,4, db10 )



重构第1～4层细节信号

d4 = wrcoef( d ,c,l, db10 ,4)

d3 = wrcoef( d ,c,l, db10 ,3)

d2 = wrcoef( d ,c,l, db10 ,2)

d1 = wrcoef( d ,c,l, db10 ,1)

-　Sampling frequency fs = 10000　 bearing outer ring fault signal fid = fopen (' bearingout.dat' , ' r' )　 failure N = 1024　xdata = fread (fid, N, ' int16' )　fclose (fid )　xdata = (xdata-mean (xdata))/std (xdata, 1)　 time-domain waveform figure (1)　plot (1: N, xdata)　xlabel (' Time t/n' )　ylabel ( ' voltage V/v' )　 db10 wavelet decomposition　4 layer one-dimensional wavelet decomposition [c, l] = wavedec (xdata, 4, ' db10' )　 1 ~ 4 reconstructed detail signal d4 = wrcoef (' d' , c, l, ' db10' , 4)　d3 = wrcoef (' d' , c, l, ' db10' , 3)　d2 = wrcoef (' d' , c, l, ' db10'　, 2)　d1 = wrcoef (' d' , c, l, ' db10' , 1)