% function s=detect(x,xx,alpha,N,w)
% 提取水印，并檢測(cè)與真實(shí)水印的相關(guān)性，返回相關(guān)系數(shù)
% x為原始音頻向量，xx為嵌入水印的音頻向量、alpha為水印幅度衰減系數(shù)、
% N為向量長(zhǎng)度，w為真實(shí)水印向量
  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)