dynidx.m

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

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);