【正文】 DCT矩陣的區(qū)域1做奇異值分解并求奇異值的平均值 str=num2str(V(k))。 %判斷最高位數(shù)的奇偶性 x=str2double(str(1))。 if xfix(x/2)*2 == 1 %生成特征水印 W(k)=1。 else W(k)=0。 end。 if Bj{k}(x2,y2)=Bj{k}(x1,y1) %提取特征水印 W1(k)=1。 else W1(k)=0。 end。 k=k+1。 end。end。L=hundun(,2,300,k1)。 %生成混沌序列W2=xor(L,W1)。 %異或運(yùn)算解密特征水印NC=sum(W.*W2)/(sqrt(sum(W.^2))*sqrt(sum(W2.^2)))。 %求歸一化相關(guān)度subplot(2,2,3)。imshow(W)subplot(2,2,4)。imshow(W2)k=0。h=1。for i=1:n:S(1) for j=1:n:S(2) if W(h)~=W2(h) B(i:i+n1,j:j+n1)=0。 k=k+1。 end。 h=h+1。 end。end。 TAF=k/(h1)。if TAF= fprintf 未篡改！\nelse if TAF=er fprintf 惡意篡改！\n else fprintf 非惡意篡改！\n end。end。subplot(2,2,2)。imshow(B)。end function L=hundun(y0,v,she,number)y(1)=1v*y0^2。for i=2:she+number y(i)=1v*y(i1)^2。end。y=y(she+1:she+number)。L=fix((sign(y)+1)/2)。 %二值化end3. 一種新的小波域盲水印算法水印置亂：function N=mix(xo,uo,f,n) %圖像置亂函數(shù)，參數(shù)為（置亂的key，原始圖像矩陣，是否灰度化圖像）if n==1 R = f(:,:,1)。 G = f(:,:,2)。 B = f(:,:,3)。 f = .*R+.*G+.*B。 %灰度化圖像endimage_size=size(f)。L=image_size(1)*image_size(2)。for i=1:L %迭代生成混沌序列 x(i)=uo*xo*(1xo)。 xo=x(i)。endk=1。for i=1:image_size(1) %按列掃描水印圖像 for j=1:image_size(2) N(k)=f(i,j)。 k=k+1。 endend[x,N]=quicksort(x,N,1,L)。 function [A,N]=quicksort(A,N,low,high) %快速排序函數(shù)，增序if low high [A,N,w]=splition2(A,N,low,high)。 [A,N]=quicksort(A,N,low,w1)。 [A,N]=quicksort(A,N,w+1,high)。end function [A,N,w]=splition2(A,N,low,high)i=low。x=A(low)。for j=low+1:high if A(j)=x i=i+1。 if i~=j。 temp=A(i)。 A(i)=A(j)。 A(j)=temp。 temp=N(i)。 N(i)=N(j)。 N(j)=temp。 clear temp。 end endendtemp=A(i)。A(i)=A(low)。A(low)=temp。temp=N(i)。N(i)=N(low)。N(low)=temp。clear temp。w=i。水印反置亂：function W = imix(xo,uo,f) %置亂逆變換函數(shù)，參數(shù)為（置亂的key，置亂后的水印圖像矩陣）image_size=size(f)。L=image_size(1)*image_size(2)。for i=1:L %迭代生成混沌序列 x(i)=uo*xo*(1xo)。 xo=x(i)。endN=1:L。[x,N]=quicksort(x,N,1,L)。for i=1:L f1(N(i))=f(i)。endk=1。for i=1:64 %按列掃描水印圖像 for j=1:64 W(i,j)=f1(k)。 k=k+1。 endend function [A,N]=quicksort(A,N,low,high) %快速排序函數(shù)，增序if low high [A,N,w]=splition2(A,N,low,high)。 [A,N]=quicksort(A,N,low,w1)。 [A,N]=quicksort(A,N,w+1,high)。end function [A,N,w]=splition2(A,N,low,high)i=low。x=A(low)。for j=low+1:high if A(j)=x i=i+1。 if i~=j。 temp=A(i)。 A(i)=A(j)。 A(j)=temp。 temp=N(i)。 N(i)=N(j)。 N(j)=temp。 clear temp。 end endendtemp=A(i)。A(i)=A(low)。A(low)=temp。temp=N(i)。N(i)=N(low)。N(low)=temp。clear temp。w=i。嵌入水?。篺unction [f,C] = sy() %水印函數(shù)q=15。 %水印圖像為二值圖像時(shí)取10，url1=39。39。 %原始圖像鏈接url2=39。39。 %水印后圖像存儲鏈接url3=39。39。 %水印圖像鏈接xo=。 %水印嵌入key1uo=。M = imread(url1)。R = M(:,:,1)。G = M(:,:,2)。B = M(:,:,3)。f = .*R+.*G+.*B。 %灰度化圖像imshow(f)title(39。水印前39。)。W = mix(xo,uo,imread(url3),0)。 %對水印置亂W = double(W)。[cA,cH,cV,cD]=dwt2(f,39。haar39。)。 %小波分解待水印圖像cA1=edge(cA,39。canny39。)。 %求近似子帶的邊緣點(diǎn)矩陣image_size = size(cA)。long=1。B={}。for i = 1:2:image_size(1) for j = 1:2:image_size(2) k = cA(i:i+1,j:j+1)。 %求分塊矩陣 k1 = cA1(i:i+1,j:j+1)。 if sum(sum(k1)) == 0 %記錄所有不含邊緣點(diǎn)的子塊 B1{1}=k。 B1{2}=mean(mean(k))。 B1{3}=[i,j]。 B{long}=B1。 long=long+1。 end endendlong=long1。W_size=size(W)。m=mean(mean(W))。B=quicksort(B,1,long)。C={}。for i=long:1:longW_size(2)+1 Wl=longi+1。 d=(B{i}{1}(1,1)+B{i}{1}(1,2))/2。 cA(B{i}{3}(1),B{i}{3}(2))=dq*(W(Wl)m)。 cA(B{i}{3}(1),B{i}{3}(2)+1)=d+q*(W(Wl)m)。 C{Wl}=B{i}{3}。 %記錄水印嵌入的塊的坐標(biāo)，key2endf=idwt2(cA,cH,cV,cD,39。haar39。)。 %逆小波變換f=uint8(f)。figureimshow(f) %水印嵌入結(jié)束title(39。水印后39。)。imwrite(f,url2) %保存圖像 function A=quicksort(A,low,high)if low high [A,w]=splition2(A,low,high)。 A=quicksort(A,low,w1)。 A=quicksort(A,w+1,high)。end function [A,w]=splition2(A,low,high)i=low。x=A{low}。for j=low+1:high if A{j}{2}=x{2} i=i+1。 if i~=j。 temp=A{i}。 A{i}=A{j}。 A{j}=temp。 clear temp。 end endendtemp=A{i}。A{i}=A{low}。A{low}=temp。clear temp。w=i。提取水印：function W = getsy(C) %提取水印函數(shù)，參數(shù)為（原圖鏈接）q=15。xo=。uo=。url1=39。39。url2=39。39。f = imread(url1)。[cA,cH,cV,cD]=dwt2(f,39。haar39。)。 %小波分解待水印圖像 K=imread(url2)。W_size1=size(K)。if length(W_size1) == 3 R = K(:,:,1)。 G = K(:,:,2)。 B = K(:,:,3)。 K = .*R+.*G+.*B。 %灰度化圖像endK=double(K)。W_size=W_size1(1)*W_size1(2)。m=mean(mean(K))。for i=1:W_size W(i)=(cA(C{i}(1),C{i}(2)+1)cA(C{i}(1),C{i}(2)))/2/q+m。endW=imix(xo,uo,W)。W=uint8(W)。W1=double(W)。if sum(sum(W1(W1~=1 amp。 W1~=0))) == 0 W=W1。endimshow(imread(url2))title(39。原始水印39。)。figureimshow(W)title(39。提取的水印39。)。4. 一種用哈達(dá)瑪矩陣進(jìn)行置亂變換的盲水印算法水印置亂：function f = mix(f) %參數(shù)為（原始圖像矩陣，置亂前是否灰度化）image_size = size(f)。H = hadamard(image_size(1))。 %生成哈達(dá)瑪矩陣H = flipud(H)。 %對哈達(dá)瑪矩陣置亂f = f*1。 %轉(zhuǎn)換邏輯值為數(shù)值f(f == 1) = 1。f(f == 0) = 1。f = f.*H。水印反置亂：function f = imix(f) %參數(shù)為（置亂后的水印圖像矩陣）image_size = size(f)。H = hadamard(image_size(1))。 %生成哈達(dá)瑪矩陣H = flipud(H)。 %對哈達(dá)瑪矩陣置亂f = f./H。f(f == 1) = 0。f(f == 1) = 1。