{'fraction_non_alphanumeric': 0.0347655384225691, 'fraction_numerical': 0.01212677706361017, 'mean_word_length': 4.516855895196507, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 13, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Determination of Inactive Powers in a Single-Phase AC Network

: Based on the development of the theory of reactive power and distortion power, starting with the works of Fryze and Budeanu, it has been found that the contradictions in the deﬁnition of the components of inactive powers are caused by errors in the introduced intermediate concepts and corresponding calculations when switching to nonlinear and non-sinusoidal AC circuits. The materials of the works of modern researchers and the numerical calculations carried out made it possible to trace the differences between reactive power and distortion power, to conﬁrm the orthogonality properties of the active, reactive power, and distortion power components. The paper deﬁnes the conditions for achieving a power balance in an AC network with nonlinear loads, compiled and tested criteria leading to the absence of distortion power in a single-phase AC network. Using the time base of the projection of the generalized vectors in vector diagrams, it is shown that compliance with the criteria for the absence of distortion power does not determine the mutual similarity of the voltage curve with the current curve for a nonlinear load. It has been found that the well-known term “distortion power” has an unfortunate wording, since this power, although it characterizes the interaction of harmonics of currents and voltages with different ordinal numbers, is not determined by the visual similarity or the degree of distortion of the load current waveforms relative to the supply voltage curve.


Introduction
At the initial stage of the power theory development, known since the beginning of the 20th century for sinusoidal AC networks, there were no contradictions in the formulations and calculated expressions [1,2]. Instead of instantaneous values of sinusoidal voltages u(t) and currents i(t), for the convenience of analysis and calculations in the well-known symbolic method, it is proposed to use their effective values U and I. The product of these values U · I determines the total or apparent power S = U · I. The reactive nature of the network z elements explains the phase shift φ between voltage u(t) = √ 2U 1 sin(ωt) and current i(t) = √ 2I 1 sin(ωt + φ). In this case, the apparent power S consists not only of active power P = U I cos(φ) capable of performing work, but also of ballast reactive power Q = U I sin(φ). The components P and Q make up the power balance S = P 2 + Q 2 , in which the proved orthogonality of P with respect to Q [1,2] is emphasized.
However, the appearance of loads with nonlinear current-voltage characteristics (CVC) led to a significant distortion of the network currents' curve shape and of supply voltages, which resulted in the formation of other inactive power components that do not perform useful work. It is obvious that the interaction of harmonics of non-sinusoidal voltage and current will create its own components of reactive power, etc. from each harmonic separately.
However, in this case, the excess of power S over power P can no longer be explained by the emergence of only reactive power Q ∑ because of all harmonics. In the scientific literature [3,4], this circumstance is explained by the appearance of another component of inactive power, which is known as the power of D distortions. However, such view is not supported by all researchers and, to this day, it causes controversy about the correct definition of the combination of Q ∑ and D, and their delineation taking into account higher harmonics.
For about 100 years, discussions have continued on the definition of the terminology for components Q ∑ and D, their calculation, and the correct interpretation of their physical nature [4][5][6][7].

Formulation of the Problem
At present, AC electrical networks do not have strictly sinusoidal voltages, and most loads have sharply nonlinear I-V characteristics, which in turn causes non-sinusoidal currents to flow through all elements of electrical networks. Therefore, an important task is to study the processes of formation of individual power components due to the action of individual harmonics of voltages and currents and to illustrate their calculations using examples of simple loads in a single-phase alternating current network.
The practical value of the research results in this article is determined by its contribution to solving the problem of calculating inactive powers in the AC network. A clear understanding of the processes of formation of components of inactive powers will make it possible to clarify their formulations and definitions from the standpoint of physical interpretation and draw up appropriate expressions for the correct calculation of the values of reactive powers and distortion powers in a non-sinusoidal AC network.
This will allow one to correctly calculate each power component in the AC network, to create new, more advanced devices for analyzing and metering electrical energy, and to successfully solve urgent energy optimization problems of redistribution and compensation of inactive power.
In addition, the proven statements of the theory of power for single-phase alternating current networks can be subsequently applied to three-phase alternating current networks, where the occurrence of asymmetric modes significantly complicates the analysis if there is no clear representation of inactive and full power [8,9].

Aspects in the Development of Inactive Capacities Theories
In the historical aspect, the works of many scientists are devoted to the definition of inactive powers in non-sinusoidal AC networks. The work of Budeanu [3] describes one of the first explanations of the mechanism of occurrence of the power imbalance, when the power is decomposed into orthogonal components P and Q, and defining Q ∑ as a set of reactive powers Q i caused by individual harmonics: where i is the ordinal numbers of the harmonics of the effective values of voltages U i and currents I i ; φ i is the shear angles between them. It is in this theory that the power balance is achieved by introducing an additional component of inactive power, mutually orthogonal with respect to the already known components and in such a way that: The value D was presented as a measure of the mutual distortion of the curves u(t) and i(t). However, it has already been noted in [10] that, in AC networks with certain combinations of higher harmonics at the level of instantaneous powers, the interaction of the components Q and D can occur and in such a way that D = 0 even in conditions of a distorted shape of the curve i(t) relative to u(t). Vice versa, D = 0 provided that the shape of the curve i(t) is similar to the shape of the curve u(t). This fact gave rise to doubts about the validity of the Budeanu theory and the search for other approaches to determining the inactive power that do not contradict the results of the experiments.
Such approaches were available in early works by Fryze [2,10,11]. In them, inactive power Q q is not divided into components Q and D or the like, but it is defined as a power imbalance: Since, by definition, Q q cannot do the job, the Fryze theory is built on dividing the instantaneous load currents into two components: active i a (t) and reactive i r (t) currents, in such a way that i(t) = i a (t) + i r (t). Assuming that the currents' components i a (t) and i r (t) are mutually orthogonal, their respective effective values I 2 a and I 2 r will define the initial current as I 2 = I 2 a + I 2 r . Hence, there is no explicit physical interpretation for power Q q in the representation of the Fryze theory, as it has been done for sinusoidal circuits.
Subsequent works [12,13], in which the development of the theory of powers is based on frequency techniques, are devoted to the establishment of the relationship between the higher harmonics of voltage and current with the components of inactive powers. This defines reduced power, and, using the ratios of the frequencies of the fundamental and higher harmonics with ordinal numbers k, the following determination of inactive power Q has been obtained: However, it, like many other similar calculated expressions in practice, had not always provided a power balance [10,14] in an AC network with a non-sinusoidal shape u(t) and i(t).
In subsequent works [15], an attempt was made to further divide the current i r (t) into orthogonal components: inactive current and a residual component. The inactive current determines the parameter of the passive reactive element (inductance or capacitance). For the harmonics with ordinal numbers n, this approach provides an expression for determining Q: The value Q in these works was obtained as a result of reducing the reactive nature of the network elements, and, from a practical point of view, it allows for determining the parameters of the compensation device for the inactive current component using passive elements. However, the origin and physical meaning of the residual component of the inactive current are not considered.
The component of current i a (t), according to the Fryze theory, was determined in subsequent works by Czarnecki [7,16] in accordance with expression (1) and the development of the provisions of the CPT (Conservative Power Theory). The inactive part was analyzed from the point of view of frequency calculation methods according to the CPC (Currents Physical Components) theory that are also using the methodology of orthogonal decomposition [17]. In this representation, reactive power is defined as the result of a large number of orthogonal components of currents I i of different frequencies and is calculated as: The total power of a linear load supplied from an AC network with a non-sinusoidal u(t), in accordance with this theory, is S = P 2 + D S 2 + Q Z 2 . The component D S has been named as the dissipated power generated by the current components depending on the frequency spectrum u(t). Therefore, only the component Q Z can be compensated by means of passive elements [7]. If the load turns out to be nonlinear, then this theory allows one to single out another component of the current i g (t), the harmonics of which are not determined by the spectrum of the periodic function, but are independently generated by the nature of the load itself [16], forming the corresponding power component D g . In this case, the total power in a single-phase network will be determined by the combination of Thus, in studies using frequency, spectral, and similar calculation methods based on the concepts of effective values of voltages and currents, no unity was found in the quantitative determination of the components of inactive powers and their physical interpretation.
In order to exclude possible errors in determining the power components, made in the course of intermediate calculations based on instantaneous values u(t) and i(t), neither frequency nor spectral, but integral techniques were used in the works of individual researchers, which do not require finding the effective values of harmonics. Thus, in the methods [13,18], the nature of the load is determined through its I-V characteristic. For this purpose, the area of the closed figure of the I-V characteristic of the load in the axes u(t) and i(t) for the period T = 2π has been determined by integration as a value [19] that is proportional to the power Q: A similar approach is described in [13], where the value Q was obtained in the course of the cyclic integration operation: It is easy to see that expressions (2)-(8) not only characterize the components of inactive capacities in different ways, but, in the case of calculations, they will lead to a significant discrepancy between the obtained results.
Therefore, in recent years, in the theory of inactive powers, the tendency of direct analysis of instantaneous values has dominated, which has become the basis of the theory of instantaneous powers. The advantage of the theory of instantaneous powers is the elimination of intermediate concepts and calculations, which, when applied to nonlinear and non-sinusoidal AC circuits, can lead to fundamental errors in definitions and calculations due to an incorrect interpretation of the physical nature of the emergence of individual power components.
The concept of considering the instantaneous values of voltages and currents, applied in the first works of Fryze [2,11], was later used in [20] as the FBD theory (Fryze, Buchholtz, Depenbrock), described in more detail in [21,22]. The theory of instantaneous powers was further developed in the concept of the pq-theory [23], which determines the instantaneous values of active p(t) and reactive q(t) power through combinations of the products of the values of projections u(t) and i(t) onto orthogonal axes. The consequence of the development of the theory of instantaneous powers was its generalization in the crossvector theory [17,24], in which the value p(t) is defined as the scalar product of vectors → u (t) and → i (t), and the value q(t) is defined as their vector product. It should be noted that the provisions of individual theories of inactive power are capable of solving a number of practical problems with some success, for example, analyzing the operating modes of electrical networks or designing technical means to compensate for reactive power [25]. However, until now, none of the presented theories have yielded an unambiguous physical interpretation or explanation of the physical meanings of all the components of inactive powers, their occurrence, and interaction for arbitrary shapes of curves u(t) and i(t).
This hinders the development of effective calculation methods and the adoption of relevant standards, in which, first of all, there would be a clear, relatively simple, and unambiguous definition of the process of occurrence of each power component for the subsequent calculation of the total inactive power in electrical networks in which nonsinusoidal voltages and currents operate. This is evidenced by the IEEE 1459-2010 Standard adopted at due time [26], in which the concept of inactive power is introduced, but the definition of power Q is limited only by the first harmonic.

Analysis of Inactive Power Components in Wave and Vector Diagrams
If voltage and current harmonics, for example, the first harmonics i = j = 1, with the same ordinal numbers interact in the AC network, then, from the point of view of the symbolic method [1], the current vector I 1 can either lag behind or ahead of the supply voltage vector U 1 by an angle φ. Characterizing these values as u(t) = U 1m sin(wt) and i(t) = I 1m sin(wt − φ), the instantaneous load power p(t) = u(t) · i(t) is conventionally divided into a constant component p = and an alternating component p ∼ : The integral sum of the components p = and p ∼ for the period T corresponding to the full reversal of the vectors I 1 and U 1 is capable of forming active power P 1 , being determined by the argument φ 1 as P 1 = 1 : In this connection, 1 The appearance of an alternating value p(t) in the combination of components p = and p ∼ over a period T characterizes the action of exchange or reactive power Q between the elements of the AC network, for example, a source and a load. As it follows from expression (9), the sign p(t) will change at a frequency of 2wt twice over the period T.
The quantity Q 1 can be traditionally expressed through the unbalance of total power S 1 = U 1 I 1 and active power P 1 , using the well-known property of orthogonality of the components Q and P as in the case of Q 2 = S 2 − P 2 , which gives the well-known expression [1,3]: As the vector I 1 lags behind the vector U 1 , the component p = decreases, and negative values p(t) appear and grow, which indicates a decrease in the component P 1 and the appearance of a growing component Q 1 .
If the harmonics of voltage U i and current I j , having different ordinal numbers i = j, interact in an AC network, then, due to the orthogonality property [27] of these harmonics with multiple frequencies j/i = n during the working period T, the work cannot be done, since, in this case, the power is: It makes no sense to further decompose such power into components p = and p ∼ . However, as it follows from expression (3), the total power in this case will be a value S ij = U i I j = 0 that is commonly called the distortion power S ij = D ij .
In this way, in an alternating current network, even with a sinusoidal voltage source, a load with a purely active character is capable of generating powers Q i and D ij , if it has a nonlinear I-V characteristic.
It is common knowledge that, if an artificial delay occurs in the opening of an electronic key in a single-phase AC network supplying a purely active load, a spectrum of higher current harmonics will form [5][6][7]. In this spectrum, only the first harmonic of the current will create a component Q 1 .
As an example, Figure 1 shows oscillograms and the spectral composition (Table 1) of the load current with R = 1 p.u. in a network with a sinusoidal voltage source of U = 100 p.u. and the frequency of f = 50 Hz when the delay in the opening of the electronic keys was about 70% or 0.007 s.
In such network with a sinusoidal voltage of the power supply, a load with a purely active character causes a phase shift in the fundamental current harmonic of −54.8 • , as if the load electrical circuit contained an inductive element, the reactive power of which can be formally calculated by expression (11). However, it is, of course, erroneous to assume that such circuits contain reactive elements that exchange power with a power source. Moreover, in separate works [15,17], it was noted that the concept of reactive power, which "characterizes the energy pumped from the source to the reactive elements of the receiver and then returned by these elements back to the source during the same period of oscillation, referred to this period", does not fully reflect the physical meaning. In such network with a sinusoidal voltage of the power supply, a load with a purely active character causes a phase shift in the fundamental current harmonic of −54.8°, as if the load electrical circuit contained an inductive element, the reactive power of which can be formally calculated by expression (11). However, it is, of course, erroneous to assume that such circuits contain reactive elements that exchange power with a power source. Moreover, in separate works [15,17], it was noted that the concept of reactive power, which "characterizes the energy pumped from the source to the reactive elements of the receiver and then returned by these elements back to the source during the same period of oscillation, referred to this period", does not fully reflect the physical meaning.  If the action of each harmonic (Table 1) is analyzed separately in this example, then it is not too difficult to see in Figure 2a that the vector of the current third harmonic 3 I , built according to the data in Table 1   If the action of each harmonic (Table 1) is analyzed separately in this example, then it is not too difficult to see in Figure 2a that the vector of the current third harmonic I 3 , built according to the data in Table 1, is three times ahead of the voltage vector per one period T.
The interaction of vectors U 1 and I 3 will create instantaneous power p 3 (t), the alternating pulsations of which (Figure 2b), due to the properties of expression (12), are distributed in such a way that the total area bounded by the curve p 3 (t) with positive values is equal to the total area bounded by this curve with negative values so that P 13 = 0. At the same time, the total power in this case is S 13 = U 1 I 3 = D 13 = 711.8 V·A. In Figure 3 waveforms explaining the power generation for the fifth (Figure 3a) and seventh (Figure 3b) harmonics are shown. The remaining higher harmonics of current interacts with voltage U 1 in a similar way; this is accompanied by the appearance of powers P 15 = 0, D 15 = 664 V·A, P 17 = 0, D 17 = 358 V·A.  In fact, power 1 Q interacts with power components ij D at the level of redistribution of instantaneous values p(t) , while the values of the components 1 Q and ij D remain unchanged [28]. This is due to the fact that one harmonic, characterizing a certain power component, as follows from the properties of the Fourier series [27], cannot be obtained or compensated for by a set of other harmonics. In this case, it is possible to correct the shape of the curve p(t) of an arbitrary AC consumer by regulating the phase shifts  In fact, power 1 Q interacts with power components ij D at the level of redistribution of instantaneous values p(t) , while the values of the components 1 Q and ij D remain unchanged [28]. This is due to the fact that one harmonic, characterizing a certain power component, as follows from the properties of the Fourier series [27], cannot be obtained or compensated for by a set of other harmonics. In this case, it is possible to correct the shape of the curve p(t) of an arbitrary AC consumer by regulating the phase shifts Thus, the actions of the instantaneous values of the components of p i (t), taking into account all higher harmonics, are completely mutually neutralized in the time intervals when the electronic switch is open (Figure 1b). This example shows that the power Q 1 = U 1 I 1 sin(φ 1 ) from the first harmonic of the current is not the prerogative of the manifestation of the properties of reactive elements in the AC network, but it can be considered as an imbalance of instantaneous values p i (t) that have caused a shift in the components p = and p ∼ .
In fact, power Q 1 interacts with power components D ij at the level of redistribution of instantaneous values p(t), while the values of the components Q 1 and D ij remain unchanged [28]. This is due to the fact that one harmonic, characterizing a certain power component, as follows from the properties of the Fourier series [27], cannot be obtained or compensated for by a set of other harmonics. In this case, it is possible to correct the shape of the curve p(t) of an arbitrary AC consumer by regulating the phase shifts of the quantities p i (t) created, for example, in accordance with expression (12) using an inactive power compensator with zero active power [29][30][31].
The set of values with orthogonal properties of harmonics is written as: In this example, there is only one voltage harmonic U 1 ; therefore,  (10) and (11) have the following values: This provides the power balance according to expression (2), since: This yields S = P 2 + Q 2 + D 2 or 1825.46 ≈ √ 712 2 + 1005 2 + 1347.4 2 . This confirms the idea in the Budeanu theory [3] that the powers P, Q, and D are mutually orthogonal, and the difference √ S 2 − P 2 determines the totality of inactive powers Q 2 + D 2 .

Analysis of Budeanu Theory Criticisms
Despite the fact that calculations similar to the above-mentioned example prove the consistency of the Budeanu's theory [3], a number of researchers, for example [4,9,10,32], pose a number of questions: (1) Why, in practice, expression (2) is not true for all measured modes of operation of consumers in the alternating current network; (2) Why distortion power D can occur in conditions when the load current curve is similar to the supply voltage curve and vice versa; D = 0 provided that the load current curve is clearly distorted relative to the supply voltage curve; The answer to the first question is determined, first of all, by the reasons that violate the property of harmonics orthogonality (12). This is possible if the parameters of the higher harmonics of the current (amplitude, phase shift) are not stationary during the operating period T, which is quite typical of AC circuits containing modern semiconductor devices. This will lead to a violation of the equality of the areas bounded by the curves p i (t) for positive and negative values, as is shown in Figure 2b. In such cases, individual current harmonics with non-stationary parameters, interacting with the mains voltage U 1 , are able to do work because, for them, P i = 0.
For example, let us assume a simple condition that per one period: one can obtain a pulsating, but positive power in Figure 2c. It is not too difficult to see that, as a result of the accepted condition, a harmonic component i 1 (t) will appear in the spectral composition of the curve i 3 (t) shown in Figure 2c, the lagging character of which will create an additional Q 1 = 31.8 V·A component relatively u 1 (t).
As a result, the active power calculated based on the first and highest harmonics will be: reactive power: distortion power: Violation of the power balance P 2 + Q 2 + D 2 = S or √ 1129 2 + 1036.8 2 + 1145.2 2 = 1913.4 > 1732.8 is just evidence of the incorrectness of calculating the effective values of the harmonics of currents or voltages under conditions of dynamic change in the parameters of the harmonic composition per one measurement period T, when the higher harmonics did not have constancy in amplitudes and phase shifts.
In this way, the power components P i , Q i , and D ij must be methodically correctly calculated based on the instantaneous power values p(t) = u(t) · i(t), which contain the most reliable information on the operating modes of the AC network.
To answer the second question, it is necessary to analyze the conditions for the appearance of D = 0-when the curves i(t) and u(t) are similar, and D = 0, when the curve i(t) is distorted relative to the curve u(t).
The condition D = 0 with the explicit presence of the spectral composition of the harmonics of the distorted curve i(t) relatively u(t) unambiguously indicates the compensation of the components in expression (13), which is possible in the case when each component D ij = 0.
Let, for example, the curves u(t) and i(t) form two harmonics of voltages U m , U n and currents I m and I n . In this case, the harmonics m and n form active powers according to expression (10): P m = U m I m cos(φ m ) and P n = U n I n cos(φ n ).
They form reactive powers according to expression (11): Proceeding from expression (2), the following can be obtained: Hence, the distortion power can be expressed as: D = U m I n (U m I n − 2U n I m cos φ m · cos φ n ) + U n I m (U n I m − 2U m I n sin φ m · sin φ n ). (14) Therefore: This is possible under the condition of simultaneous achievement of the identity in the ratios of rms values for all harmonics of voltages and currents: and the equality of phase shifts in all harmonics of voltages and currents: If the harmonics U m , U n , I m , I n that correspond to conditions (15) and (16) can be represented in the form of rotating vectors on the complex plane, then the projections of the resulting voltage U = U m + U n and current vectors I = I m + I n , unfolded in time, form (Figure 4c) a distorted curve u(t) with respect to the curve i(t). At the same time, the hodographs u(t) and i(t) of the resulting vectors are similar and mutually shifted by π/2 ( Figure 4a). Along with this, the curves u(t) and i(t) characterize the formation of purely reactive power S = Q = Q m + Q n in the AC network.   An attempt to provide the similarity of the curves u(t) and i(t) with formal preservation of the identity in their harmonic composition (Figure 4b,d) leads to the fact that, in addition to reactive power Q, distortion power D appears in the AC network, which is easily verified by numerical calculations using expressions (1)-(3), (10), (11), and (13).
In the case of a large number of harmonics in the AC network, expression (14) will take on a more cumbersome form [28,33]: However, this does not fundamentally affect the behavior of the powers P, Q, and D in the AC network.
Thus, the power D, despite its unsuccessful name, does not characterize the degree of distortion of the curve i(t) shape relative to u(t), but it indicates the presence of such interactions of harmonics of currents and voltages with different ordinal numbers, which generate additional ballast powers that are not capable of performing useful work.

1.
Based on the analysis of wave and vector diagrams, the work has shown that the inactive power of N = √ S 2 − P 2 in non-sinusoidal AC networks is not an indivisible component, but, unambiguously, it can be further divided into orthogonal components: the distortion power of the D ∑ and the reactive power of the Q ∑ , that is, which can be unambiguously calculated analytically. Taking into account the behavior of the higher harmonics of currents and voltages, the power of Q and D can be generated in an arbitrary way by active elements with nonlinear volt-ampere characteristics.

2.
It is shown that power Q is formed during the interaction of voltage and current harmonics with the same ordinal numbers Q ∑ = ∑ Q i . In addition, power D is formed during the interaction of voltage and current harmonics with different ordinal numbers, due to the orthogonal properties of which D ∑ = ∞ ∑ i=1 D 2 i . At the level of instantaneous values of p(t) = u(t) · i(t), the interaction of these powers makes it possible to correct the shape of the supply currents and voltages in the AC network. However, the average value of the powers Q ∑ and D ∑ should be zero by definition over the working period.

3.
Analytical expressions are obtained for calculating all power components in a nonsinusoidal single-phase AC network, according to which, as an illustrative example, rotating generalized vectors of voltages and currents are constructed based on individual harmonics. The time base of the projections onto the axis of these generalized vectors yields wave diagrams of the instantaneous values of the u(t) and i(t) curves, which is of practical value for solving research problems. Thus, using generalized vectors, it has been found that compliance with the criteria for the identity of the ratios of the rms values for all harmonics of voltages and currents U m U n = I m I n and the equality of phase shifts of all harmonics of voltages and currents φ m = φ n to achieve D = 0 in a non-sinusoidal network does not provide a visual similarity of curves u(t) and i(t).

4.
It is shown that the visual similarity of curves u(t) and i(t) cannot serve as a criterion for the absence of distortion power D. Power D actually does not characterize the degree of distortion of the shape of the curve i(t) relative to u(t), but indicates the presence of ballast interactions of harmonics of currents and voltages with different ordinal numbers. These, in addition to reactive power, generate additional corresponding ballast powers that are also not capable of performing useful work.

5.
The performed analysis and calculations confirm the validity of the foundations of the theory of Fryze and Budeanu regarding the concepts of active and inactive powers, orthogonality of components P i , Q, D and achieving a power balance in a network with nonlinear volt-ampere characteristics of the load, despite the absence of a strict definition of power. However, the dynamic change in the amplitudes and phase shifts of individual current and voltage harmonics characteristic of modern loads over the time period that is shorter than the operating period T complicates the calculation of the actual values of the powers Q and D. In addition, if the harmonic analysis is performed incorrectly, it can lead to calculation errors. Funding: This research was supported by the NSTU development program. This research was supported by the TPU development program.