evolution function.m

协同神经网络的程序

%序参量演化方程


%the number of order parameter
NumberOfOP = NumOfImage;
%the number of step
NumberOfStep = 400;

%prepare attention parameter
for m=1:NumberOfOP
%    lmd(m)=1;
     lmd(m)=1-abs(vmean(m)-Qmean)/max(Qmean,255-Qmean);
end
%lmd(1)=1.2;

maxL = -1000;
minL = 1000;
for k=1:NumberOfOP
   if(lmd(k)>maxL) maxL = lmd(k);end
   if(lmd(k)<minL) minL = lmd(k);end
end
   
%parameter
B=1;C=1;

for n=1:NumberOfStep
%    n
    %global field 
    D(n)=0;
    for m=1:NumberOfOP
      D(n)=D(n)+(B+C)*A(n,m)*A(n,m);
    end
    
    %evloution function
    speed = 0.1/D(n);
    %speed = 0.3;
    
    for k=1:NumberOfOP
      A(n+1,k)=speed * A(n,k)*(lmd(k)-D(n)+B*A(n,k)*A(n,k))+A(n,k);
    end
    
    %?draw the evolution progress of input vector
    mm = fix(NumberOfStep/10);
    if (0)
	 %if ((n - fix(n/mm)*mm)==0)
	   v3= zeros(1,width*hight);
	   for k=1:NumberOfOP
		   v3 = v3 + v(k,:)*A(n,k) ;
	   end
	   v3 = reshape(v3,hight,width);
	   v3 = mat2gray(v3);
   	imshow (v3);
      xlabel(n);
      pause
   end    
 end
 
 
 
%A is the order parameter,draw it;
figure
plot(A)
axis([1,NumberOfStep,-0.5,max(sqrt(minL),minA)]+0.3)
%plot(D)

figure
v3= zeros(1,width*hight);
for k=1:NumberOfOP
	v3 = v3 + v(k,:)*A(NumberOfStep,k) ;
end
v3 = reshape(v3,hight,width);
v3 = mat2gray(v3);
imshow (v3);