fcall.m

来自「用于电力系统的一个很好的分析软件」· M 代码 · 共 29 行

function fcall(a)

global DAE

if ~a.n, return, end

T = DAE.x(a.T);
x = DAE.x(a.x);
G = DAE.y(a.G);
V = DAE.y(a.vbus);
Ki = a.con(:,4);
Ti = a.con(:,5);
T1 = a.con(:,6);
Ta = a.con(:,7);
Tref = a.con(:,8);
K1 = a.con(:,10);
G_max = a.con(:,9);

DAE.f(a.T) = a.u.*(Ta - T + K1.*G.*V.^2)./T1;
DAE.f(a.x) = a.u.*Ki.*(Tref-T)./Ti;

% anti-windup limits
idx = find(x >= G_max & DAE.f(a.x) > 0);
if idx, DAE.f(a.x(idx)) = 0; end
DAE.x(a.x) = min(DAE.x(a.x),G_max);
idx = find(x <= 0 & DAE.f(a.x) < 0);
if idx, DAE.f(a.x(idx)) = 0; end
DAE.x(a.x) = max(DAE.x(a.x),0);