multil.m

来自「《面向MATLAB工具箱的神经网络理论与应用》一书中实例的所有源代码」· M 代码 · 共 26 行

function Auc=multi1(Aeec)               % 函数名与输入、输出变量名multil.m、Aeec、uc

%
      [ ps,pq] = size(Aeec);            % 求出输入的行和列数
       pn = ps/2;                       % 求出两输入变量的个数
       pe = [ eye(pn) zeros(pn)]*Aeec;  % 制作两两相乘的矩阵
       pec = [ zeros(pn) eye(pn)]*Aeec;
       pec = pec';
 pa = pe(:,1)*pec(1,:);                 % 计算矩阵的乘积
 Auc = pa(1,:);                         % 制作Auc矩阵
 for i = 2:pn
     Auc = [Auc pa(i,:)];
 end
 
 for j = 2:pq                           %  制作Aucj矩阵
        pa = pe(:,j)*pec(j,:);
        Aucj = pa(1,:);
      for i = 2:pn
          Aucj = [Aucj pa(i,:)];
      end
    Auc =[Auc;Aucj];
 end    
 
 Auc = Auc';                            % 最后的结果