% 當前延拓模式是補零
% 裝載原始圖像
load sinsin; 
% X 包含裝載的圖像
subplot(3,2,1);
image(X);
colormap(map);
title('原始圖像');
% 使用db4對X進行尺度為2的分解
[c,s] = wavedec2(X,2,'db4');
% 由系數重構第1層低頻和細節部分
% 也可以使用wrcoef2完成，或者使用下面的步驟等價實現
% Step 1: 從分解結構[c,s]中提取系數
% Step 2: 使用upcoef2進行重構
siz = s(size(s,1),:); 
ca1 = appcoef2(c,s,'db4',1); 
a1 = upcoef2('a',ca1,'db4',1,siz);
% 繪制低頻重構圖
subplot(3,2,3);
image(a1);
colormap(map);
title('低頻重構圖');
chd1 = detcoef2('h',c,s,1); 
hd1 = upcoef2('h',chd1,'db4',1,siz); 
% 繪制高頻水平部分重構圖
subplot(3,2,4);
image(hd1);
colormap(map);
title('高頻水平部分重構圖');
cvd1 = detcoef2('v',c,s,1); 
vd1 = upcoef2('v',cvd1,'db4',1,siz);
% 繪制高頻垂直部分重構圖
subplot(3,2,5);
image(vd1);
colormap(map);
title('高頻垂直部分重構圖');
cdd1 = detcoef2('d',c,s,1); 
dd1 = upcoef2('d',cdd1,'db4',1,siz);
% 繪制高頻對角部分重構圖
subplot(3,2,6);
image(dd1);
colormap(map);
title('高頻對角部分重構圖');