EXPERIMENT-1 Vinit Kumar(2K21/EE/312)
Aim: To compute the bus admittance matrix Ybus for small power system network given below using
MATLAB file.
Theory: The bus admittance matrix (Ybus) is a crucial component in power system analysis. It is used
to represent the electrical network in a simplified form and is particularly useful in solving power
flow and other analyses. The Ybus matrix relates the bus voltages and currents in a power system and
is an essential tool for studying the behaviour of the system. In a power system, the network is often
represented by buses and branches. Buses represent the connection points in the network, while
branches represent the transmission lines and other elements connecting these buses. The Ybus matrix
is derived from the admittances of the individual branches and the buses in the system. The general
form of the Ybus matrix is given by:
Ybus = Yshunt + ∑Yk
Where, Ybus is the bus admittance matrix, Yshunt represents the shunt admittances of the buses and k is
the admittance of the k-th branch. The shunt admittance of a bus is usually due to the presence of
loads or other devices connected to the bus. It is a combination of the load admittance and the shunt
capacitance of the transmission lines connected to the bus.
STEPS IN CONSTRUCTING THE Ybus BY INSPECTION:
1. The Ybus matrix is symmetric.
2. Yij, is the self-admittance (diagonal term), is equal to the sum of the primitive admittances of all
the components connected to the ith node.
3. Yij, the ijth element of the Ybus (off-diagonal element), is equal to the negative of the primitive
admittance of all components connected between nodes i and j. It is to be noted here that if more than
one component is connected in parallel between two nodes, the equivalent primitive admittance of the
components is first obtained before determining the entry in the Ybus.
Procedure:
Step-1: Enter in the branch data (i.e. R, X and B as given) in a matrix ‘D’. ‘D’ should contain five
columns and as many rows as that in the given branch data/ line data of given network. The five
columns correspond to 'From bus number', 'To bus number', R, X and B respectively. Check if data B
corresponds to Btotal or B1/2.
Step-2: Type the code given below save it in a MATLAB .m file. Save the code.
Step-3: Debug the code to generate Ybus.
Step-4: check the program results with the manually calculated values.
Code:
clc;
clear;
promptBuses = 'Enter number of buses : ';
nb = input(promptBuses) ;
z = inf (nb, nb) ;
promptArray = 'Enter the values of the array : ';
D = input (promptArray) ;
r = size (D);
rows = r(1) ;
for i=1:rows
z(D(i,1),D(i,2))=D(i,3)+j*D(i,4) ;
z(D(i,2),D(i,1))=z(D(i,1),D(i,2));
b(D(i,1),D(i,2))=D(i,5);
b(D(i,2),D(i,1))=b(D(i,1),D(i,2));
end
for i=1:nb
y0(i)=j*sum(b(i,:))/2;
end
for i=1:nb
for m=1:nb
if i~=m
Y(i,m)=-1/z(i,m);
else
Y(i,m)=y0(i)+sum(1./z(i,:));
end
end
end
disp(Y);
CIRCUIT DIAGRAM:
SINGLE LINE DIAGRAM FOR IEEE 5 BUS SYSTEM
SINGLE LINE DIAGRAM FOR IEEE 14 BUS SYSTEM
Calculations :
OBSERVATIONS:
For bus size = 5:
Line Data for IEEE 5 Bus System:
Output Ybus matrix For bus size = 14:
Line Data for IEEE 14 Bus System:
Output Ybus matrix :
CONCLUSION:
Through this experiment we have successfully learnt how to derive the Ybus matrix of any IEEE bus
system and developed a MATLAB program to derive 20 13 14 0.17093 0.34802 0 1 the same. We
also verified the results of the program with the manually calculated values.