clear?all%清除workspace中的所有變量
I=imread（‘4.bmp‘）;%讀入圖像文件
figure（1）;???????????
imshow（I）;title（‘原圖像‘）;?%顯示原來的圖像
count=imhist（I）;???%灰度頻率
[mn]=size（I）;%測矩陣的大小
h=count/（m*n）%求灰度級數
V1=zeros（2）;?????%確定對稱矩陣A
V2=zeros（2）;
U1=zeros（2256）;???%設定分區矩陣U
m1=2;??????????????%確定加權指數m
V1（1）=20;??????????%初始化聚類中心V，C=2
V1（2）=250;
l=length（h）
I=double（I）;%轉換變量的類型
diedai=1;
while?（diedai==1）?
????for?i=1:l??????%修改迭代過程中的分區矩陣
?????????if?（i-V1（1））^2>0.01&（i-V1（2））^2>0.01
?????????????tep1=（i-V1（1））^2;
?????????????tep2=（i-V1（2））^2;
?????????????U1（1i）=1.0/（1+（tep1/tep2）^（1.0/（m1-1）））;
?????????????U1（2i）=1.0/（1+（tep2/tep1）^（1.0/（m1-1）））;
?????????else
?????????????if?（i-V1（1））^2<0.01
????????????????U1（1i）=1.0;
?????????????else
????????????????U1（2i）=0.0;
?????????????end
?????????????if?（i-V1（2））^2<0.01
????????????????U1（1i）=0.0;
?????????????else
????????????????U1（