Budeanu’s Distortion Power Components Based on CPC Theory in Three-Phase Four-Wire Systems Supplied by Symmetrical Nonsinusoidal Voltage Waveforms

: Budeanu’s power theory, in its fundamental version, describes single-phase sinusoidal and nonsinusoidal systems. Over time, this elementary description has been extended to three-phase three-wire and four-wire systems, regardless of power conditions. Initially, three-phase systems were considered as three independent single-phase systems. A distinct approach was introduced by Czarnecki in his power theory (Currents’ Physical Components—CPC). The energy description and reference of the equivalent parameters of the load are comprehensive in the context of three-phase systems; Czarnecki treats such systems as a whole. This paper introduces a mathematical model to expand the basic Budeanu theory for three-phase four-wire (3-p 4-w) systems powered by symmetrical and nonsinusoidal voltage sources. The proposed approach is based on mutual elements between the fundamental Budeanu theory and the CPC theory, treating the 3-p 4-w system as a whole. In the extended Budeanu theory model, equations for the Budeanu reactive current and the Budeanu complemented reactive current are derived. The article also demonstrates their orthog-onality concerning the remaining components, indicating that each of the seven components can exist independently of the others. Furthermore, in the extended Budeanu theory, it is possible to identify which equivalent parameters of the load are responsible for the individual currents (powers) and which components are associated with the total distortion power proposed by Budeanu in 1927. All of the calculations were performed in Matlab/Simulink 2023b software.


Introduction
The exploration of energy properties in the frequency domain began in 1927 with C. I. Budeanu [1,2].This power theory, initially proposed by Budeanu, has undergone extensive development by numerous authors, evolving into increasingly sophisticated mathematical formulations [3][4][5][6][7].However, it was revealed that Budeanu's power theory misinterprets the physical phenomena occurring in systems, as demonstrated by L. S. Czarnecki in a 1987 publication [8].The conceptual error primarily concerns the interpretation of reactive power in systems powered by a nonsinusoidal but periodic voltage source.Budeanu defined reactive power in such systems as the sum of the amplitudes of the alternating components of the instantaneous power of the harmonic current and voltage.Additionally, distortion power (DB) was introduced as the sum of active power (P) and reactive power (QB) to apparent power (S), without establishing a clear connection to any physical phenomena in the system [2,[8][9][10].In the Budeanu theory, it is possible to calculate the active power of the entire system and the RMS value of the active current.
However, determining the reactance parameters of elements providing compensation and symmetry in systems supplied with sinusoidal or nonsinusoidal voltage, whether the voltage source is symmetrical or asymmetrical, was not feasible until the publication of papers [11,12].In the paper [11], a description of the current components of a single-phase load supplied with a nonsinusoidal voltage was presented with the equivalent parameters of the load.The publication [12] introduced the possibility of reactance compensation based on the extended Budeanu theory.
While the inaccuracy of the Budeanu theory was demonstrated in the publication [8], unfortunately, it continues to be utilized in universities and implemented in practical applications, such as in electricity meters or power quality parameter analyzers.
The development of the Budeanu theory, as presented by the authors, opens up opportunities for applications in various aspects of everyday life for both individual and industrial consumers.Despite the widespread adoption of LED lighting globally [40][41][42][43], it remains common practice to use the Budeanu theory when matching passive systems to compensate for reactive power, without considering the Budeanu distortion power.As highlighted in publications [11,12] and discussed in this article, a portion of the reactive power and the entire imbalance can be identified within the Budeanu distortion power.Renewable Energy Sources (RES) represent another domain where an accurate description of energy properties is crucial.Due to the possibility of single-phase and three-phase unbalanced systems, renewable sources contribute to increased load asymmetry [44][45][46][47][48] and the total harmonic distortion coefficient of voltage (THDU).Consequently, there is a growing need for a more precise depiction of the energy properties of such systems.The third rapidly developing area is wireless power transmission (WPT) [49][50][51][52][53][54].These systems utilize inductors for energy transfer and are primarily based on the power supply system available at a given location.This makes such systems sensitive to various types of asymmetry and distortion in both the supply voltage and the load itself.As is evident, a correct and consistent power theory can play a crucial role in diagnosing problems within a system accurately and cost-effectively enhancing its operational characteristics.
The article is divided into nine sections.Section 1 provides a literature review related to the Budeanu theory and the CPC theory, along with their practical applications in constructing passive and active compensators.Additionally, this section includes information about the potential application of the extended Budeanu theory in various practical scenarios, such as constructing compensators, lighting systems, wireless energy transmission, or implementing algorithms in electricity meters.Section 2 presents a mathematical model describing the CPC theory in three-phase four-wire systems, where the voltage source is symmetrical, and the waveforms generated by the source are distorted.Section 3 comprises a presentation of traditional equations used to calculate line currents in threephase systems (Ohm's law and Kirchhoff's law) and a description of the fundamental Budeanu theory.Section 4 focuses on presenting the mathematical model in the extended Budeanu theory for three-phase four-wire systems with nonsinusoidal and symmetrical power sources.This section also introduces the mutual elements in the CPC theory and the extended Budeanu theory.Section 5 discusses the second step in extending the Budeanu theory, which involves relating equivalent load parameters and the resulting currents' components.Section 6 demonstrates the orthogonality between individual components in the extended Budeanu theory.Section 7 provides a theoretical illustration, based on the mathematical model from Section 4, showcasing the steps involved in calculating the individual parameters of the load.Section 8 summarizes and discusses the previous sections, with a primary focus on Budeanu distortion power.Finally, Section 9 serves as a conclusion, summarizing the scientific results obtained and presenting potential future directions for the development or extension of the Budeanu theory.

