clear all;
load codebook16_1024;
k=256;
N=1024;
%哈德碼變換
HY=hadamard(k);
Y=[];
y=[];
%碼字
for j=1:N
    for i=1:k
        y(j,i)=I12(j,i);
    end
end
y=y';
Y1=HY*y;
Y=Y1;

%按各碼字哈德碼變換的第一維分量大小對碼書進行升序排列
for i=1:N-1
    for j=i+1:N
        if Y(1,j)<Y(1,i)
            for s=1:k
                t=Y(s,j);
                Y(s,j)=Y(s,i);
                Y(s,i)=t;
            end
        end
    end
end
Q=[];
for i=1:N
    q=i;
    for j=1:N
        if Y(1:k,i)==Y1(1:k,j)
            q=j;
            break;
        end
    end
    Q=[Q,q];
end

U=[];
R=Q;
for i=1:N
    U=[U,Y1(1:k,R(i))];
end


        
%Q=[];
%for i=1:N-1
%    q=i;
%    for j=i+1:N
%        if Y(1,j)<Y(1,q)
%            q=j;
%            for s=1:k
%                t=Y(s,q);
%                Y(s,q)=Y(s,i);
%                Y(s,i)=t;
%            end
%        end
%    end
%    Q=[Q,q];
%end
%Q=[];
%for i=1:N-1
%    q=i;
%    for j=i+1:N
%        if Y(1,j)<Y(1,i)
%            q=j;
%            for s=1:k
%                t=Y(s,j);
%                Y(s,j)=Y(s,i);
%                Y(s,i)=t;
%            end
%        end
%    end
%    Q=[Q,q];
%end