flows.m

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

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

global DAE

if ~a.n, return, end

V1 = DAE.V(a.bus1);
V2 = DAE.V(a.bus2);
cc = cos(DAE.a(a.bus1)-DAE.a(a.bus2));
den = max(sqrt(V1.^2+V2.^2-2.*V1.*V2.*cc),1e-6*ones(a.n,1));
u = DAE.x(a.vcs)./den;
Ps(a.line) = (1+u).*Ps(a.line);
Pr(a.line) = (1+u).*Pr(a.line);
Qs(a.line) = (1+u).*Qs(a.line);
Qr(a.line) = (1+u).*Qr(a.line);