模糊c均值算法

close all;
clc;
tic

%IM是输入的原图像
%IX2是分类结果
tmp=imread（‘Lenna.tif‘）;
figure（1）
imshow（tmp）;

IM=double（tmp）;

[maxXmaxY]=size（IM）;

IMM=cat（3IMIMIM）;
%初始化聚类中心（4类）
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;
    distance2=distance（::2）.*daoShu2;
    u2=1./distance2;
    distance3=distance（::3）.*daoShu2