?clear
t=0:pi/100:4*pi;
s=sin（t+pi/4）;
subplot（532）;
plot（s）;
grid;
title（‘原始信號‘）;


[cl]=wavedec（s3‘db1‘）;
grid;
ca3=appcoef（cl‘db1‘3）;%提取小波分解的低頻系數　　
cd3=detcoef（cl3）;%提取第三層的高頻系數　　　　
cd2=detcoef（cl2）;%提取第二層的高頻系數
cd1=detcoef（cl1）;%提取第一層的高頻系數　
figure（2）;
subplot（421）;
plot（ca3）;
title（‘第三層低頻系數‘）;
subplot（423）;
plot（cd1）;
title（‘第一層高頻系數‘）;
subplot（425）;
plot（cd2）;
title（‘第二層高頻系數‘）;
subplot（427）;
plot（cd3）;
title（‘第三層高頻系數‘）;
s1=waverec（cl‘db1‘）;

[cl]=wavedec（s3‘coif3‘）;