Quantifying the Nonlinear Dynamic Behavior of the DC-DC Converter via Permutation Entropy

Quantifying nonlinear dynamic behaviors, such as bifurcation and chaos, in nonlinear systems are currently being investigated. In this paper, permutation entropy is used to characterize these complex phenomena in nonlinear direct current-direct current (DC-DC) converter systems. A mode switching time sequence (MSTS), containing the information from different periodic states, is obtained in a DC-DC converter by reading the inductor current when altering the switching mode. To obtain the nonlinear characteristics of this system, the concept of permutation entropy of symbolic probability distribution properties is introduced and the structure of the chaotic system is reproduced based on the theory of phase space reconstruction. A variety of nonlinear dynamic features of the DC-DC converter are analyzed using the MSTS and permutation entropy. Finally, a current-mode-controlled buck converter is reviewed as a case to study the quantification of nonlinear phenomena using permutation entropy as one of the system parameters changes.


Background
Direct current-direct current (DC-DC) switching converters are a typical nonlinear system that enable the observation of nonlinear behaviors, such as period-doubling bifurcation, border collision bifurcation, frequency-locking phenomenon, quasi-period, and chaos [1][2][3][4][5].The study of nonlinear behaviors is helpful for improving the performance of DC-DC converters.

Formulation of the Problem of Interest for This Investigation
Generally speaking, these nonlinear phenomena can be observed by phase plane, bifurcation diagrams, and Poincaré sections in DC-DC converters [6][7][8], but these only provide qualitative analyses.To examine nonlinear phenomena in a timely manner and identify them accurately, quantitative analysis is necessary.At present, some quantitative analysis methods have been proposed to process these nonlinear features, and the largest Lyapunov exponent (LLE) is the most common [9,10].

Literature Survey
In practical applications, analyzing the diverse nonlinear phenomena in LLE is challenging due to the complex of the phase space trajectory and its solution process is complicated [7,11].Researchers [12] found the symbolic time sequence contains a large amount of DC-DC converter system state information, including nonlinear stability, bifurcation information, and degree of confusion.Both symbolic sequence analysis and entropy are utilized to quantitatively analyze DC-DC converters.
Energies 2018, 11, 2747 3 of 15 Set X = [X 1 X 2 X 3 • • • X p ] T to the state vector (with p being the numbers of variables and superscript T denoting transposition) of a DC-DC converter.So, in the ith control period, the discretized expression of X can be expressed as: where S k represents the relationship between the nitial and final terative alues in the circuit topology as the operating mode k is presented (with its expression being only related to the topology and the selected state variable).t i and t ik represent the start time of the ith period and the mode k occurs in the period, and X ik = X(t i + t i1 + t i2 + • • • + t ik ).
Suppose the system has j i (1 ≤ j i ≤ 2 m -1, where m indicates the number of switching devices) modes in the ith control period.As the system stabilizes, we take x L in the state variable X as the object of observation and sample it when the mode changes from each control period.Hence, we obtain the MSTS, i.e., x L−MSTS = {x 1 , x 11 , x 12 , . . ., x 1j 1 , x 21 , . . ., x (N−1)j (N−1) , x N1 , x N2 , . . ., x Nj N }, Considering the equation x(t i+1 ) = x(t i + t i1 + . . .+ t ij i ), the x i+1 has to be left out from Equation (2). Figure 1 shows an example of two switching devices existing in the converter system, which means it can express up to three modes, where T is the control period and t i is the time of the ith switching moment (t i −t i-1 = T).
where Sk represents the relationship between the nitial and final terative alues in the circuit topology as the operating mode k is presented (with its expression being only related to the topology and the selected state variable).ti and tik represent the start time of the ith period and the mode k occurs in the period, and 12 () Suppose the system has ji (1 ≤ ji ≤ 2 m -1, where m indicates the number of switching devices) modes in the ith control period.As the system stabilizes, we take xL in the state variable X as the object of observation and sample it when the mode changes from each control period.Hence, we obtain the MSTS, i.e., 1 -MSTS 1 11 , the xi+1 has to be left out from Equation (2).
Figure 1 shows an example of two switching devices existing in the converter system, which means it can express up to three modes, where T is the control period and ti is the time of the ith switching moment (ti−ti-1 = T).

Phase Space Reconstruction Theory
The low-dimensional coordinate system cannot reflect the complex nonlinear states, as demonstrated by Packard et al. [25].Packard et al. [25] proposed two methods of reconstructing phase space in a time sequence.The coordinate delay reconstruction method is widely used.In this method, the one-dimensional time sequence is constructed by embedding the dimension m and the delay time τ to form a new m-dimensional vector sequence.The phase space of the discrete time sequence x(n) of length N can be obtained: Afterward, the Takens embedding theorems [26] proves the rationality of phase space reconstruction and proposes that the embedding dimension of phase space should satisfy m ≥ 2D + 1, where D is the singular attractor dimension of the system and it is a geometric invariant.A dynamic system with correlation dimension D can be described through 2D + 1 independent variables at most [7].This means that the geometry of the dynamic system is completely open when m is greater than

