小波分析故障

% 小波时频分析 clc clear % 原始信号 fs=2500; t=0:1/fs:1.2; load('Sig1');s=x(:,2);load('Sig2');s1=x1(:,2);load('Sig3');s2=x2(:,2); % 连续小波变换 wavename='cmor3-3'; totalscal=256; Fc=centfrq(wavename); % 小波的中心频率 c=2*Fc*totalscal; scals=c./(1:totalscal); f=scal2frq(scals,wavename,1/fs); % 将尺度转换为频率 coefs=cwt(s,scals,wavename); coefs1=cwt(s1,scals,wavename);coefs2=cwt(s2,scals,wavename);% 求连续小波系数 figure subplot(311);imagesc(t,f,abs(coefs));set(gca,'YDir','normal');colorbar; xlabel('时间 t/s');ylabel('频率 f/Hz');title('正常信号小波时频图'); subplot(312);imagesc(t,f,abs(coefs1));set(gca,'YDir','normal');colorbar; xlabel('时间 t/s');ylabel('频率 f/Hz');title('故障信号小波时频图'); subplot(313);imagesc(t,f,abs(coefs2));set(gca,'YDir','normal');colorbar; xlabel('时间 t/s');ylabel('频率 f/Hz');title('严重故障信号小波时频图');