A New Exact Mathematical Approach for Studying Bifurcation in DCM Operated dc-dc Switching Converters

A bifurcation study for dc-dc converters operated in DCM is performed using an accurate method. When applying classical techniques significant difficulties are encountered in the calculations. For example, using the averaging method the validity of the result is limited to half the switching frequency and higher order effects are neglected Another approach is to perform a Taylor expansion of the state transition matrices. However, this is somehow also an averaging but the fact that the Taylor series is truncated leads to unacceptable inaccuracy. A new mathematical technique for discontinuous conduction mode (DCM) analysis of dc-dc switching converters is proposed in order to predict bifurcation and chaos. The proposed technique is based on exact calculation of the state transition matrices and of the Jacobian thus providing higher accuracy of the results compared to other previously reported approaches. Beside the fact the new technique allows for exact diagnosis of instability, it is also highly general, in the sense that it can be applied to any dc-dc DCM operated converter employing any type of control. The good agreement between theoretical, simulation and experimental results, with an error lower than 0.94%, confirms the validity of the proposed method.


Introduction
Power electronic switching converters are highly nonlinear systems that give arise of various phenomena when a feedback loop is closed for control.In the past years much attention was dedicated to chaos and bifurcations analysis of switching converters with emphasis on dc-dc switching topologies employing different control strategies [1,2].
As known, dc-dc switching converters are mainly used in two operating modes: continuous conduction mode (CCM) and discontinuous conduction mode (DCM).In CCM the diode currents in the off state never reach zero and this mode occurs when the inductors are relatively high, while the load current and the switching frequency are also significant.At light loads, that is at low output currents or, equivalently, at higher load resistances, small inductor values or low switching frequencies the inductor current ripple becomes comparable to its dc value and the diode current, equal to a linear combination of inductor currents, is zero for a certain time interval during transistor off time.
In [3], a current-mode controlled boost converter was subjected to study, while a sampled-data modelling and analysis of PWM dc-dc converters in both CCM and DCM is carried out in [4].A unified analysis of CCM converters employing different control techniques (e.g., one cycle control, charge control, etc.) is performed in [5][6][7][8] from a bifurcation point of view.The asynchronous switching map on studying the stability effects of circuit parameters in CCM fixed off-time controlled buck converter is performed in [9].Bifurcations due to on-load parameter fluctuations that can occur at the boundary between CCM and DCM and discontinuity mapping were investigated for a boost converter in [10], developing the theory of border collision bifurcations for this special case.
From a chaotic point of view, in [11] a DCM voltage-controlled buck converter was presented.Some important studies regarding bifurcations in DCM were reported in [12,13].A DCM buck converter with proportional control is provided in [14].In [15] a current-mode pulse train (CMPT) control technique is used with a DCM operated converter, the effects of circuit parameters on the converter dynamics are analyzed and the border collision mechanism is revealed.The first bifurcation study reported on a DCM Ćuk converter employing current mode control was performed in [16] and the onset of oscillatory inductor currents is explained, proving that oscillatory inductor currents are inherent in the Ćuk converter, even if the input voltage is a rectified one or a pure dc.
In all these studies simplified models for describing converter behavior and deriving the iterative maps are used, mainly based on averaged models [2,11].As known, these models exhibit a limited validity in the frequency domain, the upper limit being half of the switching frequency.Another approach used in [1] expands the state transition matrices in a Taylor series and only the first three terms are retained.This is somehow also an averaging that comes with its limitations.On the other side, individual studies for the buck, boost and buck-boost are performed, although all these three converters have common roots.Therefore a tool unifying the analysis for all converters would be both powerful and expected.These models can be useful in a certain range of converter parameters but outside these limits the errors may become significant.
This paper proposes a new mathematical approach that, opposed to the above mentioned ones, allows for exact iterative map derivation and exact operating point solution calculation, thus offering the exact value for the parameters that lead to bifurcation.Moreover, the applicability of the method has a high degree of generality, in the sense that it can be applied with the same formalism to any dc-dc converter.Moreover, any control technique can be employed with the converter without modifying the algorithm proposed for the analysis.Any bifurcation parameter can be studied e.g., gain k if the controller, input voltage V g , or load resistance R. The proposed algorithm can be easily implemented in Matlab™ or similar programs.

