Modeling and Harmonic Analysis of a Fractional-Order Zeta Converter

: The Zeta converter is an essential and widely used high-order converter. The current modeling studies on Zeta converters are based on the model that devices, such as capacitors and inductors, are of integer order. For this reason, this paper takes the Zeta converter as the research object and conducts an in-depth study on its fractional-order modeling. However, the existing modeling and analysis methods have high computational complexity, the analytical solutions of system variables are tedious, and it is difﬁcult to describe the ripple changes of state variables. This paper combines the principle of harmonic balance with the equivalent small parameter method (ESPM); the approximate analytic steady-state solution of the state variable can be obtained in only three iterative steps in the whole solving process. The DC components and ripples of the state variables obtained by the proposed method were compared with those obtained by the Oustaloup’s ﬁlter-based approximation method; the symbolic period results obtained by ESPM had sufﬁcient precision because they included more combinations of higher harmonics. Finally, the inﬂuence of fractional order on harmonics were analyzed. The obtained results show that the proposed method has the advantage of being less computational and easily describing changes in the ripple of the state variables. The simulation results are provided for validity of the theoretical analysis.


Introduction
Since the development of electronic converters of a high frequency and high integration of power, the scale and complexity of the systems continue to increase.Thus, the errors caused by modeling gradually cannot be ignored.Inductance and capacitance are essential energy storage components in DC-DC converters.Most modeling and analysis methods for DC-DC converters is based on integer order inductors and capacitors.The literature [1,2] shows that inductance and capacitance are both fractional order in nature, which means that errors may occur inevitably when using the integer-order model to analyze inductance and capacitance.Research on more accurate modeling and analysis methods for DC-DC converters have important theoretical significance and practical application value.
Low-order converters, such as buck, boost, and buck-boost converters, have the advantages of a small size, simple structure, and high conversion efficiency, which can easily achieve boost, step-down, and negative voltage output.They are widely used in inverter circuits and power factor correction circuits (PFC).References [3][4][5] use fractional-order calculus theory to model and analyze boost converters operating in continuous conduction mode (CCM), discontinuous conduction mode (DCM), and pseudo-continuous conduction mode (PCCM), respectively.In references [6,7], the fractional-order mathematical model and the state average model of CCM and DCM buck converters were established and analyzed, respectively, based on fractional calculus theory.The research results show that the fractional-order models of boost converters and buck converters based on fractionalorder calculus theory differ greatly from those of integer-order models.Compared with Energies 2023, 16, 3969 2 of 18 the integer-order model, the fractional-order model is more accurate and more consistent with the actual circuit essence.At present, there are mainly four definitions of fractional calculus.The research result under different definitions of fractional order calculus is also different.References [8,9] conducted a steady-state analysis of buck converters and buck-boost converters in CCM mode based on Riemann Liouville's (R-L) calculus definition and compared the results with the analysis results defined by Caputo.Compared with the Caputo definition, using the R-L definition to analyze the fractional buck converter and buck-boost converter is more accurate, and better in tracking dynamic response.However, different definitions of fractional-order calculus are not as unified as those of integer-order calculus, and it is difficult to give the expression of the analytical solution of the system according to these definitions [10,11].Therefore, numerical algorithms can be used to solve and analyze fractional-order systems.The existing modeling methods of the fraction-order converter include the state space averaging method, predictive correction method (ABM), and the Oustaloup's frequency-domain filtering method.The state-space averaging method takes the single period of the system as the research object.According to the state-space equations of the system equations in different modes, weighted average processing is performed on the linear switching function in each period, and then through the small signal disturbance processing and linearization processing, the equivalent model of the switching converter is obtained.Since the state space averaging method is "averaging" processing, and only considers the low-frequency characteristics of the system, ignoring the high-frequency dynamic characteristics of the switching converter, which will fail to correctly analyze the steady-state characteristics of the system, especially considering that the inductance and capacitance are fractional order, the inductance current ripple and capacitance voltage will be larger than the integer order.So, the error is bigger in fractional circuits.The basic idea of ABM is to establish the time model of the system, obtain the fractional differential equation of the system, then to discretize the differential of the equation, and then obtain the near value of the numerical solution through the pre-estimator, and the initial value of the next iteration can be obtained by using the corrector obtained in the calculation process.After many iterations, the numerical solution gradually converges to a stable numerical interval.The advantage of the ABM algorithm is its simple and effective nature, but due to the restriction of time complexity, the selection of initial value and the dimension of variables will affect the system's analysis of the research.The Oustaloup frequency-domain filtering method can clearly express the relationship between system state variables and circuit element parameters when analyzing fractional order systems.However, the Oustaloup's frequency-domain filtering method also has some problems, such as difficulty in modeling high-order converters, slow running speed and the impact of switching frequency on simulation results.The lower limit of fitting frequency, upper limit of fitting frequency and the order of filter also have great influence on the accuracy of fractional-order DC-DC converter modeling.The fractional calculus operation has long-term memory characteristics, which results in the numerical algorithm occupying a large memory space and the operational process being time-consuming in simulating computer software [12].References [13,14] use ESPM to model and analyze the steady-state characteristics of fractional-order buck and boost converters in CCM and DCM.Using the ESPM, the approximate steady-state analytical solution of the state variables can be obtained without considering the complex definition of fractional calculus, and the amplitudes of each harmonic in the circuit can be calculated.Compared with the numerical solution, the analytical solution obtained by ESPM is more general and can reflect the relationship between the system state variables better and more intuitively.In addition, the ESPM has the advantages of fast running speed and less memory.Now, the fractional order modeling and analysis of converters are mainly focused on low-order converters.However, high-order converters (Cuk, Sepic, Zeta) can achieve higher voltage gain, small ripple, small volume, and higher transmission efficiency, so they are widely used in wind power generation systems, fuel cells, and photovoltaic systems [15][16][17].However, the modeling and analysis methods for high-order converters Energies 2023, 16, 3969 3 of 18 are very complicated [18].The research and analysis on the modeling of Zeta converters has been achieved to some extent, but the research results are mainly based on the model that the capacitors and inductors in the circuit are of integer order.Therefore, in this paper, we take the Zeta converter as the research object and conduct an in-depth study on its fractional-order modeling.Based on the existing research, we consider the fractional-order characteristics of capacitors and inductors in the circuit and analyze the characteristics of fractional-order circuits, which not only can represent the real characteristics of the converter more accurately, but also can use the nature of fractional-order calculus to expand its application scope and further improve the performance of the converter.In this way, we can improve the theoretical basis of power electronic converter circuit research.For this reason, this paper aims to model and analyze the fractional order Zeta converter operating in CCM by using ESPM and proposes a method for obtaining approximate analytical solutions of the state variables of the inductor current and capacitor voltage.According to the obtained solutions, the DC quantity, ripple, and harmonics of the state variables can be analyzed intuitively, and then the relative performance of the DC chopper can be studied.
This paper is organized as follows.Section 2 deals with the modeling of the fractionalorder Zeta converter in CCM.Section 3 establishes an equivalent model based on the ESPM method and the approximate solution.In Section 4, simulations are performed to verify the proposed method.The effects of fractional order on the harmonic amplitude of capacitor voltage and inductor current are discussed in Section 5. Finally, conclusions are given in Section 6.

