chp6ex9.m

solve power flows by gauss seidel method ans newton Rapson

clear
basemva = 100;  accuracy = 0.001; accel = 1.8; maxiter = 100;

%        IEEE 30-BUS TEST SYSTEM (American Electric Power)
%        Bus Bus  Voltage Angle   ---Load---- -------Generator----- Static Mvar
%        No  code Mag.    Degree  MW    Mvar  MW  Mvar Qmin Qmax    +Qc/-Ql
busdata=[1   1    1      0.0     0       0.0    0.0  0.0   0   0       0
         2   0    1      0.0     20     14.64   0  0.0 0  0       0
         3   0    1      0.0     25     18.5    0.0  0.0   0   0       0
         4   0    1.0    0.0     15     9.3    0.0  0.0   0   0       0
         5   0    1.0    0.0     17    12.75    0.0  0.0 0  0       0
         6   0    1.0     0.0     20    9.69    0.0  0.0   0   0       0
         7   0    1.0     0.0    18     13.5    0.0  0.0   0   0       0
       ];

%                                        Line code
%         Bus bus   R      X     1/2 B   = 1 for lines
%         nl  nr  p.u.   p.u.   p.u.     > 1 or < 1 tr. tap at bus nl
linedata=[1   2  5.38E-02	4.86E-02	6.49E-07     1
          1   3    3.05E-02	3.75E-02	5.28E-07  1
          1   4    5.70E-02	5.15E-02	6.88E-07   1
          1   5  5.38E-02	4.86E-02	6.49E-07   1
          1   6   6.01E-02	5.42E-02	7.25E-07    1
          1   7    3.80E-02	3.43E-02	4.59E-07  1
       ];

lfybus                            % form the bus admittance matrix
lfgauss                % Load flow solution by Gauss-Seidel method
busout              % Prints the power flow solution on the screen
lineflow          % Computes and displays the line flow and losses