test.m

来自「使用分解聚类算法在IRIS数据上进行聚类分析」· M 代码 · 共 37 行

[xx1,xc] = loadtest();
%xx1=[0 0 2 2 4 4 5 6 6 7 -4 -2 -3 -3 -5 1 0 0 -1 -1 -3;6 5 5 3 4 3 1 2 1 0 3 2 2 0 2 1 -1 -2 -1 -3 -5];
x1=xx1';
x2=[0 0];
n=21;
e=zeros(0);
maxe=zeros(0);
[n1,nc1]=size(x1);
n2=n-n1;

for k=1:n-1;
e=zeros(1,n1);
for i=1:n1;
    e(1,i)=counte(i,x1,x2,n);
end
maxe(1,k)=max(max(e));         %求e值矩阵的最大值
[r,c]=find(e==maxe(1,k));      %r为最大值所在的行数，c为列数
temp=x1(c,:);
x2(k,:)=temp(:,:);
x1(c,:)=[];
n1=n1-1;
n2=n2+1;
if k>1
    if maxe(1,k)<maxe(1,k-1)
    break;
    end
end
end

disp(e);
disp(x1);
%%%disp(maxe);


        %%temp=a1(c,:);
        %%a2(k,:)=temp(:,:);