clc;
clear;
% 當前擴展模式是補零模式 (詳見dwtmode函數)
% 構造原始一維尺度信號
randn('seed',531316785) 
s = 2 + kron(ones(1,8),[1 -1]) + ((1:16).^2)/32 + 0.2*randn(1,16);
% 進行單尺度離散haar小波變換
[ca1,cd1] = dwt(s,'haar'); 
figure(1)
subplot(311); plot(s); title('原始信號'); 
subplot(323); plot(ca1); title('haar低頻系數'); 
subplot(324); plot(cd1); title('haar高頻系數');

% 進行單尺度db2離散小波變換并觀察最后系數的邊緣效果
[ca2,cd2] = dwt(s,'db2');
subplot(325); plot(ca2); title('db2低頻系數'); 
subplot(326); plot(cd2); title('db2高頻系數');axis([0 10 0 2]);
% 對于給定的小波，計算兩個相關的分解濾波器，并直接使用該濾波器計算低頻和高頻系數
[Lo_D,Hi_D] = wfilters('haar','d'); 
[ca3,cd3] = dwt(s,Lo_D,Hi_D);
figure(2)
subplot(211); plot(ca3); title('haar分解濾波器低頻系數'); 
subplot(212); plot(cd3); title('haar分解濾波器高頻系數');