%水印嵌入*********************************************************************************************（1）
%1同步信號，在音頻起始點嵌入10個幅值相同的點
[x,fs,nbits]=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）
   [x11,index]=sort(x1);
   j=fix(i/k)+1;   
   for t=1:3
      %4修改中頻系數
      x1(index(t+1))=x11(t+1)*(1+a*w(j));
      %紀錄被改變的系數和序號
      recordx(j,t)=x11(t+1);
      recordi(j,t)=index(t+1);
      %5離散余弦逆變換
      x(i:i+k-1)=idct(x1);
   end
end
%6合并得到嵌入水印后的音頻信號
x=[y;x];
%將嵌入水印后的音頻信號保存
wavwrite(x,fs,nbits,'E:\wave1\00000.wav');
%水印提取*********************************************************************************************（2）
%1同步檢測
[x2,fs,nbits]=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)>b,abs(x2(i+1)-0.02)>b),abs(x2(i+2)-0.02)>b),i<lenxx)
    i=i+1;
end
while  and(abs(x2(i+3)-0.02)<b, i<=lenxx)
   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(j,1))/recordx(j,1)-1)/a;
   w22=(x4(recordi(j,2))/recordx(j,2)-1)/a;
   w23=(x4(recordi(j,3))/recordx(j,3)-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）