load noisbloc
s= noisbloc;
ls=length(s);
plot(s);
[swa,swd] = swt(s,1,'db4');
subplot(2,1,1);plot(swa);title('低頻')
subplot(2,1,2);plot(swd);title('高頻')
A0 = iswt(swa,swd,'db4');
err=norm(s-A0)
figure(2)
subplot(3,1,1);plot(s);title('原始信號')
subplot(3,1,2);plot(A0);title('重構信號')
subplot(3,1,3);plot(s-A0);title('誤差信號')
nulcfs = zeros(size(swa));
A1 = iswt(swa,nulcfs,'db4');
D1 = iswt(nulcfs,swd,'db4');
figure(3)
subplot(2,1,1), plot(A1); 
title('重構低頻A1');
subplot(2,1,2), plot(D1); 
title('重構高頻D1');
[swa,swd] = swt(s,3,'db4');
kp = 0;
for i = 1:3
subplot(3,2,kp+1), plot(swa(i,:));
title(['低頻系數：level ',num2str(i)])
subplot(3,2,kp+2), plot(swd(i,:));
title(['高頻系數：level ',num2str(i)])
kp = kp + 2;
end
mzero = zeros(size(swd));
A = mzero;
A(3,:) = iswt(swa(3,:),mzero,'db4');  %  重構低頻分量
D = iswt(mzero,swd,'db4');  %  重構所有高頻分量
s_r=A(3,:)+D;  %  重構源信號
norm(s_r-s)    %  誤差
figure(4)
subplot(2,2,1);
plot(s);        %  源信號
title('源信號');
subplot(2,2,2)
plot(A(3,:));   %  信號的低頻部分
title('信號的低頻部分');
subplot(2,2,3)
plot(D);        %  信號的高頻部分
title('信號的高頻部分');
subplot(2,2,4)
plot(s_r);      %  重構信號
title('重構信號');