vqdataplot.m

来自「一个关于数据聚类和模式识别的程序,在生物化学,化学中因该都可以用到.希望对大家有」· M 代码 · 共 24 行

function vqDataPlot(data, center)
% vqDataPlot: Plot the result of vector quantization (used in kmeans.m and vqLBG.m)

% Find U (Copied from vqKmeans.m)
dim = size(data, 1);
dataNum = size(data, 2);
centerNum = size(center, 2);
% ====== Compute distance matrix
distMat=pairwiseSqrDistance(center, data);
% ====== Find the U (partition matrix)
[minDist, colIndex] = min(distMat);
U = zeros(size(distMat));
U(colIndex+centerNum*(0:dataNum-1)) = 1;

% Plot data
DS.input=data;
[junk, DS.output]=max(U);
dcprDataPlot(DS);

% Display the centers
clusterNum = size(center,2);
for i=1:clusterNum,
	line(center(1,i), center(2,i), 'color', 'k', 'marker', 'o', 'linestyle', 'none', 'linewidth', 2, 'markerSize', 10);
end