%装载并显示原始图像
load woman;
subplot(2,2,1);
image(X);colormap(map);
title('原始图像');
axis square;
disp('压缩前图像的大小');
whos('X')
%分解图像，提取分解结构中的第一层系数
[c,s]=wavedec2(X, 2, 'bior3.7');

ca1=appcoef2(c, s, 'bior3.7', 1);
ch1=detcoef2('h', c, s, 1);
cv1=detcoef2('v', c, s, 1);
cd1=detcoef2('d', c, s, 1);

a1=wrcoef2('a', c, s, 'bior3.7', 1);
h1=wrcoef2('h', c, s, 'bior3.7', 1);
v1=wrcoef2('v', c, s, 'bior3.7', 1);
d1= wrcoef2('d', c, s, 'bior3.7', 1);
c1=[a1, h1; v1, d1];

subplot(222);
colormap(gray);
image(c1);
axis square;
title ('分解后低频和高频信息');

ca1=appcoef2(c, s, 'bior3.7',1);
ca1=wcodemat(ca1, 440, 'mat', 0);
ca1=0.5*double(ca1);
subplot(223);
image(ca1);
colormap(gray);
axis square;
title('第一次压缩后图像(第一层低频信息)');
disp('第一次压缩图像的大小为：');
whos('ca1')


ca2=appcoef2(c, s, 'bior3.7', 2);
ca2=wcodemat(ca2, 440, 'mat', 0);
ca2=0.25*ca2;
subplot(224);
image(ca2);
colormap(gray);
axis square;
title ('第二次压缩后（第二层低频信息）');
disp('第二次压缩图像的大小为：');
whos('ca2')


在此程序的基础上根据论文中的理论将论文转为代码形式，论文在附件中有，希望matlab高手救急~酬劳方面可商议
有意者联系QQ：464423517
sxykerry@126.com