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