離散小波變換與奇異值分解的音頻信號(hào)水印算法

馮小明， 馮乃光， 汪云云

(1. 南京郵電大學(xué) 工程訓(xùn)練中心, 江蘇 南京 210003；2. 四川廣播電視大學(xué) 工程技術(shù)學(xué)院, 四川 成都 610072)

?

離散小波變換與奇異值分解的音頻信號(hào)水印算法

馮小明1， 馮乃光2， 汪云云1

(1. 南京郵電大學(xué) 工程訓(xùn)練中心, 江蘇 南京 210003；2. 四川廣播電視大學(xué) 工程技術(shù)學(xué)院, 四川 成都 610072)

針對(duì)現(xiàn)有水印算法在嵌入過程中通常會(huì)對(duì)原始音頻信號(hào)產(chǎn)生破壞的問題，設(shè)計(jì)了一種基于離散小波變換-奇異值分解(DWT-SVD)的水印嵌入方法實(shí)現(xiàn)水印信號(hào)的嵌入與提取.同時(shí)，以相關(guān)系數(shù)和誤碼率為衡量指標(biāo)，對(duì)水印算法的不可感知性、安全性及魯棒性進(jìn)行評(píng)定測(cè)試.結(jié)果表明:在經(jīng)過噪聲、濾波、剪切、壓縮等多種攻擊后，所提取的水印仍然能保持較高的清晰度，這說明DWT-SVD算法具有較強(qiáng)的抗攻擊能力，可有效保護(hù)版權(quán)人的利益.

小波變換； 水印算法； 音頻信號(hào)； 奇異值分解

在原始音頻信號(hào)中嵌入水印是一種有效的產(chǎn)品版權(quán)保護(hù)手段，但在水印嵌入過程中，往往會(huì)使原始音頻信號(hào)發(fā)生改變，使數(shù)字水印所具有的魯棒性受到破壞.基于此，尋找一種有效的音頻水印嵌入方法成為學(xué)者們研究的重點(diǎn).不可感知性和魯棒性是衡量水印算法有效性的根本指標(biāo)[1],按照對(duì)音頻信號(hào)的處理方式，水印嵌入算法可分為時(shí)域法和變換域法[2].時(shí)域法是不經(jīng)過頻域變換，直接把水印嵌入到原始音頻上，嵌入和提取過程較快，缺點(diǎn)是水印比較脆弱，抗攻擊能力較弱[3]，如信息隱藏(LSB)算法、回聲隱藏算法等[4].變換域法是先將采樣數(shù)據(jù)變換為另一個(gè)域內(nèi)的數(shù)字信號(hào)，在加入水印后再逆變?yōu)樵瓉碛蛑械臄?shù)字信號(hào)，該方法具有透明性好、魯棒性強(qiáng)等特點(diǎn)，如離散小波變換、傅里葉變換、倒譜變換等[5-7].在小波變換基礎(chǔ)上，以提高水印魯棒性為目標(biāo)，設(shè)計(jì)一種基于DWT-SVD的組合音頻水印算法.

1 基于DWT-SVD的零水印算法

1.1 水印圖像預(yù)處理

設(shè)原始水印是一個(gè)二值圖像，可以表示為

(1)

式(1)中：m1為二值圖像的寬度；m2為二值圖像的高度.因此,音頻信號(hào)為一維信號(hào)，為將水印信號(hào)嵌入音頻信號(hào)中，對(duì)二值圖像進(jìn)行降維處理，可得到二值序列，即

(2)

Logistic混沌序列具有較好的互相關(guān)性和自相關(guān)性，用其進(jìn)行水印圖像的加密可保證水印具有較好的保密性，進(jìn)行二值化處理后，可得

(3)

由此得到置亂后的水印圖像為

(4)

1.2 水印的嵌入過程

1) 將原始音頻信號(hào)分段，分成m1×m2個(gè)幀，每幀的長(zhǎng)度為1 024樣點(diǎn).2) 將三級(jí)小波變換得到的低頻分量ca3進(jìn)行降維，得到一維矩陣A0.3) 將置亂后的水印和矩陣A0分成16×16大小的塊.4) 將每個(gè)小塊進(jìn)行奇異值分解，得到奇異值矩陣Si,j和Sωi,j.5) 將Sωi,j嵌入矩陣Si,j中，完成水印信息的嵌入，Sai,j=Si,j+α×Sωi,j.式中:α為嵌入強(qiáng)度，且0<α<1.6) 將嵌入水印的矩陣進(jìn)行SVD變換，分解為含有水印的小塊Aai,j.7) 將含有水印的小塊進(jìn)行整合，將得到的矩陣進(jìn)行逆變換，最后得到含有水印的音頻信號(hào)，從而完成音頻水印的嵌入.水印嵌入流程,如圖1所示.

1.3 水印的提取過程

