clc
X=[80?10?6?2;50?1?6?4;90?6?4?6;40?5?7?3;10?1?2?4];
for?j=1:4
????zuida=max（X（:j））;
????zuixiao=min（X（:j））;
????for?i=1:5
???????
????????X（ij）=X（ij）/zuida;
????end
end

%求模糊相似矩陣?最大最小法?

R=zeros（55）;
for?i=1:5
????for?j=1:5
????????quxiao=0;
????????quda=0;
????????for?k=1:4
????????????quxiao=quxiao+min（X（ik）X（jk））;
????????????quda=quda+max（X（ik）X（jk））;
????????????R（ij）=quxiao/quda;
????????????
????????????????
????????end
????end
end


%求傳遞閉包
?a=size（R）;
B=zeros（a）;
flag=0;
while?flag==0
for?i=?1:?a
?for?j=?1:?a
?for?k=1:a
?B（?i??j?）?=?max（min（?R（?i??k）??R（?k?j）?）??B（?i??j?）?）?;%R與R內積，先取小再取大
?end
?end
end
if?B==R
?flag=1;
else
?R=B;%循環計算R傳遞閉包
end
end
R
B=reshape（R125）;
B=unique（B）;
B=sort（B）
C=size（B）;
E=C（2）;
D=zeros（55E）;
for?k=1:E
????for?i=1:5
????????for?j=1:5
????????????if?R（ij）>=B（k）
????????????????D（ijk）=1;
????????????else
????????????????D（ijk）=0;
????????????end
????????end
????end
end
D
????

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

?????文件???????1054??2015-12-23?19:28??動態聚類\dongtaijulei.m

?????目錄??????????0??2015-12-23?21:00??動態聚類

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

?????????????????1054????????????????????2