一个聚类算的实.docVIP

下午我在实验室挥汗如雨的工作，某君拿着一份东南大学的学报找到我，问我能不能在Matlab里实现资料中的那个聚类算法，我看了一下，觉得这个算法不适合使用Matlab实现，或者说我不太会在Matlab中实现这样的算法，但是这段时间一直在使用Matlab，其它语言相对生疏了，且Matlab这种解释执行的环境比较适合调试，开发速度肯定比C快多了，所以还是试着做了一下。 ????? 简略地看了一下学报中的这篇文章，算法叫做广度优先邻居搜索的聚类算法。里面有几个定义： ????? 一、数据矩阵。这是一个m行2列的矩阵，行数m=待分类的元素个数，第一列代表元素横坐标，第二列代表元素纵坐标； ????? 二、相似度矩阵simiR。这是一个归一化了的矩阵，simiR(a,b)代表元素a和b的相似度，相似度越大越相似。该矩阵是一个对称，且对角线元素为1的矩阵； ????? 三、邻居门限r。如果有两个元素a,b，且simiR(a,b)r，则这两个元素是邻居； ????? 四、类门限lambda。类A中有n个元素，有一个待分类的元素k，若类A中有超过(lambda*100)%个元素是k的邻居，则k属于类A； ????? 五、广度优先邻居搜索。即搜索某个元素的所有邻居，然后分别搜索每个邻居的所有邻居，依此； ????? 使用广度优先，正好我前几天看了一个广度优先抓取网页的Python程序，马上就能想起（参见 HYPERLINK /2007-11/119552581137716.html \t _blank 《用Python写一个小小的爬虫程序》）。 ????? 广度优先的典型算法就是使用队列：找出队首的元素的所有邻居，把它们从队尾压入，然后将队首弹出，非常方便。 ????? 如果把互为邻居的元素用无向线段相连，则所有的元素构成一个非连通图，需要多次调用广度优先搜索才能完成所有元素的分类。 ????? 下面给出代码（Matlab7.0.1）： createSimiR.m ????? 这个函数不是我写的，是某君给我的生成测试数据的代码中的函数，如果哪位看客知道这段代码的出处，请及时通知我，我会在文章中注明出处。 ????? 生成归一化相似度矩阵的函数，参数为待分类数据，返回值为归一化相似度矩阵。 function simiR=createSimiR(INPUTDATA)Xk=mean(INPUTDATA); %列均值Sk=var(INPUTDATA); %列方差numrows=size(INPUTDATA,1); %INPUTDATA行数for i=1:numrows??? INPUTDATA(i,:)=(INPUTDATA(i,:)-Xk)./Sk; %标准差标准化endfor i=1:numrows??? for j=1:2??? DATAFORCLUS(i,j)=(INPUTDATA(i,j)-min(INPUTDATA(:,j)))...??????? /(max(INPUTDATA(:,j))-min(INPUTDATA(:,j))); %极值标准化??? endend%--------相似矩阵----------%for i=1:numrows??? for j=1:numrows??????? simiR(i,j)=1-0.5*(abs(DATAFORCLUS(i,1)-DATAFORCLUS(j,1))...??????????? +abs(DATAFORCLUS(i,2)-DATAFORCLUS(j,2)));??? endend ????? 下面的代码都是我写的，转载注明出处，不为别的，如果哪位看客发现了问题或者有好的建议，能够方便地通知我。 BFSN_Algorithm.m ????? 分类算法，参数分别为G（相似度矩阵），r（邻居门限），lambda（类门限），返回值A，一个cell数组，每个元素是一个向量，包含了一个类的所有元素。 function A=BFSN_Algorithm(G,r,lambda)%广度优先搜索邻居的聚类算法实现%G为相似度矩阵%r和lambda为参数%r为邻居门限，相似度大于r即为邻居%lambda门限%未分类元素对于某类的所有元素，如果是邻居则令X(i)=1，否则为0.i为类元素的下标%对X求和并除以类元素个数，若此值大于lambda门限，则该未分类元素属于这类 %打开计时器tic%A为聚类结果A={};%k为分类计数k=1;%n为待分类元素个数n=size(G,1);%构建元素向量member=1:n;%只要member中还有未分类的元素就继续循环while numel(member)~=0??? %从member中取出一个