dynidx.m

来自「用于电力系统的一个很好的分析软件」· M 代码 · 共 14 行

function a = dynidx(a)

global DAE

if ~a.n, return, end

a.xp = zeros(a.n,1);
a.xq = zeros(a.n,1);
for i = 1:a.n
  a.xp(i) = DAE.n + 1;
  a.xq(i) = DAE.n + 2;
  DAE.n = DAE.n + 2;
end