flows.m

来自「基于PSAT 软件的多目标最优潮流计算用于中小型电力系统的分析和管理」· M 代码 · 共 20 行

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;