Phase Space Reconstruction Theory
The low-dimensional coordinate system cannot reflect the complex nonlinear states, as demonstrated by Packard et al. [25].Packard et al. [25] proposed two methods of reconstructing phase space in a time sequence.The coordinate delay reconstruction method is widely used.In this method, the one-dimensional time sequence is constructed by embedding the dimension m and the delay time τ to form a new m-dimensional vector sequence.The phase space of the discrete time sequence x(n) of length N can be obtained: Afterward, the Takens embedding theorems [26] proves the rationality of phase space reconstruction and proposes that the embedding dimension of phase space should satisfy m ≥ 2D + 1, where D is the singular attractor dimension of the system and it is a geometric invariant.A dynamic system with correlation dimension D can be described through 2D + 1 independent variables at most [7].This means that the geometry of the dynamic system is completely open when m is greater than the minimum embedding dimension 2D + 1.At this time, the geometric invariant properties of the system, such as the Lyapunov exponent, are independent of the embedding dimension m.In addition, the delay time can be obtained by the autocorrelation method.The time sequence of a state variable can reconstruct the phase space of the system, for it contains information about all variables participating in the dynamic system.According to the Takens embedding theorem, the m-dimensional state space, expressed in Equation ( 4), has the same meaning as the original state space topologically; that is, the original state variable information can be reproduced.

