Clustering by fast search and find of density(2)实现

 

matlab实现
 clear;clc;

cluster = 3;

center = 50*randn(3,2);

dcc = 10;

numOfpoint = 300;

for i = 1:numOfpoint

centerindex = randperm(3,1);

centerx = center(centerindex,1);

centery = center(centerindex,2);

theta = rand(1)*2*pi;

radius = dcc*rand(1);

point(i,1) = centerx + radius*cos(theta);

point(i,2) = centery + radius*sin(theta);

end

plot(center(:,1),center(:,2),'r*')

hold on

plot(point(:,1),point(:,2),'bo')

P = [center; point];

distance = zeros(size(P,1));

% first compute distance of each pair

for i = 1:size(P,1)

for j = i+1:size(P,1)

distance(i,j) = sqrt(sum((P(i,:)-P(j,:)).^2));

distance(j,i) = distance(i,j);

end

end

% compute local density of every point

for i = 1:size(P,1)

density(i,1) = length(find(distance(i,:)<dcc))-1;

end

%compute delta

maxdensity = find(density == max(density));

for i = 1:size(P,1)

tempindex = find(density>density(i));

if ~isempty(tempindex)

delta(i,1) = min(distance(i,tempindex));

else

delta(i,1) = max(distance(i,:));

end

end

figure

plot(density,delta,'ro')
--------------------- 
作者:cuifan0814 
来源:CSDN 
原文:https://blog.csdn.net/cuifan0814/article/details/49787691 
版权声明:本文为博主原创文章,转载请附上博文链接!

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值