fcall.m

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

function fcall(a)

global DAE Settings

if ~a.n, return, end

x = DAE.x(a.x);
y = DAE.x(a.y);
V1 = DAE.V(a.bus);
t1 = DAE.a(a.bus);
Tfv = a.con(:,13);
Tft = a.con(:,14);
k = 0.5/pi/Settings.freq;
V0 = a.dat(:,1);
t0 = a.dat(:,2);

DAE.f(a.x) = (-V1./Tfv-x)./Tfv;
DAE.f(a.y) = (-k.*(t1-t0)./Tft-y)./Tft;