chuandibibao.m

来自「模式识别-模糊聚类算法--传递闭包法」· M 代码 · 共 53 行

tic
close all;
clear all;
clc

X=[1,2,2,3,3,5,6,6,5,4;
    3.5,2.5,3.5,3,3,0.5,1.5,1.5,3,3;
    1,2,1,3,1,5,4,4,2,1;
    0,2,1,1,1,2,0,1,2,2];
[m,n]=size(X);
R=zeros(n,n);
R=X'*X;
R=R/(max(max(R)));
for I=1:n
    for J=1:n
        if I==J
            R(I,J)=1;
        end
    end
end
Q=R;
C=pingfang(R);
while R~=C
    R=C;
    C=pingfang(C);
end
R=C;

S=R;
[m,n]=size(R);
lamda=0.6;
for I=1:m
    for J=1:n
        if R(I,J)>=lamda
            R(I,J)=1;
        else
            R(I,J)=0;
        end
    end
end
R=rref(R);
sprintf('经模糊处理，以上数据可以分为%d类',rank(R))
for I=1:rank(R)
    s=[];
    for J=1:n
        if R(I,J)~=0
           s=[s,J];
        end
    end
     sprintf('第%d类为',I)
     sprintf('X(%d)',s)
end