fsclcode.m

来自「这里包含了聚类的工具箱还有很详细的文档说明」· M 代码 · 共 36 行

function [centro,nvoro] = fsclcode(X,centri,nvori,iter)

[N,d] = size(X);
[M,dc] = size(centri);

if dc ~= d
    disp('Dimensions from data and centroids is not equal?!')
    return;
end

% d = data dimension
% N = number of points

centro = centri;
nvoro = zeros(size(nvori));

alph = 0.5/iter;


[dummy,pos] = sort(rand(1,N));
%pos = 1:N;

for n = 1:N

    datum = X(pos(n),:);

    nvor = max([nvoro nvori],[],2);

    diff = nvor.*sqrt(sum((centro - ones(M,1)*datum).^2,2));
    minpos = find(diff == min(diff));
    minpos = minpos(1); % case there are two centroids, we choose the first one
    nvoro(minpos) = nvoro(minpos) + 1;

    centro(minpos,:) = (1-alph)*centro(minpos,:) + alph*datum;
end