dynidx.m

来自「这是一个很适合研究和学习用的电力系统仿真软件」· M 代码 · 共 19 行

function a = dynidx(a)

global DAE

if ~a.n, return, end

a.vcs = zeros(a.n,1);  
for i = 1:a.n
  if a.con(i,2) == 3
    a.vcs(i,1) = DAE.n + 1;
    a.vpi(i,1) = DAE.n + 2;
    DAE.n = DAE.n + 2;
  else
    a.vcs(i,1) = DAE.n + 1;
    DAE.n = DAE.n + 1;
  end
end
a.V0 = zeros(a.n,1);