Currents' Physical Components (CPC) Theory for 3-Phase 4-Wire Nonsinusoidal Symmetrical Systems
A voltage source in a power system can be a harmonics-generating source (see Figure 1).The vector of nonsinusoidal voltages supplying the LTI (Linear Time-invariant) load can be written as follows [21,22]: where n U -vector of complex RMS (CRMS) value of the voltage harmonic, and n - harmonic order.The vector of line currents can be written identically: where n I -vector of complex RMS (CRMS) value of the current harmonic.
From the active power P, the equivalent conductance of the full system can be defined, namely: ( ) The equivalent conductance determines the active current a i , the waveform of which is described by the equation: The three-phase RMS value of the active current Based on Equation (3) and the equivalent conductances for the individual harmonics, the scattered current s i is determined, the waveform of which takes the following form: ( ) The three-phase RMS value of the scattered current s i is ( ) By using the equivalent susceptances of the individual harmonics, the reactive current waveform r i can be defined: The three-phase RMS value of the reactive current r i is ( ) The overall unbalanced current u i can be expressed as ( ) The unbalanced current waveform Therefore, after the transformation, the load current i is expressed by the following equation: The three-phase RMS value of the load current i is The entire unbalanced current u i can be decomposed into three symmetrical com- ponents as follows: The decomposition of the unbalanced current into symmetrical components is based on the use of the equivalent admittance of the entire three-phase system: ( )  (18) where the generalized coefficient  is ( ) After the transformations, the current waveform of the unbalanced current of the positive sequence p u i is The three-phase RMS value of the unbalanced current of the positive sequence The waveform of the unbalanced current of the negative sequence The three-phase RMS value of the unbalanced current of the negative sequence The waveform of the unbalanced current of the zero sequence The three-phase RMS value of the unbalanced current of the zero sequence Finally, the load current can include six orthogonal components: The power equation takes the following form: where the apparent power S, calculated from the Buchholz apparent power formula [2,53], is The active power The reactive power The scattered power The unbalanced power of the positive sequence The unbalanced power of the negative sequence The unbalanced power of the zero sequence The power factor  is

Ohm's and Kirchhoff's Laws and Budeanu Theory
The line currents can be calculated based on Ohm's law: The approach to describing the energy properties of circuits comes from the description of power (initially for single-phase systems).The theory is also used to describe threephase systems.
The active power P in an individual phase of a three-phase system can be written as follows: where X denotes the successive phase from the set {R, S, T}.Identically, in Budeanu's approach, the reactive power in an individual phase B Q can be defined as "The response" to the difference between the apparent power S and the active power P and reactive power QB is the distortion power DB, which is expressed by the equation: The description of the energy properties and power components in the original Budeanu theory stops with the relationships noted above.

