Stanford 机器学习 Week8 作业:K-means Clustering and Principal Component Analysis

for i = 1:size(X,1)
dis = sum((centroids - X(i,:)) .^ 2, 2);
[t, idx(i)] = min(dis);
end

for i = 1:K
id = find(idx == i);
tot = X(id,:);
centroids(i,:) = mean(tot);
end

id = randperm(size(X,1));
id = id(1:K);
centroids = X(id);

sigma = 1/m * X' * X;
[U,S,V] = svd(sigma);

Z = X * U(:,1:K);

X_rec = Z * U(:,1:K)'