load noispol;
s=noispol;
%==========================
%用db3小波進行4層分解
[c,l]=wavedec(s,4,'db3');
subplot(5,1,1);
plot(s);
title('原始信號及各層細節信號重構圖');
Ylabel('s');
%==========================
%重構分解結構[c，l]中的高頻部分
for i=1:4
    decmp=wrcoef('d',c,l,'db3',5-i);
    subplot(5,1,i+1);
    plot(decmp);
Ylabel(['d',num2str(5-i)]);
end