Nonlinear Dynamics and Stability Analysis of a Three-Cell Flying Capacitor DC-DC Converter

: This paper presents a study of the nonlinear dynamic behavior a ﬂying capacitor four-level three-cell DC-DC buck converter. Its stability analysis is performed and its stability boundaries is determined in the multi-dimensional paramertic space. First, the switched model of the converter is presented. Then, a discrete-time controller for the converter is proposed. The controller is is responsible for both balancing the ﬂying capacitor voltages from one hand and for output current regulation. Simulation results from the switched model of the converter under the proposed controller are presented. The results show that the system may undergo bifurcation phenomena and period doubling route to chaos when some system parameters are varied. One-dimensional bifurcation diagrams are computed and used to explore the possible dynamical behavior of the system. By using Floquet theory and Filippov method to derive the monodromy matrix, the bifurcation behavior observed in the converter is accurately predicted. Based on justiﬁed and realistic approximations of the system state variables waveforms, simple and accurate expressions for these steady-state values and the monodromy matrix are derived and validated. The simple expression of the steady-state operation and the monodromy matrix allow to analytically predict the onset of instability in the system and the stability region in the parametric space is determined. Numerical simulations from the exact switched model validate the theoretical predictions.


Introduction
Power electronics converters have been widely used for industrial applications, such as motor speed variation, renewable energy technologies, and harmonic filtering [1][2][3], to adjust input voltages and currents to obtain regulated power supplies. Their structures achieve this regulation by controlling the flow of energy in the reactive storage elements and by adjusting the conduction time of the switching devices. Elementary DC-DC switching converters such as the buck, boost and buck-boost are the most commonly used. These systems include mainly switching devices (transistors and diodes), energy storage elements (inductors and capacitors) and appropriate feedback loops. Despite their apparently simple structures, numerical simulations and experimental measurements have shown that they can exhibit a plethora of nonlinear dynamic phenomena like bifurcations and chaos. The appearing of these phenomena is mainly attributed to the nonlinearity induced by the switching devices and the feedback loop [4]. These phenomena have been studied in the past by many researchers [5][6][7][8][9][10][11][12][13][14][15][16][17]. A review of the different kinds of nonlinear behaviors that could take place in such systems can be found in [18][19][20][21].
In power electronic converters, these kinds of nonlinear behaviors should be avoided and suppressed and therefore a deep understanding of the mechanism of their occurrence is necessary [22]. To analyze these behaviors, an appropriate and accurate model is required. The conventional approach used for modeling of this kind of switching systems is based on the averaging technique [1,23]. It is true that the average model reduces system complexity and can predict the low frequency behavior, but it was demonstrated that it fails to predict the nonlinear behavior of the system that could taking place at the frequency ranges close to one half switching frequency. This behavior is known in the literature as period doubling bifurcation or fast-scale instability in the nonlinear dynamics field and subharmonic oscillation in the power electronics field. With the averaging procedure, it is assumed that the steady-state behavior of the system is an equilibrium point while in reality it is a limit cycle whose stability analysis is mathematically more involved. To predict the actual behavior of the system, the exact switched model or the related discretetime model must be used. There are many works dealing with the nonlinear dynamics of switching converters most of them use a discrete-time model. While obtaining such a model for elementary switching converters is a simple task, it is not the case for advanced converter topologies such as multi-level flying capacitor converters mainly because the number of circuit configurations among which the circuit toggles is high [4]. To perform the stability analysis of the desired limit cycles of the switched model, Floquet theory can be applied without passing through the discrete-time model [24,25]. Therefore, it offers an efficient and accurate tool to perform stability analysis leading to the same results as the discrete-time modeling approach but more straightforwardly. It also separately reveals the effects of each switching on the system stability [25]. Floquet approach is based on the theory originally developed by Floquet [26] to study the stability of periodic orbits in dynamical systems with time varying periodic coefficients. These approach has been extended by Aizerman [27], and later by Filippov [28] to analyze impacting motion and stick-slip oscillations in mechanical switching systems. A good tutorial on how to use this theory is in [24].
On the other hand, flying capacitor multi-level multi-cell converters have been proposed for higher power industrial applications. These converters are based on adding and subtracting voltage of storage elements in the form of capacitors creating multi-level voltages [2]. A special characteristic of flying capacitor converters are their commutation cells that allow to operate with high supply voltages making them good candidate for high voltage applications [2]. They can be also used to reduce the stress in the switching devices and the state variables ripples by interleaving operation. Their control objective is two-fold. First, the controller must ensure a regulated output voltage/current. Second, the control circuit must ensure a voltage balancing between the different cells of the converter. This is accomplished just as in conventional switching converters by applying appropriate signals to control the time durations of the ON/OFF state of the switches. This is accomplished by using suitable error signals that represents the difference between the reference and the actual value of the state variable to be controlled. Selecting the right value of the feedback gain of these controller is crucial for ensuring stability of the system and suppressing underscored nonlinear behavior. The feedback control may be designed in the continuous time or in the discrete time domain. Discrete-time controllers have been proposed recently as substitutes of continuous-time ones due to ease in programming, insensibility to noise and the possibility of implementing complex algorithms.
The flying capacitor multi-cell mult-level converters have attracted a lot of attention from the power electronics and control communities [29][30][31][32][33][34][35][36][37][38] since their appearance in [2,3]. The reasons behind the interest for these topologies are modularity, applicability to high voltage fields and ease of balancing the flying capacitor voltage and possibility to use small reactive components [36] which could reduce the size and increase the power density. Several studies have been conducted to solve the flying capacitor balancing issues [35]. A sensorless stabilization technique is presented in [37] where a discrete-time digital platform based on FPGA was used. The use of discrete-time control techniques in the field of power electronics converters is more and more increasing. A discrete-time digital predictive control for a three-level flying capacitor buck converter with voltage balancing is presented in [38]. Discrete-time PWM control has been applied to a two-cell DC-DC converters [29][30][31][32]39]. This procedure can be adapted to any multi-cell converter with a generic number of cells [40].
Floquet theory has been successfully applied in the past to perform stability analysis of simple power electronics converters such as the conventional buck converter [25], cascaded converters [33] and interleaved converter [41] among others. The application of this approach to multi-level flying capacitor converters has not been reported for the best knowledge of the authors. It is the aim of this paper to apply this theory to predict the nonlinear dynamics and to perform stability analysis of a three-cell flying capacitor DC-DC buck converter under discrete-time controller. First, the system modeling is addressed and the nonlinear switched model is derived. Second, the different operating modes are determined according the value of the operating steady-state duty cycle. Then, a discrete time controller is proposed which will be in charge of regulating the output current and balancing the cell voltages. Subsequently, the limit cycles of the converter will be computed and the result will be compared with numerical simulations using the exact switched model. Floquet theory is applied which leading to an accurate prediction of the instability phenomena exhibited by the system. However with the exact expression of the monodromy matrix, the computational efforts are high. Therefore, some realistic approximations are used to obtain the limit cycles and the monodromy matrix analytically predicting the onset of instability in the three-cell DC-DC buck converter with a much lower computational effort. The rest of the paper is organized as follows: in Section 2, a system description and the operating modes are presented. The mathematical modeling of the power stage and the discrete-time controller are also be presented in the same section. The switched model and the switching conditions are derived. The system nonlinear dynamics is explored in Section 3 using bifurcation diagrams showing that the system may exhibit period doubling bifurcation, subharmonics and chaotic regimes when suitable system parameters are varied. Floquet theory for stability of limit cycles is presented in Section 4 and adapted for the system under study. A procedure to compute limit cycles and to analyze their stability is presented. The monodromy matrix and Floquet multipliers explaining the observed period doubling bifurcation phenomena in Section 3. A simplified approach is proposed to derive a simplified expression for the monodromy matrix and Floquet multipliers in the same section where also a validation of the theoretical predictions is provided. Finally, in Section 5, the concluding remarks of this study are given.

