calc_v_i_sum

calc_v_i_sum(Vslack, nb, nl, f, Zb, Ybf, Ybt, Yd, Sd, pv, Pg, Vg, mpopt)

calc_v_i_sum() - Solves the power flow using the current summation method.

[V, Qpv, Sf, St, Sslack, iter, success] = calc_v_i_sum(Vslack,nb,nl,f,Zb,Ybf,Ybt,Yd,Sd,pv,Pg,Vg,tol,iter_max)

Solves for bus voltages, generator reactive power, branch active and
reactive power flows and slack bus active and reactive power. The input
data consist of slack bus voltage, vector "from bus" indices, branch
impedance and shunt admittance, vector of bus shunt admittances and
load demand, as well as vectors with indicies of PV buses with their
specified voltages and active powers. It is assumed that the branches
are ordered using the principle of oriented ordering: indicies of
sending nodes are smaller then the indicies of the receiving nodes. The
branch index is equal to the index of their receiving node. Branch
addmittances are added in Yd and treated as constant admittance bus
loads. The applied method is current summation taken from:
D. Shirmohammadi, H. W. Hong, A. Semlyen and G. X. Luo,
"A compensation-based power flow method for weakly meshed distribution
and transmission networks," IEEE Transactions on Power Systems,
vol. 3, no. 2, pp. 753-762, May 1988. https://doi.org/10.1109/59.192932
and
G. X. Luo and A. Semlyen, "Efficient load flow for large weakly meshed
networks," IEEE Transactions on Power Systems, vol. 5, no. 4,
pp. 1309-1316, Nov 1990. https://doi.org/10.1109/59.99382

See also radial_pf().