fcall.m

来自「这是一个很适合研究和学习用的电力系统仿真软件」· M 代码 · 共 12 行

function fcall(a)

global DAE

if ~a.n, return, end

x = DAE.x(a.x);
V1 = DAE.V(a.bus);
iTf = 1./a.con(:,5);

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