fgamma.m

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

function g = fgamma(a)

global DAE jay

V = DAE.V(a.bus1).*exp(jay.*DAE.a(a.bus1)) - ...
    DAE.V(a.bus2).*exp(jay.*DAE.a(a.bus2));

theta = angle(V)-pi/2;

g = atan2(DAE.x(a.vq),DAE.x(a.vp))+theta-DAE.a(a.bus1);