t=0:pi/100:4*pi;
s=sin(t+3*pi/4);
plot(s);
title('原始信號s');
%====================================
%對s進行小波分解：db1  3層
[c,l]=wavedec(s,3,'db1');
%====================================
%提取小波分解的低頻系數(shù)a3
a3=appcoef(c,l,'db1',3);
%====================================
%提取小波分解的各層高頻系數(shù)
d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);	
%====================================
%繪出各系數(shù)的圖形
figure(2)
subplot(5,2,1);plot(a3);
Ylabel('a3');
subplot(5,2,3);plot(d3);
Ylabel('d3');
subplot(5,2,5);plot(d2);
Ylabel('d2');
subplot(5,2,7);plot(d1);
Ylabel('d1');
%====================================
%重構(gòu)信號s
s1=waverec(c,l,'db1');
subplot(5,2,9);plot(s1);
Ylabel('s1');
%====================================
%下面用小波‘coif3’重復(fù)上述過程
[c,l]=wavedec(s,3,'coif3');
a3=appcoef(c,l,'coif3',3);
d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);
subplot(5,2,2);plot(a3);
Ylabel('a3');
subplot(5,2,4);plot(d3);
Ylabel('d3');
subplot(5,2,6);plot(d2);
Ylabel('d2');
subplot(5,2,8);plot(d1);
Ylabel('d1');
s2=waverec(c,l,'coif3');
subplot(5,2,10);plot(s2);
Ylabel('s2');