suqlim.m

来自「基于PSAT 软件的多目标最优潮流计算用于中小型电力系统的分析和管理」· M 代码 · 共 27 行

function [qx,qn] = suqlim(a,Qmax,Qmin,bus)

if ~a.n
  qx = [];
  qn = [];
  return
end

qx = zeros(a.n, 1);
qn = zeros(a.n, 1);

for i = 1:length(bus)
  idx = find(a.bus == bus(i));
  if isempty(idx), continue, end
  qx(idx) = a.con(idx,16);
  qn(idx) = a.con(idx,17);
  maxq = sum(abs(a.con(idx,16)));
  minq = sum(abs(a.con(idx,17)));
  if maxq == 0
    qx(idx) = Qmax(i)/length(idx); 
  end
  if minq == 0
    qn(idx) = Qmin(i)/length(idx); 
  end
end