求DBSCAN算法的matlab代码！！

旖罗

各位谁有DB SCAN算法的matlab代码啊？急用啊。:@Q

风花雪月

数据挖掘中聚类算法之DBSCAN算法的matlab实现

1. function [class]=dbscan(X,k,eps)
   
2. [m,n]=size(X);
   
3. class=zeros(1,m);
   
4. clusterId = 0;
   
5. X=[class',X];
   
6. 
   
7. for i = 1:m
   
8.   clusterId = clusterId + 1;
   
9.   if X(i,1) == 0
   
10.      [X]=expandcluster(i,X,clusterId,eps,k,n);
    
11.   end  
    
12. end
    
13. 
    
14. class = X(:,1)';
    
15. 
    
16. function [X]=expandcluster(index,X,clusterId,eps,k,n)
    
17.   Di=dist(X(index,(2:n)),X(:,(2:n)));
    
18.   D = find(Di<eps);
    
19.   if length(D) < k
    
20.     X(index,1) = -1; %noise
    
21.     return;
    
22.   else
    
23.     X(index,1) = clusterId;
    
24.     D(1)=[];  
    
25.     for i=1:length(D)
    
26.       Dmatrix(i,:)=X(D(i),:);
    
27.     end
    
28.     while length(D)~=0
    
29.       Di=dist(X(D(1),(2:n)),X(:,(2:n)));
    
30.       Result = find(Di<eps);
    
31.       if length(Result) > k
    
32.         for i = 1: length(Result)   
    
33.             if X(Result(i),1) == 0|-1
    
34.                if X(Result(i),1)==0
    
35.                   D(length(D)+1)=Result(i);
    
36.                end
    
37.                X(Result(i),1)=clusterId;
    
38.             end;
    
39.          end;
    
40.        end;
    
41.         D(1)=[];
    
42.      end;
    
43.    end;
    
44. 
    
45.   
    
46. %function distan = dist(x,y)
    
47. % sum = 0;
    
48.   %nx = size(x);
    
49.   %for i = 1:nx
    
50.    %   sum = (x(i)-y(i))^2+sum;
    
51.    %end
    
52. %distan = sqrt(sum)
    
53. 
    
54. %............................................
    
55. function [Di]=dist(i,x)
    
56. 
    
57. % function: [D]=dist(i,x)
    
58. %
    
59. % Aim:
    
60. % Calculates the Euclidean distances between i and all objects in x     
    
61. %                                    
    
62. % Input:
    
63. % i - an object (1,n)
    
64. % x - data matrix (m,n); m-objects, n-variables        
    
65. %                                                                 
    
66. % Output:
    
67. % D - Euclidean distance (m,1)
    
68. 
    
69. 
    
70. [m,n]=size(x);
    
71. Di=sqrt(sum((((ones(m,1)*i)-x).^2)'));
    
72. 
    
73. if n==1
    
74.    Di=abs((ones(m,1)*i-x))';
    
75. end

复制代码

zcw20004

这个好像是算一维的吧:@L

zcw20004

再试了一次．原来多了一行注释
Error in ==> dbscan_2>dist at 61
objects in x

zcw20004

请问2楼噪声点是如何处理的:@(