Mathematical Model of the Fractional CCM Zeta Transformer
The circuit topology of the fractional Zeta converter is shown in Figure 1.In the figure, V in is the power supply voltage, δ(t) is the switching function, S T and S D are the switching tubes, R is the circuit resistance, and i L1 , i L2 , v C1 and v C2 are the inductive current and capacitor voltages of the converter, respectively.α and β are the orders of the inductor and capacitor, respectively.
higher voltage gain, small ripple, small volume, and higher transmission efficiency, so they are widely used in wind power generation systems, fuel cells, and photovoltaic systems [15][16][17].However, the modeling and analysis methods for high-order converters are very complicated [18].The research and analysis on the modeling of Zeta converters has been achieved to some extent, but the research results are mainly based on the model that the capacitors and inductors in the circuit are of integer order.Therefore, in this paper, we take the Zeta converter as the research object and conduct an in-depth study on its fractional-order modeling.Based on the existing research, we consider the fractional-order characteristics of capacitors and inductors in the circuit and analyze the characteristics of fractional-order circuits, which not only can represent the real characteristics of the converter more accurately, but also can use the nature of fractional-order calculus to expand its application scope and further improve the performance of the converter.In this way, we can improve the theoretical basis of power electronic converter circuit research.For this reason, this paper aims to model and analyze the fractional order Zeta converter operating in CCM by using ESPM and proposes a method for obtaining approximate analytical solutions of the state variables of the inductor current and capacitor voltage.According to the obtained solutions, the DC quantity, ripple, and harmonics of the state variables can be analyzed intuitively, and then the relative performance of the DC chopper can be studied.
This paper is organized as follows.Section 2 deals with the modeling of the fractional-order Zeta converter in CCM.Section 3 establishes an equivalent model based on the ESPM method and the approximate solution.In Section 4, simulations are performed to verify the proposed method.The effects of fractional order on the harmonic amplitude of capacitor voltage and inductor current are discussed in Section 5. Finally, conclusions are given in Section 6.

