Cyclostationary Approach to the Analysis of the Power in Electric Circuits under Periodic Excitations

: This paper proposes a cyclostationary based approach to power analysis carried out for electric circuits under arbitrary periodic excitation. Instantaneous power is considered to be a particular case of the two-dimensional cross correlation function (CCF) of the voltage across, and current through, an element in the electric circuit. The cyclostationary notation is used for deriving the frequency domain counterpart of CCF—voltage–current cross spectrum correlation function (CSCF). Not only does the latter exhibit the complete representation of voltage–current interaction in the element, but it can be systematically exploited for evaluating all commonly used power measures, including instantaneous power, in the form of Fourier series expansion. Simulation examples, which are given for the parallel resonant circuit excited by the periodic currents expressed as a ﬁnite sum of sinusoids and periodic train of pulses with distorted edges, numerically illustrate the components of voltage–current CSCF and the characteristics derived from it. In addition, the generalization of Tellegen’s theorem, suggested in the paper, leads to the immediate formulation of the power conservation law for each CSCF component separately.


Introduction
Modern hardware systems performing digital signal processing can easily reach high complexity, since they typically consist of many specialized devices. In turns, the majority of these devices are assembled from a variety of electronic components, in which electric currents and voltages play the primary role in the description of phenomena at the physical level. For instance, one of the typical issues solved in the electromagnetic compatibility framework [1] is the detection of potential sources of interference, which may affect the stable work of the entire device or even neighboring devices. A traditional and easy to understand model describing the property of components can be built as an electric circuit being assembled from lumped elements only or lumped elements properly combined with parts of transmission lines.
The proper choice of an appropriate model is mainly based on two points. The first one is the primary role of the component, i.e., whether it is intended to be an antenna, microstrip line, connector, frequency selective filter, or the input port of an embedded integral circuit. The second point is determined by a certain relation between the spectrum width of the processed signals and the passband (or passbands) of the analyzed component. Thus, in the case of digital signal processing, where the spectrum is considered to be concentrated in a baseband with respect to the lowest resonant frequency range, the component's electric property can often be modelled as a simplified RC circuits [2]. However, if the resonance behavior cannot be neglected, a more accurate model, exhibiting weak resonant properties, may well be introduced by means of a generalized resonant circuit, either of parallel or serial type, with relatively low Q factor, typically not exceeding the value of 2.
The rest of the paper is organized as follows. The second section presents the theoretical model of voltage-current two-dimensional cross correlation function (CCF), which appears to be a specific extension of instantaneous power. Next, the cross spectral correlation density is derived in an explicit closed form, which allows redefining the common quantities used in the power analysis directly from its components. The particular case of linear time invariant elements is also briefly explained in Section 2. In Section 3, the evaluation of the suggested characteristics is illustrated with numerical examples, where a resonant parallel circuit is excited by different periodic waveforms. The discussion in Section 4 considers the extension of the power conservation law that the components of the spectral correlation functions obey. The paper ends with conclusions.

Voltage, Current and Instantenious Power Definitions
The instantaneous power delivered to any element in a circuit is considered as a timevarying signal given by the product of the instantaneous voltage v(t) across the element and the instantaneous current i(t) through it: In this research, we focus on strictly periodic waveforms representing the voltage and currents v(t) = v(t + T), i(t) = i(t + T), which can be expressed as their Fourier series, written both in real or complex forms and i(t) = I 0 + +∞ ∑ k=1 I k cos 2π where T is the period, which is the same for the voltage and current, F = 1/T is the frequency, V k and ϕ Vk are the voltage amplitude and phase, respectively, I k and ϕ Ik are the current amplitude and phase, respectively. M (kF) V and M (kF) I denote the complex valued Fourier coefficients for the complex exponential of frequency, kF, which is written in parentheses in the superscript, in the series of the voltage and current, respectively. These coefficients can be evaluated directly as soon as the real form parameters are known where the term x takes one of the symbolic values, v or i, denoting either the voltage or the current. It is important to highlight that the right-hand sides of Equations (3) and (4) are written as sums whose summands are related to frequencies; f = kF. These frequencies can take both positive and negative values, since they are related to the integral Fourier transform: where the complex exponentials exp(j2πf ) are involved. In other words, each harmonic of nonzero physical frequency kF is represented by two complex exponentials, exp(j2πf ), of frequencies ±kF, whose complex valued amplitude Appl. Sci. 2021, 11, 9711 5 of 25 are defined using (5) as soon as physical amplitude and phase are defined. These complex exponentials provide a reasonable mathematical model, which is often more convenient than an amplitude-phase form for performing the vast majority of formal manipulation.
Strictly speaking, the Fourier transform (6) does not exist as an ordinary function of frequency f if the signal x(t) is periodic. However, the complex form of the Fourier series, written for the voltage (3) and current (4), can be Fourier transformed using (6), if one treats the formal expression of their spectra as generalized functions, or distributions [44] in frequency domain, assembled of Dirac delta functions only:

