Home > matpower7.0 > lib > runopf_w_res.m

runopf_w_res

PURPOSE ^

RUNOPF_W_RES Runs an optimal power flow with fixed zonal reserves.

SYNOPSIS ^

function [varargout] = runopf_w_res(varargin)

DESCRIPTION ^

RUNOPF_W_RES  Runs an optimal power flow with fixed zonal reserves.
   RESULTS = RUNOPF_W_RES(CASEDATA, MPOPT, FNAME, SOLVEDCASE)
   [RESULTS, SUCCESS] = RUNOPF_W_RES(CASEDATA, MPOPT, FNAME, SOLVEDCASE)

   Runs an optimal power flow with the addition of reserve requirements
   specified as a set of fixed zonal reserves. See RUNOPF for a
   description of the input and output arguments, which are the same,
   with the exception that the case file or struct CASEDATA must define
   a 'reserves' field, which is a struct with the following fields:
       zones   nrz x ng, zone(i, j) = 1, if gen j belongs to zone i
                                      0, otherwise
       req     nrz x 1, zonal reserve requirement in MW
       cost    (ng or ngr) x 1, cost of reserves in $/MW
       qty     (ng or ngr) x 1, max quantity of reserves in MW (optional)
   where nrz is the number of reserve zones and ngr is the number of
   generators belonging to at least one reserve zone and ng is the total
   number of generators.

   In addition to the normal OPF output, the RESULTS struct contains a
   new 'reserves' field with the following fields, in addition to those
   provided in the input:
       R       - ng x 1, reserves provided by each gen in MW
       Rmin    - ng x 1, lower limit on reserves provided by each gen, (MW)
       Rmax    - ng x 1, upper limit on reserves provided by each gen, (MW)
       mu.l    - ng x 1, shadow price on reserve lower limit, ($/MW)
       mu.u    - ng x 1, shadow price on reserve upper limit, ($/MW)
       mu.Pmax - ng x 1, shadow price on Pg + R <= Pmax constraint, ($/MW)
       prc     - ng x 1, reserve price for each gen equal to maximum of the
                         shadow prices on the zonal requirement constraint
                         for each zone the generator belongs to

   See T_CASE30_USERFCNS for an example case file with fixed reserves,
   and TOGGLE_RESERVES for the implementation.

   Calling syntax options:
       results = runopf_w_res(casedata);
       results = runopf_w_res(casedata, mpopt);
       results = runopf_w_res(casedata, mpopt, fname);
       results = runopf_w_res(casedata, mpopt, fname, solvedcase);
       [results, success] = runopf_w_res(...);

   Example:
       results = runopf_w_res('t_case30_userfcns');

   See also RUNOPF, TOGGLE_RESERVES, T_CASE30_USERFCNS.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [varargout] = runopf_w_res(varargin)
0002 %RUNOPF_W_RES  Runs an optimal power flow with fixed zonal reserves.
0003 %   RESULTS = RUNOPF_W_RES(CASEDATA, MPOPT, FNAME, SOLVEDCASE)
0004 %   [RESULTS, SUCCESS] = RUNOPF_W_RES(CASEDATA, MPOPT, FNAME, SOLVEDCASE)
0005 %
0006 %   Runs an optimal power flow with the addition of reserve requirements
0007 %   specified as a set of fixed zonal reserves. See RUNOPF for a
0008 %   description of the input and output arguments, which are the same,
0009 %   with the exception that the case file or struct CASEDATA must define
0010 %   a 'reserves' field, which is a struct with the following fields:
0011 %       zones   nrz x ng, zone(i, j) = 1, if gen j belongs to zone i
0012 %                                      0, otherwise
0013 %       req     nrz x 1, zonal reserve requirement in MW
0014 %       cost    (ng or ngr) x 1, cost of reserves in $/MW
0015 %       qty     (ng or ngr) x 1, max quantity of reserves in MW (optional)
0016 %   where nrz is the number of reserve zones and ngr is the number of
0017 %   generators belonging to at least one reserve zone and ng is the total
0018 %   number of generators.
0019 %
0020 %   In addition to the normal OPF output, the RESULTS struct contains a
0021 %   new 'reserves' field with the following fields, in addition to those
0022 %   provided in the input:
0023 %       R       - ng x 1, reserves provided by each gen in MW
0024 %       Rmin    - ng x 1, lower limit on reserves provided by each gen, (MW)
0025 %       Rmax    - ng x 1, upper limit on reserves provided by each gen, (MW)
0026 %       mu.l    - ng x 1, shadow price on reserve lower limit, ($/MW)
0027 %       mu.u    - ng x 1, shadow price on reserve upper limit, ($/MW)
0028 %       mu.Pmax - ng x 1, shadow price on Pg + R <= Pmax constraint, ($/MW)
0029 %       prc     - ng x 1, reserve price for each gen equal to maximum of the
0030 %                         shadow prices on the zonal requirement constraint
0031 %                         for each zone the generator belongs to
0032 %
0033 %   See T_CASE30_USERFCNS for an example case file with fixed reserves,
0034 %   and TOGGLE_RESERVES for the implementation.
0035 %
0036 %   Calling syntax options:
0037 %       results = runopf_w_res(casedata);
0038 %       results = runopf_w_res(casedata, mpopt);
0039 %       results = runopf_w_res(casedata, mpopt, fname);
0040 %       results = runopf_w_res(casedata, mpopt, fname, solvedcase);
0041 %       [results, success] = runopf_w_res(...);
0042 %
0043 %   Example:
0044 %       results = runopf_w_res('t_case30_userfcns');
0045 %
0046 %   See also RUNOPF, TOGGLE_RESERVES, T_CASE30_USERFCNS.
0047 
0048 %   MATPOWER
0049 %   Copyright (c) 2008-2016, Power Systems Engineering Research Center (PSERC)
0050 %   by Ray Zimmerman, PSERC Cornell
0051 %
0052 %   This file is part of MATPOWER.
0053 %   Covered by the 3-clause BSD License (see LICENSE file for details).
0054 %   See https://matpower.org for more info.
0055 
0056 mpc = loadcase(varargin{1});
0057 mpc = toggle_reserves(mpc, 'on');
0058 [varargout{1:nargout}] = runopf(mpc, varargin{2:nargin});
0059 
0060 if nargout > 0 && isstruct(varargout{1})
0061     varargout{1} = toggle_reserves(varargout{1}, 'off');
0062 end

Generated on Mon 24-Jun-2019 15:58:45 by m2html © 2005