Symbol Permutation Method
Entropy can quantify the periodic states and chaos by time sequence.Its concept is sourced in the second law of thermodynamics.After Shannon introduced it into information theory, Shannon entropy became an important quantitative means to measure the probability distribution of a sequence signal [22].
For any discrete sequence {x(i), i = 1, 2, . . ., n}, the corresponding probability distribution {p(i)|i = 1, 2, . . ., n } can be calculated, where p(i) is equal to the proportion of x(i) in the entire sequence.The Shannon entropy H P can be defined as H P = −∑ p(i) ln p(i).The larger the H P , the more uneven the signal distribution.This indicates that the system is highly complex, and the time sequence is difficult to predict.On the contrary, the information contained in the signal is easy to find for small H P , and the prediction can be realized with minimal error.Particularly, in the case where H P = 0, the sequence signal could be theoretically predictable, such as the step signal at t > 0.
In order to accurately quantify the nonlinear features, we first need to acquire certain statistical characteristics of the signal.In 2001, Bant et al. proposed the concept of permutation entropy based on phase space reconstruction theory, and successfully applied it to the complex quantization process of time sequences [23].Equation (5) was generated by sorting Equation (4) in ascending order: We determined the corresponding symbol permutation: For m-dimensional phase space reconstruction, we studied all n! permutations s(g) of (1, 2, . . . ,m − 1) where Then, comparing K i logically with the sequence of symbols s(g), we define p(g) as the probability that s(g) appears in K 1 ; K 2 ; . . .; K N−(M−1)τ : Energies 2018, 11, 2747 where Card is a function that represents the number of elements in a collection According to the definition of the Shannon entropy, the permutation entropy can be defined as: The whole -process -can be illustrated by the special example in Figure 2. In this graph, a sequence of period-6 was chosen as the object of analysis below.We set the embedding dimension m = 3 and the delay time τ = 2.It was assumed that the sequence is sampled from the first 13 data points onto the curve shown in Figure 2. In this case, the state space is divided into 3! (=6) parts, so there are 6 mutually different arrangement symbols.According to the definition in Equation (3), the first sample point has a three-dimensional (3D) coordinate of (1.3660,-1.0000,−0.3660).Following Equation ( 5), we found that x(i + τ) ≤ x(i + 2τ) ≤ x(i), i = 1.Then, the corresponding permutation is (−1.0000,−0.3660,1.3660)and the corresponding symbol sequence {120} is obtained by referring to Equation (6).Therefore, combined with Equation (3), the sequence (1.3660,−0.6340,−1.0000,−2.3660,−0.3660,3.0000,1.3660,−0.6340,−1.0000,−2.3660,−0.3660,3.0000,1.3660)can be converted into 9 symbol sequences for representation ({120}, {102}, {012}, {021}, {201}, {210}, {120}, {102}, {012}).
Energies 2018, 11, x FOR PEER REVIEW 5 of 15 where Card is a function that represents the number of elements in a collection According to the definition of the Shannon entropy, the permutation entropy can be defined as: The whole -process -can be illustrated by the special example in Figure 2. In this graph, a sequence of period-6 was chosen as the object of analysis below.We set the embedding dimension m = 3 and the delay time τ = 2.It was assumed that the sequence is sampled from the first 13 data points onto the curve shown in Figure 2. In this case, the state space is divided into 3! (=6) parts, so there are 6 mutually different arrangement symbols.According to the definition in Equation ( 3), the first sample point has a three-dimensional (3D) coordinate of (1.3660,-1.0000,−0.3660).Following Equation ( 5), we found that (  As analyzed above, the previous six symbol sequences are different from each other.Starting from the seventh symbol sequence, the periodic extension is performed in the order in which the first six sequences appeared.This means that all possible symbol sequences are included.The remaining other parameters are unchanged, the analysis is performed again with m = 4, and the above calculation is repeated to obtain 7 symbol sequences ({1230},{1032},{3012},{3021},{2031}, {2130},{1230}).
We considered . The permutation entropy calculated by Equation ( 8) is approximately 1.7918, provided the sampling sequence is long enough.The comparison shows that the distribution of symbol sequences in 3 m  is consistent with 4 m  .That means if the cycle count is less than m!, the m would not affect the amount of information acquired.So, the information that can be obtained is also uniform.
The permutation entropy is a method proposed according to the spatial characteristics of time.We can obtain the new state space Y by reconstructing the phase space.After sorting each mdimensional vector in Y, ( 1) Nm  sorts symbol sequence K can be obtained.For the embedding As analyzed above, the previous six symbol sequences are different from each other.Starting from the seventh symbol sequence, the periodic extension is performed in the order in which the first six sequences appeared.This means that all possible symbol sequences are included.The remaining other parameters are unchanged, the analysis is performed again with m = 4, and the above calculation is repeated to obtain 7 symbol sequences ({1230},{1032},{3012},{3021},{2031},{2130},{1230}).We considered ).The permutation entropy calculated by Equation ( 8) is approximately 1.7918, provided the sampling sequence is long enough.The comparison shows that the distribution of symbol sequences in m = 3 is consistent with m = 4.That means if the cycle count is less than m!, the m would not affect the amount of information acquired.So, the information that can be obtained is also uniform.
The permutation entropy is a method proposed according to the spatial characteristics of time.We can obtain the new state space Y by reconstructing the phase space.After sorting each m-dimensional vector in Y, N − (m − 1)τ sorts symbol sequence K can be obtained.For the embedding dimension m, the total number of corresponding order permutations is m!.To ensure the degree of uniformity of quantization, N − (m − 1)τ >> m! is generally required, and N >> m! is adopted in practical application.Furthermore, if m is small, the signal with a higher multiple periodic cannot be identified.Hence, the information volume read is correspondingly reduced, eventually resulting in a loss of information.However, if the selection of m is oversized, it often leads to N < m!.As a consequence, the identification accuracy decreases and the error trend increases.
The noise-free chaotic sequence studied by Taken's theorem of embedding indicates that there is no specific limitation on the delay time τ.τ does not take the optimal delay time only affects the Euclidean geometry of the reconstructed attractor when reconfiguring the phase space.So, it affects the calculation of the correlation dimension.However, it will not affect the reconstructive attractor's unambiguous representation of the system dynamics.In addition, a completely independent situation of the two adjacent coordinate components could occur if the value of τ is too large, in which case, the projection of the chaotic attractor's trajectory in the same direction would become irrelevant.Hence, smaller values of τ should be chosen according to the specific situation; Bandt and Pompe suggested choosing a τ value of 1 [20].Therefore, the embedding dimension m is considered to be a key parameter of spatial reconstruction, and it is commonly recommended to select m within the range of 3 to 7 for systems with less than three dimensions.
The traditional symbol sequence divides the time sequence into n regions by amplitude coarse granulation to obtain the symbol sequence.After the window with length L is selected, the block entropy (BE) of the symbol sequence is obtained and its maximum value is not more than ln(2 L ).For the permutation entropy, after the m-dimensional phase space reconstruction, the maximum value of the permutation entropy does not exceed ln(m!).That is to say, m! >> 2 L as m = L > 3. Compared with the traditional symbol sequence method, the permutation entropy can read more information and does not need to perform amplitude layering.The method of combination of the permutation entropy and MSTS is elaborated in Section 4. Combined with MSTS, we compare permutation entropy and sample entropy in Section 5. Given the logical calculation used in the method, the permutation entropy has the advantages of high robustness and fast calculation speed.

