gcall.m

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

function gcall(a)

global Bus DAE

if ~a.n, return, end

V = DAE.V(a.bus);
Dw = domega(a);

DAE.gp = DAE.gp + sparse(a.bus,1,a.con(:,2).*(V.^a.con(:,3)).* ...
                         (1+Dw).^a.con(:,4),Bus.n,1);
DAE.gq = DAE.gq + sparse(a.bus,1,a.con(:,5).*(V.^a.con(:,6)).* ...
                         (1+Dw).^a.con(:,7),Bus.n,1);