% 當前延拓模式是補零
% 裝載圖像
load sinsin;
% X包含載入的圖像
% 繪制原始圖像
figure(1);
subplot(2,2,1);
image(X);
colormap(map);
title('原始圖像');
% 使用sym5對X進行尺度為2的分解
[c,s] = wavedec2(X,2,'sym5');
% 從小波分解結構[c,s]進行尺度為1和2時的低頻重構
a1 = wrcoef2('a',c,s,'sym5',1); 
a2 = wrcoef2('a',c,s,'sym5',2);
% 繪制尺度為1時的低頻圖像
subplot(2,2,3);
image(a1);colormap(map);
title('尺度為1時的低頻圖像');
% 繪制尺度為2時的低頻圖像
subplot(2,2,4);
image(a2);colormap(map);
title('尺度為2時的低頻圖像');

% 從小波分解結構[c,s]在尺度為2時重構高頻
% 'h' 是水平方向
% 'v' 是垂直方向
% 'd' 是對角方向
hd2 = wrcoef2('h',c,s,'sym5',2); 
vd2 = wrcoef2('v',c,s,'sym5',2); 
dd2 = wrcoef2('d',c,s,'sym5',2);
% 繪制高頻圖像
figure(2);
subplot(2,2,1);
image(hd2);colormap(map);
title('尺度為2時的水平高頻圖像');
subplot(2,2,2);
image(vd2);colormap(map);
title('尺度為2時的垂直高頻圖像');
subplot(2,2,3);
image(dd2);colormap(map);
title('尺度為2時的對角高頻圖像');
% 驗證這些圖像的長度都是sX
sX = size(X)
sa1 = size(a1)
shd2 = size(hd2)