Power Stage Circuit
The studied converter is shown in Figure 1. It is based on the conventional single DC-DC buck converter but it is enriched with additional cells. The voltage stress on the switches for the conventional two-level single-cell buck converter is the total input voltage v g . The number of cells determines the number of voltage levels that the converter can handle. Therefore, for the 3-cell 4-level buck converter, the voltage stress on the switches is v g /3. The size of the inductance can be reduced by 3 2 = 9 while getting the same inductor current ripple as a conventional buck converter operating with the same switching frequency. If the inductance value is maintained the same as the one used in a conventional buck converter, the current ripple will be reduced by 3 2 = 9 in the 3-cell buck converter operating with the same switching frequency. Therefore, compared to a two-level single-cell buck converter for the same current ripple, the multi-cell buck converters have better power density because of using a smaller inductor. Figure 2 shows an illustrative inductor current waveforms for the conventional buck converter with a switching period T and inductance L and two different inductor choice for the 3-cell buck converter.

Load current error
Gain selection

Time
Switching period T Conventional buck with L 3-cell buck converter with L Each pair of switches S i and S i (i = 1, 2, 3) indicates a cell. Switches S i and diodes S i are activated in a complementary manner. When S i is ON , S i is OFF and viceversa. The voltages in the cells must be balanced to reduce the stress in these devices [2,3]. The presence of switching devices in power electronic converters produce a cyclic change of structure when defined conditions of time and variables are fulfilled. The different switch S i positions define the different circuit configurations C u 1 u 2 u 3 among which the converter toggles depending on the values of u 1 , u 2 and u 3 . The triple subscript notation that was adopted in this study correspond to the state of the different switches that is in turns related to the values the command signals u i . Namely, u i defines the state of the switches S i for i = 1, 2, 3. If u i = 1, S i is closed (ON) and if u i = 0, S i is open (OFF). For instance the state (1,1,1) for the switches S 1 , S 2 and S 3 defines configuration C 111 for which all the switches S i (i = 1, 2, 3) are closed and the state (0, 0, 0) for the same switches defines configuration C 000 for whic all the switches are open. For optimization of the state variable ripples, the driving signals are phase shifted 2π/3 [2]. The switching from one configuration to the following forms a fixed sequence as according to the operation modes. Each configuration C u 1 u 2 u 3 is described by an affine time invariant differential equations in the form ofẋ = A u 1 u 2 u 3 x + B u 1 u 2 u 3 . A u 1 u 2 u 3 and B u 1 u 2 u 3 are the system matrices and vectors during each configuration C u 1 u 2 u 3 for all the possible combinations of the values of u i for i = 1, 2, 3 for a specific mode of operation. x is the vector of state variables i L , v 1 and v 2 . The mathematical model can be obtained by applying standard Kirchhoff's voltage and current laws hence obtaining the following set of continuous-time differential equations modeling the dynamics of the converter d dt It is worth noting here that the switched model is piecewise affine and for each different configuration, the corresponding differential equation can be solved in closed form.

