fm_qlim.m

这是一个很适合研究和学习用的电力系统仿真软件

function  [Qx,Qn] = fm_qlim(maxQ,minQ,flag)
% FM_VLIM determines max and min bus reactive powers
%
% [QMAX,QMIN] = FM_QLIM(MAXQ,MINQ,FLAG)
%       MAXQ -> default max reactive power (p.u.)
%       MINQ -> default min reactive power (p.u.)
%       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://thunderbox.uwaterloo.ca/~fmilano
%
%Copyright (C) 2002-2006 Federico Milano

global Bus SW PV Supply Settings

switch flag
 
 case 'all'
  Qn = zeros(Bus.n,1);
  Qx = zeros(Bus.n,1);
  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
  idx = find(Qn == 0 & Qx == 0);
  if ~isempty(idx)
    Qn(idx) = minQ*Settings.mva;
    Qx(idx) = maxQ*Settings.mva;
  end
 
 case 'gen'
  busg = [SW.bus; PV.bus];
  Qmin = [qmin(SW); qmin(PV)];
  Qmax = [qmax(SW); qmax(PV)];
  idx = find(Qmin == 0 & Qmax == 0);
  if ~isempty(idx)
    Qmin(idx) = minQ*Settings.mva;
    Qmax(idx) = maxQ*Settings.mva;
  end
  [Qx,Qn] = suqlim(Supply,Qmax,Qmin,busg);
  [busS,idxS] = setdiff(busg,Supply.bus);
  if ~isempty(busS)
    Qn = [Qn; Qmin(idxS)];
    Qx = [Qx; Qmax(idxS)];
  end

end