gycall.m

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

function Gycall(a)

global Bus DAE Settings

if ~a.n, return, end

V = DAE.V(a.bus);
P0 = a.con(:,2);
Q0 = a.con(:,5);
ap = a.con(:,3);
aq = a.con(:,6);
bp = a.con(:,4);
bq = a.con(:,7);
dw_da = 1./a.con(:,8)/Settings.rad;

DAE.J11 = DAE.J11 + sparse(a.bus,a.bus,P0.*V.^ap.*bp.*dw_da.* ...
                           (1+a.Dw).^(bp-1),Bus.n,Bus.n);
DAE.J12 = DAE.J12 + sparse(a.bus,a.bus,P0.*ap.*V.^(ap-1).* ...
                           (1+a.Dw).^bp,Bus.n,Bus.n);
DAE.J22 = DAE.J22 + sparse(a.bus,a.bus,Q0.*aq.*V.^(aq-1).* ...
                           (1+a.Dw).^bq,Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(a.bus,a.bus,Q0.*V.^aq.*bq.*dw_da.* ...
                           (1+a.Dw).^(bq-1),Bus.n,Bus.n);