fxcall.m

电力系统分析计算程序

function Fxcall(p)

global DAE Exc

if ~p.n, return, end

z = DAE.x(p.v) > 0 & DAE.x(p.v) < p.con(:,7);

DAE.Fx = DAE.Fx + sparse(p.v,p.v,-1e-6,DAE.n,DAE.n);
DAE.Fy = DAE.Fy + sparse(p.v,p.If,z./p.con(:,2),DAE.n,DAE.m);
DAE.Gx = DAE.Gx - sparse(p.vref,p.v,z,DAE.m,DAE.n);