codeword8_512.m

基于等均值等方差的图象矢量量化Matlab程序

clear all;
load codebook8_512;
k=64;
N=512;

%码字
for j=1:N
    for i=1:k
        y1(j,i)=I13(j,i);
    end
end
y1=y1';
y=y1;

%码字均值
mj1=[];
for j=1:N
    sum1=0;
    for i=1:k
        sum1=sum1+y(i,j);
    end
    average=sum1/k;
    mj1=[mj1,average];
end
mj=mj1;

%按均值大小对码书进行升序排列
for i=1:N-1
    for j=i+1:N
        if mj(j)<mj(i)
           t=mj(j);
           mj(j)=mj(i);
           mj(i)=t;            
        end
    end
end
Q=[];
for i=1:N
    q=i;
    for j=1:N
        if mj(i)==mj1(j)
            q=j;
            break;
        end
    end
    Q=[Q,q];
end
for i=1:N
    y(1:k,i)=y1(1:k,Q(i));
end
%码字方差
vj=[];
for j=1:N
    sum2=0;
    for i=1:k
        sum2=sum2+(y(i,j)-mj(j))^2;
    end
    v1=sqrt(sum2);
    vj=[vj,v1];
end

U=[];
R=Q;
for i=1:N
    U=[U,mj1(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