DWT盲水印嵌入和提取源代码

数字水印程序，包括多分辨嵌入水印、提取水印程序，盲水印嵌入、提取程序，以及滤波攻击和JPEG攻击程序，PSNR和相似度计算公式（MATLAB）

%%%%%   盲水印嵌入算法   %%%%%
file_name=‘mark.bmp‘;
x=imread（file_name）;
p=im2bw（x0.5）;
figure（2）
imshow（p）
imwrite（p‘erzhituxiang.bmp‘‘bmp‘）
w=reshape（p164*64）;
file_name=‘lena.bmp‘;
Z=imread（file_name）;
Z_uint8=uint8（Z）;
imwrite（Z_uint8‘jiji.bmp‘‘bmp‘）;
figure（1）
imshow（Z_uint8[]）
s=double（imread（‘jiji.bmp‘））;
[cA1cH1cV1cD1]=dwt2（s‘haar‘）;
[cA2cH2cV2cD2]=dwt2（cA1‘haar‘）;

cd222=MZigZag（cA2）;
ff=zeros（14096）;
for i=1:4096
    m=round（cd222（i）.val/30）;
    t=cd222（i）.val-m*30;
 if cd222（i）.val>=0&w（i）==1
    if mod（m2）==0
        ff（i）=m*30;
    else ff（i）=m*30+30;
 end
 end
 if cd222（i）.val>=0&w（i）==0
    if mod（m2）==0
        ff（i）=m*30+30;
    else ff（i）=m*30; 
end
end
if cd222（i）.val<0&w（i）==1
    if mod（m2）==0
        ff（i）=m*30;
    else ff（i）=m*30+30;
end
end
if cd222（i）.val<0&w（i）==0
    if mod（m2）==0
        ff（i）=m*30+30;
    else ff（i）=m*30;
end
end
end
dd=invzigzag（ff）;
cax1=idwt2（ddcH2cV2cD2‘haar‘）;
water=idwt2（cax1cH1cV1cD1‘haar‘）;
water_uint8=uint8（water）;
imwrite（water_uint8‘ji1.bmp‘‘bmp‘）;
psnr=PSNR（Zwater）
cq=CQ（Zwater）
figure（3）
imshow（water[]）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-03-16 19:11  盲水印嵌入提取\
     文件        1172  2013-03-16 19:05  盲水印嵌入提取\blind_emb.m
     文件         626  2013-03-16 18:50  盲水印嵌入提取\blind_extraction.m
     文件         469  2013-03-16 18:46  盲水印嵌入提取\CQ.m
     文件         183  2013-03-16 18:49  盲水印嵌入提取\fil_code.m
     文件        1300  2013-03-16 18:47  盲水印嵌入提取\jpeg_extraction.m
     文件       66616  2007-07-04 16:54  盲水印嵌入提取\lena.bmp
     文件       12342  2011-07-09 17:23  盲水印嵌入提取\mark.bmp
     文件        1027  2013-03-16 19:11  盲水印嵌入提取\mul_emb.m
     文件        1294  2013-03-16 18:49  盲水印嵌入提取\mul_extration.m
     文件         453  2013-03-16 18:46  盲水印嵌入提取\PSNR.m