用Matlab实现的数字音频水印嵌入与提取

用Matlab实现的数字音频水印嵌入与提取

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

% 读取音频信号3.wav存到变量A中
[Afsnbits] = wavread（‘3.wav‘）;
%绘制原始声音图像 
subplot（211）; plot（A）; 
axis（[0 350000 -2 2]）;
title（‘原始声音信号波形图‘）;
% 用变量L存储音频A的长度
L = size（A）;
% 读取图像Lena.bmp存入变量M做水印信号
M = imread（‘Lena.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 < th ）
    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‘test.wav‘）;
subplot（212）; plot（G）; 
title（‘含水印的声音信号波形图‘）;
axis（[0 350000 -2 2]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2261  2009-01-09 11:44  音频水印嵌入与提取\embed.m

     文件       2046  2009-01-06 23:50  音频水印嵌入与提取\extract.m

     文件      17838  2008-12-19 22:09  音频水印嵌入与提取\Lena.bmp

     文件       3807  2009-01-09 11:22  音频水印嵌入与提取\watermark.m

     文件       2361  2009-01-09 11:25  音频水印嵌入与提取\watermark1.m

     目录          0  2009-03-11 11:24  音频水印嵌入与提取

----------- ---------  ---------- -----  ----

                28313                    6