%裝載并顯示原始圖像
load belmont2;
subplot(1,2,1);
image(X);colormap(map);title('原始圖像');
%采用默認的全局閾值
[thr,sorh,keepapp,crit]=ddencmp('cmp','wp',X);
%圖像壓縮
[Xc,treed,perf0,perfl2]=wpdencmp(X,sorh,3,'bior3.1',crit,thr,keepapp);
subplot(1,2,2);image(Xc);colormap(map);title('壓縮后的圖像');
%給出壓縮效率
disp('小波分解系數中置0的系數個數百分比：');
perf0
disp('壓縮后圖像剩余能量百分比：');
perfl2