Mathematical Model of the Fractional CCM Zeta Transformer
The circuit topology of the fractional Zeta converter is shown in Figure 1.In the figure, Vin is the power supply voltage, δ(t) is the switching function, ST and SD are the switching tubes, R is the circuit resistance, and iL1, iL2, vC1 and vC2 are the inductive current and capacitor voltages of the converter, respectively.α and β are the orders of the inductor and capacitor, respectively.
The fractional-order inductor and the voltage across the capacitor have the following relationship: The fractional-order inductor and the voltage across the capacitor have the following relationship: When both α and β are equal to 1, the object being modeled is a traditional integer order model.
When the converter operates in the continuous current mode (CCM), the switches S T and S D are controlled by the switching function δ(t).The fractional-order Zeta converter has two operating modes in one switching period T S .The inductor current and capacitor voltage are used as the state variables of the converter; fractional-order differential equations describe the state variables of the two modes of the converter: Mode 1: nT S < t ≤ (n + D)T S , where D is the on-time duty cycle and n is an integer representing a certain switching period.As shown in Figure 2, at this time δ(t) = 1, S T is turned on and S D is turned off, and the state equation is written: Working mode 1 of fractional-order CCM-Zeta converter.
SS n D T t n T , at this time δ(t) = 0, ST is turned off and SD is turned on.As shown in Figure 3, the state equation is written as: Working mode 2 of fractional-order CCM-Zeta converter.
Mode 2: (n + D)T S < t ≤ (n + 1)T S , at this time δ(t) = 0, S T is turned off and S D is turned on.As shown in Figure 3, the state equation is written as: turned on and SD is turned off, and the state equation is written: Working mode 1 of fractional-order CCM-Zeta converter.
SS n D T t n T , at this time δ(t) = 0, ST is turned off and SD is turned on.As shown in Figure 3, the state equation is written as: Working mode 2 of fractional-order CCM-Zeta converter.
In this paper, the operator symbols p α and p β are used to replace d α /dt α and d β /dt β , respectively; the nonlinear switching function δ(t) is used to represent the different modes of the fractional-order CCM Zeta converter, which is defined as: According to the mathematical model of the two modes of Equations ( 2) and ( 3) and the switching function δ(t), the differential equation of the fractional CCM Zeta converter can be obtained as Representing Formula (5) in the form of a matrix equation, we have: There are nonlinear functions 6), and the coefficient matrix The above is the state vector model of the fractional-order Zeta converter working in CCM.It can be observed that the coefficient matrix G 0 (p α1 , p α2 , p β1 , p β2 ) in the model is closely related to the order of the energy storage element of the converter.

