%裝載并顯示原始圖像
load belmont2;
subplot(1,2,1);image(X);colormap(map);title('原始圖像');
%首先利用db3小波對圖像X進行2層分解
[c,l]=wavedec2(X,2,'db3');
%全局閾值
[thr,sorh,keepapp]=ddencmp('cmp','wv',X);
%壓縮處理：對所有高頻系數進行同樣的閾值量化處理
[Xcmp,cxc,lxc,perf0,perfl2]=wdencmp('gbl',c,l,'db3',2,thr,sorh,keepapp);
%將壓縮后的圖像與原始圖像相比較
subplot(1,2,2);
image(Xcmp);colormap(map);
title('壓縮后的圖像');
%顯示相關參數
disp('小波分解系數中為0的系數個數百分比：');
perf0
disp('壓縮后保留能量百分比：');
perfl2