function printpf(baseMVA, bus, gen, branch, f, success, et, fd, mpopt)%PRINTPF   Prints power flow results.%   printpf(baseMVA, bus, gen, branch, f, success, et, fd, mpopt) prints%   powerflow results to fd (a file descriptor which defaults to STDOUT).%   mpopt is a MATPOWER options vector (see 'help mpoption' for details).%   Uses default options if this parameter is not given. The objective%   function value is given in f and the elapsed time (seconds to compute%   opf) in et.%   MATPOWER%   $Id: printpf.m,v 1.27 2007/06/26 16:16:17 ray Exp $%   by Ray Zimmerman, PSERC Cornell%   Copyright (c) 1996-2004 by Power System Engineering Research Center (PSERC)%   See http://www.pserc.cornell.edu/matpower/ for more info.%%----- initialization -----%% default argumentsif nargin < 9    mpopt = mpoption;   %% use default options    if nargin < 8        fd = 1;         %% print to stdio by default    endendif isempty(f)    isOPF = 0;      %% have only simple PF dataelse    isOPF = 1;      %% have OPF dataend%% optionsdc              = mpopt(10);        %% use DC formulation?OUT_ALL         = mpopt(32);OUT_ANY         = OUT_ALL == 1;     %% set to true if any pretty output is to be generatedOUT_SYS_SUM     = OUT_ALL == 1 | (OUT_ALL == -1 & mpopt(33));OUT_AREA_SUM    = OUT_ALL == 1 | (OUT_ALL == -1 & mpopt(34));OUT_BUS         = OUT_ALL == 1 | (OUT_ALL == -1 & mpopt(35));OUT_BRANCH      = OUT_ALL == 1 | (OUT_ALL == -1 & mpopt(36));OUT_GEN         = OUT_ALL == 1 | (OUT_ALL == -1 & mpopt(37));OUT_ANY         = OUT_ANY | (OUT_ALL == -1 & (OUT_SYS_SUM | OUT_AREA_SUM | OUT_BUS | OUT_BRANCH | OUT_GEN));if OUT_ALL == -1    OUT_ALL_LIM = mpopt(38);elseif OUT_ALL == 1    OUT_ALL_LIM = 2;else    OUT_ALL_LIM = 0;endOUT_ANY         = OUT_ANY | OUT_ALL_LIM >= 1;if OUT_ALL_LIM == -1    OUT_V_LIM       = mpopt(39);    OUT_LINE_LIM    = mpopt(40);    OUT_PG_LIM      = mpopt(41);    OUT_QG_LIM      = mpopt(42);else    OUT_V_LIM       = OUT_ALL_LIM;    OUT_LINE_LIM    = OUT_ALL_LIM;    OUT_PG_LIM      = OUT_ALL_LIM;    OUT_QG_LIM      = OUT_ALL_LIM;endOUT_ANY         = OUT_ANY | (OUT_ALL_LIM == -1 & (OUT_V_LIM | OUT_LINE_LIM | OUT_PG_LIM | OUT_QG_LIM));OUT_RAW         = mpopt(43);%% 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, PC1, PC2, QC1MIN, QC1MAX, ...    QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = 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, ...    ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;%% constantj = sqrt(-1);%% internal bus numbere2i = zeros(max(bus(:, BUS_I)), 1);     %% need internal bus numbering for a seconde2i(bus(:, BUS_I)) = [1:size(bus, 1)]';%% sizes of thingsnb = size(bus, 1);      %% number of busesnl = size(branch, 1);   %% number of branchesng = size(gen, 1);      %% number of generators%% zero out some data to make printout consistent for DC caseif dc    bus(:, [QD, BS])            = zeros(nb, 2);    gen(:, [QG, QMAX, QMIN])    = zeros(ng, 3);    branch(:, [BR_R, BR_B])     = zeros(nl, 2);end%% parametersties = find(bus(e2i(branch(:, F_BUS)), BUS_AREA) ~= bus(e2i(branch(:, T_BUS)), BUS_AREA));                        %% area inter-tiestap = ones(nl, 1);                              %% default tap ratio = 1 for linesxfmr = find(branch(:, TAP));                    %% indices of transformerstap(xfmr) = branch(xfmr, TAP);                  %% include transformer tap ratiostap = tap .* exp(j*pi/180 * branch(:, SHIFT)); %% add phase shiftersnzld = find(bus(:, PD) | bus(:, QD));sorted_areas = sort(bus(:, BUS_AREA));s_areas = sorted_areas([1; find(diff(sorted_areas))+1]);    %% area numbersna = length(s_areas);                           %% number of areasnzsh = find(bus(:, GS) | bus(:, BS));allg = find( ~isload(gen) );ong  = find( gen(:, GEN_STATUS) > 0 & ~isload(gen) );onld = find( gen(:, GEN_STATUS) > 0 &  isload(gen) );V = bus(:, VM) .* exp(sqrt(-1) * pi/180 * bus(:, VA));out = find(branch(:, BR_STATUS) == 0);          %% out-of-service branchesnout = length(out);if dc    loss = zeros(nl,1);else    loss = baseMVA * abs(V(e2i(branch(:, F_BUS))) ./ tap - V(e2i(branch(:, T_BUS)))) .^ 2 ./ ...                (branch(:, BR_R) - j * branch(:, BR_X));endfchg = abs(V(e2i(branch(:, F_BUS))) ./ tap) .^ 2 .* branch(:, BR_B) * baseMVA / 2;tchg = abs(V(e2i(branch(:, T_BUS)))       ) .^ 2 .* branch(:, BR_B) * baseMVA / 2;loss(out) = zeros(nout, 1);fchg(out) = zeros(nout, 1);tchg(out) = zeros(nout, 1);%%----- print the stuff -----if OUT_ANY    %% convergence & elapsed time    if success                     fprintf(fd, '\n****潮流計算時間 %.2f 秒***', et);    else        fprintf(fd, '\nDid not converge (%.2f seconds)\n', et);    end        %% objective function value    if isOPF        fprintf(fd, '\nObjective Function Value = %.2f $/hr', f);    endendif OUT_SYS_SUM   %% fprintf(fd, '\n================================================================================');   %% fprintf(fd, '\n|     系統概況                                                           |');   %% fprintf(fd, '\n================================================================================');   %% fprintf(fd, '\n\n基本情況                How much?              P (MW)            Q (MVAr)');   %% fprintf(fd, '\n---------------------    -------------------  -------------  -----------------');   %% fprintf(fd, '\n母線           %6d     總容量   %7.1f       %7.1f to %.1f', nb, sum(gen(allg, PMAX)), sum(gen(allg, QMIN)), sum(gen(allg, QMAX)));   %% fprintf(fd, '\n發電機         %5d     On-line Capacity     %7.1f       %7.1f to %.1f', length(allg), sum(gen(ong, PMAX)), sum(gen(ong, QMIN)), sum(gen(ong, QMAX)));    %%fprintf(fd, '\nCommitted Gens %5d     Generation (actual)  %7.1f           %7.1f', length(ong), sum(gen(ong, PG)), sum(gen(ong, QG)));   %% fprintf(fd, '\n負荷           %5d     負荷                 %7.1f           %7.1f', length(nzld)+length(onld), sum(bus(nzld, PD))-sum(gen(onld, PG)), sum(bus(nzld, QD))-sum(gen(onld, QG)));    %%fprintf(fd, '\n  Fixed        %5d       Fixed              %7.1f           %7.1f', length(nzld), sum(bus(nzld, PD)), sum(bus(nzld, QD)));    %%fprintf(fd, '\n  Dispatchable %5d       Dispatchable       %7.1f of %-7.1f%7.1f', length(onld), -sum(gen(onld, PG)), -sum(gen(onld, PMIN)), -sum(gen(onld, QG)));    %%fprintf(fd, '\nShunts         %5d     Shunt (inj)          %7.1f           %7.1f', length(nzsh), ...      %%  -sum(bus(nzsh, VM) .^ 2 .* bus(nzsh, GS)), sum(bus(nzsh, VM) .^ 2 .* bus(nzsh, BS)) );    %%fprintf(fd, '\n支路           %5d     線路損耗 (I^2 * Z)     %8.2f          %8.2f', nl, sum(real(loss)), sum(imag(loss)) );    %%fprintf(fd, '\n變壓器         %5d     支路充電功率 (inj)       -            %7.1f', length(xfmr), sum(fchg) + sum(tchg) );    %%fprintf(fd, '\nInter-ties     %5d     Total Inter-tie Flow %7.1f           %7.1f', length(ties), sum(abs(branch(ties, PF)-branch(ties, PT))) / 2, sum(abs(branch(ties, QF)-branch(ties, QT))) / 2);    %%fprintf(fd, '\nAreas          %5d', length(s_areas));