Principle of the ESPM
The ESPM is mainly based on the harmonic balance principle of solving nonlinear differential equations.By expanding the state variables of the circuit system by Fourier series, when there are enough expansion terms of the series, the sum of these series can be outstanding to approximate the periodic solution of the state variable.According to the literature [19], the generalization of the harmonic balance principle is as follows: where µ is the order of the differential operation, which can be either an integer or a non-integer; and ω is the angular frequency of the exponential function, which can be any real number.From the above formula, it can be observed that the order of the differential operation only affects the amplitude of the differentiated exponential function but has no effect on its phase.Therefore, the harmonic balance principle can be generalized to solve fractional-order nonlinear differential equations.Furthermore, the ESPM can also be improved according to this generalization.First, the state variable x is expanded to the form of the sum of the main oscillation component x 0 and the other corrections x i : The switching function δ(t) is also expanded into a similar series sum form: Energies 2023, 16, 3969 Substituting the expansions of x and δ(t) into the nonlinear functions f (x) = δ(t)x and t(x) = δ(t)e, and then combining terms with the same order ε i , f (x) and t(x) can be converted to: where f i is: According to the principle of harmonic balance, it can be assumed that x i can be expressed in the form of the sum of various harmonic components shown in Equation ( 15): where c.c represents the conjugate term, and k represents the order of the harmonic, in which τ = ωt is the normalized time.a i0 and a ik represent the DC component of x i and the amplitude of the kth harmonic component, respectively.In addition, the formula's harmonic components {E ir } are determined by the actual physical properties of the modeled object.For the DC-DC converter, due to its low-pass filtering characteristics, the main component in its state variable is DC, so set the harmonic component set {E 0 } = {0} of the main oscillation component x 0 .The harmonic components {E ir } of the other corrections x i are sequentially determined by the remainder R i generated by the previous operation.
When the converter is in a steady state, the periodic switching function δ(t) can also be expanded into a Fourier series of the form: In the formula, b 0 = 1 , where: Usually, δ 0 and δ i are taken as: Substitute the expansions of x i and δ i into f i and t i .Due to the multiplication of x i and δ i , new harmonic components are generated in f i and t i , as follows: Among them, f im represents the part of f i that is the same as the harmonic component of x i , the remainder R i+1 represents other components in f i , t im represents the part of t i that is the same as the harmonic component of x i , and the remainder H i+1 represents other components in t i components, the spectrum of x i and R i are the same as H i .Therefore, the value of the x i spectrum can be determined by R i and H i in the operation.
Substituting Equations ( 15) and ( 19) into Equation ( 6), Equation ( 20) can be deduced as: From Equation ( 21), the main oscillation component and each order correction can be obtained step by step.If the following conditions are met, the iteration will be stopped.
By solving the equation in Equation ( 20), the main oscillation component and each order correction can be obtained step by step, and finally the periodic solution of the system state variable can be expressed as: Due to the low-pass performance of the DC-DC converter, it can be concluded that the main oscillation component of the state variable is: where I 00 , I 00 , V 00 , and V 00 represent the DC components of the inductor current and capacitor voltage, respectively.From Equation (19) f i and t i expressions, it can be expressed as: Substituting Equations ( 24) and (25) into the first term of Equation ( 20), it can be expressed as: At this time, the coefficient matrix is G 00 = G 0 (0,0,0,0).According to the principle of harmonic balance, the above formula can be transformed into The analytical formula of a 00 can be obtained from Formula (27).The components of x 1 depend on R 1 and H 1 , so the set of harmonic components of the first-order correction amount x 1 is {E 1r } = {1}, and x 1 can be expressed as Where a 11 = [I 11 , I' 11 , V 11 , V' 11 ] T , bring x 0 , x 1 , δ 0 , δ 1 into f 1 and t 1 , and the remainder R 2 of f 1m and f 1 can be obtained: Since the spectrum of H 1 is {K 1r } = {1}, the spectrum of x 1 is also {1}.In the same way, selecting the same amount as the spectrum of x 1 in t 1 as t 1m , and the other components as H 2 , it can be expressed as: Substituting f 1m , t 1m , H 1 , and R 1 into the second equation of Equation ( 20), it can be reverted to: where G 01 = G 0 [(jω) α1 , (jω) α2 , (jω) β1 , (jω) β2 ], a 11 can be obtained by the matrix transformation of Formula (31).
According to the above process, a ik = [I ik , I' ik , V ik , V' ik ] can be obtained.According to the obtained main oscillation component x 0 , first-order correction amount x 1 , second-order correction amount x 2 and third-order correction amount x 3 , then the analytical expression of the state vector x can be obtained as: x = a 00 + a 20 + (a 11 + a 31 )e jτ + a 22 e j2τ + a 23 e j3τ + a 34 e j4τ + a 35 e j5τ + c.c The expressions of state variables i L1 , i L2 , v C1 , and v C2 are: where Re and Im represent the real and imaginary parts, respectively.

