Chp9ex9.m
上传用户:eighthdate
上传日期:2014-05-24
资源大小:270k
文件大小:2k
源码类别:

其他行业

开发平台:

Matlab

  1. disp('Example 9.9')
  2. basemva = 1000;  accuracy = 0.0001;  maxiter = 10;
  3. % The impedances are expressed on a 1000 MVA base.
  4. % In problems 9.7-9.9 the base is mistakenly stated as 100 MVA. 
  5. %       Bus  Bus   |V|  Ang   ---Load---      ---Gen---   Gen Mvar Injected
  6. %       No.  code p.u.  Deg   MW     Mvar   MW   Mvar   Min   Max    Mvar
  7. busdata=[1    1   1.04   0    00.0    0.0     0.0   0.0    0     0    0
  8.          2    0   1.0    0    00.0    0.0     0.0   0.0    0     0    0
  9.          3    0   1.0    0   150.0  120.0     0.0   0.0    0     0    0
  10.          4    0   1.0    0     0.0    0.0     0.0   0.0    0     0    0
  11.          5    0   1.0    0   120.0   60.0     0.0   0.0    0     0    0
  12.          6    0   1.0    0   140.0   90.0     0.0   0.0    0     0    0
  13.          7    0   1.0    0     0.0    0.0     0.0   0.0    0     0    0
  14.          8    0   1.0    0   110.0   90.0     0.0   0.0    0     0    0
  15.          9    0   1.0    0    80.0   50.0     0.0   0.0    0     0    0
  16.          10   2   1.035  0     0.0    0.0   200.0   0.0    0    180   0
  17.          11   2   1.03   0     0.0    0.0   160.0   0.0    0    120   0];
  18. %        Bus   Bus     R        X       1/2B
  19. %        No.   No.    p.u.     p.u.     p.u.
  20. linedata=[1     2     0.00     0.06     0.0000    1
  21.           2     3     0.08     0.30     0.0004    1
  22.           2     6     0.12     0.45     0.0005    1
  23.           3     4     0.10     0.40     0.0005    1
  24.           3     6     0.04     0.40     0.0005    1
  25.           4     6     0.15     0.60     0.0008    1
  26.           4     9     0.18     0.70     0.0009    1
  27.           4    10     0.00     0.08     0.0000    1
  28.           5     7     0.05     0.43     0.0003    1
  29.           6     8     0.06     0.48     0.0000    1
  30.           7     8     0.06     0.35     0.0004    1
  31.           7    11     0.00     0.10     0.0000    1
  32.           8     9     0.052    0.48     0.0000    1];
  33. %         Gen.  Ra     Xd'
  34. gendata=[ 1     0     0.20
  35.           10    0     0.15
  36.           11    0     0.25];
  37.       lfybus                              % Forms the bus admittance matrix
  38.       lfnewton               % Power flow solution by Newton-Raphson method
  39.       busout                 % Prints the power flow solution on the screen
  40.       Zbus=zbuildpi(linedata, gendata, yload)%Forms Zbus including the load
  41.       symfault(linedata, Zbus, V)    % 3-phase fault including load current