【正文】 =omark。alpha=200。R1=randn(1,8)。rand1=round(R1)。R2=randn(1,8)。rand2=round(R2)。%%%%%%% 讀取宿主圖像 %%%%%%%cimage=imread(39。39。)。figure(1)。subplot(2,3,1),imshow(cimage,[]),title(39。原始的宿主圖像39。)。[mc,nc]=size(cimage)。cimage=double(cimage)。fft0=blkproc(cimage,[8,8],39。fft239。)。%%%%%%% 嵌入水印 %%%%%%%%%fft_abs=abs(fft0)。fft_angle=angle(fft0)。fft1=fft_abs。for ii=1:mo for jj=1:no x=(ii1)*8。y=(jj1)*8。 if mark(ii,jj)==1 k=rand1。 else k=rand2。 end fft1(x+1,y+8)=fft_abs(x+1,y+8)+alpha*k(1)。 fft1(x+2,y+7)=fft_abs(x+2,y+7)+alpha*k(2)。 fft1(x+3,y+6)=fft_abs(x+3,y+6)+alpha*k(3)。 fft1(x+4,y+5)=fft_abs(x+4,y+5)+alpha*k(4)。 fft1(x+5,y+4)=fft_abs(x+5,y+4)+alpha*k(5)。 fft1(x+6,y+3)=fft_abs(x+6,y+3)+alpha*k(6)。 fft1(x+7,y+2)=fft_abs(x+7,y+2)+alpha*k(7)。 fft1(x+8,y+1)=fft_abs(x+8,y+1)+alpha*k(8)。 endendfft_abs_0=fft1。fft_abs_1=(fft_abs_0.*exp(i*fft_angle))。wimage_0=blkproc(fft_abs_1,[8,8],39。ifft239。)。wimage=abs(wimage_0)。%%%%%% 顯示嵌入水印圖象 %%%%%wimage_1=uint8(wimage)。imwrite(wimage_1,39。39。,39。bmp39。)。subplot(2,3,2),imshow(wimage,[]),title(39。嵌入水印后的圖象39。)。%%%%%%%%%%%%% 進行各種攻擊 %%%%%%%%%%%%disp(39。1加入白噪聲39。)。disp(39。2高斯低通濾波39。)。disp(39。3剪切圖像39。)。disp(39。4旋轉(zhuǎn)圖象39。)。disp(39。5直接檢測39。)。begin=input(39。請選擇攻擊（15）39。)。switch begin%%%%%%%%加入白噪聲%%%%%%%% case 1Aimage1=wimage。Wnoise=12*randn(size(Aimage1))。Aimage1=Aimage1+Wnoise。subplot(2,3,4),imshow(Aimage1,[]),title(39。加入白噪聲后的圖象39。)。att=Aimage1。imwrite(att,39。39。)。%%%%%% 高斯低通濾波 %%%%%%%% case 2Aimage2=wimage。H=fspecial(39。gaussian39。,[4,4],)。Aimage2=imfilter(Aimage2,H)。subplot(2,3,4),imshow(Aimage2,[]),title(39。高斯低通濾波后的圖象39。)。att=Aimage2。imwrite(att,39。39。)%%%%%%% 剪切攻擊 %%%%%%%%% case 3Aimage3=wimage。Aimage3(1:128,1:128)=256。subplot(2,3,4),imshow(Aimage3,[]),title(39。剪切后的圖象39。)。att=Aimage3。imwrite(att,39。39。)。%%%%%%% 旋轉(zhuǎn)攻擊 %%%%%%%%% case 4Aimage4=wimage。Aimage4=imrotate(Aimage4,10,39。bilinear39。,39。crop39。)。Aimage_4=mat2gray(Aimage4)。subplot(2,3,4),imshow(Aimage_4,[]),title(39。旋轉(zhuǎn)10度后的圖象39。)。att=Aimage_4。imwrite(att,39。39。)。%%%%%%% 沒有受到攻擊 %%%%%%% case 5subplot(2,3,4),imshow(wimage,[]),title(39。直接提取的圖象39。)。att=wimage。imwrite(att,39。39。)。end%%%%%%%%%%%%%%% 提取水印 %%%%%%%%%%%%%tmark_0=att。tmark_0=blkproc(tmark_0,[8,8],39。fft239。)。tmark_0=abs(tmark_0)。pass=zeros(1,8)。for ii=1:mo for jj=1:no x=(ii1)*8。y=(jj1)*8。 pass(1)=tmark_0(x+1,y+8)。 pass(2)=tmark_0(x+2,y+7)。 pass(3)=tmark_0(x+3,y+6)。 pass(4)=tmark_0(x+4,y+5)。 pass(5)=tmark_0(x+5,y+4)。 pass(6)=tmark_0(x+6,y+3)。 pass(7)=tmark_0(x+7,y+2)。 pass(8)=tmark_0(x+8,y+1)。 if (corr2(pass,rand1)corr2(pass,rand2)) tmark_1(ii,jj)=1。 else tmark_1(ii,jj)=0。 end endend%%%%%%%%%% 計算NC（歸一化相關(guān)系數(shù)）%%%%%%%%g_mark=double(tmark_1)。o_mark=double(omark)。[m,n]=size(g_mark)。nc_0=0。nc_1=0。nc_2=0。for ii=1:m for jj=1:n nc_0=nc_0+g_mark(ii,jj)*o_mark(ii,jj)。 nc_1=nc_1+o_mark(ii,jj)*o_mark(ii,jj)。 nc_2=nc_2+g_mark(ii,jj)*g_mark(ii,jj)。 endendNC=nc_0/sqrt(nc_1*nc_2)。%%%%%% 顯示提取水印圖象 %%%%%subplot(2,3,5),imshow(39。39。),title(39。原始水印圖象39。)。subplot(2,3,6),imshow(tmark_1,[])。name=39。提取的水印圖像39。title(strcat(num2str(name),39。NC=39。,num2str(NC)))。致 謝經(jīng)過了四年的努力，我終于完成了我的學(xué)業(yè)，并順利完成了學(xué)士畢業(yè)論文，在這期間得到了諸多老師的指導(dǎo)與幫助，還有許多親人朋友的關(guān)心與支持，在這里一并表示衷心的感謝！首先，感謝竇德召老師，他在我寫畢業(yè)論文期間給我提供了許多資料，并指導(dǎo)我寫完了畢業(yè)論文，老師嚴謹樸實的治學(xué)態(tài)度，誨人不倦的教學(xué)精神，和藹可親的處世之道給我留下了很深刻的映像，讓我受益無窮。另外，還要感謝圖書館電子資源管理室，它在我查閱資料方面也給我提供了許多幫助！感謝所有關(guān)心過我的人，謹以我的努力作為回報！