fm_pmurec.m

电力系统的psat

function [bus_pmu, bus2] = fm_pmurec(bus1,bus2,bus0)
% FM_PMUREC routine for PMU placement
%
% (...) = FM_PMUREC(...)
%
% This routine is called by FM_PMULOC
%
%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 Bus

%bus2(Bus.int(bus1)) = 1;

b = fm_busidx(bus1);
if length(b) == 1
  bus_pmu = b;
  bus2(Bus.int(b)) = 1;
  return
end

b1 = find(~bus2(Bus.int(b)));
if isempty(b1)
  bus_pmu = [];
  return
else
  b = sort(b(b1));
end

bus2(Bus.int(b)) = ones(length(b),1);

bb = [];

for k = 1:length(b)
  abb = fm_busidx(b(k));
  abb = abb(find(~bus2(Bus.int(abb))));
  bb = [bb; abb];
end

if isempty(bb); bus_pmu = []; return; end

bb = sort(bb);
num_bus = length(bb);
nodi_el = [];
for jjj = 1:num_bus-1
  nodi_el = [nodi_el; jjj + find(bb([jjj+1:num_bus]) == bb(jjj))];
end
bb(nodi_el) = [];
bb = bb(find(~bus2(Bus.int(bb))));
bus2(Bus.int(bb)) = ones(length(bb),1);

bus_pmu = sort(bb);