fxcall.m

来自「电力系统分析计算程序」· M 代码 · 共 16 行

function Fxcall(a)

global DAE

if ~a.n, return, end

iTf = 1./a.con(:,2);
iTw = 1./a.con(:,3);
k = a.u.*a.dat(:,2);

DAE.Fx = DAE.Fx + sparse(a.x,a.x,-iTf,DAE.n,DAE.n);
DAE.Fx = DAE.Fx + sparse(a.w,a.x,-a.u.*iTw,DAE.n,DAE.n);
DAE.Fx = DAE.Fx + sparse(a.w,a.w,-iTw,DAE.n,DAE.n);
DAE.Fy = DAE.Fy + sparse(a.x,a.bus,k.*iTf,DAE.n,DAE.m);
DAE.Fy = DAE.Fy + sparse(a.w,a.bus,k.*iTw,DAE.n,DAE.m);