bsvc.m

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

function B = bsvc(a)

global DAE

B = zeros(a.n,1); 

if a.ty1
  B(a.ty1) = DAE.x(a.bcv);
end

if a.ty2
  xl = a.con(a.ty2,15);
  xc = a.con(a.ty2,16);
  B(a.ty2) = (2*DAE.x(a.alpha) - sin(2*DAE.x(a.alpha)) ...
              - pi*(2-xl./xc))./(pi*xl);
end