program_13_18.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

%装载并显示原始图像
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