flows.m

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

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

global DAE

if ~a.n, return, end

% update B
a.B = btcsc(a);

V1 = DAE.V(a.bus1);
V2 = DAE.V(a.bus2);
ss = sin(DAE.a(a.bus1)-DAE.a(a.bus2));
cc = cos(DAE.a(a.bus1)-DAE.a(a.bus2));

P1 = V1.*V2.*ss.*a.B;
Ps(a.line) = Ps(a.line) + P1;
Pr(a.line) = Pr(a.line) - P1;
Qs(a.line) = Qs(a.line) + V1.*(V1-V2.*cc).*a.B;
Qr(a.line) = Qs(a.line) + V2.*(V2-V1.*cc).*a.B;