Second Order Characteristics
The normal, or nonconjugate, dyadic cross correlation function (DCCF) of a voltage and current can be defined as a function of two time variables, t 1 and t 2 : where E denotes expectation and the asterisk * in the superscript marks the complex conjugation.
As long as the signals representing the voltage and current have been bounded to be deterministic, this literally means that: and their DCCF can be written with the expectation operator omitted: Besides, the second equal sign is possible, owing to the fact that both the voltage and current are real valued functions.
Since both the voltage and current are periodic functions, their DCCF is also a periodic function, sharing the same period, T, in both its argument: Having taken the two-dimensional Fourier transform over the pair of time variables, (t 1 , t 2 ), one obtains the function of the pair of frequency variables (f 1 , f 2 ): The substitution of (11) into (13) and further factorization of its exponential kernel will lead to the spectral representation (14) which is also referred to as Loève bifrequency cross spectrum. The frequency variables (f 1 , f 2 ) are the pairwise counterparts of the time variables (t 1 , t 2 ). Some formal difficulties may be pointed out: while the periodic functions v(t) and i(t), jointly producing R vi (t 1 , t 2 ), are substituted into (13), since they do not belong to the L 1 space, integration is carried out over the whole R 2 . However, we consider that this operation can be carried out in a generalized sense, and described with rigor enough for engineering applications. As long as (14) remains valid, the formal substitution of the Fourier series (3) and (4) into it yields: where a two-dimensional Dirac delta function can be treated as a product of two onedimensional Dirac delta functions: wherever it may be considered suitable for subsequent formal transformations [45].

Cyclostationary Description
The further analysis of cyclostationary characteristics based on second order twodimensional products can be performed more effectively if the pair of independent variables forming axes spanning the time domain (t 1 , t 2 ) or frequency domain (f 1 , f 2 ) are changed for other pairs of variables. This comes after the concept of symmetric transformation, exploited by Gardner [38] and illustrated in [46] by example of a second order CS process.
The substitution is aimed at transforming the DCCF into the symmetric form of the cross correlation function (CCF) which is a function of two variables: the current time, t, and the lag, or time shift, τ. Similarly, after applying the change in the frequency variables the bifrequency spectrum (14) will be transformed into the cross spectral correlation density where α is traditionally called a cyclic frequency, whereas f is considered as a frequency exploited in the integral Fourier transform (6). After substituting (19) into (14) and formally rearranging the summands, the analytical expression of voltage-current CSCD can be expressed in the form of a double summation, involving weighted and shifted two-dimensional Dirac delta functions δ(α, f ): where ((m)) 2 denotes the modulo operation, or the remainder after m is divided by 2, and complex scalars, S (m,k) vi , being indexed by the pair of integers (m, k), are defined via the coefficients put in the complex form of the Fourier series for the waveforms of the voltage (3) and current (4): where conj() marks the complex conjugation. Practically, the introduction of the modulo operation in (21) and (22) has been performed to formally reflect the fact that the Dirac delta functions at α = mF turn out to be shifted along the f -axis by the step F/2 for any odd integer, m.
Since CSCD is a generalized function whose components are present at cyclic frequencies in, at most, a countable set, A 2 = {α|α = mF}, the application of the decomposition (16) to the term δ(α, f ) will allow rewriting Expression (21) in the form: where S mF vi ( f ) represents spectral correlation function (SCF) at the cyclic frequency α = mF: Cyclic cross-correlation function (CCCF), which is a function of lag parameter τ, taken at cyclic frequency mF is introduced as where the integration is carried out over any interval of a one period length, e.g., [0, T) or [−T/2, T/2). Conversely, CCCF at a particular cyclic frequency can be treated as a coefficient in the Fourier series expansion of CCF: At the same time, CCCF R mF vi (τ) at cyclic frequency α = mF is the inverse Fourier transform of CSCF S mF vi ( f ) at the same cyclic frequency: The substitution of the explicit representation (24) into (27) leads to the representation of CCCF in the form of the Fourier series: CCCF (28) is a periodic function of τ. However, the value of its period depends on whether m is an even or odd integer. Thus, if m is even, the period is T = 1/F, but if m is odd, the period gets the double value, which is 2T. The latter follows from the fact that the greatest common divisor of the sequence F k = [k + 1/2]F, k ∈ Z is F/2, whose reciprocal determines the period of the Fourier series (28).

Instantenious Power Representation
The direct comparison of the formal expressions written for the instantaneous power (1) and CCF of the voltage and current (18) gives an idea that the former can be considered as a particular case of the latter, supposing the lag parameter τ has been chosen equal to zero: Appl. Sci. 2021, 11, 9711 8 of 25 The Fourier series (26), possessing functional coefficients, in turn, reduces to the Fourier series with scalar coefficients: where P (mF) denotes the complex coefficients in the Fourier series, which can be evaluated as long as τ in (28) is set to zero: The Fourier series (30) highlights the fact that, in a general case, the period of instantaneous power T = 1/F is the same as the period of the voltage and current described by arbitrary periodic waveforms.

