dynidx.m

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

function a = dynidx(a)

global DAE

if ~a.n, return, end

a.x1 = zeros(a.n,1);
a.x2 = zeros(a.n,1);

for i = 1:a.n
  a.x1(i) = DAE.n + 1;
  if a.con(i,3) == 2
    a.x2(i) = DAE.n + 2;
    DAE.n = DAE.n + 2;
  else
    DAE.n = DAE.n + 1;
  end
end