1、fs=30000;%读取正常信号fid=fopen(e:1.dat,r);N=3072;x1=fread(fid,N,int16);fclose(fid);%读取故障信号fid=fopen(e:2.dat,r);x2=fread(fid,N,int16);fclose(fid);t=1:3072;%-figure(1);subplot(211);plot(t,x1);axis(1,1000,-600,600);title(正常信号);subplot(212);plot(t,x2);title(故障信号);axis(1,1000,-600,600);%-figure(2);%对正常齿轮进行2尺度
2、小波分解c,l=wavedec(x1,2,db2);%对正常齿轮进行2尺度小波分解d=wrcoef(d,c,l,db2);%对正常齿轮进行2尺度小波分解a=wrcoef(a,c,l,db2);subplot(221);plot(d);title(正常信号的高频重构);axis(1,500,-2000,2000);subplot(222);plot(a);title(正常信号的低频重构);axis(1,500,-2000,2000);%对故障齿轮进行2尺度小波分解c,l=wavedec(x2,2,db2);%对故障齿轮进行2尺度小波分解d=wrcoef(d,c,l,db2);%对故障齿轮进行2尺
3、度小波分解a=wrcoef(a,c,l,db2);subplot(223);plot(d);title(故障信号的高频重构);axis(1,500,-2000,2000);subplot(224);plot(a);title(故障信号的低频重构);axis(1,500,-2000,2000);%-c,l=wavedec(x1,2,db2);d1=wrcoef(d,c,l,db2,1);y=hilbert(d1);%对正常信号第一层细节重构信号进行希尔伯特变换y=abs(y);nfft=3072;%求正常信号波洛谱p=abs(fft(y,nfft);figure(3);subplot(211);
4、stem(0:nfft/2-1)/nfft*fs,p(1:nfft/2);axis(1000,1500,1,300000);title(正常信号第一层细节信号包络谱);xlabel(频率f/Hz);ylabel(功率谱);c,l=wavedec(x2,2,db2);d1=wrcoef(d,c,l,db2,1);%对故障信号第一层细节重构信号进行希尔伯特变换y=hilbert(d1);%求正常信号波洛谱y=abs(y);nfft=3072;p=abs(fft(y,nfft);subplot(212);stem(0:nfft/2-1)/nfft*fs,p(1:nfft/2);axis(1000,1500,1,300000);title(故障信号第一层细节信号包络谱);xlabel(频率f/Hz);ylabel(功率谱);
