基于相關標識符的頻譜水印嵌入與提取方法

張余　許金勇　柳永祥　劉忠英　陳勇

(南京電訊技術研究所,南京 210007)

?

基于相關標識符的頻譜水印嵌入與提取方法

張余許金勇柳永祥劉忠英陳勇

(南京電訊技術研究所,南京 210007)

摘要為有效標識無線電設備身份信息，快速確定無線電干擾，提出了一種基于相關標識符的頻譜水印嵌入與提取方法.該方法中發射端在標識符嵌入控制模塊的控制下將頻譜水印嵌入位置標識符嵌入到通信信號中，然后通過對含有頻譜水印嵌入位置標識符的通信信號小波分解低頻系數進行量化來嵌入頻譜水印；接收端通過分析相關標識符的自相關和互相關特性來檢測頻譜水印，對接收到的中頻信號小波分解低頻系數實施量化處理來提取頻譜水印.理論分析與仿真結果表明:提出的方法能夠在不對原始信號正常解調產生有害干擾的條件下，可有效提取頻譜水印，不需解調原始信號，可有效減少頻譜水印提取的復雜度，為有效標識和識別無線電設備身份提供了技術手段.

關鍵詞頻譜水印；相關標識符；小波分解；頻譜水印嵌入；頻譜水印提取

DOI10.13443/j.cjors.2015031701

Spectrum watermark embedding and extracting method based on correlation identifier

ZHANG YuXU JinyongLIU YongxiangLIU ZhongyingCHEN Yong

(NanjingTelecommunicationTechnologyInstitute，Nanjing210007，China)

Abstract To effectively identify the radio equipment and quickly ascertain the radio interference， this paper presents a spectrum watermark embedding and extracting method based on the correlation identifier. The correlation identifier that identified the embedding location of spectrum watermark was firstly embedded into communication signals under the control of identifier embedded control module in the transmitting terminal. Then the spectrum watermark was embedded by quantizing the low-frequency of wavelet decomposition of the communication signal with correlation identifier. The spectrum watermark was detected by analyzing these characteristic of autocorrelation and cross-correlation of correlation identifier. And it was extracted by quantizing the low-frequency of wavelet decomposition of received signals in the receiving terminal. Theory analysis and simulation results show that the proposed method can effectively extract the spectrum watermark without harmful inference to normal demodulation of communication signals， and there is no need to demodulate the communication signal， which can effectively reduce the complexity of spectrum watermark extracting. Thus， it provides a technical support for effectively identifying and ascertaining the radio equipment.

Keywordsspectrum watermark； correlation identifier； wavelet composition； spectrum watermark embedding； spectrum watermark extracting

引言

隨著信息技術的發展，人們生活與國防建設對無線電應用的依賴越來越強，無線電設備種類越來越多，數量越來越大，相互干擾日趨嚴重，無意輻射和超標輻射等現象經常發生，給無線電設備使用與管理帶來了巨大威脅，如何快速識別干擾源、查出干擾是頻譜管控急需解決的難題之一.無線電干擾源識別問題引起了學者們廣泛關注.頻譜水印技術[1]利用待傳信息中冗余空間嵌入攜帶設備身份的隱蔽信息，達到秘密信息偽裝傳遞之目的，并通過檢測隱蔽信息對輻射源進行識別.利用水印技術特點，可將無線電設備身份信息嵌入到用頻設備的信息傳輸中，在不對原始信息正常傳輸產生有害干擾條件下，對無線電設備進行身份標識，為無線電設備管理提供一種有效的身份識別途徑.頻譜水印技術來源于數字水印技術，目前還處于起步階段，但數字水印技術已取得許多成果.季秀蘭等人[2]為加強數字圖像版權保護，充分利用支持向量機非線性逼近能力計算水印嵌入強度，將水印自適應地嵌入到離散余弦變換各系數中，解決了水印嵌入強度和離散余弦變換系數的選擇難題.Naresh Kumar Trivedi等人[3]利用心理聲學聽覺模型來計算音頻信號掩蓋門限，使其嵌入的使用信息不可感知，并對水印實施擴頻調制來提高其魯棒性；Priyank Khare等人[4]為保護圖像知識產權，利用邏輯圖特征來產生混沌序列提高圖像水印的魯棒性. Aniyan， A.等人[5]針對數字圖像提出基于離散余弦變換的盲水印方法，并開展了硬件實驗工作，實驗結果顯示提出的方法能有效應對多種攻擊.Trabelsi， W. 等人[6]為在魯棒性、透明性和性能間找到均衡點，利用基于奇異值分解的頻域水印嵌入方法來嵌入水印，并插入一個雙奇異值來提高水印性能.趙軍等人[7]針對基帶信號水印嵌入算法中水印同步問題，提出了利用水印圖像清晰度評價技術識別幀同步位方法，該方法中由于在發射端基帶嵌入水印，在接收端需解調原始信號才能有效提取水印.

