% 當前延拓模式是補零
% 裝載原始圖像
load sinsin; 
% 繪制原始圖像
subplot(2,2,1);
image(X);
colormap(map);
title('原始圖像');
% X 包含裝載的圖像
% 使用db1對X進行尺度為2的分解
[c,s] = wavedec2(X,2,'db1');
sizex = size(X)
sizec = size(c)
val_s = s
% 提取尺度為2時的低頻系數
ca2 = appcoef2(c,s,'db1',2); 
sizeca2 = size(ca2)
% 繪制尺度為2時的低頻圖像
subplot(2,2,2);
image(ca2);
colormap(map);
title('尺度為2時的圖像');
% 計算尺度為1時的低頻系數
ca1 = appcoef2(c,s,'db1',1); 
sizeca1 = size(ca1)
% 繪制尺度為1時的低頻圖像
subplot(2,2,3);
image(ca1);
colormap(map);
title('尺度為1時的圖像');