二维小波变换，并显示各个子带图

clc;
clear all;
o=imread('E:\工作\科研\奇异值\image\lena.bmp');
% figure,imshow(o);
% o=[ 1 2;3 4;5 6];
o=double(o);
%进行第一层提升小波变换
% els ={'p',[-0.125 0.125],0};
% lshaarInt = liftwave('haar','int2int');
% lsnewInt = addlift(lshaarInt,els);
% [oa1,oh1,ov1,od1]=lwt2(o,lsnewInt);
% subplot(2,2,1);imshow(oa1,[])%;colormap(gray); 
% subplot(2,2,2);imshow(oh1,[]);colormap(gray); 
% subplot(2,2,3);imshow(ov1,[]);colormap(gray); 
% subplot(2,2,4);imshow(od1,[]);colormap(gray); 
map=gray;
[c,s] = wavedec2(o,3,'db1');
subplot(3,3,1);
% image(o);
imshow(o,[]);
colormap(map);
title('原始图像');
ca2=appcoef2(c,s,'db1',2);
subplot(3,3,2);
% image(ca2);
imshow(ca2,[]);
colormap(map);
title('尺度为2时低频图像');
ca1=appcoef2(c,s,'db1',1);
subplot(3,3,3);
% image(ca1);
imshow(ca1,[]);
colormap(map);
title('尺度为1时低频图像');
[chd2,cvd2,cdd2]=detcoef2('all',c,s,2);
subplot(3,3,4);
% image(chd2);
imshow(chd2,[]);
% a='尺度为2时';'高频图像的水平部分';
title(sprintf('%s \n %s','尺度为2时','高频图像的水平部分'));
subplot(3,3,5);
% image(cvd2);
imshow(cvd2,[]);
% title('尺度为2时的高频图像的垂直部分');
title(sprintf('%s \n %s','尺度为2时','高频图像的垂直部分'));
subplot(3,3,6);
% image(cdd2);
imshow(cdd2,[]);
% title('尺度为2时的高频图像的对角图像');
title(sprintf('%s \n %s','尺度为2时','高频图像的对角部分'));
[chd1,cvd1,cdd1]=detcoef2('all',c,s,1);
subplot(3,3,7);
% image(chd1);
imshow(chd1,[]);
% title('尺度为1时高频图像的水平部分');
title(strvcat('尺度为1时','高频图像的水平部分'));
subplot(3,3,8);
% image(cvd1);
imshow(cvd1,[]);
% title('尺度为1时高频图像的垂直部分');
title(strvcat('尺度为1时','高频图像的垂直部分'));
subplot(3,3,9);
% image(cdd1);
imshow(cdd1,[]);
% title('尺度为1时高频图像的对角部分');
title(strvcat('尺度为1时','高频图像的对角部分'));