gcall.m

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

function gcall(a)

global Bus DAE Settings

if ~a.n, return, end

x = DAE.x(a.x);
y = DAE.x(a.y);
V1 = DAE.V(a.bus);
t1 = DAE.a(a.bus);
Kpf = a.con(:,5);
Kpv = a.con(:,6);
alpha = a.con(:,7);
Tpv = a.con(:,8);
Kqf = a.con(:,9);
Kqv = a.con(:,10);
beta = a.con(:,11);
Tqv = a.con(:,12);
Tfv = a.con(:,13);
Tft = a.con(:,14);
k = 0.5/pi/Settings.freq;
V0 = a.dat(:,1);
t0 = a.dat(:,2);


DAE.gp = DAE.gp + sparse(a.bus,1,Kpf.*(y+k.*(t1-t0)./Tft)+Kpv.*((V1./V0).^alpha+Tpv.*(x+V1./Tfv)),Bus.n,1);
DAE.gq = DAE.gq + sparse(a.bus,1,Kqf.*(y+k.*(t1-t0)./Tft)+Kqv.*((V1./V0).^beta+Tqv.*(x+V1./Tfv)),Bus.n,1);