% 載入原始圖像
load sinsin;
% X 包含裝載的圖像
% 繪制原始圖像
subplot(2,2,1);
image(X);
colormap(map);
title('原始圖像');
% 使用sym4對X在第3層進行SWT分解
swc = swt2(X,3,'sym4');
% 第二種用法
[ca,chd,cvd,cdd] = swt2(X,3,'sym4');
% 從小波分解結構swc，重構s
a0 = iswt2(swc,'sym4');
% 繪制重構圖像1
subplot(2,2,3);
image(a0);
colormap(map);
title('重構圖像1');
% 第二種用法
a0 = iswt2(ca,chd,cvd,cdd,'sym4');
% 繪制重構圖像2
subplot(2,2,4);
image(a0);
colormap(map);
title('重構圖像2');
% 檢查重構的效果
err = max(max(abs(X-a0)))
errbis = max(max(abs(X-a0bis)))