Re:
4 - NaN - - - - -
17 36 - NaN - - - - -
18 40 - NaN - - - - -
19 42 - NaN - - - - -
20 46 0.00 NaN NaN 0.00 0.00 0.00
NaN
21 49 0.00 NaN NaN 0.00 0.00 0.00
NaN
22 54 0.00 NaN NaN 0.00 0.00 0.00
NaN
23 55 - NaN - - - - -
24 56 - NaN - - - - -
25 59 0.00 NaN NaN 0.00 0.00 0.00
NaN
26 61 0.00 NaN NaN 0.00 0.00 0.00
NaN
27 62 - NaN - - - - -
28 65 0.00 NaN NaN 0.00 0.00 0.00
NaN
29 66 0.00 NaN NaN 0.00 0.00 0.00
NaN
30 69 0.00 NaN NaN 0.00 0.00 0.00
NaN
31 70 - NaN - - - - -
32 72 - NaN - - - - -
33 73 - NaN - - - - -
34 74 - NaN - - - - -
35 76 - NaN - - - - -
36 77 - NaN - - - - -
37 80 0.00 NaN NaN 0.00 0.00 0.00
NaN
38 85 - NaN - - - - -
39 87 0.00 NaN NaN 0.00 0.00 0.00
NaN
40 89 0.00 NaN NaN 0.00 0.00 0.00
NaN
41 90 - NaN - - - - -
42 91 - NaN - - - - -
43 92 - NaN - - - - -
44 99 - NaN - - - - -
45 100 0.00 NaN NaN 0.00 0.00 0.00
NaN
46 103 0.00 NaN NaN 0.00 0.00 0.00
NaN
47 104 - NaN - - - - -
48 105 - NaN - - - - -
49 107 - NaN - - - - -
50 110 - NaN - - - - -
51 111 0.00 NaN NaN 0.00 0.00 0.00
NaN
52 112 - NaN - - - - -
53 113 - NaN - - - - -
54 116 - NaN - - - - -
557 0.00 NaN NaN 0.00 0.00 0.00
NaN
56 30 0.00 NaN NaN 0.00 0.00 0.00
NaN
57 38 0.00 NaN NaN 0.00 0.00 0.00
NaN
58 44 0.00 NaN NaN 0.00 0.00 0.00
NaN
59 51 0.00 NaN NaN 0.00 0.00 0.00
NaN
60 53 0.00 NaN NaN 0.00 0.00 0.00
NaN
61 68 0.00 NaN NaN 0.00 0.00 0.00
NaN
62 82 0.00 NaN NaN 0.00 0.00 0.00
NaN
63 83 0.00 NaN NaN 0.00 0.00 0.00
NaN
64 94 0.00 NaN NaN 0.00 0.00 0.00
NaN
65 97 0.00 NaN NaN 0.00 0.00 0.00
NaN
Total: 0.00 NaN 0.00 0.00 0.00
NaN
>>>>> Did NOT converge (1.45 seconds) <<<<<
On Mon, Feb 6, 2017 at 1:44 PM, Ammar Raza wrote:
> how to modify case30 accouding to auction case?
>
Re:
hello all,
My tutor told me to convert Case118 according to Auction case.
I modify case118 by adding dispatchable load and did runmarket() and the
programme is converging. But when I writing my program in the matrix I am
not getting the converging form.
Find the attachment.
2017-02-09 13:26 GMT+00:00 Ammar Raza :
>
> Hi
>
> This is Case118. And I am trying to solve as per auction.
>
>
> mkt.OPF = 'AC';
> mkt.auction_type = 1;
> mpc = loadcase('case118_3');
> offers.P.qty = [0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 100 150 170;
> 20 25 30;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 60 70 80;
> 70 80 100;
> 0 0 0;
> 2 2 3;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 5 7 7;
> 50 60 80;
> 10 14 20;
> 0 0 0;
> 0 0 0;
> 40 50 60;
> 40 50 70;
> 0 0 0;
> 100 130 150;
> 100 130 150;
> 120 180 200;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 120 160 180;
> 0 0 0;
> 1 1 2;
> 150 200 220;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 60 80 110;
> 10 15 15
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
> 8 10 15;
> 0 0 0;
> 0 0 0;
> 0 0 0];
>
> offers.P.prc = [0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
>20 50 60;
> 50 60;
>0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
>20 50 60;
>20 50 60;
>0 0 0;
>20 50 60;
>0 0 0;
>0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
>20 50 60;
>20 50 60;
>20 50 60;
>0 0 0;
>0 0 0;
>20 50 60;
>20 50 60;
>0 0 0;
>20 50 60;
> 20 50 60;
> 20 50 60;
>0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
>0 0 0;
>0 0 0;
>20 50 60;
>0 0 0;
>20 50 60;
>20 50 60;
>0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
>20 50 60;
>20 50 60;
>0 0 0;
> 0 0 0;
> 0 0 0;
> 0 0 0;
>20 50 60;
>0 0 0;
> 0 0 0;
> 0 0 0];
>
> bids.P.qty = [10 10 10;
> 10 10 10;
> 10 10 10;
> 10 10 10;
> 10 10 10;
> 10 10 10;
> 10 10 10;
> 10 10 10;
> 10 10 10;
> 10 10 10;
> 10 10 10];
>
> bids.P.prc = [80 60 50;
> 80 60 50;
> 80 60 50;
> 80 50 20;
> 80 60 50;
> 80 60 50;
> 80 50 20;
> 80 60 50;
> 80 50 20;
> 80 50 20;
> 80 60 50];
>
> [mpc_out, co, cb, f, dispatch, success, et] = runmarket(mpc, offers, bids,
> mkt);
>
>
>
>
>
>
> I am Getting this
>
>
> SMARTMARKET: non-convergent UOPF
>
>
> | Market Summary
> |
>
>
> Dispatch period duration: 1.00 hours
> Gen Bus Pg PriceRevenue Fix+Var Strt/Stp Total
> Earnings
> ## (MW)($/MWh) ($) Cost ($) Cost ($) Cost ($)
> ($)
> --- --- ---
Smart market program
Hello all, I am working on a smart market of a case 30 buses and I wrote a program. can you please tell me where I did wrong? thank you Kind Regard, Mohammad Ammar Raza example7.m Description: Binary data
Re: Smart market program
please check the code I have already did it. Kind Regard, Mohammad Ammar Raza On 16 February 2017 at 11:33, 严超 <1499730...@qq.com> wrote: > I think you should upload the runmarket function to check the correcness > of the code. > > > -- 原始邮件 -- > *发件人:* "Ammar Raza";; > *发送时间:* 2017年2月16日(星期四) 晚上7:29 > *收件人:* "matpower-l"; > *主题:* Smart market program > > Hello all, > > I am working on a smart market of a case 30 buses and I wrote a program. > can you please tell me where I did wrong? > > thank you > Kind Regard, > > Mohammad Ammar Raza > > > >
Re: Smart market program
can you please look at it. Kind Regard, Mohammad Ammar Raza On 16 February 2017 at 14:44, Ray Zimmerman wrote: > I don’t have the case file you are using so I can’t run your example. > Also, I don’t have time to debug your code. If you can reduce it to a > single call to runmarket() that is not behaving as you expect, I’m happy > to attempt to clear up any questions you may have. > >Ray > > > > On Feb 16, 2017, at 6:29 AM, Ammar Raza wrote: > > Hello all, > > I am working on a smart market of a case 30 buses and I wrote a program. > can you please tell me where I did wrong? > > thank you > Kind Regard, > > Mohammad Ammar Raza > > > > > > > example7.m Description: Binary data case30_1.m Description: Binary data
Re: Smart market program
Actually, I am wondering why I am getting Negative Earning in my output? Kind Regard, Mohammad Ammar Raza On 16 February 2017 at 16:32, Ray Zimmerman wrote: > I have no idea whether the results are what you expect, but it runs fine > for me. If you want me to look at it further, please reduce it to a single > call to runmarket() (with a specific set of inputs) so that I can > reproduce the error. Also, make sure to describe the result you get (what > error message if it gives you an error) and what you expected to get. > > Ray > > > On Feb 16, 2017, at 9:53 AM, Ammar Raza wrote: > > can you please look at it. > > Kind Regard, > > Mohammad Ammar Raza > > > > > On 16 February 2017 at 14:44, Ray Zimmerman wrote: > >> I don’t have the case file you are using so I can’t run your example. >> Also, I don’t have time to debug your code. If you can reduce it to a >> single call to runmarket() that is not behaving as you expect, I’m happy >> to attempt to clear up any questions you may have. >> >>Ray >> >> >> >> On Feb 16, 2017, at 6:29 AM, Ammar Raza wrote: >> >> Hello all, >> >> I am working on a smart market of a case 30 buses and I wrote a program. >> can you please tell me where I did wrong? >> >> thank you >> Kind Regard, >> >> Mohammad Ammar Raza >> >> >> >> >> >> >> > > > >
Re: Smart market program
Thank You Kind Regard, Mohammad Ammar Raza On 16 February 2017 at 16:48, Ray Zimmerman wrote: > The generator offers (in offers.P.prc) are lower than their marginal cost > (in gencost). Selling below cost results in negative earnings. > >Ray > > > > On Feb 16, 2017, at 11:37 AM, Ammar Raza wrote: > > Actually, I am wondering why I am getting Negative Earning in my output? > > Kind Regard, > > Mohammad Ammar Raza > > > > > On 16 February 2017 at 16:32, Ray Zimmerman wrote: > >> I have no idea whether the results are what you expect, but it runs fine >> for me. If you want me to look at it further, please reduce it to a single >> call to runmarket() (with a specific set of inputs) so that I can >> reproduce the error. Also, make sure to describe the result you get (what >> error message if it gives you an error) and what you expected to get. >> >> Ray >> >> >> On Feb 16, 2017, at 9:53 AM, Ammar Raza wrote: >> >> can you please look at it. >> >> Kind Regard, >> >> Mohammad Ammar Raza >> >> >> >> >> On 16 February 2017 at 14:44, Ray Zimmerman wrote: >> >>> I don’t have the case file you are using so I can’t run your example. >>> Also, I don’t have time to debug your code. If you can reduce it to a >>> single call to runmarket() that is not behaving as you expect, I’m >>> happy to attempt to clear up any questions you may have. >>> >>>Ray >>> >>> >>> >>> On Feb 16, 2017, at 6:29 AM, Ammar Raza wrote: >>> >>> Hello all, >>> >>> I am working on a smart market of a case 30 buses and I wrote a program. >>> can you please tell me where I did wrong? >>> >>> thank you >>> Kind Regard, >>> >>> Mohammad Ammar Raza >>> >>> >>> >>> >>> >>> >>> >> >> >> >> > >
Re: Smart market program
Hello everyone, Can you please tell me Why my Social Welfare is always in the negative? Kind Regard, Mohammad Ammar Raza On 16 February 2017 at 17:02, Ammar Raza wrote: > Thank You > > Kind Regard, > > Mohammad Ammar Raza > > > > > On 16 February 2017 at 16:48, Ray Zimmerman wrote: > >> The generator offers (in offers.P.prc) are lower than their marginal >> cost (in gencost). Selling below cost results in negative earnings. >> >> Ray >> >> >> >> On Feb 16, 2017, at 11:37 AM, Ammar Raza wrote: >> >> Actually, I am wondering why I am getting Negative Earning in my output? >> >> Kind Regard, >> >> Mohammad Ammar Raza >> >> >> >> >> On 16 February 2017 at 16:32, Ray Zimmerman wrote: >> >>> I have no idea whether the results are what you expect, but it runs fine >>> for me. If you want me to look at it further, please reduce it to a single >>> call to runmarket() (with a specific set of inputs) so that I can >>> reproduce the error. Also, make sure to describe the result you get (what >>> error message if it gives you an error) and what you expected to get. >>> >>> Ray >>> >>> >>> On Feb 16, 2017, at 9:53 AM, Ammar Raza wrote: >>> >>> can you please look at it. >>> >>> Kind Regard, >>> >>> Mohammad Ammar Raza >>> >>> >>> >>> >>> On 16 February 2017 at 14:44, Ray Zimmerman wrote: >>> >>>> I don’t have the case file you are using so I can’t run your example. >>>> Also, I don’t have time to debug your code. If you can reduce it to a >>>> single call to runmarket() that is not behaving as you expect, I’m >>>> happy to attempt to clear up any questions you may have. >>>> >>>>Ray >>>> >>>> >>>> >>>> On Feb 16, 2017, at 6:29 AM, Ammar Raza wrote: >>>> >>>> Hello all, >>>> >>>> I am working on a smart market of a case 30 buses and I wrote a >>>> program. can you please tell me where I did wrong? >>>> >>>> thank you >>>> Kind Regard, >>>> >>>> Mohammad Ammar Raza >>>> >>>> >>>> >>>> >>>> >>>> >>>> >>> >>> >>> >>> >> >> > untitled.fig Description: Binary data
MOST
Hello All, I am working on MOST and I have gone through the multiperiod problem example 6 and 7. I would like to add market in that problem. Can you please suggest me with an example how to developed market in them. Kind Regard, Mohammad Ammar Raza
Re: MOST
Hello, Can you please specified how we apply offers(), and bids() in the MOST. In example 6 and 7.in the MOST manual. Thankyou Kind Regard, Mohammad Ammar Raza On 23 March 2017 at 11:30, Ammar Raza wrote: > Hello All, > > I am working on MOST and I have gone through the multiperiod problem > example 6 and 7. > > I would like to add market in that problem. > > Can you please suggest me with an example how to developed market in them. > > > Kind Regard, > > Mohammad Ammar Raza > > > >
Re: MOST
hi everyone, Can you please specified me how to put Reserve data in 118 bus system in MOST. In 118 bus there is 54 gen therefore, I have to put 1's in 54 times in mpc.reserves.zones(). I did it but I am not sure it is right. %%- Reserve Data -%% %% reserve zones, element i, j is 1 if gen j is in zone i, 0 otherwise mpc.reserves.zones = [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 11 1 1 1; ]; %% reserve requirements for each zone in MW mpc.reserves.req = 100; %% reserve costs in $/MW for each gen that belongs to at least 1 zone %% (same order as gens, but skipping any gen that does not belong to any zone) mpc.reserves.cost = [ 0; 0; 5; 5; 0; 5; 5; 0; 0; 0; 5; 5; 5; 5; 5; 0; 5; 5; 5; 0; 0; 0; 0; 5; 5; 5; 0; 0; 0; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5; 5;]; %% OPTIONAL max reserve quantities for each gen that belongs to at least 1 zone %% (same order as gens, but skipping any gen that does not belong to any zone) mpc.reserves.qty = [ 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50; 50;]; Kind Regard, Mohammad Ammar Raza On 28 March 2017 at 15:55, Ray Zimmerman wrote: > If you are asking for the equivalent of runmarket() > <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/extras/smartmarket/runmarket.html> > for > MOST, I’m afraid it hasn’t been implemented. MOST assumes that gencost > has been constructed to reflect the problem you want to solve. This could > be, for example, for a market, where the offers and bids have been > converted into the corresponding gencost. If the bids and offers are > constant from period to period throughout the planning horizon, then you > could use off2case() > <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/extras/smartmarket/off2case.html> > to > do the conversion for you. If they vary from period to period, you’ll have > to figure out the appropriate gen and gencost for each period (possibly > using off2case() > <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/extras/smartmarket/off2case.html> > and then create a profile (see Section 5.1.6 in the MOST User’s Manual > <http://www.pserc.cornell.edu/matpower/docs/MOST-manual-1.0.pdf>) to > implement the period-to-period modifications. > > Ray > > > > > On Mar 24, 2017, at 11:10 AM, Ammar Raza wrote: > > Hello, > > Can you please specified how we apply offers(), and bids() in the MOST. In > example 6 and 7.in the MOST manual. > > Thankyou > > Kind Regard, > > Mohammad Ammar Raza > > > > > On 23 March 2017 at 11:30, Ammar Raza wrote: > >> Hello All, >> >> I am working on MOST and I have gone through the multiperiod problem >> example 6 and 7. >> >> I would like to add market in that problem. >> >> Can you please suggest me with an example how to developed market in them. >> >> >> Kind Regard, >> >> Mohammad Ammar Raza >> >> >> >> > >
MOST
Hi every one,
I am working on most dcopf of 118 bus system but I am getting this error.
I attached my code with this email.
Can you please specified my why I am getting this error?
Thank you,
Kind Regard
Mohammad Ammar Raza
Index exceeds matrix dimensions.
Error in get_reorder (line 29)
B = subsref(A, s);
Error in e2i_data (line 68)
newval = get_reorder(val, idx, dim);
Error in e2i_field (line 60)
e2i_data(mpc, subsref(mpc, s), ordering, dim) );
Error in toggle_reserves>userfcn_reserves_ext2int (line 135)
mpc = e2i_field(mpc, {'reserves', 'zones'}, 'gen', 2);
Error in run_userfcn (line 37)
rv = feval(userfcn.(stage)(k).fcn, rv, args);
Error in ext2int (line 211)
mpc = run_userfcn(mpc.userfcn, 'ext2int', mpc);
Error in opf (line 195)
mpc = ext2int(mpc);
Error in runopf (line 75)
[r, success] = opf(casedata, mpopt);
Error in runopf_w_res (line 58)
[varargout{1:nargout}] = runopf(mpc, varargin{2:nargin});
Error in most_ammar2 (line 39)
r1 = runopf_w_res(mpc, mpopt);
Kind Regard,
Mohammad Ammar Raza
most_ammar2.m
Description: Binary data
ex_case118.m
Description: Binary data
MOST
Hi all, I am working on bus 118 with 54 generator and 99 loads. I would like to ask you for the multi period problem i.e. 24 hours, how I will make the load profile function. Can I do like that loadprofile.value (:,1,1) = [24 x 99]; As in case3a in the MOST manual there is only one load therefore, load profile matrix is [12 x 1] loadprofile.value(:,1,1)=[ 440; 480; .. .. .. ]; Kind Regard, Mohammad Ammar Raza
Re: MOST
that means loadprofile = struct( ... 'type', 'mpcData', ... 'table', CT_TLOAD, ... 'rows', 0, ... 'col', CT_LOAD_ALL_PQ, ... 'chgtype', CT_REP, ... 'values', [] ); loadprofile.values(:, 1, :) = [ 24 x 99 ] Kind Regard, Mohammad Ammar Raza On 9 May 2017 at 14:55, Stephanie wrote: > Hi, > > I tried it a few days before, what I did is: > > loadprofile.rows=load_index; %which applies the load to the bus it locates > loadprofile.value (:,1,:) = [24 x 99]; > > And it works in my situation > > Hope it will help you > -- > Best Regards > Stephanie > > > *From:* Ammar Raza > *Date:* 2017-05-09 21:43 > *To:* MATPOWER discussion forum > *Subject:* MOST > Hi all, > > I am working on bus 118 with 54 generator and 99 loads. > > I would like to ask you for the multi period problem i.e. 24 hours, how I > will make the load profile function. Can I do like that > > loadprofile.value (:,1,1) = [24 x 99]; > > > > As in case3a in the MOST manual there is only one load therefore, load > profile matrix is [12 x 1] > > loadprofile.value(:,1,1)=[ > 440; > 480; > .. > .. > .. ]; > > > Kind Regard, > > Mohammad Ammar Raza > > > >
Re: MOST
Its means loadprofile.rows= [1 x 99 ] Thanks Stephanie for your help. Kind Regard, Mohammad Ammar Raza On 9 May 2017 at 15:25, Stephanie wrote: > Yes, and the value of rows, 0 means the load is applied to the all buses, > so you should also change the value of the rows by loadprofile.rows=load_index > (e.g. loadprofile.rows=[2 3 5]), if there are three loads located in bus 2, > 3 and 5, respectively) > > Best Regards, > Stephanie > > 2017-05-09 22:09 GMT+08:00 Ammar Raza : > >> that means >> >> loadprofile = struct( ... >> 'type', 'mpcData', ... >> 'table', CT_TLOAD, ... >> 'rows', 0, ... >> 'col', CT_LOAD_ALL_PQ, ... >> 'chgtype', CT_REP, ... >> 'values', [] ); >> >> loadprofile.values(:, 1, :) = [ 24 x 99 ] >> >> >> Kind Regard, >> >> Mohammad Ammar Raza >> >> >> >> >> On 9 May 2017 at 14:55, Stephanie wrote: >> >>> Hi, >>> >>> I tried it a few days before, what I did is: >>> >>> loadprofile.rows=load_index; %which applies the load to the bus it >>> locates >>> loadprofile.value (:,1,:) = [24 x 99]; >>> >>> And it works in my situation >>> >>> Hope it will help you >>> -- >>> Best Regards >>> Stephanie >>> >>> >>> *From:* Ammar Raza >>> *Date:* 2017-05-09 21:43 >>> *To:* MATPOWER discussion forum >>> *Subject:* MOST >>> Hi all, >>> >>> I am working on bus 118 with 54 generator and 99 loads. >>> >>> I would like to ask you for the multi period problem i.e. 24 hours, how >>> I will make the load profile function. Can I do like that >>> >>> loadprofile.value (:,1,1) = [24 x 99]; >>> >>> >>> >>> As in case3a in the MOST manual there is only one load therefore, load >>> profile matrix is [12 x 1] >>> >>> loadprofile.value(:,1,1)=[ >>> 440; >>> 480; >>> .. >>> .. >>> .. ]; >>> >>> >>> Kind Regard, >>> >>> Mohammad Ammar Raza >>> >>> >>> >>> >> >
Re: MOST
Thank you Ray
Therefore, Transmat should be like that?
transmat = cell(1, nt);
T = [ 1x 99];
[transmat{:}] = deal(T * ones(1,1));
transmat{1} = T;
Kind Regard,
Mohammad Ammar Raza
On 11 May 2017 at 13:42, Ray Zimmerman wrote:
> Just to clarify, if loadprofile.rows = 0, then loadprofile.values(:, 1,
> :) = [24 x 1] and all loads are scaled together. Otherwise, the number of
> columns in loadprofile.values(:, 1, :) must equal the length of
> loadprofile.rows and each of the specified loads is scaled separately.
>
>Ray
>
>
> On May 9, 2017, at 10:29 AM, Ammar Raza wrote:
>
> Its means
>
> loadprofile.rows= [1 x 99 ]
>
> Thanks Stephanie for your help.
>
> Kind Regard,
>
> Mohammad Ammar Raza
>
>
>
>
> On 9 May 2017 at 15:25, Stephanie wrote:
>
>> Yes, and the value of rows, 0 means the load is applied to the all buses,
>> so you should also change the value of the rows by
>> loadprofile.rows=load_index
>> (e.g. loadprofile.rows=[2 3 5]), if there are three loads located in bus 2,
>> 3 and 5, respectively)
>>
>> Best Regards,
>> Stephanie
>>
>> 2017-05-09 22:09 GMT+08:00 Ammar Raza :
>>
>>> that means
>>>
>>> loadprofile = struct( ...
>>> 'type', 'mpcData', ...
>>> 'table', CT_TLOAD, ...
>>> 'rows', 0, ...
>>> 'col', CT_LOAD_ALL_PQ, ...
>>> 'chgtype', CT_REP, ...
>>> 'values', [] );
>>>
>>> loadprofile.values(:, 1, :) = [ 24 x 99 ]
>>>
>>>
>>> Kind Regard,
>>>
>>> Mohammad Ammar Raza
>>>
>>>
>>>
>>>
>>> On 9 May 2017 at 14:55, Stephanie wrote:
>>>
>>>> Hi,
>>>>
>>>> I tried it a few days before, what I did is:
>>>>
>>>> loadprofile.rows=load_index; %which applies the load to the bus it
>>>> locates
>>>> loadprofile.value (:,1,:) = [24 x 99];
>>>>
>>>> And it works in my situation
>>>>
>>>> Hope it will help you
>>>> --
>>>> Best Regards
>>>> Stephanie
>>>>
>>>>
>>>> *From:* Ammar Raza
>>>> *Date:* 2017-05-09 21:43
>>>> *To:* MATPOWER discussion forum
>>>> *Subject:* MOST
>>>> Hi all,
>>>>
>>>> I am working on bus 118 with 54 generator and 99 loads.
>>>>
>>>> I would like to ask you for the multi period problem i.e. 24 hours, how
>>>> I will make the load profile function. Can I do like that
>>>>
>>>> loadprofile.value (:,1,1) = [24 x 99];
>>>>
>>>>
>>>>
>>>> As in case3a in the MOST manual there is only one load therefore, load
>>>> profile matrix is [12 x 1]
>>>>
>>>> loadprofile.value(:,1,1)=[
>>>> 440;
>>>> 480;
>>>> ..
>>>> ..
>>>> .. ];
>>>>
>>>>
>>>> Kind Regard,
>>>>
>>>> Mohammad Ammar Raza
>>>>
>>>>
>>>>
>>>>
>>>
>>
>
>
Re: MOST
Hi
I am trying to do example 6 of the MOST user manual with Bus 118 and I got
this error.
Error using mpopt2qpopt (line 92)
mpopt2qpopt: Sorry, no solver available for MIQP models
Error in most (line 2065)
mdo.QP.opt = mpopt2qpopt(mpopt, model, 'most');
Error in most_ammar6 (line 108)
mdo = most(mdi, mpopt);
Kind Regard,
Mohammad Ammar Raza
On 15 May 2017 at 20:14, Ray Zimmerman wrote:
> In a deterministic problem, transmat can simply be a scalar with the
> number of periods. See the description of transmat in the help for
> loadmd()
> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/most/loadmd.html>
> .
>
>Ray
>
>
>
> On May 15, 2017, at 7:37 AM, Ammar Raza wrote:
>
> Thank you Ray
>
> Therefore, Transmat should be like that?
>
> transmat = cell(1, nt);
> T = [ 1x 99];
> [transmat{:}] = deal(T * ones(1,1));
> transmat{1} = T;
>
> Kind Regard,
>
> Mohammad Ammar Raza
>
>
>
>
> On 11 May 2017 at 13:42, Ray Zimmerman wrote:
>
>> Just to clarify, if loadprofile.rows = 0, then loadprofile.values(:, 1,
>> :) = [24 x 1] and all loads are scaled together. Otherwise, the number
>> of columns in loadprofile.values(:, 1, :) must equal the length of
>> loadprofile.rows and each of the specified loads is scaled separately.
>>
>>Ray
>>
>>
>> On May 9, 2017, at 10:29 AM, Ammar Raza wrote:
>>
>> Its means
>>
>> loadprofile.rows= [1 x 99 ]
>>
>> Thanks Stephanie for your help.
>>
>> Kind Regard,
>>
>> Mohammad Ammar Raza
>>
>>
>>
>>
>> On 9 May 2017 at 15:25, Stephanie wrote:
>>
>>> Yes, and the value of rows, 0 means the load is applied to the all
>>> buses, so you should also change the value of the rows by
>>> loadprofile.rows=load_index
>>> (e.g. loadprofile.rows=[2 3 5]), if there are three loads located in bus 2,
>>> 3 and 5, respectively)
>>>
>>> Best Regards,
>>> Stephanie
>>>
>>> 2017-05-09 22:09 GMT+08:00 Ammar Raza :
>>>
>>>> that means
>>>>
>>>> loadprofile = struct( ...
>>>> 'type', 'mpcData', ...
>>>> 'table', CT_TLOAD, ...
>>>> 'rows', 0, ...
>>>> 'col', CT_LOAD_ALL_PQ, ...
>>>> 'chgtype', CT_REP, ...
>>>> 'values', [] );
>>>>
>>>> loadprofile.values(:, 1, :) = [ 24 x 99 ]
>>>>
>>>>
>>>> Kind Regard,
>>>>
>>>> Mohammad Ammar Raza
>>>>
>>>>
>>>>
>>>>
>>>> On 9 May 2017 at 14:55, Stephanie wrote:
>>>>
>>>>> Hi,
>>>>>
>>>>> I tried it a few days before, what I did is:
>>>>>
>>>>> loadprofile.rows=load_index; %which applies the load to the bus it
>>>>> locates
>>>>> loadprofile.value (:,1,:) = [24 x 99];
>>>>>
>>>>> And it works in my situation
>>>>>
>>>>> Hope it will help you
>>>>> --
>>>>> Best Regards
>>>>> Stephanie
>>>>>
>>>>>
>>>>> *From:* Ammar Raza
>>>>> *Date:* 2017-05-09 21:43
>>>>> *To:* MATPOWER discussion forum
>>>>> *Subject:* MOST
>>>>> Hi all,
>>>>>
>>>>> I am working on bus 118 with 54 generator and 99 loads.
>>>>>
>>>>> I would like to ask you for the multi period problem i.e. 24 hours,
>>>>> how I will make the load profile function. Can I do like that
>>>>>
>>>>> loadprofile.value (:,1,1) = [24 x 99];
>>>>>
>>>>>
>>>>>
>>>>> As in case3a in the MOST manual there is only one load therefore, load
>>>>> profile matrix is [12 x 1]
>>>>>
>>>>> loadprofile.value(:,1,1)=[
>>>>> 440;
>>>>> 480;
>>>>> ..
>>>>> ..
>>>>> .. ];
>>>>>
>>>>>
>>>>> Kind Regard,
>>>>>
>>>>> Mohammad Ammar Raza
>>>>>
>>>>>
>>>>>
>>>>>
>>>>
>>>
>>
>>
>
>
Re: MOST
Hi Ray,
Thank you for your email.
Kind Regard,
Mohammad Ammar Raza
On 17 May 2017 at 14:51, Ray Zimmerman wrote:
> Hi Mohammad,
>
> A unit commitment problem with quadratic generator costs or ramping costs
> requires a solver that can handle Mixed Integer Quadratic Programming
> (MIQP) problems. The only MATPOWER-compatible solvers that handle MIQP are
> Gurobi, CPLEX and MOSEK, as stated in the System Requirements section of
> the MOST User’s Manual.
>
> So, you will either need to install one of these solvers, or modify your
> problem so that in does not include any quadratic costs.
>
>Ray
>
>
> On May 16, 2017, at 5:30 PM, Ammar Raza wrote:
>
> Hi
>
> I am trying to do example 6 of the MOST user manual with Bus 118 and I got
> this error.
>
> Error using mpopt2qpopt (line 92)
> mpopt2qpopt: Sorry, no solver available for MIQP models
>
> Error in most (line 2065)
> mdo.QP.opt = mpopt2qpopt(mpopt, model, 'most');
>
> Error in most_ammar6 (line 108)
> mdo = most(mdi, mpopt);
>
>
> Kind Regard,
>
> Mohammad Ammar Raza
>
>
>
>
> On 15 May 2017 at 20:14, Ray Zimmerman wrote:
>
>> In a deterministic problem, transmat can simply be a scalar with the
>> number of periods. See the description of transmat in the help for
>> loadmd()
>> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/most/loadmd.html>
>> .
>>
>>Ray
>>
>>
>>
>> On May 15, 2017, at 7:37 AM, Ammar Raza wrote:
>>
>> Thank you Ray
>>
>> Therefore, Transmat should be like that?
>>
>> transmat = cell(1, nt);
>> T = [ 1x 99];
>> [transmat{:}] = deal(T * ones(1,1));
>> transmat{1} = T;
>>
>> Kind Regard,
>>
>> Mohammad Ammar Raza
>>
>>
>>
>>
>> On 11 May 2017 at 13:42, Ray Zimmerman wrote:
>>
>>> Just to clarify, if loadprofile.rows = 0, then loadprofile.values(:,
>>> 1, :) = [24 x 1] and all loads are scaled together. Otherwise, the
>>> number of columns in loadprofile.values(:, 1, :) must equal the length
>>> of loadprofile.rows and each of the specified loads is scaled
>>> separately.
>>>
>>>Ray
>>>
>>>
>>> On May 9, 2017, at 10:29 AM, Ammar Raza wrote:
>>>
>>> Its means
>>>
>>> loadprofile.rows= [1 x 99 ]
>>>
>>> Thanks Stephanie for your help.
>>>
>>> Kind Regard,
>>>
>>> Mohammad Ammar Raza
>>>
>>>
>>>
>>>
>>> On 9 May 2017 at 15:25, Stephanie wrote:
>>>
>>>> Yes, and the value of rows, 0 means the load is applied to the all
>>>> buses, so you should also change the value of the rows by
>>>> loadprofile.rows=load_index
>>>> (e.g. loadprofile.rows=[2 3 5]), if there are three loads located in bus 2,
>>>> 3 and 5, respectively)
>>>>
>>>> Best Regards,
>>>> Stephanie
>>>>
>>>> 2017-05-09 22:09 GMT+08:00 Ammar Raza :
>>>>
>>>>> that means
>>>>>
>>>>> loadprofile = struct( ...
>>>>> 'type', 'mpcData', ...
>>>>> 'table', CT_TLOAD, ...
>>>>> 'rows', 0, ...
>>>>> 'col', CT_LOAD_ALL_PQ, ...
>>>>> 'chgtype', CT_REP, ...
>>>>> 'values', [] );
>>>>>
>>>>> loadprofile.values(:, 1, :) = [ 24 x 99 ]
>>>>>
>>>>>
>>>>> Kind Regard,
>>>>>
>>>>> Mohammad Ammar Raza
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On 9 May 2017 at 14:55, Stephanie wrote:
>>>>>
>>>>>> Hi,
>>>>>>
>>>>>> I tried it a few days before, what I did is:
>>>>>>
>>>>>> loadprofile.rows=load_index; %which applies the load to the bus it
>>>>>> locates
>>>>>> loadprofile.value (:,1,:) = [24 x 99];
>>>>>>
>>>>>> And it works in my situation
>>>>>>
>>>>>> Hope it will help you
>>>>>> --
>>>>>> Best Regards
>>>>>> Stephanie
>>>>>>
>>>>>>
>>>>>> *From:* Ammar Raza
>>>>>> *Date:* 2017-05-09 21:43
>>>>>> *To:* MATPOWER discussion forum
>>>>>> *Subject:* MOST
>>>>>> Hi all,
>>>>>>
>>>>>> I am working on bus 118 with 54 generator and 99 loads.
>>>>>>
>>>>>> I would like to ask you for the multi period problem i.e. 24 hours,
>>>>>> how I will make the load profile function. Can I do like that
>>>>>>
>>>>>> loadprofile.value (:,1,1) = [24 x 99];
>>>>>>
>>>>>>
>>>>>>
>>>>>> As in case3a in the MOST manual there is only one load therefore,
>>>>>> load profile matrix is [12 x 1]
>>>>>>
>>>>>> loadprofile.value(:,1,1)=[
>>>>>> 440;
>>>>>> 480;
>>>>>> ..
>>>>>> ..
>>>>>> .. ];
>>>>>>
>>>>>>
>>>>>> Kind Regard,
>>>>>>
>>>>>> Mohammad Ammar Raza
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>
>>>
>>>
>>
>>
>
>
Re: MOST
hi,
I am wondering in the MATLAB R2015a version my code is running fine but in
MATLAB R2016a version I am getting error of the same code. The error is
Error using miqps_ot (line 346)
miqps_ot: EXITFLAG from price computation stage = 0
Error in miqps_matpower (line 220)
miqps_ot(H, c, A, l, u, xmin, xmax, x0, vtype, opt);
Error in most (line 2072)
mdo.QP.lambda ] = miqps_matpower( mdi.QP.H, mdi.QP.C, ...
Error in ammar1 (line 111)
mdo = most(mdi, mpopt);
Kind Regard,
Mohammad Ammar Raza
On 22 May 2017 at 15:04, Ammar Raza wrote:
> Hi Ray,
>
> Thank you for your email.
>
> Kind Regard,
>
> Mohammad Ammar Raza
>
>
>
>
> On 17 May 2017 at 14:51, Ray Zimmerman wrote:
>
>> Hi Mohammad,
>>
>> A unit commitment problem with quadratic generator costs or ramping costs
>> requires a solver that can handle Mixed Integer Quadratic Programming
>> (MIQP) problems. The only MATPOWER-compatible solvers that handle MIQP are
>> Gurobi, CPLEX and MOSEK, as stated in the System Requirements section of
>> the MOST User’s Manual.
>>
>> So, you will either need to install one of these solvers, or modify your
>> problem so that in does not include any quadratic costs.
>>
>>Ray
>>
>>
>> On May 16, 2017, at 5:30 PM, Ammar Raza wrote:
>>
>> Hi
>>
>> I am trying to do example 6 of the MOST user manual with Bus 118 and I
>> got this error.
>>
>> Error using mpopt2qpopt (line 92)
>> mpopt2qpopt: Sorry, no solver available for MIQP models
>>
>> Error in most (line 2065)
>> mdo.QP.opt = mpopt2qpopt(mpopt, model, 'most');
>>
>> Error in most_ammar6 (line 108)
>> mdo = most(mdi, mpopt);
>>
>>
>> Kind Regard,
>>
>> Mohammad Ammar Raza
>>
>>
>>
>>
>> On 15 May 2017 at 20:14, Ray Zimmerman wrote:
>>
>>> In a deterministic problem, transmat can simply be a scalar with the
>>> number of periods. See the description of transmat in the help for
>>> loadmd()
>>> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/most/loadmd.html>
>>> .
>>>
>>>Ray
>>>
>>>
>>>
>>> On May 15, 2017, at 7:37 AM, Ammar Raza wrote:
>>>
>>> Thank you Ray
>>>
>>> Therefore, Transmat should be like that?
>>>
>>> transmat = cell(1, nt);
>>> T = [ 1x 99];
>>> [transmat{:}] = deal(T * ones(1,1));
>>> transmat{1} = T;
>>>
>>> Kind Regard,
>>>
>>> Mohammad Ammar Raza
>>>
>>>
>>>
>>>
>>> On 11 May 2017 at 13:42, Ray Zimmerman wrote:
>>>
>>>> Just to clarify, if loadprofile.rows = 0, then loadprofile.values(:,
>>>> 1, :) = [24 x 1] and all loads are scaled together. Otherwise, the
>>>> number of columns in loadprofile.values(:, 1, :) must equal the length
>>>> of loadprofile.rows and each of the specified loads is scaled
>>>> separately.
>>>>
>>>>Ray
>>>>
>>>>
>>>> On May 9, 2017, at 10:29 AM, Ammar Raza wrote:
>>>>
>>>> Its means
>>>>
>>>> loadprofile.rows= [1 x 99 ]
>>>>
>>>> Thanks Stephanie for your help.
>>>>
>>>> Kind Regard,
>>>>
>>>> Mohammad Ammar Raza
>>>>
>>>>
>>>>
>>>>
>>>> On 9 May 2017 at 15:25, Stephanie wrote:
>>>>
>>>>> Yes, and the value of rows, 0 means the load is applied to the all
>>>>> buses, so you should also change the value of the rows by
>>>>> loadprofile.rows=load_index
>>>>> (e.g. loadprofile.rows=[2 3 5]), if there are three loads located in bus
>>>>> 2,
>>>>> 3 and 5, respectively)
>>>>>
>>>>> Best Regards,
>>>>> Stephanie
>>>>>
>>>>> 2017-05-09 22:09 GMT+08:00 Ammar Raza :
>>>>>
>>>>>> that means
>>>>>>
>>>>>> loadprofile = struct( ...
>>>>>> 'type', 'mpcData', ...
>>>>>> 'table', CT_TLOAD, ...
>>>>>> 'rows', 0, ...
>>>>>> 'col', CT_LOAD_ALL_PQ, ...
>>>>>> 'chgtype', CT_REP, ...
>>>>>> 'values', [] );
>>>>>>
>>>>>> loadprofile.values(:, 1, :) = [ 24 x 99 ]
>>>>>>
>>>>>>
>>>>>> Kind Regard,
>>>>>>
>>>>>> Mohammad Ammar Raza
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> On 9 May 2017 at 14:55, Stephanie wrote:
>>>>>>
>>>>>>> Hi,
>>>>>>>
>>>>>>> I tried it a few days before, what I did is:
>>>>>>>
>>>>>>> loadprofile.rows=load_index; %which applies the load to the bus it
>>>>>>> locates
>>>>>>> loadprofile.value (:,1,:) = [24 x 99];
>>>>>>>
>>>>>>> And it works in my situation
>>>>>>>
>>>>>>> Hope it will help you
>>>>>>> --
>>>>>>> Best Regards
>>>>>>> Stephanie
>>>>>>>
>>>>>>>
>>>>>>> *From:* Ammar Raza
>>>>>>> *Date:* 2017-05-09 21:43
>>>>>>> *To:* MATPOWER discussion forum
>>>>>>> *Subject:* MOST
>>>>>>> Hi all,
>>>>>>>
>>>>>>> I am working on bus 118 with 54 generator and 99 loads.
>>>>>>>
>>>>>>> I would like to ask you for the multi period problem i.e. 24 hours,
>>>>>>> how I will make the load profile function. Can I do like that
>>>>>>>
>>>>>>> loadprofile.value (:,1,1) = [24 x 99];
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> As in case3a in the MOST manual there is only one load therefore,
>>>>>>> load profile matrix is [12 x 1]
>>>>>>>
>>>>>>> loadprofile.value(:,1,1)=[
>>>>>>> 440;
>>>>>>> 480;
>>>>>>> ..
>>>>>>> ..
>>>>>>> .. ];
>>>>>>>
>>>>>>>
>>>>>>> Kind Regard,
>>>>>>>
>>>>>>> Mohammad Ammar Raza
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>
