That is a MATLAB bug in their intlinprog() function. I would recommend using Gurobi or CPLEX for running MOST. I prefer Gurobi <https://user.gurobi.com/download/licenses/free-academic>. It is free for academic use.
Ray > On Mar 18, 2019, at 2:20 PM, Carlos Ferrandon Cervantes > <[email protected]> wrote: > > Thank you Doctor for your responses. I've been testing the software with > different values of reserve. I am using the 24 RTS bus. There is a point when > I am increasing the minimum reserve in the system when I get this error: > > Error using intlinprog (line 142) > INTLINPROG encountered an internal error that caused the solution to lose > integer and/or linear constraint > feasibility. We are sorry for the inconvenience. > > Please contact technical support for assistance with your problem, quoting > the code "-2_4". > > Error in miqps_ot (line 281) > intlinprog(c, intcon, Ai, bi, Ae, be, xmin, xmax, ot_opt); > > Error in miqps_matpower (line 220) > miqps_ot(H, c, A, l, u, xmin, xmax, x0, vtype, opt); > > Error in most (line 2154) > mdo.QP.lambda ] = miqps_matpower( mdi.QP.H, mdi.QP.C, ... > > Error in most_ex6_uc_DCFULL (line 117) > mdo = most(mdi, mpopt); > > I've noticed that the lower bound compared to the Optimal objective value is > somewhat far from each other. I dont know what this error means. > > Thank you as usual in advance, > > Carlos > > El mié., 13 mar. 2019 a las 19:42, Ray Zimmerman (<[email protected] > <mailto:[email protected]>>) escribió: > 1. For a deterministic case, you want to use the default values for > contingency constraint parameters in the xGen data, since they are not > relevant when there is a single state per period. Technically, in the > implementation they are independent, so it is possible to include both, but > the contingency constraints make no sense and have no impact on the > deterministic case. See the footnote on p. 35 of the MOST User’s Manual > <http://www.pserc.cornell.edu/matpower/docs/MOST-manual-1.0.1.pdf>. > > 2. Yes. See equations (7.2)–(7.5) in the MATPOWER User’s Manual > <http://www.pserc.cornell.edu/matpower/docs/MATPOWER-manual-7.0b1.pdf> for > the formulation. > > Ray > > > >> On Mar 13, 2019, at 12:37 PM, Carlos Ferrandon Cervantes >> <[email protected] <mailto:[email protected]>> wrote: >> >> Thank you Doctor for your response. I coded it that way just before calling >> most(mdi, mpopt) and it worked well and I've compared the results with fixed >> reserves and without. Now I have the following questions: >> >> 1. About the reserves, for a UC optimisation I had been using loadxgendata(' >> ',mpc) but now after adding mdi.FixedReserves(t,1,1)= mpc.reserves I >> wonder which one has the priority in the problem. Does one deactivate the >> other? Or what is the action here from the program. >> >> 2. Also, does the zonal value of reserve mpc.reserves.req is the minimum >> reserve that all generators allocated to certain zone in mpc.reserves.zones >> have to meet for the system? >> >> Thank you so much in advance, >> >> Carlos Ferrandon >> >> >> El lun., 11 mar. 2019 a las 15:48, Ray Zimmerman (<[email protected] >> <mailto:[email protected]>>) escribió: >> Hi Carlos, >> >> Sorry the documentation for this isn’t better. For fixed zonal reserves, >> MOST uses the basic structure used by MATPOWER and described in section >> 7.6.1 of the MATPOWER User’s Manual >> <http://www.pserc.cornell.edu/matpower/docs/MATPOWER-manual-7.0b1.pdf>. You >> will need to essentially specify the mpc.reserves structure for each period >> for MOST. After using loadmd() to create a MOST input data structure mdi, >> simply assign the corresponding to mpc.reserves to mdi.FixedReserves(t, 1, >> 1) for each period t. >> >> Hope this helps, >> >> Ray >> >> >>> On Mar 10, 2019, at 11:57 AM, Carlos Ferrandon Cervantes >>> <[email protected] <mailto:[email protected]>> wrote: >>> >>> Hello Doctor: >>> >>> I'm running a unit commitment (UC) simulation and I'm looking to set a >>> minimum variable value of MW of reserve for the whole system per each time >>> step (t) of the simulation. Assuming it is a 24hr simulation there could be >>> 24 different values of minimum MW of reserve connected to the system. Do >>> you have an idea how this can be performed in MOST? >>> >>> Thank you very much in advance >>> >>> Carlos Ferrandon >>> >>> >>> >>> -- >>> Carlos Ferrandon >> >> >> >> -- >> Carlos Ferrandon > > > > -- > Carlos Ferrandon
