www.gusucode.com > 基于lingo求所以解，对潮流计算求出所有解 > matpower4.1/extras/se/outputpfsoln.m
function outputpfsoln(baseMVA, bus, gen, branch, converged, et, type_solver, iterNum)
%OUTPUTPFSOLN Output power flow solution.
% MATPOWER
% $Id: outputpfsoln.m,v 1.4 2010/04/26 19:45:26 ray Exp $
% by Rui Bo
% and Ray Zimmerman, PSERC Cornell
% Copyright (c) 1996-2010 by Power System Engineering Research Center (PSERC)
% Copyright (c) 2009-2010 by Rui Bo
%
% This file is part of MATPOWER.
% See http://www.pserc.cornell.edu/matpower/ for more info.
%
% MATPOWER is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published
% by the Free Software Foundation, either version 3 of the License,
% or (at your option) any later version.
%
% MATPOWER is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with MATPOWER. If not, see <http://www.gnu.org/licenses/>.
%
% Additional permission under GNU GPL version 3 section 7
%
% If you modify MATPOWER, or any covered work, to interface with
% other modules (such as MATLAB code and MEX-files) available in a
% MATLAB(R) or comparable environment containing parts covered
% under other licensing terms, the licensors of MATPOWER grant
% you additional permission to convey the resulting work.
%% define named indices into bus, gen, branch matrices
[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, ...
GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;
[F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, ...
RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST] = idx_brch;
fd = 1; % output to screen
%% sizes of things
nb = size(bus, 1); %% number of buses
nl = size(branch, 1); %% number of branches
ng = size(gen, 1); %% number of generators
%% parameters
ong = find( gen(:, GEN_STATUS) > 0);
nzld = find(bus(:, PD) | bus(:, QD));
%% calculate losses
loss = branch(:, PF) + j*branch(:, QF) + branch(:, PT) + j*branch(:, QT);
%% ---output case and solver information
fprintf(fd, '\n\n');
if type_solver == 1 % newton's method
fprintf(fd, 'Newton''s method is chosen to solve Power Flow.\n');
elseif type_solver == 2 % decoupled method
fprintf(fd, 'Decoupled method is chosen to solve Power Flow.\n');
else
fprintf('Error: unknow ''type_solver''.\n');
pause
end
if converged
fprintf(fd, '\nConverged in %.2f seconds\n', et);
else
fprintf(fd, '\nDid not converge (%.2f seconds)\n', et);
end
fprintf(fd, '\n[iteration number]: %d\n', iterNum);
%% ---output generation information
fprintf(fd, '\n================================================================================');
fprintf(fd, '\n| Generator Data |');
fprintf(fd, '\n================================================================================');
fprintf(fd, '\n Gen Bus Status Pg Qg ');
fprintf(fd, '\n # # (MW) (MVAr) ');
fprintf(fd, '\n---- ----- ------ -------- --------');
for k = 1:length(ong)
i = ong(k);
fprintf(fd, '\n%3d %6d %2d ', i, gen(i, GEN_BUS), gen(i, GEN_STATUS));
if gen(i, GEN_STATUS) > 0 & (gen(i, PG) | gen(i, QG))
fprintf(fd, '%10.2f%10.2f', gen(i, PG), gen(i, QG));
else
fprintf(fd, ' - - ');
end
end
fprintf(fd, '\n -------- --------');
fprintf(fd, '\n Total: %9.2f%10.2f', sum(gen(ong, PG)), sum(gen(ong, QG)));
fprintf(fd, '\n');
%% ---output bus information
fprintf(fd, '\n================================================================================');
fprintf(fd, '\n| Bus Data |');
fprintf(fd, '\n================================================================================');
fprintf(fd, '\n Bus Voltage Generation Load ');
fprintf(fd, '\n # Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr)');
fprintf(fd, '\n----- ------- -------- -------- -------- -------- --------');
for i = 1:nb
fprintf(fd, '\n%5d%7.3f%9.3f', bus(i, [BUS_I, VM, VA]));
g = find(gen(:, GEN_STATUS) > 0 & gen(:, GEN_BUS) == bus(i, BUS_I));
if ~isempty(g)
fprintf(fd, '%10.2f%10.2f', sum(gen(g, PG)), sum(gen(g, QG)));
else
fprintf(fd, ' - - ');
end
if bus(i, PD) | bus(i, QD)
fprintf(fd, '%10.2f%10.2f ', bus(i, [PD, QD]));
else
fprintf(fd, ' - - ');
end
end
fprintf(fd, '\n -------- -------- -------- --------');
fprintf(fd, '\n Total: %9.2f %9.2f %9.2f %9.2f', ...
sum(gen(ong, PG)), sum(gen(ong, QG)), ...
sum(bus(nzld, PD)), ...
sum(bus(nzld, QD)));
fprintf(fd, '\n');
%% ---output bus information
fprintf(fd, '\n================================================================================');
fprintf(fd, '\n| Branch Data |');
fprintf(fd, '\n================================================================================');
fprintf(fd, '\nBrnch From To From Bus Injection To Bus Injection Loss (I^2 * Z) ');
fprintf(fd, '\n # Bus Bus P (MW) Q (MVAr) P (MW) Q (MVAr) P (MW) Q (MVAr)');
fprintf(fd, '\n----- ----- ----- -------- -------- -------- -------- -------- --------');
fprintf(fd, '\n%4d%7d%7d%10.2f%10.2f%10.2f%10.2f%10.3f%10.2f', ...
[ [1:nl]', branch(:, [F_BUS, T_BUS]), ...
branch(:, [PF, QF]), branch(:, [PT, QT]), ...
real(loss), imag(loss) ...
]');
fprintf(fd, '\n -------- --------');
fprintf(fd, '\n Total:%10.3f%10.2f', ...
sum(real(loss)), sum(imag(loss)));
fprintf(fd, '\n');