%裝載并顯示原始圖像
load flujet;
subplot(1,2,1);
image(X);
colormap(map);
title('原始圖像');

%首先利用db3小波對(duì)圖像X進(jìn)行2層分解
[c,l]=wavedec2(X,2,'db3');

%全局閾值
[thr,sorh,keepapp]=ddencmp('cmp','wv',X);

%壓縮處理：對(duì)所有高頻系數(shù)進(jìn)行同樣的閾值量化處理
[Xcmp,cxc,lxc,perf0,perfl2]=wdencmp('gbl',c,l,'db3',2,thr,sorh,keepapp);

%將壓縮后的圖像于原始圖像相比較
subplot(1,2,2);
image(Xcmp);colormap(map);
title('壓縮后的圖像');

%顯示相關(guān)參數(shù)
disp('小波分解系數(shù)中為0的系數(shù)個(gè)數(shù)百分比：');
perf0
disp('壓縮后保留能量百分比：');
perfl2