Average Power
The coefficient at m = 0 in (30) is of particular interest due to the fact that P (0) represents the constant, or DC component, of the power p(t), also referred to as average power P av , which can be evaluated by means of time averaging over one period: Active power is traditionally interpreted as the rate of energy exchange in the long run. Provided the passive sign convention is assumed, an element of the circuit will be called a load if P av > 0. In contrast, if P av < 0, the element is an energy source. The remaining case, P av = 0, means the element is neither an absorber nor a supplier in the long run. However, it is important to notice that P av = 0 does not imply p(t) = 0.
The average power expression via CSCD components is followed immediately from (31): Equation (22), remaining valid for any m, can be specialized to the particular case where m = 0: S As long as both the voltage and current remain real valued functions, the coefficients in their expansions, (3) and (4), exhibit Hermitian symmetry: This property also establishes the Hermitian symmetry of the weights of Dirac delta functions at zero cyclic frequency: Using the well known identity held for any complex A: A + A * = 2 eA = 2 eA * , the average power can be rewritten as follows: where e marks the real part of the complex number, as well as m is marking the imaginary part below.

Average Power Distribution
Each component, which can be denoted as P vi , is treated as a part of the active power associated with frequency mF, as it appears to be contributing to the total average power separately. The consideration lying behind this idea is the following. Suppose one applied a band pass filter tuned so as to let only the nth harmonic in the Fourier series, representing the voltage (3) or current (4), pass through it. Then, according to (34), there will be only the pair made of S vi . However, as the voltage and current are real valued and, therefore, symmetry (36) holds, those components will inevitably satisfy the equation: The distribution of the components contributing to the average power, or furthermore, the extraction of their real (active) or imaginary (reactive) parts, among the discrete set of frequencies A 1 = {α|α = kF}, appears to be useful in the power analysis of the circuit, since it helps identify the main frequency channels of the energy absorption taking place in the circuit and measure their intensities.

Relation to Apparent Power
Apparent power is a quantity measured as a product of voltage and current effective values: The effective value of the voltage is typically considered as its root mean square (RMS) value, which can be expressed via the amplitudes or, alternatively, the coefficients of a complex Fourier series (3): Similarly, the effective value of the current can be evaluated using the parameters of its Fourier series (4), as: As soon as the substitution of (40) and (41) into (39) is made and rearranging the pairwise summand in accordance with equation (22) is performed, the apparent power can be expressed as: Appl. Sci. 2021, 11, 9711 10 of 25 Then, the dimensionless power factor can be evaluated as the ratio:

Relation to 1D Cross Correlation Function
The one-dimensional (conventional) cross correlation function (CCF1) can be obtained as CCCF taken at zero cyclic frequency: The frequency domain counterpart of CCF1 is known to be cross spectral density (CSD), which is, in turn, the CSCF at zero cyclic frequency: This consists of Dirac delta functions only since the voltage and current are periodic waveforms.
In contrast to the one-dimensional autocorrelation function (ACF1), which could have been constructed for the voltage waveform R v (τ) itself, the CCF1 is not obliged to exhibit the even symmetry as a function of lag shift τ. Similarly, CSD should not be thought of as a real valued function, or possessing real valued coefficients at Dirac delta functions, in the case of a generalized function, whereas power spectral densities (PSD) that could have been derived for the voltage waveform S v (f ) would be real valued. However, even-odd decomposition is applicable to the real valued CCF1 and, according to the properties of the Fourier transform, it relates to the real-imaginary decomposition made for CSD: where the vertical arrows designate the pairs "signal-spectrum," interconnected by the Fourier transform. The decomposition of the CSCF can be moved down to the weights of Dirac delta functions in its representation:

Linear Passive Elements and Causality
If an element of an electric circuit is linear, the interrelation between the voltage, v(t), across it and the current, i(t), through it can be established using the model of the linear time invariant (LTI) system. Such a system can be completely described in a time domain by means of its impulse response (IR). Impedance impulse response, z(t), describes the LTI system whose input is the current and output is the voltage: where * denotes the convolution expanded after the second equals sign. If the voltage is considered to be the system input and, consequently, the current is its output, the LTI system models admittance with impulse response y(t): If the LTI system under consideration belongs to the class of stable systems [47], it can be equivalently described by frequency response (FR), which is the Fourier transform of the system impulse response. Thus, the link between the voltage (3) and current (4) spectra can be established in its multiplicative form as: where Z(f ) and Y(f ) are the FR of the impedance and admittance, respectively. The relation (50), together with the properties of multiplying a generalized function by an ordinary one [44], gives the immediate relation between the coefficient of the complex form Fourier series for the voltage (3) and current (4): As soon as FRs are known, the link between CSCD (21) and the equivalent characteristic of the voltage and current can be established in explicit form: where S v (α, f ) and S i (α, f ) are the spectral correlation densities (SCD) of the voltage and current, respectively. Relation (52) allows writing the immediate representation of the CSCF at cyclic frequency, α, via the SCF of a voltage and current, denoted, respectively, as S In turn, the components of CSCD, i.e., the weights of Dirac delta functions, can be expressed in the closed form: where S (m,k) v and S (m,k) i can be evaluated in accordance with (22), assuming that both of their multipliers are chosen from the same signal representation.
In particular, at zero cyclic frequency α = 0, (54) gets reduced to: Furthermore, the average power components distributed over frequency are: where R( f ) = eZ( f ) and G( f ) = eY( f ) represent, respectively, the element resistance and conductance as the frequency dependent function.

