www.pudn.com > dct+ss.rar > 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)