Matlab實(shí)現(xiàn)簡單擴(kuò)頻語音水印算法詳解
一、實(shí)驗(yàn)背景
1.實(shí)驗(yàn)?zāi)康?/h3>
了解擴(kuò)頻通信原理，掌握擴(kuò)頻水印算法的基本原理，設(shè)計(jì)并實(shí)現(xiàn)一種基于音頻的擴(kuò)頻水印算法，了解參數(shù)對(duì)擴(kuò)頻水印算法性能的影響。
2.實(shí)驗(yàn)環(huán)境
(1) Windows 11 操作系統(tǒng)；
(2) Matlab R2020b 科學(xué)計(jì)算軟件；
(3) WAV音頻文件。
3.原理簡介
①擴(kuò)頻基本原理
擴(kuò)頻是一種能在高噪聲環(huán)境下可靠傳輸數(shù)據(jù)的重要通信技術(shù)，其基本原理是：信號(hào)在大于所需的帶寬內(nèi)進(jìn)行傳輸，數(shù)據(jù)的帶寬擴(kuò)展是通過一個(gè)與數(shù)據(jù)獨(dú)立的碼字完成的，并且在接收端需要該碼字的一個(gè)同步接收，以進(jìn)行解擴(kuò)和數(shù)據(jù)恢復(fù)。
②擴(kuò)頻通信的特點(diǎn)
占據(jù)頻帶很寬，每個(gè)頻段上的能量很低；
即使幾個(gè)頻段的信號(hào)丟失，仍可恢復(fù)信號(hào)；
利用相互正交的擴(kuò)頻碼，可以利用這個(gè)優(yōu)點(diǎn)設(shè)計(jì)水印算法。
③實(shí)驗(yàn)算法
本例中設(shè)計(jì)一種簡單的算法：利用正交的PN序列代表0、1信號(hào)，并將其疊加到信號(hào)DCT域。提取水印時(shí)，利用PN序列的正交性可以較為準(zhǔn)確地恢復(fù)水印。
二、基礎(chǔ)知識(shí)
1.PN序列
PN序列(Pseudo-noise Sequence)，又稱偽噪聲序列，這類序列具有類似隨機(jī)噪聲的一些統(tǒng)計(jì)特性，但和真正的隨機(jī)信號(hào)不同，它可以重復(fù)產(chǎn)生和處理，故稱作偽隨機(jī)噪聲序列。
PN序列一般用于擴(kuò)展信號(hào)頻譜。PN序列的擴(kuò)頻是指用一個(gè) 序列去乘以一個(gè)信息符號(hào)，序列碼片的時(shí)間遠(yuǎn)小于信息符號(hào)的時(shí)間，由信號(hào)的時(shí)間與頻譜的關(guān)系，我們可以知道擴(kuò)頻后的序列的頻譜是展寬的。
由于PN序列的相關(guān)性很低，只有在發(fā)送的PN序列和接收的PN序列相同，并且其碼片同步時(shí)才能得到一個(gè)相關(guān)峰。當(dāng)發(fā)射時(shí)，信號(hào)的功率是低于噪聲的功率的，如果不知道PN序列，則較難得知碼片的信息。因此，擴(kuò)頻通信具有一定的加密性。
2.時(shí)域到頻域變換的原因
時(shí)域是信號(hào)在時(shí)間軸隨時(shí)間變化的總體概括，頻域是把時(shí)域波形的表達(dá)式做傅立葉變化，得到復(fù)頻域的表達(dá)式，所畫出的波形就是頻譜圖。是描述頻率變化和幅度變化的關(guān)系。將時(shí)域變換到頻域是為了做頻譜分析，即在時(shí)域中我們研究信號(hào)的時(shí)間特性，在頻域中我們研究信號(hào)的頻率特性，由于信號(hào)往往在頻域比在時(shí)域更加簡單、直觀，所以大部分信號(hào)分析的工作是在頻域進(jìn)行的。
3.三種時(shí)域到頻域變換的區(qū)別
FTP: 離散傅里葉變換，對(duì)于連續(xù)信號(hào)的靜態(tài)變換;
DWT:離散小波變換，對(duì)短時(shí)間(瞬間)信號(hào)的動(dòng)態(tài)變換;
DCT:離散余弦變換，對(duì)連續(xù)信號(hào)的動(dòng)態(tài)變換。
三、算法源碼
1.PN產(chǎn)生函數(shù)
function out =pn_gen(g,init,shift)
format=1;
out_len = 0;
in_len= 0;
out=[];
%check parameter fomat, ether g2 =[100000101]or g1 = [820]
tp = max(g);
if tp==1
format = 2;% format of parameter
in_len = length(g) -1;
else
fommat = 1;
in_len= g(1);
end
out_len= 2^in_len-1; % length of output
out = zeros(1, out_len);
for n = 1:out_len
out(n) = init(in_len);
if fommat==1
tp=0;
for m=2:length(g)
tp = mod((tp+init(g(m) + 1)), 2);%caculate new init(1)
tp = mod((tp + init(in_len- g(m))),2); %caculate new init(1)
end
else
tp= init .* g(2 : (in_len+1));
tp = mod(sun(tp), 2);
end
init=[tp init(1 : (in_len-1))];
end
for n = (shift- 1):-1:0
out = [out(2 : out_len),out(1)];
end
2.隱藏算法
function o = hide_ds(fragment, data, s, atten, pn0, pn1)
[row, col] = size(s);
if(row> col)
s=s';
end
i =1;
n = min( floor(length(s) / fragment), length(data));
o = s;
len = length(pn0);
base = fragment -len + 1;
for i=1 : n
st =(i- 1) *fragment+1;
ed = i*fragment;
tmp = dct(s(st:ed));
attenl = atten *max(abs(tmp));
if data(i) == 1
tmp(base:fragment)=tmp( base:fragment)+ attenl*pn1;
else
tmp(base:fragment) = tmp(base:fragment) + attenl*pn0;
end
o(st : ed) = idct( tmp);
end
3.提取算法
function out = dh_ds(fragment, in, pn0, pnl)
[row, col] = size(in);
if(row> col)
in = in';
end
i=1;
len = floor(length(in) / fragment);
out =[];
len_pn = length( pn0);% length of pn
base = fragment - len_pn+1;
for i = 1:len
st = (i-1) *fragment +1;
ed = i *fragment;
p = dct(in(st: ed));
t0 = sum( p(base:fragment).* pn0);
t1 = sum( p(base:fragment).* pnl);
if t1>t0
out(i) = 1;
else
out(i) = 0;
end
end
4.測試腳本
% 1 select cover audio
[fname, pname] = uigetfile(' *.wav','Select cover audio');
sourcename = strcat(pname, fname) ;
s = audioread( sourcename)';
s_len = length(s);
% 2 generate msg to be embedded
frag = 256;
msg_len = floor(s_len / frag);
msg = randsrc(1, msg_len, [0 1]);
% 3 generate PN
degree = 7;
pn0 = 2*pn_gen([degree 6 0], [zeros(1, degree - 1) 1],0) - 1;
pn1 = 2*pn_gen([degree 6 0],[zeros(1, degree - 1) 1],1) - 1;
%4 embed msg
atten = 0.005;
bld = hide_ds(frag, msg, s, atten, pn0, pn1);
% 5 save the stegoed-audio
audiowrite('hide.wav', bld, 8e+3);
% 6 select stegoed-audio
[fname, pname] = uigetfile('*.wav','Select stegoed-audio');
sourcename = strcat(pname, fname);
steg = audioread(sourcename)' ;
% 7 extract msg
out = dh_ds(frag, steg, pn0, pn1);
% 8 compute ebr
fid = 1;
ebr = sum( abs(msg -out)) / s_len;
fprintf(fid, 'ebr;%f\n', ebr);
四、運(yùn)行測試
1.無攻擊（誤碼率0.000976）:
① 運(yùn)行test.m,選擇提前準(zhǔn)備好的wav文件（載體音頻）：
② 得到攜帶水印的音頻hide.wave, 可利用音頻處理軟件對(duì)音頻進(jìn)行攻擊，觀察攻擊后水印的恢復(fù)情況。
③最后提取水印，并計(jì)算誤碼率
2.AU格式轉(zhuǎn)換（誤碼率0.001921）:
（wav->mp3->wav）
3.壓縮與解壓縮（誤碼率0.002029）:
到此這篇關(guān)于Matlab實(shí)現(xiàn)簡單擴(kuò)頻語音水印算法詳解的文章就介紹到這了,更多相關(guān)Matlab 擴(kuò)頻語音水印算法內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
這篇文章主要為大家詳細(xì)介紹了C++中的運(yùn)算符重載的相關(guān)資料，文中的示例代碼講解詳細(xì)，對(duì)我們學(xué)習(xí)C++有一定幫助，需要的可以參考一下2022-09-09
這篇文章主要為大家詳細(xì)介紹了C語言實(shí)現(xiàn)猜數(shù)字游戲，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2019-11-11
VS2019開發(fā)簡單的C/C++動(dòng)態(tài)鏈接庫并進(jìn)行調(diào)用的實(shí)現(xiàn)
這篇文章主要介紹了VS2019開發(fā)簡單的C/C++動(dòng)態(tài)鏈接庫并進(jìn)行調(diào)用的實(shí)現(xiàn)，文中通過示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-03-03
C++程序內(nèi)存棧區(qū)與堆區(qū)模型案例分析
一直以來總是對(duì)這個(gè)問題的認(rèn)識(shí)比較朦朧，我相信很多朋友也是這樣的，總是聽到內(nèi)存一會(huì)在棧上分配，一會(huì)又在堆上分配，那么它們之間到底是怎么的區(qū)別呢，讓我們一起來看看2022-03-03
這篇文章主要介紹了C語言實(shí)現(xiàn)的PNPoly算法代碼例子,PNPoly算法j是判斷一個(gè)坐標(biāo)點(diǎn)是否在不規(guī)則多邊形內(nèi)部的算法,需要的朋友可以參考下2014-07-07
本篇文章就來整理一下 Qt 中使用線程的五種方式，方便后期回顧。前面兩種比較簡單，一筆帶過了，主要介紹后面三種，感興趣的朋友跟隨小編一起看看吧2021-09-09
這篇文章主要介紹了C語言的字符空間與非字符空間,本文給大家介紹的非常詳細(xì)，具有參考借鑒價(jià)值，需要的朋友可以參考下,希望能給你帶來幫助2021-08-08
Matlab實(shí)現(xiàn)二維散點(diǎn)主方向直方圖的繪制詳解
這篇文章主要為大家詳細(xì)介紹了如何利用Matlab實(shí)現(xiàn)二維散點(diǎn)主方向直方圖的繪制，文中的示例代碼講解詳細(xì)，對(duì)我們學(xué)習(xí)Matlab有一定幫助，需要的可以參考一下2022-09-09
c++中關(guān)于int、long、long?long等取值范圍
這篇文章主要介紹了c++中關(guān)于int、long、long?long等取值范圍,具有很好的參考價(jià)值,希望對(duì)大家有所幫助,如有錯(cuò)誤或未考慮完全的地方,望不吝賜教2024-02-02
今天我們學(xué)習(xí)的是鏈棧，也就是說棧的鏈?zhǔn)浇Y(jié)構(gòu)，我們運(yùn)用順序鏈的方式來實(shí)現(xiàn)。首先呢，鏈棧是不存在存儲(chǔ)空間滿的情況的，所以可以說它是個(gè)無底洞，然而我們之前學(xué)的順序棧是有額定空間的2022-06-06