test.m

LSK算法源码,欢迎大家交流

%function mark()
m = 256; L = 8;
mark = 13*ones(1,m*m);mark(3),mark(4)
i = m/(2^L) * m/(2^L);
e = log2(m/(2^L));
for j = 1:i
    mark(j) = 0;
end;i,j
j = 1;
while i < m*m
    %for j = 1:e
        mark(i+j) = e;
        mark(2*i+j) = e;
        mark(3*i+j) = e;
    %end;
    i = 4*i;
    e = e+1;
end;

% map = mapping(m);
% temp(map) = mark;
% mark = temp;
% Mark = ones(m,m);
% for i = 1:m
%     for j = 1:m
%         Mark(i,j) = mark((i-1)*m+j);
%     end;
% end;
% 
% Mark