xfirst.m

来自「电力系统分析计算程序」· M 代码 · 共 15 行

function xfirst(a)

global DAE

if ~a.n, return, end

DAE.x(a.alpha) = 0;
DAE.x(a.Pm) = 0;
idx = find(~a.con(:,7));
if ~isempty(idx)
  warn(a,idx,'Measurement time constant Tm cannot be 0. Tm = 1e-3 will be used.')
  a.con(idx,7) = 1e-3;
end