gycall.m

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

function Gycall(a)

global Bus DAE

if ~a.n, return, end

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

DAE.J12 = DAE.J12 + sparse(a.bus,a.bus, ...
                           2.*Plz.*V1./V1_0.^2+Pli./V1_0,Bus.n,Bus.n);
DAE.J22 = DAE.J22 + sparse(a.bus,a.bus, ...
                           2.*Qlz.*V1./V1_0.^2+Qli./V1_0+Kv./Tf,Bus.n,Bus.n);