現有數字水印研究成果中大多要求接收端解調出原始信號，而日常監測中難以做到，其通用性和適用性較差.雖然CN 102780502 A國家發明專利[8]公開了一種為無線電臺施加水印的方法，不需解調原始信息，但為了不影響原始信號解調，發射振蕩器相位變化不能太大，接收端對細微特征的捕捉很困難，接收機實現難度大.因此，文章提出了一種基于相關標識符的頻譜水印嵌入與提取方法.該方法利用相關碼作為頻譜水印嵌入位置的標識，可減少接收端頻譜水印存在性檢測的復雜度，有利于頻譜水印嵌入起始位置的精確定位；選擇原始中頻信號小波分解低頻系數作為頻譜水印的嵌入位置，可提高頻譜水印的魯棒性；利用量化準則來嵌入與提取頻譜水印信號，在頻譜水印提取時不需解調原始無線電信號，可提高其通用性，減少接收端頻譜水印提取的復雜度.

1基于相關標識符的頻譜水印嵌入方法

為使己方接收設備能容易地發現和提取頻譜水印而第三方難以提取，在嵌入頻譜水印時，首先選擇一個自相關性強、互相關性弱的相關碼來標識頻譜水印的嵌入，然后對含有標識信息的通信信號進行小波分解，并通過對小波分解系數的量化來實現頻譜水印的嵌入，在確保不影響原始通信信號正常傳輸的同時，將標識無線電設備身份的頻譜水印隱蔽地傳遞出去.圖1給出了基于相關標識符的頻譜水印嵌入方法的實現過程，主要包括標識符生成與嵌入、頻譜水印生成、頻譜水印嵌入等過程，圖中字母或變量的含義見下文.

圖1　基于相關標識符的頻譜水印嵌入的實現過程

1.1標識符生成與嵌入

假設頻譜水印嵌入位置的標識符序列為b(n)，經基帶調制和脈沖成型后的基帶標識符信號z(t)可表示為

(1)

式中，g3(t)，g4(t)分別為二進制序列1和0對應的脈沖成型波形.為便于分析，通常選用存在g3(t)=-g4(t)關系的波形來成型.經載波調制后的中頻標識符信號zc(t)可表示為

zc(t)=z(t)cos(2πfct+φz).

(2)

式中: fc為無線電設備發射端載波頻率； φz為頻譜水印嵌入位置標識符載波調制后的初始相位.

為確保頻譜水印嵌入位置標識符的嵌入不擾亂無線電設備信息的傳輸，標識符只能嵌入到通信信號的有用信息序列后、同步序列前，因此，標識符嵌入必須在標識符嵌入控制模塊控制下實現嵌入，標識符嵌入控制模塊將記錄無線電設備信號有用信息和同步碼的位置，并將此信息傳遞給信號合成器來控制頻譜水印嵌入位置標識符的嵌入.

信號合成器在標識符嵌入控制模塊的控制下，將頻譜水印嵌入位置標識符周期性嵌入到原始信號序列的上一幀有用序列后和下一幀同步序列前，得到含有頻譜水印嵌入位置標識符后的信號y(l)，可表示為

(3)

式中， yi(l)為第i段含有頻譜水印嵌入位置標識符的通信信號，可表示為

yi(l)=

(4)

式中: tz0為嵌入頻譜水印嵌入位置標識符的起始時間； tz為嵌入頻譜水印嵌入位置標識符的持續時間； txl為嵌入頻譜水印的周期，必須保證頻譜水印所有信息在此周期內能夠完全被嵌入，且包括M幀原始信號數據，若每幀數據進行載波調制后的持續時間為tf，則txl=M×tf.tz0，tz，txl和tf的單位均為s.

