dynidx.m

电力系统分析计算程序

function a = dynidx(a)

global DAE

if ~a.n, return, end

a.theta_p = zeros(a.n,1);
a.omega_m = zeros(a.n,1);
a.iqs = zeros(a.n,1);
a.idc = zeros(a.n,1);
for i = 1:a.n
  a.omega_m(i) = DAE.n + 1;
  a.theta_p(i) = DAE.n + 2;
  a.iqs(i) = DAE.n + 3;
  a.idc(i) = DAE.n + 4;
  DAE.n = DAE.n + 4;
end

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