Simulation Results
The typical load of electrical or electronic systems can be modelled by the use of circuits that consists of element sets, such as resistors, capacitors and inductors, connecting to each other. The resonance phenomenon can occur in such circuits under some conditions. The circuit becomes resonant and has selecting properties. The study of the power distribution of such circuits is of particular interest.
The parallel resonant circuit with feeding current source, which is presented in Figure 1a, was chosen for numerical simulation. The following parameters of the circuit were established: R 1 = 20 Ω, R 2 = 125 Ω, L = 0.8 mH, C = 320 nF. Those parameters make it possible to determine the following circuit characteristics: Q factor Q = 1.25, resonant frequency f r = 10 MHz, characteristic impedance ρ = 50 Ω.
cuits that consists of element sets, such as resistors, capacitors and inductors, connecting to each other. The resonance phenomenon can occur in such circuits under some conditions. The circuit becomes resonant and has selecting properties. The study of the power distribution of such circuits is of particular interest.
The parallel resonant circuit with feeding current source, which is presented in Figure 1a, was chosen for numerical simulation. The following parameters of the circuit were established: R1 = 20 Ω, R2 = 125 Ω, L = 0.8 mH, C = 320 nF. Those parameters make it possible to determine the following circuit characteristics: Q factor Q = 1.25, resonant frequency fr = 10 MHz, characteristic impedance ρ = 50 Ω. The FR describing the input impedance of the circuit is defined as: The input impedance of the circuit versus frequency is presented in Figure 1b, where its absolute value and phase are combined in the same double sided plot.

Finite Sum of Harmonics Current Excitation
The periodic current i(t), consisting of three sinusoids and a dc component with period T = 0.1 μs, is considered as the first example: where the frequency is F = 1/T = 10 MHz. The current, i(t), voltage, v(t), and power, p(t), signals are depicted in Figure 2. They are periodic functions with the same period, T. The FR describing the input impedance of the circuit is defined as: The input impedance of the circuit versus frequency is presented in Figure 1b, where its absolute value and phase are combined in the same double sided plot.

Finite Sum of Harmonics Current Excitation
The periodic current i(t), consisting of three sinusoids and a dc component with period T = 0.1 µs, is considered as the first example: where the frequency is F = 1/T = 10 MHz. The current, i(t), voltage, v(t), and power, p(t), signals are depicted in Figure 2. They are periodic functions with the same period, T. The spectra of current, i(t), and voltage, v(t), are presented in Figure 3a,b ,respectively, as the weights of the Dirac delta functions representing Fourier coefficients (5) and measured in [A] and [V], respectively. The colors of the circles designate the phase of each Fourier coefficient, according to the color bar, ranging from −π to π. The same visualizing technique will be used for revealing the phase information in other figures. The spectra of current, i(t), and voltage, v(t), are presented in Figure 3a,b, respectively, as the weights of the Dirac delta functions representing Fourier coefficients (5) and mea-sured in [A] and [V], respectively. The colors of the circles designate the phase of each Fourier coefficient, according to the color bar, ranging from −π to π. The same visualizing technique will be used for revealing the phase information in other figures. The spectra of current, i(t), and voltage, v(t), are presented in Figure 3a,b ,respectively, as the weights of the Dirac delta functions representing Fourier coefficients (5) and measured in [A] and [V], respectively. The colors of the circles designate the phase of each Fourier coefficient, according to the color bar, ranging from −π to π. The same visualizing technique will be used for revealing the phase information in other figures. It can be seen in Figure 3b that, for the voltage spectra, the maximum in absolute value weights of the delta function are at frequencies ±10 MHz, which correspond to the circuit resonant frequency, fr, and the fundamental frequency, F, of considered signals. Figure 4 shows the weights of two-dimensional Dirac delta functions, forming the bifrequency spectrum, Svi(f1, f2), of the current, i(t), and the voltage, v(t), in accordance with (15). These delta functions are arranged with step 10 MHz, which corresponds to the fundamental frequency, F, of signals i(t), v(t) and p(t). It can be seen that the section of Svi(f1, f2) at f1 = ±10 MHz have maximum weights. It can be seen in Figure 3b that, for the voltage spectra, the maximum in absolute value weights of the delta function are at frequencies ±10 MHz, which correspond to the circuit resonant frequency, f r , and the fundamental frequency, F, of considered signals. Figure 4 shows the weights of two-dimensional Dirac delta functions, forming the bifrequency spectrum, S vi (f 1 , f 2 ), of the current, i(t), and the voltage, v(t), in accordance with (15). These delta functions are arranged with step 10 MHz, which corresponds to the fundamental frequency, F, of signals i(t), v(t) and p(t). It can be seen that the section of S vi (f 1 , f 2 ) at f 1 = ±10 MHz have maximum weights.  The absolute values, in μW, and phases, in radians, of the delta function weights forming the bifrequency spectrum, Svi(f1, f2), are presented in Table 1, where each cell corresponds to a pair of frequencies, f1 and f2, defined in the rows and columns, respectively.  The absolute values, in µW, and phases, in radians, of the delta function weights forming the bifrequency spectrum, S vi (f 1 , f 2 ), are presented in Table 1, where each cell corresponds to a pair of frequencies, f 1 and f 2 , defined in the rows and columns, respectively. The alternative representation of the bifrequency spectrum is the diagram shown in Figure 5a, where each delta-function is depicted as a circle whose area is proportional to the absolute value of S vi (f 1 , f 2 ). This way of presenting a 3D image of S vi (f 1 ,f 2 ), on plane, can be convenient for the visualization and rapid estimation of the relative values of S vi (f 1 ,f 2 ) components. We preferred to use a diagrammatic representation in subsequent typical cases. The frequencies f 1 , f 2 range from −30 to 30 MHz.   Figure 6 shows the distribution of the real and imaginary parts of average power. It can be seen that the real part contributes, in average power, more than the imaginary part. The weights of the delta function of the real part are highest at frequencies ±10 MHz, which is the circuit resonant frequency, fr, and the fundamental frequency, F, of considered signals.    Figure 6 shows the distribution of the real and imaginary parts of average power. It can be seen that the real part contributes, in average power, more than the imaginary part. The weights of the delta function of the real part are highest at frequencies ±10 MHz, which is the circuit resonant frequency, f r , and the fundamental frequency, F, of considered signals.
functions are presented. The cyclic frequency, α, ranges from −60 to 60 MHz, the frequency, f, ranges from −30 to 30 MHz. Figure 6 shows the distribution of the real and imaginary parts of average power. It can be seen that the real part contributes, in average power, more than the imaginary part. The weights of the delta function of the real part are highest at frequencies ±10 MHz, which is the circuit resonant frequency, fr, and the fundamental frequency, F, of considered signals.    These power quantities can be verified using the traditional method. The periodic voltage v(t), which is depicted in Figure 1a, can be found analytically using frequency characteristic (57). Its formal representation, which is written with 2-digit precision for the sake of compactness, is as follows: The effective values of the voltage Veff and current Ieff, evaluated using (40) and (41), are, respectively, 153.0825 mV and 3.5355 mA, whose product yields the apparent power by definition (39): Papp = 541.2 μVA.
As soon as the voltage and current are known in closed form, the average power Pav can also be evaluated analytically using (32) where the substitution p(t) = v(t)i(t) is made and the integration limits are chosen from 0 to T = 1/F. This yields Pav = 386.5 μW, which is exactly the same value as found above.