Combination of MSTS and Permutation Entropy
In common DC-DC switching converters that include both an inductor and a capacitor, the trend in the current variable or voltage variable can reflect the operating mode of the system.For a switched capacitor DC-DC converter [27], since it does not contain an inductor, the capacitor voltage can be used as an object for extracting the MSTS.Peak current single loop control is characterized by the inability to stabilize the output voltage and produces a rich nonlinear behavior.We used the current variable as the observation signal for the common converter (Figure 1).The procedure of K as m = 3, τ = 1 can be observed in Figure 3. !Nm  is adopted in practical application.Furthermore, if m is small, the signal with a higher multiple periodic cannot be identified.Hence, the information volume read is correspondingly reduced, eventually resulting in a loss of information.However, if the selection of m is oversized, it often leads to !Nm  .As a consequence, the identification accuracy decreases and the error trend increases.
The noise-free chaotic sequence studied by Taken's theorem of embedding indicates that there is no specific limitation on the delay time τ.τ does not take the optimal delay time only affects the Euclidean geometry of the reconstructed attractor when reconfiguring the phase space.So, it affects the calculation of the correlation dimension.However, it will not affect the reconstructive attractor's unambiguous representation of the system dynamics.In addition, a completely independent situation of the two adjacent coordinate components could occur if the value of τ is too large, in which case, the projection of the chaotic attractor's trajectory in the same direction would become irrelevant.Hence, smaller values of τ should be chosen according to the specific situation; Bandt and Pompe suggested choosing a τ value of 1 [20].Therefore, the embedding dimension m is considered to be a key parameter of spatial reconstruction, and it is commonly recommended to select m within the range of 3 to 7 for systems with less than three dimensions.
The traditional symbol sequence divides the time sequence into n regions by amplitude coarse granulation to obtain the symbol sequence.After the window with length L is selected, the block entropy (BE) of the symbol sequence is obtained and its maximum value is not more than ln(2 ) L .For the permutation entropy, after the m-dimensional phase space reconstruction, the maximum value of the permutation entropy does not exceed ln( !) m .That is to say, !2 L m  as 3 mL .Compared with the traditional symbol sequence method, the permutation entropy can read more information and does not need to perform amplitude layering.The method of combination of the permutation entropy and MSTS is elaborated in Section 4. Combined with MSTS, we compare permutation entropy and sample entropy in Section 5. Given the logical calculation used in the method, the permutation entropy has the advantages of high robustness and fast calculation speed.

Combination of MSTS and Permutation Entropy
In common DC-DC switching converters that include both an inductor and a capacitor, the trend in the current variable or voltage variable can reflect the operating mode of the system.For a switched capacitor DC-DC converter [27], since it does not contain an inductor, the capacitor voltage can be used as an object for extracting the MSTS.Peak current single loop control is characterized by the inability to stabilize the output voltage and produces a rich nonlinear behavior.We used the current variable as the observation signal for the common converter (Figure 1).The procedure of K as m = 3, τ = 1 can be observed in Figure 3. Combined with the above analysis, the permutation entropy has the following characteristics when applied to a DC-DC converter: Combined with the above analysis, the permutation entropy has the following characteristics when applied to a DC-DC converter: (1) Assuming the signal is random, such as white noise, permutation entropy is expected to be close to the maximum value ln(m!).Considering that there are at least two operation modes in period-1, the recurring number of the MSTS is greater than or equal to 2. Hence, the minimum value of PE is ln(2) = 0.6931.(2) When the system does not enter the chaotic state by changing m, the acquired permutation entropy remains uniform.This means that increasing or decreasing m as appropriate does not affect the value of the permutation entropy.(3) DCM and CCM can be further divided into different working patterns according to the different switching numbers of modes.Border collision occurs when the working pattern combination changes.The complexity of the systems in DCM working patterns is usually greater than that of the CCM working pattern.Therefore, the more DCM patterns the system contains, the higher the value of the permutation entropy.After entering chaos, the permutation sequence K shows a certain disorder.Because the chaotic motion is pseudo-random, the permutation entropy obtained in the chaotic state is not able to reach the maximum value.

