matlab 三维dbscan代码

以下是MATLAB中实现三维DBSCAN算法的示例代码：

function [IDX, isnoise]=DBSCAN(X,epsilon,MinPts)
    C=0;
    n=size(X,1);    
    D=pdist2(X,X);    
    visited=false(n,1);    
    IDX=zeros(n,1);
    for i=1:n
        if ~visited(i)            
            visited(i)=true;
            Neighbors=RegionQuery(i);
            if numel(Neighbors)<MinPts
                IDX(i)=-1;             
            else
                C=C+1;               
                ExpandCluster(i,Neighbors,C);
            end
        end
    end
    
    isnoise=IDX==-1;
    IDX(isnoise)=0;
    
    function ExpandCluster(i,Neighbors,C)
        IDX(i)=C;
        k = 1;
        while true
            j = Neighbors(k);
            if ~visited(j)
                visited(j)=true;
                Neighbors2=RegionQuery(j);
                if numel(Neighbors2)>=MinPts
                    Neighbors=[Neighbors Neighbors2]; %#ok
                end
            end
            if IDX(j)==0
                IDX(j)=C;
            end
            
            k = k + 1;
            if k > numel(Neighbors)
                break;
            end
        end
    end
 
    function Neighbors=RegionQuery(i)
        Neighbors=find(D(i,:)<=epsilon);
    end
 
end

该函数接受三个输入参数：X是一个n行3列的矩阵，表示三维空间中的点集；epsilon是邻域半径；MinPts是最小点数阈值。输出参数包括数据点的簇分配结果和噪声点的布尔向量。在函数内部，使用RegionQuery函数来获取每个点的邻域，使用ExpandCluster函数扩展簇，最终生成一个簇分配结果。