嵌入頻譜水印嵌入位置標識符將增加系統的開銷，若頻譜水印嵌入位置標識符嵌入過于頻繁，頻譜水印嵌入位置標識符序列長度過長，增加開銷將更多，系統將無法忍受，設計合適的頻譜水印嵌入標識符是保證頻譜水印透明度和頻譜水印檢測魯棒性的關鍵.為便于頻譜水印嵌入位置標識符的檢測，在此選用自相關性強、互相關性弱的16bits的m序列作為頻譜水印嵌入位置標識符，在頻譜水印提取端充分利用其自相關特性來檢測頻譜水印的嵌入位置，既可降低其檢測復雜度，也可保證增加的開銷在系統忍受的范圍之內.

1.2頻譜水印生成

通信信號不存在像音頻信號那樣具有人類聽覺系統的掩蔽效應[9]，受噪聲影響較大，若采用圖像形式的水印來標識其頻譜身份，將會影響通信信號本身的特征，因此，在此采用編號序列作為頻譜水印來標識無線電設備身份.

頻譜水印主要由無線電設備型號和編號等信息組成，如編號為3 201的超短波調頻電臺的頻譜水印可用12個字符96bits的序列V來表示.為增加頻譜水印檢測的魯棒性，提高頻譜水印檢測的正確率，用擴頻增益為P的擴頻碼對頻譜水印序列進行擴頻調制，得到擴頻調制后的水印序列為W，擴頻調制后的序列長度為96Pbits.

1.3頻譜水印嵌入

頻譜水印嵌入包括分段通信信號采樣、小波分解、頻譜水印信號嵌入和小波重構等過程.

1) 分段通信信號采樣

以采樣率fs對第i段含有頻譜水印嵌入位置標識符的通信信號進行采樣，得到第i段離散通信信號yi(k)，可表示為

yi(k)=

(5)

式中，[u]表示對變量u向下取整； [u]表示對變量u向上取整.用矢量形式可表示為

yi=[zi，xi].

(6)

式中:zi=z(m)，[((i-1)(tz+txl)+tz0)fs]≤m<[((i-1)(tz+txl)+tz+tz0)fs]；xi=x(m)，[((i-1)(tz+txl)+tz+tz0)fs]≤m<[(i(tz+txl)+tz0)fs].

2) 小波分解

對第i段含有頻譜水印嵌入位置標識符的通信信號yi進行n層小波分解，得到第n層小波分解的低頻系數Ca，n和各層小波分解的高頻系數Cd，n、…、Cd，2、Cd，1，其中第n層小波分解低頻系數的長度為[(tz+txl)fs]/2n，各層小波分解高頻系數的長度分別為[(tz+txl)fs]/2n、…、[(tz+txl)fs]/22、[(tz+txl)fs]/2，且第n層小波分解低頻系數的長度必須滿足大于經擴頻調制后的頻譜水印的序列長度，即[(tz+txl)fs]/2n≥96P.

3) 頻譜水印信號嵌入

通信信號可通過小波分解分成近似分量和細節分量，近似分量集中了信號的大部分能量，可表征信號的絕大多數信號特征，而細節分量僅表示信號的細微特征，受噪聲影響較大，容易受濾波、除噪等信號處理的影響[10-11].因此，選擇原始中頻信號小波變換的近似分量作為水印嵌入的對象.

由于量化水印方法[12]是根據不同的水印信息用不同量化器去量化原始信號，不像常見數字水印中利用的加性準則[13]、乘性準則[13-14]和指數準則[13]的水印嵌入方法那樣在水印提取端要求知道原始信號或原始信號變換域的變換系數，因此，以嵌入系數α為量化步長對各小波分解低頻系數Ca，n進行量化，并根據頻譜水印信號與量化結果間的關系嵌入水印.

利用量化準則可得到小波分解低頻系數的量化結果，即

Cq，n(m)=[|Ca，n(m)|/α+1/2].

(7)

式中:Cq，n(m)為第n層中頻信號的小波分解第m個低頻系數的量化結果；Ca，n(m)為第n層中頻信號的小波分解第m個低頻系數； |h|表示對h的取模運算.

為便于頻譜水印檢測，在此以分段中頻信號小波分解的第一個低頻系數為起點嵌入水印信號.

當m≤96P時，若2mod(Cq，n(m)，2)-1=w(m)，則

Cb，n(m)=sign(Ca，n(m))Cq，n(m)α.