Pulse Train Current Excitation
The second example of signal i(t) to consider is a periodic train of rectangular pulses followed with period T = 0.1 μs, with exponentially smoothed leading and trailing edges. The current i(t) can be defined as: The following quantitative power characteristics were calculated: the average power (33) P av = 386.5 µW, apparent power (39) P app = 541.2 µVA, power factor (43) PF = 0.71.
These power quantities can be verified using the traditional method. The periodic voltage v(t), which is depicted in Figure 1a, can be found analytically using frequency characteristic (57). Its formal representation, which is written with 2-digit precision for the sake of compactness, is as follows: The effective values of the voltage V eff and current I eff , evaluated using (40) and (41), are, respectively, 153.0825 mV and 3.5355 mA, whose product yields the apparent power by definition (39): P app = 541.2 µVA.
As soon as the voltage and current are known in closed form, the average power P av can also be evaluated analytically using (32) where the substitution p(t) = v(t)i(t) is made and the integration limits are chosen from 0 to T = 1/F. This yields P av = 386.5 µW, which is exactly the same value as found above.

Pulse Train Current Excitation
The second example of signal i(t) to consider is a periodic train of rectangular pulses followed with period T = 0.1 µs, with exponentially smoothed leading and trailing edges. The current i(t) can be defined as: where where u(t) is Heaviside function, and the time constant ∆ = 1/α is equal to one tenth of the width of each pulse. The current, i(t), voltage, v(t), and power, p(t), signals are depicted in Figure 8. They are periodic functions with the same period, T. The spectra of current i(t) and voltage v(t) are presented in Figures 9 and 10, respectively. The spectrum (Figure 9) of the current (60) contains more harmonics" in comparison with the spectrum (Figure 3a) of the current (58), discussed above. However, the voltage spectra in Figures 3b and 10 consist of a similar number of meaningful delta functions because the absolute values of delta functions related the components higher than the third become negligible. For the case of the voltage excited by the pulse current, this can be explained by the selective properties of resonant circuits. It can be seen in Figure 10 that delta functions having the maximum weights by absolute value are at frequencies ±10 MHz, which correspond to the circuit resonant frequency, fr, and the fundamental frequency, F, of the considered signals.   Figure 11 shows the diagram of CSCD Svi(α,f) found in accordance with (21). It has more components than the diagram in Figure 5b because of the richer spectral composition of the rectangular pulse train current (60). The spectra of current i(t) and voltage v(t) are presented in Figures 9 and 10, respectively. The spectrum (Figure 9) of the current (60) contains more harmonics" in comparison with the spectrum (Figure 3a) of the current (58), discussed above. However, the voltage spectra in Figures 3b and 10 consist of a similar number of meaningful delta functions because the absolute values of delta functions related the components higher than the third become negligible. For the case of the voltage excited by the pulse current, this can be explained by the selective properties of resonant circuits. It can be seen in Figure 10 that delta functions having the maximum weights by absolute value are at frequencies ±10 MHz, which correspond to the circuit resonant frequency, f r , and the fundamental frequency, F, of the considered signals. The spectra of current i(t) and voltage v(t) are presented in Figures 9 and 10, respectively. The spectrum (Figure 9) of the current (60) contains more harmonics" in comparison with the spectrum (Figure 3a) of the current (58), discussed above. However, the voltage spectra in Figures 3b and 10 consist of a similar number of meaningful delta functions because the absolute values of delta functions related the components higher than the third become negligible. For the case of the voltage excited by the pulse current, this can be explained by the selective properties of resonant circuits. It can be seen in Figure 10 that delta functions having the maximum weights by absolute value are at frequencies ±10 MHz, which correspond to the circuit resonant frequency, fr, and the fundamental frequency, F, of the considered signals.   Figure 11 shows the diagram of CSCD Svi(α,f) found in accordance with (21). It has more components than the diagram in Figure 5b because of the richer spectral composition of the rectangular pulse train current (60). The spectra of current i(t) and voltage v(t) are presented in Figures 9 and 10, respectively. The spectrum (Figure 9) of the current (60) contains more harmonics" in comparison with the spectrum (Figure 3a) of the current (58), discussed above. However, the voltage spectra in Figures 3b and 10 consist of a similar number of meaningful delta functions because the absolute values of delta functions related the components higher than the third become negligible. For the case of the voltage excited by the pulse current, this can be explained by the selective properties of resonant circuits. It can be seen in Figure 10 that delta functions having the maximum weights by absolute value are at frequencies ±10 MHz, which correspond to the circuit resonant frequency, fr, and the fundamental frequency, F, of the considered signals.   Figure 11 shows the diagram of CSCD Svi(α,f) found in accordance with (21). It has more components than the diagram in Figure 5b because of the richer spectral composition of the rectangular pulse train current (60).  Figure 11 shows the diagram of CSCD S vi (α,f ) found in accordance with (21). It has more components than the diagram in Figure 5b because of the richer spectral composition of the rectangular pulse train current (60).           Figure 13 shows the spectrum of instantaneous power p(t). The delta functions forming the spectrum of p(t) are at the same frequencies as the delta functions forming the spectrum of instantaneous current, i(t), and voltage, v(t), signals.  Figure 13 shows the spectrum of instantaneous power p(t). The delta functions forming the spectrum of p(t) are at the same frequencies as the delta functions forming the spectrum of instantaneous current, i(t), and voltage, v(t), signals. The following quantitative power characteristics were calculated: the average power (33) P av = 199.4 µW, apparent power (39) P app = 264.7 µVA, power factor (43) PF = 0.75.

