Based on what I see below, it looks like the difference is from shunts. Notice
there are some non-zeros in the BS column of the bus matrix.
Ray
>> mpc=loadcase('case300’);
>> mpopt=mpoption('OPF_ALG_DC',200,'OUT_BUS',0,'OUT_BRANCH',0,'OUT_ALL_LIM',0,'VERBOSE',2);
>> r=rundcopf(mpc,mpopt);
MATPOWER Version 4.1, 14-Dec-2011 -- DC Optimal Power Flow
MATLAB Interior Point Solver -- MIPS, Version 1.0, 07-Feb-2011
it objective step size feascond gradcond compcond costcond
---- ------------ --------- ------------ ------------ ------------ ------------
0 474427.12 0.119951 7750 6288.56 0
1 706292.47 23.586 5.03366e-16 0.914459 1900.17 0.488726
2 706811.13 2.8981 9.94131e-16 8.36768e-13 245.27 0.000734348
3 706424.94 1.4629 7.36521e-16 0.00385424 86.8407 0.000546389
4 706332.29 0.67175 5.0377e-16 1.13509e-12 12.8299 0.000131144
5 706301.52 0.42584 5.94015e-16 7.43969e-13 2.32315 4.35643e-05
6 706294.23 0.21314 3.77564e-16 6.39677e-13 0.49385 1.03216e-05
7 706292.64 0.10513 3.77521e-16 9.14906e-13 0.113122 2.26163e-06
8 706292.34 0.049447 3.14584e-16 5.36245e-13 0.0255032 4.16802e-07
9 706292.31 0.019494 5.03323e-16 3.28603e-13 0.00481366 5.00004e-08
10 706292.3 0.0044363 7.92828e-16 2.24905e-13 0.000608401 2.12219e-09
11 706292.3 0.00028709 6.2915e-16 5.39439e-13 6.14459e-05 9.09311e-12
12 706292.3 2.9243e-06 3.7749e-16 8.97475e-13 6.14465e-06 1.31861e-15
13 706292.3 1.7543e-07 3.08677e-16 7.94211e-13 6.14465e-07 1.64826e-16
Converged!
Converged in 0.88 seconds
Objective Function Value = 706292.30 $/hr
================================================================================
| System Summary |
================================================================================
How many? How much? P (MW) Q (MVAr)
--------------------- ------------------- ------------- -----------------
Buses 300 Total Gen Capacity 32678.4 0.0 to 0.0
Generators 69 On-line Capacity 32678.4 0.0 to 0.0
Committed Gens 69 Generation (actual) 23527.1 0.0
Loads 199 Load 23525.8 0.0
Fixed 199 Fixed 23525.8 0.0
Dispatchable 0 Dispatchable -0.0 of -0.0 -0.0
Shunts 17 Shunt (inj) -1.2 0.0
Branches 411 Losses (I^2 * Z) 0.00 0.00
Transformers 107 Branch Charging (inj) - 0.0
Inter-ties 0 Total Inter-tie Flow 0.0 0.0
Areas 1
Minimum Maximum
------------------------- --------------------------------
Voltage Magnitude 0.929 p.u. @ bus 9033 1.073 p.u. @ bus 149
Voltage Angle -17.48 deg @ bus 528 59.16 deg @ bus 7166
Lambda P 40.03 $/MWh @ bus 7166 40.03 $/MWh @ bus 7049
Lambda Q 0.00 $/MWh @ bus 1 0.00 $/MWh @ bus 1
On Mar 6, 2014, at 11:31 AM, Victor Hugo Hinojosa M. <[email protected]>
wrote:
> Dear Dr. Zimmerman,
> I run the DCOPF algorithm to solve the IEEE 300-bus system considering the
> default technical data for the system. Using the MIPS algorithm, the solution
> converges in 13 iterations with an objective function of 706292.30 $/hr. When
> I analyze the solution, the sum of the power generation (23527.15 MW) is
> bigger than the load power (23525.85 MW). The commands used are the following:
> mpc=case300; mpopt=mpoption; mpopt=mpoption('OPF_ALG',200,'OUT_ALL',0);
> salida_matpower=dcopf(mpc,mpopt);
> I’d like to know why there’s a difference of 1.3 MW
> (sum(salida_matpower.gen(:,2))-sum(mpc.bus(:,3))
> )
> Regards,
> Vh
