dynidx.m

来自「基于PSAT 软件的多目标最优潮流计算用于中小型电力系统的分析和管理」· M 代码 · 共 28 行

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