Dear Stephanie,
Please see the Gurobi documentation for the meaning of the -4 exit code.
According to the Status Codes on p. 621 of
http://www.gurobi.com/documentation/7.0/refman.pdf, it looks like it means
infeasible or unbounded. My guess is that the problem is somehow infeasible.
There are many possible reasons it could be infeasible, either because the
problem for an individual hour is infeasible (e.g. insufficient generation to
meet demand) or the hour-to-hour constraints make it infeasible (e.g. ramp
rates too small).
I would probably begin by trying to construct the 24 individual single hour
problems to see if they solve. If so, try relaxing ramp rates.
Also, it seems as though you figured out how to set the load profile, but just
in case, here’s the response I wrote to the question you sent me off-list the
other day, in case it’s useful …
You definitely do not need to call getprofiles() 20 times. You simply need to
define a single profile that modifies all of the loads as desired. The load
profile defined in ex_load_profile.m is already set to do this ('rows' = 0,
'col' = CT_LOAD_ALL_PQ). It uses a target quantity ('chgtype' = CT_REP) rather
than a direct scale factor ('chgtype' = CT_REL). To understand the meaning of
the various fields in a profile, you need to become familiar, first, with
apply_changes()
<http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/apply_changes.html>
which is used to implement profiles, and then with apply_profile()
<http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/most/apply_profile.html>
(see also idx_ct
<http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/idx_ct.html> and
idx_profile
<http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/most/idx_profile.html>).
More information is available in Section 9.3.5 in the MATPOWER User’s Manual
<http://www.pserc.cornell.edu/matpower/docs/MATPOWER-manual-6.0.pdf> and
Section 5.1.6 in the MOST User’s Manual
<http://www.pserc.cornell.edu/matpower/docs/MOST-manual-1.0.pdf>.
Hope this helps,
Ray
> On Apr 25, 2017, at 1:22 PM, 陈丽霞 <[email protected]> wrote:
>
> Dear professor:
>
> I try to use matpower to calculate the Unit Commitment problem of a load
> profile of 24h in IEEE30 bus. And I refer to the example of 3 buses in
> most_ex6_uc.m, and change the case file to case30 and write the code in the
> following way to set the load profile:
>
> profiles = getprofiles('load_profile');
>
> profiles.rows=load_index';
>
> profiles.values(:,1,:)=PLOAD’;
>
> where load_index is the bus index which loads are in, and PLOAD is the load
> of each bus at each hour, and I just run the Base:No Network situation. But
> it turns out that:
>
> MOST: MIQP solver 'GUROBI' failed with exit flag = -4
>
> You can query the workspace to debug.
>
> When finished, type the word "return" to continue.
>
> I don’t know how this kind of error will appear, has anyone ever met this
> problem? Is there anything else need to be set?
>
> Attached are my codes, I would be most appreciated if you can reply. Thanks
> in advance.
>
>
>
> Best wishes
>
> Sincerely, yours
>
> Stephanie
>
> <2017.txt><load_profile.m><xgd_uc.m><2017.txt>
