Scribd
Upload
12 views

Ex No:4 Load Flow Analysis by Gauss-Siedel Method

The document describes a load flow analysis using the Gauss-Siedel method. It provides sample input data for a 5 bus system including bus details, line impedances, real and reactive power values. It then forms the y-bus matrix and performs one iteration of the Gauss-Siedel load flow algorithm, providing the final calculated bus voltages as output.

Uploaded by

Syed Asif
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
12 views

Ex No:4 Load Flow Analysis by Gauss-Siedel Method

The document describes a load flow analysis using the Gauss-Siedel method. It provides sample input data for a 5 bus system including bus details, line impedances, real and reactive power values. It then forms the y-bus matrix and performs one iteration of the Gauss-Siedel load flow algorithm, providing the final calculated bus voltages as output.

Uploaded by

Syed Asif
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
You are on page 1/ 3

EX NO:4

LOAD FLOW ANALYSIS BY GAUSS-SIEDEL METHOD

% to form the y-bus matrix


disp(' THE INPUT DATAS ARE AS FOLLOWS');
nb=input('THE NUMBER OF BUSES\n');
nl=input('THE NUMBER OF LINES\n');
sb=input('THE NUMBER OF STARTING BUSES\n');
eb=input('THE NUMBER OF ENDING BUSES\n');
sli=input('ENTER THE DETAILS OF SERIES LINE IMPEDANCE\n');
lca=input('ENTER THE DETAIL OF LINE IMPEDANCE\n');
% execution starts here
y=zeros(nb,nb);
for m=1:nb
for n=1:nb
if m==n||m~=n
if sli(m,n)==0
n=n+1;
else
y(m,m)=y(m,m)+sli(m,n)^-1+lca(m,n);
y(m,n)=-sli(m,n)^-1;
y(n,m)=y(m,n);
end
end
end
end
%output of y-bus
disp(' THE RESULT OF Y BUS MATRIX is ');
ybus=y
m=input('ENTER THE NUMBER OF PV BUSES');
for i=1:m
v(i+1)=input('ENTER THE PV BUS VOLTAGE');
qmax(i)=input('ENTER MAXIMUM LIMIT OF Q FOR PV BUS');
qmin(i)=input('ENTER MINIMUM LIMIT OF Q FOR PV BUS');
end
p=input('ENTER THE REAL POWERS');
q=input('ENTER THE REACTIVE POWERS');
v(1)=input('ENTER THE SLACK BUS VOLTAGE');
itrn=input('ENTER NUMBER OF ITERATION');
for i=m+2:nb
v(i)=1;
end
for i=2:m+1
yv=0;
for j=1:nb
yv=yv+y(i,j)*v(j);
end

yv=(conj(v(i))*yv);
q(i)=-imag(yv);
if q(i)<qmin(i-1)
q(i)=qmin(i-1);
else
q(i)>qmax(i-1)
q(i)=qmax(i-1)
end
yv=0;
for j=1:nb
if i~=j
yv=yv+y(i,j)*v(j);
end
end
v1=(complex(p(i),-q(i))/(conj(v(i)))-yv)/y(i,i);
del(i)=angle(v1);
[re,im]=pol2cart(del(i),v(i));
v(i)=complex(re,im);
end
for i=m+2:nb
yv=0;
for j=1:nb
if i~=j
yv=yv+y(i,j)*v(j);
end
end
v(i)=(complex(-p(i),q(i))/(conj(v(i)))-yv)/y(i,i);
end
for k=2:itrn
for i=2:nb-1
yv=0;
for j=1:nb
if i~=j
yv=y(i,j)*v(j);
end
end
v(i)=(complex(-p(i).q(i))/(conj(v(i)))-yv)/y(i,i);
end
end
disp(['the final value of voltage after the iteration',num2str(itrn),'are']);
disp(v);

OUTPUT:
THE INPUT DATAS ARE AS FOLLOWS
THE NUMBER OF BUSES
5

THE NUMBER OF LINES


7
THE NUMBER OF STARTING BUSES
[1 1 2 2 2 3 4]
THE NUMBER OF ENDING BUSES
[2 3 3 4 5 4 5]
ENTER THE DETAILS OF SERIES LINE IMPEDANCE
[0 0.2+0.05j 0.8+0.24j 0 0
0.2+0.05j 0 0.06+0.18j 0.02+0.08j 0.04+0.12j
0.8+0.24j 0.06+0.18j 0 0.01+0.03j 0
0 0.02+0.08j 0.01+0.03j 0 0.08+0.24j
0 0.04+0.12j 0 0.08+0.24j 0]
ENTER THE DETAIL OF LINE IMPEDANCE
[0 0.06j 0.05j 0 0
0.06j 0 0.04j 0.04j 0.03j
0.05j 0.04j 0 0.02j 0
0 0.04j 0.02j 0 0.07j
0 0.03j 0 0.07j 0]
the result of y bus matrix is
y bus =
5.8527 - 1.4105i -4.7059 + 1.1765i -1.1468 + 0.3440i
0
0
-4.7059 + 1.1765i 11.8137 -25.2712i -1.6667 + 5.0000i -2.9412 +11.7647i -2.5000 + 7.5000i
-1.1468 + 0.3440i -1.6667 + 5.0000i 12.8135 -35.2340i -10.0000 +30.0000i
0
0
-2.9412 +11.7647i -10.0000 +30.0000i 14.1912 -45.3847i -1.2500 + 3.7500i
0
-2.5000 + 7.5000i
0
-1.2500 + 3.7500i 3.7500 -11.1500i
ENTER THE NUMBER OF PV BUSES 1
ENTER THE PV BUS VOLTAGE 1
ENTER MAXIMUM LIMIT OF Q FOR PV BUS 1.5
ENTER MINIMUM LIMIT OF Q FOR PV BUS -1.0
ENTER THE REAL POWERS [0 0.6 0.45 0.4 0.6]
ENTER THE REACTIVE POWERS [0 0 0.15 0.05 0.1]
ENTER THE SLACK BUS VOLTAGE 1.06
ENTER NUMBER OF ITERATION 1
THE FINAL VALUE OF VOLTAGE AFTER THE ITERATION1ARE
1.0600

1.0000 + 0.0021i 0.9960 - 0.0091i 0.9966 - 0.0141i 0.9826 - 0.0517i

You might also like