By your instruction, the problem is smoothed. Thank you very much!
But I found an interesting thing, the PF result of 3.0b3 is slightly
different from that of 3.0.0 or 3.1b2. The attachment is my test system and
results of different versions. The only difference is the Q of bus 1. I set
the Qmax and Qmin of bus 1 to zero so that I think the Q calculated should
also be zero. But it is not zero by using 3.0b3 while it is zero by using
3.0.0 or 3.1b2.
By the way, I want to ask another question. If the initial values of Vg
should be 1 or other?
2007/4/27, Ray Zimmerman <[email protected]>:
When upgrading you should completely remove all files from the previous
version from your Matlab path. Then add all of the files from your new
version to your Matlab path.
If 'isload' is not defined it means you do not have the new version
installed correctly ... isload.m is one of the files in the new version
and should be in your path.
--
Ray Zimmerman
Senior Research Associate
428-B Phillips Hall, Cornell University, Ithaca, NY 14853
phone: (607) 255-9645
On Apr 27, 2007, at 1:27 AM, Jun Wu wrote:
Dear Ray Zimmerman,
Is the higher version compatible with the lower version well? I found some
problems. I used 3.0b3 before, which is lack of the new feature you've
referred. So I copy all files of 3.1b2 to replace them. But the same
program can not work with the error--Undefined function or method 'isload'
for input arguments of type 'double'. I don't know how update it simply
Yours sincerely,
Wu Jun
Unit: State Key Lab. of Power Systems,
Dep. of Electrical Engineering,
Tsinghua University, P.R. China.
Email: [email protected]
2007-04-27
2007/4/27, Ray Zimmerman <[email protected]>:
>
> On Apr 26, 2007, at 11:02 AM, Jun Wu wrote:
>
> Thanks for your attention and suggestion. I think reading the manual
> more will help me. But I never saw it before. Is this the latest feature for
> version 3.1b2? Does 3.0b3 has this setting?
>
>
>
> This feature was added in 3.0b4 and was included in the final 3.0.0release (
> 3.0b4 was a beta version). It is also included in the beta releases for
> 3.1.
>
> By the way, what do below settings mean?
>
> 38 - OUT_ALL_LIM, -1 control constraint info output
> [ -1 - individual flags control what constraint info
> prints]
> [ 0 - no constraint info (overrides individual flags)
> ]
> [ 1 - binding constraint info (overrides individual
> flags)]
> [ 2 - all constraint info (overrides individual flags)
> ]
> 39 - OUT_V_LIM, 1 control output of voltage limit info
>
> [ 0 - don't print
> ]
> [ 1 - print binding constraints only
> ]
> [ 2 - print all constraints
> ]
> [ (same options for OUT_LINE_LIM, OUT_PG_LIM, OUT_QG_LIM)
> ]
> 40 - OUT_LINE_LIM, 1 control output of line limit info
> 41 - OUT_PG_LIM, 1 control output of gen P limit info
> 42 - OUT_QG_LIM, 1 control output of gen Q limit info
>
>
>
> These options control the constraint sections of the output for the OPF
> solver only.
>
>
> Ray
>
> 2007/4/26, Ray Zimmerman <[email protected]>:
> >
> > By default, the power flow solver does not take into account the
> > voltage or reactive power limits. However [1] ...
> >
> > For the AC power flow solvers, if the ENFORCE_Q_LIMS option is set to
> > true (default is false), then if any generator reactive power limit is
> > violated after running the AC power flow, the corresponding bus is con-
> > verted to a PQ bus, with the reactive output set to the limit, and the case
> > is re-run. The voltage magnitude at the bus will deviate from the specified
> > value in order to satisfy the reactive power limit. If the generator at the
> > reference bus is reaches a reactive power limit and the bus is converted to
> > a PQ bus, the first re- maining PV bus will be used as the slack bus for the
> > next iteration. This may result in the real power out- put at this generator
> > being slightly off from the specified values.
> >
> >
> > [1] Taken from p. 10 of the User's Manual
(http://www.pserc.cornell.edu/matpower/manual.pdf
> > ).
> >
> > --
> > Ray Zimmerman
> > Senior Research Associate
> > 428-B Phillips Hall, Cornell University, Ithaca, NY 14853
> > phone: (607) 255-9645
> >
> >
> > On Apr 26, 2007, at 5:23 AM, Wu Jun wrote:
> >
> > Hi all,
> >
> >
> > From the case30.m in Matpower, Vmax and Vmin for buses are listed,
> > also the Qmax and Qmin for generators. I want to know how they work during
> > the process of calculating the power flow. If one of them are violated, how
> > does Matpower treat it? Does Matpower keep the original outputs or change
> > the calculation condition automatically then put out a result without
> > breaking any limit?
> >
> >
> > --------------
> > Yours sincerely,
> > Wu Jun
> > Unit: State Key Lab. of Power Systems,
> > Dep. of Electrical Engineering,
> > Tsinghua University, P.R. China.
> > Email: [email protected]
> > 2007-04-26
> >
> >
> >
> >
> >
> >
> >
>
>
>
Newton's method power flow converged in 4 iterations.
Converged in 0.11 seconds
================================================================================
| System Summary |
================================================================================
How many? How much? P (MW) Q (MVAr)
--------------------- -------------------- -------- -----------------
Buses 30 Total Gen Capacity 890.5 -75.0 to +257.5
Generators 6 On-line Capacity 890.5 -75.0 to +257.5
Committed Gens 6 Generation (current) 290.7 124.0
Loads 21 Load 283.4 126.2
Branches 41 Losses (I^2 * Z) 7.33 32.46
Transformers 4 Branch Charging (inj) - 34.4
Areas 3 Shunt (inj) 0.0 0.2
Inter-ties 7 Total Inter-tie Flow 105.2 39.7
Minimum Maximum
------------------------- --------------------------------
Voltage Magnitude 0.985 p.u. @ bus 26 1.088 p.u. @ bus 13
Voltage Angle -12.66 deg @ bus 30 0.00 deg @ bus 1
P Losses (I^2*R) - 1.52 MW @ line 2-5
Q Losses (I^2*X) - 6.38 MVAr @ line 2-5
================================================================================
| Bus Data |
================================================================================
Bus Voltage Generation Load
# Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr)
--- ------- -------- -------- -------- -------- --------
1 1.050 0.000 138.77 1.06 - -
2 1.034 -2.749 57.56 16.11 21.70 12.70
3 1.024 -4.602 - - 2.40 1.20
4 1.018 -5.510 - - 7.60 1.60
5 1.000 -8.990 24.56 22.25 94.20 19.00
6 1.011 -6.332 - - - -
7 0.999 -7.957 - - 22.80 10.90
8 1.012 -6.334 35.00 32.27 30.00 30.00
9 1.031 -7.881 - - - -
10 1.016 -9.788 - - 5.80 2.00
11 1.065 -5.934 17.93 17.61 - -
12 1.044 -9.336 - - 11.20 7.50
13 1.088 -8.142 16.91 34.74 - -
14 1.027 -10.237 - - 6.20 1.60
15 1.020 -10.272 - - 8.20 2.50
16 1.025 -9.800 - - 3.50 1.80
17 1.014 -10.006 - - 9.00 5.80
18 1.007 -10.826 - - 3.20 0.90
19 1.002 -10.957 - - 9.50 3.40
20 1.005 -10.723 - - 2.20 0.70
21 1.004 -10.270 - - 17.50 11.20
22 1.005 -10.260 - - - -
23 1.006 -10.620 - - 3.20 1.60
24 0.995 -10.727 - - 8.70 6.70
25 1.003 -10.774 - - - -
26 0.985 -11.205 - - 3.50 2.30
27 1.017 -10.522 - - - -
28 1.006 -6.745 - - - -
29 0.997 -11.767 - - 2.40 0.90
30 0.986 -12.661 - - 10.60 1.90
-------- -------- -------- --------
Total: 290.73 124.04 283.40 126.20
================================================================================
| Branch Data |
================================================================================
Brnch From To From Bus Injection To Bus Injection Loss (I^2 * Z)
# Bus Bus P (MW) Q (MVAr) P (MW) Q (MVAr) P (MW) Q (MVAr)
----- ----- ----- -------- -------- -------- -------- -------- --------
1 1 2 91.00 -1.54 -89.56 0.13 1.442 4.32
2 1 3 47.77 2.60 -46.83 -3.12 0.945 3.87
3 2 4 29.39 -1.33 -28.93 -1.13 0.461 1.40
4 3 4 44.43 1.92 -44.18 -2.08 0.249 0.72
5 2 5 58.28 4.43 -56.76 -2.38 1.519 6.38
6 2 6 37.76 0.18 -36.98 -1.73 0.778 2.36
7 4 6 37.46 5.92 -37.29 -6.27 0.166 0.58
8 5 7 -12.88 5.63 12.98 -7.42 0.097 0.24
9 6 7 36.12 2.82 -35.78 -3.48 0.345 1.06
10 6 8 -0.75 -3.50 0.75 2.58 0.001 0.00
11 6 9 13.76 -1.98 -13.76 2.36 0.000 0.38
12 6 10 11.57 6.64 -11.57 -5.74 0.000 0.90
13 9 11 -17.93 -16.45 17.93 17.61 0.000 1.16
14 9 10 31.69 14.09 -31.69 -12.85 0.000 1.24
15 4 12 28.04 -4.31 -28.04 6.25 0.000 1.94
16 12 13 -16.91 -32.98 16.91 34.74 0.000 1.76
17 12 14 8.09 3.12 -8.00 -2.95 0.085 0.18
18 12 15 18.26 9.71 -18.00 -9.20 0.260 0.51
19 12 16 7.40 6.40 -7.32 -6.22 0.083 0.17
20 14 15 1.80 1.35 -1.79 -1.34 0.011 0.01
21 16 17 3.82 4.42 -3.79 -4.36 0.027 0.06
22 15 18 6.21 3.40 -6.16 -3.29 0.052 0.11
23 18 19 2.96 2.39 -2.95 -2.37 0.009 0.02
24 19 20 -6.55 -1.03 6.56 1.06 0.015 0.03
25 10 20 8.84 1.92 -8.76 -1.76 0.074 0.17
26 10 17 5.21 1.46 -5.21 -1.44 0.009 0.02
27 10 21 15.79 9.28 -15.68 -9.04 0.113 0.24
28 10 22 7.62 4.12 -7.56 -4.01 0.053 0.11
29 21 22 -1.82 -2.16 1.82 2.16 0.001 0.00
30 15 23 5.38 4.64 -5.33 -4.54 0.048 0.10
31 22 24 5.74 1.85 -5.70 -1.78 0.041 0.06
32 23 24 2.13 2.94 -2.12 -2.90 0.017 0.04
33 24 25 -0.88 -1.98 0.89 1.99 0.009 0.02
34 25 26 3.55 2.37 -3.50 -2.30 0.046 0.07
35 25 27 -4.44 -4.36 4.48 4.44 0.042 0.08
36 28 27 17.77 9.22 -17.77 -7.78 0.000 1.44
37 27 29 6.19 1.67 -6.10 -1.51 0.087 0.17
38 27 30 7.09 1.67 -6.93 -1.36 0.164 0.31
39 29 30 3.70 0.61 -3.67 -0.54 0.034 0.06
40 8 28 4.25 -0.31 -4.23 -4.01 0.013 0.04
41 6 28 13.57 4.01 -13.53 -5.21 0.034 0.12
-------- --------
Total: 7.331 32.46Newton's method power flow converged in 4 iterations.
Converged in 0.02 seconds
================================================================================
| System Summary |
================================================================================
How many? How much? P (MW) Q (MVAr)
--------------------- ------------------- ------------- -----------------
Buses 30 Total Gen Capacity 890.5 -75.0 to 257.5
Generators 6 On-line Capacity 890.5 -75.0 to 257.5
Committed Gens 6 Generation (actual) 290.7 123.0
Loads 21 Load 283.4 126.2
Fixed 21 Fixed 283.4 126.2
Dispatchable 0 Dispatchable 0.0 of 0.0 0.0
Shunts 2 Shunt (inj) 0.0 0.2
Branches 41 Losses (I^2 * Z) 7.33 32.46
Transformers 4 Branch Charging (inj) - 34.4
Inter-ties 7 Total Inter-tie Flow 105.2 39.7
Areas 3
Minimum Maximum
------------------------- --------------------------------
Voltage Magnitude 0.985 p.u. @ bus 26 1.088 p.u. @ bus 13
Voltage Angle -12.66 deg @ bus 30 0.00 deg @ bus 1
P Losses (I^2*R) - 1.52 MW @ line 2-5
Q Losses (I^2*X) - 6.38 MVAr @ line 2-5
================================================================================
| Bus Data |
================================================================================
Bus Voltage Generation Load
# Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr)
----- ------- -------- -------- -------- -------- --------
1 1.050 0.000 138.77 0.00 - -
2 1.034 -2.749 57.56 16.11 21.70 12.70
3 1.024 -4.602 - - 2.40 1.20
4 1.018 -5.510 - - 7.60 1.60
5 1.000 -8.990 24.56 22.25 94.20 19.00
6 1.011 -6.332 - - - -
7 0.999 -7.957 - - 22.80 10.90
8 1.012 -6.334 35.00 32.27 30.00 30.00
9 1.031 -7.881 - - - -
10 1.016 -9.788 - - 5.80 2.00
11 1.065 -5.934 17.93 17.61 - -
12 1.044 -9.336 - - 11.20 7.50
13 1.088 -8.142 16.91 34.74 - -
14 1.027 -10.237 - - 6.20 1.60
15 1.020 -10.272 - - 8.20 2.50
16 1.025 -9.800 - - 3.50 1.80
17 1.014 -10.006 - - 9.00 5.80
18 1.007 -10.826 - - 3.20 0.90
19 1.002 -10.957 - - 9.50 3.40
20 1.005 -10.723 - - 2.20 0.70
21 1.004 -10.270 - - 17.50 11.20
22 1.005 -10.260 - - - -
23 1.006 -10.620 - - 3.20 1.60
24 0.995 -10.727 - - 8.70 6.70
25 1.003 -10.774 - - - -
26 0.985 -11.205 - - 3.50 2.30
27 1.017 -10.522 - - - -
28 1.006 -6.745 - - - -
29 0.997 -11.767 - - 2.40 0.90
30 0.986 -12.661 - - 10.60 1.90
-------- -------- -------- --------
Total: 290.73 122.98 283.40 126.20
================================================================================
| Branch Data |
================================================================================
Brnch From To From Bus Injection To Bus Injection Loss (I^2 * Z)
# Bus Bus P (MW) Q (MVAr) P (MW) Q (MVAr) P (MW) Q (MVAr)
----- ----- ----- -------- -------- -------- -------- -------- --------
1 1 2 91.00 -1.54 -89.56 0.13 1.442 4.32
2 1 3 47.77 2.60 -46.83 -3.12 0.945 3.87
3 2 4 29.39 -1.33 -28.93 -1.13 0.461 1.40
4 3 4 44.43 1.92 -44.18 -2.08 0.249 0.72
5 2 5 58.28 4.43 -56.76 -2.38 1.519 6.38
6 2 6 37.76 0.18 -36.98 -1.73 0.778 2.36
7 4 6 37.46 5.92 -37.29 -6.27 0.166 0.58
8 5 7 -12.88 5.63 12.98 -7.42 0.097 0.24
9 6 7 36.12 2.82 -35.78 -3.48 0.345 1.06
10 6 8 -0.75 -3.50 0.75 2.58 0.001 0.00
11 6 9 13.76 -1.98 -13.76 2.36 0.000 0.38
12 6 10 11.57 6.64 -11.57 -5.74 0.000 0.90
13 9 11 -17.93 -16.45 17.93 17.61 0.000 1.16
14 9 10 31.69 14.09 -31.69 -12.85 0.000 1.24
15 4 12 28.04 -4.31 -28.04 6.25 0.000 1.94
16 12 13 -16.91 -32.98 16.91 34.74 0.000 1.76
17 12 14 8.09 3.12 -8.00 -2.95 0.085 0.18
18 12 15 18.26 9.71 -18.00 -9.20 0.260 0.51
19 12 16 7.40 6.40 -7.32 -6.22 0.083 0.17
20 14 15 1.80 1.35 -1.79 -1.34 0.011 0.01
21 16 17 3.82 4.42 -3.79 -4.36 0.027 0.06
22 15 18 6.21 3.40 -6.16 -3.29 0.052 0.11
23 18 19 2.96 2.39 -2.95 -2.37 0.009 0.02
24 19 20 -6.55 -1.03 6.56 1.06 0.015 0.03
25 10 20 8.84 1.92 -8.76 -1.76 0.074 0.17
26 10 17 5.21 1.46 -5.21 -1.44 0.009 0.02
27 10 21 15.79 9.28 -15.68 -9.04 0.113 0.24
28 10 22 7.62 4.12 -7.56 -4.01 0.053 0.11
29 21 22 -1.82 -2.16 1.82 2.16 0.001 0.00
30 15 23 5.38 4.64 -5.33 -4.54 0.048 0.10
31 22 24 5.74 1.85 -5.70 -1.78 0.041 0.06
32 23 24 2.13 2.94 -2.12 -2.90 0.017 0.04
33 24 25 -0.88 -1.98 0.89 1.99 0.009 0.02
34 25 26 3.55 2.37 -3.50 -2.30 0.046 0.07
35 25 27 -4.44 -4.36 4.48 4.44 0.042 0.08
36 28 27 17.77 9.22 -17.77 -7.78 0.000 1.44
37 27 29 6.19 1.67 -6.10 -1.51 0.087 0.17
38 27 30 7.09 1.67 -6.93 -1.36 0.164 0.31
39 29 30 3.70 0.61 -3.67 -0.54 0.034 0.06
40 8 28 4.25 -0.31 -4.23 -4.01 0.013 0.04
41 6 28 13.57 4.01 -13.53 -5.21 0.034 0.12
-------- --------
Total: 7.331 32.46function [baseMVA, bus, gen, branch, areas, gencost] = case30
%CASE30 Power flow data for 30 bus, 6 generator case.
% Please see 'help caseformat' for details on the case file format.
%
% Based on data from ...
% Alsac, O. & Stott, B., "Optimal Load Flow with Steady State Security",
% IEEE Transactions on Power Apparatus and Systems, Vol. PAS 93, No. 3,
% 1974, pp. 745-751.
% ... with branch parameters rounded to nearest 0.01, shunt values divided
% by 100 and shunt on bus 10 moved to bus 5, load at bus 5 zeroed out.
% Generator locations, costs and limits and bus areas were taken from ...
% Ferrero, R.W., Shahidehpour, S.M., Ramesh, V.C., "Transaction analysis
% in deregulated power systems using game theory", IEEE Transactions on
% Power Systems, Vol. 12, No. 3, Aug 1997, pp. 1340-1347.
% Generator Q limits were derived from Alsac & Stott, using their Pmax
% capacities. V limits and line |S| limits taken from Alsac & Stott.
% MATPOWER
% $Id: case30.m,v 1.7 2004/09/21 02:41:49 ray Exp $
%%----- Power Flow Data -----%%
%% system MVA base
baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va
baseKV zone Vmax Vmin
bus = [
1 3 0 0 0 0 1 1 0
110 1 1.1 0.97;
2 2 21.7 12.7 0 0 1 1 0
110 1 1.1 0.97;
3 1 2.4 1.2 0 0 1 1 0
110 1 1.1 0.97;
4 1 7.6 1.6 0 0 1 1 0
110 1 1.1 0.97;
5 1 94.2 19 0 0 1 1 0
110 1 1.1 0.97;
6 1 0 0 0 0 1 1 0
110 1 1.1 0.97;
7 1 22.8 10.9 0 0 1 1 0
110 1 1.1 0.97;
8 1 30 30 0 0 1 1 0
110 1 1.1 0.97;
9 1 0 0 0 0 1 1 0
110 1 1.1 0.97;
10 1 5.8 2 0 0.19 3 1 0
110 1 1.1 0.97;
11 1 0 0 0 0 1 1 0
110 1 1.1 0.97;
12 1 11.2 7.5 0 0 2 1 0
110 1 1.1 0.97;
13 2 0 0 0 0 2 1 0
110 1 1.1 0.97;
14 1 6.2 1.6 0 0 2 1 0
110 1 1.1 0.97;
15 1 8.2 2.5 0 0 2 1 0
110 1 1.1 0.97;
16 1 3.5 1.8 0 0 2 1 0
110 1 1.1 0.97;
17 1 9 5.8 0 0 2 1 0
110 1 1.1 0.97;
18 1 3.2 0.9 0 0 2 1 0
110 1 1.1 0.97;
19 1 9.5 3.4 0 0 2 1 0
110 1 1.1 0.97;
20 1 2.2 0.7 0 0 2 1 0
110 1 1.1 0.97;
21 1 17.5 11.2 0 0 3 1 0
110 1 1.1 0.97;
22 2 0 0 0 0 3 1 0
110 1 1.1 0.97;
23 2 3.2 1.6 0 0 2 1 0
110 1 1.1 0.97;
24 1 8.7 6.7 0 0.04 3 1 0
110 1 1.1 0.97;
25 1 0 0 0 0 3 1 0
110 1 1.1 0.97;
26 1 3.5 2.3 0 0 3 1 0
110 1 1.1 0.97;
27 2 0 0 0 0 3 1 0
110 1 1.1 0.97;
28 1 0 0 0 0 1 1 0
110 1 1.1 0.97;
29 1 2.4 0.9 0 0 3 1 0
110 1 1.1 0.97;
30 1 10.6 1.9 0 0 3 1 0
110 1 1.1 0.97;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax
Pmin
gen = [
1 138.53 -2.58 0 0 1.05 100 1 238.53
0;
2 57.56 2.43 60 -20 1.0338 100 1 157.56
0;
5 24.56 22.25 62.5 -15 1.0058 100 1 124.56
0;
8 35 32.27 50 -15 1.023 100 1 135
0;
11 17.93 17.61 40 -10 1.0913 100 1 117.93
0;
13 16.91 24.96 45 -15 1.0883 100 1 116.91
0;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio
angle status
branch = [
1 2 0.0192 0.0575 0.0528 130 130 130 0
0 1;
1 3 0.0452 0.1852 0.0408 130 130 130 0
0 1;
2 4 0.057 0.1737 0.0368 65 65 65 0
0 1;
3 4 0.0132 0.0379 0.0084 130 130 130 0
0 1;
2 5 0.0472 0.1983 0.0418 130 130 130 0
0 1;
2 6 0.0581 0.1763 0.0374 65 65 65 0
0 1;
4 6 0.0119 0.0414 0.009 90 90 90 0
0 1;
5 7 0.0460 0.1160 0.0204 70 70 70 0
0 1;
6 7 0.0267 0.082 0.017 130 130 130 0
0 1;
6 8 0.012 0.042 0.009 32 32 32 0
0 1;
6 9 0 0.208 0 65 65 65 0.9847
0 1;
6 10 0 0.556 0 32 32 32 0.9629
0 1;
9 11 0 0.208 0 65 65 65 0
0 1;
9 10 0 0.11 0 65 65 65 0
0 1;
4 12 0 0.256 0 65 65 65 0.9873
0 1;
12 13 0 0.14 0 65 65 65 0
0 1;
12 14 0.1231 0.2559 0 32 32 32 0
0 1;
12 15 0.0662 0.1304 0 32 32 32 0
0 1;
12 16 0.0945 0.1987 0 32 32 32 0
0 1;
14 15 0.221 0.1997 0 16 16 16 0
0 1;
16 17 0.0824 0.1932 0 16 16 16 0
0 1;
15 18 0.107 0.2185 0 16 16 16 0
0 1;
18 19 0.0639 0.1292 0 16 16 16 0
0 1;
19 20 0.034 0.068 0 32 32 32 0
0 1;
10 20 0.0936 0.209 0 32 32 32 0
0 1;
10 17 0.0324 0.0845 0 32 32 32 0
0 1;
10 21 0.0348 0.0749 0 32 32 32 0
0 1;
10 22 0.0727 0.1499 0 32 32 32 0
0 1;
21 22 0.0116 0.0236 0 32 32 32 0
0 1;
15 23 0.1 0.202 0 16 16 16 0
0 1;
22 24 0.115 0.179 0 16 16 16 0
0 1;
23 24 0.132 0.27 0 16 16 16 0
0 1;
24 25 0.1885 0.3292 0 16 16 16 0
0 1;
25 26 0.2554 0.38 0 16 16 16 0
0 1;
25 27 0.1093 0.2087 0 16 16 16 0
0 1;
28 27 0 0.396 0 65 65 65 0.9581
0 1;
27 29 0.2198 0.4153 0 16 16 16 0
0 1;
27 30 0.3202 0.6027 0 16 16 16 0
0 1;
29 30 0.2399 0.4533 0 16 16 16 0
0 1;
8 28 0.0636 0.2 0.0428 32 32 32 0
0 1;
6 28 0.0169 0.0599 0.013 32 32 32 0
0 1;
];
%%----- OPF Data -----%%
%% area data
areas = [
1 8;
2 23;
3 26;
];
%% generator cost data
% 1 startup shutdown n x0 y0 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
gencost = [
2 0 0 3 0.02 2 0;
2 0 0 3 0.0175 1.75 0;
2 0 0 3 0.0625 1 0;
2 0 0 3 0.00834 3.25 0;
2 0 0 3 0.025 3 0;
2 0 0 3 0.025 3 0;
];
return;