Trasmission Line Example
The third example is a lossless transmission line (Figure 14a) with a load, which coincides with the load of the circuit in Figure 1a.

Trasmission Line Example
The third example is a lossless transmission line (Figure 14a) with a load, which coincides with the load of the circuit in Figure 1a. The input impedance of transmission line can be evaluated as [48]: where c is the speed of wave propagation, l is the length of transmission line, Z is defined by the formula (1), and Z0 is characteristic impedance. For numerical simulation, the following parameters were chosen: c = 3 × 10 8 m/s, Z0 = 50 Ω, l = 12 m. Consider the current source, which is the same as in the second example, defined by formulas (60) and (61). The input impedance versus frequency is presented in Figure 14b, where its absolute value and phase are combined in the same double sided plot.
The current, i(t), voltage, v(t), and power, p(t), signals are depicted in Figure 15. They are periodic functions with period T. The input impedance of transmission line can be evaluated as [48]: where c is the speed of wave propagation, l is the length of transmission line, Z is defined by the formula (1), and Z 0 is characteristic impedance. For numerical simulation, the following parameters were chosen: c = 3 × 10 8 m/s, Z 0 = 50 Ω, l = 12 m. Consider the current source, which is the same as in the second example, defined by formulas (60) and (61). The input impedance versus frequency is presented in Figure 14b, where its absolute value and phase are combined in the same double sided plot.
The current, i(t), voltage, v(t), and power, p(t), signals are depicted in Figure 15. They are periodic functions with period T. by the formula (1), and Z0 is characteristic impedance. For numerical simulation, the following parameters were chosen: c = 3 × 10 8 m/s, Z0 = 50 Ω, l = 12 m.
Consider the current source, which is the same as in the second example, defined by formulas (60) and (61). The input impedance versus frequency is presented in Figure 14b, where its absolute value and phase are combined in the same double sided plot.
The current, i(t), voltage, v(t), and power, p(t), signals are depicted in Figure 15. They are periodic functions with period T.  The spectrum of current i(t) is the same as in Figure 9, the spectrum of voltage v(t) is presented in Figure 16. It can be seen that the maximum by absolute weight value delta function are at frequencies ±20 MHz. This can be explained by the higher absolute value of the input impedance (62) at frequencies ±20 MHz, compared to ±10 MHz. The spectrum of current i(t) is the same as in Figure 9, the spectrum of voltage v(t) is presented in Figure 16. It can be seen that the maximum by absolute weight value delta function are at frequencies ±20 MHz. This can be explained by the higher absolute value of the input impedance (62) at frequencies ±20 MHz, compared to ±10 MHz.  Figure 17 shows the diagram of CSCD Svi(α,f) (21). It has even more components than diagram in Figure 11 because more spectral components of |V(f)| in Figure 16 have sufficient value, compared to |V(f)| in Figure 10.    Figure 17 shows the diagram of CSCD S vi (α,f ) (21). It has even more components than diagram in Figure 11 because more spectral components of |V(f )| in Figure 16 have sufficient value, compared to |V(f )| in Figure 10. The spectrum of current i(t) is the same as in Figure 9, the spectrum of voltage v(t) is presented in Figure 16. It can be seen that the maximum by absolute weight value delta function are at frequencies ±20 MHz. This can be explained by the higher absolute value of the input impedance (62) at frequencies ±20 MHz, compared to ±10 MHz.  Figure 17 shows the diagram of CSCD Svi(α,f) (21). It has even more components than diagram in Figure 11 because more spectral components of |V(f)| in Figure 16 have sufficient value, compared to |V(f)| in Figure 10.  Figure 18 shows the distribution of the real and imaginary parts of average power. The weights of the delta function of the real part are maximum at frequencies ±10 MHz, which correspond to fundamental frequency, F, of the considered signals. The imaginary Figure 17. CSCD for the case of transmission line and pulse train current. Figure 18 shows the distribution of the real and imaginary parts of average power. The weights of the delta function of the real part are maximum at frequencies ±10 MHz, which correspond to fundamental frequency, F, of the considered signals. The imaginary part of delta function weights is maximum by its absolute value at frequencies ±30 MHz.  Figure 19 shows the spectrum of p(t). The delta functions forming the spectrum of instantaneous power p(t) are at the same frequencies as the delta functions forming the spectrum of instantaneous current i(t) and voltage v(t). The first five spectral components are higher with respect to the dc component than took place in the case of the resonant circuit without transmission line ( Figure 13).

