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