fm_qlim.m

电力系统分析计算程序

function  [Qx,Qn] = fm_qlim(flag)
% FM_VLIM determines max and min bus reactive powers
%
% [QMAX,QMIN] = FM_QLIM(FLAG)
%       FLAG -> 'all' full bus vector
%               'gen' vector of generator buses
%       QMAX -> vector of max bus reactive powers (p.u.)
%       QMIN -> vector of min bus reactive powers (p.u.)
%
%Author:    Federico Milano
%Date:      27-Dic-2005
%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 Bus SW PV Supply Settings

switch flag

 case 'all'

  Qn = getzeros(Bus);
  Qx = getzeros(Bus);
  if PV.n
    Qn(PV.bus) = qmin(PV);
    Qx(PV.bus) = qmax(PV);
  end
  if SW.n
    Qn(SW.bus) = qmin(SW);
    Qx(SW.bus) = qmax(SW);
  end

 case 'gen'

  busg = [SW.bus; PV.bus];
  Qmin = [qmin(SW); qmin(PV)];
  Qmax = [qmax(SW); qmax(PV)];
  [Qx,Qn] = suqlim(Supply,Qmax,Qmin,busg);
  if Settings.octave
    busS = setdiff(busg,Supply.bus);
    [dummy,idxS] = ismember(busS,Supply.bus);
  else    
    [busS,idxS] = setdiff(busg,Supply.bus);
  end
  if ~isempty(busS)
    Qn = [Qn; Qmin(idxS)];
    Qx = [Qx; Qmax(idxS)];
  end

end