pqlim.m

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

function [a, check] = pqlim(a,vlim,sp,lambda,one)  

global DAE

check = 0;

if ~vlim, return, end
if ~a.n, return, end

vmax_idx = find(DAE.V(a.bus) > a.con(:,6) & a.vmax);
vmin_idx = find(DAE.V(a.bus) < a.con(:,7) & a.vmin);
if ~isempty(vmin_idx)
  a.vmin(vmin_idx(1)) = 0;
  fm_disp([sp,'PQ load at bus #',fvar(a.bus(vmin_idx(1)),4), ...
           ' reached V_min at lambda = ',fvar(lambda-one,9)])
  check = 1;
  return
end
if ~isempty(vmax_idx)
  a.vmax(vmax_idx(1)) = 0;
  fm_disp([sp,'PQ load at bus #',fvar(a.bus(vmax_idx(1)),4), ...
           ' reached V_max at lambda = ',fvar(lambda-one,9)])
  check = 1;
end