clc;?
M=8;

omark=double（imread（‘C:\Users\micheal\Desktop\信息理論課設\new?code&picture\3.tif‘））/255;
mo=size（omark1）;
no=size（omark2）;
%?以下生成水印信息?%
mark=wextend（‘2D‘‘sym‘omarkM）;?%///////////////////////////////////
[mono]=size（mark）;
alpha=30;
rand1=randn（18）;
rand2=randn（18）;
%%%%%%?讀取宿主圖像?%%%%%%%%
cimage=imread（‘C:\Users\micheal\Desktop\信息理論課設\new?code&picture\2.tif‘）;%////////////////////////////////F:\1\nwu.tif???F:\1\Cameraman.tif
figure（1）;
subplot（231）imshow（cimage[]）title（‘直接提取的圖像‘）;
[mcnc]=size（cimage）;
cda0=blkproc（cimage[88]‘dct2‘）;?
%%%%%%%?嵌入水印?%%%%%%%%%
cda0=wextend（‘2D‘‘sym‘cda0M）;%////////////////////////////////////////
cda1=cda0;??
for?i=M+1:mo-M?
for?j=M+1:no-M?
?????x=i;y=j;
????????if?mark（ij）==1
????????k=rand1;
????????else
????????k=rand2;?
????????end
????cda1（x+1y+8）=cda0（x+1y+8）+alpha*k（1）;
????cda1（x+2y+7）=cda0（x+2y+7）+alpha*k（2）;
????cda1（x+3y+6）=cda0（x+3y+6）+alpha*k（3）;
????cda1（x+4y+5）=cda0（x+4y+5）+alpha*k（4）;
????cda1（x+5y+4）=cda0（x+5y+4）+alpha*k（5）;
????cda1（x+6y+3）=cda0（x+6y+3）+alpha*k（6）;
????cda1（x+7y+2）=cda0（x+7y+2）+alpha*k（7）;
????cda1（x+8y+1）=cda0（x+8y+1）+alpha*k（8）;
????end
end
%%%%%%?嵌入水印后圖象?%%%%%
cda2=cda1（M+1:mo-MM+1:no-M）;
wimage=blkproc（cda2[88]‘idct2‘）;?
wimage_1=uint8（wimage）;
imwrite（wimage_1‘withmark.bmp‘‘bmp‘）;
subplot（232）imshow（wimage[]）title（‘嵌入水印后圖象‘）;
%%%%%%%%%%%%%%?進行攻擊測試?%%%%%%%%%%%
disp（‘1-->加入白噪聲‘）;
disp（‘2-->高斯低通濾波‘）;
disp（‘3-->剪切圖像‘）;
disp（‘4-->旋轉攻擊‘）;
disp（‘5-->直接檢測‘）
begin=input（‘請選擇攻擊（1-5）：‘）
switch?begin
%%%%%%%?加入白噪聲?%%%%%%%%
????case?1
Aimage1=wimage;
Wnoise=20*randn（size（Aimage1））;
Aimage1=Aimage1+Wnoise;
subplot（234）imshow（Aimage1[]）title（‘加入白噪聲后的圖象‘）;
att=Aimage1;
imwrite（att‘whitenoiseimage.bmp‘）;
%%%%%%%?高斯低通濾波?%%%%%%%
????case?2
Aimage2=wimage;
H=fspecial（‘gaussian‘[44]0.5）;
Aimage2=imfilter（Aimage2H）;
subplot（234）imshow（Aimage2[]）title（‘高斯低通濾波后的圖象‘）;
att=Aimage2;
imwrite（att‘gaussianimage.bmp‘）
%%%%%%%%?剪切攻擊?%%%%%%%%
????case?3
Aimage3=wimage;
Aimage3（1:1281:128）=256;
subplot（234）imshow（Aimage3[]）title（‘剪切后的圖象‘）;
att=Aimage3;
imwrite（att‘cutpartimage.bmp‘）;
%%%%%%%%?旋轉攻擊?%%%%%%%%
????case?4
Aimage4=wimage;
Aimage4=imrotate（Aimage410‘bilinear‘‘crop‘）;
Aimage_4=mat2gray（Aimage4）;
subplot（234）imshow（Aimage_4[]）title（‘旋轉10?度后的圖象‘）;
att=Aimage_4;
imwrite（att‘rotatedimage.bmp‘）;
%%%%%%%?沒有受到攻擊?%%%%%%%
????case?5
subplot（234）imshow（wimage[]）title（‘直接提取的圖像‘）;
att=wimage;
imwrite（att‘directimage.bmp‘）;
end
%%%%%%%%%%%%%?提取水印?%%%%%%%%%%%%%%
tmark_0=wimage;

