>From the Introduction section in the MOST User’s Manual ><http://www.pserc.cornell.edu/matpower/docs/MOST-manual-1.0.pdf> …
While the problem formulation is general and incorporates a full nonlinear AC network model, the current implementation is limited to DC power flow modeling of the network. Some work has been done on an AC implementation, but it is not yet ready for release. Maybe someday … but currently, MOST is limited to a DC model of the network. — Ray > On Nov 19, 2017, at 12:42 PM, Mohammed Alhajri <[email protected]> wrote: > > Hello all > > i wanted to compare between AC OPF results and DC OPF results using MOST > > My code for AC OPF is: > > ++++++++++++++++++++++++++++++++++++++++++++ > > define_constants; > verbose = 1; > mpopt = mpoption('verbose', verbose); > mpopt = mpoption(mpopt, 'out.gen', 1); > mpopt = mpoption(mpopt, 'model', 'AC'); > mpopt = mpoption(mpopt, 'opf.ac.solver', 'MIPS'); > mpopt = mpoption(mpopt, 'most.solver', 'MIPS'); > if ~verbose > mpopt = mpoption(mpopt, 'out.all', 0); > end > > casefile = 'ex_case3b'; > mpc = loadcase(casefile); > xgd = loadxgendata('ex_xgd_ramp', mpc); > profiles = getprofiles('ex_load_profile'); > nt = size(profiles(1).values, 1); % number of periods > > mdi = loadmd(mpc, nt, xgd, [], [], profiles); > mdo = most(mdi, mpopt); > EPg1 = mdo.results.ExpectedDispatch; % active generation > Elam1 = mdo.results.GenPrices; % nodal energy price > if verbose > most_summary(mdo); > end > > > ++++++++++++++++++++++++++++++++++ > > For DC OPF > > ++++++++++++++++++++++++++++++++++ > > %% set up options > define_constants; > verbose = 1; > mpopt = mpoption('verbose', verbose); > mpopt = mpoption(mpopt, 'out.gen', 1); > mpopt = mpoption(mpopt, 'model', 'DC'); > mpopt = mpoption(mpopt, 'opf.dc.solver', 'MIPS'); > mpopt = mpoption(mpopt, 'most.solver', 'MIPS'); > mpopt = mpoption(mpopt, 'most.dc_model', 1); > if ~verbose > mpopt = mpoption(mpopt, 'out.all', 0); > end > > casefile = 'ex_case3b'; > mpc = loadcase(casefile); > xgd = loadxgendata('ex_xgd_ramp', mpc); > profiles = getprofiles('ex_load_profile'); > nt = size(profiles(1).values, 1); % number of periods > > %%----- Multiperiod DC OPF (w/ramp) ----- > mdi = loadmd(mpc, nt, xgd, [], [], profiles); > mdo = most(mdi, mpopt); > EPg1 = mdo.results.ExpectedDispatch; % active generation > Elam1 = mdo.results.GenPrices; % nodal energy price > if verbose > most_summary(mdo); > end > > +++++++++++++++++++++++++++++++ > > > i suppose to get in the results, isn't it? > > the results are attached, i don't find any difference between them! > <DC OPF.pdf><AC OPF.pdf>
