Home > matpower4.1 > bustypes.m

bustypes

PURPOSE ^

BUSTYPES Builds index lists for each type of bus (REF, PV, PQ).

SYNOPSIS ^

function [ref, pv, pq] = bustypes(bus, gen)

DESCRIPTION ^

BUSTYPES   Builds index lists for each type of bus (REF, PV, PQ).
   [REF, PV, PQ] = BUSTYPES(BUS, GEN)
   Generators with "out-of-service" status are treated as PQ buses with
   zero generation (regardless of Pg/Qg values in gen). Expects BUS and
   GEN have been converted to use internal consecutive bus numbering.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [ref, pv, pq] = bustypes(bus, gen)
0002 %BUSTYPES   Builds index lists for each type of bus (REF, PV, PQ).
0003 %   [REF, PV, PQ] = BUSTYPES(BUS, GEN)
0004 %   Generators with "out-of-service" status are treated as PQ buses with
0005 %   zero generation (regardless of Pg/Qg values in gen). Expects BUS and
0006 %   GEN have been converted to use internal consecutive bus numbering.
0007 
0008 %   MATPOWER
0009 %   $Id: bustypes.m,v 1.13 2010/04/26 19:45:25 ray Exp $
0010 %   by Ray Zimmerman, PSERC Cornell
0011 %   Copyright (c) 1996-2010 by Power System Engineering Research Center (PSERC)
0012 %
0013 %   This file is part of MATPOWER.
0014 %   See http://www.pserc.cornell.edu/matpower/ for more info.
0015 %
0016 %   MATPOWER is free software: you can redistribute it and/or modify
0017 %   it under the terms of the GNU General Public License as published
0018 %   by the Free Software Foundation, either version 3 of the License,
0019 %   or (at your option) any later version.
0020 %
0021 %   MATPOWER is distributed in the hope that it will be useful,
0022 %   but WITHOUT ANY WARRANTY; without even the implied warranty of
0023 %   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
0024 %   GNU General Public License for more details.
0025 %
0026 %   You should have received a copy of the GNU General Public License
0027 %   along with MATPOWER. If not, see <http://www.gnu.org/licenses/>.
0028 %
0029 %   Additional permission under GNU GPL version 3 section 7
0030 %
0031 %   If you modify MATPOWER, or any covered work, to interface with
0032 %   other modules (such as MATLAB code and MEX-files) available in a
0033 %   MATLAB(R) or comparable environment containing parts covered
0034 %   under other licensing terms, the licensors of MATPOWER grant
0035 %   you additional permission to convey the resulting work.
0036 
0037 %% constants
0038 [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
0039     VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
0040 [GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
0041     MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ...
0042     QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
0043 
0044 %% get generator status
0045 % bus_gen_status = zeros(size(bus, 1), 1);
0046 % bus_gen_status(gen(:, GEN_BUS)) = gen(:, GEN_STATUS) > 0;
0047 nb = size(bus, 1);
0048 ng = size(gen, 1);
0049 Cg = sparse(gen(:, GEN_BUS), (1:ng)', gen(:, GEN_STATUS) > 0, nb, ng);  %% gen connection matrix
0050                                         %% element i, j is 1 if, generator j at bus i is ON
0051 bus_gen_status = Cg * ones(ng, 1);      %% number of generators at each bus that are ON
0052 
0053 
0054 %% form index lists for slack, PV, and PQ buses
0055 ref = find(bus(:, BUS_TYPE) == REF & bus_gen_status);   %% reference bus index
0056 pv  = find(bus(:, BUS_TYPE) == PV  & bus_gen_status);   %% PV bus indices
0057 pq  = find(bus(:, BUS_TYPE) == PQ | ~bus_gen_status);   %% PQ bus indices
0058 
0059 %% pick a new reference bus if for some reason there is none (may have been shut down)
0060 if isempty(ref)
0061     ref = pv(1);                %% use the first PV bus
0062     pv = pv(2:length(pv));      %% take it off PV list
0063 end

Generated on Mon 26-Jan-2015 15:00:13 by m2html © 2005