qianru_3.m

一些程序

function [a]=qianru_1(x,N,w,s)
%x为待嵌入水印的小波系数
%N为嵌入位数
%w为水印信息，长度为N
%s为嵌入强度系数
x_length=length(x);

%找出从大到小第N大的值
xx=abs(x);
for i=1:N        
    j=x_length;
    while j>=i+1
        if xx(j)>xx(j-1)
            xxx=xx(j-1);
            xx(j-1)=xx(j);
            xx(j)=xxx;
        end
        j=j-1;       
    end
end
xxxx=xx(j);

d=find(abs(x)>=xxxx);

i=1;
a=x;

while i<=N
    if w(i)==0
        a(d(i))=((a(d(i)))-s);
    else
        a(d(i))=((a(d(i)))+s);
    end
    i=i+1;
end