kmeans_evaluation1.m

来自「kmeans算法源程序 共6个文件 很好用的」· M 代码 · 共 28 行

function [J]= kmeans_evaluation1(y, class_label, k, c);
%
%
%
[nr, nc] = size(y);

for i = 1:1:k
    t = 1;
    for j = 1:1:nr
        if i == class_label(j)
            x(t,1:c) = y(j,1:c);
            t = t+1;
        end
    end
    
    [xr,xc] = size(x);
    d(i) = 0;
    for j = 1:1:xr-1
        for p = j+1:1:xr
            d(i) = d(i) + (x(j,:)-x(p,:))*(x(j,:)-x(p,:))';
        end
    end
    clear x
end
J = sum(d)/nr;