State-Space Description of the Discontinuous Conduction Mode
In the most general case, DCM is defined as the operation mode in which at least one asynchronous switching with respect to the control signals occurs.The passive switches, namely the diodes, are responsible for this type of operation.The main waveforms of a converter in DCM are depicted in Figure 1, where q(t) is the transistor switching function, d n is the duty cycle in the nth period, h n T is the length of the second topological state when the active switch is off and the diode is forward biased, and the third topological state is when both semiconductors are off.T is the switching period.Compared to CCM, in DCM a third topological state is introduced, corresponding to the hatched zone in Figure 1.Obviously, from Figure 1 it results that the lengths of the three topological states are d n T, h n T and (1 − d n − h n )T for the on state, off state and the third state respectively.
In this paper the following notations have been adopted.Indexes of a discrete variable z obtained by sampling the continuous time variable z at the time moments nT are written as subscripts.For example z n = z(nT).The values of the continuous time variable at some time moment is denoted using rounded brackets, for example z(nT + d n T) is the value of the continuous variable z at the time moment nT + d n T. The state-space approach is used in this paper for converter description and analysis.As in each topological state a different linear time invariant (LTI) circuit is configured, a dc-dc converter is obviously a piecewise linear system.The set of LTI equations that govern converter operation differ from one topological state to another.They can be written as:

dx A x t B u t t nT n d T dt dx A x t B u t t n d T n d h T dt dx A x t B u t t n d h T n T dt
where n is a positive integer, x(t) is the state vector, u(t) is the input vector containing the voltage sources and Ai, Bi, 1, 3 i = , are the state matrices corresponding to the three topological states.In these conditions it is known [1,3] that converter operation can be described by the following discrete state equation: where the transition matrices involved in (2) are: ; ; ;

d T A d T A d T A d h T A d h T A h T A d T T A d h T A T A
When matrices Ai are nonsingular, then 1 ψ , 2 ψ , 3 ψ can be written in a simplified form as: ) The set of LTI equations that govern converter operation differ from one topological state to another.They can be written as: where n is a positive integer, x(t) is the state vector, u(t) is the input vector containing the voltage sources and A i , B i , i = 1, 3, are the state matrices corresponding to the three topological states.In these conditions it is known [1,3] that converter operation can be described by the following discrete state equation: where the transition matrices involved in (2) are: When matrices A i are nonsingular, then Ψ 1 , Ψ 2 , Ψ 3 can be written in a simplified form as: where I is the unity matrix of converter order.Even when singular, matrices Ψ 1 , Ψ 2 , Ψ 3 can be easily evaluated from (3) using symbolic calculation.It is also important to remark that if B i is a zero matrix, Energies 2018, 11, 663 4 of 25 then the corresponding Ψ i is also zero.This is the situation in the boost converter where B 3 = 0 0 t .More general, in any converter with the transistor connected in series with the supply voltage, matrices B i are zero in all topological states when the transistor is off.
It is known that because of its periodic nature, the steady state solution of the discrete model is constant and therefore the duty cycles d n and h n are constant, denoted by D and H and the input vector u n is also constant, U. Consequently, the steady-state transition matrices in (3) will be denoted by capitals as below: In these conditions the steady-state constant state vector X can be easily found from (2) imposing x n + 1 = x n = X, from which one obtains:

