【正文】 o=abs_block(ii,jj)。abs_block(ii,jj)=abs_block(ii,jj)*(1+k*pn_sequence_zero(ll))。 abs_block(blocksizeii+1,blocksizejj+1)=abs_block(blocksizeii+1…,blocksizejj+1)+abs_block(ii,jj)abs_block_o。 ll=ll+1。 end end end else for ii=1:blocksize for jj=1:blocksize if (filter_m(ii,jj)==1) abs_block_o=abs_block(ii,jj)。abs_block(ii,jj)=abs_block(ii,jj)*(1+k*pn_sequence_one(ll))。abs_block(blocksizeii+1,blocksizejj+1)=abs_block(blocksizeii+1…,blocksizejj+1)…+abs_block(ii,jj)abs_block_o。 ll=ll+1。 end end end end % 進(jìn)行傅立葉逆變換abs_block=fftshift(abs_block)。 watermarked_image(y:y+blocksize1,x:x+blocksize1)=…abs(ifft2(abs_block.*exp(i*angle_block)))。 % 移動到下一塊 if (x+blocksize)=Nc x=1。 y=y+blocksize。 else x=x+blocksize。 end waitbar(kk/length(message_vector),h)。endclose(h)。% 轉(zhuǎn)換為uint8，并寫入watermarked_image_int=uint8(watermarked_image)。imwrite(watermarked_image_int,39。39。,39。bmp39。)。% 計算運行時間elapsed_time=cputimestart_time,% 計算psnrpsnr=psnr(cover_object,watermarked_image),% 顯示水印，嵌入水印圖像與原始圖像figure(1)subplot(1,2,1)。imshow(message,[])。title(39。水印39。)。subplot(1,2,2)imshow(outImg,[])。title(39。置亂水印39。)。figure(2)subplot(1,2,1)。imshow(watermarked_image,[])name=39。嵌入水印圖像39。title(strcat(num2str(name),39。k=39。,num2str(k),39。PSNR=39。,num2str(psnr)))。subplot(1,2,2)。imshow(cover_object,[])。title(39。原始圖像39。)。PSNR(峰值信噪比)函數(shù)PSNR(峰值信噪比)衡量嵌入水印后圖像的失真程度%計算PSNR(峰值信噪比)function dPSNR=psnr(ImageA,ImageB)if (size(ImageA,1)~=size(ImageB,1))or(size(ImageA,2)~=size(ImageB,2)) error(39。ImageAImageB39。)。 dPSNR=0。 return。endImageA=double(ImageA)。ImageB=double(ImageB)。M=size(ImageA,1)。N=size(ImageA,2)。d=0 。for i=1:M for j=1:N d=d+(ImageA(i,j)ImageB(i,j)).^2。 endenddPSNR=10*log10((M*N*max(max(ImageA.^2)))/d)。returnNC(歸一化相關(guān)系數(shù))計算NC(歸一化相關(guān)系數(shù))衡量提取水印與原始水印的相似程度%計算nc(歸一化相關(guān)系數(shù))function dNC=nc(ImageA,ImageB)if (size(ImageA,1)~=size(ImageB,1))or(size(ImageA,2)~=size(ImageB,2)) error(39。ImageAImageB39。)。 dNC=0。 return。endImageA=double(ImageA)。ImageB=double(ImageB)。M=size(ImageA,1)。N=size(ImageA,2)。d1=0 。d2=0。d3=0。for i=1:M for j=1:N d1=d1+ImageA(i,j)*ImageB(i,j)。 d2=d2+ImageA(i,j)*ImageA(i,j)。 d3=d3+ImageB(i,j)*ImageB(i,j)。 endenddNC=d1/(sqrt(d2)*sqrt(d3))。return水印提取代碼%提取源碼clcclear all。%保存開始時間start_time=cputime。iTimes=20。 %置亂次數(shù)blocksize=8。 % 設(shè)置塊的大小filter_m=[ 1,1,1,1,1,1,1,1。 % 濾波矩陣 1,1,1,1,1,1,1,1。 1,1,0,0,0,0,1,1。 1,1,0,0,0,0,1,1。 0,0,0,0,0,0,0,0。 0,0,0,0,0,0,0,0。 0,0,0,0,0,0,0,0。 0,0,0,0,0,0,0,0 ]。% 讀入嵌入水印圖像file_name=39。39。watermarked_image=double(imread(file_name))。% 嵌入水印圖像矩陣的行數(shù)與列數(shù)Mw=size(watermarked_image,1)。 Nw=size(watermarked_image,2)。 % 最大可嵌入信息量max_message=Mw*Nw/(blocksize^2)。% 讀入原始水印file_name=39。39。orig_watermark=double(imread(file_name))。% 原始水印矩陣的行數(shù)與列數(shù)Mo=size(orig_watermark,1)。 No=size(orig_watermark,2)。 %置隨機(jī)數(shù)發(fā)生器的狀態(tài)為1100key=1100。rand(39。state39。,key)。% 產(chǎn)生偽隨機(jī)序列pn_sequence_zero=round(2*(rand(1,sum(sum(filter_m)))))。pn_sequence_one=round(2*(rand(1,sum(sum(filter_m)))))。 % 將圖像分塊x=1。 y=1。h=waitbar(0,39。提取水印，請等待39。)。for (kk = 1:max_message) % 傅立葉變換fft_block_w=fft2(watermarked_image(y:y+blocksize1,x:x+blocksize1))。 abs_block_w=abs(fftshift(fft_block_w))。 ll=1。 for ii=1:blocksize for jj=1:blocksize if (filter_m(ii,jj)==1) sequence(ll)=abs_block_w(ii,jj)。 ll=ll+1。 end end end % 計算sequence與pn_sequence_zero和pn_sequence_one的相關(guān)系數(shù) correlation_zero(kk)=corr2(pn_sequence_zero,sequence)。 correlation_one(kk)=corr2(pn_sequence_one,sequence)。 % 移動到下一塊 if (x+blocksize) = Nw x=1。 y=y+blocksize。 else x=x+blocksize。 end waitbar(kk/max_message,h)。endclose(h)。%如果correlation_zerocorrelation_one，那么message_vector=0，反之為1for (kk=1:Mo*No) if correlation_zero(kk)correlation_one(kk) %if (correlation(kk) mean(correlation(1:Mo*No))) message_vector(kk)=0。 else message_vector(kk)=1。 endend% Arnold置亂tempImg=reshape(message_vector(1:Mo*No),Mo,No)。message_arnold=tempImg。for n=1:iTimes % 次數(shù) for u=1:Mo for v=1:No temp=tempImg(u,v)。 ax=mod((u1)+(v1),Mo)+1。 ay=mod((u1)+2*(v1),No)+1。 outImg(ax,ay)=temp。 end end tempImg=outImg。endmessage=outImg。% 計算運行時間elapsed_time=cputimestart_time,%計算NC（歸一化相關(guān)系數(shù)）NC=nc(message,orig_watermark)% 顯示提取水印與原始水印figure(3)subplot(1,2,1)。imshow(message,[])。name=39。提取水印39。title(strcat(num2str(name),39。 NC=39。,num2str(NC)))。subplot(1,2,2)。imshow(orig_watermark,[])title(39。原始水印39。)。附錄C 離散余弦變換(DCT)的數(shù)字水印嵌入提取源碼clcs=512。 %原圖像長度N=32。 %水印圖像長度K=8。 D=zeros(s)。 %zeros(n)產(chǎn)生nn的全0方陣E=。I=imread(39。39。)。subplot(2,2,1)。imshow(I)。title(39。原始圖像39。)。I=double(I)/512。 %把I轉(zhuǎn)換成double型%讀入水印J=imread(39。39。)。subplot(2,2,2)。imshow(J)。title(39。水印圖像39。)。%對水印縮放J=double(imresize(J,[64,64]))。subplot(2,2,3)。imshow(J)。title(39?？s放水印圖像39。)。%嵌入水印for p=1:384/K for q=1:384/K x=(p1)*K+1。 y=(q1)*K+1。 I_dct=I(x:x+K1,y:y+K1)。 I_dct1=dct2(I_dct)。 if J(p,q)==0 alfa=1。 else alfa=1。 end I_dct2=I_dct1+alfa*E。 I_dct=idct2(I_dct2)。 D(x:x+K1,y:y+K1)=I_dct。 end