0001 function t_totcost(quiet)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014 if nargin < 1
0015 quiet = 0;
0016 end
0017
0018 n_tests = 22;
0019
0020 t_begin(n_tests, quiet);
0021
0022
0023 [GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
0024 MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ...
0025 QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
0026 [PW_LINEAR, POLYNOMIAL, MODEL, STARTUP, SHUTDOWN, NCOST, COST] = idx_cost;
0027
0028
0029
0030
0031 gencost = [
0032 2 0 0 3 0.01 0.1 1 0 0 0 0 0;
0033 2 0 0 5 0.0006 0.005 0.04 0.3 2 0 0 0;
0034 1 0 0 4 0 0 10 200 20 600 30 1200;
0035 1 0 0 4 -30 -2400 -20 -1800 -10 -1000 0 0;
0036 ];
0037
0038 t = 'totcost - quadratic';
0039 t_is(totcost(gencost, [0;0;0;0]), [1;2;0;0], 8, t);
0040 t_is(totcost(gencost, [1;0;0;0]), [1.11;2;0;0], 8, t);
0041 t_is(totcost(gencost, [2;0;0;0]), [1.24;2;0;0], 8, t);
0042
0043 t = 'totcost - 4th order polynomial';
0044 t_is(totcost(gencost, [0;0;0;0]), [1;2; 0;0], 8, t);
0045 t_is(totcost(gencost, [0;1;0;0]), [1;2.3456;0;0], 8, t);
0046 t_is(totcost(gencost, [0;2;0;0]), [1;2.8096;0;0], 8, t);
0047
0048 t = 'totcost - pwl (gen)';
0049 t_is(totcost(gencost, [0;0;-10;0 ]), [1;2;-200;0], 8, t);
0050 t_is(totcost(gencost, [0;0;5;0 ]), [1;2;100;0], 8, t);
0051 t_is(totcost(gencost, [0;0;10;0]), [1;2;200;0], 8, t);
0052 t_is(totcost(gencost, [0;0;15;0]), [1;2;400;0], 8, t);
0053 t_is(totcost(gencost, [0;0;20;0]), [1;2;600;0], 8, t);
0054 t_is(totcost(gencost, [0;0;25;0]), [1;2;900;0], 8, t);
0055 t_is(totcost(gencost, [0;0;30;0]), [1;2;1200;0], 8, t);
0056 t_is(totcost(gencost, [0;0;35;0]), [1;2;1500;0], 8, t);
0057
0058 t = 'totcost - pwl (load)';
0059 t_is(totcost(gencost, [0;0;0;10 ]), [1;2;0;1000], 8, t);
0060 t_is(totcost(gencost, [0;0;0;-5 ]), [1;2;0;-500], 8, t);
0061 t_is(totcost(gencost, [0;0;0;-10]), [1;2;0;-1000], 8, t);
0062 t_is(totcost(gencost, [0;0;0;-15]), [1;2;0;-1400], 8, t);
0063 t_is(totcost(gencost, [0;0;0;-20]), [1;2;0;-1800], 8, t);
0064 t_is(totcost(gencost, [0;0;0;-25]), [1;2;0;-2100], 8, t);
0065 t_is(totcost(gencost, [0;0;0;-30]), [1;2;0;-2400], 8, t);
0066 t_is(totcost(gencost, [0;0;0;-35]), [1;2;0;-2700], 8, t);
0067
0068 t_end;