fxcall.m

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

function Fxcall(a)

global DAE

if ~a.n, return, end

iTf = a.u./a.con(:,5);
Kv = a.u.*a.con(:,12);

DAE.Fx = DAE.Fx + sparse(a.x,a.x,-iTf,DAE.n,DAE.n);
DAE.Fy = DAE.Fy + sparse(a.x,a.vbus,-iTf.*iTf,DAE.n,DAE.m);
DAE.Gx = DAE.Gx + sparse(a.vbus,a.x,Kv,DAE.m,DAE.n);