som_demo1.m

patten regnization source从1-14章能运行

clear
clf,figure(1)
N=100; % 阿 努矾胶磐狼 钎夯狼 肮荐
N2=N+N; B1=ceil(N/2); B2=N+B1;
eta=0.2;
means=[0.7  -0.8      
       0.7  -0.8];
var= [0.2  0.2];

x=datagen([N N],[means;var]); % x : 2N by 2
x=randomize(x);% 青狼 蔼阑 罚待窍霸 集扁
disp('11俺狼 春快繁阑 荤侩窍咯 急屈利栏肺 硅凯阑 父甸绢 1俊辑 11鳖瘤 扼骇傅.');
ncenter=11; % 荤侩且 努矾胶磐 春繁狼 肮荐
w=rand(ncenter,2)-0.5*ones(ncenter,2);  % 檬扁 春快繁篮 漂隆 傍埃惑俊 烙狼肺 困摹矫糯.

subplot(121),plot(x(:,1),x(:,2),'r.',w(:,1),w(:,2),'*-')
axis([-2 2 -2 2])
title('檬扁拳')
i=1; iter=1; converge=0;
while converge==0,
   dn=ones(ncenter,1)*x(i,:)-w;
   ddn=sum((dn.*dn)')'; % ddn: ncenter by 1
   [tmp,istar]=min(ddn);
   if istar==1,
      w([istar:istar+1],:)=w([istar:istar+1],:)+eta*(ones(2,1)*x(i,:)-w([istar:istar+1],:));
   elseif istar==ncenter,
      w([istar-1:istar],:)=w([istar-1:istar],:)+eta*(ones(2,1)*x(i,:)-w([istar-1:istar],:));
   else
      w([istar-1:istar+1],:)=w([istar-1:istar+1],:)+eta*(ones(3,1)*x(i,:)-w([istar-1:istar+1],:));
   end
   subplot(122),
   plot(x(:,1),x(:,2),'r.',x(i,1),x(i,2),'o',w(:,1),w(:,2),'*-')
   title(['Iteration = ' num2str(iter)])
   drawnow
   i=rem(i+1,N2); iter=iter+1;
   if i==0,
      x=randomize(x);  % 努矾胶磐狼 鉴辑甫 罚待窍霸 集绰促 
      i=1;             % 犁硅摹等 涝仿栏肺 矫累 
   end
   if rem(iter,50)==0,
      eta=eta*0.9;
      %if isempty(converge), converge=0; end
      if iter >=200, converge=1; end
   end
end
subplot(122),plot(x(:,1),x(:,2),'r.',w(:,1),w(:,2),'*-')
text(w(1,1),w(1,2)+0.2,'1'), text(w(ncenter,1),w(ncenter,2)+0.2,int2str(ncenter))
title(['Iteration = ' num2str(iter)])
axis([-2 2 -2 2])