【正文】 水印算法[J]．計算機安全，2009，16（1）：6567．[4] 劉宏斌，王麗君，王克剛，陳升．小波域中的魯棒性盲水印算法[J]．電子技術(shù)應用，2008，27（3）：122125．[5] 強世東，強英．基于小波域的自適應圖象盲水印算法[J]．計算機工程與應用，2005，41（24）：4749．[6] 王友衛(wèi)，申鉉京，陳海鵬．基于Slant變換的半脆弱水印算法[J]．北京工業(yè)大學學報，2013，39（1）：7680．[7] Emir Ganic, Ahmet M．A DFTBased SemiBlind Multiple Watermarking Scheme For Images[C]．Eskicioglu．Bedford Avenue CUNY Brooklyn College, 2009: 112134．[8] 何冰．基于SVD和Radon變換的抗旋轉(zhuǎn)攻擊盲水印算法[J]．計算機工程與應用，2012，48（20）：200205．[9] 林愛英，吳莉莉，蘇金，鄭寶周．一種用哈達瑪矩陣進行置亂變換的盲水印算法[J]．合肥工業(yè)大學學報(自然科學版)，2010，33（9）：13511354．[10] 劉珍慧，馮桂，鄭劍雄．一種新的小波域盲水印算法[J]．信息安全與通信保密，2006，13（9）：126127．[11] 葉天語．基于子塊區(qū)域分割和自嵌入技術(shù)的全盲多功能圖像水印算法[J]．通信學報，2013，34（3）：148156．致 謝致 謝這篇論文的完成過程充滿了艱辛，但也飽含著溫暖．在論文的寫作過程中我遇到了無數(shù)的困難和障礙，是老師和同學的無私幫助讓我度過了難關(guān)．尤其要強烈感謝我的論文指導老師—常志文老師，他對我進行了無私的指導和幫助，不厭其煩的幫助我進行論文的修改和改進，是一位耐心負責的好老師！同時要感謝李鵬松老師、朱秀麗老師、劉洪偉老師和張中華老師，他們耐心細致地幫助我進行論文的后期修改，使我的論文更加完美，他們都是敬業(yè)負責的好老師，好園??！感謝各位在編程過程中幫我解決技術(shù)性問題的同學，謝謝你們！感謝我的舍友，在論文的撰寫和排版過程中提供熱情的幫助！感謝這篇論文所涉及到的各位學者，本文引用了各位學者的研究文獻，如果沒有你們的研究成果的幫助和啟發(fā)，我將很難完成本篇論文的寫作，謝謝各位！在此向幫助和指導過我的各位表示最誠摯的謝意！由于我的學術(shù)水平有限，所寫論文難免有不足之處，懇請各位老師和學友批評和指正！ 劉世江2015年6月附 錄附 錄1. 基于SVD和Radon變換的抗旋轉(zhuǎn)攻擊盲水印算法水印置亂：function f=mix(f,n) %參數(shù)為（原始圖像矩陣，置亂次數(shù)）image_size=size(f)。for x=1:n for i=1:image_size(1) for j=1:image_size(2) if i j a=[1,1。1,0]*[i。j]+(image_size39。+[1。1])。 g(a(1),a(2))=f(i,j)。 else a=[1,1。1,0]*[i。j]+([0。image_size(2)]+[1。1])。 g(a(1),a(2))=f(i,j)。 end end end f=g。end水印反置亂：function f=imix(M,n) %參數(shù)為（置亂后的水印圖像矩陣，置亂次數(shù)）f = M。image_size=size(f)。for x=1:n for i=1:image_size(1) for j=1:image_size(2) if i + j = image_size(1)+1 a=[0,1。1,1]*[i。j]+(image_size39。+[1。2])。 g(a(1),a(2))=f(i,j)。 else a=[0,1。1,1]*[i。j]+(image_size39。+[1。1]).*[1。2]。 g(a(1),a(2))=f(i,j)。 end end end f=g。end嵌入水印：function f = sy(url,url2,n) %水印函數(shù)f=imread(url)。imshow(f)title(39。水印前39。)。[cA,cH,cV,cD]=dwt2(f,39。haar39。)。[U,S,V]=svd(cA)。W1=mix(imread(39。39。),n)。b=2。Y=round(S/b)。image_size=size(S)。for i=1:image_size(1) for j=1:image_size(2) if mod((Y(i,j)+W1(i,j)),2)==1 K(i,j)=(Y(i,j))*b。 else K(i,j)=(Y(i,j)+)*b。 end endend[U1,S1,V1]=svd(K)。A1=U*S1*V39。f=idwt2(A1,cH,cV,cD,39。haar39。)。figureimshow(uint8(f))title(39。水印后39。)。imwrite(uint8(f),url2) %保存圖像提取水?。篺unction [shuiyin,T] = getsy(url1,url2,n) %參數(shù)為（原圖鏈接，水印后圖像鏈接，水印置亂次數(shù)）f=imread(url1)。imshow(f)title(39。水印前39。)。[cA,cH,cV,cD]=dwt2(f,39。haar39。)。 %二維小波分解[U,S,V]=svd(cA)。 %奇異值分解W1=mix(imread(39。39。),n,0)。 %對水印置亂b=2。 %量化步長Y=round(S/b)。 %計算量化值image_size=size(S)。for i=1:image_size(1) for j=1:image_size(2) if mod((Y(i,j)+W1(i,j)),2)==1 K(i,j)=(Y(i,j))*b。 else K(i,j)=(Y(i,j)+)*b。 end endend[U1,S1,V1]=svd(K)。 %對K進行奇異值分解A1=U*S1*V39。 %合成新的低頻子帶f=idwt2(A1,cH,cV,cD,39。haar39。)。 %用新的低頻子帶替換原低頻子帶，進行逆離散小波變換figureimshow(uint8(f)) %水印嵌入結(jié)束title(39。水印后39。)。%imwrite(uint8(f),url2) %保存圖像 %前面的步驟是為了得到U1和V1，即參考向量f=imread(url2)。 %讀取水印后的圖像[A1,cH,cV,cD]=dwt2(f,39。haar39。)。 %提取有水印的低頻子帶[U,S1,V]=svd(A1)。 %提取嵌入的對角矩陣K=U1*S1*V139。 %重新合成KY=floor(K/b)。 %逆推出量化值image_size=size(Y)。for i=1:image_size(1) %逆推出置亂后的水印 for j=1:image_size(2) if mod(Y(i,j),2)==1 W1(i,j)=1。 else W1(i,j)=0。 end endendshuiyin=imix(W1,n)。 %進行置亂反變換得到原始水印figureimshow(imread(39。39。))title(39。原始水印39。)。figureimshow(shuiyin)title(39。提取的水印39。)。T=sum(sum(abs(imread(39。39。)shuiyin)))。 %計算重構(gòu)的水印圖像和原始水印圖像的相似度if T = 150 fprintf(39。有水?。n39。)。else fprintf(39。無水印！\n39。)。end2. 基于子塊區(qū)域分割和自嵌入技術(shù)的全盲多功能圖像水印算法嵌入水?。篺unction B=sy(n,m,p) %n為分割份數(shù)，m為子矩陣區(qū)域1行列數(shù)，p為水印嵌入強度，測試參數(shù)B=sy(32,3,20)。M=imread(39。39。)。B=.*M(:,:,1)+.*M(:,:,2)+.*M(:,:,3)。 %圖像灰度化subplot(1,2,1)。imshow(B)。S=size(B)。k=1。for i=1:n:S(1) for j=1:n:S(2) Bj{k}=dct2(B(i:i+n1,j:j+n1))。 %分割矩陣并做離散余弦變換 V(k)=mean(svd(Bj{k}(1:m,1:m)))。 %多每個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。 k=k+1。 end。end。L=hundun(,2,300,k1)。 %生成混沌序列W1=xor(L,W)。 %異或運算加密特征水印 x1=5。y1=1。x2=4。y2=1。 %嵌入水印位置for i=1:k1 %嵌入水印 f1=Bj{i}(x1,y1)Bj{i}(x2,y2)。 if W1(i)==0 amp。amp。 f1p Bj{i}(x1,y1)=Bj{i}(x1,y1)+(pf1)/2。 Bj{i}(x2,y2)=Bj{i}(x2,y2)(pf1)/2。 end。 f2=Bj{i}(x2,y2)Bj{i}(x1,y1)。 if W1(i)==1 amp。amp。 f2p Bj{i}(x1,y1)=Bj{i}(x1,y1)(pf2)/2。 Bj{i}(x2,y2)=Bj{i}(x2,y2)+(pf2)/2。 end。end。k=1。for i=1:n:S(1) for j=1:n:S(2) B(i:i+n1,j:j+n1)=idct2(Bj{k})。 %做逆離散余弦變換后覆蓋生成圖像 k=k+1。 end。end。subplot(1,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)。 %二值化end提取水印：function [NC,TAF]=getsy(n,m,p,er,T) %n為分割份數(shù)，m為子矩陣區(qū)域1行列數(shù)，p為水印嵌入強度，er為判斷篡改類型的閥值，測試參數(shù)[p,q]=getsy(32,3,20,1)。x1=5。y1=1。x2=4。y2=1。 %提取水印位置B=sy(n,m,p)。 if T==1 M=imread(39。39。)。 H=.*M(:,:,1)+.*M(:,:,2)+.*M(:,:,3)。 %圖像灰度化 B(189:220,165:196)=B(189:220,165:196)+H*50。 %惡意篡改 %B(1:220,1:196)=255。 %非惡意篡改endfigureimshow(B)figuresubplot(2,2,1)。imshow(B)。S=size(B)。k=1。for i=1:n:S(1) for j=1:n:S(2) Bj{k}=dct2(B(i:i+n1,j:j+n1))。 %分割矩陣并做離散余弦變換 V(k)=mean(svd(Bj{k}(1:m,1:m)))。 %多每個