Scribd
Upload
161 views

Loadf Flow Studies Using Gauss Siedel

For a 4 bus power system with generators at all buses and loads at buses 2 and 3, the document provides the bus admittance matrix and uses the Gauss-Siedel method in MATLAB to calculate the voltages and angles at each bus over iterations until the solution converges. The MATLAB code performs the load flow calculation using the Gauss-Siedel method and outputs the final voltages after convergence within 0.00001 p.u.

Uploaded by

durgendra
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
161 views

Loadf Flow Studies Using Gauss Siedel

For a 4 bus power system with generators at all buses and loads at buses 2 and 3, the document provides the bus admittance matrix and uses the Gauss-Siedel method in MATLAB to calculate the voltages and angles at each bus over iterations until the solution converges. The MATLAB code performs the load flow calculation using the Gauss-Siedel method and outputs the final voltages after convergence within 0.00001 p.u.

Uploaded by

durgendra
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/ 6

Load Flow Studies using Gauss Siedel

method
Aim: To find load flow solution of the given power system using GaussSeidel method theoretically for one iteration and obtain full solution
using MATLAB.
Apparatus: Matlab R2007 or above versions
Algorithm:
Problem:
1. For the sample power system shown below, the generators are
connected at all the four buses, while loads are at buses 2 and 3.
Values of real and reactive powers are listed in the table. All buses
other than the slack are PQ type. Assuming a flat voltage start, find
the voltages and bus angles at the three buses at the end of first GS
iteration.
2. Inputdata:

Bu
1s
2
3
4

Pi , pu
0.5
1.0
0.3

Qi , pu
0.5
0.2
0.1

Vi , pu
1.04 0 0
-

Remarks
Slack bus
PQ bus
PQ bus
PQ bus

For the above system, the bus admittance matrix is


Y=[ 3-j*9
-2+j*6

-2+j*6

-1+j*3

3.666-j*11 -0.666+j*2

-1+j*3 -0.666+j*2

3.666-j*11

-2+j*6

-1+j*3

0
-1+j*3
-2+j*6
3-j*9]

Write a MATLAB program to solve the load flow e q u a ti o n s of


the above sample power system by using Gauss-Seidal method.
MATLAB PROGRAM:
clc;
n=4;
V=[1.04 1 1 1];
Y=[3-j*9
-2+j*6
-1+j*3
0
-2+j*6
3.666-j*11
-0.666+j*2
-1+j*3
-1+j*3
-0.666+j*2
3.666-j*11
-2+j*6
0
-1+j*3
-2+j*6
3-j*9 ]
diff=10;
noofiter=1;
Vprev=V;
while(diff>0.00001|noofiter==1),
abs(V);
abs(Vprev);
Vprev=V;
P=[inf 0.5 -1 0.3];
Q=[inf -0.2 0.5 -0.1];
S=[inf+j*inf (0.5-j*0.2) (-1.0+j*0.5) (0.3-j*0.1)];
for i=2:n,
sumyv=0;
for k=1:n,
if (i~=k)
sumyv=sumyv+Y(i,k)*V(k);
end
end
V(i)=(1/Y(i,i))*((P(i)-j*Q(i))/conj(V(i))-sumyv);
end

V
diff=max(abs(abs(V(2:n))-abs(Vprev(2:n))));
noofiter=noofiter+1;
end
Output:
V=
1.0400 + 0.0000i 1.0191 + 0.0464i 1.0280 - 0.0870i 1.0250 - 0.0092i
1.0400 + 0.0000i 1.0290 + 0.0269i 1.0352 - 0.0934i 1.0334 - 0.0208i
1.0400 + 0.0000i 1.0335 + 0.0223i 1.0401 - 0.0999i 1.0385 - 0.0269i
1.0400 + 0.0000i 1.0360 + 0.0193i 1.0427 - 0.1034i 1.0413 - 0.0304i
1.0400 + 0.0000i 1.0374 + 0.0176i 1.0442 - 0.1053i 1.0429 - 0.0324i
1.0400 + 0.0000i 1.0382 + 0.0167i 1.0450 - 0.1064i 1.0437 - 0.0335i
1.0400 + 0.0000i 1.0386 + 0.0161i 1.0455 - 0.1070i 1.0442 - 0.0341i
1.0400 + 0.0000i 1.0388 + 0.0158i 1.0457 - 0.1074i 1.0445 - 0.0344i
1.0400 + 0.0000i 1.0390 + 0.0157i 1.0459 - 0.1076i 1.0446 - 0.0346i
1.0400 + 0.0000i 1.0390 + 0.0156i 1.0460 - 0.1077i 1.0447 - 0.0347i
1.0400 + 0.0000i 1.0391 + 0.0155i 1.0460 - 0.1078i 1.0448 - 0.0348i
1.0400 + 0.0000i 1.0391 + 0.0155i 1.0460 - 0.1078i 1.0448 - 0.0348i
1.0400 + 0.0000i 1.0391 + 0.0155i 1.0460 - 0.1078i 1.0448 - 0.0348i
1.0400 + 0.0000i 1.0391 + 0.0155i 1.0461 - 0.1078i 1.0448 - 0.0349i
Result: Load flow solution of the given power system using GaussSeidel method was obtained using using MATLAB.

You might also like