gcall.m

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

function gcall(a)

global Bus DAE

if ~a.n, return, end

x = DAE.x(a.x);
V1 = DAE.V(a.bus);
Tf = a.con(:,5);
Plz = a.con(:,6);
Pli = a.con(:,7);
Plp = a.con(:,8);
Qlz = a.con(:,9);
Qli = a.con(:,10);
Qlp = a.con(:,11);
Kv = a.con(:,12);
V1_0 = a.dat(:,1);

DAE.gp = DAE.gp + sparse(a.bus,1,Plz.*(V1./V1_0).^2 + ...
                         Pli.*(V1./V1_0)+Plp,Bus.n,1);
DAE.gq = DAE.gq + sparse(a.bus,1,Qlz.*(V1./V1_0).^2 + ...
                         Qli.*(V1./V1_0)+Qlp+Kv.*(x+V1./Tf),Bus.n,1);