(8)

式中: mod(h，2)表示對h取模2的余數； sign(h)表示取h的符號.若2mod(Cq，n(m)，2)-1≠w(m)，但如Cq，n(m)=[|Cq，n(m)|/α]，則

Cb，n(m)=sign(Ca，n(m))(Cq，n(m)+1)α.

(9)

若2mod(Cq，n(m)，2)-1≠w(m)，但如Cq，n(m)≠[|Cq，n(m)|/α]，則

Cb，n(m)=sign(Ca，n(m))(Cq，n(m)-1)α.

(10)

當m>96P，Cq，n(m)=Ca，n(m).

4) 中頻信號重構

2基于相關標識符的頻譜水印檢測與提取方法

基于相關標識符的頻譜水印檢測與提取的實現過程如圖2所示，主要包括頻譜水印信號檢測與嵌入位置定位和頻譜水印提取，圖中字母和變量的含義見下文.

圖2　基于相關標識符的頻譜水印檢測與提取的實現過程

2.1頻譜水印信號檢測與嵌入位置定位

(11)

(12)

式中，ti為滑動窗第i次滑動的起始時間.然后求h(t)的功率譜密度S(f)，即

S(f)=|FFT(h(t))|2.

(13)

只要通過搜索S(f)的相關峰值來確定是否存在頻譜水印信號，若有相關峰存在，則停止搜索，否則繼續滑動窗口，進行相乘運算和功率譜運算，直到出現較大的峰值.

由于頻譜水印提取過程與頻譜水印嵌入過程類似，在提取頻譜水印時必須精確定位頻譜水印，否則難以準確提取頻譜水印，因此，精確定位頻譜水印嵌入標識符是能否成功提取頻譜水印的關鍵.當檢測到頻譜水印存在后，利用現有載頻同步算法精確估計頻譜水印嵌入標識符的相位，實現頻譜水印嵌入位置的精確定位.

2.2頻譜水印提取

根據頻譜水印嵌入標識符定位結果，對接收到的中頻信號實施分段和采樣，獲得第j段離散中頻信號rj(n).假設頻譜水印檢測端已知分段長度tz+txl和采樣率fs，對于半合作狀態的頻譜水印檢測而言，此假設是合理的.

Uj={uj(k)|，0≤k≤P}，

(14)

(15)

根據以上方法，對接收端獲得的d段中頻信號進行相同處理，提取頻譜水印，最后根據大數定理，對恢復的頻譜水印序列實施大數判決，獲得頻譜水印序列綜合判決結果V′.

3實驗結果分析

為驗證算法的性能，在此對頻譜水印進行了透明性、魯棒性、水印提取正確概率和嵌入量仿真分析實驗.其中透明性用于衡量所嵌入的水印信息不能影響原始通信信號的質量，在此，利用信噪比來衡量通信信號的客觀質量評價指標，它將嵌入的水印信息看作是添加到原始通信信號里的噪聲，用來評價嵌入水印后的通信信號與原始通信信號間的失真度；魯棒性用來表示頻譜水印算法抗外界干擾的能力，可以利用提取出的頻譜水印信息與原始頻譜水印之間的相似性來衡量頻譜水印算法的魯棒性；正確提取概率用來表示正確提取頻譜水印的可能性，定義為在滿足頻譜水印的透明性和魯棒性要求情況下，在實施N次提取頻譜水印操作中，正確提取頻譜水印的次數；嵌入量用于衡量在原始通信信號中嵌入頻譜水印的強度.由于頻譜水印的長度一定，嵌入量(用嵌入系數表示，嵌入系數越大，嵌入量越大)直接表現為頻譜水印周期嵌入的次數，若嵌入次數越多，即嵌入量越大，越有利于頻譜水印的檢測與提取，但通信信號的透明性越差，反之，則嵌入量越小，通信信號的透明性越好，但越不利于頻譜水印的檢測與提取.

3.1頻譜水印檢測性能分析

實驗中使用的通信信號幀結構由同步幀、有用信息幀和截止幀組成，頻譜水印嵌入位置標識符在標識符嵌入控制模塊控制下周期性地嵌入到通信信號截止幀和同步幀之間，其比特速率為10kbps，小波分解時采樣率為2MHz，用長度為96bits的二進制序列表示頻譜水印，利用長度為P的m序列對其進行擴頻調制.頻譜水印嵌入位置標識符采用長度為16bits的m序列.圖3給出了m序列長度為63bits，頻譜水印嵌入系數為0.5，收到5段中頻信號，每段中頻信號持續時間為0.1s的某段中頻信號嵌入頻譜水印前后的時域比較圖.

