Hi Yang, It’s not clear whether you simply want to use the fixed reserves extension (see Section 7.6.1 in the MATPOWER User’s Manual<https://matpower.org/docs/MATPOWER-manual-7.0.pdf>) or you want to implement your own OPF extension. In the first case you can simply use runopf_w_res(). In the second, please see Chapter 7 and the code for the various extensions described in Section 7.6.
Ray On Dec 25, 2019, at 8:11 PM, Yang <[email protected]<mailto:[email protected]>> wrote: Dear sir, I want to do the extensible OPF, and I read paper MATPOWER’s Extensible Optimal Power Flow Architecture. The objective function is min f(x)+fu(x,z),and fu=sum(ci*ri), following the code: Ar = [I I]; om = add_vars(om, 'R', ng, [], Rmin, Rmax); om = add_constraints(om, 'Pg_plus_R', Ar, [], Pmax, {'Pg', 'R'}); om = add_constraints(om, 'Rreq', Az, Rreq, [], {'R'}); om = add_costs(om, 'Rcost', struct('N',I,'Cw',Rcost), {'R'}); I add these into .m file, as for add_vars if I should change [], Rmin, Rmax into data I need, and do the same with other codes. I don't know if it's right. Best regards. Yang.