Application Example
In this section, the peak current control buck DC-DC converter is analyzed using the above method.By quantifying its complexity, the behaviors of the period-doubling bifurcation, border collision bifurcation, and chaos present in the converter are quantitatively described.The circuit model is shown in Figure 4. Relevant parameters were selected as follows: E = 20 V, R = 19 Ω, L = 3.3 mH, C = 1000 uF, T = 400 us, and the variation range of the control variable I ref was controlled within 0.12-1.32.As shown in Figure 4, there are two switching devices in the circuit: the fully controlled device V and the uncontrollable diode VD.Combined with the principle of the buck converter, we obtained three working modes of the circuit, as shown in Figures 5 and 6.These figures depict the switching mode (V,VD) = (on,off )→(V,VD) = (off,on)→(V,VD) = (off,off ) of each switching period.The inductor current is x L > 0 during mode-1 and mode-2 and x L = 0 during mode-3.In Figure 6, Ib 1 and Ib 2 are the borderlines in the buck converter with peak current control.Analyzing the characteristics of the buck converter and using the precise discretization method to model the system, we determined the work pattern of each control cycle.After multiple iterations, the inductor current x L was numerically sampled under the overall steady state, and MSTS was obtained.Here, m = 4 and τ = 1 were selected.Figure 7 shows the permutation entropy for the MSTS as a function of the reference current for different nonlinear feature.For the data of the chaotic state part, since the irregularity was relatively strong, cycle extension was performed to homogenize the order.
Energies 2018, 11, x FOR PEER REVIEW 7 of 15 (1) Assuming the signal is random, such as white noise, permutation entropy is expected to be close to the maximum value ln( !) m .Considering that there are at least two operation modes in period-1, the recurring number of the MSTS is greater than or equal to 2. Hence, the minimum value of PE is ln(2) 0.6931  .
(2) When the system does not enter the chaotic state by changing m,the acquired permutation entropy remains uniform.This means that increasing or decreasing m as appropriate does not affect the value of the permutation entropy.
(3) DCM and CCM can be further divided into different working patterns according to the different switching numbers of modes.Border collision occurs when the working pattern combination changes.The complexity of the systems in DCM working patterns is usually greater than that of the CCM working pattern.Therefore, the more DCM patterns the system contains, the higher the value of the permutation entropy.After entering chaos, the permutation sequence K shows a certain disorder.Because the chaotic motion is pseudo-random, the permutation entropy obtained in the chaotic state is not able to reach the maximum value.

Application Example
In this section, the peak current control buck DC-DC converter is analyzed using the above method.By quantifying its complexity, the behaviors of the period-doubling bifurcation, border collision bifurcation, and chaos present in the converter are quantitatively described.The circuit model is shown in Figure 4. Relevant parameters were selected as follows: E = 20 V, R = 19 Ω, L = 3.3 mH, C = 1000 uF, T = 400 us, and the variation range of the control variable Iref was controlled within 0.12-1.32.As shown in Figure 4, there are two switching devices in the circuit: the fully controlled device V and the uncontrollable diode VD.Combined with the principle of the buck converter, we obtained three working modes of the circuit, as shown in Figures 5 and 6.These figures depict the switching mode (V,VD) = (on,off)→(V,VD) = (off,on)→(V,VD) = (off,off) of each switching period.The inductor current is xL ＞ 0 during mode-1 and mode-2 and xL = 0 during mode-3.In Figure 6, Ib1 and Ib2 are the borderlines in the buck converter with peak current control.Analyzing the characteristics of the buck converter and using the precise discretization method to model the system, we determined the work pattern of each control cycle.After multiple iterations, the inductor current xL was numerically sampled under the overall steady state, and MSTS was obtained.Here, m = 4 and τ = 1 were selected.Figure 7 shows the permutation entropy for the MSTS as a function of the reference current for different nonlinear feature.For the data of the chaotic state part, since the irregularity was relatively strong, cycle extension was performed to homogenize the order.As shown in Figure 8, the inductor current xL is sampled at the operating mode interval and the bifurcation diagram of the system with Iref as a changing parameter can be obtained.
By comparing Figures 7 and 8, the dynamic behavior of the buck converter in the peak current control can be clearly analyzed.Similarly, each state of the border collision bifurcation and the standard bifurcation are represented in Table 1.In this table, S and MSTS-PE indicate the permutation symbol of the MSTS and the corresponding permutation entropy, respectively.As shown in Figure 8, the inductor current xL is sampled at the operating mode interval and the bifurcation diagram of the system with Iref as a changing parameter can be obtained.
By comparing Figures 7 and 8, the dynamic behavior of the buck converter in the peak current control can be clearly analyzed.Similarly, each state of the border collision bifurcation and the standard bifurcation are represented in Table 1.In this table, S and MSTS-PE indicate the permutation symbol of the MSTS and the corresponding permutation entropy, respectively.As shown in Figure 8, the inductor current x L is sampled at the operating mode interval and the bifurcation diagram of the system with I ref as a changing parameter can be obtained.

