绝对值减数法的模糊聚类分析的matlab程序解释

求模糊相似矩阵

a=zeros(1,N) %创建空矩阵加入excel里的数据

X=[unnamed] %输入样本矩阵

I=ones(31,31); %I是元素全为1矩阵

for i=1:31

for j=1:31

Y(i,j)=(abs(X(i,1)- X(j,1))+abs(X(i,2)- X(j,2))+ abs(X(i,3)- X(j,3))); %绝对值减数标定算法

end

end

R=round((I-0.005*Y)*100)/100 %得到模糊相似矩阵

计算R的传递闭包

B=zeros(31);

flag=0; %设置标志

C=R;

while flag==0

for i=1:31

for j=1:31

for k=1:31

B(i,j)=max(min(C(i,k),C(k,j)),B(i,j)); %计算传递闭包，R与R内积，先取小再取大

end

end

end

if B==C

flag=1;

else

C=B; %循环计算R传递闭包

end

end

B

t = B %求出传递闭包

计算K截矩阵

L=unique(B)';

D=zeros(31);

for m=length(L):-1:1

k=L(m);

for i=1:31

for j=1:31

if B(i,j)>=k

D(i,j)=1;

else D(i,j)=0; %求截距阵，当bij≥λ时，bij(λ) =1；当bij＜λ时，bij(λ)=0 end

end

end

fprintf('当分类系数k=：\n');

disp(L(m));

fprintf('所得截距阵为：\n');

disp(D);

end