loss2bus

LOSS2BUS   Accumulates branch real power losses at downstream buses.

   BUS_LOSS = LOSS2BUS(MPC)

   Takes a solved AC power flow case as input and returns an
   NB x 1 vector containing all branch active power losses
   accumulated to the bus at the downstream end of the branch.

0001 function bus_loss = loss2bus(mpc)
0002 %LOSS2BUS   Accumulates branch real power losses at downstream buses.
0003 %
0004 %   BUS_LOSS = LOSS2BUS(MPC)
0005 %
0006 %   Takes a solved AC power flow case as input and returns an
0007 %   NB x 1 vector containing all branch active power losses
0008 %   accumulated to the bus at the downstream end of the branch.
0009 
0010 %   MATPOWER
0011 %   Copyright (c) 2014-2016, Power Systems Engineering Research Center (PSERC)
0012 %   by Ray Zimmerman, PSERC Cornell
0013 %
0014 %   This file is part of MATPOWER Extras.
0015 %   Covered by the 3-clause BSD License (see LICENSE file for details).
0016 %   See https://github.com/MATPOWER/matpower-extras for more info.
0017 
0018 define_constants;
0019     
0020 %% create external to internal bus map
0021 nb  = size(mpc.bus, 1);     %% number of buses
0022 nl  = size(mpc.branch, 1);  %% number of branches
0023 mb  = max(abs([mpc.bus(:, BUS_I); mpc.gen(:, GEN_BUS); ...
0024             mpc.branch(:, F_BUS); mpc.branch(:, T_BUS)]));
0025 e2i = sparse(mpc.bus(:, BUS_I), ones(nb, 1), 1:nb, mb, 1);
0026 
0027 %% assign losses to downstream buses
0028 loss = mpc.branch(:, PF) + mpc.branch(:, PT);
0029 bus_loss = zeros(nb, 1);
0030 for j = 1:nl    %% need to use loop to accumulate for multiple lines per bus
0031     if mpc.branch(j, PF) >= mpc.branch(j, PT)
0032         b = T_BUS;
0033     else
0034         b = F_BUS;
0035     end
0036     ib = e2i(mpc.branch(j, b));
0037     bus_loss(ib) = bus_loss(ib) + loss(j);
0038 end