Development of the Budeanu Theory for 3-Phase 4-Wire Nonsinusoidal Symmetrical Systems-Step I
Noting that the active power P in the CPC theory and the Budeanu theory are identical, the equivalent conductance of the overall system can be calculated (3), and thus: ( ) From this, it is possible to determine the active current in the Budeanu theory, which will have an identical waveform and an identical three-phase RMS value as in the CPC theory: From the Budeanu reactive power relation (41), the equivalent susceptance eB B of the whole system can be found, namely: Based on this, the Budeanu reactive current rB i can be defined, the waveform of which takes the form: The three-phase RMS value of the Budeanu reactive current It should be noted that the reactive current in CPC theory and the Budeanu reactive current are currents that are different* with respect to three-phase RMS value and waveform, namely: * is the situation in which both currents are identical.The reactive currents in the two approaches are identical if and only if the susceptance does not change with the harmonic order for both the capacitive and inductive nature of the load -condition in the equation (80).
To the first step of the extension of the Budeanu theory, the waveform of the distortion current must be added.However, it is not determined from the equivalent parameters of the load but only by the difference in the individual currents: ( ) The three-phase RMS value of the Budeanu distortion current dB Point 2: The scattered current from the CPC theory is kept and defined as in Equations ( 6) and (7).
Point 3: The unbalanced current of the positive sequence in the CPC theory remains and is defined as in Equations ( 20) and (21).
Point 4: The unbalanced current of the negative sequence in the CPC theory remains and is defined as in Equations ( 23) and (23).
Point 5: The unbalanced current of the zero sequence in the CPC theory remains and is defined as in Equations ( 24) and (25).
Point 6: Look more closely at the Budeanu reactive current-relations ( 45) and (46).The value of the Budeanu reactive current (waveform and three-phase RMS value) is less than the value of the reactive current in the CPC theory.In order to find the missing value of the current in the Budeanu theory and to relate it to the equivalent parameter of the load (to give a physical interpretation), it is necessary to look at the equation for the scattered current in the CPC theory, which is the difference between the conductances of the individual harmonics and the equivalent conductance of the whole system, since From the susceptance of the individual harmonics and the equivalent susceptance of the entire system, we can determine the waveform of the Budeanu complemented reactive current crB i , which is ( ) The three-phase RMS value of the Budeanu complemented reactive current crB i is ( ) In summary, the load current described by the extended Budeanu theory (a combination of the Budeanu theory and the Currents' Physical Components theory), can have seven components: These components are orthogonal; therefore: The power equation takes the following form: where the apparent power S, calculated from the Buchholz apparent power relation, is The active power The Budeanu reactive power The Budeanu complemented reactive power The unbalanced power of the zero sequence The power factor  is

Orthogonality of the Components in the Extended Budeanu Theory
Vectors are mutually orthogonal when the scalar product of the three-phase quantities [2,21,22] is zero: If the three-phase quantities are represented in the form of complex RMS values, namely: then their scalar product is equal to All components in the original CPC theory are mutually orthogonal [21].The question of the orthogonality of the individual components with respect to the Budeanu reactive current and the Budeanu complemented reactive current remains.( ) 10. Orthogonality of the Budeanu complemented reactive current and the unbalanced current of the negative sequence: 11. Orthogonality of the Budeanu complemented reactive current and the unbalanced current of the zero sequence:

Theoretical Illustration
The system in Figure 2 is supplied from a nonsinusoidal symmetrical voltage source.The distorted voltage source, apart from the fundamental harmonic (230 V, 50 Hz), also generates harmonics of higher orders, i.e., n = 3, 5, 7 (3rd-15 V, 5th-25 V, 7th-10 V).The three-phase load is formed by linear elements in the form of resistance, inductive reactance and capacitive reactance, which were selected randomly, albeit with the load asymmetry.The system scheme is shown in Figure 2. The resistance, inductive reactance and capacitive reactance values for the fundamental harmonic are listed in Table 1.The waveform of the nonsinusoidal voltage at the R, S, T terminals is shown in Figure 3.The three-phase RMS value of the nonsinusoidal supply voltage shown in Figure 3 is presents the complex values of the admittances in the individual lines for each harmonic.

Harmonic Order Line R in [S] Line S in [S] Line T in [S]
n = 1 .e− Based on Ohm's law and Kirchhoff's law, the complex values of the harmonic currents in the individual lines were calculated.Table 3 presents the resulting complex values of the nonsinusoidal line currents.

Harmonic Order Line R in [A] Line S in [A] Line T in [A]
n = 1   The apparent power (28), obtained by multiplying the three-phase RMS value of the supply voltage [2] by the three-phase RMS value of the load current (39), is  5 presents the equivalent parameters in the form of equivalent conductance and susceptance of the individual harmonics of a three-phase load.15)-( 18) described in CPC theory.Based on relations ( 5), ( 7), ( 9), ( 21), (23), and ( 25), the three-phase RMS values of the currents' components described by the CPC theory were calculated.The three-phase RMS values of the load currents' components for the individual harmonics described by the CPC theory are presented in Table 7. Table 8 shows the complex values of the currents' components and the total current of a three-phase load for the individual harmonics following the CPC theory..e Based on Table 8 and relations ( 4), ( 6), ( 8), ( 20), (22), and ( 24), Figures 5-10 show the waveforms of the individual load currents' components described in the CPC theory.In a system with an asymmetrical structure of a three-phase linear load, there is a current flow in the neutral line, the RMS value of which is associated with the three-phase RMS value of the unbalanced current of the zero sequence multiplied by the value of 3 , or use relation (37), and is equal to Summing the values of the waveforms of the currents' components shown in Figures 5-10, we find the value of the waveform of the line currents based on the CPC theory.This waveform of line currents is presented in Figure 11.5), ( 7), ( 21), ( 23), ( 25), (47), and (53), the three-phase RMS values of the currents' components described in the extended Budeanu theory were calculated.The three-phase RMS values of the load currents' components for the individual harmonics defined in the extended Budeanu theory are presented in Table 9. .e Based on Table 10 and relations ( 46) and ( 52), the waveforms of the two currents' components of the three-phase load described by the extended Budeanu theory are shown in Figures 12 and 13.The waveforms of the active currents (Figure 5), the scattered currents (Figure 6), the three unbalanced currents of the positive sequence (Figure 7), the negative sequence (Figure 8) and the zero sequence (Figure 9) are equal in both approaches.