Operating Mode under Study
The three-cell buck converter has different operating modes that depend on the value of the steady-state duty cycle D of the switches S i . It should be noted here that D is the same for all the switches S i . This study was carried out with the condition 1 3 < D < 2 3 . The same procedures can be followed for other ranges of D. The driving signals sequence is given by (u 1 , u 2 , u 3 ) = (1, 0, 1), (1, 0, 0), (1, 1, 0), (0, 1, 0), (0, 1, 1), (0, 0, 1) as depicted in Figure 3. It develops a periodic orbit with six configurations per switching cycle. The sequence is shown as follows Switching period T The current paths in each circuit configuration are highlighted in Figure 4. The matrices A u 1 u 2 u 3 , the vectors B u 1 u 2 u 3 for each topology and a brief description of the behavior in the converter elements for each case are shown in Table 1.

Discrete-Time Controller
Power electronic converters are generally used to provide a desired voltage/current for the load. This can be accomplished through suitable feedback loops. For flying capacitor multi-cell converters, the controller task it to regulate the voltage/current of the load and also to balance the flying capacitor voltages to effectively reduce the switches voltage stresses. It can be demonstrated that if the duty cycles for all the cells is the same, a natural balancing phenomenon takes place without using the flying capacitor voltages in the feedback loop [34]. However, the transient response speed in this case depends on the circuit parameters and could not undesirably very slow. To avoid this problem, the voltages of the cell are also used in the feedback loop. The controller must be able to regulate the current i L and balance the input voltage between the different cells. The expression of the duty cycles in terms of the state variables error can be fixed according to the control objectives. A purely proportional control will be used here with the aim to minimize the error between the current and voltages and their references values. As the system contains three controlled switches, three different duty cycles must be generated and their corresponding driving signals must be applied for the three switches. The error between the controlled outputs and their references are appropriately combined and used in a discrete time modulation strategy. In particular, the inductor current (which is also the load current) i L is controlled to a desired current reference i ref and the flying capacitor voltages v 1 and v 2 are controlled to 2/3 and 1/3 of the input voltage respectively. The charge and discharge of the capacitors depend on the flow direction of the current across it. The discrete-time duty cycles are given by where κ i , κ 1  are defined as the ratio between the fraction of the switching period T during which the switch S i is closed (ON state) and this period and therefore they must be constrained between 0 and 1. This means that they can be saturated to 0 or to 1 if the expressions in (3)-(5) give a value outside the interval (0, 1). Therefore, a saturation function must be used to limit the value of d i [n] in this interval. The final expressions for the duty cycles are where sat(·) is the saturation function defined as follows The driving signals u i are appropriately phase-shifted in such a way that the signal u 1 is high and the switch S 1 is closed (ON) during the time interval (nT, (n + d 1 [n])T), the signal u 2 is high and the switch S 2 is closed (ON) during the time interval (n +

Nonlinear Behavior from Time-Domain Numerical Simulations
Our concern in this section is to present some graphic results to support and complement the study. The fixed parameter values used in this study are depicted in Table 2. The parameters κ i and i ref are varied. A powerful tool to have a wide panoramic view of the different dynamical behavior that the system can have when parameters are varied is a bifurcation diagram. In order to explore the dynamics of the system, bifurcation diagrams for the system are plotted by considering the current feedback gain κ i as bifurcation parameter which is varied within the range κ i ∈ (0.04, 0.  Figure 5 as an example. From this figure, it is possible to see that the system exhibits a stable periodic behavior for values of κ i less than a critical value κ i,cri ≈ 0.058 A −1 . At this critical value of the feedback gain, the system period-1 limit cycle loses its stability and the attractor of the system evolves on a period-2 limit cycle, i.e., the system undergoes a period doubling bifurcation and a period-2 limit cycle emerges. By further increasing the feedback gain κ i , the resulting period-2 limit cycle undergoes a period doubling resulting in a period-4 limit cycle. For higher values of κ i , the system undergoes more bifurcations phenomena leading eventually to chaotic behavior. Nonlinear behavior manifested by different periodicity parametric windows and chaos can also be observed.  In order to check the results obtained from the previous bifurcation diagrams, timedomain simulations have been carried out by using the exact switched circuit-level model implemented in PSIM© software. The steady-state time-domain waveforms are depicted in Figure 6 for different values of κ i . It is worth noting that in the top left panel of Figure 6, the flying capacitor voltages oscillate at the switching frequency f s = 40 kHz, the inductor current oscillate at the double of the switching frequency 2 f s = 80 kHz due to interleaving effect and all the duty cycles are the same and constant in steady-state. In the top right panel of Figure 6, the flying capacitor voltages oscillate at one half the switching frequency f s 2 = 20 kHz, the inductor current oscillate at the switching frequency f s = 40 kHz and all the duty cycles are jumping up and down between two different values one of them is saturated to its lower level 0. In the bottom left panel of Figure 6, the flying capacitor voltages oscillate at one fourth the switching frequency  Figure 6, all the state variables and duty cycles are chaotic and saturated upper and the lower values of the duty cycles are reached. It is also worth to note the increase in the current ripple amplitude after losing the stability of the fundamental period-1 limit cycle which could have a damaging effect on the system efficiency and stresses in the switching devices and which also could shorten their lifetime. Therefore, the stability analysis of this limit cycle is of concern both from a theoretical and from a practical point of view.  Floquet theory for limit cycle stability analysis will be applied to predict the onset of instability in the three-cell buck converter. The theoretical background of this theory adapted for the three-cell buck converter is presented in the following section.

Limit Cycles Computation
In steady state, for 1 3 < D < 2 3 , there will be six time instants at which a switching from a configuration C u 1 u 2 u 3 to another will take place (see Figure 3). These are t 0 = nT, t 2 = (n + 1 3 )T and t 4 = (n + 2 3 )T. The other remaining switching instants are t 1 = n + (D − 2 3 )T, t 3 = (n + D)T and t 5 = n + (D − 1 3 )T. The steady-state duty cycle D can be obtained by imposing volt-second balance on the inductor current [1]. This leads to D = Ri ref /v g . If 1 3 < D < 2 3 , during one switching cycle, the three-cell buck converter changes its structure among the six different configurations shown in Figure 4 and the the system will have a limit cycle using six different affine models. In steady-state, the limit cycle x(t) starts with initial conditions x(0), at the start of the switching cycle and evolves according to the switching sequence in (2) as shown below where Φ u 1 u 2 u 3 and Ψ u 1 u 2 u 3 are given by It is worth to note that the limit cycle can be determined once its initial state x(0) is obtained. This initial state x(0) can be computed by imposing T-periodicity hence getting where Φ and Ψ are given by the following expressions The stability properties for the limit cycle represented by the vector x(0) can be revealed by using the monodromy matrix which will be determined below.

Accurate Stability Analysis Using Floquet Theory
As explained before, the switched model of the system under study is time periodic with the switching period T determining the periodicity of the desired limit cycle whose stability can be performed by using Floquet theory combined with Filippov method [24,25]. The main tool for studying the stability of limit cycles using this theory is the monodromy matrix M. This matrix is defined such that the dynamics close to a limit cycle can be expressed as followsx (t + T) = Mx(t) ∀t (10) where the overhat stands for small signal variations. The eigenvalues of M are called Floquet multipliers and they determine the amount of contraction or expansion close to a trajectory and therefore they can be used to determine the local stability of limit cycles. The limit cycle will be stable if all the Floquet multipliers have modulus less than 1. For piecewice affine systems like the power converter considered in this study, the monodromy matrix can be constructed from the product of the state transition matrices corresponding to each sub-interval and the corresponding saltation matrix [24][25][26][27][28]. For the case of the three-cell buck converter of this study, the monodromy matrix M can be written as the product of the state transition matrices Φ u 1 u 2 u 3 and saltation matrices at switching instants between their respective time intervals. Namely, the monodromy matrix can be expressed as follows Figure 7 shows an illustration of how to construct the monodromy matrix using the state transition and the saltation matrices. The state transition matrices from one switching instant to another are just exponential matrices of the different A u 1 u 2 u 3 matrices in the piecewise affine modelẋ = A u 1 u 2 u 3 x + B u 1 u 2 u 3 . The saltation matrices can be obtained by first determining the conditions governing the switching in the previous model. These switching conditions are given by the following expressions where ρ is the vector of the reference signals that can be expressed by and K 1 , K 3 and K are the gradients corresponding to the switching functions h i , i = 1, 3, 5. The saltation matrices S 001→101 , S 101→100 , S 100→110 , S 110→010 , S 010→011 and S 011→001 can be obtained using the Filippov method [24,25,28] and their expression are where I is a unitary matrix with appropriate dimension. As stated before, due to interleaving effect, the three switching instants 0, t 2 = (D + 1 3 )T and t 4 = (D + 2 3 )T are given in a fixed pattern and therefore one has x(t) x(t + T ) x(τ 1 ) This implies that Therefore, expression (11) for the monodromy matrix becomes (29) To obtain the gradients ∂h i (x, t)/∂t for i = 1, 3, 5, the expression of x(0) must be obtained in terms of t for t 5 ≤ t < T. To get this expression, let us write x(t) for t 5 ≤ t < T in terms of x(0) Therefore, the sampled state x(0) can be expressed as follows where matrix η is given by With the expression of x(0) in terms x(t), the gradient vectors become as follows Finally, the partial time derivative ∂h i (x, t)/∂t for i = 1, 3, 5 are With the previous expressions for all the terms needed for constructing the monodromy matrix, this can be derived and the stability analysis of the limit cycles can be performed. It is worth noting here that the previous theoretical approach can be applied to any kind of limit cycles. These refer to the different periodic orbits under steady state operation. This section considers the period-1 limit cycle which is the only desired behavior for this kind of converters in any industrial application.
On the other hand, although the obtained monodromy matrix is accurate for predicting the onset of instability, it is far from being simple and cannot be applied directly for system design. In the coming subsection an approximate expression for this matrix will be derived. It will allow showing insight into how the different system parameters influence the system dynamics.

Approximate Expression of the Monodromy Matrix and Stability Region in the Parametric Space
The computational efforts with the previous approach is very high and it also precludes obtaining simple expressions for the stability boundary in the parametric space. In this section, some realistic assumptions will be made in order to obtain simple expressions for the stability borders. These expressions will be validated from numerical simulations using the exact semi-analytical approach detailed previously. A comparison will be performed between the results obtained from the exact approach and those obtained from the simplified approximate approach. The purpose of this comparison is to give a further support to the feasibility of the using the approximate expression of the monodromy matrix for design purposes.
In practice, the switching period is selected much smaller than the time constants of the different converter configurations. Under this condition, it is possible to approximate the matrix exponential as follows e At ≈ I + At. By using this simplified expression, we can obtain an approximate expression for the limit cycles, saltation matrices and finally monodromy matrix. The approximate expression of the limit cycle evaluated at the start of the switching period is given by In order to make a validation of the simplified expression of the x(0) let us compare the result from both the exact and the approximate approach. Namely, the validity of (39) is first demonstrated before obtaining the approximate expression of the monodromy matrix. Figure 8 shows the evolution of the voltage and the current coordinates of x(0) as function of the feedback gain κ i for two different values of the current reference i ref using the exact expression (9) and the approximate one (39). A remarkable agreement can be clearly observed which demonstrates that for design-purposes, the approximate expression (39) can be used without having to use the state exact transition matrices whose computation could be very time consuming in a repetitive simulation process. This also indicates that very little error expected if the stability analysis is performed using an approximate expression for the monodromy matrix.
With the approximated vector x(0) given in (39) and the simplified expressions of the state transition and the saltation matrices, the expressions of the monodromy matrix evaluated at the previous limit cycle can be explicitly expressed as follows  Table 3 shows a comparative analysis of the results from the exact and the approximate monodromy matrix. This table shows the limit cycle represented by their initial values x(0), the corresponding Floquet multiplies, the exact and the approximate monodromy matrix M and the stability of limit cycles. The values of the Floquet multipliers are obtained by solving the characteristic equation det(M − λI) = 0 for both cases. As can be observed, the matching between the results is remarkable.  Figure 9 depicts the Floquet multipliers loci when the parameter κ i is varied. The loci was obtained from the exact and the approximate expression of the monodromy matrix. The loci from both the approximate and the exact expressions of the monodromy matrix shows that the critical value of the feedback gain at which one Floquet multipliers crosses the unit circle is 0.058 which in a perfect agreement with the numerical simulations performed from the switched model and presented in Section 3. For 0.04 < κ i < 0.058 (A −1 ), all the Floquet multipliers are inside the unit disk hence indicating that the limit cycle is stable. However, at κ i ≈ 0.058 A −1 , one Floquet multiplier crosses the unit circle from the point (−1,0) on the complex plane and the system undergoes a flip bifurcation. For κ i > 0.058, one of the Floquet multiplier is outside the unit disk being real and negative, the period-1 limit cycle becomes unstable. This explains subhamronic oscillations exhibited by the time waveforms and the bifurcation diagrams in Section 3.   Now that an approximate expression for the monodormy matrix is available in closed form, the Floquet multipliers can be found analytically. The approximate expressions for the Floquet multipliers in terms of system parameters are given by The approximate closed-form expressions of the Floquet multipliers can be used to determine the stability region of the desired limit cycle in the parameter space. A sufficient condition for stability is that all the Floquet multipliers µ i lie inside the unitary circle. For practical design reasons, the feedback coefficients are restricted within R + . These limit cycles will be stable whenever the following conditions are fulfilled for the feedback gains κ i , κ 1 and κ 2 The first constraint (44) when evaluated for the set of parameter values of Table 2 gives a critical value of the feedback gain κ i ≈ 0.058 A −1 which is in perfect agreement with the numerical simulations in Section 3. Expressions (44)-(46) can be used to get 2-dimensional diagrams that would help in setting the design limits for the parameters since as they show the safe range of the parameters that would lead to a stable system. The stability boundaries can be plotted in a suitable projection of the parametric space. The parameters can be selected in this space for obtaining a desired response for the system under step type change of a certain reference. For instance, based on a small-signal analysis, the feedback gains that leads to a deadbeat response [42], characterized by the fact that all Floquet multipliers are zero, are given by the following expression The intersection of the different curves representing (47)-(49) would lead to deadbeat response. Figure 10 shows the constraints (44), (47) and (48) in the (κ i , κ 1 ) parameter space. The stable and the unstable region as well as the set of parameter values leading to a deadbeat response are shown in the same figure. The system parameters are selected on the intersection point, and since the system is of order three, the system will settle within three switching periods after a disturbance takes place. The stable and the unstable region as well as the set of parameter values leading to a deadbeat response are shown in the same figure. The left panel of Figure 11 shows the Floquet multipliers locus when the parameters are selected at the deadbeat response point in Figure 10. These multipliers are obtained from both the exact and the approximate monodromy matrix. As can predicted by the theoretical analysis, all the multipliers are located at the origin of the complex plane. The matching between the results in this particular case is also remarkable. The right panel of Figure 11 shows the corresponding system response to a 5% step change in the input voltage v g with the feedback. It can be observed that the duration of the transient after a perturbation in the input voltage takes place is three periods agreed with the theoretical prediction.  Figure 11. Deadbeat response of the system. Left: The Floquet multipliers are at the origin of the comlpex plane. Right: The system settling time is three switching periods.

Conclusions
Power electronics converters are intrinsically complex nonlinear dynamical systems. Under parameter changes, undesired nonlinear behavior in the form bifurcations and subharmonic oscillation may take place. Understanding the nonlinear behavior and determining the effect of parameters in these systems could help in their design process. The first challenge to understand and to study these systems is obtaining an accurate model that can predict the behavior of such complex systems. Like other conventional power converters, multi-cell flying capacitor converters may also exhibit undesired nonlinear behavior in the form of bifurcation phenomena that could take place for some parameter values. First an accurate model has been used for numerically simulating the system behavior and predicting the nonlinear behavior of a three-cell flying capacitor buck converter under a discrete-time controller. In particular, bifurcation diagrams for the system under this controller has been computed. Floquet theory is a powerful tool to analyze and predict nonlinear behavior in this kind of systems. We have illustrated the use of the theory by studying a three-cell flying capacitor buck converter under discrete-time controller. First we have analyzed the switched nonlinear model of the system. Using Floquet theory, the expression of the monodromy matrix has been used to reveal the effect of suitable parameters on the system. However, the exact expression of this matrix gives little insight out the impact of each system parameter. Under some conditions, the monodromy matrix can be approximated by a simple expression whose validity has been demonstrated by numerical simulations. The simplified expression was used leading to very similar results than the ones obtained by the exact expression of the monodromy matrix. In particular, the critical value of parameters at which instability takes place using the exact monodromy matrix model were shown to be very similar to the ones obtained by an approximate expression of this matrix. The advantage of this latter is that it hows insight into how the different system parameters influence its dynamical behavior. We used the expression of this matrix to determine analytically the limit cycles, their stability status and the onset of their bifurcation. Using the simplified expression of the monodromy matrix, the stability boundaries can be easily plotted in terms of suitable parameters such as the feedback gains, input voltage, load resistance, switching frequency, filter inductance and flying capacitor capacitance values. Finally, a deadbeat-based tuning of the controller gains was used that has lead to a response with minimum settling time after a step change in system parameters. Therefore, the deep understanding of the nonlinear dynamics of the converters has allowed a design of the system with an optimum dynamic response as illustrated by simulation results from the nonlinear switched model of the converter.