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