Analysis and Discussion of the Obtained Results
What remains to be solved are the components of the distortion current and thus the distortion power described by Budeanu.In three-phase four-wire systems supplied from a symmetrical nonsinusoidal voltage source and a three-phase asymmetrical load formed by resistance, inductance and capacitance in any configuration in each line, the distortion power defined by Budeanu is expressed by five components, namely: the scattered power, the Budeanu complemented reactive power and the three unbalanced powers of the positive, negative and zero sequence.Similarly, the five currents' components relate to the Budeanu distortion current.
Figure 15 shows the conjunction of currents' components in the extended Budeanu theory and the CPC theory.Figure 16, on the other hand, shows an identical conjunction of the components associated with the Budeanu distortion power.In accordance with the data in Table 11, the Budeanu distortion current waveform described by the extended Budeanu theory is shown in Figure 15.
In the analyzed example, the Budeanu distortion current, the waveform of which is shown in Figure 17, consists mainly of the zero sequence unbalanced current u z I .Hence, the typical 120° shift of the individual phases is not observed.It should also be noted that the load in Figure 2, by suitably matching the imbalance in the elements in the individual lines, results in the Budeanu distortion power consisting of 5 components.At this point, it should be mentioned that the components of the Budeanu distortion power do not necessarily have to be five.It is possible to find situations in which the Budeanu distortion power consists of one, two, three or four powers' components.
As presented in the theoretical illustration, the understanding of distortion current allows for its interpretation based on the number of components.If the distortion current has only one component, it indicates that the three-phase load is symmetrical, consisting solely of L or C elements (it is worth mentioning a case that is not necessarily unrealistic but deviates from theoretical assumptions, namely a load composed of resistive elements that change with the frequency of harmonics).When there are two components in the distortion current, the three-phase load is symmetrical and consists of RC or RL elements.The distortion current also has only two components when the load is composed of reactive elements in the form of an inductor L or a capacitor C, whose parameters do not change with the harmonic order (denying an increase or decrease in susceptance with an increase in frequency).In this case, specific harmonics occur, namely for harmonics with the positive sequence (n = 3k + 1), the two components of distortion current are the unbalanced currents of negative sequence and zero sequence.For harmonics with the negative sequence (n = 3k − 1), the two components of the distortion current are the unbalanced currents of positive sequence and zero sequence, and for harmonics with zero sequence (n = 3k), the two components are the unbalanced currents of positive sequence and negative sequence.A distortion current composed of three components exists when the threephase load is asymmetrical and consists of reactive elements in the form of an inductor L or a capacitor C, with parameters unchanged with the harmonic order, and at least two distinct sequence harmonics.The scenario with three components in the distortion current also occurs when the three-phase load is asymmetrical and composed of R, L, or C elements, with only specific harmonics present.Specifically, for harmonics with the positive sequence (n = 3k + 1), the three components of distortion current are the unbalanced currents of negative sequence, zero sequence and the Budeanu complemented reactive current.For harmonics with the negative sequence (n = 3k − 1), the three components are the unbalanced currents of positive sequence, zero sequence and Budeanu complemented reactive current, and for harmonics with zero sequence (n = 3k), the three components are the unbalanced currents of positive sequence, negative sequence and Budeanu complemented reactive current.The case of four components in the distortion current occurs when the load is asymmetrical and consists of R, L, C, or variable R for at least two distinct sequence harmonics.The four components are also present when the three-phase load is asymmetrical and composed of RC or RL elements, with only specific harmonics occurring.Specifically, for harmonics with the positive sequence (n = 3k + 1), the four components of distortion current are the unbalanced currents of negative sequence, zero sequence, the Budeanu complemented reactive current and scattered current.For harmonics with the negative sequence (n = 3k − 1), the four components are the unbalanced currents of positive sequence zero sequence, Budeanu complemented reactive current and scattered current and for harmonics with zero sequence (n = 3k), the four components are the unbalanced currents of positive sequence, negative sequence, Budeanu complemented reactive current and scattered current.
As can be observed, Budeanu distortion current (distortion power) occurs when the load admittance varies with the harmonic order, i.e., As observed from the analysis of specific cases describing distortion current components, this is possible due to the extension of Budeanu's theory with the inclusion of Budeanu complemented reactive current described in this article for three-phase four-wire systems powered by a nonsinusoidal but symmetrical voltage source.Further development of the extended Budeanu theory should consider the possibility of constructing and synthesizing passive compensation systems and expanding the mathematical model for systems with asymmetrical power supply.