Design Equation of Zeta Converter
Applying Kirchhoff's voltage law on the Zeta converter circuit for the first and second mode, the equations are derived below.The ripple of the current through the energy transferring (input) inductor can be expressed as, The output inductor current ripple can be expressed as, The capacitor ripple voltages ∆V C1 & ∆V C2 can be derived from the Kirchhoff's current law for first and second mode as where The simulation parameters of Figure 1 are shown in Table 1.

DC Components and Ripple Analysis
The fractional-order CCM Zeta converter is modeled on the Oustaloup's filter-based approximation method [20], and the simulation model is shown in Figure 4.
The simulation parameters of Figure 1 are shown in Table 1.

DC Components and Ripple Analysis
The fractional-order CCM Zeta converter is modeled on the Oustaloup's filter-based approximation method [20], and the simulation model is shown in Figure 4. Using the Oustaloup's filter-based approximation method, the fractional-order inductor and capacitor in Figure 4 are replaced by the equivalent circuit of the fractionalresistance chain.The fractional-resistance chain of fractional-order inductor is constructed, as shown in Figure 5a.The fractional-resistance chain of fractional-order capacitor is shown in Figure 5b.Using the Oustaloup's filter-based approximation method, the fractional-order inductor and capacitor in Figure 4 are replaced by the equivalent circuit of the fractional-resistance chain.The fractional-resistance chain of fractional-order inductor is constructed, as shown in Figure 5a.The fractional-resistance chain of fractional-order capacitor is shown in Figure 5b.The DC components of the state variables obtained by the Oustaloup's filter-based approximation method and ESPM are shown in Table 2.It can be seen from Table 2 that the difference between the results of the inductor current and capacitor voltage DC components obtained by ESPM and the Oustaloup's filter-based approximation method is relatively small, and the changing trends are the same.
The waveform of inductor current and capacitor voltage of different orders obtained by the Oustaloup's filter-based approximation method and ESPM are shown in Figure 6 and Figure 7, respectively.The figure's green solid line and the red dotted line in the figure are the curves obtained by the Oustaloup's filter-based approximation method and ESPM, respectively.The DC components of the state variables obtained by the Oustaloup's filter-based approximation method and ESPM are shown in Table 2.It can be seen from Table 2 that the difference between the results of the inductor current and capacitor voltage DC components obtained by ESPM and the Oustaloup's filter-based approximation method is relatively small, and the changing trends are the same.
It can be seen from Figures 6 and 7 that the values of the inductor current and capa itor voltage are consistent with the trends shown in Table 2. Furthermore, the steady-sta ripples from the Oustaloup's approximation method and the ESPM are compared, whe the green solid lines represent the results from the Oustaloup's approximation metho and the red dotted lines represent the results from the ESPM.The harmonic magnitude are also order dependent.Specifically, the harmonic magnitudes of iL and vC increase wi the decreases of α and β, respectively.Waveforms coming from these two methods sho good agreement with each other, which proves the accuracy and feasibility of ESPM mod eling of the fractional CCM Zeta converter.
It can be seen from Figures 6 and 7 that the values of the inductor current and capacitor voltage are consistent with the trends shown in Table 2. Furthermore, the steady-state ripples from the Oustaloup's approximation method and the ESPM are compared, where the green solid lines represent the results from the Oustaloup's approximation method, and the red dotted lines represent the results from the ESPM.The harmonic magnitudes are also order dependent.Specifically, the harmonic magnitudes of i L and v C increase with the decreases of α and β, respectively.Waveforms coming from these two methods show good agreement with each other, which proves the accuracy and feasibility of ESPM modeling of the fractional CCM Zeta converter.
The inductance current increment of the circuit model can be measured from Figures 6  and 7, and the corresponding theoretical values can be obtained according to the formula for ripple given in Equation ( 46), and the comparison between the theoretical and simulated values is shown in Table 3.The ripple error percentage is shown in Table 4. Combining Tables 2-4 it can be seen that the DC components of the inductor current and capacitor voltage are basically not affected by the fractional order in a steady state.The fractional order mainly affects the ripple of the inductor current and capacitor voltage, and the ripple amplitude increases sharply as the order decreases.Analysis of the data in Table 3 shows that there is a certain error between the simulated and theoretical values, but the error is small and stays in a reasonable range.The comparison between the theoretical and simulated values fully verifies the correctness of the fractional order model and the theoretical derivation.

