% 當前延拓模式是補零
% 裝載原始圖像
load sinsin;
% X 包含裝載的圖像
% 繪制原始圖像
subplot(1,2,1);
image(X);
colormap(map);
title('原始圖像');
sX = size(X);
% 使用db4對X進行單尺度分解
[cA1,cH1,cV1,cD1] = dwt2(X,'db4');
% 利用尺度為1時的系數直接重構
A0 = idwt2(cA1,cH1,cV1,cD1,'db4',sX);
% 繪制重構圖像
subplot(1,2,2);
image(A0);
title('重構圖像');
% 檢查重構效果
max(max(abs(X-A0)))