detect.m

改程序实现了数字水印的功能

% function s=detect(x,xx,alpha,N,w)
% 提取水印，并检测与真实水印的相关性，返回相关系数
% x为原始音频向量，xx为嵌入水印的音频向量、alpha为水印幅度衰减系数、
% N为向量长度，w为真实水印向量
  N=400
  alpha=0.2
  load('G:\matlab\work\DCT+ss\wmk400.mat')
   x=wavread('01')
   y=wavread('01wA')
 
 DCTI=dct(x)
 yy=y
 DCTJ=dct(yy)
 [Xdown,ind]=sort(abs(DCTI),'descend')
 Ind=ind(1:N)

 v=randn(N,1)
  for i=1:N
           v(i)=(DCTJ(Ind(i))/DCTI(Ind(i))-1)/alpha
           v(i)=v(i)+4.2016
  end
  meanv=mean(v)
  
    
  s=v'*w/sqrt(v'*v)