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