dynidx.m

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

function a = dynidx(a)

global DAE

if ~a.n, return, end

a.bcv = zeros(a.n,1);
a.alpha = zeros(a.n,1);
a.vm = zeros(a.n,1);
type = a.con(:,5);
for i = 1:a.n
  if type(i) == 1
    a.bcv(i) = DAE.n + 1;
    DAE.n = DAE.n + 1;
  elseif type(i) == 2
    a.alpha(i) = DAE.n + 1;
    a.vm(i) = DAE.n + 2;
    DAE.n = DAE.n + 2;
  end
end
a.bcv = a.bcv(find(a.bcv));
a.alpha = a.alpha(find(a.alpha));
a.vm = a.vm(find(a.vm));
a.Be = zeros(a.n,1);
a.Vref = zeros(a.n,1);