main.m

ISODATA算法

function main()
%@project:ISODATA算法
%@author:fantasy
%@date:2006.5.4

%@ISODATA算法主函数
%产生测试数据
clear;
Dim=2;
c=4;
Nc=2;
Selta_n=2;
Selta_s=1;
Selta_D=8;
L=1;
I=20;

Patterns=TestData(Dim,c);
Classes=ISODATA(Patterns,c,Nc,Selta_n,Selta_s,Selta_D,L,I);
subplot(1,2,1),plot(Patterns(1,:),Patterns(2,:),'.');

subplot(1,2,2),
hold on;

Color={'r.','b.','g.','y.'};
for iClass=1:length(Classes)
    Patterns=getPatterns(Classes(iClass));
    plot(Patterns(1,:),Patterns(2,:),Color{iClass});
end


%*************************************************************
%产生测试数据
function X=TestData(n,c)
Z=[randn(c,n)*20];%模式中心
PatternNum=10;%每类模式数
X=[];
for(i=1:c)
    w=[];
    for(j=1:PatternNum)
        w=[w;Z(i,:)];
    end
    X=[X;randn(PatternNum,n)+w];
end
X=X';