a) It would take some time for me to dig into the details (Carlos, maybe you remember more off the top of your head), but it has to do with the fact that the constant term (“no-load cost" is only included in the cost if the unit is committed, so it is actually removed from the part that is a function of Pg and handled separately. I believe, if I am not mistaken, the constant term for generator i in period t, scenario j, contingency k, can be found in mdo.flow(t,j,k).mpc.fixed_gencost(i). And I believe the probability-weighted "no-load" cost for generator i in period t can be found in mdo.UC.c00(i,t).
b) It is possible. Simply change the col field of the profile from CT_TAREALOAD to CT_TLOAD. See Table 9-2 in the MATPOWER User’s Manual <http://www.pserc.cornell.edu/matpower/docs/MATPOWER-manual-6.0.pdf>, in the description of MATPOWER’s apply_changes() <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/apply_changes.html> function, which is used to implement the profiles in MOST. Hope this helps, Ray > On Mar 19, 2018, at 8:09 PM, Bainan Xia <[email protected]> wrote: > > Hi guys, > > I’ve got following two questions when solving multi-period UC problem. > > a) Given the result structure ‘mdo’, I’m calculating the total > generation cost for time period t in following way: > mpc = mdo.flow(t).mpc; > total_cost = sum(totcost(mpc.gencost, mpc.gen(:, PG))); > However, I found the total cost is much less than my expectation. Then I > compared mdo.flow(t).mpc.gencost and mpc.gencost and found all the constant > terms (in my case, cost curves are quadratic and last column defines the > constant term) are zero in the previous matrix, which is different from the > base case mpc structure. I’m wondering what causes the issue. > > b) Previously, from the example cases (ex_load_profile), I’ve learned > how to assign the time series load profile for each area/load zone. Now, I > would like to assign specific time series load profile for each bus and > cannot find an example doing so. I’m wondering whether there is a way to do > that. Otherwise, I may use the naïve method: making the #area/#load zone the > same as the #bus, so that I could assign the load profile in the old way. > > Thanks for the help! > > Best, > Bainan
