pqdata.m

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

function data = pqdata(a,idx,qlim,sp,lambda,one)

if ~a.n, return, end

idx = idx(1);
data = [a.con(idx,[1 2 3 10 6 8 9]), 0];
data(1,4) = -a.pg(idx);  

switch qlim
 case 'qmax'
  data(1,5) = -a.con(idx,6);  
  fm_disp([sp,'SW gen. at bus #', fvar(a.bus(idx),4), ...
           ' reached Q_max at lambda = ',fvar(lambda-one,9)])
 case 'qmin'
  data(1,5) = -a.con(idx,7);  
  fm_disp([sp,'SW gen. at bus #', fvar(a.bus(Qmin_idx(1)),4), ...
           ' reached Q_min at lambda = ',fvar(lambda-one,9)])
end