newtonpf_I_polar

newtonpf_I_polar(Ybus, Sbus, V0, ref, pv, pq, mpopt)

newtonpf_I_polar() - Solves power flow using full Newton’s method (current/cartesian).

[V, CONVERGED, I] = NEWTONPF_I_POLAR(YBUS, SBUS, V0, REF, PV, PQ, MPOPT)

Solves for bus voltages using a full Newton-Raphson method, using nodal
current balance equations and polar coordinate representation of
voltages, given the following inputs:
    YBUS  - full system admittance matrix (for all buses)
    SBUS  - handle to function that returns the complex bus power
            injection vector (for all buses), given the bus voltage
            magnitude vector (for all buses)
    V0    - initial vector of complex bus voltages
    REF   - bus index of reference bus (voltage ang reference & gen slack)
    PV    - vector of bus indices for PV buses
    PQ    - vector of bus indices for PQ buses
    MPOPT - (optional) MATPOWER option struct, used to set the
            termination tolerance, maximum number of iterations, and
            output options (see MPOPTION for details).

The bus voltage vector contains the set point for generator
(including ref bus) buses, and the reference angle of the swing
bus, as well as an initial guess for remaining magnitudes and
angles.

Returns the final complex voltages, a flag which indicates whether it
converged or not, and the number of iterations performed.

See also runpf(), newtonpf(), newtonpf_S_cart(), newtonpf_I_cart().