动态聚类MATLAB 代码

根据矩阵计算聚类中心 并分类，再计算聚类中心，迭代出最终结果

clc
X=[80 10 6 2;50 1 6 4;90 6 4 6;40 5 7 3;10 1 2 4];
for j=1:4
    zuida=max（X（:j））;
    zuixiao=min（X（:j））;
    for i=1:5
       
        X（ij）=X（ij）/zuida;
    end
end

%求模糊相似矩阵 最大最小法 

R=zeros（55）;
for i=1:5
    for j=1:5
        quxiao=0;
        quda=0;
        for k=1:4
            quxiao=quxiao+min（X（ik）X（jk））;
            quda=quda+max（X（ik）X（jk））;
            R（ij）=quxiao/quda;
            
                
        end
    end
end


%求传递闭包
 a=size（R）;
B=zeros（a）;
flag=0;
while flag==0
for i= 1: a
 for j= 1: a
 for k=1:a
 B（ i  j ） = max（min（ R（ i  k）  R（ k j） ）  B（ i  j ） ） ;%R与R内积，先取小再取大
 end
 end
end
if B==R
 flag=1;
else
 R=B;%循环计算R传递闭包
end
end
R
B=reshape（R125）;
B=unique（B）;
B=sort（B）
C=size（B）;
E=C（2）;
D=zeros（55E）;
for k=1:E
    for i=1:5
        for j=1:5
            if R（ij）>=B（k）
                D（ijk）=1;
            else
                D（ijk）=0;
            end
        end
    end
end
D
    

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

     文件       1054  2015-12-23 19:28  动态聚类\dongtaijulei.m

     目录          0  2015-12-23 21:00  动态聚类

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

                 1054                    2