c均值聚类matlab算法和程序

c均值聚类的matlab程序，里面含有源程序，图片，算法介绍，还包括一份报告，详细介绍了算法的过程，希望读者都能看懂。

function[IX2]=fcm（IM）;
%IM是输入的原图象
%IX2是分类结果
[maxXmaxY]=size（IM）;

IMM=cat（3IMIMIM）;
%初始化聚类中心（3类）
cc1=8;
cc2=100;
cc3=200;


ttFcm=0;
while（ttFcm<15）
    ttFcm=ttFcm+1;
    c1=repmat（cc1maxXmaxY）;
    c2=repmat（cc2maxXmaxY）;
    c3=repmat（cc3maxXmaxY）;
    c=cat（3c1c2c3）;
    ree=repmat（0.000001maxXmaxY）;
    ree1=cat（3reereeree）;
    distance=IMM-c;
    distance=distance.*distance+ree1;
    daoShu=1./distance;
    daoShu2=daoShu（;;1）+daoShu（;;2）+daoShu（;;3）;
    
    
    
    
    %计算隶属度u
    distance1=distance（;;1）.*daoShu2;
    u1=1./distance1;
    distance1=distance（;;2）.*daoShu2;
    u2=1./distance2;
    distance1=distance（;;3）.*daoShu2;
    u3=1./distance3;
    
    %计算聚类中心z
    ccc1=sum（sum（u1.*u1.*IM））/sum（sum（u1.*u1））;
    ccc2=sum（sum（u2.*u2.*IM））/sum（sum（u2.*u2））;
    ccc3=sum（sum（u3.*u3.*IM））/sum（sum（u3.*u3））;
    
    tmpMatrix=[abs（cc1-ccc1）/cc1abs（cc2-ccc2）/cc2abs（cc3-ccc3）/cc3];
    pp=cat（3u1u2u3）;
    
    for i=1:maxX
    for j=1:maxY
        if max（pp（ij））==u1（ij）
            IX2（ij）=1;
        elseif max（pp（ij:））==u2（ij）
            IX2（ij）=2;
        else
            IX2（ij）=3;
        end
    end
end
%判断结束条件
   if max（tmpMatrix）<0.0001
       break;
   else
       cc1=ccc1;
       cc2=ccc2;
       cc3=ccc3;
   end
   
   for i=1:maxX
       for j=1:maxY
           if IX2（ij）==3
               IMMM（ij）=240;
           elseif IX2（ij）==2
               IMMM（ij）=130;
           else
               IMMM（ij）=20;
           end        
       end
   end
   %显示每次聚类分割结果
   figure（1）;
   imshow（uint8（IMMM））;
   
   end
   
   for i=1:maxX
       for j=1:maxY
           if IX2（ij）==3
               IMMM（ij）=240;
           elseif IX2（ij）==2
               IMMM（ij）=130;
           else
               IMMM（ij）=20;
           end
       end
   end
   
   %显示最终分类结果
   IMMM=uint8（IMMM）;
   figure（2）;
   imshow（IMMM）;
   
   end
   
   

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

     文件     263222  2007-12-19 09:48  实验1\1998.09.bmp

     文件       2083  2010-03-31 18:17  实验1\fcm.m

     文件       2195  2010-03-31 19:30  实验1\fcm3.m

     文件       2630  2010-03-31 19:46  实验1\fcm4.asv

     文件       2629  2010-04-13 16:52  实验1\fcm4.m

     文件      81765  2010-03-31 20:15  实验1\fig.fig

     文件     419837  2010-03-31 20:15  实验1\fig3.fig

     文件      81784  2010-03-31 20:15  实验1\figure2.fig

     文件        975  2010-03-31 14:59  实验1\shang.m

     文件       1029  2010-03-31 18:18  实验1\shiyanyi2.m

     文件     157025  2009-06-03 08:22  实验1\vvvvv.JPG

     文件     907776  2010-06-15 12:01  实验1\实验1报告.doc

     目录          0  2010-06-15 12:01  实验1

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

              1922950                    13