load woman;
%打開原始圖像
imshow(X,map);
% X包含原始圖像信息. 
% map包含打開的色圖
nbcol = size(map,1); 
[cA1,cH1,cV1,cD1] = dwt2(X,'db1');
% 對圖像X執行單層分解，小波為db1
cod_X = wcodemat(X,nbcol); 
cod_cA1 = wcodemat(cA1,nbcol); 
cod_cH1 = wcodemat(cH1,nbcol); 
cod_cV1 = wcodemat(cV1,nbcol); 
cod_cD1 = wcodemat(cD1,nbcol); 
%圖像編碼
figure 
subplot(221)
imshow(cod_cA1,map)
title('近似細節系數')
subplot(222)
imshow(cod_cH1,map)
title('水平細節系數')
subplot(223)
imshow(cod_cV1,map)
title('垂直細節系數')
subplot(224)
imshow(cod_cD1,map)
title('對角細節系數')

Y = idwt2(cA1,cH1,cV1,cD1,'db1',nbcol);
% 對分解的細節系數執行單層重構，小波為db1
figure;
imshow(Y,map);
title('重構后的圖像')