flows.m

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

function [Pr,Ps,Qr,Qs] = flows(a,Pr,Ps,Qr,Qs)

global DAE

if ~a.n, return, end

g = fgamma(a);
V1 = DAE.V(a.bus1);
V2 = DAE.V(a.bus2);
theta = DAE.a(a.bus1)-DAE.a(a.bus2)+g;
ss = sin(theta);
cc = cos(theta);
c1 = sqrt(DAE.x(a.vp).^2+DAE.x(a.vq).^2).*a.y;
Ps(a.line) = Ps(a.line) + c1.*V2.*ss;
Pr(a.line) = Pr(a.line) - c1.*V2.*ss;
Qs(a.line) = Qs(a.line) + c1.*V1.*cos(g)-DAE.x(a.iq).*V1;
Qr(a.line) = Qr(a.line) - c1.*V2.*cc;