fcall.m

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

function fcall(a)

global DAE

if ~a.n, return, end

V = DAE.y(a.vbus);
P0 = a.dat(:,1);
Q0 = a.dat(:,2);
V0 = a.dat(:,3);

DAE.f(a.xp) = -a.u.*(DAE.x(a.xp)./a.con(:,5) + ...
                     P0.*(V./V0).^a.con(:,7) - ...
                     P0.*(V./V0).^a.con(:,8));
DAE.f(a.xq) = -a.u.*(DAE.x(a.xq)./a.con(:,6) + ...
                     Q0.*(V./V0).^a.con(:,9) - ...
                     Q0.*(V./V0).^a.con(:,10));