(a) 嵌入頻譜水印前的中頻信號

(b) 嵌入頻譜水印后的中頻信號

(c) 嵌入頻譜水印前后中的中頻信號幅度差值圖3　某段中頻信號嵌入頻譜水印信號前后的時域比較圖

通過對比嵌入水印前后中頻信號波形圖可知，嵌入頻譜水印的中頻信號在幅度上有細微差別，但并不會改變其主要特征和分布趨勢，且根據計算可知，嵌入水印前后的信噪比為22.9dB，對原通信信號的解調影響不大.

3.2頻譜水印嵌入位置標識符檢測性能分析

為檢驗頻譜水印嵌入位置標識符檢測性能，圖4給出了頻譜水印嵌入標識符定位性能曲線.由圖4(a)可知，嵌入系數對頻譜水印嵌入標識符定位性能幾乎沒有影響.而從圖4(b)可知，當嵌入系數一定時，信噪比對其影響較大，在信噪比大于6dB時，頻譜水印嵌入標識符位置檢測錯誤率就可低于10-3.圖4(c)給出了頻譜水印嵌入標識符定位精度(用頻譜水印嵌入標識符檢測獲得的位置偏離真實位置的符號數)對頻譜水印提取正確率影響性能曲線，由大量實驗分析可知，當頻譜水印嵌入標識符定位精度小于等于1個碼元長度，提出的算法能夠正確地提取頻譜水印，反之則難以正確提取頻譜水印.

(a) 信噪比為10 dB時，嵌入系數對標識符定位錯誤率的影響

(b) 嵌入系數為0.5，信噪比對標識符定位錯誤率的影響

(c) 信噪比為10 dB時，嵌入系數為0.5，標識符定位精度對頻譜水印提取正確率的影響圖4　相關標識符定位性能曲線

圖5　中頻信號分段長度對嵌入水印前后的信噪比

3.3頻譜水印提取性能分析

為分析嵌入系數對頻譜水印算法透明性的影響，針對分段信號的不同持續時間，仿真分析了頻譜水印算法的透明性，如圖5所示.圖5中嵌入品質表示嵌入頻譜水印前后的信噪比，m序列長度為63bits，針對5段中頻信號實施頻譜水印嵌入操作(即累積5次)，為確保頻譜水印能完全嵌入到每段中頻信號中去，其持續時間必須大于0.05s.由圖5可知，隨著每段中頻信號長度增加，其嵌入品質逐漸增大，增加幅度隨著每段中頻信號持續時間的增加而降低，當每段中頻信號持續時間大于某一值后，嵌入品質增加越來越不明顯.同時，頻譜水印嵌入系數對嵌入品質影響也較大，如當每段中頻信號持續時間為0.4s時，若嵌入系數為0.4，則嵌入品質為31.477 5dB；若嵌入系數為0.5，則嵌入品質為28.9dB；若嵌入系數為0.6，則嵌入品質為27.382 4dB.

為分析頻譜水印提取的性能，圖6、7給出了頻譜水印提取的性能曲線.其中，圖6中將提取的頻譜水印與原始頻譜水印間的相似系數為1的情況視為接收端能夠正確提取頻譜水印，每段中頻信號持續時間為0.1s.

由圖6(a)可知，當利用5段中頻信號累積來提取頻譜水印時，頻譜水印嵌入系數大于等于0.4時，能完全正確地提取水印，頻譜水印嵌入系數小于等于0.3時，幾乎完全不能正確提取水印；當利用10段中頻信號累積來提取頻譜水印時，頻譜水印嵌入系數大于等于0.4時，能完全正確提取水印，頻譜水印嵌入系數大于等于0.35時，其提取正確率可達98.2%，頻譜水印嵌入系數小于等于0.3時，幾乎完全不能正確提取水印；當利用15段中頻信號累積來提取頻譜水印時，頻譜水印嵌入系數大于等于0.3時，能完全正確提取水印，頻譜水印嵌入系數小于等于0.2時，幾乎完全不能正確提取水印.因此，嵌入系數對頻譜水印提取正確率的影響大，而累積次數對頻譜水印提取正確率有一定影響，但不是很明顯.