Conclusions
As presented in this paper, it is possible to develop Budeanu's theory to a form in which the currents' components of a three-phase four-wire load supplied from a symmetrical and nonsinusoidal voltage source can be defined.With the development of Budeanu's theory, it is possible to accurately characterize the Budeanu distortion current by the components of the load current and not by the difference in complex RMS values (CRMS) vectors.This knowledge makes it possible to recognize the nature of the distortion power in the analyzed cases, which is not only related to the imbalance of the load but also to the variation of active or reactive powers.
The mathematical model in the extended Budeanu theory describing the equivalent parameters of a three-phase four-wire load is capable of distinguishing the orthogonal components of the current, and this makes it possible to determine the parameters of a passive compensator both with ideal parameters and a compensator minimizing reactive current and unbalanced current, which contributes to reducing the currents' components associated with the reactance elements of the load (reactive powers and unbalanced powers).
on Kirchhoff's first law, the current in the neutral line is three-phase RMS value of the load current is

) 5 .
Development of the Budeanu Theory for 3-Phase 4-Wire Nonsinusoidal Symmetrical Systems-Step II Assumptions in developing the Budeanu theory for three-phase four-wire systems supplied with nonsinusoidal symmetrical voltage: Point 1: The active currents in the two theories are equal to each other and are represented by identical equivalent load parameters-the equivalent conductance of the load.
− -describes the scattered current, so eB bn BB − -describes the missing value of the Budeanu reactive current(51)

1 . 2 . 3 . 4 . 5 .
Orthogonality of the Budeanu reactive current and the active current: Orthogonality of the Budeanu reactive current and the scattered current: Orthogonality of the Budeanu complemented reactive current and the active current: Orthogonality of the Budeanu complemented reactive current and the scattered current: Orthogonality of the Budeanu complemented reactive current and the Budeanu reactive current:

73) 6 . 7 . 8 .
Orthogonality of the Budeanu reactive current and the unbalanced current of the positive sequence: Orthogonality of the Budeanu reactive current and the unbalanced current of the negative sequence: Orthogonality of the Budeanu reactive current and the unbalanced current of the zero sequence:

