gycall.m

来自「这是一个很适合研究和学习用的电力系统仿真软件」· M 代码 · 共 25 行

function Gycall(a)

global Bus DAE Settings

if ~a.n, return, end

V1 = DAE.V(a.bus);
Kpf = a.con(:,5);
Kpv = a.con(:,6);
alpha = a.con(:,7);
Tpv = a.con(:,8);
Kqf = a.con(:,9);
Kqv = a.con(:,10);
beta = a.con(:,11);
Tqv = a.con(:,12);
Tfv = a.con(:,13);
Tft = a.con(:,14);
k = 0.5/pi/Settings.freq;
V0 = a.dat(:,1);

DAE.J12 = DAE.J12 + sparse(a.bus,a.bus,Kpv.*((V1./V0).^alpha.*alpha./V1+Tpv./Tfv),Bus.n,Bus.n);
DAE.J11 = DAE.J11 + sparse(a.bus,a.bus,Kpf.*k./Tft,Bus.n,Bus.n);
DAE.J22 = DAE.J22 + sparse(a.bus,a.bus,Kqv.*((V1./V0).^beta.*beta./V1+Tqv./Tfv),Bus.n,Bus.n);
DAE.J21 = DAE.J21 + sparse(a.bus,a.bus,Kqf.*k./Tft,Bus.n,Bus.n);