load freqbrk;
s=freqbrk;
[c,l]=wavedec(s,6,'db5');
subplot(8,1,1);
plot(s);
title('使用db5小波分解6層：s=a6+d6+d5+d4+d3+d2+d1');
Ylabel('s');
a6=wrcoef('a',c,l,'db5',6);
subplot(8,1,2);
plot(a6);
Ylabel('a6');
d6=wrcoef('d',c,l,'db5',6);
subplot(8,1,3);
plot(d6);
Ylabel('d6');
d5=wrcoef('d',c,l,'db5',5);
subplot(8,1,4);
plot(d5);
Ylabel('d5');
d4=wrcoef('d',c,l,'db5',4);
subplot(8,1,5);
plot(d4);
Ylabel('d4');
d3=wrcoef('d',c,l,'db5',3);
subplot(8,1,6);
plot(d3);
Ylabel('d3');
d2=wrcoef('d',c,l,'db5',2);
subplot(8,1,7);
plot(d2);
Ylabel('d2');
d1=wrcoef('d',c,l,'db5',1);
subplot(8,1,8);
plot(d1);
Ylabel('d1');