9 .
Orthogonality of the Budeanu complemented reactive current and the unbalanced current of the positive sequence:

Figure 2 .
Figure 2. Scheme of a three-phase four-wire system supplied with a nonsinusoidal symmetrical voltage with an asymmetrical linear load.

Figure 3 .
Figure 3.The waveform of nonsinusoidal voltage values at the terminals of a three-phase linear load.

Figure 4
Figure4shows the waveform of line currents resulting from supplying a three-phase load with a nonsinusoidal symmetrical voltage.

Figure 4 .i
Figure 4. Waveform of the line currents obtained from Ohm's law and Kirchhoff's law.The three-phase RMS value of the nonsinusoidal current of the load shown in Figure 4 is 87 336 A .= i Table 4 presents the three-phase active powers (40) and reactive powers (41) of the individual harmonics and the powers of the entire system.
Based on expression (3), the equivalent conductance of the system in Figure2 is

Figure 5 .
Figure 5.The waveform of active currents based on the CPC theory.

Figure 6 .
Figure 6.The waveform of scattered currents based on the CPC theory.

Figure 7 .
Figure 7.The waveform of reactive currents based on the CPC theory.

Figure 8 .
Figure 8.The waveform of unbalanced currents of the positive sequence based on the CPC theory.

Figure 9 .
Figure 9.The waveform of unbalanced currents of the negative sequence based on the CPC theory.

Figure 10 .
Figure 10.The waveform of unbalanced currents of the zero sequence based on the CPC theory.

Figure 11 .
Figure 11.The waveform of the line currents based on the CPC theory.As can be noted, the waveforms in Figures 4 and 11 are similar.According to the extended Budeanu theory and relation (45), the equivalent susceptance of the system from Figure 2 is are similar.According to the extended Budeanu theory and relation (45), the equivalent susceptance of the system from Figure 2 is e 0 1011 S B. = Based on relations (

Figure 12 .
Figure 12.The waveform of the Budeanu reactive currents based on the extended Budeanu theory.

Figure 13 .
Figure 13.The waveform of the Budeanu complemented reactive currents based on the extended Budeanu theory.By summing the values of the currents' components shown in Figures 5, 6, 8-10, 12 and 13, we obtain the value of the waveform of the line currents based on the extended Budeanu theory.This waveform of line currents is presented in Figure 14.
, 8-10, 12 and 13, we obtain the value of the waveform of the line currents based on the extended Budeanu theory.This waveform of line currents is presented in Figure14.

Figure 14 .
Figure 14.The waveform of the line currents based on the extended Budeanu theory.As can be seen, the waveforms in Figures 4, 11 and 14 are exactly the same.

Figure 15 .
Figure 15.Connection of the 3-p 4-w currents' components of the load in the extended Budeanu and the CPC theories.

Figure 16 .
Figure 16.Connection of the 3-p 4-w powers' components of the load in the extended Budeanu and the CPC theories.Table 11 presents the complex values of the currents' components and the Budeanu distortion current for the individual harmonics based on the extended Budeanu theory.

Figure 17 .
Figure 17.The waveform of the Budeanu distortion currents based on the extended Budeanu theory.
where n and h denote the orders of harmonics in the set can be presented in the form:

Table 1 .
List of resistance, inductive reactance and capacitive reactance values for the fundamental harmonic.

Table 4 .
List of three-phase active and reactive powers values of individual harmonics and the total three-phase active and reactive powers of the entire system.

Table 5 .
List of conductance and susceptance values of individual harmonics.

Table 6
presents the unbalanced admittances and equivalent admittance values obtained from relations (

Table 6 .
List of unbalanced admittances and equivalent admittance values in CPC theory.

Table 7 .
List of three-phase RMS values of currents' components for individual harmonics described by the CPC theory.

Table 8 .
List of values of complex currents' components and total current for individual harmonics based on the CPC theory.

Table 9 .
List of three-phase RMS values of the currents' components for the individual harmonics described by the extended Budeanu theory.

Table 10
presents the complex values of the currents' components and the total current of a three-phase load for the individual harmonics based on the extended Budeanu theory.

Table 10 .
List of values of complex currents' components and total current for individual harmonics based on the extended Budeanu theory.

Table 11 .
List of complex values of currents' components and distortion current for individual harmonics based on the extended Budeanu theory.