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