even_k.m

来自「神经网络VC++代码 人工神经网络原理及仿真实例.」· M 代码 · 共 32 行

function c=even_k(x,c)
[np,ni]=size(x);
[ni,m]=size(c);
y=zeros(1,np);
q=zeros(1,m);

for i=1:np
   p=x(i,:)'*ones(1,m);
   q=sum((p-c).^2);
   j=find(q==min(q)); 
   [l,n]=size(j);
   if n==1     
      y(i)=j;
   else
      y(i)=min(j);
   end
end   %确定x(i,:) (i =1,2,…,np)属于哪一类

for j=1:m
   z=find(y==j);
   [l,n]=size(z);
   if n==1
      c(:,j)=x(z,:)';
   else
         u=[];
        for i=1:n
            u=[u;x(z(i),:)];
        end
         c(:,j)=(sum(u)/n)';  
     end
  end         %确定经过聚类后新的中心