Iref/A S MSTS-PE State
By comparing Figures 7 and 8, the dynamic behavior of the buck converter in the peak current control can be clearly analyzed.Similarly, each state of the border collision bifurcation and the standard bifurcation are represented in Table 1.In this table, S and MSTS-PE indicate the permutation symbol of the MSTS and the corresponding permutation entropy, respectively.1.As Iref ∈ [0.1200,0.2778], the system works in the F3 pattern of period-1.In this pattern, there are three modes in each period of the system, which means the system is in the DCM state.Then, the first border collision bifurcation occurs in the system at Iref = 0.2779 A. Because the inductor current in the steady state collides with Ib2, the operating mode changes and the working patterns of the converter system change from F3 to F2. 2. During Iref ∈[0.2779,0.8295], the system still works in the period-1 state, but the operating mode is F2.At this time, the system has two modes in each period.When Iref = 0.8296, the system generates period doubling bifurcation and enters period-2.Since there is no border collision bifurcation, the working patterns are extended to F2F2. 3. When Iref ∈[0.8296,1.1946], the system works in period-2.We found that two bounder collisions of the system occur in this area because the inductance current of the system collides with the borderline of Ib2 and Ib1 at Iref = 0.9007 and Iref = 1.1578A, respectively.Hence, the operating pattern of the converter system changes from F2F2 at Iref ∈ [0.8296,0.9006] to F2F3 at Iref ∈ [0.9007,1.1577]and to F1F3 at Iref ∈[1.1578,1.1946].From the entropy value of three period-2 orbits shown in Figure 7, the permutation entropy value of period-2 is greater than the corresponding entropy value of any working pattern in the entirety of period-1.The sample entropy, used for assessing the complexity of the time sequence shown in Figure 9, is the negative logarithm of the probability of if two sets of simultaneous data points of length m have a distance smaller  1.As Iref ∈ [0.1200,0.2778], the system works in the F3 pattern of period-1.In this pattern, there are three modes in each period of the system, which means the system is in the DCM state.Then, the first border collision bifurcation occurs in the system at Iref = 0.2779 A. Because the inductor current in the steady state collides with Ib2, the operating mode changes and the working patterns of the converter system change from F3 to F2. 2. During Iref ∈[0.2779,0.8295], the system still works in the period-1 state, but the operating mode is F2.At this time, the system has two modes in each period.When Iref = 0.8296, the system generates period doubling bifurcation and enters period-2.Since there is no border collision bifurcation, the working patterns are extended to F2F2. 3. When Iref ∈[0.8296,1.1946], the system works in period-2.We found that two bounder collisions of the system occur in this area because the inductance current of the system collides with the borderline of Ib2 and Ib1 at Iref = 0.9007 and Iref = 1.1578A, respectively.Hence, the operating pattern of the converter system changes from F2F2 at Iref ∈ [0.8296,0.9006] to F2F3 at Iref ∈ [0.9007,1.1577]and to F1F3 at Iref ∈[1.1578,1.1946].From the entropy value of three period-2 orbits shown in Figure 7, the permutation entropy value of period-2 is greater than the corresponding entropy value of any working pattern in the entirety of period-1.The sample entropy, used for assessing the complexity of the time sequence shown in Figure 9, is the negative logarithm of the probability of if two sets of simultaneous data points of length m have a distance smaller 1.
As I ref ∈ [0.1200,0.2778], the system works in the F 3 pattern of period-1.In this pattern, there are three modes in each period of the system, which means the system is in the DCM state.Then, the first border collision bifurcation occurs in the system at I ref = 0.2779 A. Because the inductor current in the steady state collides with Ib 2 , the operating mode changes and the working patterns of the converter system change from F 3 to F 2 .

2.
During I ref ∈[0.2779,0.8295], the system still works in the period-1 state, but the operating mode is F 2 .At this time, the system has two modes in each period.When I ref = 0.8296, the system generates period doubling bifurcation and enters period-2.Since there is no border collision bifurcation, the working patterns are extended to F 2 F 2 .

