newnode.m

这个程序是采用模块化的方式对潮流进行求解

   %首先依据bus确定系统的PQ节点、PV节点、平衡节点数分别为nPQ、nPV和nSW，
   %并分别将PQ节点、PV节点、平衡节点的数据分别放入数组PQ、PV和SW中
   function newnode(bus,line)
   global G B V delta P Q  n m nl l my bus line nodenum
   [nb,mb]=size(bus);
   [nl,ml]=size(line);
  nSW = 0;                 % number of swing bus counter
  nPV = 0;                 % number of PV bus counter
  nPQ = 0;                 % number of PQ bus counter
  for i = 1:nb,              % nb为总节点数
    type= bus(i,6);
    if type == 3,
       nSW = nSW + 1;     % increment swing bus counter
       SW(nSW,:)=bus(i,:);
    elseif type == 2,
        nPV = nPV +1;      % increment PV bus counter
       PV(nPV,:)=bus(i,:);
    else
        nPQ = nPQ + 1;     % increment PQ bus counter
       PQ(nPQ,:)=bus(i,:);
    end
  end; 
  
  n=nb;   
  m=nPQ;
  %按照PQ节点、PV节点、平衡节点的次序形成bus数据，并建立节点编号的对应关系： 
  bus=[PQ;PV;SW];
  newbus=[1:nb]';
  nodenum=[newbus bus(:,1)];
  bus(:,1)=newbus; 

  %依据新的节点编号修改line数据：
  for i=1:nl
      for j=1:2
          for k=1:nb
              if line(i,j)==nodenum(k,2)
                 line(i,j)=nodenum(k,1);
                 break
              end
          end
      end
  end