Gap Statistic算法 Matlab代码

总共三个文件

gap_stat.m

clear all;close all;
SampleNum=30;                                %样本观测点的数量
B=1000;                                      %参考数据集的数量
MaxK=10;                                     %最大的聚类数
u=[0 0;4 4;4 -4];                            %各类样本的均值向量
sigma=[1 0;0 1];                             %协方差矩阵


SampleSet=[];
[uM,~]=size(u);
for i=1:uM                                   %生成样本数据
    SampleSet=[SampleSet;mvnrnd(u(i,:),sigma,fix(SampleNum/uM))];
end
figure(1);
plot(SampleSet(:,1),SampleSet(:,2),'.');     %画二维样本分布图
title('样本数据');


Wk=log(CompuWk(SampleSet,MaxK));             %计算样本的Wk


interval=minmax(SampleSet');                 %样本各维度的取值区间
[SampleM,SampleN]=size(SampleSet);
RefSet=zeros(SampleM,SampleN,B);
for b=1:B