load  noissin;
x = noissin;
figure(1);
subplot(611);
plot(x);
ylabel('x');
% 使用db5小波對信號進行5層分解
[c,l]=wavedec(x,5,'db5');
for i=1:5
    % 對分解的第5層到第1層的低頻系數近行重構
    a = wrcoef('a',c,l,'db5',6-i);
    subplot(6,1,i+1);
    plot(a);
    ylabel(['a',num2str(6-i)]);
end
figure(2);
subplot(611);
plot(x);
ylabel('x');
for i=1:5
    % 對分解的第5層到第1層的高頻系數近行重構
    d = wrcoef('d',c,l,'db5',6-i);
    subplot(6,1,i+1);
    plot(d);
    ylabel(['d',num2str(6-i)]);
end