【Matlab源码】分布式电源接入对配电网的影响.rar

function printpf(baseMVA, bus, gen, branch, f, success, et, fd, mpopt) %PRINTPF 打印潮流计算的信息 % % printpf(results) % printpf(results, fd) % printpf(results, fd, mpopt) % printpf(baseMVA, bus, gen, branch, f, success, et) % printpf(baseMVA, bus, gen, branch, f, success, et, fd) % printpf(baseMVA, bus, gen, branch, f, success, et, fd, mpopt) %%----- initialization ----- %% default arguments if isstruct(baseMVA) have_results_struct = 1; results = baseMVA; if nargin < 3 || isempty(gen) mpopt = mpoption; %% use default options else mpopt = gen; end if mpopt(32) == 0 && mpopt(43) == 0 %% OUT_ALL or OUT_RAW return; %% nothin' to see here, bail out now end if nargin < 2 || isempty(bus) fd = 1; %% print to stdio by default else fd = bus; end [baseMVA, bus, gen, branch, success, et] = ... deal(results.baseMVA, results.bus, results.gen, results.branch, ... results.success, results.et); if isfield(results, 'f') && ~isempty(results.f) f = results.f; else f = []; end else have_results_struct = 0; if nargin < 9 mpopt = mpoption; %% use default options if nargin < 8 fd = 1; %% print to stdio by default end end if mpopt(32) == 0 && mpopt(43) == 0 %% OUT_ALL or OUT_RAW return; %% nothin' to see here, bail out now end end isOPF = ~isempty(f); %% FALSE -> only simple PF data, TRUE -> OPF data %% options isDC = mpopt(10); %% use DC formulation? OUT_ALL = mpopt(32); OUT_ANY = OUT_ALL == 1; %% set to true if any pretty output is to be generated OUT_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; end OUT_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; end OUT_ANY = OUT_ANY || (OUT_ALL_LIM == -1 && (OUT_V_LIM || OUT_LINE_LIM || OUT_PG_LIM || OUT_QG_LIM)); OUT_RAW = mpopt(43); ptol = 1e-6; %% tolerance for displaying shadow prices %% 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; %% internal bus number i2e = bus(:, BUS_I); e2i = sparse(max(i2e), 1); e2i(i2e) = (1:size(bus, 1))'; %% sizes of things nb = size(bus, 1); %% number of buses nl = size(branch, 1); %% number of branches ng = size(gen, 1); %% number of generators %% zero out some data to make printout consistent for DC case if isDC bus(:, [QD, BS]) = zeros(nb, 2); gen(:, [QG, QMAX, QMIN]) = zeros(ng, 3); branch(:, [BR_R, BR_B]) = zeros(nl, 2); end %% parameters ties = find(bus(e2i(branch(:, F_BUS)), BUS_AREA) ~= bus(e2i(branch(:, T_BUS)), BUS_AREA)); %% area inter-ties tap = ones(nl, 1); %% default tap ratio = 1 for lines xfmr = find(branch(:, TAP)); %% indices of transformers tap(xfmr) = branch(xfmr, TAP); %% include transformer tap ratios tap = tap .* exp(1j*pi/180 * branch(:, SHIFT)); %% add phase shifters nzld = find(bus(:, PD) | bus(:, QD)); sorted_areas = sort(bus(:, BUS_AREA)); s_areas = sorted_areas([1; find(diff(sorted_areas))+1]); %% area numbers nzsh = 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 branches nout = length(out); if isDC loss = zeros(nl,1); else loss = baseMVA * abs(V(e2i(branch(:, F_BUS))) ./ tap - V(e2i(branch(:, T_BUS)))) .^ 2 ./ ... (branch(:, BR_R) - 1j * branch(:, BR_X)); end fchg = 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, '\nConverged in %.2f seconds', 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); end end if OUT_SYS_SUM fprintf(fd, '\n================================================================================'); fprintf(fd, '\n| System Summary |'); fprintf(fd, '\n================================================================================'); fprintf(fd, '\n\nHow many? How much? P (MW) Q (MVAr)'); fprintf(fd, '\n--------------------- ------------------- ------------- -----------------'); fprintf(fd, '\nBuses %6d Total Gen Capacity %7.1f %7.1f to %.1f', nb, sum(gen(allg, PMAX)), sum(gen(allg, QMIN)), sum(gen(allg, QMAX))); fprintf(fd, '\nGenerators %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, '\nLoads %5d Load %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, '\nBranches %5d Losses (I^2 * Z) %8.2f %8.2f', nl, sum(real(loss)), sum(imag(loss)) ); fprintf(fd, '\nTransformers %5d Branch Charging (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)); fprintf(fd, '\n'); fprintf(fd, '\n Minimum