windup.m

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

function windup(a)

if ~a.n, return, end

global DAE

tyb = find(a.con(:,2) > 3);
if isempty(tyb), return, end
vamax = a.con(:,16);
fm_windup(a.va(tyb),vamax(tyb),0)

S2 = a.con(:,22);
S2 = (((DAE.x(a.omega)-1) < 0) | S2) & S2 >= 0;
idx = find(S2(tyb));
if isempty(idx), return, end
vathr = a.con(:,17);
fm_windup(a.va(tyb(idx)),vathr(tyb(idx)),0)