Fractional-Order Zeta Converter CCM Discriminant
From the method of reference [21], the expressions for the fractional order inductance currents ∆i L1 and ∆i L2 can be obtained as: In Equation ( 46), Г(α) is the gamma function.Additionally, V C1 = −V C2 , so by taking L = L 1 = L 2 , α = α 1 = α 2 , the fractional order inductor L 1 α1 and L 2 α2 current ripple calculation formula can be expressed uniformly as follows.
Equation ( 47) shows that when the fractional-order inductors L 1 α1 and L 2 α2 have equal inductance values and orders, the inductor current ripple is calculated by the same formula, and the inductor current ripple is related to its order.The curve of inductor current ripple ∆i L with order α can be obtained by using Matlab software, as shown in Figure 8.As the order decreases, the inductor current ripple ∆i L increases sharply, and the lower the order, the faster the ripple increases.For a fractional-order Zeta converter to operate in CCM, the current flowing through the converter diode SD must always be greater than 0 in the ( )  The diode current iD is always greater than 0. Equivalently, the DC component of the inductor current should always be greater than half of its current ripple amplitude, yielding the following expression: For a fractional-order Zeta converter to operate in CCM, the current flowing through the converter diode S D must always be greater than 0 in the (n + D)T S < t ≤ (n + 1)T S interval.In the (n + D)T S < t ≤ (n + 1)T S interval, the current flowing through the diode is equal to the sum of the two inductor currents.
The diode current i D is always greater than 0. Equivalently, the DC component of the inductor current should always be greater than half of its current ripple amplitude, yielding the following expression: Substituting Equation (46) into Equation ( 49) and simplifying it can be expressed in the following form: The effect of order α 1 and α 2 of the fractional-order inductor on the converter operation mode can be obtained from Equation (50), as shown in Figure 9.
the converter diode SD must always be greater than 0 in the ( ) The diode current iD is always greater than 0. Equivalently, the DC component of the inductor current should always be greater than half of its current ripple amplitude, yielding the following expression: Substituting Equation (46) into Equation ( 49) and simplifying it can be expressed in the following form: The effect of order α1 and α2 of the fractional-order inductor on the converter operation mode can be obtained from Equation (50), as shown in Figure 9.  Figure 9a shows the three-dimensional diagram of the current relative to the order α 1 -α 2 plane change of Equation (50).When this current is greater than 0, it indicates that the converter is operating in CCM mode; when this current is less than 0, it indicates that the converter is operating in DCM mode; when this current is equal to 0, it indicates that it is operating in the critical state of CCM and DCM; the boundary lines of the two operating modes are shown in Figure 9b.From Figure 9b, the intersection coordinates (0.71, 1) and (1, 0.71) of the dividing line and the edge of α 1 -α 2 plane can be obtained.Only when the inductance order is greater than 0.71 order can the converter run in CCM mode.

