search

MATPOWER / mp-opt-model

MP-Opt-Model
Other
9 stars 3 forks source link

matpower 7.1, most, ex_case3b.m #3

Closed ajaythakur01 closed 3 years ago

ajaythakur01 commented 3 years ago

In ex_case3b.m, while run this program , there is some errors which is not predictable, please provide me some description about the given errors :- Error using miqps_ot (line 347) miqps_ot: EXITFLAG from price computation stage = -2

Error in miqps_master (line 218) [x, f, eflag, output, lambda] = ...

Error in opt_model/solve (line 189) [x, f, eflag, output, lambda] = ...

Error in most (line 2014) [mdo.QP.x, mdo.QP.f, mdo.QP.exitflag, mdo.QP.output, mdo.QP.lambda ] = ...

Error in TestProgram_MPUC_1 (line 21) mdo = most(mdi,mpopt);

ajaythakur01 commented 3 years ago

I am working on Deterministic Unit Commitment, after running this program there are some errors. this is my program: close all; clc; casefile = 'ex_case3b'; mpc = loadcase(casefile); xgd = loadxgendata('ex_xgd_uc',mpc); [iwind,mpc,xgd] = addwind('ex_wind_uc',mpc,xgd); profiles = getprofiles('ex_wind_profile_d',iwind); profiles = getprofiles('ex_load_profile',profiles); nt=size(profiles(1).values,1); %number of periods mpc_full = mpc; %save for later xgd_full = xgd; %save for later STARTUP = 2; SHUTDOWN = 3; mpc.gencost(:,[STARTUP SHUTDOWN])= 0; %remove startup/shutdown costs xgd.MinUp(2) = 1; %remove min up-time constraint xgd.PositiveLoadFollowReserveQuantity(3)= 250; %remove ramp reserve xgd.PositiveLoadFollowReservePrice(3)= 1e-6; %constraint and costs xgd.NegativeLoadFollowReservePrice(3)= 1e-6; mpopt = mpoption('out.all', 0, 'verbose', 2); mpopt = mpoption(mpopt,'most.dc_model',1); %use DC network model(default) mdi = loadmd(mpc,nt,xgd,[],[],profiles); mdo = most(mdi,mpopt); ms = most_summary(mdo); %print results,depending on 'vebrose'option mpc.gencost(2,[STARTUP SHUTDOWN])=[200 200]; mpc.gencost(3,[STARTUP SHUTDOWN])=[3000 600]; %equivalent to doing: mpc = mpc_full mdi = loadmd(mpc,nt,xgd,[],[],profiles); mdo = most(mdi,mpopt); xgd.MinUp(2) = 3; mdi = loadmd(mpc,nt,xgd,[],[],profiles); mdo = most(mdi,mpopt); xgd.PositiveLoadFollowReserveQuatity(3)= 100; %restore ramp reserve xgd.PositiveLoadFollowReservePrice(3)= 10; %constraint and costs xgd.NegativeLoadFollowReservePrice(3)= 10; %equivalent to doing:xgd = xgd_full; mdi = loadmd(mpc,nt,xgd,[],[],profiles); %mpopt = mpoption(mpopt,'most.storage.cyclic',1); [iess,mpc,xgd,sd] = addstorage('ex_storage',mpc,xgd); mdo = most(mdi,mpopt);

Errors are: Error using miqps_ot (line 347) miqps_ot: EXITFLAG from price computation stage = -2

Error in miqps_master (line 218) [x, f, eflag, output, lambda] = ...

Error in opt_model/solve (line 189) [x, f, eflag, output, lambda] = ...

Error in most (line 2014) [mdo.QP.x, mdo.QP.f, mdo.QP.exitflag, mdo.QP.output, mdo.QP.lambda ] = ...

Error in TestProgram_MPUC_1 (line 21) mdo = most(mdi,mpopt);

rdzman commented 3 years ago

See MATPOWER/most#21.