由圖6(b)可知，當嵌入系數為0.5時，嵌入水印前后信噪比為22.904dB.由圖可以看出，在累積次數一定時，隨著信噪比的提高，頻譜水印提取正確率逐步提高，當信噪比高于一定值時，其正確率可達到100%.如累積5次，信噪比大于等于8dB時，其頻譜水印提取正確率為100%；累積10次，信噪比大于等于7dB時，其頻譜水印提取正確率為98.5%；累積15次，信噪比大于等于7dB時，其頻譜水印提取正確率為100%.因此，在嵌入系數一定時，頻譜水印提取正確率不但與信噪比有關，還與累積次數有關，當累積次數和信噪比達到某一個值后，其頻譜水印提取正確率可達100%.

(a) 信噪比為10 dB

(b) 嵌入系數為0.5圖6　頻譜水印提取正確率性能曲線

圖7(a)給出嵌入系數為0.5，每段中頻信號持續時間為0.1s，頻譜水印采用不同長度的擴頻碼情況下，信噪比對頻譜水印提取正確率的影響曲線.由圖可知:當頻譜水印采用的擴頻碼長度為63，信噪比大于等于8dB時，能夠完全正確提取水印信息；當頻譜水印采用的擴頻碼長度為127，信噪比大于等于7dB時，能夠完全正確提取水印信息.因此，頻譜水印采用的擴頻碼長度對頻譜水印提取有一定的影響，擴頻碼長度越長，其頻譜水印提取算法越能在更低的信噪比條件下正確提取頻譜水印.

(a) 嵌入系數為0.5

(b) 嵌入系數為0.5圖7　頻譜水印提取正確率隨信噪比變化曲線

圖7(b)給出了在不同相似系數要求下，信噪比對頻譜水印提取正確率的影響曲線，圖中嵌入系數為0.5，每段中頻信號持續時間為0.1s，累積次數為5次.若相似度要求不高時，頻譜水印提取算法能夠在信噪比較低的情況下，提取出頻譜水印.如相似系數大于0.7時，可在信噪比為5dB以上的噪聲環境中提取頻譜水印；相似系數大于0.9時，可在信噪比為6dB以上的噪聲環境中提取頻譜水印.

4結論

針對無線電設備頻譜身份識別問題，研究提出了一種基于相關標識符的頻譜水印嵌入與提取方法.理論分析和仿真實驗證明:提出的頻譜水印嵌入與提取方法的有效性，能夠在不影響原始通信信號的正確解調的基礎上，正確地提取頻譜水印，即在嵌入系數為0.5條件下，提出方法能夠保證其嵌入頻譜水印后的信噪比大于20dB，頻譜水印的嵌入不會對原始信號解調造成有明顯影響，頻譜水印嵌入透明性好；在信道信噪比大于7dB時，頻譜水印嵌入標識符定位錯誤率小于0.000 3，頻譜水印提取正確率能夠達到100%，頻譜水印具有一定的魯棒性.這些研究成果將為有效標識和識別無線電設備身份、快速定位無線電干擾提供技術手段.

參考文獻

[1]劉忠英， 許金勇, 趙杭生. 基于頻譜信號的頻譜水印技術研究[J]. 河北科技大學學報, 2011, 32(S1): 116-118.

LIUZY,XUJY,ZHAOHS.Spectrumwatermarkingtechniquebasedonspectrumsignal[J].JournalofHebeiUniversityofScienceandTechnology, 2011, 32(S1):116-118. (inChinese).

[2] 季秀蘭. 基于離散余弦變換和支持向量機的數字圖像水印算法[J]. 科技通報, 2012, 28(11):153-182.

JIXL.Digitalimagewatermarkingalgorithmbasedondiscretecosinetransformandsupportvectormachines[J].Bulletinofscienceandtechnology, 2012, 28(11):153-182. (inChinese).

[3]TRIVEDINK,TRIVEDIMC.Analgorithmicdigitalaudiowatermarkinginperceptualdomainusingdirectsequencespreadspectrum[C]//2014FourthInternationalConferenceonCommunicationSystemsandNetworkTechnologies.Bhopal, 2014: 805-809.

