inputdata_net2.m

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

function [SB,NodeNum,para_line,para_trans2,para_trans3,para_capac,para_bus]=InputData_Net2
% InputData_Net2 records the parameter of line,two/three windings transformer ,capacitor
% and bus.
% In this program,1 represents power transistion line,2 represents two
% winding transformer,n1 is higher voltage side and n2 is lower voltage 
% side of two winding transformer,3 represents three winding transformer,
% 4 represents capacitor,5/6/7 respectively represents PQ/PV/UA bus

SB=1000;
NodeNum=9;
% the input format of line's parameter is defined as follows:
% para_line=[1,n1,n2, r_unit,  x_unit,   b_unit,    Len, UB]
para_line  =[1 2 3    0.053    0.043     2.94e-6    50  220;...
             1 2 6    0.053    0.043     2.94e-6    60  220;...
             1 3 6    0.053    0.043     2.94e-6    30  220;...
             1 4 7    0.053    0.043     2.94e-6    80  110];

step1=2.5/100;
Umin1=220*(1-2*step1); Umax1=220*(1+2*step1);Uset1=220*(1-2*step1);
step2=1.25/100;
Umin2=110*(1-8*step2); Umax2=110*(1+8*step2);Uset2=110*(1+8*step2);
% the input format of two winding transformer parameter is defined as follows:
% para_trans2=[2,n1,n2, Pk,     Uk,   P0,     I0,  Sn,    [U1  U2],     [UB1 UB2],[Umin   Umax],step]
para_trans2  =[2  2  1  665.2   14   189.6  0.299  246  [Uset1  15.75]  [220  15] [Umin1 Umax1],step1;...
               2  7  8  153.08 10.67 28.68  0.16   40   [Uset2  10.5]   [110  10] [Umin2 Umax2],step2];    


step3=1.25/100;
Umin3=220*(1-8*step3); Umax3=220*(1+8*step3);Uset3=220*(1+8*step3);
% the input format of three winding transformer parameter is defined as follows:
% para_trans3=[3,n1,n2,n3,n0,[Pk12 Pk23 Pk31], [Uk12 Uk23 Uk31],P0,     I0,  [Sn1 Sn2 Sn3], [U1 U2 U3],    [UB1 UB2 UB3],[Umin  Umax],step]
para_trans3  =[3 3  4  5  9  [638.9 125 163.1] [12.9  8.2  23] 74.75 0.04444 [180 180 90] [Uset3 121 10.5] [220 110 10],[Umin3 Umax3],step3];

% the input format of capacitor parameter is defined as follows:
% para_capac=[4,n,Qcr,[Qcrmin Qcrmax],step]
step=[];
Qcrmin=[]; Qcrmax=[];
para_capac=[];

% the input format of bus's parameter is defined as follows:
% para_bus=[type, n, PG,  QG,   PD,  QD,  U,   theta,   [Umin       Umax],          UB ]
para_bus  =[  6   1  150   0    0    0   15.75   0    [15.75*0.95 15.75*1.05]       15 ;...
              5   2   0    0   150  100  220     0    [220*1        220*1.05]      220 ;...
              5   3   0    0    0    0   220     0    [220*0.95     220*1.05]      220 ;...
              5   4   0    0    0    0   121     0    [121*1        121*1.05]      110 ;...
              5   5   0    0   15    5   10.5    0    [10*1          10*1.05]       10 ;...
              7   6   0    0    0    0   220     0    [220*0.95     220*1.05]      220 ;...
              5   7   0    0    0    0   110     0    [110*0.95     110*1.05]      110 ;...
              5   8   0    0  100   80   10.5    0    [10*1          10*1.05]       10 ;...
              5   9   0    0    0    0   10.5    0    [121*1        121*1.05]       10 ];
   
% the input format of bus's parameter is defined as follows:
% para_gener=[8,n,[PGmin PGmax],[QGmin QGmax],a,b,c]