embed.m

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

%function   embed(x,alpha,N,w)
%在向量x中嵌入水印，并将含有水印的音频向量写入02.wav
%alpha 为幅度衰减系数、N为向量长度、w为水印向量
    N=1000
    alpha=0.2
    load('G:\matlab\work\DCT+ss\wmk.mat')
    x=wavread('04')
   
    
    DCTI=dct(x)
    [Xdown,index]=sort(abs(DCTI),'descend')
    Index=index(1:N)
   
    for i=1:N
    DCTI(Index(i))=DCTI(Index(i))*(1+alpha*w(i))
    end
    xx=idct(DCTI)
    wavwrite(xx,44100,'04w')