The New Method for Determining Bifurcation Parameter
Because the control of a dc-dc converter is of discrete nature, stability analysis can be performed if an iterative map of the form: is available.Here p bif is a generic bifurcation parameter.The proposed approach is based on the small-signal assumption, requiring the determination of the Jacobian matrix J = ∂ f ∂x n from (7) in a certain steady-state operating point (OP) {u n = U, d n = D, x n = X} and the calculation of the eigenvalues (characteristic multipliers) from the well-known equation: However, in the most general case, the discrete state equation of a dc-dc converter is not of the form (7), but it can be written as: Clearly the form of Equation ( 9) differs from Equation ( 7) because additional variables d n and h n appear.Comparing (2) to ( 9) it easily results that function g in ( 9) is: As previously mentioned, function g depends not only on x n and u n , but also on d n and h n .This is due to the fact that matrices Φ i and Ψ i , i = {1,2,3}, depend on d n and h n , as (3) reveals.It is clear that (10) is not an iterative map of the form (7), as the additional discrete variables d n and h n appear.However, the Jacobian can be still evaluated as will be described below, eliminating parameters d n and h n in the small-signal model, as they depend on x n and u n .In order to achieve this let us first express the state vector at the end of the first and second topological states in terms of the state vector at the beginning of the corresponding topological state.Because in each of these states the circuit is LTI, this can be easily performed and one obtains: x(nT In order to eliminate d n and h n from ( 9), two additional relationships are needed, both relating d n , x n , u n , and possibly p bif .In the most general case, these equations are of the form In any dc-dc converter Equation ( 13) is the control equation, (e.g., proportional voltage control, current mode control, charge control, etc.).Generally, function F 1 is the result of a comparison between converter variables and reference magnitudes and F 1 is not very complicated for most of control techniques.In this paper proportional voltage control is adopted as an example, but the reader can easily apply the principle to any type of control.In proportional voltage-mode control, illustrated in Figure 2 for a boost converter, the difference between the reference voltage V ref and the output voltage v c is amplified and then compared to a rising saw tooth signal v saw .When amplified difference signal is higher than the saw tooth the transistor is switched on, otherwise it will be off.This is the so called natural sampling PWM control, but also uniform sampling can be used, that is when the output voltage is sampled at the beginning of each period and the value of the sample is provided to the control circuit.From the schematic in Figure 2 it results that transistor off switching occurs when: where V U , V L are the peak and the valley value of the sawtooth signal and k is the gain.From (15) the discrete duty cycle is: where the row vector V exv extracts the capacitor voltage from state vector.As function F 1 requires x n as argument, x(nT + d n T) has to be expressed in terms of x n .This is achieved substituting x(nT + d n T) from ( 11) in (16), resulting in: From ( 17) function F 1 is easily identified as: Function F 2 that appears in (14) typically describes the operation mode of the converter.In DCM it states that inductor current equals zero at the end of the second topological state.Mathematically this can be written as: Energies 2018, 11, 663 6 of 25 where the vector V exi is a row vector that extracts the inductor current from the state vector.Obviously, from ( 19) and ( 14) it follows that and intermediate form for of F 2 is given by: The final form of F 2 is found substituting x((n + d n + h n )T) from successively substituting (11) in (12) and the result in (20), obtaining that: Concluding, in order to obtain an equation of the form (7) we have to eliminate d n and h n from the system of equations: where g, F 1 and F 2 are given by ( 10), ( 18) and ( 21 , , , Concluding, in order to obtain an equation of the form (7) we have to eliminate dn and hn from the system of equations: ) , , , , , , , , 0 , , , 0 where g, F1 and F2 are given by ( 10), ( 18) and ( 21 In most of the cases it is impossible to provide an explicit expression of type (7) for xn + 1 from (22).But as long as stability is of concern, local stability can be still exactly predicted.First it has to be stressed out that pbif is a parameter that influences the operating point and for local stability analysis it is a constant value.In the following, the symbol ^ will be used to denote the small-signal perturbations.Linearizing the three equations in (22) around the constant steady-state operating point defined by {U, D, H, X} and taking into account that 0 bif p ∧ = because pbif is constant, one obtains: where in general the symbol  In most of the cases it is impossible to provide an explicit expression of type ( 7) for x n + 1 from (22).But as long as stability is of concern, local stability can be still exactly predicted.First it has to be stressed out that p bif is a parameter that influences the operating point and for local stability analysis it is a constant value.In the following, the symbol ˆwill be used to denote the small-signal perturbations.Linearizing the three equations in (22) around the constant steady-state operating point defined by {U, D, H, X} and taking into account that ∧ p bi f = 0 because p bif is constant, one obtains: From (24) the Jacobian J is easily identified as: Equation ( 25) is the final result that will be used in the analysis, where g F 1 and F 2 are defined by ( 9), ( 13) and ( 14) respectively in the general cases.For DCM operation the expressions of these functions are given by ( 10), ( 18) and ( 21) respectively.Examining (25) it can be seen that the partial derivatives of the functions g, F 1 and F 2 with respect to x n , d n and h n need be calculated.Also the exact steady-state operation point has to be calculated because the linearization process requires these derivatives to be evaluated in the steady-state operating point.This will be presented and performed in the next paragraph.

Partial Derivatives of Functions g, F 1 , F 2 and Exact Steady-State Operating Point Calculation
For deriving the partial derivatives of function g evaluated in the operating point, Equation (10) is used.After some simple algebra the final result is: The partial derivatives of F 1 evaluated in the operating point are straight from (18): Partial derivatives of F 2 in the operating point are directly found from (21).The result is: Energies 2018, 11, 663 In order to complete the Jacobian determination, the exact operating point needs to be calculated.This can be achieved using (6) and taking into account that in steady state the discrete model becomes a dc one, that is x n + 1 = x n = X, d n = D, h n = H and the steady-state transition matrices are those given by (5).It is obvious that X depends on D and H.The relative length D is found from (18) written for steady state, with d n = D. Thus a transcendent equation for D results in the form: After solving this equation, with D known, H is obtained in a similar way from (21) in steady state.The equation for finding H is: In conclusion, first D is found from (35), then with D known H is found from (36) and finally X is given by (6).In other words, Equations ( 35), ( 36) and ( 6) constitute a nonlinear system that has to be numerically solved in order to obtain the steady state solution {D, H, X}, with U known.

Bifurcation Analysis Using the Proposed Mathematical Model
As an example how the theory developed can be applied, a Matlab™ bifurcation analysis is performed for a boost dc-dc converter employing proportional voltage-mode control, similar to that in Figure 2. As Figure 2 shows, parasitics for C and L are ignored although they can be any time included, as their presence only modify the state matrices without any influence on the algorithm.
The authors also performed simulations and experiments on high order converters, namely on the fourth order Sepic, Zeta and excellent agreement between theory, simulations and experiments was found.For the sake of simplicity and for space reasons only the results related to the Boost converter will be presented in the paper.
The bifurcation parameter p bif was chosen to be the amplifier gain k.Converter parameters are: V g = 16 V; V ref = 22 V; C = 220 µF; T = 333.33µs; R = 78 Ω; L = 1209 µH; Transistor on resistance is R on = 0.2 Ω and diode forward voltage is V D = 0.4 V; Saw tooth parameters are V U = 3.5 V; V L = 0.7 V.
For the boost converter under study, choosing the state vector as , the state matrices are: The program flow chart is presented in Figure 3.The Matlab TM code both for the main program BOOST_DCM.m and for the function f_system.mthat is used with the fsolve command for calculating the dc operating point is provided in the Appendix A.
In Table 1 the evolution of the characteristic multipliers values as the bifurcation parameter k is progressively increased is presented.Four decimal places are retained as this is sufficient to provide the accuracy of the results with an at least 0.01%.
The increment for the bifurcation parameter was 1 × 10 6 .It can be remarked that one eigenvalue is always zero and this is explained by the fact that in DCM in any switching cycle the inductor current starts and ends at zero [10].
It can be remarked that at k = 1.1589only one negative characteristic multiplier exits the unit circle.It is known that if a negative real eigenvalue moves out of the unit circle, while other eigenvalues are Energies 2018, 11, 663 9 of 25 within the unit circle, the converter exhibits period-doubling bifurcation [2,10,11].This prognosis will be confirmed by circuit simulation and experiments in the next paragraphs.In Table 1 the evolution of the characteristic multipliers values as the bifurcation parameter k is progressively increased is presented.Four decimal places are retained as this is sufficient to provide the accuracy of the results with an at least 0.01%.
The increment for the bifurcation parameter was 1 × 10 6 .It can be remarked that one eigenvalue is always zero and this is explained by the fact that in DCM in any switching cycle the inductor current starts and ends at zero [10].
It can be remarked that at k = 1.1589only one negative characteristic multiplier exits the unit circle.It is known that if a negative real eigenvalue moves out of the unit circle, while other eigenvalues are within the unit circle, the converter exhibits period-doubling bifurcation [2,10,11].This prognosis will be confirmed by circuit simulation and experiments in the next paragraphs.

Validation through Circuit Simulation
In order to validate the theory developed and implemented in Matlab™, a circuit simulation is performed using the Caspoc™ package [17].For correct results, the variation of the bifurcation parameter has to be very slow compared to converter dynamics, such that the operation to consist of a sequence of quasi steady-states, otherwise significant transients will occur and the bifurcation diagram will be altered.In our case k is varied from 1.1 to 1.2 in 1 s.The simulation schematic is presented in Figure 4.The bifurcation diagram is obtained in SCOPE3 displaying the values of k on the x axis and the sampled output voltage on the y axis.

Validation through Circuit Simulation
In order to validate the theory developed and implemented in Matlab™, a circuit simulation is performed using the Caspoc™ package [17].For correct results, the variation of the bifurcation parameter has to be very slow compared to converter dynamics, such that the operation to consist of a sequence of quasi steady-states, otherwise significant transients will occur and the bifurcation diagram will be altered.In our case k is varied from 1.1 to 1.2 in 1 s.The simulation schematic is presented in Figure 4.The bifurcation diagram is obtained in SCOPE3 displaying the values of k on the x axis and the sampled output voltage on the y axis.
The simulated bifurcation diagram obtained in SCOPE3 is enlarged presented in Figure 5, with a magnified image around the bifurcation point in Figure 6.The theoretical value of 1.159, marked with the red circle on the graph is the bifurcation threshold.In the small table below the Figure 6 the values in the column denoted by y represent voltage of the power supply and the output voltage values respectively, when the bifurcation appears.From Figure 6 it can be seen that the first bifurcation with period doubling occurs at k = 1.159, which is an excellent agreement with the Matlab™ prognosis of 1.1589, with 0.092% relative error.In Figure 7a, a time representation of the inductor current in steady-state operation at a value of k = 1.1, lower than the bifurcation value is presented.In Figure 7b the phase portrait of the system for the same value of k is depicted.From Figure 6 it can be seen that the first bifurcation with period doubling occurs at k = 1.159, which is an excellent agreement with the Matlab™ prognosis of 1.1589, with 0.092% relative error.In Figure 7a, a time representation of the inductor current in steady-state operation at a value of k = 1.1, lower than the bifurcation value is presented.In Figure 7b the phase portrait of the system for the same value of k is depicted.From Figure 6 it can be seen that the first bifurcation with period doubling occurs at k = 1.159, which is an excellent agreement with the Matlab™ prognosis of 1.1589, with 0.092% relative error.In Figure 7a, a time representation of the inductor current in steady-state operation at a value of k = 1.1, lower than the bifurcation value is presented.In Figure 7b the phase portrait of the system for the same value of k is depicted.Both representations in Figure 7a,b confirm that stable operation is achieved.
In Figure 7c, the inductor current waveform for a value of k = 1.2, a value that exceeds the bifurcation threshold, being located in the range where the bifurcation diagram prognoses a period doubling is enfaced.The corresponding phase portrait is shown in Figure 8. Figures 7c and 8, reveal that period doubling occurs.
It can be remarked that still increasing k after the first bifurcation (for example for k = 1.215), period 4 bifurcation occurs.This bifurcation can also be numerically determined calculating the Jacobian of the function f(f), where f is the function in relationship (6).Obviously, the calculations are more complex, but the approach is the same as the one used for determining the period 2 bifurcation.Further increase of the bifurcation parameter causes the system to burst into chaos.Both representations in Figure 7a,b confirm that stable operation is achieved.In Figure 7c, the inductor current waveform for a value of k = 1.2, a value that exceeds the bifurcation threshold, being located in the range where the bifurcation diagram prognoses a period doubling is enfaced.The corresponding phase portrait is shown in Figure 8. Figures 7c and 8, reveal that period doubling occurs.

Experimental Results
The experimental setup schematic is presented in Figure 9.The differential amplifier consisting of IC1 and resistors R1 to R4 impose the value of the gain k.With the typical condition of the differential amplifier fulfilled, that is It can be remarked that still increasing k after the first bifurcation (for example for k = 1.215), period 4 bifurcation occurs.This bifurcation can also be numerically determined calculating the Jacobian of the function f (f ), where f is the function in relationship (6).Obviously, the calculations are more complex, but the approach is the same as the one used for determining the period 2 bifurcation.Further increase of the bifurcation parameter causes the system to burst into chaos.

Experimental Results
The experimental setup schematic is presented in Figure 9.The differential amplifier consisting of IC1 and resistors R 1 to R 4 impose the value of the gain k.With the typical condition of the differential amplifier fulfilled, that is gain k is given by For the values of the resistors used, k range is between 1.08 and 2.16.The saw tooth waveform is the voltage across capacitor C that together with potentiometer P2 imposes the operation frequency of the SG3524 represented by IC0 circuit and the switching frequency of the boost converter.The two outputs of the SG3524 circuit are connected in parallel, providing the short set pulses of the flip-flop.The pulse width is adjusted by potentiometer P1.IC2 is used to amplify the pulses such that to cross the threshold voltage of the flip-flop.The reset signal for the flip-flop is provided by the comparator represented by IC3.The switching frequency is 3 kHz, equal to that used in the simulation.
In Figure 10 system behavior for k = 1.08 is presented.The theory developed predicts stable operation both in Matlab™ and Caspoc™.This is experimentally confirmed by the waveforms in Figure 10, where all signals exhibit the switching period, with the typical DCM shape of the inductor current.The phase portrait is depicted in Figure 11, additionally confirming stable operation.Then the gain factor k is then progressively increased and the inductor current is carefully monitored.When two consecutive periods differ it is decided that the system enters bifurcation and the corresponding value of k is retained.Figure 12 presents the inductor current and output voltage of the converter when the system goes from stable operation to bifurcation.It can be seen that two consecutive periods are not identical as the two inductor current peak values begin to slightly differ.The value of parameter k corresponding to this situation is k = 1.17, while the Matlab™ and Caspoc™ programs both predict that bifurcation appears at k = 1.16.Obviously, the relative error between the simulation and experiment is very low, equal to 0.94%.The error sources come from the fact that nonideal devices parameters, such as R on and V D are not exactly known as only the catalog values were used and the series resistance of the inductor and capacitor were neglected.Other error sources are of course given by the equipment used but these are lower compared to those mentioned above.Thus the experimental result accurately validates the theoretical method proposed.Further increasing the gain will reveal a more evident bifurcation.For example, in Figure 13 inductor current and output voltage are depicted for k = 1.20.The period of the system is clear half of the switching period (period 2 subharmonic) and the phase portrait enfaced in Figure 14 is typical for a period 2 bifurcation system.Further increasing the gain will reveal a more evident bifurcation.For example, in Figure 13 inductor current and output voltage are depicted for k = 1.20.The period of the system is clear half of the switching period (period 2 subharmonic) and the phase portrait enfaced in Figure 14 is typical for a period 2 bifurcation system.
A higher value of the gain will lead the converter into chaos.This situation is presented in Figure 15, both the waveforms and the phase portrait exhibiting the chaotic behavior.A higher value of the gain will lead the converter into chaos.This situation is presented in Figure 15, both the waveforms and the phase portrait exhibiting the chaotic behavior.

Input Voltage as a Bifurcation Parameter
In practice also the input voltage and the load resistance could affect the stability and the bifurcation behavior of a dc-dc.For this reason and in order to additionally prove the validity of the proposed method, a bifurcation analysis was carried out in the situation when the bifurcation   In Figure 18 the behavior for an input voltage of 17.2 V, higher than the bifurcation threshold is depicted.The bifurcation phenomenon is obvious.Further increase of the input voltage will drive the system into chaos, as Figure 19 presents for an input voltage equal to 17.8 V.In Figure 18 the behavior for an input voltage of 17.2 V, higher than the bifurcation threshold is depicted.The bifurcation phenomenon is obvious.Further increase of the input voltage will drive the system into chaos, as Figure 19 presents for an input voltage equal to 17.8 V.    Hence it can be concluded that the proposed method is valid, powerful and with a high degree of generality.Hence it can be concluded that the proposed method is valid, powerful and with a high degree of generality.

Conclusions
A new mathematical method for analyzing and predicting bifurcation behavior of DCM operated dc-dc converters is proposed.Compared to other methods this new approach performs exact calculation, without truncating any state matrix and thus correctly analyzing all the waveforms that in most of the approaches are approximated to be piecewise linear.On the other side, the new method is general, in the sense that it can be applied to any converter topology employing different types of control by simply changing only the matrices that describe the state-space model and one constraint equation.The proposed algorithm can be easily absorbed by any mathematical program that is able to calculate the matrix exponential.The method provides excellent accuracy, better than any other previous reported approach.The experimental results confirmed the perfect agreement between theory and practice.The proposed technique can also be used to exactly predict the parameter values that determine the operation at the border between DCM and CCM and thus to practically provide exact prediction of border collision.The method is different from those in [1,2,9,12,15] and can be also successfully applied both to uniform sampling and to different control techniques [15] employing CCM or other discontinuous modes, such as discontinuous capacitor voltage mode (DCVM) [18,19].

Figure 1 .
Figure 1.Transistor switching function, inductor voltage and inductor current in a DCM operated dc-dc converter.

Figure 1 .
Figure 1.Transistor switching function, inductor voltage and inductor current in a DCM operated dc-dc converter.
the derivative of function z with respect to the variable r evaluated in the operating point.
the symbol ∂z ∂r OP denotes the derivative of function z with respect to the variable r evaluated in the operating point.As the model described by (23) is LTI, from its last two equations ∧ d n and ∧ h n can be expressed in terms of ∧ x n and ∧ u n , leading to the proposed mathematical model:

Energies 2018 ,
11, x FOR PEER REVIEW 9 of 26 Initialize pbif =k and its increment incr Define the state vectors and matrices Determine the steay-state solution {U, D, H, X} using fsolve command Calculate the Jacobian and the characteristic multipliers l pbif = pbif + incr 1 : pbif , D, H, X, l S T O P Declare the symbolic variables Initialize converter circuit and control Declare the symbolic variables Initialize converter circuit and control Evaluate the partial derivatives:

Figure 3 .
Figure 3. Flow chart of the Matlab™ program for determining the critical bifurcation parameter value.

Figure 3 .
Figure 3. Flow chart of the Matlab™ program for determining the critical bifurcation parameter value.

Figure 4 .
Figure 4. Caspoc™ schematic for the proportional voltage-mode controlled boost converter.Figure 4. Caspoc™ schematic for the proportional voltage-mode controlled boost converter.

Figure 6 .
Figure 6.Magnified image around the bifurcation point.Clearly bifurcation appears at k = 1.159.The table below provides the instantaneous values of the simulated waveforms, together with their minimum and maximum values.

Figure 5 .
Figure 5.The simulated bifurcation diagram of the boost DCM converter employing proportional voltage-mode control.

Energies 2018 , 26 Figure 5 .
Figure 5.The simulated bifurcation diagram of the boost DCM converter employing proportional voltage-mode control.

Figure 6 .
Figure 6.Magnified image around the bifurcation point.Clearly bifurcation appears at k = 1.159.The table below provides the instantaneous values of the simulated waveforms, together with their minimum and maximum values.

Figure 6 .
Figure 6.Magnified image around the bifurcation point.Clearly bifurcation appears at k = 1.159.The table below provides the instantaneous values of the simulated waveforms, together with their minimum and maximum values.

Figure 7 .
Figure 7. (a) The inductor current waveform for k = 1.1.The horizontal dotted line denotes the zero level and the solid line the inductor current.(b) Phase portrait for stable operation, k = 1.1.(c) Inductor current for parameter k = 1.2, higher than bifurcation threshold value of 1.1589.Unstable operation with period doubling (period 2 subharmonic) is obvious.

Figure 9 .
Figure 9.Schematic of the experimental boost DCM converter employing proportional voltage-mode control.

Figure 9 .
Figure 9. Schematic of the experimental boost DCM converter employing proportional voltage-mode control.

Figure 10 .
Figure 10.Set signal and reset signals of IC4, inductor current and output voltage waveform (this up to down order) for k = 1.08.Stable operation can be remarked.

Figure 10 . 26 Figure 11 .
Figure 10.Set signal and reset signals of IC4, inductor current and output voltage waveform (this up to down order) for k = 1.08.Stable operation can be remarked.Energies 2018, 11, x FOR PEER REVIEW 16 of 26

Figure 12 .
Figure 12.Inductor current and output voltage waveforms for k = 1.17.Notice that the peak inductor current values start to slightly differ in two consecutive periods and the period of the inductor current doubles compared to stable operation.

Figure 12 .
Figure 12.Inductor current and output voltage waveforms for k = 1.17.Notice that the peak inductor current values start to slightly differ in two consecutive periods and the period of the inductor current doubles compared to stable operation.

Figure 13 .
Figure 13.Inductor current and output voltage waveforms for k = 1.20.Bifurcation with period 2 operation is evident.

Figure 13 .
Figure 13.Inductor current and output voltage waveforms for k = 1.20.Bifurcation with period 2 operation is evident.

Figure 13 .
Figure 13.Inductor current and output voltage waveforms for k = 1.20.Bifurcation with period 2 operation is evident.

Figure 15 .
Figure 15.Chaotic operation of the converter corresponding to k = 1.50.

Figure 19 .
Figure 19.Chaotic operation of the converter corresponding to Vg = 17.80 V.

Figure 19 .
Figure 19.Chaotic operation of the converter corresponding to V g = 17.80 V.

Table 1 .
Characteristic multipliers of the boost converter for increasing values of the gain k.The bold value in the first column denotes the threshold when bifurcation occurs.

Table 1 .
Characteristic multipliers of the boost converter for increasing values of the gain k.The bold value in the first column denotes the threshold when bifurcation occurs.