fm_genstatus.m

来自「用于电力系统的一个很好的分析软件」· M 代码 · 共 48 行

function  u = fm_genstatus(idx)
% FM_GENSTATUS finds and remove static generators
%
% CHECK = FM_GENSTATUS(IDX)
%       IDX   = bus index where to look for generators
%       U = 0 -> generator off-line (or no generator found)
%       U = 1 -> generator on-line
%
%Author:    Federico Milano
%Date:      24-Aug-2007
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://www.uclm.es/area/gsee/Web/Federico
%
% Copyright (C) 2002-2008 Federico Milano

global SW PV PQ

u = zeros(length(idx),1);

for i = 1:length(idx)
  k = idx(i);
  if k <= 0, continue, end
  idx_sw = findbus(SW,k);
  idx_pv = findbus(PV,k);
  idx_pq = findgen(PQ,k);

  if isempty(idx_sw)
    u_sw = 0;
  else
    u_sw = SW.u(idx_sw);
  end

  if isempty(idx_pv)
    u_pv = 0;
  else
    u_pv = PV.u(idx_pv);
  end

  if isempty(idx_pq)
    u_pq = 0;
  else
    u_pq = PQ.u(idx_pq);
  end

  u(i) = u_sw | u_pv | u_pq;
end