0001 function t_opf_dc_mips(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 (IPOPT): ';
0062 mpopt = mpoption('out.all', 0, 'verbose', verbose);
0063 mpopt = mpoption(mpopt, 'opf.dc.solver', 'IPOPT');
0064
0065
0066 ib_data = [1:BUS_AREA BASE_KV:VMIN];
0067 ib_voltage = [VM VA];
0068 ib_lam = [LAM_P LAM_Q];
0069 ib_mu = [MU_VMAX MU_VMIN];
0070 ig_data = [GEN_BUS QMAX QMIN MBASE:APF];
0071 ig_disp = [PG QG VG];
0072 ig_mu = (MU_PMAX:MU_QMIN);
0073 ibr_data = (1:ANGMAX);
0074 ibr_flow = (PF:QT);
0075 ibr_mu = [MU_SF MU_ST];
0076 ibr_angmu = [MU_ANGMIN MU_ANGMAX];
0077
0078
0079 load soln9_dcopf;
0080
0081
0082 t = t0;
0083 [baseMVA, bus, gen, gencost, branch, f, success, et] = rundcopf(casefile, mpopt);
0084 t_ok(success, [t 'success']);
0085 t_is(f, f_soln, 3, [t 'f']);
0086 t_is( bus(:,ib_data ), bus_soln(:,ib_data ), 10, [t 'bus data']);
0087 t_is( bus(:,ib_voltage), bus_soln(:,ib_voltage), 3, [t 'bus voltage']);
0088 t_is( bus(:,ib_lam ), bus_soln(:,ib_lam ), 3, [t 'bus lambda']);
0089 t_is( bus(:,ib_mu ), bus_soln(:,ib_mu ), 2, [t 'bus mu']);
0090 t_is( gen(:,ig_data ), gen_soln(:,ig_data ), 10, [t 'gen data']);
0091 t_is( gen(:,ig_disp ), gen_soln(:,ig_disp ), 3, [t 'gen dispatch']);
0092 t_is( gen(:,ig_mu ), gen_soln(:,ig_mu ), 3, [t 'gen mu']);
0093 t_is(branch(:,ibr_data ), branch_soln(:,ibr_data ), 10, [t 'branch data']);
0094 t_is(branch(:,ibr_flow ), branch_soln(:,ibr_flow ), 3, [t 'branch flow']);
0095 t_is(branch(:,ibr_mu ), branch_soln(:,ibr_mu ), 2, [t 'branch mu']);
0096
0097
0098
0099
0100
0101
0102 mpc = loadcase(casefile);
0103 mpc.A = sparse([1;1;1;2;2;2],[10;11;13;11;12;14],[-1;1;-1;1;-1;-1],2,14);
0104 mpc.u = [0; 0];
0105 mpc.l = [-Inf; -Inf];
0106 mpc.zl = [0; 0];
0107
0108 mpc.N = sparse([1;2], [13;14], [1;1], 2, 14);
0109 mpc.fparm = ones(2,1) * [1 0 0 1];
0110 mpc.H = sparse(2,2);
0111 mpc.Cw = [1000;1];
0112
0113 t = [t0 'w/extra constraints & costs 1 : '];
0114 [r, success] = rundcopf(mpc, mpopt);
0115 t_ok(success, [t 'success']);
0116 t_is(r.gen(1, PG), 116.15974, 4, [t 'Pg1 = 116.15974']);
0117 t_is(r.gen(2, PG), 116.15974, 4, [t 'Pg2 = 116.15974']);
0118 t_is(r.var.val.z, [0; 0.3348], 4, [t 'user vars']);
0119 t_is(r.cost.usr, 0.3348, 3, [t 'user costs']);
0120
0121
0122 mpc = loadcase(casefile);
0123 mpc.A = sparse([1;1;1;2;2;2],[19;20;25;20;21;26],[-1;1;-1;1;-1;-1],2,26);
0124 mpc.u = [0; 0];
0125 mpc.l = [-Inf; -Inf];
0126 mpc.zl = [0; 0];
0127
0128 mpc.N = sparse([1;2], [25;26], [1;1], 2, 26);
0129 mpc.fparm = ones(2,1) * [1 0 0 1];
0130 mpc.H = sparse(2,2);
0131 mpc.Cw = [1000;1];
0132
0133 t = [t0 'w/extra constraints & costs 2 : '];
0134 [r, success] = rundcopf(mpc, mpopt);
0135 t_ok(success, [t 'success']);
0136 t_is(r.gen(1, PG), 116.15974, 4, [t 'Pg1 = 116.15974']);
0137 t_is(r.gen(2, PG), 116.15974, 4, [t 'Pg2 = 116.15974']);
0138 t_is(r.var.val.z, [0; 0.3348], 4, [t 'user vars']);
0139 t_is(r.cost.usr, 0.3348, 3, [t 'user costs']);
0140
0141 t = [t0 'infeasible : '];
0142
0143 mpc = loadcase(casefile);
0144 mpc.A = sparse([1;1], [10;11], [1;1], 1, 14);
0145 mpc.u = Inf;
0146 mpc.l = 600;
0147 [r, success] = rundcopf(mpc, mpopt);
0148 t_ok(~success, [t 'no success']);
0149
0150 if have_fcn('octave')
0151 warning(s1.state, 'Octave:load-file-in-path');
0152 end
0153
0154 t_end;