Psa 1
Psa 1
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 Yk 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.
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:
1 1 2 0.004 0.0533 0
5 4 5 0.006 0.08 0
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
the same. We also verified the results of the program with the manually
calculated values.