Research Article  Open Access
SteadyState Response of Periodically Switched Linear Circuits via Augmented TimeInvariant Nodal Analysis
Abstract
We focus on the simulation of periodically switched linear circuits. The basic notation and theoretical framework are presented, with emphasis on the differences between the linear timeinvariant and the timevarying cases. For this important class of circuits and sources defined by periodic signals, the computation of their steadystate response is carried out via the solution of an augmented timeinvariant MNA equation in the frequencydomain. The proposed method is based on the expansion of the unknown voltages and currents in terms of Fourier series and on the automatic generation of augmented equivalents of the circuit components. The above equivalents along with the information on circuit topology allow creating, via circuit inspection, a timeinvariant MNA equation, the solution of which provides the coefficients of both the time and the frequencydomain responses of the circuit. Analytical and numerical examples are used to stress the generality and benefits of the proposed approach.
1. Introduction
Nowadays, a number of alternative methods and tools for the analysis of timevarying circuits and systems that exhibit a periodic switching behavior have been developed and are currently available to designers. They are used to predict both the time and the frequencydomain responses of electrical and electronic circuits that are characterized by the possibly periodic activity of internal switches during the initial design phase. Wellknown examples are the power converters that are massively used to energize almost any equipment or appliance.
A seminal contribution to circuit analysis has been proposed by Zadeh and Desoer in [1, 2] where the basic definitions and mathematical tools are presented as the natural extension of the results available for the linear timeinvariant (LTI) case [3]. In the 70s, Liou developed a closedform solution for the complete frequencydomain response of periodically switched linear circuits [4]. Based on the above results, Strom and Signell [5] extended the previous works to more complex networks with an arbitrary number of switches. The common background of the above contributions is the analysis of a circuit by means of a stateequation formalism, intended to compute the socalled generalized transfer functions and impulse responses of timevarying circuits. These approaches, however, require a strong mathematical development and heavy numerical computations. Yuan and Opal [6, 7] introduced an innovative simulation approach via the classical modified nodal analysis (MNA) tool. In the latter contribution, the circuits are frozen and solved in all the different operating states of the switches. The dynamical effects of switches on the circuit solution are then taken into account via a suitable coupling involving the dynamical elements within the net. More recently, other techniques such as Floquet theory [8], truncated harmonic balance [9], and timeinvariant multifrequency (TIMFS) model [10] have been proposed in order to study and to simulate linear periodically timevarying systems both in frequency and in timedomains. Without loss of generality, the readers should also refer to [11, 12] where additional properties of timevarying systems are available and to [13–15] where applications are presented.
The present contribution extends the results of previously published papers in this field and provides an alternative approach to the steadystate analysis of periodically switched linear (PSL) circuits. The proposed method is based on the generation of augmented timeinvariant characteristics of the PSL elements directly in frequencydomain via the harmonic balance technique and on the computation of steadystate response of a PSL circuit by means of the standard solution of an augmented linear timeinvariant nodal equation generated from circuit inspection only. The proposed solution dramatically reduces the mathematical complexity and offers a modular approach to circuit analysis, leading to accurate results with remarkable simulation speedup. The paper collects a brief overview of the essential theory needed to handle PSL circuits in order to provide a selfcontained document allowing the readers to readily apply the proposed methodology to their own application problems. To this aim, simple analytical examples and topological intuition are used to demonstrate the generality of the method. Its feasibility and strength are instead provided by a real application example involving a switching power converter in [16]. The proposed approach is used to predict and to understand the conducted emission of switching devices, improving the stateoftheart results [17, 18]. In the above assessment, the conducted disturbances are quantified by means of the frequency spectrum of the steadystate response of the current supplying the converter.
The remaining part of the paper is organized as follows. Section 2 outlines the basic notation and definitions of LTI systems. Their generalization to the case of timevarying circuits is included in Section 3. Section 4 focuses on the important class of periodically switched circuits and illustrates the basic constitutive operators and rules governing their behavior. The core of the proposed method is presented in Sections 5 and 6 where augmented characteristic equations of resistive PSL elements are derived and used to generate a timeinvariant MNA equation describing a generic PSL circuit. Specific emphasis is given to a topological approach for the construction of the advocated augmented nodal equations in matrix form. Section 7 collects the results of the application of the proposed solution to two example circuits, with additional details on the accuracy and efficiency of the method. Summary and conclusions are given in Section 8.
2. Linear TimeInvariant Systems
This section summarizes the standard symbols and definitions used in this paper. A generic singleinput singleoutput LTI system [3] is defined by the following compact notation: where and are the input and the output signals, respectively, and is the operator or rule accounting for the system behavior.
For a system defined by (1), the evolution of the output signal can be computed from the impulse response by means of the convolution integral that writes where , being the Dirac delta function. Additionally, the wellknown counterpart of (2) in the frequencydomain writes where is the frequencydomain network function obtained via the Fourier transform of the impulse response ; that is,
3. Linear TimeVarying Circuits
This section collects the essential theoretical results needed to extend the characterization of LTI systems to the more complex case of timevarying circuits. The readers should refer to [1] for a more exhaustive and indepth treatment. Here, the discussion is aimed at presenting the notation and the main mathematical tools allowing handling these systems.
For the sake of illustration, the example of Figure 1 is considered. It shows a simple timevarying circuit consisting of two resistors, one capacitor, and one switch (Figure 1(a)) and its system representation in terms of interconnected blocks (Figure 1(b)). From the above example, clearly, the output of the system depends on the position of the switch . For a fixed state of , the output can be defined by either or and the standard convolution integral given by (2) can be used to compute the system response. However, to account for a possible timevarying activity of the system, (2) must be replaced by where the generalized impulse response is defined by (where in a causal linear timevarying system for ) according to [1, 11]. It is relevant to remark that, in the previous equation, plays the role of the observation time and represents the socalled excitation time accounting for the position of the delta function feeding the system. In other words, the excitation time is used to explore the timevarying nature of the system under different operating conditions.
(a)
(b)
The dependence of the generalized impulse response on the additional variable requires the definition of the bifrequency transfer function that can be computed by the impulse response via the following integral:
The above equation extends the Fourier transform of the classical impulse response to the twodimensional case (see [1] for additional details). Similarly, the generalization of (3) becomes
From the previous result, the bifrequency transfer function can be considered as a map between the input and the output , where the variables and are the input and output frequencies, respectively.
Equation (7) suggests that a sinusoidal excitation with angular frequency feeding the system produces an output response characterized by a possibly richer frequencydomain behavior (i.e., turns out to be a function of ). On the contrary, it is well known that LTI systems operating in the sinusoidal steadystate at produce responses characterized by the same input frequency. In the latter case, would consist of two delta functions located at .
In this framework, an important role is played by the timevarying transfer function: It relates the timedomain response of the system to the corresponding frequencydomain excitation , where the variable needs to be replaced by , leading to
Based on [11], that provides a summary of direct and inverse transformations involving the above defined operators; the bifrequency transfer can be alternatively obtained from (8) via
The above equations constitute the basic set of definitions needed to solve both analytically and numerically timevarying circuits. This overview also highlights the main differences between timeinvariant and timevarying systems. In the latter case, the relations among the input and output variables in both time and frequencydomain involve more complex operators rather than simple products and convolutions. The mathematics is unavoidably more complex. However, the framework is well established and allows us to effectively employ the readily available results to real application examples.
4. Periodically Switched Linear Circuits
This section deals with the analysis of PSL circuits by means of the generalized transfer functions described in the previous section. A PSL circuit can be seen as a particular case of a timevarying linear system, where the circuit configuration changes periodically in time due to the effect of the switches. The example of Figure 1 becomes a PSL circuit when the position of the switch is periodically changed between the close and the open states. It is important to remark that this class of circuits represents an important subset of timevarying systems with a large number of applications in different engineering domains.
For PSL circuits, the timevarying transfer function (8) turns out to be a periodic function of that can be expanded in Fourier series, leading to [6] where is the fundamental angular frequency of the switch ( being the characteristic period) and are new objects called aliasing transfer functions, defined as
Substituting (11) into (10) gives
Hence, the output in frequencydomain can finally be computed via (7) and writes The above equation is the generalization of (3) for linear periodically switched systems. This can be appreciated by observing that (14) is equivalent to (3) when .
In order to better understand the behavior of a PSL circuit, it is important to specialize the above equations for the case of a cisoidal excitation . The frequencydomain response of the system is obtained by substituting the frequencydomain description of the input signal into (14), leading to
The above equation, that is the Fourier transform of a signal described by its Fourier series, can be rewritten in timedomain by means of a complex sum of trigonometric functions with angular frequencies :
From the previous results it is clear that PSL circuits are completely described in both frequency and timedomains by the aliasing transfer functions (see [5, 11] for a more general definition and additional details). It is relevant to remark that the complexity of the analytical solution of a PSL circuit is reasonable and affordable for the case or resistive switched circuits only (as the example of Figure 1 with the capacitor replaced by a shortcircuit). In the more general case of dynamical circuits, however, the complexity increases and numerical methods or specialized techniques are needed to compute the solution of the system (e.g., see [4–6]).
In order to overcome the above limitations, the next two sections develop an alternative and effective approach to circuit analysis. A generic circuit is seen as the interconnection of linear timeinvariant elements (as resistors, capacitors, inductors, …) and PSL resistive elements (as the twoterminal element of Figure 2). Thereby, the circuit is solved by suitably reinterpreting the characteristics of the circuit elements within a new framework involving augmented voltage and current variables. The proposed solution allows solving any switching dynamical circuit that can be decomposed into the interconnection of the above two classes of elements.
(a)
(b)
5. PSL Resistive Elements
In this section, the constitutive relation of PSL resistive circuit elements and their equivalent characteristics expressed in terms of augmented timeinvariant relations are derived. For conciseness, the discussion and the proposed procedure are based on the twoterminal switching element of Figure 2. The method, however, is general and can be readily applied to the multiport case. The above extension is briefly outlined hereafter in this section and is applied to solve an example problem involving a twoport PSL resistive block in Section 7.
5.1. Constitutive Relations
In the scheme of Figure 2, the PSL element is driven by an ideal voltage source to allow the computation of its generalized admittance representation, that is, the constitutive relation involving the port voltage and the current as the input and output variables, respectively. The resistive switching element is the same of Figure 1 where the switch is characterized by a periodic switching activity between the open and the close positions with a characteristic period . The switch is assumed to be closed in the first half part of the period (e.g., ) and open in the remaining part. A cisoidal signal is chosen as the voltage excitation in this example.
Since the switching block does not include dynamical elements, the port current can be alternatively interpreted as the juxtaposition in time of the current response computed from the analysis of the circuit of Figure 2 at each operating state of the switch (i.e., open and close). Hence, the current writes where and are the equivalent conductances of the PSL twoterminal element when the switch is close or open, respectively, and is the window function defined by
The direct application of the Fourier transform to (17) leads, via a relatively long but straightforward manipulation, to where and are the coefficients of the Fourier series expansion of the periodic function :
It is important to notice that (19) is equivalent to the general equation (14), where a singlefrequency excitation is considered ( being the Fourier transform of the input cisoidal signal). Equations (14) and (19) allow defining the aliasing transfer functions as
In this example, the functions do not depend on the variable since the circuit in each operating state is resistive and its solution involves instantaneous relations only. According to (13), the generalized bifrequency admittance can be written as
The readers should refer to the appendix that collects additional details on the behavior of the example circuit of Figure 2 with the aim of further stressing the key features of PSL circuits.
It is relevant to notice that the proposed procedure can be applied to multiport resistive elements as well, provided that their port admittance characteristic can be computed as for the example twoterminal element of Figure 2. In practice, voltage sources are applied to the different ports of the multiport and the current responses are interpreted as the juxtaposition in time of the port current responses computed for the different operating states of the switches, leading to matrix relations similar to (19).
5.2. Augmented Characteristics
Equation (19) suggests that the frequencydomain steadystate response of a PSL circuit with independent voltage or current sources defined by cisoidal excitations with angular frequency can be suitably expressed in terms of an infinite sum of delta functions [11]. For practical reasons, this infinite series is then truncated. For example, for the case of the circuit of Figure 2, where are the Fourier coefficients of the th harmonic of the current through the PSL element and is the total number of (both positive and negative) harmonics considered, being the total number of terms defining the expansion (23). A similar relation holds for the voltage .
The previous interpretation of the electrical variables along with (7) (with and ) and (21) leads to the following generalized Ohm’s law describing the PSL resistive block: for .
It is important to notice that the periodic behavior of the system and the convenient interpretation of voltages and currents in terms of series expansions allow replacing the integral operator with the products of two sums. From the above equation it is also clear that voltage and current variables are defined by truncated Fourier expansions, which take into account harmonics, from to .
Equation (24), along with the representation of the current by means of (23), can be represented via the following augmented characteristic equation of the twoterminal PSL element in matrix form: where the voltage and current variables and are replaced by the vectors and . The above vectors, that have dimension , are filled in by the harmonic coefficients defining the steadystate response of the corresponding variables and play the role of the new port variables of a possibly augmented circuit element. Clearly, this is the same underlying interpretation of the socalled harmonic balance that has been used along with other methods for the solution of PSL circuits in [19–22]. In the latter papers, however, the proposed simulation technique involves more complex and cumbersome computations. In our proposal, instead, the complexity of the circuit solution is dramatically reduced by the interpretation of the characteristics of PSL blocks according to (25).
Equation (25), that turns out to be a linear timeinvariant representation of the PSL element, can be written in a more compact form as where is the matrix of (25).
For conciseness, the above procedure has been detailed for the case of a twoterminal element only. Its extension to the more general multiport case is straightforward and requires to follow the same steps above.
6. Circuit Solution via Augmented MNA
This section adopts the classical MNA approach for circuit analysis to the case of PSL circuits. For the sake of illustration, the derivation is based on the example circuit of Figure 3, where a nonideal voltage source with inductive series impedance is connected to the switching resistive block of Figure 2. In the schematic of Figure 3, the PSL block is represented by its generalized bifrequency admittance defined by means of (22) and the ideal voltage source is defined by the cisoidal excitation .
The discussion starts with the wellknown stamp arising from the analysis of the schematic of Figure 3 where a simple admittance (e.g., ) replaces the PSL block. In this case, the MNA equation in matrix form writes
When the impedance is substituted by , the unknowns, that is, , , and , need to be replaced by the corresponding vectors , , and collecting the coefficients of the harmonic series expansion of the nodal unknowns.
The augmented characteristic of the PSL twoterminal element defined by (25) and (26) and the companion relations that can be readily obtained for the classical LTI elements allow replacing the original MNA equation with where all the entries of the matrix are submatrices with dimension , where , , and . On the right hand side of the MNA equation, the excitation term is composed of vectors of dimension defined as and .
Clearly, (28) turns out to be times larger than the corresponding MNA equation of a circuit with the PSL element replaced by a linear timeinvariant impedance. However, the proposed extended matrix, that belongs to the same class of (27), can be readily solved via simple linear inversion. What is more important, the steadystate response of the PSL circuit can be computed by reinterpreting the coefficients of the new voltage and current unknowns by means of (23). Equation (29) defines the unwrapped and complete extended MNA stamp for the example circuit of Figure 3 with .
In sum, the proposed simulation method is valid for PSL circuits with both constant sources (a constant is seen as a cisoidal signal with null angular frequency) and any arbitrary periodic excitation, provided that the latter signal is expended in Fourier series. Also, the solution of the augmented MNA equation like (29) provides the direct computation of the steadystate response of the circuit.
For the sake of illustration, the method has been derived for the case of a timevarying circuit with one PSL resistive twoterminal element only. The method, however, is general and can be suitably applied to arbitrary circuits with one or more multiport resistive PSL elements, provided that the switching elements can be described by means of an admittance representation. It is also useful to remark that the proposed method is still valid if the resistive switching element is described by its impedance representation. The same procedure suggested above for the generation of the augmented characteristic (26) of the example twoterminal element of Figure 1 can be followed. As an example, Section 7 includes the results obtained by considering a circuit with a twoport switching resistive element and briefly outlines the procedure for the computation of its augmented characteristic. The solution of the proposed augmented MNA problem introduces benefits in terms of robustness and efficiency with respect to the classical ODE piecewise approach in timedomain.
7. Numerical Results
This section collects the numerical results obtained by applying the proposed method to the solution of the two example circuits. The responses of the circuits computed by means of MATLAB and the standard ordinary differential equation (ODE) integration routines are assumed as the reference curves in this study. The reference curves are compared with the solution obtained by means of the linear inversion of the augmented matrix equation like (29).
The first validation deals with the same example circuit of Figure 3. For this test case, Figures 4 and 5 show both the time and frequencydomain steadystate responses of the current and voltage , respectively, thus highlighting the good accuracy of the proposed method in reproducing the reference responses. The timedomain comparison is carried out by considering the series approximation of the unknowns with an increasing number of the expansion orders . For plot readability, a maximum order of has been considered since a larger order generates responses overlapping the reference curve very well and does not allow appreciating the differences among the approximated responses:
In order to provide a better quantitative comparison, Table 1 collects the main figures on the accuracy and efficiency of the ODEbased and of the proposed MNAbased methods implemented in MALAB environment. The accuracy is assessed by computing the relative mean square error between the reference and the predicted responses. Also, the overhead of the proposed method in computing the aliasing admittances is included in the total CPUtime. The numbers in the table confirm the strengths of the proposed approach that has a number of advantages over alternative stateoftheart methods. Mainly, it allows computing the solution of the nodal variables of a periodically switched circuit in both frequency and timedomains by means of the solution of a single linear problem. The coefficients in the unknown vector of (29) are the numbers multiplying either the complex trigonometric functions defining the timedomain responses or the deltas defining the different terms of the corresponding frequencydomain spectra.

