calc_v_y_sum

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

calc_v_y_sum() - Solves the power flow using the admittance summation method.

[V, Qpv, Sf, St, Sslack, iter, success] = calc_v_y_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 admittance summation taken from:
Dragoslav Rajičić, Rubin Taleski, Two novel methods for radial and
weakly meshed network analysis, Electric Power Systems Research,
Volume 48, Issue 2, 15 December 1998, Pages 79-87
https://doi.org/10.1016/S0378-7796(98)00067-4