Thank you Ray. I tried to run the optimal power flow again by restarting the MATLAB and setting the "Working directory" again. Now, it works fine. Thank you so much.
On Tue, May 31, 2016 at 7:23 PM, Ray Zimmerman <[email protected]> wrote: > I am not sure why you are only seeing 2 generators. I was not able to > duplicate your result. I copied and pasted the data from your e-mail into a > new case file and get the result below (different solutions for the cases > with and without reactive costs) … > > Ray > > > >> mpopt = mpoption('opf.ac.solver', 'MIPS', 'out.suppress_detail', 1); > >> mpc = loadcase('case5Stagg') > > mpc = > > version: '2' > baseMVA: 100 > bus: [5x13 double] > gen: [3x21 double] > branch: [7x13 double] > gencost: [6x7 double] > > >> r = runopf(mpc, mpopt) > > MATPOWER Version 5.1, 20-Mar-2015 -- AC Optimal Power Flow > MATLAB Interior Point Solver -- MIPS, Version 1.2, 20-Mar-2015 > (using built-in linear solver) > Converged! > > Converged in 0.41 seconds > Objective Function Value = 2021.47 $/hr > > ================================================================================ > | System Summary > | > > ================================================================================ > > How many? How much? P (MW) Q (MVAr) > --------------------- > ------------------- ------------- ----------------- > Buses 5 Total Gen Capacity 250.0 -200.0 to > 250.0 > Generators 3 On-line Capacity 250.0 -200.0 to > 250.0 > Committed Gens 3 Generation (actual) 168.7 76.2 > Loads 4 Load 165.0 79.7 > Fixed 4 Fixed 165.0 79.7 > Dispatchable 0 Dispatchable -0.0 of -0.0 -0.0 > Shunts 0 Shunt (inj) -0.0 0.0 > Branches 7 Losses (I^2 * Z) 3.74 11.23 > Transformers 0 Branch Charging (inj) - 14.7 > Inter-ties 0 Total Inter-tie Flow 0.0 0.0 > Areas 1 > > Minimum Maximum > > ------------------------- -------------------------------- > Voltage Magnitude 0.991 p.u. @ bus 5 1.050 p.u. @ bus 1 > Voltage Angle -4.74 deg @ bus 5 0.00 deg @ bus 1 > P Losses (I^2*R) - 1.34 MW @ line 2-5 > Q Losses (I^2*X) - 4.01 MVAr @ line 2-5 > Lambda P 14.45 $/MWh @ bus 1 15.42 $/MWh @ bus 5 > Lambda Q 0.05 $/MWh @ bus 1 0.40 $/MWh @ bus 5 > > r = > > version: '2' > baseMVA: 100 > bus: [5x17 double] > gen: [3x25 double] > branch: [7x21 double] > gencost: [6x7 double] > order: [1x1 struct] > om: [1x1 opf_model] > x: [16x1 double] > mu: [1x1 struct] > f: 2.0215e+03 > var: [1x1 struct] > nln: [1x1 struct] > et: 0.4122 > success: 1 > raw: [1x1 struct] > > >> mpc1 = mpc; > >> mpc1.gencost = mpc1.gencost(1:3, :); > >> r = runopf(mpc1, mpopt) > > MATPOWER Version 5.1, 20-Mar-2015 -- AC Optimal Power Flow > MATLAB Interior Point Solver -- MIPS, Version 1.2, 20-Mar-2015 > (using built-in linear solver) > Converged! > > Converged in 0.08 seconds > Objective Function Value = 2012.97 $/hr > > ================================================================================ > | System Summary > | > > ================================================================================ > > How many? How much? P (MW) Q (MVAr) > --------------------- > ------------------- ------------- ----------------- > Buses 5 Total Gen Capacity 250.0 -200.0 to > 250.0 > Generators 3 On-line Capacity 250.0 -200.0 to > 250.0 > Committed Gens 3 Generation (actual) 168.7 75.8 > Loads 4 Load 165.0 79.7 > Fixed 4 Fixed 165.0 79.7 > Dispatchable 0 Dispatchable -0.0 of -0.0 -0.0 > Shunts 0 Shunt (inj) -0.0 0.0 > Branches 7 Losses (I^2 * Z) 3.69 11.07 > Transformers 0 Branch Charging (inj) - 14.9 > Inter-ties 0 Total Inter-tie Flow 0.0 0.0 > Areas 1 > > Minimum Maximum > > ------------------------- -------------------------------- > Voltage Magnitude 1.001 p.u. @ bus 5 1.050 p.u. @ bus 1 > Voltage Angle -4.84 deg @ bus 5 0.00 deg @ bus 1 > P Losses (I^2*R) - 1.32 MW @ line 2-5 > Q Losses (I^2*X) - 3.95 MVAr @ line 2-5 > Lambda P 14.45 $/MWh @ bus 1 15.38 $/MWh @ bus 5 > Lambda Q -0.00 $/MWh @ bus 1 0.25 $/MWh @ bus 5 > > r = > > version: '2' > baseMVA: 100 > bus: [5x17 double] > gen: [3x25 double] > branch: [7x21 double] > gencost: [3x7 double] > order: [1x1 struct] > om: [1x1 opf_model] > x: [16x1 double] > mu: [1x1 struct] > f: 2.0130e+03 > var: [1x1 struct] > nln: [1x1 struct] > et: 0.0820 > success: 1 > raw: [1x1 struct] > > > > > > On May 31, 2016, at 2:52 AM, Saranya A <[email protected]> wrote: > > Hi, > I am using MATPOWER 5.1. > > I am trying to run opf for a 5 bus system mentioned in a book by "Stagg > and El-Abiad". > > Following is my input system data: > > function mpc = case5Stagg > %CASE5Stagg This represents the 5 bus system in Stagg El-Abiad book > > %% MATPOWER Case Format : Version 2 > mpc.version = '2'; > > %%----- Power Flow Data -----%% > %% system MVA base > mpc.baseMVA = 100; > > %% bus data > % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin > mpc.bus = [ > 1 3 0 0 0 0 1 1.05 0 230 1 1.05 0.95; > 2 2 20 9.7 0 0 1 1 0 230 1 1.05 0.95; > 3 1 45 22 0 0 1 1 0 230 1 1.05 0.95; > 4 2 40 19 0 0 1 1 0 230 1 1.05 0.95; > 5 1 60 29 0 0 1 1 0 230 1 1.05 0.95; > ]; > > %% generator data > % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max > Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf > mpc.gen = [ > 1 100 20 100 -100 1 100 1 125 20 0 0 0 0 0 0 0 0 0 0 0; > 2 80 20 100 -100 1 100 1 125 20 0 0 0 0 0 0 0 0 0 0 0; > 4 0 50 50 0 1 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0; > ]; > > %% branch data > % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax > mpc.branch = [ > 1 2 0.02 0.06 0.03 180 180 180 0 0 1 -360 360; > 1 3 0.08 0.24 0.025 180 180 180 0 0 1 -360 360; > 2 3 0.06 0.18 0.02 180 180 180 0 0 1 -360 360; > 2 4 0.06 0.18 0.02 180 180 180 0 0 1 -360 360; > 2 5 0.04 0.12 0.01 180 180 180 0 0 1 -360 360; > 3 4 0.01 0.03 0.01 180 180 180 0 0 1 -360 360; > 4 5 0.08 0.24 0.025 180 180 180 0 0 1 -360 360; > ]; > > %%----- OPF Data -----%% > %% generator cost data > % 1 startup shutdown n x1 y1 ... xn yn > % 2 startup shutdown n c(n-1) ... c0 > mpc.gencost = [ > 2 0 0 3 0.0420 7.50 75; > 2 0 0 3 0.0420 7.50 75; > 2 0 0 3 0 0 0; > 2 0 0 3 0.0043 -0.0383 0.2992; > 2 0 0 3 0.0043 -0.0383 0.2992; > 2 0 0 3 0 0.1324 0; > ]; > > This system has 2 generators at bus 1 and 2 and a capacitor at bus 4. > > I am also including cost of generation of reactive power from the > generators as well as the capacitors(refer mpc.gencost). > > But after I 'runopf' , the opf results are the same with/without reactive > power generation cost. > > results= > > version: '2' > baseMVA: 100 > bus: [5x17 double] > gen: [2x25 double] > branch: [7x21 double] > gencost: [2x7 double] > order: [1x1 struct] > om: [1x1 opf_model] > x: [14x1 double] > mu: [1x1 struct] > f: 2.0196e+03 > var: [1x1 struct] > nln: [1x1 struct] > et: 0.0930 > success: 1 > raw: [1x1 struct] > > You can see that the gencost is only a 2x7 double. Why is the program not > taking the other rows into account? > > Please help me out. > Thank you. > > >
