MAKEAPQ Construct linear constraints for generator capability curves. [APQH, UBPQH, APQL, UBPQL, DATA] = MAKEAPQ(BASEMVA, GEN) Constructs the parameters for the following linear constraints implementing trapezoidal generator capability curves, where Pg and Qg are the real and reactive generator injections. APQH * [Pg; Qg] <= UBPQH APQL * [Pg; Qg] <= UBPQL DATA constains additional information as shown below. Example: [Apqh, ubpqh, Apql, ubpql, data] = makeApq(baseMVA, gen); data.h [QC1MAX-QC2MAX, PC2-PC1] data.l [QC2MIN-QC1MIN, PC1-PC2] data.ipqh indices of gens with general PQ cap curves (upper) data.ipql indices of gens with general PQ cap curves (lower)
0001 function [Apqh, ubpqh, Apql, ubpql, data] = makeApq(baseMVA, gen) 0002 %MAKEAPQ Construct linear constraints for generator capability curves. 0003 % [APQH, UBPQH, APQL, UBPQL, DATA] = MAKEAPQ(BASEMVA, GEN) 0004 % 0005 % Constructs the parameters for the following linear constraints 0006 % implementing trapezoidal generator capability curves, where 0007 % Pg and Qg are the real and reactive generator injections. 0008 % 0009 % APQH * [Pg; Qg] <= UBPQH 0010 % APQL * [Pg; Qg] <= UBPQL 0011 % 0012 % DATA constains additional information as shown below. 0013 % 0014 % Example: 0015 % [Apqh, ubpqh, Apql, ubpql, data] = makeApq(baseMVA, gen); 0016 % 0017 % data.h [QC1MAX-QC2MAX, PC2-PC1] 0018 % data.l [QC2MIN-QC1MIN, PC1-PC2] 0019 % data.ipqh indices of gens with general PQ cap curves (upper) 0020 % data.ipql indices of gens with general PQ cap curves (lower) 0021 0022 % MATPOWER 0023 % $Id: makeApq.m 1635 2010-04-26 19:45:26Z ray $ 0024 % by Ray Zimmerman, PSERC Cornell 0025 % and Carlos E. Murillo-Sanchez, PSERC Cornell & Universidad Autonoma de Manizales 0026 % Copyright (c) 1996-2010 by Power System Engineering Research Center (PSERC) 0027 % 0028 % This file is part of MATPOWER. 0029 % See http://www.pserc.cornell.edu/matpower/ for more info. 0030 % 0031 % MATPOWER is free software: you can redistribute it and/or modify 0032 % it under the terms of the GNU General Public License as published 0033 % by the Free Software Foundation, either version 3 of the License, 0034 % or (at your option) any later version. 0035 % 0036 % MATPOWER is distributed in the hope that it will be useful, 0037 % but WITHOUT ANY WARRANTY; without even the implied warranty of 0038 % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 0039 % GNU General Public License for more details. 0040 % 0041 % You should have received a copy of the GNU General Public License 0042 % along with MATPOWER. If not, see <http://www.gnu.org/licenses/>. 0043 % 0044 % Additional permission under GNU GPL version 3 section 7 0045 % 0046 % If you modify MATPOWER, or any covered work, to interface with 0047 % other modules (such as MATLAB code and MEX-files) available in a 0048 % MATLAB(R) or comparable environment containing parts covered 0049 % under other licensing terms, the licensors of MATPOWER grant 0050 % you additional permission to convey the resulting work. 0051 0052 %% define named indices into data matrices 0053 [GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ... 0054 MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ... 0055 QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen; 0056 0057 %% data dimensions 0058 ng = size(gen, 1); %% number of dispatchable injections 0059 0060 %% which generators require additional linear constraints 0061 %% (in addition to simple box constraints) on (Pg,Qg) to correctly 0062 %% model their PQ capability curves 0063 ipqh = find( hasPQcap(gen, 'U') ); 0064 ipql = find( hasPQcap(gen, 'L') ); 0065 npqh = size(ipqh, 1); %% number of general PQ capability curves (upper) 0066 npql = size(ipql, 1); %% number of general PQ capability curves (lower) 0067 0068 %% make Apqh if there is a need to add general PQ capability curves; 0069 %% use normalized coefficient rows so multipliers have right scaling 0070 %% in $$/pu 0071 if npqh > 0 0072 data.h = [gen(ipqh,QC1MAX)-gen(ipqh,QC2MAX), gen(ipqh,PC2)-gen(ipqh,PC1)]; 0073 ubpqh = data.h(:, 1) .* gen(ipqh,PC1) + data.h(:, 2) .* gen(ipqh,QC1MAX); 0074 for i=1:npqh, 0075 tmp = norm(data.h(i,:)); 0076 data.h(i,:) = data.h(i, :) / tmp; 0077 ubpqh(i) = ubpqh(i) / tmp; 0078 end 0079 Apqh = sparse([1:npqh, 1:npqh]', [ipqh; ipqh+ng], ... 0080 data.h(:), npqh, 2*ng); 0081 ubpqh = ubpqh / baseMVA; 0082 else 0083 data.h = []; 0084 Apqh = sparse(0, 2*ng); 0085 ubpqh = []; 0086 end 0087 0088 %% similarly Apql 0089 if npql > 0 0090 data.l = [gen(ipql,QC2MIN)-gen(ipql,QC1MIN), gen(ipql,PC1)-gen(ipql,PC2)]; 0091 ubpql= data.l(:, 1) .* gen(ipql,PC1) + data.l(:, 2) .* gen(ipql,QC1MIN) ; 0092 for i=1:npql, 0093 tmp = norm(data.l(i, : )); 0094 data.l(i, :) = data.l(i, :) / tmp; 0095 ubpql(i) = ubpql(i) / tmp; 0096 end 0097 Apql = sparse([1:npql, 1:npql]', [ipql; ipql+ng], ... 0098 data.l(:), npql, 2*ng); 0099 ubpql = ubpql / baseMVA; 0100 else 0101 data.l = []; 0102 Apql = sparse(0, 2*ng); 0103 ubpql = []; 0104 end 0105 0106 data.ipql = ipql; 0107 data.ipqh = ipqh;