DBSCAN 聚类 自手动调节查看 matlab

clc 
clear all
n = 200;
a = linspace(0,8*pi,n/2);
u = [5*cos(a)+5 10*cos(a)+5]'+1*rand(n,1);
v = [5*sin(a)+5 10*sin(a)+5]'+1*rand(n,1);
Data=[u v];
CoreIndex=[];
radius=2;
NUMS=3;
for i=1:size(Data,1)
    Count=0;
    for j=1:size(Data,1)
        if sqrt((Data(i,1)-Data(j,1)).^2+(Data(i,2)-Data(j,2)).^2)<radius
            Count=Count+1;
        end
    end
    if Count>=NUMS
        CoreIndex=[CoreIndex i];
    end
end
%% 改编i的值可以得到不用的聚类，每次只能聚集一种类别，自己可以添加代码完成
i=150
    Set=[];
    for j=1:size(Data,1) 
    if sqrt((Data(CoreIndex(i),1)-Data(j,1))^2+(Data(CoreIndex(i),2)-Data(j,2))^2)<radius
    Set=[Set j];
    end
    end
    while 1
        TempIn=size(Set,2);
    for k=1:TempIn
        for n=1:size(Data,1)
        if sqrt((Data(Set(k),1)-Data(n,1))^2+(Data(Set(k),2)-Data(n,2))^2)<radius
            flag=0;
            for ii=1:size(Set,2)
            if Set(ii)==n
                flag=1;
            end
            end
            if flag ==0
            Set=[Set n];
            end  
        end
        end
    end
    if TempIn==size(Set,2)
            break;
    end
    end
plot(Data(:,1),Data(:,2),'b*');
for i=1:size(Set,2)
   hold on
    plot(Data(Set(i),1),Data(Set(i),2),'ro')
end