FORMATION OF NETWORK MATRICES Z BUS
EX.NO:2
DATE:
AIM
To develop a OCTAVE program to obtain bus impedance matrix Z-bus of the given
power system
SOFTWARE REQUIRED
OCTAVE
THEORY
Z-bus matrix is an important matrix used in different kinds of power system study
such as short circuit study, load flow study etc. In short circuit analysis the generator
uses transformer impedance must be taken into account. In quality analysis the two-
short element are neglected by forming the z-bus matrix which is used to compute the
voltage distribution factor. This can be largely obtained by reversing the y-bus formed
by resection method or by analytical method.
Taking inverse of the y-bus for large system in time conditioning managing
modification in the system requires whole process to be repeated to voltage changes in
the system. In such cases z-bus computed to z-bus solving algorithm.
ALGORITHM
Step 1: Zhis added from a new bus to the reference bus (ie. a new branch is added and
the dimension of Z-bus goes up by one). This is type-l modification.
Step 2: Zo is added from a new bus to an old bus (ie., a new branch is added and
The dimension of Z-bus goes up by one). This is type-2 modification.
Step 3: Z, connects an old bus to the reference branch (ie.a new loop is formed but the
dimension of z-bus does not change). This is type-3 modification
Step 4: Z connects two old buses (ie.,new loop is formed but the dimension of Z-bus,
does not change). This is type-4 modification.
Step 5: Z connects two new buses (Zbus remains unaffected in this case). This
situation can be avoided by suitable numbering of buses.
Step 6 Print the z-bus matrix.
PROBLEM
Find the ZBus for the three bus network shown in figure
PROGRAM:
OCTAVE Program for formation of Z bus matrix:
clear
zprimary=[
1100.25
22101
33 10.1
4200.25
5230.11
[elements coloumns)=size(zprimary)
zbus=[]
currentbusno=0
for count=1 elements,
[rows cols=size(zbus)
from=zprimary(count,2)
to-zprimary(count,3)
value=zprimary(count,4)
newbus=max(from, to) ref=min(from to)
if newbus currentbusno&ref==0
zbus=[zbus zeros(rows,1)
zeros(1.cols) value]
currentbusno=newbus
continue