Analysis of Harmonic Components in Different Orders
The steady-state approximate analytical solution of the fractional-order CCM Zeta converter is obtained through the analytical modeling method in Section 2. It is observed that both the DC component and the ripple of the energy storage element are related to the order of the energy storage element in the study.The specific change rule is as follows: when the order of the energy storage element decreases, the ripple amplitude of the state variable increases; in the contrary, the ripple amplitude of the state variable decreases.When the inductor and capacitor are in different orders, the harmonic amplitudes of the converter inductor current and capacitor voltage are obtained, as shown in Figures 10  and 11, respectively.that both the DC component and the ripple of the energy storage element are related to the order of the energy storage element in the study.The specific change rule is as follows: when the order of the energy storage element decreases, the ripple amplitude of the state variable increases; in the contrary, the ripple amplitude of the state variable decreases.When the inductor and capacitor are in different orders, the harmonic amplitudes of the converter inductor current and capacitor voltage are obtained, as shown in Figures 10 and  11, respectively.(0.95,0.95,1,1) (0.95,0.95,0.95,0.95)(0.9,0.9,0.95,0.95)(0.9,0.9,0.9,0.9)(0.85,0.85,0.85,0.85)

RMS values of capacitor voltage(V)
Fundamental Harmonic 2nd Harmonic 3rd Harmonic 4th Harmonic 5th Harmonic From Figures 10 and 11, as the order of the inductor and capacitor elements decreases, the amplitudes of the fundamental, second, third, fourth, and fifth harmonics of the inductor current and capacitor voltage increase.Since the ripple of the state variable is jointly affected by the harmonics, the changing trend of the harmonics in the state variable is consistent with the changing trend of the ripple amplitude.The change in the inductance order greatly influences the inductor current, while the change in the capacitance order has a relatively small influence on the inductor current.
The modified termination criterion of the analytical solution proposed in Equation ( 21) is verified by the change in each harmonic amplitude.The condition of the discriminant stability were analyzed with respect to the harmonic components.The effective values of the harmonic components of the converter inductor current and capacitor voltage, as a percentage of the fundamental component, were obtained, as shown in Figure 12.As can be seen from Figure 12, the value of the modified termination discriminant condition increases as the fractional order of the inductor and capacitor decreases.(21) is verified by the change in each harmonic amplitude.The condition of the discriminant stability were analyzed with respect to the harmonic components.The effective values of the harmonic components of the converter inductor current and capacitor voltage, as a percentage of the fundamental component, were obtained, as shown in Figure 12.As can be seen from Figure 12, the value of the modified termination discriminant condition increases as the fractional order of the inductor and capacitor decreases.

Conclusions
Based on the extended harmonic balance principle and the idea of disturbance, ESPM was used to model the fractional CCM zeta converter.The steady-state analytical expression was obtained.According to the steady-state analytical expression, the fractional-order CCM Zeta converter was modeled in MATLAB and compared with the Oustaloup's filter-based approximation method.Furthermore, the harmonic amplitude of each order of the state variables of the fractional CCM Zeta converter was obtained, and the influence of the order of the energy storage element on the harmonic components of each order was analyzed.The obtained results show that: (1) ESPM can avoid the discussion of the applicability of several fractional calculus definitions to the upper and lower limits of the integration under different circumstances, and overcome the problem that it is difficult for the fractional system to obtain specific expressions.The obtained solutions conclude practical physical significance, and the analysis results are consistent with those obtained by the Oustaloup's filter-based approximation method.(2) The amplitude of each harmonic of the fractional converter is related to the order of the inductance and capacitance components.With all other parameters unchanged, when the fractional order of inductance and capacitance decreases, the amplitude of 18.00% (0.95,0.95,1,1) (0.95,0.95,0.95,0.95)(0.9,0.9,0.95,0.95)(0.9,0.9,0.9,0.9)(0.85,0.85,0.85,0.85) Inductance current harmonic component as a percentage of DC component(%) Fundamental Harmonic 2nd Harmonic 3rd Harmonic

