matlab silhouette函数,相当于Matlab的聚类质量函数？

我在下面展示了一个在MATLAB和Python/Numpy中的示例silhouette实现(请记住，我对MATLAB比较流利)：

1) MATLAB语言function s = mySilhouette(X, IDX)

%# X : matrix of size N-by-p, data where rows are instances

%# IDX: vector of size N, cluster index of each instance (starting from 1)

%# s : vector of size N, silhouette score value of each instance

N = size(X,1); %# number of instances

K = numel(unique(IDX)); %# number of clusters

%# compute pairwise distance matrix

D = squareform( pdist(X,'euclidean').^2 );

%# indices belonging to each cluster

kIndices = accumarray(IDX, 1:N, [K 1], @(x){sort(x)});

%# compute a,b,s for each instance

%# a(i): average distance from i to all other data within the same cluster.

%# b(i): lowest average dist from i to the data of another single cluster

a = zeros(N,1);

b = zeros(