FCM图像分割算法MATTLABA实现

function [IX2]=fcm（IM）;
%IM是输入的源图像
%IX2是分类结果
IM=imread（‘trees.tif‘）;
%只取RGB的第一个值
IM=IM（::1）;
figure（1）
imshow（uint8（IM））
[maxXmaxY]=size（IM）;
IM=double（IM）;
%IMM是4维矩阵，由IM连接而成
IMM=cat（4IMIMIMIM）;
%初始化聚类中心（4类）
cc1=8;
cc2=80;
cc3=160;
cc4=230;
ttFcm=0;
while（ttFcm<15）
    ttFcm=ttFcm+1;
    c1=repmat（cc1maxXmaxY）;
    c2=repmat（cc2maxXmaxY）;
    c3=repmat（cc3maxXmaxY）;
    c4=repmat（cc4maxXmaxY）;
    c=cat（4c1c2c3c4）;
    ree=repmat（0.000001maxXmaxY）;
    ree1=cat（4reereereeree）;
    
    distance=IMM-c;
    distance=distance.*distance+ree1;
    
    daoshu=1./distance;
    
    daoshu2=daoshu（::1）+daoshu（::2）+daoshu（::3）+daoshu（::4）;
    %计算隶属度u
        distance1=distance（::1）.*daoshu2;
        u1=1./distance1;
        d