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