Next Article in Journal
Study on Development Characteristics of Partial Discharge in Oil-Pressboard Insulation under Constant DC Voltage
Previous Article in Journal
A Review on Flame Stabilization Technologies for UAV Engine Micro-Meso Scale Combustors: Progress and Challenges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Harmonic Analysis of a Fractional-Order Zeta Converter

Guangxi Key Laboratory of Power System Optimization and Energy-Saving Technology, School of Electrical Engineering, Guangxi University, Nanning 530004, China
*
Author to whom correspondence should be addressed.
Energies 2023, 16(9), 3969; https://doi.org/10.3390/en16093969
Submission received: 16 March 2023 / Revised: 26 April 2023 / Accepted: 27 April 2023 / Published: 8 May 2023

Abstract

:
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 difficult 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 filter-based approximation method; the symbolic period results obtained by ESPM had sufficient precision because they included more combinations of higher harmonics. Finally, the influence 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.

1. 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 fractional-order calculus theory differ greatly from those of integer-order models. Compared with 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 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.

2. 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:
v L ( t ) = L d α i L d t α i C ( t ) = C d β v C d t β
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 ST and SD are controlled by the switching function δ(t). The fractional-order Zeta converter has two operating modes in one switching period TS. 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: n T 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, ST is turned on and SD is turned off, and the state equation is written:
d α 1 i L 1 d t α 1 = v i n L 1 d α 2 i L 2 d t α 2 = v i n v C 1 v C 2 L 2 d β 1 v C 1 d t β 1 = i L 2 C 1 d β 2 v c 2 d t β 2 = i L 2 C 2 v c 2 R C 2
Mode 2: ( n + D ) T S < t ( n + 1 ) T S , 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:
d α 1 i L 1 d t α 1 = v C 1 L 1 d α 2 i L 2 d t α 2 = v c 2 L 2 d β 1 v C 1 d t β 1 = i L 1 C 1 d β 2 v C 2 d t β 2 = i L 2 C 2 v c 2 R C 2
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:
δ ( t ) = 1 , n T S < t ( n + D ) T S 0 , ( n + D ) T S < t ( n + 1 ) T S
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
p α 1 i L 1 v C 1 L 1 ( 1 δ ( t ) ) V i n L 1 δ ( t ) = 0 P α 2 i L 2 + v C 1 L 2 δ ( t ) + v C 2 L 2 V i n L 2 δ ( t ) = 0 P β 1 v C 1 i L 2 C 1 δ ( t ) + i L 1 C 1 ( 1 δ ( t ) ) = 0 P β 2 v C 2 i L 2 C 2 + v C 2 R C 2 = 0
Representing Formula (5) in the form of a matrix equation, we have:
G 0 ( p α 1 , p α 2 , p β 1 , p β 2 ) x + G 1 f ( x ) G 2 t ( x ) = 0
There are nonlinear functions f(x) = δ(t)x, t(x) = δ(t)e, and e = V i n L 1 V i n L 2 0 0 T in Formula (6), and the coefficient matrix G0(pα1, pα2, pβ1, pβ2), G1, G2 are, respectively,
G 0 ( p α 1 , p α 2 , p β 1 , p β 2 ) = p α 1 0 1 L 1 0 0 p α 2 0 1 L 2 1 C 1 0 p β 1 0 0 1 C 2 0 p β 2 + 1 R C 2
G 1 = 0 0 1 L 1 0 0 0 1 L 2 0 1 C 1 1 C 1 0 0 0 0 0 0 , G 2 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1
The above is the state vector model of the fractional-order Zeta converter working in CCM. It can be observed that the coefficient matrix G0(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.

3. Equivalent Model of Converter Based on ESPM

3.1. 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:
d μ e ω τ d t μ = ω μ e ω τ
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 x0 and the other corrections xi:
x = x 0 + i = 1 n ε i x i
The switching function δ(t) is also expanded into a similar series sum form:
δ ( t ) = δ 0 + i = 1 n ε i δ i
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:
f ( x ) = f 0 + ε f 1 + ε 2 f 2 + t ( x ) = t 0 + ε t 1 + ε 2 t 2 +
where fi is:
f 0 = δ 0 x 0 f 1 = δ 0 x 1 + δ 1 x 0 f 2 = δ 0 x 2 + δ 1 x 2 + δ 2 x 0 f 3 = δ 0 x 3 + δ 1 x 2 + δ 2 x 1 + δ 3 x 0
ti is:
t 0 = δ 0 e t 1 = δ 1 e t 2 = δ 2 e t 3 = δ 3 e
According to the principle of harmonic balance, it can be assumed that xi can be expressed in the form of the sum of various harmonic components shown in Equation (15):
x i = k E i r x i k = a i 0 + ( a i k e j k τ + c . c )
where c.c represents the conjugate term, and k represents the order of the harmonic, in which τ = ωt is the normalized time. ai0 and aik represent the DC component of xi and the amplitude of the kth harmonic component, respectively. In addition, the formula’s harmonic components {Eir} 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 {E0} = {0} of the main oscillation component x0. The harmonic components {Eir} of the other corrections xi are sequentially determined by the remainder Ri 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:
δ ( t ) = b 0 + k = 1 ( b k e j k τ + c . c )
In the formula, b 0 = 1 T S 0 T S δ ( t ) d t = D , b k = 1 2 ( α k j β k ) , where:
α k = 2 T S 0 T S δ ( t ) cos ( k ω t ) d t = sin ( 2 k π D ) k π β k = 2 T S 0 T S δ ( t ) sin ( k ω t ) d t = 1 cos ( 2 k π D ) k π
Usually, δ0 and δi are taken as:
δ 0 = b 0 + b 1 e j τ + c . c δ i = b 2 i e j 2 i τ + b 2 i + 1 e j ( 2 i + 1 ) τ + c . c
Substitute the expansions of xi and δi into fi and ti. Due to the multiplication of xi and δi, new harmonic components are generated in fi and ti, as follows:
f i = f i m + φ R i + 1 t i = t i m + φ H i + 1
Among them, fim represents the part of fi that is the same as the harmonic component of xi, the remainder Ri+1 represents other components in fi, tim represents the part of ti that is the same as the harmonic component of xi, and the remainder Hi+1 represents other components in ti components, the spectrum of xi and Ri are the same as Hi. Therefore, the value of the xi spectrum can be determined by Ri and Hi in the operation.
Substituting Equations (15) and (19) into Equation (6), Equation (20) can be deduced as:
G 0 ( p α 1 , p α 2 , p β 1 , p β 2 ) x 0 + G 1 f 0 m - G 2 t 0 m = 0 G 0 ( p α 1 , p α 2 , p β 1 , p β 2 ) x 1 + G 1 ( f 1 m + R 1 ) - G 2 ( t 1 m + H 1 ) = 0 G 0 ( p α 1 , p α 2 , p β 1 , p β 2 ) x 2 + G 1 ( f 2 m + R 2 ) - G 2 ( t 2 m + H 2 ) = 0 G 0 ( p α 1 , p α 2 , p β 1 , p β 2 ) x 3 + G 1 ( f 3 m + R 3 ) - G 2 ( t 3 m + H 3 ) = 0
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.
T o l e r a n c e = a i k 2 a 00 2 1 %
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:
x = x 0 + x 1 + x 2 + x 3

3.2. Steady-State Analytical Solution of Fractional-Order CCM Zeta Converter

3.2.1. Solving the Main Oscillation Component x0

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:
x 0 = a 00 = I 00 , I 00 , V 00 , V 00 T
where I00, I′00, V00, and V′00 represent the DC components of the inductor current and capacitor voltage, respectively. From Equation (19) fi and ti expressions, it can be expressed as:
f 0 m = a 00 b 0 t 0 m = a 00 e
R 1 = a 00 ( b 1 e j τ + c . c ) H 1 = b 1 e e j τ + c . c
Substituting Equations (24) and (25) into the first term of Equation (20), it can be expressed as:
G 00 a 00 + G 1 D a 00 G 2 D e = 0
At this time, the coefficient matrix is G00 = G0(0,0,0,0). According to the principle of harmonic balance, the above formula can be transformed into
0 0 1 L 1 0 0 0 0 1 L 2 1 C 1 0 0 0 0 1 C 2 0 1 R C 2 + D 0 0 1 L 1 0 0 0 1 L 2 0 1 C 1 1 C 1 0 0 0 0 0 0 I 00 I 00 V 00 V 00 = D 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 V i n L 1 V i n L 2 0 0
The analytical formula of a00 can be obtained from Formula (27).

3.2.2. Solving the First-Order Correction Amount x1

The components of x1 depend on R1 and H1, so the set of harmonic components of the first-order correction amount x1 is {E1r} = {1}, and x1 can be expressed as
x 1 = a 11 e j τ + c . c
Where a11 = [I11, I’11, V11, V’11]T, bring x0, x1, δ0, δ1 into f1 and t1, and the remainder R2 of f1m and f1 can be obtained:
f 1 m = b 0 a 11 e j τ + c . c R 2 = b 1 a ¯ 11 + b ¯ 1 a 11 + ( b 1 a 11 + b 2 a 00 ) e j 2 τ + b 3 a 00 e j 3 τ + c . c
Since the spectrum of H1 is {K1r} = {1}, the spectrum of x1 is also {1}. In the same way, selecting the same amount as the spectrum of x1 in t1 as t1m, and the other components as H2, it can be expressed as:
t 1 m = 0 H 2 = b 2 e e j τ + b 3 e e j τ + c . c
Substituting f1m, t1m, H1, and R1 into the second equation of Equation (20), it can be reverted to:
( G 01 + G 1 b 0 ) a 11 = ( G 1 b 1 a 00 + G 2 b 1 e )
where G01 = G0[()α1, ()α2, ()β1, ()β2], a11 can be obtained by the matrix transformation of Formula (31).

3.2.3. Solve the Second-Order Correction Amount x2

According to R2 and H2, it is known that the set of harmonic components of the second-order correction amount x2 is {E2r} = {0, 2, 3}; then, the expression of the second-order correction amount is:
x 2 = a 20 + a 22 e j 2 τ + a 23 e j 3 τ + c . c
The DC component in the second-order correction is a20 = [I20, I’20, V20, V’20]T, the second harmonic is a22 = [I22, I’22, V22, V’22]T, and the third harmonic is a23 = [I23, I’23, V23, V’23]T. Substituting x0, x1, x2, δ0, δ1, and δ2 into f2 and t2, it can be expressed as:
f 2 m = b 0 a 20 + ( b 0 a 22 + b ¯ 1 a 23 + b 3 a ¯ 11 ) e j 2 τ + ( b 0 a 23 + b 1 a 22 + b 2 a 11 ) e j 3 τ + c . c R 3 = ( b 1 a 20 + b ¯ 1 a 22 + b 2 a ¯ 11 ) e j τ + ( b 1 a 23 + b 3 a 11 + b 4 a 00 ) e j 4 τ + b 5 a 00 e j 5 τ + c . c
t 2 m = 0 H 3 = b 4 e e j 4 τ + b 5 e e j 5 τ + c . c
Bringing x2, R2, f2m, H2, and t2m into the third formula of Formula (20), it can be expressed as:
( G 00 + G 1 b 0 ) a 20 = G 1 ( b 1 a ¯ 11 + b ¯ 1 a 11 ) ( G 02 + G 1 b 0 ) a 22 = G 1 ( b 1 a 11 + b 2 a 00 + b 3 a ¯ 11 ) + G 2 b 2 e ( G 03 + G 1 b 0 ) a 23 = G 1 ( b 1 a 22 + b 2 a 11 + b 3 a 00 ) + G 2 b 3 e
The matrix G0k = [(jkω)α1, (jkω)α2, (jkω)β1, (jkω)β2] (k = 0, 2, 3) in Equation (35), a20, a22 and a23 can be obtained by matrix transformation.

3.2.4. Solving the Third-Order Correction Amount x3

According to R3 and H3, it can be seen that the set of harmonic components of the third-order correction amount x3 is {E3r} = {1, 4, 5}; then, the expression of the third-order correction amount is:
x 3 = a 31 e j τ + a 34 e j 4 τ + a 35 e j 5 τ + c . c
The fundamental harmonic, fourth harmonic, and fifth harmonic in the third-order correction amount are, respectively, a31 = [I31, I’31, V31, V’31]T, a34 = [I34, I’34, V34, V’34]T, and a35 = [I35, I’35, V35, V’35]T. Substituting x0, x1, x2, x3, δ0, δ1, δ2, and δ3 into f3 and t3, it can be reverted to:
f 3 m = ( b 0 a 31 + b 3 a ¯ 22 ) e j τ + ( b 0 a 34 + b ¯ 1 a 35 + b 2 a 22 + b 5 a ¯ 11 ) e j 4 τ + ( b 0 a 35 + b 1 a 34 + b 2 a 23 + b 3 a 22 + b 4 a 11 ) e j 5 τ + c . c R 4 = b ¯ 1 a 31 + b 2 a ¯ 22 + b 3 a ¯ 23 + ( b 1 a 31 + b 2 a 20 ) e j 2 τ + ( b ¯ 1 a 34 + b 3 a 20 + b 4 a ¯ 11 ) e j 3 τ + c . c
t 3 m = 0
Bringing x3, R3, f3m, H3, and t3m into the fourth formula of Formula (20), it can be expressed as:
( G 01 + G 1 b 0 ) a 31 = G 1 ( b 1 a 20 + b 3 a ¯ 22 + b ¯ 1 a 22 + b 2 a ¯ 11 ) ( G 04 + G 1 b 0 ) a 34 = G 1 ( b 1 a 23 + b 2 a 22 + b 3 a 11 + b 4 a 00 + b 5 a ¯ 11 ) + G 2 b 4 e ( G 05 + G 1 b 0 ) a 35 = G 1 ( b 0 a 35 + b 1 a 34 + b 2 a 23 + b 3 a 22 + b 4 a 11 + b 5 a 00 ) + G 2 b 5 e
In the coefficient matrix G0k = G1[(jkω)α1, (jkω)α2, (jkω)β1, (jkω)β2] (k = 1, 4, 5) in Formula (39), a31, a34, and a35 can all be passed through matrix transformation.
According to the above process, aik = [Iik, I’ik, Vik, V’ik] can be obtained. According to the obtained main oscillation component x0, first-order correction amount x1, second-order correction amount x2 and third-order correction amount x3, 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 j 2 τ + a 23 e j 3 τ + a 34 e j 4 τ + a 35 e j 5 τ + c . c
The expressions of state variables iL1, iL2, vC1, and vC2 are:
i L 1 ( I 00 + I 20 ) + 2 ( Re I 11 cos ω t Im I 11 sin ω t + Re I 22 cos 2 ω t Im I 22 sin 2 ω t + Re I 23 cos 3 ω t Im I 23 sin 3 ω t + Re I 31 cos ω t Im I 31 sin ω t + Re I 34 cos 4 ω t Im I 34 sin 4 ω t + Re I 35 cos 5 ω t Im I 35 sin 5 ω t ) i L 2 ( I 00 + I 20 ) + 2 ( Re I 11 cos ω t Im I 11 sin ω t + Re I 22 cos 2 ω t Im I 22 sin 2 ω t + Re I 23 cos 3 ω t Im I 23 sin 3 ω t + Re I 31 cos ω t Im I 31 sin ω t + Re I 34 cos 4 ω t Im I 34 sin 4 ω t + Re I 35 cos 5 ω t Im I 35 sin 5 ω t ) v C 1 ( V 00 + V 20 ) + 2 ( Re V 11 cos ω t Im V 11 sin ω t + Re V 22 cos 2 ω t Im V 22 sin 2 ω t + Re V 23 cos 3 ω t Im V 23 sin 3 ω t + Re V 31 cos ω t Im V 31 sin ω t + Re V 34 cos 4 ω t Im V 34 sin 4 ω t + Re V 35 cos 5 ω t Im V 35 sin 5 ω t ) v C 2 ( V 00 + V 20 ) + 2 ( Re V 11 cos ω t Im V 11 sin ω t + Re V 22 cos 2 ω t Im V 22 sin 2 ω t + Re V 23 cos 3 ω t Im V 23 sin 3 ω t + Re V 31 cos ω t Im V 31 sin ω t + Re V 34 cos 4 ω t Im V 34 sin 4 ω t + Re V 35 cos 5 ω t Im V 35 sin 5 ω t )
where Re and Im represent the real and imaginary parts, respectively.

4. Simulation Comparison and Validation of Different Methods

4.1. 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,
Δ I L 1 = D V i n f s L 1 , L 1 = D V i n f s Δ I L 1
The output inductor current ripple can be expressed as,
Δ I L 2 = D V i n f s L 1 , L 2 = D V i n f s Δ I L 2
The capacitor ripple voltages ΔVC1 & ΔVC2 can be derived from the Kirchhoff’s current law for first and second mode as
Δ V C 1 = 1 C 1 0 T o n i C 1 d t , C 1 = D V 0 f s R L Δ V C 1
Δ V C 2 = 1 C 2 0 T / 2 Δ I L 2 4 d t , C 2 = D V i n 8 f s 2 L 2 Δ V C 2
where V0 = VC2.
The simulation parameters of Figure 1 are shown in Table 1.

4.2. 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 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.
It can be seen from Figure 6 and Figure 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 iL and vC 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 Figure 6 and Figure 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 Table 2, Table 3 and Table 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.

4.3. Fractional-Order Zeta Converter CCM Discriminant

From the method of reference [21], the expressions for the fractional order inductance currents ΔiL1 and ΔiL2 can be obtained as:
Δ i L 1 = V i n ( D T ) α 1 L 1 α 1 Γ ( α 1 ) Δ i L 2 = ( V i n V C 1 V C 2 ) ( D T ) α 2 L 2 α 2 Γ ( α 2 )
In Equation (46), Г(α) is the gamma function. Additionally, VC1 = −VC2, so by taking L = L1 = L2, α = α1 = α2, the fractional order inductor L1α1 and L2α2 current ripple calculation formula can be expressed uniformly as follows.
Δ i L = Δ i L 1 = Δ i L 2 = V i n ( D T ) α L α Γ ( α )
Equation (47) shows that when the fractional-order inductors L1α1 and L2α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 ΔiL with order α can be obtained by using Matlab software, as shown in Figure 8. As the order decreases, the inductor current ripple ΔiL 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 ( 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.
i S D = i L 1 + i L 2
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:
( I L 1 1 2 Δ i L 1 ) + ( I L 2 1 2 Δ i L 2 ) > 0
Substituting Equation (46) into Equation (49) and simplifying it can be expressed in the following form:
( I L 1 + I L 2 ) V i n 2 [ ( D T ) α 1 L 1 α 1 Γ ( α 1 ) + ( D T ) α 2 L 2 α 2 Γ ( α 2 ) ] > 0
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.

5. 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 Figure 10 and Figure 11, respectively.
From Figure 10 and Figure 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.

6. 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.

Author Contributions

Supervision, L.X.; writing—original draft preparation, D.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article. Further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Jonscher, A.K. Dielectric Relaxation in Solids, Chelsea; Dielectrics Pub.: London, UK, 1983. [Google Scholar]
  2. Westerlund, S. Causal Consulting: Dead Matter Has Memory; Casual Consulting: Kalmar, Sweden, 2002. [Google Scholar]
  3. Chen, X. Research on Analytical Modeling Methods and Nonlinear Dynamic Characteristics of Fractional-Order DC/DC Converters. Ph.D. Thesis, South China University of Technology, Guangzhou, China, 2018. [Google Scholar]
  4. Kianpoor, N.; Yousefi, M.; Bayati, N. Fractional Order Modelling of DC-DC Boost Converters. In Proceedings of the 2019 IEEE 28th International Symposium on Industrial Electronics (ISIE), Vancouver, BC, Canada, 12–14 June 2019. [Google Scholar]
  5. Tan, C.; Liang, Z.S. Fractional order modeling and analysis of Boost converter in inductor current pseudo-continuous mode. Acta Phys. Sin. 2017, 63, 58–67. [Google Scholar]
  6. Yang, R.C.; Liao, X.Z. Modeling and analysis of fractional order Buck converter using Caputo–Fabrizio derivative. Energy Rep. 2020, 6, 440–445. [Google Scholar] [CrossRef]
  7. Wang, F.; Ma, X.K. Modeling and analysis of the fractional order buck converter in DCM operation by using fractional calculus and the circuit-Averaging Technique. J. Power Electron. 2013, 13, 1008–1015. [Google Scholar] [CrossRef]
  8. Wei, Z.H.; Zhang, B.; Jiang, Y.W. Analysis and Modeling of Fractional-Order Buck Converter Based on Riemann-Liouville Derivative. IEEE Access 2019, 7, 162768–162777. [Google Scholar] [CrossRef]
  9. Xie, L.L.; Liu, Z.P.; Zhang, B. A Modeling and Analysis Method for CCM Fractional Order Buck-Boost Converter by Using R–L Fractional Definition. J. Electr. Eng. Technol. 2020, 15, 1651–1661. [Google Scholar] [CrossRef]
  10. Xie, L.L.; Liu, Z.P.; Ning, K.Z. Fractional-Order Adaptive Sliding Mode Control for Fractional-Order Buck-Boost Converters. J. Electr. Eng. Technol. 2022, 17, 1693–1704. [Google Scholar] [CrossRef]
  11. Sarwar, S. On the Existence and Stability of Variable Order Caputo Type Fractional Differential Equations. Fractal Fract. 2022, 6, 51. [Google Scholar] [CrossRef]
  12. Wang, Z.Q.; Cao, J.Y. Numerical Solution of Fractional-Order Differential Integration Equations and Their Error Analysis; Xi’an Jiaotong University Press: Xi’an, China, 2015. [Google Scholar]
  13. Chen, X.; Chen, Y.F.; Zhang, B. A Modeling and Analysis Method for Fractional-Order DC-DC Converters. IEEE Trans. Power Electron. 2017, 32, 7034–7044. [Google Scholar] [CrossRef]
  14. Chen, Y.F.; Chen, X.; Hu, J. A symbolic analysis method for fractional-order boost converter in discontinuous conduction mode. In Proceedings of the IECON 2017—43rd Annual Conference of the IEEE Industrial Electronics Society, Beijing, China, 29 October–1 November 2017; IEEE: New York, NY, USA, 2017. [Google Scholar]
  15. Xu, J.; Li, X.; Liu, H. Fractional-Order Modeling and Analysis of a Three-Phase Voltage Source PWM Rectifier. IEEE Access 2020, 8, 13507–13515. [Google Scholar] [CrossRef]
  16. Sharma, M.; Rajpurohit, B.S.; Agnihotri, S.; Rathore, A.K. Development of Fractional Order Modeling of Voltage Source Converters. IEEE Access 2020, 8, 142652–142662. [Google Scholar] [CrossRef]
  17. Chen, B.; Li, C.; Wilson, B.; Huang, Y. Fractional modeling and analysis of coupled MR damping system. IEEE/CAA J. Autom. Sin. 2016, 3, 288–294. [Google Scholar]
  18. Chen, Y.F.; Xiao, W.X.; Zhang, B. Nonlinear Modeling and Harmonic Analysis of Magnetic Resonant WPT System Based on Equivalent Small Parameter Method. IEEE Trans. Ind. Electron. 2019, 66, 6604–6612. [Google Scholar] [CrossRef]
  19. Tseng, C.-C.; Lee, S.-L. Computation of fractional derivatives using Fourier transform and digital FIR differentiator. Signal Process. 2000, 80, 151–159. [Google Scholar] [CrossRef]
  20. Mohamed, A.; Rachid, M. Numerical Simulations of Fractional Systems: An Overview of Existing Methods and Improvements. Nonlinear Dyn. 2004, 38, 117–131. [Google Scholar]
  21. Xie, L.L.; Ning, K.Z.; Qin, L. Fractional Order Modeling and Simulation Analysis of Cuk Converter. Comput. Simul. 2022, 39, 313–320. [Google Scholar]
Figure 1. Fractional-order Zeta converter.
Figure 1. Fractional-order Zeta converter.
Energies 16 03969 g001
Figure 2. Working mode 1 of fractional-order CCM-Zeta converter.
Figure 2. Working mode 1 of fractional-order CCM-Zeta converter.
Energies 16 03969 g002
Figure 3. Working mode 2 of fractional-order CCM-Zeta converter.
Figure 3. Working mode 2 of fractional-order CCM-Zeta converter.
Energies 16 03969 g003
Figure 4. Simulation model of fractional-order CCM Zeta converter.
Figure 4. Simulation model of fractional-order CCM Zeta converter.
Energies 16 03969 g004
Figure 5. Equivalent circuit of fractional element. (a) Fractional-order inductor; (b) fractional-order capacitor.
Figure 5. Equivalent circuit of fractional element. (a) Fractional-order inductor; (b) fractional-order capacitor.
Energies 16 03969 g005
Figure 6. Inductance current waveform of different orders: (a) (α1, α2, β1, β2) = (1, 1, 1, 1); (b) (α1, α2, β1, β2) = (0.95, 0.95, 1, 1); (c) (α1, α2, β1, β2) = (0.95, 0.95, 0.95, 0.95); (d) (α1, α2, β1, β2) = (0.9, 0.9, 0.95, 0.95); (e) (α1, α2, β1, β2) = (0.9, 0.9, 0.9, 0.9); (f) (α1, α2, β1, β2) = (0.85, 0.85, 0.85, 0.85).
Figure 6. Inductance current waveform of different orders: (a) (α1, α2, β1, β2) = (1, 1, 1, 1); (b) (α1, α2, β1, β2) = (0.95, 0.95, 1, 1); (c) (α1, α2, β1, β2) = (0.95, 0.95, 0.95, 0.95); (d) (α1, α2, β1, β2) = (0.9, 0.9, 0.95, 0.95); (e) (α1, α2, β1, β2) = (0.9, 0.9, 0.9, 0.9); (f) (α1, α2, β1, β2) = (0.85, 0.85, 0.85, 0.85).
Energies 16 03969 g006
Figure 7. Capacitance voltage waveform of different orders: (a) (α1, α2, β1, β2) = (1, 1, 1, 1); (b) (α1, α2, β1, β2) = (0.95, 0.95, 1, 1); (c) (α1, α2, β1, β2) = (0.95, 0.95, 0.95, 0.95); (d) (α1, α2, β1, β2) = (0.9, 0.9, 0.95, 0.95); (e) (α1, α2, β1, β2) = (0.9, 0.9, 0.9, 0.9); (f) (α1, α2, β1, β2) = (0.85, 0.85, 0.85, 0.85).
Figure 7. Capacitance voltage waveform of different orders: (a) (α1, α2, β1, β2) = (1, 1, 1, 1); (b) (α1, α2, β1, β2) = (0.95, 0.95, 1, 1); (c) (α1, α2, β1, β2) = (0.95, 0.95, 0.95, 0.95); (d) (α1, α2, β1, β2) = (0.9, 0.9, 0.95, 0.95); (e) (α1, α2, β1, β2) = (0.9, 0.9, 0.9, 0.9); (f) (α1, α2, β1, β2) = (0.85, 0.85, 0.85, 0.85).
Energies 16 03969 g007aEnergies 16 03969 g007b
Figure 8. The relationship between inductance current ripple and inductance order.
Figure 8. The relationship between inductance current ripple and inductance order.
Energies 16 03969 g008
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.
Energies 16 03969 g009
Figure 10. Harmonic components of inductor currents with different orders in the CCM region.
Figure 10. Harmonic components of inductor currents with different orders in the CCM region.
Energies 16 03969 g010
Figure 11. Harmonic components of capacitor voltages with different orders in the CCM region.
Figure 11. Harmonic components of capacitor voltages with different orders in the CCM region.
Energies 16 03969 g011
Figure 12. Change in effective value of harmonic component in converter state variable with order.
Figure 12. Change in effective value of harmonic component in converter state variable with order.
Energies 16 03969 g012
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParametersValues
Vin12 V
R10 Ω
D0.4
fs25,000 Hz
L12 mH
L22 mH
C110−5 F
C210−5 F
Table 2. Comparison of DC components of state variables.
Table 2. Comparison of DC components of state variables.
Order (α1, α2, β1, β2)ESPM (a00 + a20)Oustaloup’s Method
(0.85, 0.85, 0.85, 0.85)(0.5574, 0.7515, −7.5145, 7.5145)(0.6040, 0.7536, −7.5270, 7.530)
(0.9, 0.9, 0.9, 0.9)(0.5417, 0.7836, −7.8359, 7.8359)(0.5676, 0.7826, −7.8250, 7.8230)
(0.9, 0.9, 0.95, 0.95)(0.5494, 0.7948, −7.9478, 7.9478)(0.5813, 0.7954, −7.950, 7.950)
(0.95, 0.95, 0.95, 0.95)(0.5354, 0.7955, −7.9546, 7.9546)(0.5596, 0.7956, −7.9530, 7.9530)
(0.95, 0.95, 1, 1)(0.5383, 0.7995, −7.9955, 7.9955)(0.5625, 0.7998, −7.9940, 7.9950)
(1, 1, 1, 1)(0.5330, 0.7998, −7.9975, 7.9975)(0.5497, 0.8009, −7.990, 7.9970)
Table 3. Comparison of the theoretical value and the simulation value.
Table 3. Comparison of the theoretical value and the simulation value.
Order (α1, α2, β1, β2)ΔiLTheoretical ValueESPMOustaloup
(1, 1, 1, 1)ΔiL1/A0.0960.0870.09501
ΔiL2/A0.0960.08830.0949
(0.95, 0.95, 1, 1)ΔiL1/A0.17020.18120.1858
ΔiL2/A0.17020.17230.1851
(0.95, 0.95, 0.95, 0.95)ΔiL1/A0.17020.18530.1874
ΔiL2/A0.17020.16780.1839
(0.9, 0.9, 0.95, 0.95)ΔiL1/A0.30120.3130.3435
ΔiL2/A0.30120.29680.3358
(0.9, 0.9, 0.9, 0.9)ΔiL1/A0.30120.31320.3435
ΔiL2/A0.30120.29180.3201
(0.85, 0.85, 0.85, 0.85)ΔiL1/A0.5320.57160.5825
ΔiL2/A0.5320.52460.5637
Table 4. Ripple error percentage.
Table 4. Ripple error percentage.
Order (α1, α2, β1, β2)ΔiLESPMOustaloup
(1, 1, 1, 1)ΔiL1/A9.37%1.03%
ΔiL2/A8.02%1.14%
(0.95, 0.95, 1, 1)ΔiL1/A6.46%9.16%
ΔiL2/A1.23%8.75%
(0.95, 0.95, 0.95, 0.95)ΔiL1/A8.87%10.1%
ΔiL2/A1.41%8.05%
(0.9, 0.9, 0.95, 0.95)ΔiL1/A3.91%14.04%
ΔiL2/A1.46%11.48%
(0.9, 0.9, 0.9, 0.9)ΔiL1/A3.98%14.04%
ΔiL2/A3.12%6.27%
(0.85, 0.85, 0.85, 0.85)ΔiL1/A7.44%9.49%
ΔiL2/A1.39%5.95%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xie, L.; Wan, D. Modeling and Harmonic Analysis of a Fractional-Order Zeta Converter. Energies 2023, 16, 3969. https://doi.org/10.3390/en16093969

AMA Style

Xie L, Wan D. Modeling and Harmonic Analysis of a Fractional-Order Zeta Converter. Energies. 2023; 16(9):3969. https://doi.org/10.3390/en16093969

Chicago/Turabian Style

Xie, Lingling, and Di Wan. 2023. "Modeling and Harmonic Analysis of a Fractional-Order Zeta Converter" Energies 16, no. 9: 3969. https://doi.org/10.3390/en16093969

APA Style

Xie, L., & Wan, D. (2023). Modeling and Harmonic Analysis of a Fractional-Order Zeta Converter. Energies, 16(9), 3969. https://doi.org/10.3390/en16093969

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop