Home > matpower7.0 > lib > update_mupq.m

update_mupq

PURPOSE ^

UPDATE_MUPQ Updates values of generator limit shadow prices.

SYNOPSIS ^

function gen = update_mupq(baseMVA, gen, mu_PQh, mu_PQl, data)

DESCRIPTION ^

UPDATE_MUPQ  Updates values of generator limit shadow prices.
   GEN = UPDATE_MUPQ(BASEMVA, GEN, MU_PQH, MU_PQL, DATA)

   Updates the values of MU_PMIN, MU_PMAX, MU_QMIN, MU_QMAX based
   on any shadow prices on the sloped portions of the generator
   capability curve constraints.

   MU_PQH - shadow prices on upper sloped portion of capability curves
   MU_PQL - shadow prices on lower sloped portion of capability curves
   DATA   - "data" struct returned by MAKEAPQ

   See also MAKEAPQ.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function gen = update_mupq(baseMVA, gen, mu_PQh, mu_PQl, data)
0002 %UPDATE_MUPQ  Updates values of generator limit shadow prices.
0003 %   GEN = UPDATE_MUPQ(BASEMVA, GEN, MU_PQH, MU_PQL, DATA)
0004 %
0005 %   Updates the values of MU_PMIN, MU_PMAX, MU_QMIN, MU_QMAX based
0006 %   on any shadow prices on the sloped portions of the generator
0007 %   capability curve constraints.
0008 %
0009 %   MU_PQH - shadow prices on upper sloped portion of capability curves
0010 %   MU_PQL - shadow prices on lower sloped portion of capability curves
0011 %   DATA   - "data" struct returned by MAKEAPQ
0012 %
0013 %   See also MAKEAPQ.
0014 
0015 %   MATPOWER
0016 %   Copyright (c) 1996-2016, Power Systems Engineering Research Center (PSERC)
0017 %   by Ray Zimmerman, PSERC Cornell
0018 %   and Carlos E. Murillo-Sanchez, PSERC Cornell & Universidad Nacional de Colombia
0019 %
0020 %   This file is part of MATPOWER.
0021 %   Covered by the 3-clause BSD License (see LICENSE file for details).
0022 %   See https://matpower.org for more info.
0023 
0024 %% extract the constraint parameters
0025 [ipqh, ipql, Apqhdata, Apqldata] = ...
0026     deal(data.ipqh, data.ipql, data.h, data.l);
0027 
0028 %% define named indices into data matrices
0029 [GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
0030     MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ...
0031     QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
0032 
0033 %% combine original limit multipliers into single value
0034 muP = gen(:, MU_PMAX) - gen(:, MU_PMIN);
0035 muQ = gen(:, MU_QMAX) - gen(:, MU_QMIN);
0036 
0037 %% add P and Q components of multipliers on upper sloped constraint
0038 if ~isempty(ipqh)
0039     muP(ipqh) = muP(ipqh) - mu_PQh .* Apqhdata(:,1)/baseMVA;
0040     muQ(ipqh) = muQ(ipqh) - mu_PQh .* Apqhdata(:,2)/baseMVA;
0041 end
0042 
0043 %% add P and Q components of multipliers on lower sloped constraint
0044 if ~isempty(ipql)
0045     muP(ipql) = muP(ipql) - mu_PQl .* Apqldata(:,1)/baseMVA;
0046     muQ(ipql) = muQ(ipql) - mu_PQl .* Apqldata(:,2)/baseMVA;
0047 end
0048 
0049 %% split back into upper and lower multipliers based on sign
0050 gen(:, MU_PMAX) = (muP > 0) .*  muP;
0051 gen(:, MU_PMIN) = (muP < 0) .* -muP;
0052 gen(:, MU_QMAX) = (muQ > 0) .*  muQ;
0053 gen(:, MU_QMIN) = (muQ < 0) .* -muQ;

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