小波变换的源程序DWT

基于小波变换的源程序，简单易懂,实用性强

function t1
%数字水印的嵌入与提取
A=imread（‘1.tif‘）;
B=imread（‘2.tif‘）;
A=rgb2gray（A）;
B=rgb2gray（B）;


K=64;

for i=1:4
     for j=1:4
         x=（i-1）*K+1;
         y=（j-1）*K+1;
         block=A（x:x+K-1y:y+K-1）;
         [ca1ch1cv1cd1]=dwt2（block‘db1‘）;
         [ca2ch2cv2cd2]=dwt2（B‘db1‘）;
         ca3=ca1+（ca2+rand（32））/1000;
         p=idwt2（ca3ch1cv1cd1‘db1‘size（block））; 
         for q=1:64;
             for o=1:64
                  M（（i-1）*64+q（j-1）*64+o）=p（qo）;
              end
          end
     end
end 
     

      
      

 %display original image and watermark
figure（1）
imshow（A）;
title（‘原始图象‘）;

%display  watermark
figure（2）
imshow（B）;
title（‘需嵌入的水印‘）;

% display watermarked image
figure（3）
imshow（uint8（M））;
title（‘嵌入水印后的图象‘）;


%数字水印的提取
ca4=1000*（ca3-ca1）-rand（32）;
figure（4）;
k=idwt2（ca4ch2cv2cd2‘db1‘size（B））;
imshow（uint8（k））;
title（‘提取出来的数字水印图象‘）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-05-24 18:02  基于dwt的数字水印源代码\
     文件      208036  2006-05-21 13:17  基于dwt的数字水印源代码\1.tif
     文件       20044  2006-05-18 14:58  基于dwt的数字水印源代码\2.tif
     文件      196662  2005-05-12 00:53  基于dwt的数字水印源代码\pic.bmp
     文件         936  2006-05-24 12:45  基于dwt的数字水印源代码\t1.m
     文件       20044  2006-05-18 14:58  基于dwt的数字水印源代码\watermark.bmp
     文件        2102  2005-05-28 15:06  基于dwt的数字水印源代码\watermark1.bmp