Initiation

 Subroutine Initiation read in the input full model data and convert the
 data format to be good to generate the full model admittance matrix.
 
 [NFROM,NTO,BraNum,LineB,ShuntB,BCIRC,BusNum,NUMB,SelfB,mpc,ExBus,newbusnum,oldbusnum] = Initiation(mpc,ExBus)

 INPUT DATA:
   mpc, struct, full model data in MATPOWER case format
   ExBus, n*1 vector, includes indices of external buses

 OUTPUT DATA
   NFROM, 1*n array, includes indices of all from end buses
   NTO, 1*n array, includes indices of all to end buses
   BraNum, 1*n array, includes indices of all branches
   LineB, 1*n array, includes line admittance of branches
   ShuntB, 1*n array, includes line shunt admittance of branches
   mpc, struct, full model in MATPOWER case format
   ExBus, n*1 vector, includes external bus indices
   newbusnum, 1*n array, indices of buses in internal bus numbering
   oldbusnum, 1*n array, indices of buses in original bus numbering
 NOTE:
   All output bus indices are in internal bus numbering.

0001 function [NFROM,NTO,BraNum,LineB,ShuntB,BCIRC,BusNum,NUMB,SelfB,mpc,ExBus,newbusnum,oldbusnum] = Initiation(mpc,ExBus)
0002 % Subroutine Initiation read in the input full model data and convert the
0003 % data format to be good to generate the full model admittance matrix.
0004 %
0005 % [NFROM,NTO,BraNum,LineB,ShuntB,BCIRC,BusNum,NUMB,SelfB,mpc,ExBus,newbusnum,oldbusnum] = Initiation(mpc,ExBus)
0006 %
0007 % INPUT DATA:
0008 %   mpc, struct, full model data in MATPOWER case format
0009 %   ExBus, n*1 vector, includes indices of external buses
0010 %
0011 % OUTPUT DATA
0012 %   NFROM, 1*n array, includes indices of all from end buses
0013 %   NTO, 1*n array, includes indices of all to end buses
0014 %   BraNum, 1*n array, includes indices of all branches
0015 %   LineB, 1*n array, includes line admittance of branches
0016 %   ShuntB, 1*n array, includes line shunt admittance of branches
0017 %   mpc, struct, full model in MATPOWER case format
0018 %   ExBus, n*1 vector, includes external bus indices
0019 %   newbusnum, 1*n array, indices of buses in internal bus numbering
0020 %   oldbusnum, 1*n array, indices of buses in original bus numbering
0021 % NOTE:
0022 %   All output bus indices are in internal bus numbering.
0023 
0024 %   MATPOWER
0025 %   Copyright (c) 2014-2016, Power Systems Engineering Research Center (PSERC)
0026 %   by Yujia Zhu, PSERC ASU
0027 %
0028 %   This file is part of MATPOWER/mx-reduction.
0029 %   Covered by the 3-clause BSD License (see LICENSE file for details).
0030 %   See https://github.com/MATPOWER/mx-reduction/ for more info.
0031 
0032 mpc.bus = sortrows(mpc.bus,1); % sort the buses
0033 oldbusnum = mpc.bus(:,1);
0034 newbusnum = [1:size(mpc.bus,1)]';
0035 mpc.bus(:,1) = newbusnum;
0036 mpc.branch(:,1) = interp1(oldbusnum,newbusnum,mpc.branch(:,1)); % change the branch terminal bus number
0037 mpc.branch(:,2) = interp1(oldbusnum,newbusnum,mpc.branch(:,2)); % change the branch terminal bus number
0038 mpc.gen(:,1)=interp1(oldbusnum,newbusnum,mpc.gen(:,1)); % change the generator bus number
0039 ExBus = interp1(oldbusnum,newbusnum,ExBus')';
0040 % bus data
0041 NUMB = newbusnum;
0042 BusNum = size(mpc.bus,1);
0043 SelfB = mpc.bus(:,6)./mpc.baseMVA;
0044 % branch data
0045 % status = mpc.branch(:,11);
0046 % delete all out of service branches
0047 % mpc.branch(status==0,:)=[];
0048 % status(status==0)=[];
0049 BraNum = size(mpc.branch,1);
0050 NFROM = mpc.branch(:,1);
0051 NTO = mpc.branch(:,2);
0052 LineB = 1./mpc.branch(:,4); % calculate the branch susceptance (b)
0053 ShuntB = mpc.branch(:,5)/2; % branch shunts;
0054 BCIRC = GenerateBCIRC(mpc.branch);
0055 % update SelfB
0056 for i = 1:BraNum
0057     SelfB(NFROM(i))=SelfB(NFROM(i))+LineB(i)+ShuntB(i);
0058     SelfB(NTO(i))=SelfB(NTO(i))+LineB(i)+ShuntB(i);
0059 end
0060 
0061 
0062 end