gycall.m

电力系统分析计算程序

function Gycall(a)

global DAE Settings

if ~a.n, return, end

V = DAE.y(a.vbus);
P0 = a.u.*a.con(:,2);
Q0 = a.u.*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)/(2*pi*Settings.freq);

DAE.Gy = DAE.Gy + sparse(a.bus,a.bus,P0.*V.^ap.*bp.*dw_da.* ...
                           (1+a.Dw).^(bp-1),DAE.m,DAE.m);
DAE.Gy = DAE.Gy + sparse(a.bus,a.vbus,P0.*ap.*V.^(ap-1).* ...
                           (1+a.Dw).^bp,DAE.m,DAE.m);
DAE.Gy = DAE.Gy + sparse(a.vbus,a.vbus,Q0.*aq.*V.^(aq-1).* ...
                           (1+a.Dw).^bq,DAE.m,DAE.m);
DAE.Gy = DAE.Gy + sparse(a.vbus,a.bus,Q0.*V.^aq.*bq.*dw_da.* ...
                           (1+a.Dw).^(bq-1),DAE.m,DAE.m);