fxcall.m

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

function Fxcall(a)

global DAE Settings

if ~a.n, return, end

V = DAE.y(a.vbus);
P0 = a.u.*a.con(:,2);
Q0 = a.u.*a.con(:,5);
ap = a.con(:,3);
aq = a.con(:,6);
bp = a.con(:,4);
bq = a.con(:,7);

DAE.Fx = DAE.Fx + sparse(a.x,a.x,-a.u./a.con(:,8),DAE.n,DAE.n);
DAE.Fy = DAE.Fy + sparse(a.x,a.bus,-a.u./a.con(:,8)./ ...
                         a.con(:,8)/(2*pi*Settings.freq),DAE.n,DAE.m);
DAE.Gx = DAE.Gx + sparse(a.bus,a.x,P0.*V.^ap.*bp.* ...
                         (1+a.Dw).^(bp-1),DAE.m,DAE.n);
DAE.Gx = DAE.Gx + sparse(a.vbus,a.x,Q0.*V.^aq.*bq.* ...
                         (1+a.Dw).^(bq-1),DAE.m,DAE.n);