3.
When I ref ∈[0.8296,1.1946], the system works in period-2.We found that two bounder collisions of the system occur in this area because the inductance current of the system collides with the borderline of Ib From the entropy value of three period-2 orbits shown in Figure 7, the permutation entropy value of period-2 is greater than the corresponding entropy value of any working pattern in the entirety of period-1.The sample entropy, used for assessing the complexity of the time sequence shown in Figure 10, is the negative logarithm of the probability of if two sets of simultaneous data points of length m have a distance smaller than r, then two sets of simultaneous data points of length m + 1 also have a distance smaller than r [28].However, with the change in sample entropy, we could not find the bifurcation information of the system under period-1 and period-2 states.4.
As I ref ∈[1.1947,1.2624], the system works in the period-4 state with the working pattern F 1 F 2 F 1 F 3 , and there are seven modes in this state.Then, the fifth time border collision bifurcation of a fixed point and Ib 1 curve collision occur once more in the system at I ref = 1.2625 A. At this time, LLE shown in Figure 9 is greater than 0 and the mode switching shows disorder, which means that the system is beginning enter the chaotic state.When I ref ≥ 1.2625 A, a number of border collisions occur in the chaotic state, as depicted in Figure 7.Then, the value of permutation entropy is no longer stable, exhibiting similar characteristics to the LLE, and its corresponding value is greater than any period.The characteristics described are consistent with Bandt C et al. [20].
Moreover, the fourth border collision bifurcation occurs at I ref = 1.1947 A. Since the system has unstable characteristic roots at the same time, the system generates period-doubling bifurcation and enters the period-4 state.The working pattern is extended while F 2 switches to F 3 .
border collisions occur in the chaotic state, as depicted in Figure 7.Then, the value of permutation entropy is no longer stable, exhibiting similar characteristics to the LLE, and its corresponding value is greater than any period.The characteristics described are consistent with Bandt C et al. [20].
Moreover, the fourth border collision bifurcation occurs at Iref = 1.1947 A. Since the system has unstable characteristic roots at the same time, the system generates period-doubling bifurcation and enters the period-4 state.The working pattern is extended while F2 switches to F3.  than r, then two sets of simultaneous data points of length m + 1 have a distance smaller than r [28].However, with the change in sample entropy, we could not find the bifurcation information of the system under period-1 and period-2 states.4. As Iref ∈ [1.1947,1.2624], the system works in the period-4 state with the working pattern F1F2F1F3, and there are seven modes in this state.Then, the fifth time border collision bifurcation of a fixed point and Ib1 curve collision occur once more in the system at Iref = 1.2625 A. At this time, LLE shown in Figure 10 is greater than 0 and the mode switching shows disorder, which means that the system is beginning enter the chaotic state.When Iref ≥ 1.2625 A, a number of border collisions occur in the chaotic state, as depicted in Figure 7.Then, the value of permutation entropy is no longer stable, exhibiting similar characteristics to the LLE, and its corresponding value is greater than any period.The characteristics described are consistent with Bandt C et al. [20].
Moreover, the fourth border collision bifurcation occurs at Iref = 1.1947 A. Since the system has unstable characteristic roots at the same time, the system generates period-doubling bifurcation and enters the period-4 state.The working pattern is extended while F2 switches to F3.To further explain permutation entropy, the time-domain simulating waveforms with phase plane, and the experimental results at the typical values of I ref, are illustrated in Figures 11 and 12, respectively.We found that the values of the permutation entropy correspond to the complexity of the working waveforms.In summary, the presented MSTS-PE can quantify the complexity of the operation states of the buck converter, and the established MSTS can characterize the bifurcation phenomena of the system.To further explain permutation entropy, the time-domain simulating waveforms with phase plane, and the experimental results at the typical values of Iref, are illustrated in Figures 11 and 12, respectively.We found that the values of the permutation entropy correspond to the complexity of the working waveforms.In summary, the presented MSTS-PE can quantify the complexity of the operation states of the buck converter, and the established MSTS can characterize the bifurcation phenomena of the system.

Discussion and Conclusions
Phase space reconstruction is a method of reconstructing attractors based on finite data to study the dynamic behavior of a system.Its basic idea is that the evolution of any component in the system is determined by its interactions with other components.In DC-DC converters, the dynamic information of one component is implicit in any related component; that is, the original dynamic system model can be reconstructed with one variable observation in the system.
The quantization process of the sequence information based on the permutation entropy, and comparison between permutation entropy and the traditional symbol sequence method, are described in Section 3. In addition, the permutation entropy and sample entropy of buck converter are shown in Figures 7 and 9.The permutation entropy seems to be a better indicator in some circumstances because it has high robust performance and anti-interference ability for analysis of the mixed signal.Permutation entropy overcomes shortcomings, like the poor relative consistency of sample entropy and time consuming calculation of multi-scale entropy, which is easily affected by the non-stationarity and outliers of time sequences.For a better description of nonlinear dynamic behaviors of a system from a time sequence, MSTS was proposed to describe the dynamic features of a DC-DC converter.In addition, distinguished from joint entropy, permutation entropy combined with MSTS only needs one sequence to identify the system features.This concept is also consistent with the concept of phase space reconstruction.
Based on the theory of phase space reconstruction, the permutation entropy of the DC-DC converter under the determined parameters is obtained by sorting each group of m-dimensional vectors.In this paper, we found that the higher the frequency of mode switching, the more complex the obtained MSTS.This means that the higher the complexity of the corresponding system, the larger the permutation entropy value.With the fluctuation in the permutation entropy value of different control parameters, the nonlinear dynamic behaviors of a DC-DC converter can be better quantified.

