Chameleon变色龙层次聚类算法实现

Chameleon变色龙层次聚类算法实现

clear; 
clc; 
 
%产生随机点 
img=zeros（100200）; 
data=rand（1100）; 
x=zeros（1100）; 
y=zeros（1100）; 
for i=1:100 
    loc=ceil（data（1i）*20000）; %x为一个包含小数点的数组ceil（x）就是把x中所有小数点后面的去掉 整数位上加1
    x（1i）=ceil（loc/200）; 
    y（1i）=mod（loc200）; 
    if （y（1i）==0） 
        y（1i）=200; 
    end 
    img（x（1i）y（1i））=255; 
end 
%image（img）; 
 
%构造距离矩阵 
processed=1000000; 
for i=1:100 
    for j=i+1:100 
       dist（ij）=（x（1i）-x（1j））*（x（1i）-x（1j））+（y（1i）-y（1j））*（y（1i）-y（1j））; 
       dist（ji）=processed; 
   end 
   dist（ii）=processed; 
end 
 
for i=1:100 
    category（1i）=i; 
end 
num_category=100; 
%反复挑选距离最近的两类合并 
while （num_category>7） 
    mindist=1000000; 
 
    for i=1:100 
        for j=i+1:100 
            if （dist（ij）                mindist=dist（ij）; 
                minp1=i; 
                minp2=j; 
            end 
        end 
    end 
    dist（minp1minp2）=processed; 
    dist（minp2minp1）=processed; 
    for i=1:100 
        if （category（1i） == minp2） 
            category（1i）=minp1; 
        end 
    end 
    for i=1:100 
        if （dist（minp2i）            dist（minp1i）=dist（minp2i）; 
        end 
        dist（minp2i）=processed; 
        dist（iminp2）=processed; 
    end 
    num_category=num_category-1; 
end 
 
%整理类别标记 
for i=1:100; 
    labeled（1i）=0; 
end 
categories=0; 
for i=1:100 
    kind=category（1i）; 
    if （labeled（1kind）==0） 
        categories=categories+1; 
        for j=i:100 
            if（category（1j）==kind） 
                category（1j）=categories; 
            end 
        end 
        labeled（1categories）=1; 
    end 
end 
 
%显示 
x1=zeros（1100）; 
y1=zeros（1100）; 
x2=zeros（1100）; 
y2=zeros（1100）; 
x3=zeros（1100）; 
y3=zeros（1100）; 
x4=zeros（1100）; 
y4=zeros（1100）; 
x5=zeros（1100）; 
y5=zeros（1100）; 
x6=zeros（1100）; 
y6=zeros（1100）; 
x7=zeros（1100）; 
y7=zeros（1100）; 
 
for i=1:100 
    if（category（1i）==1） 
        x1（1i）=x（1i）; 
        y1（1i）=y（1i）; 
    elseif （category（1i）==2） 
        x2（1i）=x（1i）; 
        y2（1i）=y（1i）; 
    elseif （category（1i）==3） 
        x3（1i）=x（1i）; 
        y3（1i）=y（1i）; 
    elseif （category（1i）==4） 
        x4（1i）=x（1i）; 
        y4（1i）=y（1i）; 
    elseif （category（1i）==5） 
        x5（1i）=x（1i）; 
        y5（1i）=y（1i）; 
    elseif （category（1i）==6） 
        x6（1i）=x（1i）; 
        y6（1i）=y（1i）; 
    elseif （category（1i）==7） 
        x7（1i）=x（1i）; 
        y7（1i）=y（1i）; 
    end 
end 
plot（x1y1‘b.‘x2y2‘g.‘x3y3‘r.‘x4y4‘c.‘x5y5‘m.‘x6y6‘y.‘x7y7‘k.‘）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    1838080  2014-04-28 21:02  Chameleon算法.ppt

     文件       2750  2014-04-29 22:03  chameleon.m

----------- ---------  ---------- -----  ----

              1840830                    2