水印提取過程即嵌入水印的逆過程.1) 將嵌入水印信號(hào)的音頻進(jìn)行三級(jí)小波變換.2) 低頻分量ca3升維，得到二維矩陣A1.3) 將矩陣A1分成16×16大小的塊，每個(gè)小塊中都含有水印信號(hào).4) 將音頻塊奇異值分解，得到矩陣Sai,j,Uai,j和Vai,j.5) 為得到嵌入水印奇異值矩陣，對(duì)水印嵌入公式進(jìn)行逆運(yùn)算，Sωi,j=(Sai,j-Si,j)/α.6) 將含水印的奇異值矩陣反變換，得到水印小塊Aωi,j.7) 將每個(gè)小塊整合為16×16的矩陣，實(shí)現(xiàn)水印的提取過程.水印提取流程，如圖2所示.

圖1 音頻水印信號(hào)的嵌入流程 圖2 音頻水印信號(hào)的提取流程Fig.1 Embedding process of audio watermark signal Fig.2 Extraction process of audio watermark signal

2 仿真實(shí)驗(yàn)與結(jié)果分析

2.1 不可感知性測(cè)試

為了測(cè)試嵌入水印對(duì)原始音頻信號(hào)的影響，采用Matlab分析軟件進(jìn)行聽覺感知仿真，選用單聲道、16 bit原始音頻文件，采樣頻率為44.1 kHz，音頻格式為WAV.水印圖像采用16 px×16 px的二值圖像.采用Logistic混沌序列加密，嵌入強(qiáng)度α為0.2，根據(jù)原始音頻和水印信息的長(zhǎng)度進(jìn)行采樣分幀，仿真結(jié)果,如圖3所示.由圖3可知：原始音頻信號(hào)在水印圖像嵌入的前后，音頻信號(hào)的波形幅度沒有發(fā)生明顯的變化.經(jīng)計(jì)算嵌入水印前后音頻信號(hào)信噪比為35.036 7 dB.該數(shù)值在聽覺上很難感覺到音頻水印的存在，說明水印算法具有較好的不可感知性.

(a) 嵌入水印前的音頻信號(hào) (b) 嵌入水印后的音頻信號(hào) (c) 嵌入水印前后的波形差圖3 水印嵌入前后音頻信號(hào)波形對(duì)比Fig.3 Comparison of audio signal waveform before and after watermark embedding

2.2 安全性測(cè)試

為檢驗(yàn)Logistic混沌序列加密的可靠性與安全性，對(duì)嵌入的水印信號(hào)進(jìn)行保密性測(cè)試，分別采用正確的和錯(cuò)誤的密鑰對(duì)水印圖像進(jìn)行提取，分別與原始水印圖像對(duì)比，其結(jié)果如圖4所示.由圖4可知：當(dāng)采用正確密鑰提取水印時(shí)，得到的水印圖像非常清晰，與原始水印相比，基本沒有發(fā)生任何改變.由此可見，采用Logistic混沌序列對(duì)音頻水印信號(hào)進(jìn)行加密和置亂具有較好的安全性，滿足保密性要求.

(a) 原始水印 (b) 正確密匙下提取的水印 (c) 錯(cuò)誤密匙下提取的水印 圖4 水印圖像的安全性測(cè)試Fig.4 Security testing of watermark image

2.3 魯棒性測(cè)試

為了測(cè)試算法魯棒性，分別對(duì)嵌入水印的古典音樂和搖滾音樂信號(hào)進(jìn)行攻擊實(shí)驗(yàn)，具體攻擊內(nèi)容如下.1) 高斯噪聲.在待測(cè)試的音頻信號(hào)中加入均值為0的正態(tài)分布高斯噪聲.2) 重采樣.改變待測(cè)試音頻的采樣頻率，由44.1 kHz變?yōu)?2.05 kHz.3) 低通濾波.采用2階低通濾波器對(duì)待測(cè)試音頻進(jìn)行濾波，濾波器的截止頻率為4 Hz.4) 重量化.先將待測(cè)試音頻從16 bit量化為8 bit,再量化為16 bit.5) 剪切.將待測(cè)試音頻隨機(jī)裁剪掉2 s.6) 小波去噪.利用一維小波分析方法對(duì)信號(hào)進(jìn)行降噪處理.7) MP3壓縮.將待測(cè)試音頻壓縮到56 Kb·s-1.攻擊完成后，對(duì)水印圖像進(jìn)行提取，為進(jìn)一步分析水印的清晰程度，分別采用相關(guān)系數(shù)和誤碼率作為評(píng)價(jià)指標(biāo)對(duì)其進(jìn)行判定，相關(guān)系數(shù)為

(5)

誤碼率RBE是檢驗(yàn)音頻水印信號(hào)抗攻擊能力的另一個(gè)指標(biāo)，其表達(dá)式為

(6)

式(6)中：Ebit為錯(cuò)誤的比特?cái)?shù)量；Tbit為總的比特?cái)?shù)量.

受到攻擊后提取的水印圖像及其屬性值，如表1所示.由表1可知：受到攻擊后的音頻信號(hào)提取出來的水印圖像大多清晰可見，基本不受各攻擊操作的影響，只有MP3壓縮水印效果稍顯模糊，但水印圖像的內(nèi)容仍然可以很容易分辨出來，這說明所采用的水印算法抗攻擊能力強(qiáng)、可行有效.

