Load Flow Problem
Problem Statement: Solve the load flow problem below in Matlab, or similar program, using the Gauss-Seidel Method.
- Two graphs. One showing voltage magnitude 
at all buses vs. iteration number (k). A second showing voltage angle 
at all buses vs. iteration number (k).
- Your code. Please comment as much as possible.
- Any handwritten work you have. Again, please comment as much as possible to show me what you did.
The table below provides all of the values you need. Use a flat start initial condition.
|
all values in per unit (pu)
|
|
Component
|
|V|
|
theta
|
Pinj
|
Qinj
|
status
|
Gseries
|
Bseries
|
Bshunt
|
|
Bus #1 (infinite bus)
|
1.006
|
0.000
|
0.035
|
0.005
|
1
|
n/a
|
n/a
|
n/a
|
|
Bus #2 (voltage-controlled)
|
0.955
|
0.000
|
0.551
|
0.277
|
1
|
n/a
|
n/a
|
n/a
|
|
Bus #3 (load / PQ)
|
1.000
|
0.000
|
-0.205
|
0.145
|
1
|
n/a
|
n/a
|
n/a
|
|
X-Line 1-2 (pi)
|
n/a
|
n/a
|
n/a
|
n/a
|
1
|
0.365
|
-0.401
|
0.000
|
|
X-Line 2-3 (pi)
|
n/a
|
n/a
|
n/a
|
n/a
|
1
|
0.027
|
-0.322
|
0.000
|
|
X-Line 1-3 (pi)
|
n/a
|
n/a
|
n/a
|
n/a
|
1
|
0.328
|
-0.949
|
0.000
|
|
|
|
|
|
|
Y = G + jB
|