aa.zip_小波分解_小波变换

clc;close all;clear all; X = csvread('tek0001CH1.csv',14,1); N=1000; t=1:N; figure(1); subplot(211); plot(X,'linewidth',2); xlabel('样本序号 N'); ylabel('幅值A'); %叠加信号 [c,l]=wavedec(X,7,'db5');%重构第1-7层逼近系数 a7=wrcoef('a',c,l,'db5',7); a6=wrcoef('a',c,l,'db5',6); a5=wrcoef('a',c,l,'db5',5); a4=wrcoef('a',c,l,'db5',4); a3=wrcoef('a',c,l,'db5',3); a2=wrcoef('a',c,l,'db5',2); a1=wrcoef('a',c,l,'db5',1);%显示逼近系数 figure(3) subplot(711) plot(a7,'linewidth',2); ylabel('a7'); subplot(712) plot(a6,'linewidth',2); ylabel('a6'); subplot(713) plot(a5,'linewidth',2); ylabel('a5'); subplot(714) plot(a4,'linewidth',2); ylabel('a4'); subplot(715) plot(a3,'linewidth',2); ylabel('a3'); subplot(716) plot(a2,'linewidth',2); ylabel('a2'); subplot(717) plot(a1,'linewidth',2); ylabel('a1'); xlabel('样本序号 N');%重构第1-7层细节系数 d7=wrcoef('d',c,l,'db5',7); d6=wrcoef('d',c,l,'db5',6); d5=wrcoef('d',c,l,'db5',5); d4=wrcoef('d',c,l,'db5',4); d3=wrcoef('d',c,l,'db5',3); d2=wrcoef('d',c,l,'db5',2); d1=wrcoef('d',c,l,'db5',1); %显示细节系数 figure(4) subplot(711) plot(d7,'linewidth',2); ylabel('d7'); subplot(712) plot(d6,'linewidth',2); ylabel('d6'); subplot(713) plot(d5,'linewidth',2); ylabel('d5'); subplot(714) plot(d4,'linewidth',2); ylabel('d4'); subplot(715) plot(d3,'linewidth',2); ylabel('d3'); subplot(716) plot(d2,'linewidth',2); ylabel('d2'); subplot(717) plot(d1,'linewidth',2); ylabel('d1'); xlabel('样本序号 N'); ylabel('幅值A'); %叠加信号 %%%%%%%%% 绘制时间域小波波形 %%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %cparam=2*fc*totalscal; %%scal=cparam./a; %coefs=cwt(s,scal,wavename,1/fs); %f=scal2frq(scal,wavename,1/fs); %subplot(212) %imagesc(t,f,abs(coefs)); %colormap(jet) %colorbar %xlabel('时间 t/s') %ylabel('频率 f/Hz') %title('小波时频图')