%裝載并顯示原始圖像
load flujet;
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