dynidx.m

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

function a = dynidx(a)

global DAE

if ~a.n, return, end

a.vp = a.vp(find(a.vp));
a.vq = a.vq(find(a.vq));
a.iq = a.iq(find(a.iq));
for i = 1:a.n
  a.vp(i) = DAE.n + 1;
  a.vq(i) = DAE.n + 2;     
  a.iq(i) = DAE.n + 3;
  DAE.n = DAE.n + 3;
end
a.gamma = zeros(a.n,1);
a.vp0 = zeros(a.n,1);
a.vq0 = zeros(a.n,1);
a.Vref = zeros(a.n,1);