Matlab有关DCT域音频水印嵌入程序-31767689audiowatermark.rar

Matlab有关DCT域音频水印嵌入程序-31767689audiowatermark.rar
DCT域的音频水印嵌入的程序
Figure35.jpg Matlab有关DCT域音频水印嵌入程序

%水印嵌入*********************************************************************************************（1）
%1同步信号，在音频起始点嵌入10个幅值相同的点
[xfsnbits]=wavread（‘E:\wave1\0000.wav‘）;
lenx=length（x）;
w=wavread（‘E:\wave1\0001.wav‘）;
y=[0.02; 0.02 ;0.02; 0.02 ;0.02; 0.02 ;0.02; 0.02 ;0.02 ;0.02];
%input m
m=1000;
k=fix（lenx/10/m）*10;
w1=w（1:m）;
%input a
a=50;
%2分段处理
for i=1:k:m*k
   %3离散余弦变换
    x1=dct（x（i:i+k-1））;
   %排序找出中频系数mid 用公式mid‘=mid（1+a*w）
   [x11index]=sort（x1）;
   j=fix（i/k）+1;   
   for t=1:3
      %4修改中频系数
      x1（index（t+1））=x11（t+1）*（1+a*w（j））;
      %纪录被改变的系数和序号
      recordx（jt）=x11（t+1）;
      recordi（jt）=index（t+1）;
      %5离散余弦逆变换
      x（i:i+k-1）=idct（x1）;
   end
end
%6合并得到嵌入水印后的音频信号
x=[y;x];
%将嵌入水印后的音频信号保存
wavwrite（xfsnbits‘E:\wave1\00000.wav‘）;
%水印提取*********************************************************************************************（2）
%1同步检测
[x2fsnbits]=wavread（‘E:\wave1\00000.wav‘）;
lenxx=length（x2）-m*k;
i=1;
%input b
b=0.001;
while and（or（or（abs（x2（i）-0.02）>babs（x2（i+1）-0.02）>b）abs（x2（i+2）-0.02）>b）i    i=i+1;
end
while  and（abs（x2（i+3）-0.02）   i=i+1;
end
%判断是否有同步检测信号
if i>lenxx
  error（‘没有同步检测信号 ‘）;
end
%2含有水印的音频信号做分段处理
x3=x2（（i+3）:（i+2+m*k））;
length（x3）
%3原始音频信号做分段处理
for i=1:k:m*k  
%4含有水印的音频信号做离散余弦变换
   x4=dct（x3（i:i+k-1））;
   j=fix（i/k）+1;
%5水印提取
   w21=（x4（recordi（j1））/recordx（j1）-1）/a;
   w22=（x4（recordi（j2））/recordx（j2）-1）/a;
   w23=（x4（recordi（j3））/recordx（j3）-1）/a;
   w2（j）=（w21+w22+w23）/3;
end
%6求相关系数
cox1=（w1‘*w2‘）/（norm（w1）*norm（w2））
%水印攻击*********************************************************************************************（3）
%1加噪音，    转（2）
%2滤波   ，  转（2）
%3剪切，     转（2）
%4重量化，   转（2）
%5重采样，   转（2）
%6MPEG压缩，转（2）




   
   

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

     文件       2033  2002-12-01 14:13  31767689audiowatermark\watermark.m

     文件       3543  2009-02-01 16:19  31767689audiowatermark\使用帮助：新手必看.htm

     文件        183  2009-02-02 09:47  31767689audiowatermark\Matlab中文论坛--助努力的人完成毕业设计.url

     目录          0  2009-04-16 17:36  31767689audiowatermark

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

                 5977                    5