Not without writing your own code to do it. Check the MATPOWER-L archives for
multi-period OPF. For example …
https://www.mail-archive.com/[email protected]/msg00702.html
<https://www.mail-archive.com/[email protected]/msg00702.html>
— Ray
> On Nov 21, 2017, at 11:34 AM, Mohammed Alhajri <[email protected]> wrote:
>
> Actually our target is to do an economic dispatch for a given network in
> presence of Solar PV Power
>
> We want to do the ED for multiple periods,24hr, and we want to use AC OPF to
> compare the simulation results with the practical results
>
>
> is there any way to do this task?
>
> بتاريخ ٢٠١٧/١١/٢١ ٨:١٧ م، كتب "Mohammed Alhajri" <[email protected]
> <mailto:[email protected]>>:
> Thank you Ray for your clarification
>
> بتاريخ ٢٠١٧/١١/٢١ ٧:٠٦ م، كتب "Ray Zimmerman" <[email protected]
> <mailto:[email protected]>>:
> 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]
>> <mailto:[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>
>
