dynidx.m

这是一个很适合研究和学习用的电力系统仿真软件

function a = dynidx(a)

global DAE

if ~a.n, return, end

a.ty1 = find(a.con(:,2) == 1);
a.ty2 = find(a.con(:,2) == 2);

n1 = length(a.ty1);
n2 = length(a.ty2);

a.tg1 = zeros(n1,1);
a.tg2 = zeros(n1,1);
a.tg3 = zeros(n1,1);
a.tg  = zeros(n2,1);

for i = 1:n1
  a.tg1(i) = DAE.n + 1;
  a.tg2(i) = DAE.n + 2;
  a.tg3(i) = DAE.n + 3;
  DAE.n = DAE.n + 3;
end
if n2
  a.tg = DAE.n + [1:n2]';
  DAE.n = DAE.n + n2;
end