Discussion on Conservation Law
From its definition (18), CCF suggests performing the multiplication of time shifted versions of the voltage and current. We are going to show, in this section, that this can lead us to important properties linking the components of CSCD (22) over all the elements of which any circuit consists.
Let us consider an electric circuit assembled of N elements, whose currents and voltages are denoted by in(t) and vn(t), respectively. Tellegen's theorem [10] states that, whenever the currents obey the equation set constituting Kirchhoff's current law and the voltages obey the equation set constituting Kirchhoff's voltage law (KVL), the instantaneous powers pn(t) = vn(t)in(t) will submit to the balance equation:  Figure 19 shows the spectrum of p(t). The delta functions forming the spectrum of instantaneous power p(t) are at the same frequencies as the delta functions forming the spectrum of instantaneous current i(t) and voltage v(t). The first five spectral components are higher with respect to the dc component than took place in the case of the resonant circuit without transmission line ( Figure 13).  Figure 19 shows the spectrum of p(t). The delta functions forming the spectrum of instantaneous power p(t) are at the same frequencies as the delta functions forming the spectrum of instantaneous current i(t) and voltage v(t). The first five spectral components are higher with respect to the dc component than took place in the case of the resonant circuit without transmission line ( Figure 13).

Discussion on Conservation Law
From its definition (18), CCF suggests performing the multiplication of time shifted versions of the voltage and current. We are going to show, in this section, that this can lead us to important properties linking the components of CSCD (22) over all the elements of which any circuit consists.
Let us consider an electric circuit assembled of N elements, whose currents and voltages are denoted by in(t) and vn(t), respectively. Tellegen's theorem [10] states that, whenever the currents obey the equation set constituting Kirchhoff's current law and the voltages obey the equation set constituting Kirchhoff's voltage law (KVL), the instantaneous powers pn(t) = vn(t)in(t) will submit to the balance equation: The following quantitative power characteristics were calculated: the average power (33) P av = 297.1 µW, apparent power (39) P app = 537.2 µVA, power factor (43) PF = 0.55.

