tiqu_3.m

来自「一些程序」· M 代码 · 共 29 行

%从小波系数中提取水印
function a=tiqu_3(x1,x2,N)
%x1为原始的未嵌水印的小波系数
%x2为已嵌入水印的小波系数
%N为嵌入位数
%s为量化系数
x1_length=length(x1);
%找出从大到小第N大的值
xx=abs(x1);
% a=ones(1,N);
for i=1:N        
    j=x1_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(x1)>=xxxx);

for i=1:N
    a(i)=(x1(d(i))<x2(d(i)));         
end