tmark_0=wextend（‘2D‘‘sym‘tmark_0M）;%///////////////////////////////////
tmark_0=blkproc（tmark_0[88]‘dct2‘）;
tmark_1=tmark_0;%//////////////////////////////
pass=zeros（18）;
for?i=M+1:mo-M
????for?j=M+1:no-M
?????????x=i;y=j;
??????
????????pass（1）=tmark_0（x+1y+8）;
????????pass（2）=tmark_0（x+2y+7）;
????????pass（3）=tmark_0（x+3y+6）;
????????pa

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????6885600??2013-07-02?15:53??數字水印算法.docx

?????文件??????33951??2013-07-01?14:13??new?code&picture\2.jpg

?????文件??????66024??2013-07-02?00:24??new?code&picture\2.tif

?????文件???????8208??2013-07-02?00:23??new?code&picture\3.jpg

?????文件??????14632??2013-07-02?00:24??new?code&picture\3.tif

?????文件??????66614??2013-07-02?13:38??new?code&picture\cutpartimage.bmp

?????文件??????66614??2013-07-02?13:35??new?code&picture\gaussianimage.bmp

?????文件?????130130??2013-07-02?13:43??new?code&picture\QQ截圖20130702134307.jpg

?????文件??????66614??2013-07-02?13:42??new?code&picture\rotatedimage.bmp

?????文件???????6231??2013-07-02?00:28??new?code&picture\Untitled2.m

?????文件???????4321??2013-07-02?13:32??new?code&picture\Untitled8.m

?????文件???????1403??2013-07-01?22:58??new?code&picture\Untitled9.m

?????文件??????66614??2013-07-02?13:32??new?code&picture\whitenoiseimage.bmp

?????文件??????66614??2013-07-02?13:42??new?code&picture\withmark.bmp

?????文件?????132291??2013-07-02?13:39??new?code&picture\剪切.jpg

?????文件??????74632??2013-07-01?22:54??new?code&picture\新建文件夾\airplane.jpg

?????文件??????66614??2013-07-02?01:03??new?code&picture\新建文件夾\cutpartimage.bmp

?????文件??????66614??2013-07-02?01:04??new?code&picture\新建文件夾\directimage.bmp

?????文件??????66614??2013-07-02?01:02??new?code&picture\新建文件夾\gaussianimage.bmp

?????文件??????67471??2013-07-01?22:40??new?code&picture\新建文件夾\linux.jpg

?????文件??????18675??2013-07-01?22:48??new?code&picture\新建文件夾\linux2.jpg

?????文件???????5755??2013-07-01?21:44??new?code&picture\新建文件夾\mark2.jpg

?????文件??????66614??2013-07-02?01:04??new?code&picture\新建文件夾\rotatedimage.bmp

?????文件??????26659??2013-07-01?22:33??new?code&picture\新建文件夾\sea.jpg

?????文件???????8208??2013-07-01?22:34??new?code&picture\新建文件夾\simplepanda.jpg

?????文件??????66614??2013-07-02?00:28??new?code&picture\新建文件夾\whitenoiseimage.bmp

?????文件??????66614??2013-07-02?01:04??new?code&picture\新建文件夾\withmark.bmp

?????文件??????18375??2013-07-01?22:42??new?code&picture\新建文件夾\write.jpg

?????文件??????38868??2013-07-02?01:03??new?code&picture\新建文件夾\剪切圖像.jpg

?????文件??????39382??2013-07-02?01:02??new?code&picture\新建文件夾\加入白噪聲.jpg

............此處省略11個文件信息