Hello *Kazi,
* Yah you are right. This is happening due to the xx matrix is empty,
because of the Pmin and Pmax of 15th generator is zero (Synch condenser).
While making cost function in 5 blocks, this takes empty matrix. This can be
temporarily avoided by putting very small value for Pmax of 15 th generator
in the "mpc.gen" matrix.
Dr Ray Zimmerman may consider this error, and add small checking for
synchronous condenser, with if condition in *case2off at 67 line *to solve
this issue permanently in the future version. Bye.
On Wed, Sep 29, 2010 at 6:35 PM, Kazi Hasan <[email protected]> wrote:
> Hello,
>
>
>
> runmarket.m is considering ‘case9’ as default data file.
>
> All other cases e.g. case14, case30, case39, case57, case118 and case300
> are working fine in runmarket.
>
>
>
> BUT WITH THE ‘case24_ieee_rts’ IT SHOWS THE FOLLOWING ERROR . . . .
>
>
>
> *??? Subscripted assignment dimension mismatch.*
>
> * *
>
> *Error in ==> case2off at 67*
>
> * gencost(ig, COST:2:(COST + 2*(npts-1) )) = xx;*
>
> * *
>
> *Error in ==> runmarket at 124*
>
> * [q, p] = case2off(mpc.gen, mpc.gencost);*
>
>
>
> Can anyone please advise me, what change I need to make in
> ‘case24_ieee_rts’ to run with runmarket.m.
>
>
>
> Thank you.
>
>
>
> Kazi
>
>
>
> *Kazi Nazmul Hasan*
>
> PhD Candidate
>
> Power and Energy Systems Research Group
>
> (Queensland Geothermal Energy Centre of Excellence)
>
> School of Information Technology & Electrical Engineering
>
> *The University of Queensland*
>
> Brisbane, QLD 4072
>
> *Australia*
>
> Phone: +61 7 3365 1654
>
>
>
>
>
>
>
--
-------------------------------------------------------------
WITH REGARDS,
V.RAVIKUMAR PANDI,
-------------------------------------------------------------
HAVE A NICE DAY
-------------------------------------------------------------
function mpc = case24_ieee_rts
%CASE24_IEEE_RTS Power flow data for the IEEE RELIABILITY TEST SYSTEM.
% Please see CASEFORMAT for details on the case file format.
%
% This system data is from the IEEE RELIABILITY TEST SYSTEM, see
%
% IEEE Reliability Test System Task Force of the Applications of
% Probability Methods Subcommittee, "IEEE reliability test system,"
% IEEE Transactions on Power Apparatus and Systems, Vol. 98, No. 6,
% Nov./Dec. 1979, pp. 2047-2054.
%
% IEEE Reliability Test System Task Force of Applications of
% Probability Methods Subcommittee, "IEEE reliability test system-96,"
% IEEE Transactions on Power Systems, Vol. 14, No. 3, Aug. 1999,
% pp. 1010-1020.
%
% Cost data is from Web site run by Georgia Tech Power Systems Control
% and Automation Laboratory:
%
% http://pscal.ece.gatech.edu/testsys/index.html
%
% MATPOWER case file data provided by Bruce Wollenberg.
% MATPOWER
% $Id: case24_ieee_rts.m,v 1.4 2010/03/10 18:08:13 ray Exp $
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va
baseKV zone Vmax Vmin
mpc.bus = [
1 2 108 22 0 0 1 1 0
138 1 1.05 0.95;
2 2 97 20 0 0 1 1 0
138 1 1.05 0.95;
3 1 180 37 0 0 1 1 0
138 1 1.05 0.95;
4 1 74 15 0 0 1 1 0
138 1 1.05 0.95;
5 1 71 14 0 0 1 1 0
138 1 1.05 0.95;
6 1 136 28 0 -100 2 1 0
138 1 1.05 0.95;
7 2 125 25 0 0 2 1 0
138 1 1.05 0.95;
8 1 171 35 0 0 2 1 0
138 1 1.05 0.95;
9 1 175 36 0 0 1 1 0
138 1 1.05 0.95;
10 1 195 40 0 0 2 1 0
138 1 1.05 0.95;
11 1 0 0 0 0 3 1 0
230 1 1.05 0.95;
12 1 0 0 0 0 3 1 0
230 1 1.05 0.95;
13 3 265 54 0 0 3 1 0
230 1 1.05 0.95;
14 2 194 39 0 0 3 1 0
230 1 1.05 0.95;
15 2 317 64 0 0 4 1 0
230 1 1.05 0.95;
16 2 100 20 0 0 4 1 0
230 1 1.05 0.95;
17 1 0 0 0 0 4 1 0
230 1 1.05 0.95;
18 2 333 68 0 0 4 1 0
230 1 1.05 0.95;
19 1 181 37 0 0 3 1 0
230 1 1.05 0.95;
20 1 128 26 0 0 3 1 0
230 1 1.05 0.95;
21 2 0 0 0 0 4 1 0
230 1 1.05 0.95;
22 2 0 0 0 0 4 1 0
230 1 1.05 0.95;
23 2 0 0 0 0 3 1 0
230 1 1.05 0.95;
24 1 0 0 0 0 4 1 0
230 1 1.05 0.95;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax
Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10
ramp_30 ramp_q apf % Unit Code
mpc.gen = [
1 10 0 10 0 1.035 100 1 20
16 0 0 0 0 0 0 0 0 0
0 0; % U20
1 10 0 10 0 1.035 100 1 20
16 0 0 0 0 0 0 0 0 0
0 0; % U20
1 76 0 30 -25 1.035 100 1 76
15.2 0 0 0 0 0 0 0 0 0
0 0; % U76
1 76 0 30 -25 1.035 100 1 76
15.2 0 0 0 0 0 0 0 0 0
0 0; % U76
2 10 0 10 0 1.035 100 1 20
16 0 0 0 0 0 0 0 0 0
0 0; % U20
2 10 0 10 0 1.035 100 1 20
16 0 0 0 0 0 0 0 0 0
0 0; % U20
2 76 0 30 -25 1.035 100 1 76
15.2 0 0 0 0 0 0 0 0 0
0 0; % U76
2 76 0 30 -25 1.035 100 1 76
15.2 0 0 0 0 0 0 0 0 0
0 0; % U76
7 80 0 60 0 1.025 100 1 100
25 0 0 0 0 0 0 0 0 0
0 0; % U100
7 80 0 60 0 1.025 100 1 100
25 0 0 0 0 0 0 0 0 0
0 0; % U100
7 80 0 60 0 1.025 100 1 100
25 0 0 0 0 0 0 0 0 0
0 0; % U100
13 95.1 0 80 0 1.02 100 1 197
69 0 0 0 0 0 0 0 0 0
0 0; % U197
13 95.1 0 80 0 1.02 100 1 197
69 0 0 0 0 0 0 0 0 0
0 0; % U197
13 95.1 0 80 0 1.02 100 1 197
69 0 0 0 0 0 0 0 0 0
0 0; % U197
14 0 35.3 200 -50 0.98 100 1 1e-10
0 0 0 0 0 0 0 0 0 0
0 0; % SynCond
15 12 0 6 0 1.014 100 1 12
2.4 0 0 0 0 0 0 0 0 0
0 0; % U12
15 12 0 6 0 1.014 100 1 12
2.4 0 0 0 0 0 0 0 0 0
0 0; % U12
15 12 0 6 0 1.014 100 1 12
2.4 0 0 0 0 0 0 0 0 0
0 0; % U12
15 12 0 6 0 1.014 100 1 12
2.4 0 0 0 0 0 0 0 0 0
0 0; % U12
15 12 0 6 0 1.014 100 1 12
2.4 0 0 0 0 0 0 0 0 0
0 0; % U12
15 155 0 80 -50 1.014 100 1 155
54.3 0 0 0 0 0 0 0 0 0
0 0; % U155
16 155 0 80 -50 1.017 100 1 155
54.3 0 0 0 0 0 0 0 0 0
0 0; % U155
18 400 0 200 -50 1.05 100 1 400
100 0 0 0 0 0 0 0 0 0
0 0; % U400
21 400 0 200 -50 1.05 100 1 400
100 0 0 0 0 0 0 0 0 0
0 0; % U400
22 50 0 16 -10 1.05 100 1 50
10 0 0 0 0 0 0 0 0 0
0 0; % U50
22 50 0 16 -10 1.05 100 1 50
10 0 0 0 0 0 0 0 0 0
0 0; % U50
22 50 0 16 -10 1.05 100 1 50
10 0 0 0 0 0 0 0 0 0
0 0; % U50
22 50 0 16 -10 1.05 100 1 50
10 0 0 0 0 0 0 0 0 0
0 0; % U50
22 50 0 16 -10 1.05 100 1 50
10 0 0 0 0 0 0 0 0 0
0 0; % U50
22 50 0 16 -10 1.05 100 1 50
10 0 0 0 0 0 0 0 0 0
0 0; % U50
23 155 0 80 -50 1.05 100 1 155
54.3 0 0 0 0 0 0 0 0 0
0 0; % U155
23 155 0 80 -50 1.05 100 1 155
54.3 0 0 0 0 0 0 0 0 0
0 0; % U155
23 350 0 150 -25 1.05 100 1 350
140 0 0 0 0 0 0 0 0 0
0 0; % U350
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio
angle status angmin angmax
mpc.branch = [
1 2 0.0026 0.0139 0.4611 175 250 200 0
0 1 -360 360;
1 3 0.0546 0.2112 0.0572 175 208 220 0
0 1 -360 360;
1 5 0.0218 0.0845 0.0229 175 208 220 0
0 1 -360 360;
2 4 0.0328 0.1267 0.0343 175 208 220 0
0 1 -360 360;
2 6 0.0497 0.192 0.052 175 208 220 0
0 1 -360 360;
3 9 0.0308 0.119 0.0322 175 208 220 0
0 1 -360 360;
3 24 0.0023 0.0839 0 400 510 600 1.03
0 1 -360 360;
4 9 0.0268 0.1037 0.0281 175 208 220 0
0 1 -360 360;
5 10 0.0228 0.0883 0.0239 175 208 220 0
0 1 -360 360;
6 10 0.0139 0.0605 2.459 175 193 200 0
0 1 -360 360;
7 8 0.0159 0.0614 0.0166 175 208 220 0
0 1 -360 360;
8 9 0.0427 0.1651 0.0447 175 208 220 0
0 1 -360 360;
8 10 0.0427 0.1651 0.0447 175 208 220 0
0 1 -360 360;
9 11 0.0023 0.0839 0 400 510 600 1.03
0 1 -360 360;
9 12 0.0023 0.0839 0 400 510 600 1.03
0 1 -360 360;
10 11 0.0023 0.0839 0 400 510 600 1.02
0 1 -360 360;
10 12 0.0023 0.0839 0 400 510 600 1.02
0 1 -360 360;
11 13 0.0061 0.0476 0.0999 500 600 625 0
0 1 -360 360;
11 14 0.0054 0.0418 0.0879 500 625 625 0
0 1 -360 360;
12 13 0.0061 0.0476 0.0999 500 625 625 0
0 1 -360 360;
12 23 0.0124 0.0966 0.203 500 625 625 0
0 1 -360 360;
13 23 0.0111 0.0865 0.1818 500 625 625 0
0 1 -360 360;
14 16 0.005 0.0389 0.0818 500 625 625 0
0 1 -360 360;
15 16 0.0022 0.0173 0.0364 500 600 625 0
0 1 -360 360;
15 21 0.0063 0.049 0.103 500 600 625 0
0 1 -360 360;
15 21 0.0063 0.049 0.103 500 600 625 0
0 1 -360 360;
15 24 0.0067 0.0519 0.1091 500 600 625 0
0 1 -360 360;
16 17 0.0033 0.0259 0.0545 500 600 625 0
0 1 -360 360;
16 19 0.003 0.0231 0.0485 500 600 625 0
0 1 -360 360;
17 18 0.0018 0.0144 0.0303 500 600 625 0
0 1 -360 360;
17 22 0.0135 0.1053 0.2212 500 600 625 0
0 1 -360 360;
18 21 0.0033 0.0259 0.0545 500 600 625 0
0 1 -360 360;
18 21 0.0033 0.0259 0.0545 500 600 625 0
0 1 -360 360;
19 20 0.0051 0.0396 0.0833 500 600 625 0
0 1 -360 360;
19 20 0.0051 0.0396 0.0833 500 600 625 0
0 1 -360 360;
20 23 0.0028 0.0216 0.0455 500 600 625 0
0 1 -360 360;
20 23 0.0028 0.0216 0.0455 500 600 625 0
0 1 -360 360;
21 22 0.0087 0.0678 0.1424 500 600 625 0
0 1 -360 360;
];
%%----- OPF Data -----%%
%% area data
% area refbus
mpc.areas = [
1 1;
2 3;
3 8;
4 6;
];
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [ %
bus Pmin Pmax Qmin Qmax Unit Code
2 1500 0 3 0 130 400.6849; %
1 16 20 0 10 U20
2 1500 0 3 0 130 400.6849; %
1 16 20 0 10 U20
2 1500 0 3 0.014142 16.0811 212.3076;
% 1 15.2 76 -25 30 U76
2 1500 0 3 0.014142 16.0811 212.3076;
% 1 15.2 76 -25 30 U76
2 1500 0 3 0 130 400.6849; %
2 16 20 0 10 U20
2 1500 0 3 0 130 400.6849; %
2 16 20 0 10 U20
2 1500 0 3 0.014142 16.0811 212.3076;
% 2 15.2 76 -25 30 U76
2 1500 0 3 0.014142 16.0811 212.3076;
% 2 15.2 76 -25 30 U76
2 1500 0 3 0.052672 43.6615 781.521;
% 7 25 100 0 60 U100
2 1500 0 3 0.052672 43.6615 781.521;
% 7 25 100 0 60 U100
2 1500 0 3 0.052672 43.6615 781.521;
% 7 25 100 0 60 U100
2 1500 0 3 0.00717 48.5804 832.7575; %
13 69 197 0 80 U197
2 1500 0 3 0.00717 48.5804 832.7575; %
13 69 197 0 80 U197
2 1500 0 3 0.00717 48.5804 832.7575; %
13 69 197 0 80 U197
2 1500 0 3 0 0 0; % 14
SynCond
2 1500 0 3 0.328412 56.564 86.3852;
% 15 2.4 12 0 6 U12
2 1500 0 3 0.328412 56.564 86.3852;
% 15 2.4 12 0 6 U12
2 1500 0 3 0.328412 56.564 86.3852;
% 15 2.4 12 0 6 U12
2 1500 0 3 0.328412 56.564 86.3852;
% 15 2.4 12 0 6 U12
2 1500 0 3 0.328412 56.564 86.3852;
% 15 2.4 12 0 6 U12
2 1500 0 3 0.008342 12.3883 382.2391;
% 15 54.3 155 -50 80 U155
2 1500 0 3 0.008342 12.3883 382.2391;
% 16 54.3 155 -50 80 U155
2 1500 0 3 0.000213 4.4231 395.3749;
% 18 100 400 -50 200 U400
2 1500 0 3 0.000213 4.4231 395.3749;
% 21 100 400 -50 200 U400
2 1500 0 3 0 0.001 0.001; % 22
10 50 -10 16 U50
2 1500 0 3 0 0.001 0.001; % 22
10 50 -10 16 U50
2 1500 0 3 0 0.001 0.001; % 22
10 50 -10 16 U50
2 1500 0 3 0 0.001 0.001; % 22
10 50 -10 16 U50
2 1500 0 3 0 0.001 0.001; % 22
10 50 -10 16 U50
2 1500 0 3 0 0.001 0.001; % 22
10 50 -10 16 U50
2 1500 0 3 0.008342 12.3883 382.2391;
% 23 54.3 155 -50 80 U155
2 1500 0 3 0.008342 12.3883 382.2391;
% 23 54.3 155 -50 80 U155
2 1500 0 3 0.004895 11.8495 665.1094;
% 23 140 350 -25 150 U350
];