As a second and more realistic test case, the circuit of Figure 6 is considered. It is composed of a twoport PSL resistive block (highlighted by means of a gray frame) and standard linear timeinvariant twoterminal elements. The topology of this second example is similar to the ones occurring in a number of applications (e.g., the DCDC boost converter is defined by a similar circuit with a different PSL block involving two switching components, i.e., the diode and the MOS transistor). As already done in Section 5 for the twoterminal PSL element of Figure 3, the characteristic of the twoport resistive PSL block needs to be suitably replaced by an approximated characteristic involving the expansion of its port voltages and currents (e.g., see (24) and (26)). Specifically, the integral equations governing the port behavior of the PSL block, can be suitably replaced, via the same procedure of Section 5, by where the interpretation of the variables and matrices is the same as the one introduced in the previous section. From the above equation, it is clear that the PSL block is now suitable to be plugged into the augmented MNA equation describing the circuit of Figure 6, leading to where is a submatrix with dimensions and . The remaining entries of the matrix and the nodal unknowns are the same already defined in (29).
Figure 7 shows the reference and the predicted responses of the voltage and the current of this second test circuit. The curves in the figure highlight that the proposed method produces the steadystate responses of a PSL circuit and does not include the initial transient observed in the reference response of the circuit calculated via an ODEbased solver with null initial conditions of the dynamical components. As far as the performance of the proposed method is concerned, the trend given by the numbers of Table 1 is also confirmed for this second example (see Table 2 for detailed information on the accuracy and performance of the method for an increasing number of the expansion orders ). This second example allows claiming the generality of the proposed approach that can be applied to circuits with multiterminal PSL blocks without modifying the basic rules introduced in the previous section for the generation of the augmented MNA equation.

