clc;
clear;
% 載入一維原始信號
load noisbloc; s = noisbloc;
% 畫出原始信號
subplot(4,1,1);plot(s);title('原始信號');
% 使用db1進行在第3層SWT分解.
[swa,swd] = swt(s,3,'db1');
% 畫出從1～3層離散小波變換的低頻系數和高頻系數
subplot(4,2,3);plot(swa(3,:));title('第3層低頻系數');
subplot(4,2,5);plot(swa(2,:));title('第2層低頻系數');
subplot(4,2,7);plot(swa(1,:));title('第1層低頻系數');
subplot(4,2,4);plot(swd(3,:));title('第3層高頻系數');
subplot(4,2,6);plot(swd(2,:));title('第2層高頻系數');
subplot(4,2,8);plot(swd(1,:));title('第1層高頻系數');