function?[IX2]=fcm（IM）;
%IM是輸入的源圖像
%IX2是分類結(jié)果
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