(a) Current
(b) Voltage
8. Conclusions
This paper addressed the simulation of the steadystate behavior of an important class of timevarying circuits consisting of linear timeinvariant (possibly dynamical) circuit elements, periodically switched linear resistive multiport elements, and ideal sources defined by periodic excitations. The proposed solution extends the results of previously published papers in this field and is based on the generation of an augmented timeinvariant MNA equation governing the circuit behavior. The new nodal unknowns in the advocated MNA equation correspond to the harmonic coefficients defining the expansions of the unknown voltage and current variables. A topological approach is considered, with emphasis on an intuitive physical based interpretation of the PSL elements in the network. The proposed solution is proven to offer a modular approach to circuit analysis, leading to accurate results with good simulation speedups. Also, it is sufficiently general to handle a wide class of circuit topologies occurring in real application problems. The feasibility and strengths of the method are demonstrated on simple analytical examples and two validation test circuits.
Appendix
This appendix collects some useful remarks on the general behavior of a PSL circuit to highlight its features and to provide a graphical interpretation of the operators involved in the circuit analysis. The discussion is based on the circuit of Figure 2.
As an example, Figure 8 shows the magnitude of defined by (22) and computed for Ω and kHz. To improve the readability of the plot, the axis is discretized and four samples are considered only. This figure allows getting a practical interpretation of the transfer function of a timevarying system. Due to the term in (13), the cut of the twodimensional surface for a fixed value turns out to be given by the superposition of delta functions occurring at frequencies (). Clearly, the position of such deltas on the axis are shifted to the right (or to the left) when the values of the input frequencies increase (or decrease).
Figure 9 shows the timedomain current response of the circuit of Figure 2 to a sinusoidal excitation , kHz, and its corresponding frequencydomain spectrum. The above figure highlights that the steadystate response of a periodically timevarying system contains an infinite number of components at frequencies .
(a)
(b)
The example discussed in this appendix clearly illustrates that a linear timevarying circuit excited by a sinusoidal signal produces new harmonics even if the circuit does not contain active elements (i.e., the circuit of Figure 2 has an ideal voltage source connected to a passive twoterminal switching resistive element).
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 L. A. Zadeh, “Frequency analysis of variable networks,” Proceedings of the IRE, vol. 38, no. 3, pp. 291–299, 1950. View at: Google Scholar
 L. A. Zadeh and C. A. Desoer, Linear System Theory: The State Space Approach, McGrawHill, New York, NY, USA, 1963.
 A. Papoulis, Signal Analysis, McGrawHill, Auckland, New Zealand, 1977.
 M.L. Liou, “Exact analysis of linear circuits containing periodically operatd wwthes with applications,” IEEE Transactions on Circuit Theory, vol. 19, no. 2, pp. 146–154, 1972. View at: Google Scholar
 T. Strom and S. Signell, “Analysis of periodically switched linear circuits,” IEEE Transactions on Circuits and Systems, vol. 24, no. 10, pp. 531–541, 1977. View at: Publisher Site  Google Scholar  MathSciNet
 F. Yuan and A. Opal, Computer Methods for Analysis of MixedMode Switching Circuits, Kluwer Academic, New York, NY, USA, 2004.
 F. Yuan and A. Opal, “Sensitivity analysis of periodically switched linear circuits using an adjoint network technique,” in Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS '99), pp. V331–V334, June 1999. View at: Google Scholar
 J. A. Á. Martín, J. R. Melgoza, and J. J. R. Pasaye, “Exact steady state analysis in power converters using Floquet decomposition,” in Proceedings of the North American Power Symposium (NAPS '11), pp. 1–7, Boston, Mass, USA, August 2011. View at: Publisher Site  Google Scholar
 J. Kovar, Z. Kolka, and D. Biolek, “Comparison of averaging and harmonic balance methods for switched DCDC converters,” in Proceedings of the 17th International Conference on Mixed Design of Integrated Circuits and Systems (MIXDES '10), pp. 402–407, Warsaw, Poland, June 2010. View at: Google Scholar
 H. Behjati, L. Niu, A. Davoudi, and P. L. Chapman, “Alternative timeinvariant multifrequency modeling of PWM DCDC converters,” IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 60, no. 11, pp. 3069–3079, 2013. View at: Publisher Site  Google Scholar
 T. A. C. M. Claasen and W. F. G. Mecklenbrauker, “On stationary linear timevarying systems,” IEEE Transactions on Circuits and Systems, vol. 29, no. 3, pp. 169–184, 1982. View at: Publisher Site  Google Scholar
 H. Sandberg, E. Mollerstedt, and Bernhardsson, “Frequencydomain analysis of linear timeperiodic systems,” IEEE Transactions on Automatic Control, vol. 50, no. 12, pp. 1971–1983, 2005. View at: Publisher Site  Google Scholar
 P. M. Mellacheruvu, S. M. Mahajan, C. L. Carnal, and J. J. Biernacki, “Application of fourier and laplace transform techniques for modeling an electric circuit having timevarying components,” in Proceedings of the 38th Southeastern Symposium on System Theory (SSST '06), pp. 147–151, March 2006. View at: Google Scholar
 J. Lerdworatawe and W. Namgoong, “Generalized linear periodic timevarying analysis for noise reduction in an active mixer,” IEEE Journal of SolidState Circuits, vol. 42, no. 6, pp. 1339–1351, 2007. View at: Publisher Site  Google Scholar
 M. L. Liou, Y.L. Kuo, and C. F. Lee, “A tutorial on computeraided analysis of switchedcapacitor circuits,” Proceedings of the IEEE, vol. 17, no. 8, pp. 987–1005, 1983. View at: Google Scholar
 R. Trinchero, I. S. Stievano, and F. G. Canavero, “Steadystate analysis of switching power converters via augmented timeinvariant equivalents'',” IEEE Transactions on Power Electronics, vol. 29, no. 11, pp. 5657–5661, 2014. View at: Publisher Site  Google Scholar
 H. Bishnoi, P. Mattavelli, R. Burgos, and D. Boroyevich, “EMI behavioral models of DCfed threephase motor drive systems,” IEEE Transactions on Power Electronics, vol. 29, no. 9, pp. 4633–4645, 2014. View at: Publisher Site  Google Scholar
 I. Cadirci, B. Saka, and Y. Eristiren, “Practical EMIfilterdesign procedure for highpower highfrequency SMPS according to MILSTD 461,” IEE Proceedings—Electric Power Applications, vol. 152, no. 4, pp. 775–782, 2005. View at: Google Scholar
 M. Okumura, H. Tanimoto, T. Itakura, and T. Sugawara, “Numerical noise analysis for nonlinear circuits with a periodic large signal excitation including cyclostationary noise sources,” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 40, no. 9, pp. 581–590, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Dobrovolny, G. Vandersteen, P. Wambacq, and S. Donnay, “Analysis and whitebox modeling of weakly nonlinear timevarying circuits,” in Proceedings of the Design, Automation and Test in Europe Conference and Exhibition, pp. 624–629, 2003. View at: Publisher Site  Google Scholar
 J. Roychowdhury, “Reducedorder modeling of timevarying systems,” IEEE Transactions on Circuits and Systems II: Analog and Digital Signal Processing, vol. 46, no. 10, pp. 1273–1288, 1999. View at: Publisher Site  Google Scholar
 T. Mei, H. Thornquist, E. Keiter, and S. Hutchinson, “Structure preserving reducedorder modeling of linear periodic timevarying systems,” in Proceedings of the IEEE/ACM International Conference on ComputerAided Design (ICCAD '11), pp. 361–366, San Jose, Calif, USA, November 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Riccardo Trinchero et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.