fxcall.m

这是一个很适合研究和学习用的电力系统仿真软件

function Fxcall(a)

global DAE Bus Settings

if ~a.n, return, end

V = DAE.V(a.bus);
P0 = a.con(:,2);
Q0 = 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,-1./a.con(:,8),DAE.n,DAE.n);
DAE.Fy = DAE.Fy + sparse(a.x,a.bus,-1./a.con(:,8)./ ...
                         a.con(:,8)/Settings.rad,DAE.n,2*Bus.n);
DAE.Gx = DAE.Gx + sparse(a.bus,a.x,P0.*V.^ap.*bp.* ...
                         (1+a.Dw).^(bp-1),2*Bus.n,DAE.n);
DAE.Gx = DAE.Gx + sparse(a.bus+Bus.n,a.x,Q0.*V.^aq.*bq.* ...
                         (1+a.Dw).^(bq-1),2*Bus.n,DAE.n);