matlab 密度聚类中心点,密度峰值聚类算法MATLAB程序

密度峰值聚类算法MATLAB程序

数据见：MATLAB中“fitgmdist”的用法及其GMM聚类算法，保存为gauss_data.txt文件，数据最后一列是类标签。

1. MATLAB程序

clear all

close all

data_load=dlmread('gauss_data.txt');

[num,dim]=size(data_load); %数据最后一列是类标签

data=data_load(:,1:dim-1); %去掉标签的数据

mdist=pdist(data); %两两行之间距离

A=tril(ones(num))-eye(num);

[x,y]=find(A~=0);

% Column 1: id of element i, Column 2: id of element j', Column 3: dist(i,j)'

xx=[x y mdist'];

ND=max(xx(:,2));

NL=max(xx(:,1));

if (NL>ND)

ND=NL;

end

N=size(xx,1);

for i=1:ND

for j=1:ND

dist(i,j)=0;

end

end

for i=1:N

ii=xx(i,1);

jj=xx(i,2);

dist(ii,jj)=xx(i,3);

dist(jj,ii)=xx(i,3);

end

percent=2.0;

fprintf('average percentage of neighbours (hard coded): %5.6f\n', percent);

position=round(N*percent/100);

sda=sort(xx(:,3));

dc=sda(position);

fprintf('Computing Rho with gaussian kernel of radius: %12.6f\n', dc);

for i=1:ND

rho(i)=0.;

end

%

% Gaussian kernel

%

fo