flows.m

用于电力系统的一个很好的分析软件

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

global DAE

if nargin == 7
  varargout{1} = [varargin{1}; a.bus1];
  varargout{2} = [varargin{2}; a.bus2];
end

if ~a.n, return, end

Idc = a.u.*DAE.x(a.Idc);
Sr = a.u.*DAE.y(a.Sr);
Si = a.u.*DAE.y(a.Si);
Vrdc = DAE.y(a.Vrdc);
Vidc = DAE.y(a.Vidc);

VISn = a.dat(:,8);

P1 = VISn.*Vrdc.*Idc;
P2 = VISn.*Vidc.*Idc;
Q1 = sqrt(Sr.*Sr-P1.*P1);
Q2 = sqrt(Si.*Si-P2.*P2);

Ps = [Ps; P1];
Pr = [Pr;-P2];
Qs = [Qs; Q1];
Qr = [Qr; Q2];