SOM Neural Nework

clear all;
clc;

X=[1 2;
   3 2;
   2 3;
   4 1;
   9 9;
   9 8;
   10 9;
   9 10;
   1 8;
   2 9;
   3 8;
   4 10]

W=[1 2;
   10 9;
   1 8]

[m, n]=size(X);
[m1, n1]=size(W);

alfa=0.1;
decalfa=0.1;
maksiterasi=130;
iterasi=0;

while (iterasi <= maksiterasi)
    iterasi=iterasi+1;
    for i=1:m,
        for j=1:m1,
            d(j)=0;
            for k=1:n,
                d(j)=d(j)+(X(i,k)-W(j,k))^2;
            end,
            d(j)=sqrt(d(j));
            [jmin idx]=min(d);
           
            W(idx,:)=W(idx,:)+alfa*(X(i,:)-W(idx,:));
        end,
    end,
    alfa=alfa-alfa*decalfa;
    disp(strcat('Iterasi = ',int2str(iterasi),' Alfa = ',num2str(alfa)));
   
end,
W
H=[];
for i=1:m,
    for j=1:m1,
        d(j)=0;
        for k=1:n,
            d(j)=d(j)+alfa*(X(i,k)-W(j,k))^2;
        end,
        d(j)=sqrt(d(j));
        [jmin idx]=min(d);
    end,
    H=[H idx];
end,
H