%BY?YY??模糊聚類算法?2008.3.4
%采用數(shù)量積求相似矩陣?
%輸入樣本矩陣
X=input（‘please?input?X:‘）;
[NM]=size（X）;
%N是樣本數(shù)?M是特征數(shù)
%N個(gè)樣本M個(gè)指標(biāo)的平均值矩陣和標(biāo)準(zhǔn)差矩陣
P=zeros（M1）;
for?j=1:M
????for?i=1:N
????????P（j）=P（j）+X（ij）;
????end
????P（j）=P（j）/N;
end
P1=zeros（M1）;
for?j=1:M
????for?i=1:N
????????P1（j）=P1（j）+（X（ij）-P（j））*（X（ij）-P（j））;
????end
????P1（j）=P1（j）/N;
????P1（j）=sqrt（P1（j））;
end
%標(biāo)準(zhǔn)化將數(shù)據(jù)歸一化
for?i=1:N
????for?j=1:M
????????X1（ij）=（X（ij）-P（j））/P1（j）;
????end
end
for?j=1:M
???A（j）=min（X1（:j））;%最小
???B（j）=max（X1（:j））;%最大
end?
for?i=1:N
????for?j=1:M
????????X2（ij）=（X1（ij）-A（j））/（B（j）-A（j））;
????end?
end?
X2
%以上都是數(shù)據(jù)標(biāo)準(zhǔn)化
%下面建立模糊相似矩陣
R1=zeros（NN）;
for?i=1:N
????for?k=1:N
????????if?i==k?;
????????else?
????????????
???????????for?j=1:M
????????????R1（ik）=?R1（ik）+X2（ij）*X2（kj）;
???????????end
????????end
????end
end
m=max（R1）;
m=max（m）;
R=zeros（NN）;
for?i=1:N
????for?j=1:N
????????if?i==j
????????????R（ij）=1;
????????else
????????????R（ij）=R1（ij）/m;
????????end
????end
end
R
%建立模糊等價(jià)關(guān)系矩陣
flag=0;
C=R;
T=zeros（NN）;
while?flag==0
????for?i=1:N
????????for?j=1:N
????????????for?k=1:N
????????????T（ij）=max（min（C（ik）C（kj））T（ij））;%計(jì)算傳遞閉包?不能少T（ij）這項(xiàng)?否則T（i，i）將不再為1?取大取小合運(yùn)算
????????????end
????????end
????end
????if?T==C
????????flag=1;
????else
????????C=T;
????end
end
T?%輸出傳遞閉包
%取不同的置信水平q，q需要設(shè)定
q=0.25;
for?i=1:N
????for?j=1:N
????????if?T（ij）>=q
????????????T1（ij）=1;
????????else
????????????T1（ij）=0;
????????end
????end
end
T1?%T1為一個(gè)P截矩陣?如有對(duì)應(yīng)的列向量相等，則是同一類
%輸出結(jié)果
F=zeros（NN）;
k=1;
for?i=1:N
????e=1;
????F（ke）=i;
????for?j=1:N
????????if?i==j;
????????else
????????????if?isequal（T1（i:）T1（j:））
???????????????e=e+1;
???????????????F（ke）=j;
????????????end
????????end
????end
????k=k+1;
end
F
????????????



????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????2013??2008-03-11?15:40??FCM.m

?????文件???????3703??2008-03-06?01:41??fcm.C

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

?????????????????5934????????????????????3