clc;clear;close?all;
%?水印嵌入程序

%?讀取音頻信號audio.wav存到變量A中
[Afsnbits]?=?wavread（‘login.wav‘）;
%繪制原始聲音圖像?
subplot（211）;?plot（A）;?
axis（[0?350000?-2?2]）;
title（‘原始聲音信號波形圖‘）;
%?用變量L存儲音頻A的長度
L?=?size（A）;
%?讀取圖像Lena.bmp存入變量M做水印信號
M?=?imread（‘mark.bmp‘）;
%?將圖像Lena.bmp轉化為二值圖并存入變量BW
BW?=?im2bw（M）;

%?計算水印矩陣大小
[M1M2]?=?size（BW）;
%?M12為中間變量，避免每次都計算M1*M2
M12?=?M1*M2;?
%?降維，將水印信息得到的一維序列存入序列C中
C?=?reshape（BW1M12）;
n?=?M12;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?對水印信號進行擴頻處理效果不是很好
%?擴頻系數為2
n?=?M12*2;
M?=?zeros（n1）;
%?產生密鑰序列M
for?k?=?1?:?n
????if?mod（k4）?==?0
????????M（k）?=?1;
????else?
????????M（k）?=?0;
????end
????%?水印信號序列分別按位與密鑰異或
????l?=?ceil（k/2）;
????S（k）?=?bitxor（C（l）M（k））;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%?嵌入一個水印信息需要的音頻數據為N
N?=?10;?
length?=?n*10;
%?將原始音頻信號分解為Ae和Ar兩部分
i?=?1?:?length;?
j?=?[1];
%?取矩陣A的l到length行構建矩陣Ae
Ae?=?A（ij）;
%?Ae（ij）
i?=?length+1?:?L;
%?取矩陣A的length到L行第一列構建矩陣Ar
Ar?=?A（ij）;


%?建立元胞B，每個音頻數據段Ae（m）是B的一個元素
k?=?1;
%?建立M1?x?M2行l列的元胞
B?=?cell（n1）;
th?=?n*N;
%?當k小于Ae的長度時，矩陣Ae每10行作為一個音頻數據段存入元胞B中
while?（?k?????i?=?k?:?k+9;
????m?=?（k+9）/10;
????B{m1}?=?Ae（ij）;
????k?=?k+10;
end


%?建立元胞D，并將元胞B中離散余弦變換了的元素存入元胞D中
D?=?cell（n1）;
for?i?=?1?:?n
????D{i1}?=?dct（B{i1}）;
end

%?建立元胞E，儲存嵌有水印信號的中頻系數
E?=?cell（n1）;
E?=?D;
for?i?=?1?:?n
%?????%?將水印信息C嵌入音頻中
%?????E{i1}（3）?=?（D{i1}（3））*（1+2*C（i））;
????%?將擴頻了的水印信息S嵌入音頻中
????E{i1}（3）?=?（D{i1}（3））*（1+2*S（i））;
end

%?建立元胞F
F?=?cell（n1）;
%?將元胞E中離散余弦反變換了的元素存入元胞F中
for?i?=?1?:?n
????F{i1}?=?idct（E{i1}）;
end

%?將所有分離的段合并創建到一維矩陣G中
G?=?F{11};
%?i從3到M1*M2依次取值
for?i?=?2?:?n
????%?將元胞G和元胞F第i組元素合并創建到一維矩陣G
????G=[G;?F{i1}];
end
%?將元胞G和矩陣Ar合并創建一維矩陣G
G?=[G;?Ar];
wavwrite（Gfsnbits‘mark_audio.wav‘）;
subplot（212）;?plot（G）;?
title（‘含水印的聲音信號波形圖‘）;
axis（[0?350000?-2?2]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2275??2010-05-09?14:36??音頻\1embed.m

?????文件???????2056??2010-05-09?14:36??音頻\1extract.m

?????文件???????2965??2010-05-09?14:56??音頻\2dwt?audio.m

?????文件???????1284??2010-05-09?14:58??音頻\3DCT_Audio.M

?????文件???????1187??2010-05-09?16:13??音頻\4dct+SS_audio.m

?????文件???????2365??2010-05-09?16:15??音頻\6audio_watermark1.m

?????文件???????2265??2010-05-09?16:16??音頻\7音頻水印嵌入與提取\embed.m

?????文件???????2046??2009-01-06?23:50??音頻\7音頻水印嵌入與提取\extract.m

?????文件??????17838??2008-12-19?22:09??音頻\7音頻水印嵌入與提取\Lena.bmp

?????文件??????88266??2010-04-16?11:19??音頻\7音頻水印嵌入與提取\login.wav

?????文件???????5174??2010-05-08?19:12??音頻\7音頻水印嵌入與提取\mark.bmp

?????文件?????????44??2010-05-09?16:17??音頻\7音頻水印嵌入與提取\test.wav

?????文件???????3793??2010-05-09?16:18??音頻\7音頻水印嵌入與提取\watermark.m

?????文件???????2361??2009-01-09?11:25??音頻\7音頻水印嵌入與提取\watermark1.m

?????文件???????5003??2009-05-09?09:32??音頻\8DWT.txt

?????文件???????1841??2010-05-08?23:01??音頻\9音頻水印程序.txt

?????文件??????88266??2010-04-16?11:19??音頻\login.wav

?????文件???????5174??2010-05-08?19:12??音頻\mark.bmp

?????文件??????88252??2010-05-09?14:55??音頻\modified?login.wav

?????文件???????2848??2010-05-09?16:41??音頻\復件?embed.m

?????文件???????2056??2010-05-09?16:50??音頻\復件?extract.m

?????目錄??????????0??2010-05-09?16:17??音頻\7音頻水印嵌入與提取

?????目錄??????????0??2010-05-19?10:17??音頻

-----------?---------??----------?-----??----

???????????????327359????????????????????23