為了進(jìn)一步驗(yàn)證DWT-SVD算法的有效性，與單純采用小波變換作為水印嵌入算法進(jìn)行魯棒性對(duì)比測(cè)試，測(cè)試過程中，選取相同的古典音樂作為原始音頻信號(hào)，采用同一格式的圖像作為水印.水印嵌入后，分別對(duì)兩個(gè)含水印的音頻信號(hào)進(jìn)行魯棒性測(cè)試，并采取相關(guān)系數(shù)及誤碼率進(jìn)行評(píng)定，結(jié)果如表2所示.由表2可知:與只采用小波變換水印算法相比，采用小波變換與奇異值分解相結(jié)合的水印算法各項(xiàng)性能都大幅提升，對(duì)各種攻擊的抵抗能力更強(qiáng).

表1 受到攻擊后提取的水印圖像及其屬性值

表2 魯棒性測(cè)試結(jié)果

3 結(jié)束語

針對(duì)嵌入過程容易造成原始信號(hào)損壞的問題，提出一種DWT-SVD組合音頻水印算法.該方法結(jié)合了小波變換多分辨率與奇異值分解抗干擾能力強(qiáng)的特點(diǎn)，嵌入過程不改變?cè)妓⌒畔?，同時(shí)基于Logistic混沌序列對(duì)二值水印信息進(jìn)行加密，以滿足水印的安全性要求.測(cè)試結(jié)果表明，DWT-SVD水印算法可大幅提高水印的相關(guān)系數(shù)，有效降低水印信息的誤碼率，對(duì)數(shù)字音頻產(chǎn)品的版權(quán)具有一定的保護(hù)作用.

[1] 楊榆,雷敏,鈕心忻,等.DWT-DCT-QR結(jié)合的音頻盲水印算法[J].重慶大學(xué)學(xué)報(bào),2012,35(8):62-66.

[2] 吳倪.基于DCT域的數(shù)字水印算法及其在版權(quán)保護(hù)中的應(yīng)用[D].長(zhǎng)沙：中南大學(xué),2013：12-15.

[3] 雷敏.音頻數(shù)字水印與隱寫分析算法研究[D].北京：北京郵電大學(xué),2011：23-27.

[4] 林曉丹.一種基于支持向量機(jī)的數(shù)字音頻認(rèn)證方法[J].華僑大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,32(3):153-155.

[5] 于曉敏,于曉坤,耿蕊.基于離散小波變換的盲水印研究與仿真[J].計(jì)算機(jī)仿真,2011,28(8):287-290.

[6] 劉素楠,方志宏,李淑芝.采用混沌加密和重復(fù)碼的復(fù)倒譜音頻水印算法[J].計(jì)算機(jī)工程與應(yīng)用,2010,46(4):79-82.

[7] 董斌.基于變換域的音頻水印研究及實(shí)現(xiàn)[D].長(zhǎng)春：吉林大學(xué),2011：17-21.

[8] 全笑梅,張鴻賓.基于小波包域聽覺感知分析的統(tǒng)計(jì)音頻水印算法[J].電子學(xué)報(bào),2007,35(4):673-678.

[9] 彭維,高健,孫瑞鵬,等.基于DCT 系數(shù)比較的音頻水印算法[J].計(jì)算機(jī)應(yīng)用與軟件,2014,31(11):158-160.

[10] 馮祥斌,陳永紅.應(yīng)用P-Fibonacci加密的模糊自適應(yīng)水印算法[J].華僑大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,35(5):227-292.

(責(zé)任編輯： 錢筠 英文審校： 吳逢鐵)

Watermarking Algorithm of Audio Signal Based on Discrete Wavelet Transform and Singular Value Decomposition

FENG Xiaoming1, FENG Naiguang2, WANG Yunyun1

(1. Engineering Training Center, Nangjing University of Posts and Telecommunications, Nangjing 210003, China;2. College of Engineering and Technology, Sichuan Radio and Television University, Chengdu 610073, China)

The existing watermarking algorithm in the embedding process usually issues the destruction of the original audio signal, and design a watermark embedding method based on the discrete wavelet transform-singular value decomposition (DWT-SVD) to realize the watermark embedding and extraction. At the same time, with the correlation coefficient and bit error rate as the measuring index, the non-sense, security and robustness of the watermarking algorithm are evaluated and tested. Results show that after noise, filtering, cropping, compression multiple attacking, the extracted watermark can still maintain a high clarity, which indicates that DWT-SVD algorithm has strong anti-attack ability, can effectively protect the copyright owner′s interests.

wavelet transform; watermarking algorithm; audio signal; singular value decomposition

10.11830/ISSN.1000-5013.201606022

2016-10-13

馮小明(1960-)，女，高級(jí)工程師，主要從事計(jì)算機(jī)技術(shù)的研究.E-mail:671003753@qq.com.

國家自然科學(xué)基金青年基金資助項(xiàng)目(61300165)

TP 751.1

A

1000-5013(2016)06-0770-04