0001 function t_opf_dc_mips_sc(quiet)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033 if nargin < 1
0034 quiet = 0;
0035 end
0036
0037 num_tests = 23;
0038
0039 t_begin(num_tests, quiet);
0040
0041 [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
0042 VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
0043 [GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
0044 MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ...
0045 QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
0046 [F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, ...
0047 TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ...
0048 ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;
0049
0050 casefile = 't_case9_opf';
0051 if quiet
0052 verbose = 0;
0053 else
0054 verbose = 0;
0055 end
0056 if have_fcn('octave')
0057 s1 = warning('query', 'Octave:load-file-in-path');
0058 warning('off', 'Octave:load-file-in-path');
0059 end
0060
0061 t0 = 'DC OPF (MIPS-sc): ';
0062 mpopt = mpoption('out.all', 0, 'verbose', verbose);
0063 mpopt = mpoption(mpopt, 'opf.dc.solver', 'MIPS', 'mips.step_control', 1);
0064
0065 s2 = warning('query', 'MATLAB:singularMatrix');
0066
0067
0068 ib_data = [1:BUS_AREA BASE_KV:VMIN];
0069 ib_voltage = [VM VA];
0070 ib_lam = [LAM_P LAM_Q];
0071 ib_mu = [MU_VMAX MU_VMIN];
0072 ig_data = [GEN_BUS QMAX QMIN MBASE:APF];
0073 ig_disp = [PG QG VG];
0074 ig_mu = (MU_PMAX:MU_QMIN);
0075 ibr_data = (1:ANGMAX);
0076 ibr_flow = (PF:QT);
0077 ibr_mu = [MU_SF MU_ST];
0078 ibr_angmu = [MU_ANGMIN MU_ANGMAX];
0079
0080
0081 load soln9_dcopf;
0082
0083
0084 t = t0;
0085 [baseMVA, bus, gen, gencost, branch, f, success, et] = rundcopf(casefile, mpopt);
0086 t_ok(success, [t 'success']);
0087 t_is(f, f_soln, 3, [t 'f']);
0088 t_is( bus(:,ib_data ), bus_soln(:,ib_data ), 10, [t 'bus data']);
0089 t_is( bus(:,ib_voltage), bus_soln(:,ib_voltage), 3, [t 'bus voltage']);
0090 t_is( bus(:,ib_lam ), bus_soln(:,ib_lam ), 3, [t 'bus lambda']);
0091 t_is( bus(:,ib_mu ), bus_soln(:,ib_mu ), 2, [t 'bus mu']);
0092 t_is( gen(:,ig_data ), gen_soln(:,ig_data ), 10, [t 'gen data']);
0093 t_is( gen(:,ig_disp ), gen_soln(:,ig_disp ), 3, [t 'gen dispatch']);
0094 t_is( gen(:,ig_mu ), gen_soln(:,ig_mu ), 3, [t 'gen mu']);
0095 t_is(branch(:,ibr_data ), branch_soln(:,ibr_data ), 10, [t 'branch data']);
0096 t_is(branch(:,ibr_flow ), branch_soln(:,ibr_flow ), 3, [t 'branch flow']);
0097 t_is(branch(:,ibr_mu ), branch_soln(:,ibr_mu ), 2, [t 'branch mu']);
0098
0099
0100
0101
0102
0103
0104 mpc = loadcase(casefile);
0105 mpc.A = sparse([1;1;1;2;2;2],[10;11;13;11;12;14],[-1;1;-1;1;-1;-1],2,14);
0106 mpc.u = [0; 0];
0107 mpc.l = [-Inf; -Inf];
0108 mpc.zl = [0; 0];
0109
0110 mpc.N = sparse([1;2], [13;14], [1;1], 2, 14);
0111 mpc.fparm = ones(2,1) * [1 0 0 1];
0112 mpc.H = sparse(2,2);
0113 mpc.Cw = [1000;1];
0114
0115 t = [t0 'w/extra constraints & costs 1 : '];
0116 [r, success] = rundcopf(mpc, mpopt);
0117 t_ok(success, [t 'success']);
0118 t_is(r.gen(1, PG), 116.15974, 4, [t 'Pg1 = 116.15974']);
0119 t_is(r.gen(2, PG), 116.15974, 4, [t 'Pg2 = 116.15974']);
0120 t_is(r.var.val.z, [0; 0.3348], 4, [t 'user vars']);
0121 t_is(r.cost.usr, 0.3348, 3, [t 'user costs']);
0122
0123
0124 mpc = loadcase(casefile);
0125 mpc.A = sparse([1;1;1;2;2;2],[19;20;25;20;21;26],[-1;1;-1;1;-1;-1],2,26);
0126 mpc.u = [0; 0];
0127 mpc.l = [-Inf; -Inf];
0128 mpc.zl = [0; 0];
0129
0130 mpc.N = sparse([1;2], [25;26], [1;1], 2, 26);
0131 mpc.fparm = ones(2,1) * [1 0 0 1];
0132 mpc.H = sparse(2,2);
0133 mpc.Cw = [1000;1];
0134
0135 t = [t0 'w/extra constraints & costs 2 : '];
0136 [r, success] = rundcopf(mpc, mpopt);
0137 t_ok(success, [t 'success']);
0138 t_is(r.gen(1, PG), 116.15974, 4, [t 'Pg1 = 116.15974']);
0139 t_is(r.gen(2, PG), 116.15974, 4, [t 'Pg2 = 116.15974']);
0140 t_is(r.var.val.z, [0; 0.3348], 4, [t 'user vars']);
0141 t_is(r.cost.usr, 0.3348, 3, [t 'user costs']);
0142
0143 t = [t0 'infeasible : '];
0144 warning('off', 'MATLAB:singularMatrix');
0145
0146 mpc = loadcase(casefile);
0147 mpc.A = sparse([1;1], [10;11], [1;1], 1, 14);
0148 mpc.u = Inf;
0149 mpc.l = 600;
0150 [r, success] = rundcopf(mpc, mpopt);
0151 t_ok(~success, [t 'no success']);
0152
0153 if have_fcn('octave')
0154 warning(s1.state, 'Octave:load-file-in-path');
0155 end
0156 warning(s2.state, 'MATLAB:singularMatrix');
0157
0158 t_end;