Discussion and Conclusions
Phase space reconstruction is a method of reconstructing attractors based on finite data to study the dynamic behavior of a system.Its basic idea is that the evolution of any component in the system is determined by its interactions with other components.In DC-DC converters, the dynamic information of one component is implicit in any related component; that is, the original dynamic system model can be reconstructed with one variable observation in the system.
The quantization process of the sequence information based on the permutation entropy, and comparison between permutation entropy and the traditional symbol sequence method, are described in Section 3. In addition, the permutation entropy and sample entropy of buck converter are shown in Figures 7 and 10.The permutation entropy seems to be a better indicator in some circumstances because it has high robust performance and anti-interference ability for analysis of the mixed signal.Permutation entropy overcomes shortcomings, like the poor relative consistency of sample entropy and time consuming calculation of multi-scale entropy, which is easily affected by the non-stationarity and outliers of time sequences.For a better description of nonlinear dynamic behaviors of a system from a time sequence, MSTS was proposed to describe the dynamic features of a DC-DC converter.In addition, distinguished from joint entropy, permutation entropy combined with MSTS only needs one sequence to identify the system features.This concept is also consistent with the concept of phase space reconstruction.
Based on the theory of phase space reconstruction, the permutation entropy of the DC-DC converter under the determined parameters is obtained by sorting each group of m-dimensional vectors.In this paper, we found that the higher the frequency of mode switching, the more complex the obtained MSTS.This means that the higher the complexity of the corresponding system, the larger the permutation entropy value.With the fluctuation in the permutation entropy value of different control parameters, the nonlinear dynamic behaviors of a DC-DC converter can be better quantified.

Figure 1 .
Figure 1.Sampling at the mode switching time to build a mode switching time sequence (MSTS.)

Figure 1 .
Figure 1.Sampling at the mode switching time to build a mode switching time sequence (MSTS).

Figure 2 .
Figure 2. Procedure to obtain the permutation symbols from the sample period-6 signal as delay time τ = 2.

Figure 2 .
Figure 2. Procedure to obtain the permutation symbols from the sample period-6 signal as delay time τ = 2.

Figure 4 .
Figure 4. Buck converter with peak current control.Figure 4. Buck converter with peak current control.

Figure 4 .
Figure 4. Buck converter with peak current control.Figure 4. Buck converter with peak current control.

Figure 6 .
Figure6.(a) F 1 working pattern means that only mode-1 exists in one period; (b) F 2 means that mode-1 is performed first in the period, then the switch state changes once and enters into mode-2; and (c) F 3 represents the mode in the period change from mode-1 to mode-2 to mode-3.

Figure 7 .
Figure 7. Permutation entropy of buck converter with Iref as the changing parameter.

Figure 8 .
Figure 8. Bifurcation diagram of a buck converter with Iref as the changing parameter.

Figure 7 .Figure 7 .
Figure 7. Permutation entropy of buck converter with I ref as the changing parameter.

Figure 8 .
Figure 8. Bifurcation diagram of a buck converter with Iref as the changing parameter.

Figure 8 .
Figure 8. Bifurcation diagram of a buck converter with I ref as the changing parameter.

2 and
Ib 1 at I ref = 0.9007 and I ref = 1.1578A, respectively.Hence, the operating pattern of the converter system changes from F 2 F 2 at I ref ∈[0.8296,0.9006] to F 2 F 3 at I ref ∈[0.9007,1.1577]and to F 1 F 3 at I ref ∈[1.1578,1.1946].

Figure 9 .Figure 10 .
Figure 9. Sample Entropy of buck converter with Iref as the change parameter

Figure 9 .
Figure 9. Largest Lyapunov exponent (LLE) of a buck converter with I ref as the changing parameter.

Figure 9 .
Figure 9. Sample Entropy of buck converter with Iref as the change parameter

Figure 10 .
Figure 10.Sample Entropy of buck converter with I ref as the change parameter

Figure 11 .Figure 12 .
Figure 11.Waveforms of inductor current x L and phase plane at (a) I ref = 0.20 A; (b) I ref = 0.50 A; (c) I ref = 0.85 A; (d) I ref = 1.00 A; (e) I ref = 1.17 A; (f) I ref = 1.20 A; and (g) I ref = 1.27 A. In (a) to (g), x C and x L represent the inductor current and capacitor voltage, respectively, and the red marker dot corresponds the mode switching moment.

Figure 12 .
Figure 12.Experimental waveforms of inductor current x L and switching signal at (a) I ref = 0.20 A; (b) I ref = 0.50 A; (c) I ref = 0.85 A; (d) I ref = 1.00 A; (e) I ref = 1.17 A; (f) I ref = 1.20 A; and (g) I ref = 1.27 A.

Table 1 .
Permutation sequence and permutation entropy of the buck converter as reference current Iref changes.

Table 1 .
Permutation sequence and permutation entropy of the buck converter as reference current I ref changes.