fm_busidx.m

电力系统的psat

function  b = fm_busidx(bus_no)
% FM_BUSIDX find input and output line at a certain bus
%
% IDX = FM_BUSIDX(BUS_NO)
%       BUS_NO bus number
%       IDX = [line_from_#, line_to_#]
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://thunderbox.uwaterloo.ca/~fmilano
%
% Copyright (C) 2002-2005 Federico Milano
%
% This toolbox is free software; you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation; either version 2.0 of the License, or
% (at your option) any later version.
%
% This toolbox is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANDABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
% General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this toolbox; if not, write to the Free Software
% Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307,
% USA.

global Line

% indice delle correnti misurate:  I_idx => [#corrente, #linea, from bus, to bus, sign]

I_idx_from = find(Line.con(:,1) == bus_no);
I_idx_to = find(Line.con(:,2) == bus_no);
n_current = length(I_idx_from) + length(I_idx_to);
b = [Line.con(I_idx_from,2); Line.con(I_idx_to,1)];