Files
2026-07-13 13:30:25 +08:00

16 lines
321 B
Matlab
Executable File

function mu = kseeds(X, k)
% Perform kmeans++ seeding
% Input:
% X: d x n data matrix
% k: number of seeds
% Output:
% mu: d x k seeds
% Written by Mo Chen (sth4nth@gmail.com).
n = size(X,2);
D = inf(1,n);
mu = X(:,ceil(n*rand));
for i = 2:k
D = min(D,sum((X-mu(:,i-1)).^2,1));
mu(:,i) = X(:,randp(D));
end