0001 function t_modcost(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 = 162;
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 gencost0 = [
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
0039 gencost = modcost(gencost0, 5, 'SCALE_F');
0040
0041 t = 'modcost SCALE_F - quadratic';
0042 t_is(totcost(gencost, [0;0;0;0])/5, [1;2;0;0], 8, t);
0043 t_is(totcost(gencost, [1;0;0;0])/5, [1.11;2;0;0], 8, t);
0044 t_is(totcost(gencost, [2;0;0;0])/5, [1.24;2;0;0], 8, t);
0045
0046 t = 'modcost SCALE_F - 4th order polynomial';
0047 t_is(totcost(gencost, [0;0;0;0])/5, [1;2; 0;0], 8, t);
0048 t_is(totcost(gencost, [0;1;0;0])/5, [1;2.3456;0;0], 8, t);
0049 t_is(totcost(gencost, [0;2;0;0])/5, [1;2.8096;0;0], 8, t);
0050
0051 t = 'modcost SCALE_F - pwl (gen)';
0052 t_is(totcost(gencost, [0;0;5;0 ])/5, [1;2;100;0], 8, t);
0053 t_is(totcost(gencost, [0;0;10;0])/5, [1;2;200;0], 8, t);
0054 t_is(totcost(gencost, [0;0;15;0])/5, [1;2;400;0], 8, t);
0055 t_is(totcost(gencost, [0;0;20;0])/5, [1;2;600;0], 8, t);
0056 t_is(totcost(gencost, [0;0;25;0])/5, [1;2;900;0], 8, t);
0057 t_is(totcost(gencost, [0;0;30;0])/5, [1;2;1200;0], 8, t);
0058 t_is(totcost(gencost, [0;0;35;0])/5, [1;2;1500;0], 8, t);
0059
0060 t = 'modcost SCALE_F - pwl (load)';
0061 t_is(totcost(gencost, [0;0;0;-5 ])/5, [1;2;0;-500], 8, t);
0062 t_is(totcost(gencost, [0;0;0;-10])/5, [1;2;0;-1000], 8, t);
0063 t_is(totcost(gencost, [0;0;0;-15])/5, [1;2;0;-1400], 8, t);
0064 t_is(totcost(gencost, [0;0;0;-20])/5, [1;2;0;-1800], 8, t);
0065 t_is(totcost(gencost, [0;0;0;-25])/5, [1;2;0;-2100], 8, t);
0066 t_is(totcost(gencost, [0;0;0;-30])/5, [1;2;0;-2400], 8, t);
0067 t_is(totcost(gencost, [0;0;0;-35])/5, [1;2;0;-2700], 8, t);
0068
0069
0070 gencost = modcost(gencost0, 2, 'SCALE_X');
0071
0072 t = 'modcost SCALE_X - quadratic';
0073 t_is(totcost(gencost, [0;0;0;0]*2), [1;2;0;0], 8, t);
0074 t_is(totcost(gencost, [1;0;0;0]*2), [1.11;2;0;0], 8, t);
0075 t_is(totcost(gencost, [2;0;0;0]*2), [1.24;2;0;0], 8, t);
0076
0077 t = 'modcost SCALE_X - 4th order polynomial';
0078 t_is(totcost(gencost, [0;0;0;0]*2), [1;2; 0;0], 8, t);
0079 t_is(totcost(gencost, [0;1;0;0]*2), [1;2.3456;0;0], 8, t);
0080 t_is(totcost(gencost, [0;2;0;0]*2), [1;2.8096;0;0], 8, t);
0081
0082 t = 'modcost SCALE_X - pwl (gen)';
0083 t_is(totcost(gencost, [0;0;5;0 ]*2), [1;2;100;0], 8, t);
0084 t_is(totcost(gencost, [0;0;10;0]*2), [1;2;200;0], 8, t);
0085 t_is(totcost(gencost, [0;0;15;0]*2), [1;2;400;0], 8, t);
0086 t_is(totcost(gencost, [0;0;20;0]*2), [1;2;600;0], 8, t);
0087 t_is(totcost(gencost, [0;0;25;0]*2), [1;2;900;0], 8, t);
0088 t_is(totcost(gencost, [0;0;30;0]*2), [1;2;1200;0], 8, t);
0089 t_is(totcost(gencost, [0;0;35;0]*2), [1;2;1500;0], 8, t);
0090
0091 t = 'modcost SCALE_X - pwl (load)';
0092 t_is(totcost(gencost, [0;0;0;-5 ]*2), [1;2;0;-500], 8, t);
0093 t_is(totcost(gencost, [0;0;0;-10]*2), [1;2;0;-1000], 8, t);
0094 t_is(totcost(gencost, [0;0;0;-15]*2), [1;2;0;-1400], 8, t);
0095 t_is(totcost(gencost, [0;0;0;-20]*2), [1;2;0;-1800], 8, t);
0096 t_is(totcost(gencost, [0;0;0;-25]*2), [1;2;0;-2100], 8, t);
0097 t_is(totcost(gencost, [0;0;0;-30]*2), [1;2;0;-2400], 8, t);
0098 t_is(totcost(gencost, [0;0;0;-35]*2), [1;2;0;-2700], 8, t);
0099
0100
0101 gencost = modcost(gencost0, 3, 'SHIFT_F');
0102
0103 t = 'modcost SHIFT_F - quadratic';
0104 t_is(totcost(gencost, [0;0;0;0])-3, [1;2;0;0], 8, t);
0105 t_is(totcost(gencost, [1;0;0;0])-3, [1.11;2;0;0], 8, t);
0106 t_is(totcost(gencost, [2;0;0;0])-3, [1.24;2;0;0], 8, t);
0107
0108 t = 'modcost SHIFT_F - 4th order polynomial';
0109 t_is(totcost(gencost, [0;0;0;0])-3, [1;2; 0;0], 8, t);
0110 t_is(totcost(gencost, [0;1;0;0])-3, [1;2.3456;0;0], 8, t);
0111 t_is(totcost(gencost, [0;2;0;0])-3, [1;2.8096;0;0], 8, t);
0112
0113 t = 'modcost SHIFT_F - pwl (gen)';
0114 t_is(totcost(gencost, [0;0;5;0 ])-3, [1;2;100;0], 8, t);
0115 t_is(totcost(gencost, [0;0;10;0])-3, [1;2;200;0], 8, t);
0116 t_is(totcost(gencost, [0;0;15;0])-3, [1;2;400;0], 8, t);
0117 t_is(totcost(gencost, [0;0;20;0])-3, [1;2;600;0], 8, t);
0118 t_is(totcost(gencost, [0;0;25;0])-3, [1;2;900;0], 8, t);
0119 t_is(totcost(gencost, [0;0;30;0])-3, [1;2;1200;0], 8, t);
0120 t_is(totcost(gencost, [0;0;35;0])-3, [1;2;1500;0], 8, t);
0121
0122 t = 'modcost SHIFT_F - pwl (load)';
0123 t_is(totcost(gencost, [0;0;0;-5 ])-3, [1;2;0;-500], 8, t);
0124 t_is(totcost(gencost, [0;0;0;-10])-3, [1;2;0;-1000], 8, t);
0125 t_is(totcost(gencost, [0;0;0;-15])-3, [1;2;0;-1400], 8, t);
0126 t_is(totcost(gencost, [0;0;0;-20])-3, [1;2;0;-1800], 8, t);
0127 t_is(totcost(gencost, [0;0;0;-25])-3, [1;2;0;-2100], 8, t);
0128 t_is(totcost(gencost, [0;0;0;-30])-3, [1;2;0;-2400], 8, t);
0129 t_is(totcost(gencost, [0;0;0;-35])-3, [1;2;0;-2700], 8, t);
0130
0131
0132 gencost = modcost(gencost0, -4, 'SHIFT_X');
0133
0134 t = 'modcost SHIFT_X - quadratic';
0135 t_is(totcost(gencost, [0;0;0;0]-4), [1;2;0;0], 8, t);
0136 t_is(totcost(gencost, [1;0;0;0]-4), [1.11;2;0;0], 8, t);
0137 t_is(totcost(gencost, [2;0;0;0]-4), [1.24;2;0;0], 8, t);
0138
0139 t = 'modcost SHIFT_X - 4th order polynomial';
0140 t_is(totcost(gencost, [0;0;0;0]-4), [1;2; 0;0], 8, t);
0141 t_is(totcost(gencost, [0;1;0;0]-4), [1;2.3456;0;0], 8, t);
0142 t_is(totcost(gencost, [0;2;0;0]-4), [1;2.8096;0;0], 8, t);
0143
0144 t = 'modcost SHIFT_X - pwl (gen)';
0145 t_is(totcost(gencost, [0;0;5;0 ]-4), [1;2;100;0], 8, t);
0146 t_is(totcost(gencost, [0;0;10;0]-4), [1;2;200;0], 8, t);
0147 t_is(totcost(gencost, [0;0;15;0]-4), [1;2;400;0], 8, t);
0148 t_is(totcost(gencost, [0;0;20;0]-4), [1;2;600;0], 8, t);
0149 t_is(totcost(gencost, [0;0;25;0]-4), [1;2;900;0], 8, t);
0150 t_is(totcost(gencost, [0;0;30;0]-4), [1;2;1200;0], 8, t);
0151 t_is(totcost(gencost, [0;0;35;0]-4), [1;2;1500;0], 8, t);
0152
0153 t = 'modcost SHIFT_X - pwl (load)';
0154 t_is(totcost(gencost, [0;0;0;-5 ]-4), [1;2;0;-500], 8, t);
0155 t_is(totcost(gencost, [0;0;0;-10]-4), [1;2;0;-1000], 8, t);
0156 t_is(totcost(gencost, [0;0;0;-15]-4), [1;2;0;-1400], 8, t);
0157 t_is(totcost(gencost, [0;0;0;-20]-4), [1;2;0;-1800], 8, t);
0158 t_is(totcost(gencost, [0;0;0;-25]-4), [1;2;0;-2100], 8, t);
0159 t_is(totcost(gencost, [0;0;0;-30]-4), [1;2;0;-2400], 8, t);
0160 t_is(totcost(gencost, [0;0;0;-35]-4), [1;2;0;-2700], 8, t);
0161
0162 t = 'modcost empty gencost';
0163 gencost = modcost([], 7);
0164 t_ok(isempty(gencost), t);
0165
0166
0167 alpha = [10; 9; 8; 7];
0168 gencost = modcost(gencost0, alpha, 'SCALE_F');
0169
0170 t = 'modcost vector SCALE_F - quadratic';
0171 t_is(totcost(gencost, [0;0;0;0])./alpha, [1;2;0;0], 8, t);
0172 t_is(totcost(gencost, [1;0;0;0])./alpha, [1.11;2;0;0], 8, t);
0173 t_is(totcost(gencost, [2;0;0;0])./alpha, [1.24;2;0;0], 8, t);
0174
0175 t = 'modcost vector SCALE_F - 4th order polynomial';
0176 t_is(totcost(gencost, [0;0;0;0])./alpha, [1;2; 0;0], 8, t);
0177 t_is(totcost(gencost, [0;1;0;0])./alpha, [1;2.3456;0;0], 8, t);
0178 t_is(totcost(gencost, [0;2;0;0])./alpha, [1;2.8096;0;0], 8, t);
0179
0180 t = 'modcost vector SCALE_F - pwl (gen)';
0181 t_is(totcost(gencost, [0;0;5;0 ])./alpha, [1;2;100;0], 8, t);
0182 t_is(totcost(gencost, [0;0;10;0])./alpha, [1;2;200;0], 8, t);
0183 t_is(totcost(gencost, [0;0;15;0])./alpha, [1;2;400;0], 8, t);
0184 t_is(totcost(gencost, [0;0;20;0])./alpha, [1;2;600;0], 8, t);
0185 t_is(totcost(gencost, [0;0;25;0])./alpha, [1;2;900;0], 8, t);
0186 t_is(totcost(gencost, [0;0;30;0])./alpha, [1;2;1200;0], 8, t);
0187 t_is(totcost(gencost, [0;0;35;0])./alpha, [1;2;1500;0], 8, t);
0188
0189 t = 'modcost vector SCALE_F - pwl (load)';
0190 t_is(totcost(gencost, [0;0;0;-5 ])./alpha, [1;2;0;-500], 8, t);
0191 t_is(totcost(gencost, [0;0;0;-10])./alpha, [1;2;0;-1000], 8, t);
0192 t_is(totcost(gencost, [0;0;0;-15])./alpha, [1;2;0;-1400], 8, t);
0193 t_is(totcost(gencost, [0;0;0;-20])./alpha, [1;2;0;-1800], 8, t);
0194 t_is(totcost(gencost, [0;0;0;-25])./alpha, [1;2;0;-2100], 8, t);
0195 t_is(totcost(gencost, [0;0;0;-30])./alpha, [1;2;0;-2400], 8, t);
0196 t_is(totcost(gencost, [0;0;0;-35])./alpha, [1;2;0;-2700], 8, t);
0197
0198
0199 gencost = modcost(gencost0, alpha, 'SCALE_X');
0200
0201 t = 'modcost vector SCALE_X - quadratic';
0202 t_is(totcost(gencost, [0;0;0;0].*alpha), [1;2;0;0], 8, t);
0203 t_is(totcost(gencost, [1;0;0;0].*alpha), [1.11;2;0;0], 8, t);
0204 t_is(totcost(gencost, [2;0;0;0].*alpha), [1.24;2;0;0], 8, t);
0205
0206 t = 'modcost vector SCALE_X - 4th order polynomial';
0207 t_is(totcost(gencost, [0;0;0;0].*alpha), [1;2; 0;0], 8, t);
0208 t_is(totcost(gencost, [0;1;0;0].*alpha), [1;2.3456;0;0], 8, t);
0209 t_is(totcost(gencost, [0;2;0;0].*alpha), [1;2.8096;0;0], 8, t);
0210
0211 t = 'modcost vector SCALE_X - pwl (gen)';
0212 t_is(totcost(gencost, [0;0;5;0 ].*alpha), [1;2;100;0], 8, t);
0213 t_is(totcost(gencost, [0;0;10;0].*alpha), [1;2;200;0], 8, t);
0214 t_is(totcost(gencost, [0;0;15;0].*alpha), [1;2;400;0], 8, t);
0215 t_is(totcost(gencost, [0;0;20;0].*alpha), [1;2;600;0], 8, t);
0216 t_is(totcost(gencost, [0;0;25;0].*alpha), [1;2;900;0], 8, t);
0217 t_is(totcost(gencost, [0;0;30;0].*alpha), [1;2;1200;0], 8, t);
0218 t_is(totcost(gencost, [0;0;35;0].*alpha), [1;2;1500;0], 8, t);
0219
0220 t = 'modcost vector SCALE_X - pwl (load)';
0221 t_is(totcost(gencost, [0;0;0;-5 ].*alpha), [1;2;0;-500], 8, t);
0222 t_is(totcost(gencost, [0;0;0;-10].*alpha), [1;2;0;-1000], 8, t);
0223 t_is(totcost(gencost, [0;0;0;-15].*alpha), [1;2;0;-1400], 8, t);
0224 t_is(totcost(gencost, [0;0;0;-20].*alpha), [1;2;0;-1800], 8, t);
0225 t_is(totcost(gencost, [0;0;0;-25].*alpha), [1;2;0;-2100], 8, t);
0226 t_is(totcost(gencost, [0;0;0;-30].*alpha), [1;2;0;-2400], 8, t);
0227 t_is(totcost(gencost, [0;0;0;-35].*alpha), [1;2;0;-2700], 8, t);
0228
0229
0230 gencost = modcost(gencost0, alpha, 'SHIFT_F');
0231
0232 t = 'modcost vector SHIFT_F - quadratic';
0233 t_is(totcost(gencost, [0;0;0;0])-alpha, [1;2;0;0], 8, t);
0234 t_is(totcost(gencost, [1;0;0;0])-alpha, [1.11;2;0;0], 8, t);
0235 t_is(totcost(gencost, [2;0;0;0])-alpha, [1.24;2;0;0], 8, t);
0236
0237 t = 'modcost vector SHIFT_F - 4th order polynomial';
0238 t_is(totcost(gencost, [0;0;0;0])-alpha, [1;2; 0;0], 8, t);
0239 t_is(totcost(gencost, [0;1;0;0])-alpha, [1;2.3456;0;0], 8, t);
0240 t_is(totcost(gencost, [0;2;0;0])-alpha, [1;2.8096;0;0], 8, t);
0241
0242 t = 'modcost vector SHIFT_F - pwl (gen)';
0243 t_is(totcost(gencost, [0;0;5;0 ])-alpha, [1;2;100;0], 8, t);
0244 t_is(totcost(gencost, [0;0;10;0])-alpha, [1;2;200;0], 8, t);
0245 t_is(totcost(gencost, [0;0;15;0])-alpha, [1;2;400;0], 8, t);
0246 t_is(totcost(gencost, [0;0;20;0])-alpha, [1;2;600;0], 8, t);
0247 t_is(totcost(gencost, [0;0;25;0])-alpha, [1;2;900;0], 8, t);
0248 t_is(totcost(gencost, [0;0;30;0])-alpha, [1;2;1200;0], 8, t);
0249 t_is(totcost(gencost, [0;0;35;0])-alpha, [1;2;1500;0], 8, t);
0250
0251 t = 'modcost vector SHIFT_F - pwl (load)';
0252 t_is(totcost(gencost, [0;0;0;-5 ])-alpha, [1;2;0;-500], 8, t);
0253 t_is(totcost(gencost, [0;0;0;-10])-alpha, [1;2;0;-1000], 8, t);
0254 t_is(totcost(gencost, [0;0;0;-15])-alpha, [1;2;0;-1400], 8, t);
0255 t_is(totcost(gencost, [0;0;0;-20])-alpha, [1;2;0;-1800], 8, t);
0256 t_is(totcost(gencost, [0;0;0;-25])-alpha, [1;2;0;-2100], 8, t);
0257 t_is(totcost(gencost, [0;0;0;-30])-alpha, [1;2;0;-2400], 8, t);
0258 t_is(totcost(gencost, [0;0;0;-35])-alpha, [1;2;0;-2700], 8, t);
0259
0260
0261 gencost = modcost(gencost0, -alpha, 'SHIFT_X');
0262
0263 t = 'modcost vector SHIFT_X - quadratic';
0264 t_is(totcost(gencost, [0;0;0;0]-alpha), [1;2;0;0], 8, t);
0265 t_is(totcost(gencost, [1;0;0;0]-alpha), [1.11;2;0;0], 8, t);
0266 t_is(totcost(gencost, [2;0;0;0]-alpha), [1.24;2;0;0], 8, t);
0267
0268 t = 'modcost vector SHIFT_X - 4th order polynomial';
0269 t_is(totcost(gencost, [0;0;0;0]-alpha), [1;2; 0;0], 8, t);
0270 t_is(totcost(gencost, [0;1;0;0]-alpha), [1;2.3456;0;0], 8, t);
0271 t_is(totcost(gencost, [0;2;0;0]-alpha), [1;2.8096;0;0], 8, t);
0272
0273 t = 'modcost vector SHIFT_X - pwl (gen)';
0274 t_is(totcost(gencost, [0;0;5;0 ]-alpha), [1;2;100;0], 8, t);
0275 t_is(totcost(gencost, [0;0;10;0]-alpha), [1;2;200;0], 8, t);
0276 t_is(totcost(gencost, [0;0;15;0]-alpha), [1;2;400;0], 8, t);
0277 t_is(totcost(gencost, [0;0;20;0]-alpha), [1;2;600;0], 8, t);
0278 t_is(totcost(gencost, [0;0;25;0]-alpha), [1;2;900;0], 8, t);
0279 t_is(totcost(gencost, [0;0;30;0]-alpha), [1;2;1200;0], 8, t);
0280 t_is(totcost(gencost, [0;0;35;0]-alpha), [1;2;1500;0], 8, t);
0281
0282 t = 'modcost vector SHIFT_X - pwl (load)';
0283 t_is(totcost(gencost, [0;0;0;-5 ]-alpha), [1;2;0;-500], 8, t);
0284 t_is(totcost(gencost, [0;0;0;-10]-alpha), [1;2;0;-1000], 8, t);
0285 t_is(totcost(gencost, [0;0;0;-15]-alpha), [1;2;0;-1400], 8, t);
0286 t_is(totcost(gencost, [0;0;0;-20]-alpha), [1;2;0;-1800], 8, t);
0287 t_is(totcost(gencost, [0;0;0;-25]-alpha), [1;2;0;-2100], 8, t);
0288 t_is(totcost(gencost, [0;0;0;-30]-alpha), [1;2;0;-2400], 8, t);
0289 t_is(totcost(gencost, [0;0;0;-35]-alpha), [1;2;0;-2700], 8, t);
0290
0291 t = 'modcost vector empty gencost';
0292 gencost = modcost([], alpha);
0293 t_ok(isempty(gencost), t);
0294
0295 t_end;