Power Systems Analysis Lab
BEEE306P
Digital Assignment - 2
Simulation of Y-bus matrix
Navya Kukreja
22BEE0093
Aim
To calculate the elements yij of the admittance-bus matrix from the given bus
network.
Theory
The nodal admittance matrix is given by the equation (for a three-bus system):
I1 Y11 Y12 Y13 V1
I2 = Y21 Y22 Y23 · V2
I3 Y31 Y32 Y33 V3
then, the nodal admittance matrix is
Y11 Y12 Y13
Y21 Y22 Y23
Y31 Y32 Y33
The diagonal elements are dependent on the given bus system and calculated
accordingly
1
Theoretical calculations
Given the system,
1 2 0.15 0.6 0.02 1
1 3 0.10 0.4 0.02 1
1 4 0.15 0.6 0.025 1
2 3 0.05 0.2 0.01 1
3 4 0.05 0.2 0.01 1
We have,
y12 = Z112 = 1
0.15+j0.6 = 0.39 − j1.56
y13 = Z113 = 1
0.10+j0.4 = 0.58 − j2.32
y14 = Z114 = 1
0.15+j0.6 = 0.39 − j1.56
y23 = Z123 = 1
0.05+j0.2 = 1.76 − j4.7
and diagonal elements are given by,
y11 = y12 + y13 + y14 + j 0.02 0.02 0.025
2 + j 2 + j 2 = 1.36 − j5.43
y22 = y21 + y23 + j 0.01 0.02
2 + j 2 = 1.56 − j6.25
y33 = y12 + y13 + y14 + j 0.02 0.02 0.025
2 + j 2 + j 2 = 1.36 − j5.43
y44 = y12 + y13 + y14 + j 0.02 0.02 0.025
2 + j 2 + j 2 = 1.36 − j5.43
Algorithm
1. User enters the line data matrix, which contains the start and end buses of
each transmission line, resistance (R), reactance (X), half-line charging suscep-
tance (Bc ), and transformation ratio (a).
2. The number of buses and transmission lines is determined. Using this in-
formation, an nbus × nbus zero matrix is initialized as Ybus , which will be populated
2
with admittance values representing the network’s electrical properties.
3. The off-diagonal elements of Ybus are updated by iterating through the trans-
mission lines. For each line, the mutual admittance between the connected buses
is calculated and added using the formula Ybus (nl , nr ) = −y/a. Symmetry is en-
sured by setting Ybus (nr , nl ) = Ybus (nl , nr ).
4. The diagonal elements of Ybus are updated by summing the contributions
from all connected lines for each bus. These include terms from the sending end,
y/a2 + Bc , and the receiving end, y + Bc . The total is added to the diagonal entry
Ybus (n, n).
5. Once all elements are calculated, the fully populated bus admittance matrix
Ybus is displayed. This matrix is essential for power system analysis and represents
the admittance relationships between buses in the network.
3
MATLAB code
clc
clear all
linedata = [1 2 0.15 0.6 0.02 1
1 3 0.10 0.4 0.02 1
1 4 0.15 0.6 0.025 1
2 3 0.05 0.2 0.01 1
3 4 0.05 0.2 0.01 1];
% Ybus formation
j = sqrt ( -1) ;
i = sqrt ( -1) ;
nl = linedata (: ,1) ;
nr = linedata (: ,2) ;
R = linedata (: ,3) ;
X = linedata (: ,4) ;
Bc = j * linedata (: ,5) ;
a = linedata (: ,6) ;
nbr = length ( linedata (: ,1) ) ;
nbus = max ( max ( nl ) , max ( nr ) ) ;
Z = R + j*X;
y = ones ( nbr , 1) ./ Z ;
for n = 1:1: nbr
if a ( n ) <= 0
a(n) = 1
else end
Ybus = zeros ( nbus , nbus ) ;
for k = 1:1: nbus
Ybus ( nl ( k ) , nr ( k ) ) = Ybus ( nl ( k ) , nr ( k ) ) - y ( k ) / a ( k ) ;
Ybus ( nr ( k ) , nl ( k ) ) = Ybus ( nl ( k ) , nr ( k ) ) ;
end
end
for n =1: nbr
for k = 1: nbr
if nl ( k ) == n
Ybus (n , n ) = Ybus (n , n ) + y ( k ) /( a ( k ) ^2) + Bc ( k ) ;
elseif nr ( k ) == n
Ybus (n , n ) = Ybus (n , n ) + y ( k ) + Bc ( k ) ;
else
end
end
4
end
disp (" The Ybus matrix is : ") ;
disp ( Ybus ) ;
Results
Output from given input parameters are,
”The Ybus matrix is:
Columns 1 through 3
1.3725 − 5.4252i −0.3922 + 1.5686i −0.5882 + 2.3529i
−0.3922 + 1.5686i 1.5686 − 6.2445i −1.1765 + 4.7059i
−0.5882 + 2.3529i −1.1765 + 4.7059i 2.9412 − 11.7247i
−0.3922 + 1.5686i 0.0000 + 0.0000i 0.0000 + 0.0000i
Column 4
−0.3922 + 1.5686i
0.0000 + 0.0000i
0.0000 + 0.0000i
1.5686 − 6.2395i
Summary
We have successfully calculated the admittance matrix for the given bus system
and output is consistent with theoretical calculations.
