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