You need to add the capability curves to the gen matrix (columns PC1 through
QC2MAX, see section 5.4.3 in the User's Manual). Your Q bids also need to have
either zero or 9 rows (see the error message).
--
Ray Zimmerman
Senior Research Associate
419A Warren Hall, Cornell University, Ithaca, NY 14853
phone: (607) 255-9645
On Feb 12, 2012, at 6:52 PM, Carlos Gonzalez Almeida wrote:
> Dear Dr. Zimmerman,
>
> In the case "'t_auction_case'", I want to have reactive offers and bids as
> well as capability curve limits. What changes I have to apply to the gencost
> and gen matrices?
> I faced with following error when I ran the following code.
>
> clear all;
> mpc = loadcase('t_auction_case');
> mkt.OPF = 'AC';
> mkt.auction_type = 0;
> offers.P.qty = [ ...
> 12 24 24;
> 12 24 24;
> 12 24 24;
> 12 24 24;
> 12 24 24;
> 12 24 24 ];
> offers.P.prc = [ ...
> 20 50 60;
> 20 40 70;
> 20 42 80;
> 20 44 90;
> 20 46 75;
> 20 48 60 ];
> bids.P.qty = [ ...
> 10 10 10;
> 10 10 10;
> 10 10 10 ];
> bids.P.prc = [ ...
> 100 70 60;
> 100 50 20;
> 100 60 50 ];
>
> %++++++++Reactive bids and offers++++++++
> offers.Q.qty = [ ...
> 12 24 24;
> 12 24 24;
> 12 24 24;
> 12 24 24;
> 12 24 24;
> 12 24 24 ];
> offers.Q.prc = [ ...
> 20 50 60;
> 20 40 70;
> 20 42 80;
> 20 44 90;
> 20 46 75;
> 20 48 60 ];
> bids.Q.qty = [ ...
> 10 10 10;
> 10 10 10;
> 10 10 10 ];
> bids.Q.prc = [ ...
> 100 70 60;
> 100 50 20;
> 100 60 50 ];
> %++++++++++++++++++++++++
>
> [mpc_out, co, cb, f, dispatch, success, et] = runmarket(mpc, offers, bids,
> mkt);
>
> Error
> ??? Error using ==> off2case>idx_vecs at 459
> number of reactive power bids must be zero or match the total number of rows
> in gen
>
> Error in ==> off2case at 102
> [idxPo, idxPb, idxQo, idxQb] = idx_vecs(offers, bids, G, L, haveQ);
>
> Error in ==> smartmkt at 74
> [gen, genoffer] = off2case(mpc.gen, mpc.gencost, offers, bids, mkt.lim);
>
> Error in ==> runmarket at 141
> [co, cb, r, dispatch, success] = smartmkt(mpc, offers, bids, mkt, mpopt);
>
> Regards
>
> Carlos
>