Discussion on Conservation Law
From its definition (18), CCF suggests performing the multiplication of time shifted versions of the voltage and current. We are going to show, in this section, that this can lead us to important properties linking the components of CSCD (22) over all the elements of which any circuit consists.
Let us consider an electric circuit assembled of N elements, whose currents and voltages are denoted by i n (t) and v n (t), respectively. Tellegen's theorem [10] states that, whenever the currents obey the equation set constituting Kirchhoff's current law and the voltages obey the equation set constituting Kirchhoff's voltage law (KVL), the instantaneous powers p n (t) = v n (t)i n (t) will submit to the balance equation: The equation expressing KVL in the lth loop looks like the linear combination of the voltages: where λ ln is a coefficient taking one of three possible values: 1, −1 or 0. If the element has not been included in the lth loop, λ ln = 0. If the nth element is a part of the loop, the mutual directions of the element voltage and the loop trace determine λ ln = 1 if they coincide, and λ ln = −1 otherwise. The structure of equation set describing KVL establishes the fact that they are invariant to time shift θ applied to all voltages, v n (t + θ), simultaneously. Similarly, the equation for KCL in the pth node is the linear combination of the currents: where β pn is a coefficient taking values out of the set {1, −1, 0}. If the nth element is not adjacent to the pth loop, β pn = 0. If the current associated with the nth element is leaving the node, β pn = 1, and if the current is entering the node, β pn = −1. Provided the time shift σ applied to all currents, i n (t + σ), simultaneously, the KCL equation set remains valid. More importantly, time shifts applied to voltages and currents do not need to be equal to each other. Returning to the genuine Tellegen's theorem, we may conclude that: where the additional index, n, written in the subscript of R vi n (t, τ), highlights the possession of CCF by the nth element of the circuit. Equation (66) provides us with a generalized version of Tellegen's theorem in a time domain. However, the series made of Dirac delta functions in (21) forms the frequency domain counterpart of orthogonal basis for formal expansion of CCF into double Fourier series: In other words, provided the complex exponential term is rewritten into a concise expression: the orthogonality condition will hold: where δ pq denotes the Kronecker delta, which is equal to 1 if p = q and 0 otherwise. Having applied the component wise orthogonality property to relation (66), we come to a new version of the conservation law: which means that each component of (22) obeys the conservation law separately.
The direct yet useful corollary of conservation law (70) is the fact that an appropriately chosen function, depending on a subset containing some of the components, will also sub-mit to the conservation law. In particular, a sum of individually transformed components can assemble a function of such a type: where g() defines a one variable function, satisfying a linearity condition, and V ⊂ Z 2 determines the set of pairs (m, k) chosen for performing the summation. For instance, the average power (33) and the coefficients of Fourier series (31), which is expressing the instantaneous power, can be treated this way if the function g(x) = x and V = {M} × Z, where scalar M is the coefficient number: it will be equal to zero for the average power. Another example is combination: which can be fairly treated as the reactive power, according to Budeanu's theory [15]. It also obeys to the conservation law, as it follows from (70), with g(x) = m(x) and V = {0} × N.
On the other hand, some other transformations may not agree with componentwise Equation (70). Thus, it should not be surprising that the apparent power (42), which actually does not obey the conservation law, as it has not got any reasons, rooted in (70), to be governed by it.

Conclusions
The approach proposed in this paper to power analysis in electric circuits under arbitrary periodic excitation suggests a systematic framework, which is mainly based on two ideas. The first is expressing the voltage across, and current through, an element by their complex valued Fourier series; this allows one to consider them as first order CS processes. The second idea consists in applying the set of conventional CS tools, such as CCF and SCF, to the voltage-current pair in the element under analysis.
The main advantages of the proposed method are generally based on its methodology. The method reveals the elementary, so called "atomic", power components, which still obey the conservation law. Proper combinations of these components, including nonlinear, as for apparent power by (42), can be used for representing any quantity, characterizing power. As it is shown in this paper, the description based on CS properties perfectly matches the nature of instantaneous power as a product of the voltage and current related to the same element of electric circuit, or one port.
On the one hand, such an approach may look overly complicated since there exist alternative approaches to the same problem based on real valued, or amplitude and phase representation. The introduction of complex numbers for describing power components may be looking excessive, since complex Fourier series are more common for researchers dealing with signal processing tasks than for researchers working in the fields of electric machines and power delivery.
At the current development stage, the proposed method was implemented as a specialized software utility used for postprocessing the data obtained as the output of an applied software performing circuit analysis in steady state mode. One of the classes of the problems typically solved by means of this software is the numerical optimization of the power factor, which is to be achieved by the allowed variation of the load, e.g., by adding extra reactance or changing the length of the transmission line. However, the classical optimization approach is usually performed for a single harmonic, whereas the other harmonics are considered primarily as the distortion. In contrast, the proposed method utilizes the compact component wise representation of the full instantaneous power in the element of the circuit excited by an arbitrary periodic waveform. This paves the way to more accurate optimization procedures, taking into account higher harmonics of the signal.
Although, from an applied mathematician's point of view, CCF can be treated as an extended model of the time varying instantaneous power, in which the auxiliary lag parameter has been introduced, the physical meaning of CCF may still require deeper clarification. Thus, although each of the components given by (22) obeys conservation law (70) separately, they will contribute to the Fourier series representing the instantaneous power (30) only in aggregated form (31). Other aggregations are carried out to form the average and reactive powers. In addition, the question of the realizability of the circuit, which can implement the product of the time-translated voltage and currents, is omitted in this paper, leaving this issue open.
Future development of the CS framework proposed in this paper can be carried out in several directions. The first one consists in its application to the circuits containing nonlinear or time varying elements. As soon as the simplified model of a nonlinear element is involved, the relation between the voltage and current in such an element can be described by its volt-ampere characteristics. In contrast, the time varying behavior of an element can be described by linear periodically time-varying, which are perfectly covered by the CS theory. In both cases, instantaneous power still remains periodic and keeps its representative ability as the quantitative measure of the voltage-current interaction.
Another generalization of the proposed theory can be performed if the circuit is excited by voltages and currents whose sources are modeled as random processes. Moreover, promising results are expected for cases of wide sense stationary and wide sense cyclostationary random processes. However, since the voltage and current cease to be deterministic, the expectation starts taking on an important part in the power representation. Not only can this procedure be carried out in a probabilistic manner using ensemble averaging, but it can also be performed in the infinite time limiting process, according to the concepts developed within fraction of time theory [49]. Those may also become the promising topics for further research.
Finally, after proper generalization, the proposed approach can be used for describing the power in a multiport, where voltages and currents can be measured in dedicated ports rather than separate elements. Furthermore, starting from voltage-current interaction, the approach may be used in electromagnetics for describing power produced by appropriate cross product electric and magnetic fields expressed as component vectors in three dimensional space.