dynidx.m

用于电力系统的一个很好的分析软件

function a = dynidx(a)

global DAE Syn

if ~a.n, return, end

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

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

for i = 1:a.n
  switch a.con(i,2)
   case 1
    a.tg1(i) = DAE.n + 1;
    a.tg2(i) = DAE.n + 2;
    a.tg3(i) = DAE.n + 3;
    DAE.n = DAE.n + 3;
   case 2
    a.tg(i) = DAE.n + 1;
    DAE.n = DAE.n + 1;
  end
end

a.wref = DAE.m + [1:a.n]';
DAE.m = DAE.m + a.n;

a.pm = Syn.pm(a.syn);