[4]KHAREP,VERMAAK,SRIVASTAVAVK.Digitalimagewatermarkingschemeinwaveletdomainusingchaoticencryption[C]//StudentsConferenceonEngineeringandSystems.Allahabad, 2014:1-4.

[5]ANIYUANA,DEEPAJ.Hardwareimplementationofarobustwatermarkingtechniquefordigitalimages[C]//IEEERecentAdvancesinIntelligentComputationalSystems.Trivandrum, 2013: 293-298.

[6]TRABELSIW,SELMIMH.Multi-signaturerobustvideowatermarking[C]//1stInternationalConferenceonAdvancedTechnologiesforSignalandImageProcessing.Sousse, 2014: 158-163.

[7] 趙軍, 余松, 葛云露, 等. 射頻水印系統中的同步算法[J]. 電波科學學報, 2012, 27(3): 465-469.

ZHAOJ,YUS,GEYL,etal.WatermarkingsynchronizationinRFwatermarkingsystem[J].Chinesejournalofradioscience, 2012, 27(3): 465-469. (inChinese)

[8] 劉奕彤. 一種為無線電臺施加水印的方法: 中國, 201110121146.8[P]. 2011-05-11.

[9] 劉芳. 基于離散小波變換的音頻數字水印研究[D]. 北京: 北京化工大學, 2010.

LIUF.Researchonaudiodigitalwatermarkingbasedondiscretewavelettransform[D].Beijing:BeijingUniversityofChemicalTechnology, 2010. (inChinese).

[10]吳彥華, 鐘子發. 基于小波變換的K-L變換去噪應用研究[J]. 電波科學學報, 2005, 20(6):745-748+757.

WUYH,ZHONGZF.ApplicationofK-Ltransformbasedonwavelettransforminreductionofnoise[J].Chinesejournalofradioscience, 2005, 20(6):745-748+757. (inChinese).

[11]肖艷軍, 李建勛. 基于小波變換的信號濾波在探地雷達中的應用[J]. 電波科學學報, 2006, 21(1): 140-145+149.

XIAOYJ,LIJX.Signalfilteringbasedonwavelettransformanditsapplicationingroundpenetratingradar[J].Chinesejournalofradioscience, 2006, 21(1): 140-145+149. (inChinese).

[12]張文婕. 小波域自適應音頻水印算法研究[D]. 成都: 西南交通大學, 2010.

ZHANGWJ.Studyonadaptiveaudiowatermarkingalgorithmsindomainwavelet[D].Chengdu:SouthwestJiaotongUniversity, 2010. (inChinese).

[13]夏嵬. 抗幾何攻擊的魯棒數字圖像水印算法研究[D]. 武漢: 華中科技大學, 2010.

XIAW.Researchongeometricallyrobustdigitalimagewatermarkingalgorithm[D].Wuhan:HuazhongUniversityofScienceandTechnology, 2010. (inChinese).

[14]黃兒松, 劉金華, 文汝紅. 基于小波視覺模型的乘性水印算法[J]. 計算機應用, 2011, 31(8):2165-2173.

HUANGES,LIUJH,WENRH.Multiplicativewatermarkingalgorithmbasedonwaveletvisualmodel[J].Journalofcomputerapplications, 2011, 31(8): 2165-2173. (inChinese).

張余(1983-),男,四川人,南京電訊技術研究所工程師,碩士,主要從事電磁頻譜技術方向的研究.

許金勇(1976-),男,江西人,南京電訊技術研究所工程師,博士,主要從事電磁頻譜技術方向的研究.

柳永祥(1974-),男,湖北人,南京電訊技術研究所高級工程師,碩士,主要從事電磁頻譜技術方向的研究.

作者簡介

中圖分類號TN014

文獻標志碼A

文章編號1005-0388(2016)01-0185-08

收稿日期：2015-03-17

張余, 許金勇, 柳永祥, 等. 基于相關標識符的頻譜水印嵌入與提取方法[J]. 電波科學學報,2016,31(1):185-192. DOI: 10.13443/j.cjors.2015031701

ZHANG Y, XU J Y, LIU Y X, et al. Spectrum watermark embedding and extracting method based on correlation identifier [J]. Chinese journal of radio science,2016,31(1):185-192. (in Chinese). DOI: 10.13443/j.cjors.2015031701

資助項目: 國家自然科學基金(No.61102092；61471395)

聯系人： 張余 E-mail: zhyu63@163.com