fxcall.m

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

function Fxcall(a)

global DAE Bus Settings

if ~a.n, return, end

Kpf = a.con(:,5);
Kpv = a.con(:,6);
Tpv = a.con(:,8);
Kqf = a.con(:,9);
Kqv = a.con(:,10);
Tpv = a.con(:,8);
Tqv = a.con(:,12);
Tfv = a.con(:,13);
Tft = a.con(:,14);
k = 0.5/pi/Settings.freq;

DAE.Fx = DAE.Fx + sparse(a.x,a.x,-1./Tfv,DAE.n,DAE.n);
DAE.Fx = DAE.Fx + sparse(a.y,a.y,-1./Tft,DAE.n,DAE.n);
DAE.Fy = DAE.Fy + sparse(a.x,a.bus+Bus.n,-1./Tfv./Tfv,DAE.n,2*Bus.n);
DAE.Fy = DAE.Fy + sparse(a.y,a.bus,-k./Tft./Tft,DAE.n,2*Bus.n);
DAE.Gx = DAE.Gx + sparse(a.bus,a.x,Kpv.*Tpv,2*Bus.n,DAE.n);
DAE.Gx = DAE.Gx + sparse(a.bus,a.y,Kpf,2*Bus.n,DAE.n);
DAE.Gx = DAE.Gx + sparse(a.bus+Bus.n,a.x,Kqv.*Tqv,2*Bus.n,DAE.n);
DAE.Gx = DAE.Gx + sparse(a.bus+Bus.n,a.y,Kqf,2*Bus.n,DAE.n);