clear;
close?all;

aaa=imread（‘chepai.jpg‘）;
date=rgb2gray（aaa）;
[mn]=size（date）;
subplot（221）imshow（date）;title（‘?圖一?原圖像‘）
subplot（222）imhist（date）;title（‘?圖二?原圖像的灰度直方圖‘）
[countx]=imhist（date）;
count=count/sum（count）;
date=double（date）;
v=zeros（2）;
u=zeros（2562）;
v（1）=fix（rand（1）*256）;
v（2）=fix（rand（1）*256）;????????????????????%兩個隨機聚類中心
v1=zeros（2）;
while?（（v1（1）-v（1））^2+（v1（2）-v（2））^2）>0.001
????v1（1）=v（1）;
????v1（2）=v（2）;
????for?i=1:256
???????j=i-1;
?????if?（j-v（1））^2<0.001
????????u（i1）=1.0;
????????u（i2）=0.0;?????????????%屬于接近Ｖ１的歸為一類。
?????else?if（j-v（2））^2<0.001
??????????u（i1）=0.0;
??????????u（i2）=1.0;%屬于接近Ｖ２的歸為一類。　
?????????else%否則，更新隸屬度。　
????????????????tep=0.0;
????????????????for?j1=1:2?
????????????????????tep=tep+1/（k（jj）-2*k（jv（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1719??2010-06-07?10:51??fcm.m

?????文件?????????41??2010-06-07?10:53??k.m

-----------?---------??----------?-----??----

?????????????????1760????????????????????2