Sorry, I meant the GS column. Ray
On Mar 7, 2014, at 9:00 AM, Ray Zimmerman <[email protected]> wrote: > 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
