A method for automatically determining an impedance (Z x_real, Z x_imag) of a plurality of electrical components associated with multiple power phases, the method comprising: obtaining quadrature voltage values (Vxa, Vxβ) and quadrature current values (Ixα, Ixβ) associated with an electrical component of the plurality of electric components at a frequency of interest (ω base); using at least one processor (32, 172), computing at least one impedance value (Z x_real, Z x_imag) of the electrical component at the frequency of interest (ω base) based on the quadrature voltage and current values (Vxa, Vxβ, Ixα, Ixβ), characterized in that: computing the at least one impedance value (Z x_real, Z x_imag) of the electrical component at the frequency of interest (ω base) comprises: using the at least one processor (32, 172), computing a voltage magnitude value (Vmag) as a square root of a sum of squares of the quadrature voltage values (Vxa, Vxβ); using the at least one processor (32, 172), computing the voltage vector angle (θ) based on one of the quadrature voltage values (Vxa, Vxβ) and the voltage magnitude value (Vmag); using the at least one processor (32, 172), computing a current magnitude value (Imag) as a square root of a sum of squares of the quadrature current values (Ixα, Ixβ); using the at least one processor (32, 172), computing estimated d-axis and q-axis current component values (I e d, I e q) at least partially according to the quadrature current values (Ixα, Ixβ) and a voltage vector angle (θ); using the at least one processor (32, 172), computing an estimated impedance angle (δ) at least partially according to the estimated d-axis and q-axis current component values (I e d, I e q); and using the at least one processor (32, 172), computing (26, 27) the at least one impedance value (Z x_real, Z x_imag) of the electrical component at least partially according to magnitudes (Vmag, Imag) of the quadrature voltage and current values (Vxa, Vxβ, Ixα, Ixβ) and the estimated impedance angle (δ). obtaining quadrature voltage values (Vxa, Vxβ) and quadrature current values (Ixα, Ixβ) associated with an electrical component of the plurality of electric components at a frequency of interest (ω base); using at least one processor (32, 172), computing at least one impedance value (Z x_real, Z x_imag) of the electrical component at the frequency of interest (ω base) based on the quadrature voltage and current values (Vxa, Vxβ, Ixα, Ixβ), characterized in that: computing the at least one impedance value (Z x_real, Z x_imag) of the electrical component at the frequency of interest (ω base) comprises: using the at least one processor (32, 172), computing a voltage magnitude value (Vmag) as a square root of a sum of squares of the quadrature voltage values (Vxa, Vxβ); using the at least one processor (32, 172), computing the voltage vector angle (θ) based on one of the quadrature voltage values (Vxa, Vxβ) and the voltage magnitude value (Vmag); using the at least one processor (32, 172), computing a current magnitude value (Imag) as a square root of a sum of squares of the quadrature current values (Ixα, Ixβ); using the at least one processor (32, 172), computing estimated d-axis and q-axis current component values (I e d, I e q) at least partially according to the quadrature current values (Ixα, Ixβ) and a voltage vector angle (θ); using the at least one processor (32, 172), computing an estimated impedance angle (δ) at least partially according to the estimated d-axis and q-axis current component values (I e d, I e q); and using the at least one processor (32, 172), computing (26, 27) the at least one impedance value (Z x_real, Z x_imag) of the electrical component at least partially according to magnitudes (Vmag, Imag) of the quadrature voltage and current values (Vxa, Vxβ, Ixα, Ixβ) and the estimated impedance angle (δ). using the at least one processor (32, 172), computing a voltage magnitude value (Vmag) as a square root of a sum of squares of the quadrature voltage values (Vxa, Vxβ); using the at least one processor (32, 172), computing the voltage vector angle (θ) based on one of the quadrature voltage values (Vxa, Vxβ) and the voltage magnitude value (Vmag); using the at least one processor (32, 172), computing a current magnitude value (Imag) as a square root of a sum of squares of the quadrature current values (Ixα, Ixβ); using the at least one processor (32, 172), computing estimated d-axis and q-axis current component values (I e d, I e q) at least partially according to the quadrature current values (Ixα, Ixβ) and a voltage vector angle (θ); using the at least one processor (32, 172), computing an estimated impedance angle (δ) at least partially according to the estimated d-axis and q-axis current component values (I e d, I e q); and using the at least one processor (32, 172), computing (26, 27) the at least one impedance value (Z x_real, Z x_imag) of the electrical component at least partially according to magnitudes (Vmag, Imag) of the quadrature voltage and current values (Vxa, Vxβ, Ixα, Ixβ) and the estimated impedance angle (δ).