Scribd
Upload
64 views

Lab-01 Formation of Y-Bus: Date-21/07/2014 Programme Code

The document describes a MATLAB code that calculates the bus admittance matrix (Y-bus) for a power network. It takes in the line data including the from and to bus numbers, resistance, reactance, and charging admittance for each line. It then constructs the Y-bus matrix with the appropriate admittance values based on the line data. The code is then tested on a 4 bus sample network and the resulting Y-bus matrix is displayed.

Uploaded by

Gaurav Singh
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
64 views

Lab-01 Formation of Y-Bus: Date-21/07/2014 Programme Code

The document describes a MATLAB code that calculates the bus admittance matrix (Y-bus) for a power network. It takes in the line data including the from and to bus numbers, resistance, reactance, and charging admittance for each line. It then constructs the Y-bus matrix with the appropriate admittance values based on the line data. The code is then tested on a 4 bus sample network and the resulting Y-bus matrix is displayed.

Uploaded by

Gaurav Singh
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Lab-01

Formation of Y-bus
Date-21/07/2014

Programme code:
clc
clear all
nbranch=input('Enter total no. of lines in network');
disp('Enter line data')
for p=1:1:nbranch
fb=input('From bus no. :');
tb=input('To bus no. :');
r=input('Resistance of line =');
x=input('Reactance of line =');
B=input('Line Charging Admittance"b/2"=');
T=input('Tap=');
z = r + 1i*x;
y = 1./z;
Ldata(p,:)=[fb tb r x B T y];
end
fb = Ldata(:,1);
tb = Ldata(:,2);
r = Ldata(:,3);
x = Ldata(:,4);
b = Ldata(:,5);
a = Ldata(:,6);
y = Ldata(:,7);
b = 1i*b;
nbus = max(max(fb),max(tb));
Y = zeros(nbus,nbus);
for k=1:nbranch
Y(fb(k),tb(k)) = Y(fb(k),tb(k))-y(k);
Y(tb(k),fb(k)) = Y(fb(k),tb(k));
end
for m =1:nbus
for n =1:nbranch
if fb(n) == m
Y(m,m) = Y(m,m) + y(n) + b(n);
elseif tb(n) == m
Y(m,m) = Y(m,m) + y(n) + b(n); end
end
end

Data :

enter total no. of lines in network 4


enter line data
from bus no. :1
to bus no. :2
resistance of line =0.01008
reactance of line =0.0504
line charging admittance"b/2"=0.05125
tap=0
from bus no. :1
to bus no. :3
resistance of line =0.00744
reactance of line =0.0372
line charging admittance"b/2"=0.03875
tap=0
from bus no. :2
to bus no. :4
resistance of line =0.00744
reactance of line =0.0372
line charging admittance"b/2"=0.03875
tap=0
from bus no. :3
to bus no. :4
resistance of line =0.01272
reactance of line =0.0636
line charging admittance"b/2"=0.06375
tap=0
>> Y

Result:

Y=
8.9852 -44.8360i -3.8156 +19.0781i -5.1696 +25.8478i
0
-3.8156 +19.0781i 8.9852 -44.8360i
0
-5.1696 +25.8478i
-5.1696 +25.8478i
0
8.1933 -40.8638i -3.0237 +15.1185i
0
-5.1696 +25.8478i -3.0237 +15.1185i 8.1933 -40.8638i
>>

You might also like