% 當前延拓模式是補零
% 裝載原始圖像
load sinsin; 
% X 包含原始圖像
% 繪制原始圖像
subplot(3,3,1);
image(X);
colormap(map);
title('原始圖像');

% 使用db1對X進行尺度為2的分解
[c,s] = wavedec2(X,2,'db1');
sizex = size(X)
sizec = size(c)
val_s = s 
% 對于各個方向，從小波分解結構[c,s]中提取尺度為2時的高頻系數
[chd2,cvd2,cdd2] = detcoef2('all',c,s,2); 
sizecd2 = size(chd2)
% 繪制尺度為2時高頻圖像
subplot(3,3,4);
image(chd2);
title('尺度為2時高頻圖像的水平部分');
subplot(3,3,5);
image(cvd2);
title('尺度為2時高頻圖像的垂直部分');
subplot(3,3,6);
image(cdd2);
title('尺度為2時高頻圖像的對角部分');
% 對于各個方向，從小波分解結構[c,s]中提取尺度為1時的高頻系數
[chd1,cvd1,cdd1] = detcoef2('all',c,s,1); 
sizecd1 = size(chd1)
% 繪制尺度為1時高頻圖像
subplot(3,3,7);
image(chd1);
title('尺度為1時高頻圖像的水平部分');
subplot(3,3,8);
image(cvd1);
title('尺度為1時高頻圖像的垂直部分');
subplot(3,3,9);
image(cdd1);
title('尺度為1時高頻圖像的對角部分');