pfsoln.m

粒子群优化算法(PSO)是一种进化计算技术(evolutionary computation)

function [bus, gen, branch] = pfsoln(baseMVA, bus0, gen0, branch0, Ybus, Yf, Yt, V, ref, pv, pq);
%PFSOLN  Updates bus, gen, branch data structures to match power flow soln.
%   [bus, gen, branch] = pfsoln(baseMVA, bus0, gen0, branch0, ...
%                                   Ybus, Yf, Yt, V, ref, pv, pq)

%   MATPOWER Version 2.0
%   by Ray Zimmerman, PSERC Cornell    12/19/97
%   Copyright (c) 1996, 1997 by Power System Engineering Research Center (PSERC)
%   See http://www.pserc.cornell.edu/ for more info.

%% constants
j = sqrt(-1);
nl = size(branch0, 1);		%% number of lines

%% ----- define named indices into bus, gen, branch matrices
PQ=1; PV=2; REF=3; NONE=4; BUS_I=1; BUS_TYPE=2; PD=3; QD=4; GS=5; BS=6; BUS_AREA=7;
VM=8; VA=9; BASE_KV=10;	ZONE=11; VMAX=12; VMIN=13; LAM_P=14; LAM_Q=15; MU_VMAX=16; MU_VMIN=17; 

GEN_BUS=1; PG=2; QG=3; QMAX=4; QMIN=5; VG=6; MBASE=7; GEN_STATUS=8; PMAX=9; PMIN=10;
MU_PMAX=11; MU_PMIN=12; MU_QMAX=13; MU_QMIN=14; % MU: KT multiplier

F_BUS=1;  T_BUS=2;  BR_R=3;   BR_X=4; BR_B=5; % standard notation (in input)
RATE_A=6; RATE_B=7; RATE_C=8; TAP=9;  SHIFT=10; BR_STATUS=11; % standard notation (in input)
PF=12; QF=13; PT=14; QT=15; MU_SF=16; MU_ST=17; % MU_SF: idx of MU on MVA limit at "f" bus (in opf)

%% initialize return values
bus		= bus0;
gen		= gen0;
branch	= branch0;
 
%%----- update bus voltages -----
bus(:, VM) = abs(V);
bus(:, VA) = angle(V) * 180 / pi;

%%----- update Qg for all gens and Pg for swing bus -----
%% generator info
on = find(gen(:, GEN_STATUS));				%% which generators are on?
gbus = gen(on, GEN_BUS);					%% what buses are they at?
refgen = find(gen(on, GEN_BUS) == ref);		%% which is the reference gen?

%% compute total injected bus powers
%% This is slow in Matlab 5 ...
% Sg = V(gbus) .* conj(Ybus(gbus, :) * V);
%% ... so we do this instead ...
temp = Ybus.';
Sg = V(gbus) .* conj(temp(:, gbus).' * V);

%% update Qg for all generators
gen(:, QG) = zeros(size(gen, 1), 1);				%% zero out all Qg
gen(on, QG) = imag(Sg) * baseMVA + bus(gbus, QD);	%% inj Q + local Qd

%% update Pg for swing bus
gen(on(refgen), PG) = real(Sg(refgen)) * baseMVA + bus(ref, PD);	%% inj P + local Pd

%%----- update/compute branch power flows -----
Sf = V(branch(:, F_BUS)) .* conj(Yf * V) * baseMVA;	%% complex power at "from" bus
St = V(branch(:, T_BUS)) .* conj(Yt * V) * baseMVA;	%% complex power injected at "to" bus
branch(:, [PF, QF, PT, QT]) = [real(Sf) imag(Sf) real(St) imag(St)];

return;