fxcall.m

来自「这是一个很适合研究和学习用的电力系统仿真软件」· M 代码 · 共 13 行

function Fxcall(a)

global DAE Bus

if ~a.n, return, end

iTf = 1./a.con(:,5);
Kv = 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.bus+Bus.n,-iTf.*iTf,DAE.n,2*Bus.n);
DAE.Gx = DAE.Gx + sparse(a.bus+Bus.n,a.x,Kv,2*Bus.n,DAE.n);