% 裝載原始圖像
load belmont2; 
nbc = size(map,1);
% 使用sym4小波對圖像進行3層分解
wname = 'sym4'; lev = 3;
[c,s] = wavedec2(X,lev,wname);
% 圖像壓縮，使用推薦參數選擇各層不同的閾值
alpha = 1.5; m = 2.7*prod(s(1,:));
[thr,nkeep] = wdcbm2(c,s,alpha,m)
% 使用上面的閾值和硬閾值處理方式進行圖像壓縮
[xd,cxd,sxd,perf0,perfl2] =wdencmp('lvd',c,s,wname,lev,thr,'h');
% 畫出原始圖像和壓縮后的圖像
colormap(pink(nbc));
subplot(121), image(wcodemat(X,nbc)),
title('原始圖像')
subplot(122), image(wcodemat(xd,nbc)),
title('壓縮后的圖像')
disp('恢復百分比: ');
num2str(perfl2)
disp('壓縮百分比: ');
num2str(perf0)