fcall.m

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

function fcall(a)

global DAE

if ~a.n, return, end

x = DAE.x(a.x);
V1 = DAE.y(a.vbus);
iTf = a.u./a.con(:,5);

DAE.f(a.x) = -(V1.*iTf+x).*iTf;