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

其他行业

开发平台:

Matlab

  1. clear
  2. basemva = 100;  accuracy = 0.001; maxiter = 12;
  3. %        IEEE 30-BUS TEST SYSTEM (American Electric Power)
  4. %        Bus Bus  Voltage Angle   ---Load---- -------Generator----- Injected
  5. %        No  code Mag.    Degree  MW    Mvar  MW  Mvar Qmin Qmax     Mvar
  6. busdata=[1   1    1.06    0.0     0.0   0.0    0.0  0.0   0   0       0
  7.          2   2    1.0     0.0   21.70  12.7   40.0 48.822 0   0       0
  8.          3   0    1.0     0.0     2.4   1.2    0.0  0.0   0   0       0
  9.          4   0    1.06    0.0     7.6   1.6    0.0  0.0   0   0       0
  10.          5   2    1.0     0.0    94.2  19.0    0.0 35.975 0   0       0
  11.          6   0    1.0     0.0     0.0   0.0    0.0  0.0   0   0       0
  12.          7   0    1.0     0.0    22.8  10.9    0.0  0.0   0   0       0
  13.          8   2    1.0     0.0    30.0  30.0    0.0 30.826 0   0       0
  14.          9   0    1.0     0.0     0.0   0.0    0.0  0.0   0   0       0
  15.         10   0    1.0     0.0     5.8   2.0    0.0  0.0   0   0      19
  16.         11   2    1.0     0.0     0.0   0.0    0.0 16.119 0   0       0
  17.         12   0    1.0     0       11.2  7.5    0    0     0   0       0
  18.         13   2    1.0     0        0    0.0    0   10.423 0   0       0
  19.         14   0    1       0       6.2   1.6    0    0     0   0       0
  20.         15   0    1       0       8.2   2.5    0    0     0   0       0
  21.         16   0    1       0       3.5   1.8    0    0     0   0       0
  22.         17   0    1       0       9.0   5.8    0    0     0   0       0
  23.         18   0    1       0       3.2   0.9    0    0     0   0       0
  24.         19   0    1       0       9.5   3.4    0    0     0   0       0
  25.         20   0    1       0       2.2   0.7    0    0     0   0       0
  26.         21   0    1       0      17.5  11.2    0    0     0   0       0
  27.         22   0    1       0       0     0.0    0    0     0   0       0
  28.         23   0    1       0       3.2   1.6    0    0     0   0       0
  29.         24   0    1       0       8.7   6.7    0    0     0   0      4.3
  30.         25   0    1       0       0     0.0    0    0     0   0       0
  31.         26   0    1       0       3.5   2.3    0    0     0   0       0
  32.         27   0    1       0       0     0.0    0    0     0   0       0
  33.         28   0    1       0       0     0.0    0    0     0   0       0
  34.         29   0    1       0       2.4   0.9    0    0     0   0       0
  35.         30   0    1       0      10.6   1.9    0    0     0   0       0];
  36. %                                        Line code
  37. %         Bus bus   R      X     1/2 B   = 1 for lines
  38. %         nl  nr  p.u.   p.u.   p.u.     > 1 or < 1 tr. tap at bus nl
  39. linedata=[1   2   0.0192   0.0575   0.02640    1
  40.           1   3   0.0452   0.1852   0.02040    1
  41.           2   4   0.0570   0.1737   0.01840    1
  42.           3   4   0.0132   0.0379   0.00420    1
  43.           2   5   0.0472   0.1983   0.02090    1
  44.           2   6   0.0581   0.1763   0.01870    1
  45.           4   6   0.0119   0.0414   0.00450    1
  46.           5   7   0.0460   0.1160   0.01020    1
  47.           6   7   0.0267   0.0820   0.00850    1
  48.           6   8   0.0120   0.0420   0.00450    1
  49.           6   9   0.0      0.2080   0.0    0.978
  50.           6  10   0         .5560   0      0.969
  51.           9  11   0         .2080   0          1
  52.           9  10   0         .1100   0          1
  53.           4  12   0         .2560   0      0.932
  54.          12  13   0         .1400   0          1
  55.          12  14    .1231    .2559   0          1
  56.          12  15    .0662    .1304   0          1
  57.          12  16    .0945    .1987   0          1
  58.          14  15    .2210    .1997   0          1
  59.          16  17    .0824    .1923   0          1
  60.          15  18    .1073    .2185   0          1
  61.          18  19    .0639    .1292   0          1
  62.          19  20    .0340    .0680   0          1
  63.          10  20    .0936    .2090   0          1
  64.          10  17    .0324    .0845   0          1
  65.          10  21    .0348    .0749   0          1
  66.          10  22    .0727    .1499   0          1
  67.          21  22    .0116    .0236   0          1
  68.          15  23    .1000    .2020   0          1
  69.          22  24    .1150    .1790   0          1
  70.          23  24    .1320    .2700   0          1
  71.          24  25    .1885    .3292   0          1
  72.          25  26    .2544    .3800   0          1
  73.          25  27    .1093    .2087   0          1
  74.          28  27     0       .3960   0      0.968
  75.          27  29    .2198    .4153   0          1
  76.          27  30    .3202    .6027   0          1
  77.          29  30    .2399    .4533   0          1
  78.           8  28    .0636    .2000   0.0214     1
  79.           6  28    .0169    .0599   0.065      1];
  80. lfybus                            % form the bus admittance matrix
  81. perturb          % Load flow solution power perturbation technique
  82. busout              % Prints the power flow solution on the screen
  83. lineflow          % Computes and displays the line flow and losses