本帖最后由 馨菲smile 于 2015-7-9 16:43 编辑
clear
filedir1 = 'F:\huhu\01\';
for Num1 = 1:1:116
filename1 = strcat(filedir1,num2str(Num1),'-1.mat');
A =[];
t = [];
A = load (filename1);
X(Num1,:)=A.AB(:,2);
end
t(:,1)=A.AB(:,1);
K = 4;
dataSet = load('X');
[row,col] = size(dataSet);
% 存储质心矩阵
centSet = zeros(K,col);
% 随机初始化质心
for i= 1:col
minV = min(dataSet(:,i));
rangV = max(dataSet(:,i)) - minV;
centSet(:,i) = repmat(minV,[K,1]) + rangV*rand(K,1);
end
% 用于存储每个点被分配的cluster以及到质心的距离
clusterAssment = zeros(row,2);
clusterChange = true;
while clusterChange
clusterChange = false;
% 计算每个点应该被分配的cluster
for i = 1:row
% 这部分可能可以优化
minDist = 10000;
minIndex = 0;
for j = 1:K
distCal = distEclud(dataSet(i,:) , centSet(j,:));
if (distCal < minDist)
minDist = distCal;
minIndex = j;
end
end
if minIndex ~= clusterAssment(i,1)
clusterChange = true;
end
clusterAssment(i,1) = minIndex;
clusterAssment(i,2) = minDist;
end
% 更新每个cluster 的质心
for j = 1:K
simpleCluster = find(clusterAssment(:,1) == j);
centSet(j,:) = mean(dataSet(simpleCluster',:));
end
end
figure
%scatter(dataSet(:,1),dataSet(:,2),5)
for i = 1:K
pointCluster = find(clusterAssment(:,1) == i);
scatter(dataSet(pointCluster,1),dataSet(pointCluster,2),5)
hold on
end
%hold on
scatter(centSet(:,1),centSet(:,2),300,'+')
hold off
end
% 计算欧式距离
function dist = distEclud(vecA,vecB)
dist = sqrt(sum(power((vecA-vecB),2)));
end
我的01文件夹里是.m文件,运行总会显示‘Undefined function or variable 'K'.’
不知道是哪里出了错,请哪位老师指导下。多谢多谢!