Dear Prof. Zimmerman,
After replacing the new smartmkt code I get an error as follows for the
following code:
define_constants;
mpc0 = loadcase('case6ww');
mpc = load2disp(mpc0);
mpc.branch(:, RATE_A) = 0; % disable flow limits
offers.P.qty = [200; 150; 180];
offers.P.prc = [30; 40; 50];
bids.P.qty = [70; 70; 70];
bids.P.prc = [100; 100; 100];
mkt.OPF = 'DC';
[r, co, cb] = runmarket(mpc, offers, bids, mkt);
co.P.prc
Shutting down generator 3.
??? Index exceeds matrix dimensions.
Error in ==> smartmkt at 206
qgcost = mpc.gencost(ng+(1:ng), :);
Error in ==> runmarket at 141
[co, cb, r, dispatch, success] = smartmkt(mpc, offers, bids, mkt, mpopt);
What is the reason?
Best Regards
D. Xu
On Thu, Mar 1, 2012 at 17:33, Ray Zimmerman <[email protected]> wrote:
> Sorry, the problem may have been due to a bug in smartmkt.m. I realized
> based on an error experienced by someone else that runmarket had never been
> tested with reactive power *costs*, only with reactive bids and offers.
> I've attached a new version of smartmkt.m that should hopefully fix that
> bug.
>
>
>
>
> --
> Ray Zimmerman
> Senior Research Associate
> 419A Warren Hall, Cornell University, Ithaca, NY 14853
> phone: (607) 255-9645
>
>
>
>
> On Feb 29, 2012, at 7:06 PM, Dailan Xu wrote:
>
> Dear Prof. Zimmerman,
>
> I have trusted you and I have sent you the m.files also I don't want
> distribute these files and I sent them only to you.
>
> I tried your suggestions but again it doesn't work. I don't want you to
> debug the errors. The main objective is this as I want to have active and
> reactive offers and bids, therefore, want to maximize social welfare and
> consequently want to have reactive costs. I am not sure either it is
> necessary or not.
>
> It is my pleasure working with you.
>
> Best Regards
>
> Dailan
>
> On Wed, Feb 29, 2012 at 17:21, Ray Zimmerman <[email protected]> wrote:
>
>> You are setting your generator costs to be piecewise linear costs with 3
>> points defining the curve (6 more columns) but then you only provide 3 more
>> columns (which look like they are polynomial coefficients).
>>
>> But honestly, I really can't afford to keep providing this level of
>> debugging support for you. If you need further assistance of this type, we
>> can discuss a consulting arrangement and I'll be happy to work with you.
>> You can also feel free to post on the MATPOWER-L and if anyone else has
>> time to help that will be up to them.
>>
>> --
>> Ray Zimmerman
>> Senior Research Associate
>> 419A Warren Hall, Cornell University, Ithaca, NY 14853
>> phone: (607) 255-9645
>>
>>
>>
>> On Feb 29, 2012, at 10:46 AM, Dailan Xu wrote:
>>
>> Dear Prof. Zimmerman,
>>
>> I want to use reactive cost as objective function, i.e both and active
>> and reactive costs. When I do runmarket also runopf I get the following
>> error. The files are attached.
>>
>> Error in ==> totcost at 53
>> pp = mkpp(x(i, 1:gencost(i, NCOST))', [(y(i,j2) - y(i,j1)) ./
>> (x(i,j2) -
>> x(i,j1)); y(i,j1)]');
>>
>> Error in ==> smartmkt at 191
>> fcost = mkt.t * totcost(mpc.gencost, zeros(ng, 1) ); %% fixed costs
>>
>> Error in ==> runmarket at 141
>> [co, cb, r, dispatch, success] = smartmkt(mpc, offers, bids, mkt, mpopt);
>>
>> Error in ==> ANM_9 at 299
>> [r, co, cb, f, dispatch, success, et] =
>> runmarket(mpc,offers,bids,mkt,mpoption);
>>
>>
>> Best Regards
>>
>> Dailan
>>
>> <ANM_9.m><case84_newdispatch6_10.m>
>>
>>
>>
>
>
>