Conclusions
Based on the extended harmonic balance principle and the idea of disturbance, ESPM was used to model the fractional CCM zeta converter.The steady-state analytical expression was obtained.According to the steady-state analytical expression, the fractional-order CCM Zeta converter was modeled in MATLAB and compared with the Oustaloup's filter-based approximation method.Furthermore, the harmonic amplitude of each order of the state variables of the fractional CCM Zeta converter was obtained, and the influence of the order of the energy storage element on the harmonic components of each order was analyzed.The obtained results show that: (1) ESPM can avoid the discussion of the applicability of several fractional calculus definitions to the upper and lower limits of the integration under different circumstances, and overcome the problem that it is difficult for the fractional system to obtain specific expressions.The obtained solutions conclude practical physical significance, and the analysis results are consistent with those obtained by the Oustaloup's filter-based approximation method.(2) The amplitude of each harmonic of the fractional converter is related to the order of the inductance and capacitance components.With all other parameters unchanged, when the fractional order of inductance and capacitance decreases, the amplitude of the harmonic components of each order in the state variable increases, increasing the amplitude of the inductor current and capacitance voltage ripple of the fractional order converter.(3) Compared with the numerical simulation method, the proposed method can better describe the change in the state variable ripple.The computational complexity is significantly reduced, the simulation speed is fast, and the memory consumption is small.
The ESPM is a fast algorithm for solving the steady-state periodic solution of Zeta circuits.The algorithm is a symbolic algorithm which overcomes the drawback of too many variable symbols in similar symbolic algorithms and uses matrix operations to make the whole solution process intuitive and clear.However, the equivalent small parameter method also has some problems; for example, the modeling of a high-order converter is complicated, and the converter with low-order-of-energy storage components needs to derive the high-order correction.This paper uses ESPM to model and analyze the fractional-order Zeta converter in CCM mode, and deduces the boundary conditions of the CCM mode and DCM mode.The modeling method proposed in this paper can be easily extended and applied to the fractional-order Zeta converter in DCM mode and other converter circuits.

Figure 4 .
Figure 4. Simulation model of fractional-order CCM Zeta converter.

Figure 4 .
Figure 4. Simulation model of fractional-order CCM Zeta converter.

Figure 8 .
Figure 8.The relationship between inductance current ripple and inductance order.
interval, the current flowing through the diode is equal to the sum of the two inductor currents.

Figure 8 .
Figure 8.The relationship between inductance current ripple and inductance order.
interval, the current flowing through the diode is equal to the sum of the two inductor currents.

Figure 9 .
Figure 9. Influence of inductor order on operation mode of the converter.(a) Three-dimensional graph.(b) Boundary between CCM and DCM.Figure 9. Influence of inductor order on operation mode of the converter.(a) Three-dimensional graph.(b) Boundary between CCM and DCM.

Figure 9 .
Figure 9. Influence of inductor order on operation mode of the converter.(a) Three-dimensional graph.(b) Boundary between CCM and DCM.Figure 9. Influence of inductor order on operation mode of the converter.(a) Three-dimensional graph.(b) Boundary between CCM and DCM.

Figure 10 .
Figure 10.Harmonic components of inductor currents with different orders in the CCM region. 0

Figure 11 .
Figure 11.Harmonic components of capacitor voltages with different orders in the CCM region.

Figure 12 .
Figure 12.Change in effective value of harmonic component in converter state variable with order.

Figure 12 .
Figure 12.Change in effective value of harmonic component in converter state variable with order.

Table 2 .
Comparison of DC components of state variables.

Table 2 .
Comparison of DC components of state variables.

Table 3 .
Comparison of the theoretical value and the simulation value.

Table 3 .
Comparison of the theoretical value and the simulation value.
Harmonic components of inductor currents with different orders in the CCM region.theorder of the energy storage element in the study.The specific change rule is as follows: when the order of the energy storage element decreases, the ripple amplitude of the state variable increases; in the contrary, the ripple amplitude of the state variable decreases.When the inductor and capacitor are in different orders, the harmonic amplitudes of the converter inductor current and capacitor voltage are obtained, as shown inFigures 10 and  11, respectively.