schema_art2_f1.m

ART2网络的MATLAB程序

function [u,p]=Schema_ART2_F1(I,F2_feedback,a,b,theta,e,m,NumNeurons)
    
    u = zeros(1,m);                  %【0 0 0 0 0】
    p = zeros(1,m);                  %【 0 0 0 0 0】同 u
    x = zeros(1,m);                  %【0 0 0 0 0】
    w = zeros(1,m);                  %【0 0 0 0 0】
    q = zeros(1,m);                  %【 0 0 0 0 0】同 u
    v = zeros(1,m);                   %【 0 0 0 0 0】同 u
  
    
       
       temp = ones(1,m);                            %对每一个新的输入，temp清为【1，1，1，1，1】
       
       %这个运算的稳定结果是u全部接近1？？？？？——》连续两次计算相差不大
       while sum(abs(temp-u)) >= 0.001        %u 不是全部为1， 则进行while内运算
          y = zeros(NumNeurons,1);               %假定一开始，F2层输出为0
          temp = u;                                       %保留原来的u？

          w = I + a * u;                                   % 从 u和输入I,获得w
          p = u+F2_feedback;                        %p为什么不考虑F2的返回？？？
          x = w / ( e + norm(w,2) );                 % 从w获得x
          q = p / ( e + norm(p,2) );                  %从p获得q
         
          for j = 1:m
             % 反映 f(x(k))
             if x(j) < theta
                x(j) = 0;
             else
                x(j) = x(j);
             end
             % 反映 f(q(k))
             if q(j) < theta
                q(j) = 0;
             else
                q(j) = q(j);
             end
          end     
          
          v = x + b * q;
          u = v / ( e + norm(v,2) );                 %从v 获得 u
       end
       %到此获得稳定的u, v等
end