【正文】 subplot(2,3,6),imshow(tmark_1,[])。 else tmark_1(ii,jj)=0。tmark_0=blkproc(tmark_0,[8,8],39。Aimage_4=mat2gray(Aimage4)。39。加入白噪聲后的圖象39。disp(39。wimage=abs(wimage_0)。y=(jj1)*8。%%%%%%% 讀取宿主圖像 %%%%%%%cimage=imread(39。title(strcat(num2str(name),39。 else tmark_1(i,j)=0。end%%%%%%%%%%%%% 提取水印 %%%%%%%%%%%%%%tmark_0=att。Aimage4=imrotate(Aimage4,10,39。H=fspecial(39。bmp39。 else k=rand2。alpha=30。4. 將含水印的系數(shù)序列與水印信號進行相關(guān)計算，得到歸一化相關(guān)系數(shù)。計算機只能處理離散信號，且在數(shù)字圖象處理中，輸入圖像和輸出圖像通常都是二維的，一般表示成二維數(shù)字矩陣，因此這里直接討論二維離散傅立葉變換（DFT）和二維快速傅立葉變換（FFT）。用嵌入水印的系數(shù)替代原來位置的系數(shù)，再進行逆DCT（Inverse Discrete Cosine Transform,IDCT）變換，就得到了含有水印的圖像X。 MATLAB函數(shù)介紹在研究水印技術(shù)時，我們處理的圖像數(shù)據(jù)是二維信號，而聲音是一維信號，所以在這里，我門只簡單介紹與水印有關(guān)的函數(shù)。若與不相關(guān)，的概率等于具有Gaussian分布的隨機變量超過其均值T倍方差的概率。由于水印信號被安排在了最底位上，是不可見的，基于同樣的原因，可以輕易的被移去，所以它的魯棒性很差。一般而言，非盲水印的魯棒性較強，而其應(yīng)用受到存儲成本的限制。Ruanaidh和Pun利用傅立葉域?qū)θ中缘男D(zhuǎn)，平移和縮放變換具有不變性的特點，將水印嵌入到傅立葉域來達到對這些攻擊的魯棒性。1. 數(shù)據(jù)庫標識有時一些文件中提示數(shù)據(jù)的標識信息往往比文件本身更重要或者一些音像文件需要將說明注釋（如字幕等）與音像本身結(jié)合起來，這就可以通過數(shù)字水印技術(shù)加以解決。 數(shù)字水印的基本特點數(shù)字水印是加在數(shù)字圖象、音頻或視頻中的微弱信號，這個信號是人們能夠建立產(chǎn)品所有權(quán)、辨認購買者或提供數(shù)字產(chǎn)品的一些額外信息。數(shù)字水印技術(shù)便是二十世紀九十年代初出現(xiàn)的為解決這些問題的一門嶄新的技術(shù)，也是近幾十年來國內(nèi)外專家和學者研究的一大熱點。一方面，它將促進多媒體技術(shù)、網(wǎng)絡(luò)技術(shù)、通信技術(shù)、信號處理技術(shù)等多門新興技術(shù)的有機結(jié)合，促進多媒體網(wǎng)絡(luò)的進一步繁榮；另一方面，它將有助于多媒體信息版權(quán)保護及其沖突問題的解決。宿主信號S是指被嵌入水印的信號（原始信號）。2. 水印檢測算法水印檢測是上述嵌入過程的逆過程，以下詳細介紹一下水印相關(guān)檢測技術(shù)。5. 語言簡單，內(nèi)涵豐富MATLAB語言中最基本最重要的成分是函數(shù)，同一函數(shù)名，不同數(shù)目的輸入變量及不同數(shù)目的輸出變量，代表著不同的含義?；贒CT域的數(shù)字水印算法，一般而言有兩種算法。水印檢測提取過程如圖5圖5 DCT域數(shù)字水印檢測提取框圖水印檢測提取步驟如下：1. 對含有水印圖像進行88的塊DCT變換。嵌入規(guī)則為當水印圖像元素為‘0’時，按照式子（其中alfa為尺度因子，mark為隨機序列，為原圖像的分塊DFT系數(shù)），將一個隨機序列與原始圖像塊的幅度譜對應(yīng)元素進行乘性疊加；為‘1’時，用另一個偽隨機序列與幅度譜對應(yīng)元素進行乘性疊加。 clc。cda0=blkproc(cimage,[8,8],39。 endend%%%%%% 嵌入水印后圖象 %%%%%%wimage=blkproc(cda1,[8,8],39。Aimage1=Aimage1+Wnoise。Aimage3(1:128,1:128)=256。39。 pass(3)=tmark_0(x+3,y+6)。 nc_2=nc_2+g_mark(i,j)*g_mark(i,j)。no=size(omark,2)。fft0=blkproc(cimage,[8,8],39。 fft1(x+7,y+2)=fft_abs(x+7,y+2)+alpha*k(7)。)。begin=input(39。,[4,4],)。)。直接提取的圖象39。 pass(3)=tmark_0(x+3,y+6)。 nc_2=nc_2+g_mark(ii,jj)*g_mark(ii,jj)。另外，還要感謝圖書館電子資源管理室，它在我查閱資料方面也給我提供了許多幫助！感謝所有關(guān)心過我的人，謹以我的努力作為回報！ 30 。for ii=1:m for jj=1:n nc_0=nc_0+g_mark(ii,jj)*o_mark(ii,jj)。 pass(1)=tmark_0(x+1,y+8)。)。imwrite(att,39。H=fspecial(39。5直接檢測39。subplot(2,3,2),imshow(wimage,[]),title(39。 fft1(x+5,y+4)=fft_abs(x+5,y+4)+alpha*k(5)。[mc,nc]=size(cimage)。))/255。for i=1:m for j=1:n nc_0=nc_0+g_mark(i,j)*o_mark(i,j)。 pass(1)=tmark_0(x+1,y+8)。att=Aimage_4。39。disp(‘5直接檢測’)begin=input(‘請選擇攻擊（15）：’)switch begin%%%%%%% 加入白噪聲 %%%%%%%% case 1Aimage1=wimage。 cda1(x+7,y+2)=cda0(x+7,y+2)+alpha*k(7)。subplot(2,3,1),imshow(cimage,[]),title(‘原始的宿主圖像’)?？傊@是一門嶄新而富有挑戰(zhàn)的學科，而且，隨著網(wǎng)絡(luò)的發(fā)展其前景也是非常樂觀的，它需要各個領(lǐng)域的專家聯(lián)合攻關(guān)共同努力才能取得理想的效果。3. 將每個88的DFT系數(shù)矩陣，從每一塊的中頻段取出個系數(shù)。 基于DCT域數(shù)字水印檢測提取算法DCT域水印的驗證過程（即包括水印的檢測和提取過程），就是將待測圖像Y（注意此時的圖像也許已經(jīng)遭受了圖象處理或惡意攻擊）進行DCT運算；然后按著嵌入水印時選取的DCT系數(shù)，選取其含有水印的系數(shù)，進行水印提取，然后利用相關(guān)檢測法判斷水印是否存在。MATLAB中的許多函數(shù)可以直接用來做攻擊測試：1）旋轉(zhuǎn)：rotate()可以對圖像進行任意角度的旋轉(zhuǎn)；2）剪裁：imcrop()可以按精確定位的各點坐標進行剪裁；3）濾波：filter()和filter2()可實現(xiàn)對一維信號和二維信號的濾波；4）加各種噪聲：imnoise()可以對圖像加入各種噪聲，如白噪聲、椒鹽噪聲等。4. 移植性好，開放性好MATLAB是用C語言編寫的，而C語言的可移植性很好。另外基于DFT和DWT算法與上述算法具有相似的原理。數(shù)字水印系統(tǒng)的一般模型如圖1所示：圖1 數(shù)字水印系統(tǒng)基本模型水印嵌入器的輸入量有三個：水印信號M，宿主信號S和密鑰K。數(shù)字水印技術(shù)是近十年來橫跨信號處理、數(shù)字通信、密碼學、計算機網(wǎng)絡(luò)技術(shù)等多學科的新興技術(shù)，具有潛在的應(yīng)用市場和良好的應(yīng)用前景。媒體偵破技術(shù)就可以指出圖像的哪部分被篡改了，指出被插入到原圖像中的新對象等等。而一些具有特殊意義的數(shù)字信息，如涉及司法訴訟、政府機要等信息，更是遭到了不法分子地惡意攻擊和隨意篡改等，這一系列問題給當今科學家?guī)砹司薮筇魬?zhàn)。3. 系統(tǒng)升級日常生活中常常涉及到舊裝系統(tǒng)升級情況，這可以通過將“增強層”嵌入到所發(fā)送的數(shù)據(jù)中來給傳統(tǒng)的信號發(fā)射系統(tǒng)升級。2000年，黃繼武、易開樣等人還提出了一種DCT域數(shù)字水印算法：首先把圖像分成88的不重疊像素塊，經(jīng)過分塊DCT變換后，得到有DCT系數(shù)組成的頻率塊，然后隨機選取一些頻率塊，水印信號嵌入到由密鑰控制選擇的一些DCT系數(shù)中。根據(jù)內(nèi)容劃分水印又可以分為有意義水印和無意義水印。紋理映射編碼方法則是將數(shù)字信息隱藏于數(shù)字圖象的任意紋理部分，該算法對于濾波、壓縮和旋轉(zhuǎn)等操作具有抵抗能力。 MATLAB研究數(shù)字水印的優(yōu)點MATLAB語言有不同于其他高級語言