matlab音频水印算法研究+代码

%加入白噪声的音频水印程序
clear;
[y ,fs] = audioread('mei.wav');   %读入原始音频文件
[c,l]=wavedec(y,3,'db4');               %三级小波分解
ca3=appcoef(c,l,'db4',3);%提取3级近似系数，绘制原始信号和近似系数。
cd3=detcoef(c,l,3);%提取 3 所指定的级别上的细节系数
cd2=detcoef(c,l,2);%提取 2 所指定的级别上的细节系数
cd1=detcoef(c,l,1);%提取 1 所指定的级别上的细系数节
x=ca3;    %提取近似低频系数
len=length(y);  %计算音频数据长度
x1=x;%x重新赋值给x1
s=max(abs(x))*0.2;%对近似系数的绝对值，求其中的最大值，再乘以0.2
i=find(abs(x)>s);lx=length(x(i));         %找出大于最大值0.2倍的序列
figure;%创建图窗窗口
subplot(2,2,1);%将当前图窗划分为 2×2 网格，并在1指定的位置创建坐标区
plot(ca3);                            %画出低频系数图
title('低频系数图形');%标题
subplot(2,2,2);%将当前图窗划分为 2×2 网格，并在2指定的位置创建坐标区
plot(cd3);%画出第3级的细节系数。
title('cd3');%标题
subplot(2,2,3);%将当前图窗划分为 2×2 网格，并在3指定的位置创建坐标区
plot(cd2);%画出第2级的细节系数。
axis([0 10e4 -0.5 0.5]);%设置坐标轴范围和纵横比
title('cd2');%标题cd2
subplot(2,2,4);%将当前图窗划分为 2×2 网格，并在3指定的位置创建坐标区
plot(cd1);%画出cd1
title('cd1');%标题cd1
randn('seed',10);      %产生随机高斯序列
mark=randn(1,lx);%创建一个由随机数组成的 1×lx 向量。
ss=mark;%mark赋值给ss
rr=ss*0.1;                       %设置水印嵌入强度
x(i)=x(i).*(1+2*rr');                   %嵌入水印
c1=[x',cd3',cd2',cd1'];  %创建矩阵c1
s1=waverec(c1,l,'db4');%重建信号
file1='已加水印.wav';
dd=length(s1);                     %调整s1的长度,使之可以分成两列
s11=s1;                            %s1的值不能改变,因为后面还需要用到
if rem(dd,2)==1                    %如果s1是奇数,则去掉最后一个数,将新数组定义为s11
    s11=s1(1:dd-1);               
end
ee=reshape(s11,[],2);                %将s1调整成2列的数组
audiowrite(file1,y,fs);%将加水印的音频写入file1.wav
figure;%创建图窗
subplot(3,1,1);%创建3*1的网格，把该图放在第一个位置
plot(y);                   %画出原信号图
axis([0 18e4 -2 2]);%设置坐标轴范围和纵横比
title('原信号的图');     %标题               
subplot(3,1,2);plot(ss);    %创建3*1的网格，把该图放在第二个位置
title('水印图');
subplot(3,1,3);plot(s1);                  %画出嵌入了水印的信号图
title('加入了水印的声音信号')%加标题
kk.wave = wavread('已加水印'); 
yc=kk.wave;
dy=length(y);
if rem(dy,2)==1                         %如果y为奇数,处理同s1,但是因为后面不需要用到y,所以不必定义一个新数组
    y=y(1:dd-1);
end
y1=reshape(y,[],2);                       %调整数组y的维数,使之可以和yc做运算
fz=sum(y1.*y1);  %计算嵌入了水印的信号的信噪比,".*"用于实现对应元素的乘法

结果：

matlab群：912166339