order parameter evolution equation.m

协同神经网络的程序

%序参量演化方程


%the number of order parameter
NumberOfOP = 3;
%the number of step
NumberOfStep = 15;

%prepare order parameter
% this is a group of 'good' sample
A(1,1)=0.65;
A(1,2)=0.25;
A(1,3)=0.4;

%prepare order parameter
%A(1,1)=0.7;
%A(1,2)=0.5;
%A(1,3)=0.3;
for k=1:NumberOfOP
   A(1,k)=rand;
end
%prepare attention parameter
for m=1:NumberOfOP
    lmd(m)=1;
end

%parameter
B=1;C=1;

 for n=1:NumberOfStep
    %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 = 1/D(n);
    %speed = 1;
    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
 end
 
%A is the order parameter,draw it;
plot(A)
axis([1,NumberOfStep,0,1.1])
%plot(D)