I'm so sorry that the picture can't be seen in the system.
For Part One, the first objective funtion is min.
sum[(ai+bi*Pi+ci*Pi^2)+(di*ri)], and the second is
sum[(ai+bi*Pi+ci*Pi^2)+(αi+βi*Pi+γi*Pi^2+μi*exp(ηi*Pi))].
For Part Two, the objective function is
sum{sum{ti*[(ai+bi*Pi+ci*Pi^2)+(αi+βi*Pi+γi*Pi^2+μi*exp(ηi*Pi))]}}.
Thanks again.
At 2020-01-03 01:36:43, "Ray Daniel Zimmerman" <[email protected]> wrote:
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) 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]> 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.
