addgen2mpc

ADDGEN2MPC Adds a set of generators to a MATPOWER case struct.

   [NEW_MPC, IDX] = ADDGEN2MPC(MPC, GEN, GENCOST, GEN_TYPE)

   Inserts a set of generators to a MATPOWER case struct by adding the
   info of the new generators (in inputs GEN, GENCOST and GEN_TYPE) to the
   bottom of the corresponding tables (fields of MPC). Dimensions must
   be consistent.

   Inputs:
       MPC : standard MATPOWER case struct, with the following additional
             fields:
               .genfuel : (optional) cell array of strings with fuel type
                          (filled with 'unknown' if missing)
               .i<type> : vector of generator indices for units of the
                          specified fuel type
       GEN: standard MATPOWER generator matrix for units to be added
       GENCOST: standard MATPOWER generator cost matrix for units to be added
       GEN_TYPE: string or cell array of names of fuel or type of new gens

   Outputs:
       NEW_MPC : the new MATPOWER case with the additional generators
                 appended to GEN, GENCOST, GENFUEL and additional field
                   .i<type> : vector of generator indices for units of the
                              specified fuel type
       IDX : generator indices of the newly added generators

0001 function [mpco, NewGenIdx] = addgen2mpc(mpci, gen, gencost, fuel)
0002 %ADDGEN2MPC Adds a set of generators to a MATPOWER case struct.
0003 %
0004 %   [NEW_MPC, IDX] = ADDGEN2MPC(MPC, GEN, GENCOST, GEN_TYPE)
0005 %
0006 %   Inserts a set of generators to a MATPOWER case struct by adding the
0007 %   info of the new generators (in inputs GEN, GENCOST and GEN_TYPE) to the
0008 %   bottom of the corresponding tables (fields of MPC). Dimensions must
0009 %   be consistent.
0010 %
0011 %   Inputs:
0012 %       MPC : standard MATPOWER case struct, with the following additional
0013 %             fields:
0014 %               .genfuel : (optional) cell array of strings with fuel type
0015 %                          (filled with 'unknown' if missing)
0016 %               .i<type> : vector of generator indices for units of the
0017 %                          specified fuel type
0018 %       GEN: standard MATPOWER generator matrix for units to be added
0019 %       GENCOST: standard MATPOWER generator cost matrix for units to be added
0020 %       GEN_TYPE: string or cell array of names of fuel or type of new gens
0021 %
0022 %   Outputs:
0023 %       NEW_MPC : the new MATPOWER case with the additional generators
0024 %                 appended to GEN, GENCOST, GENFUEL and additional field
0025 %                   .i<type> : vector of generator indices for units of the
0026 %                              specified fuel type
0027 %       IDX : generator indices of the newly added generators
0028 
0029 %   MOST
0030 %   Copyright (c) 2013-2016, Power Systems Engineering Research Center (PSERC)
0031 %   by Daniel Munoz-Alvarez and Ray Zimmerman, PSERC Cornell
0032 %
0033 %   This file is part of MOST.
0034 %   Covered by the 3-clause BSD License (see LICENSE file for details).
0035 %   See http://www.pserc.cornell.edu/matpower/ for more info.
0036 
0037 [ng, ncg]     = size(mpci.gen);
0038 [nr, nc]      = size(gen);
0039 [ng2, ncgc]   = size(mpci.gencost);
0040 [nr2, ncgc2]  = size(gencost);
0041 if nr ~= nr2
0042     error('addgen2mpc: number of rows in GEN (%d) and GENCOST (%d) must agree', nr, nr2);
0043 end
0044 if ng ~= ng2
0045     error('addgen2mpc: number of rows in MPC.GEN (%d) and MPC.GENCOST (%d) must agree', ng, ng2);
0046 end
0047 
0048 %% append rows to GEN
0049 mpco = mpci;
0050 if nc <= ncg
0051     mpco.gen = [mpco.gen ; gen zeros(nr, ncg-nc)];
0052 else
0053     error('addgen2mpc: number of columns of GEN (%d) exceeds number of columns of MPC.GEN (%d)', nc, ncg);
0054 end
0055 
0056 %% append rows to GENCOST
0057 dim_gencost = size(gencost);
0058 dim_gencost_cur = size(mpco.gencost);
0059 if dim_gencost(2) < dim_gencost_cur(2)
0060     mpco.gencost = [mpco.gencost ; [gencost zeros(dim_gencost(1),dim_gencost_cur(2)-dim_gencost(2))]];
0061 elseif dim_gencost(2) > dim_gencost_cur(2)
0062     mpco.gencost = [mpco.gencost zeros(dim_gencost_cur(1),dim_gencost(2)-dim_gencost_cur(2))];
0063     mpco.gencost = [mpco.gencost ; gencost];
0064 else
0065     mpco.gencost = [mpco.gencost ; gencost];
0066 end
0067 
0068 %% initialize 'genfuel'  and 'i<type>' fields if missing
0069 if ~isfield(mpco, 'genfuel')
0070     mpco.genfuel = mat2cell(repmat('unknown', ng, 1), ones(ng,1), 7);
0071 end
0072 if ~isfield(mpco,['i' fuel])
0073     mpco.(['i' fuel]) = [];
0074 end
0075 
0076 %% append 'genfuel' and 'i<type>' fields
0077 if iscell(fuel)
0078     for i = dim_gencost(1)
0079         mpco.genfuel = [mpco.genfuel ; fuel{i}];
0080         mpco.(['i' fuel{i}]) = [ mpco.(['i' fuel{i}]) ; dim_gencost_cur(1) + i ];
0081     end
0082 elseif ischar(fuel)
0083     for i = 1:dim_gencost(1)
0084         mpco.genfuel = [ mpco.genfuel ; fuel ];
0085         mpco.(['i' fuel]) = [ mpco.(['i' fuel]) ; dim_gencost_cur(1) + i ];
0086     end
0087 else
0088     error('addgen2mpc: GEN_TYPE must be a string (or cell array of strings) indicating the fuel type of the new generators');
0089 end
0090 
0091 NewGenIdx = ( ng + 1 : size(mpco.gen,1) )';