initbus.m

％ 该Matlab程序基于牛顿－拉夫逊算法

function [PQ,PV,UA]=InitBus(para_bus,SB)
% This program was used for initialize the bus which input format is defined as follows:
% para_bus=[type,n,PG,  QG,  PD,  QD,  U,   theta,   [Umin Umax],           UB   ]

k1=1;k2=1;k3=1;
for n=1:size(para_bus,1)
    if para_bus(n,1)==5%PQ bus
       PQ(1,k1)=para_bus(n,2);
       PQ(2,k1)=(para_bus(n,5)-para_bus(n,3)+(para_bus(n,6)-para_bus(n,4))*i)/SB;
       k1=k1+1;
    end
    if para_bus(n,1)==6%PV bus
       PV(1,k2)=para_bus(n,2);
       PV(2,k2)=(para_bus(n,5)-para_bus(n,3))/SB;
       PV(3,k2)=para_bus(n,7)/para_bus(n,11);
       k2=k2+1;
    end
    if para_bus(n,1)==7%UA bus
       UA(1,k3)=para_bus(n,2);
       UA(2,k3)=para_bus(n,7)/para_bus(n,11);
       k3=k3+1;
    end
end