Next Article in Journal
On the Impacts of the Global Sea Level Dynamics
Next Article in Special Issue
Application of Tilt Integral Derivative for Efficient Speed Control and Operation of BLDC Motor Drive for Electric Vehicles
Previous Article in Journal
Research on Application of Fractional Calculus Operator in Image Underlying Processing
Previous Article in Special Issue
Adaptive Residual Useful Life Prediction for the Insulated-Gate Bipolar Transistors with Pulse-Width Modulation Based on Multiple Modes and Transfer Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional-Order LCL Filters: Principle, Frequency Characteristics, and Their Analysis

1
College of Electrical Engineering, Guangxi University, Nanning 530004, China
2
College of Mechanical and Electrical and Quality Technology Engineering, Nanning University, Nanning 530200, China
3
Electric Power Research Institute of Guangxi Power Grid Corporation, Nanning 530023, China
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2024, 8(1), 38; https://doi.org/10.3390/fractalfract8010038
Submission received: 3 December 2023 / Revised: 26 December 2023 / Accepted: 2 January 2024 / Published: 5 January 2024

Abstract

:
The fractional-order LCL filter, composed of two fractional-order inductors and one fractional-order capacitor, is a novel fractional-order π-type circuit introduced in recent years. Based on mathematical modeling, this article comprehensively studies the principles and frequency characteristics of fractional-order LCL filters. Five critical properties are derived and rigorously demonstrated. One of the most significant findings is that we identify the necessary and sufficient condition for resonance in fractional-order LCL filters when the sum of the orders of the fractional-order inductors and the fractional-order capacitor is equal to 2, which provides a theoretical foundation for effectively avoiding resonance in fractional-order LCL filters. The correctness of our theoretical derivation and analysis was confirmed through digital simulations. This study reveals that fractional-order LCL filters exhibit more versatile operational characteristics than traditional integer-order LCL filters, paving the way for broader application prospects.

1. Introduction

The development of fractional calculus and the research on the implementation methods of fractional-order elements have opened up new approaches for the design, modeling, and analysis of circuits [1,2,3,4,5,6,7,8,9,10,11]. Fractional-order circuits have become an emerging interdisciplinary research area and have drawn much interest [5].
Fractional-order inductors and fractional-order capacitors [12,13] are the essential building elements of fractional-order circuits. The mathematical models of a fractional-order inductor and a fractional-order capacitor can be expressed as [6,7,8].
{ L d α i L d t α = u L , C d β u C d t β = i C ,
where L is the inductance of a fractional-order inductor, and C is the capacitance of a fractional-order inductor. α and β are the order of a fractional-order inductor and a fractional-order capacitor, respectively. α and β satisfy 0 < α , β < 2 . Particularly, when α = 1 and β = 1 , a fractional-order inductor and a fractional-order capacitor are traditional inductors and capacitors, respectively. u L is the inductor voltage, i L is the inductor current, u C is the capacitor voltage, and i C is the capacitor current.
From (1), we know that in addition to the traditional parameter (such as L or C), the fractional-order element introduces an additional parameter, namely order, which makes the element and system more complex and more difficult to analyze, on the one hand, but, on the other hand, it increases the flexibility and degree of freedom and brings in immense versatility towards design and application. For more information about the operating characteristics of fractional-order inductors and fractional-order capacitors, please refer to [3,6,7,8,12,13].
In recent years, many realizations and fabrication methods of fractional-order inductors and fractional-order capacitors have been reported [14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30], and the orders of some capacitors, such as super-capacitors, were verified as possibly being much further away from 1 by experiments [27]. At the same time, fractional-order inductors and fractional-order capacitors have gradually been applied to various traditional integer-order circuits, and counterpart fractional-order circuits are constructed, such as fractional-order chaotic circuits [3], fractional-order sinusoidal oscillators [4], fractional-order RLC/RC/RL circuits [6,7,10,11], fractional-order power electronic converters [8,9,31,32,33,34,35,36], etc.
The traditional integer-order LCL filter [37,38,39], which consists of two inductors, L1 and L2, and one capacitor, C, is a pivotal circuit extensively utilized in power electronic converters, such as in a grid-connected inverter, as shown in Figure 1a.
Figure 1b displays the Bode plot of an integer-order LCL filter, where ωrp denotes resonant frequency. When ω < ωrp and ω > ωrp, the asymptotic slopes of the logarithmic magnitude-frequency characteristic curve of an integer-order LCL filter are −20 dB/dec and −60 dB/dec, respectively [39]. The integer-order LCL filter plays a key role in filtering high-order harmonics in the AC side voltage and improving the quality of the current injected into the grid [37,38,39]. Compared with L filters, integer-order LCL filters can achieve a better high-frequency harmonic attenuation effect. However, the integer-order LCL filter has an inherent resonant peak, which will destroy the stability of the system when the system is not designed and controlled properly [37,38,39].
A fractional-order capacitor is introduced to replace the integer-order capacitor in an integer-order LCL filter to form an LCαL filter [40]. Theoretical analysis and numerical simulations show that the LCαL filter provides a wider bandwidth to mitigate higher-order resonant frequencies than its integer-order counterpart. The circuit and mathematical model of a fractional-order LCL filter are preliminarily given in [41]. However, the detailed theoretical derivation and analysis of the frequency characteristics are missing.
This study employs Caputo’s fractional calculus [1] to investigate the mathematical modeling and frequency characteristic analysis of fractional-order LCL filters. The aim is to expand the LCL filters from the integer-order domain to the fractional-order domain, laying a theoretical groundwork for the engineering application of fractional-order LCL filters. The main contributions and originality of this research are summarized as follows:
(1)
This paper pioneered a method for the theoretical analysis of the frequency characteristics of fractional-order LCL filters, summarized their five critical properties, and systematically revealed their principles and frequency characteristics.
(2)
It is found that the necessary and sufficient condition for resonance in the magnitude-frequency characteristic curve of fractional-order LCL filters is that the sum of the orders of the fractional-order inductor and the fractional-order capacitor is equal to 2. This provides a theoretical basis for effectively avoiding resonance in fractional-order LCL filters.
(3)
This paper fills the gap in the research on the frequency characteristics of general fractional systems with 2 α + β -order (where α , β ( 0 , 2 ) ).
The rest of this paper is organized as follows: the circuit and mathematical models of a fractional-order LCL filter are given in Section 2. In Section 3, the frequency characteristics of a fractional-order LCL filter are studied in depth, and five important properties of a fractional-order LCL filter are systematically presented. In Section 4, digital simulations are performed to verify the correctness of the theoretical derivation and analysis. Finally, conclusions are drawn in Section 5.

2. The Circuit and Mathematical Models of a Fractional-Order LCL Filter

By introducing the fractional-order inductors and the fractional-order capacitor to replace the inductors and the capacitor in an integer-order LCL filter, respectively, the main circuit of a fractional-order LCL filter can be obtained, as shown in Figure 2a, where the fractional-order LCL filter is used for the AC side filtering of a grid-connected inverter. L 1 denotes the inductance of the input side fractional-order inductor in units of H / sec ( 1 α 1 ) , and its order is α 1 ; L 2 denotes the inductance of the output side fractional-order inductor in units of H / sec ( 1 α 2 ) , and its order is α 2 ; C denotes the capacitance of the fractional-order capacitor in units of F / sec ( 1 β ) , and its order is β . All the value ranges of α 1 , α 2 , and β are (0, 2). Obviously, a fractional-order LCL filter is the general form of an LCL filter, and an integer-order LCL filter is just the special case of a fractional-order LCL filter when α 1 = α 2 = β = 1 . In addition, u dc is the DC voltage provided by the renewable energy generation or energy storage unit; u i , u g , and u C are the input voltage, the output voltage, and the fractional-order capacitor’s terminal voltage of a fractional-order LCL filter, respectively. i 1 , i 2 , and i C are the input current, the output current, and the fractional-order capacitor’s current of a fractional-order LCL filter, respectively. When a fractional-order LCL filter is used for the AC side filtering of a grid-connected inverter, u i , u g , and i 2 are also known as the AC side voltage of the inverter, the power grid voltage, and the current injected into the grid, respectively.
According to the Kirchhoff voltage law and the Kirchhoff current law, the time domain mathematical model of a fractional-order LCL filter can be expressed as
{ L 1 d α 1 i 1 d t α 1 = u i u C , C d β u C d t β = i 1 i 2 , L 2 d α 2 i 2 d t α 2 = u C u g .
By performing the Laplace transform on (2), the s-domain expression of a fractional-order LCL filter can be obtained as
{ L 1 s α 1 i 1 ( s ) = u i ( s ) u C ( s ) , C s β u C ( s ) = i 1 ( s ) i 2 ( s ) , L 2 s α 2 i 2 ( s ) = u C ( s ) u g ( s ) .
According to (3), we obtain the model structure of a fractional-order LCL filter, as shown in Figure 2b.
According to (3) and Figure 2b, the output current expression of a fractional-order LCL filter can be obtained as
i 2 ( s ) = 1 L 1 L 2 C s α 1 + α 2 + β + L 1 s α 1 + L 2 s α 2 u i ( s ) L 1 C s α 1 + β + 1 L 1 L 2 C s α 1 + α 2 + β + L 1 s α 1 + L 2 s α 2 u g ( s ) .
In order to simplify the analysis, let α 1 = α 2 = α ; then, (4) can be simplified as
i 2 ( s ) = 1 L 1 L 2 C s 2 α + β + ( L 1 + L 2 ) s α u i ( s ) L 1 C s α + β + 1 L 1 L 2 C s 2 α + β + ( L 1 + L 2 ) s α u g ( s ) .
Generally, u g ( s ) is regarded as the disturbance of the system, so the transfer function from the input voltage u i ( s ) to the output current i 2 ( s ) can be obtained; that is, the transfer function of a fractional-order LCL filter is
G gi ( s ) = i 2 ( s ) u i ( s ) = 1 L 1 L 2 C s α 1 s α + β + A .
where
A = ( L 1 + L 2 ) / ( L 1 L 2 C ) ,   α + β ( 0 , 4 ) .
From (6), we know that a fractional-order LCL filter is a fractional-order system with order 2 α + β , where α , β ( 0 , 2 ) . Until now, no study has systematically introduced the frequency characteristics of such fractional-order systems.
According to (6), the frequency domain model of a fractional-order LCL filter can be expressed as
G gi ( j ω ) = 1 L 1 L 2 C ( j ω ) α 1 ( j ω ) α + β + A .
Substituting the Euler equation
( j ω ) α = ω α e j α π / 2 = ω α cos α π 2 + j ω α sin α π 2
into (8), the frequency domain expression of a fractional-order LCL filter can be further expressed as
G gi ( j ω ) = 1 L 1 L 2 C ω α ( cos α π 2 + j sin α π 2 ) 1 ω α + β cos ( α + β ) π 2 + A + j ω α + β sin ( α + β ) π 2 .
Based on (9), the expressions of the magnitude-frequency characteristics and the phase-frequency characteristics of a fractional-order LCL filter can be obtained as
| G gi ( j ω ) | = 1 / ( L 1 L 2 C ω α ) ( ω α + β + A cos ( α + β ) π 2 ) 2 + A 2 sin 2 ( α + β ) π 2 .
G gi ( j ω ) = arctan ( tan π α 2 ) arctan ω α + β sin ( α + β ) π 2 ω α + β cos ( α + β ) π 2 + A .

3. The Frequency Characteristics and Analysis of a Fractional-Order LCL Filter

3.1. The Resonance Characteristics of a Fractional-Order LCL Filter

In (10), when α + β ( 0 , 1 ] [ 3 , 4 ) , since cos [ ( α + β ) π / 2 ] 0 , with the angular frequency ω increasing from 0, the denominator of | G gi ( j ω ) | increases monotonically; that is, | G gi ( j ω ) | decreases monotonically. In this case, there is no resonance in the magnitude-frequency characteristics of a fractional-order LCL filter.
When α + β ( 1 , 3 ) , cos [ ( α + β ) π / 2 ] < 0 , let us define an intermediate variable
ω r = [ A cos ( α + β ) π / 2 ] 1 α + β ,   α + β ( 1 , 3 ) ,
namely, ω r α + β = A cos [ ( α + β ) π / 2 ] .
Substituting the expression of ω r into (10), we obtain
| G gi ( j ω ) | = 1 / ( L 1 L 2 C ω α ) ( ω α + β ω r α + β ) 2 + A 2 sin 2 ( α + β ) π 2 . α + β ( 1 , 3 )
When the angular frequency satisfies
ω = ω r = [ A cos ( α + β ) π 2 ] 1 α + β ,   α + β ( 1 , 3 ) ,
| G gi ( j ω ) | can be expressed as
| G gi ( j ω r ) | = 1 / ( L 1 L 2 C ω r α ) A | sin [ ( α + β ) π / 2 ] | ,   α + β ( 1 , 3 ) .
According to (15), we know that when α + β = 2 , due to A sin [ ( α + β ) π / 2 ] = 0 , we have | G gi ( j ω r ) | = . This means that the magnitude-frequency characteristics of a fractional-order LCL filter exhibit a resonance peak at ω = ω r . Conversely, if the magnitude-frequency characteristics of a fractional-order LCL filter have a resonance peak when ω = ω r , that is, | G gi ( j ω r ) | = , it follows that α + β = 2 must be satisfied. For ease of distinction, we denote the ω r when α + β = 2 as ω rp , which represents the resonant frequency.
From the above analysis, we discern that the fractional-order LCL filter possesses the following resonance property:
Property 1. 
When α 1 = α 2 = α , the necessary and sufficient condition for a fractional-order LCL filter to have a resonance peak is
α + β = 2 .
Substituting α + β = 2 into (12), the expression of the resonant frequency ω rp can be obtained as follows:
ω rp = A = ( L 1 + L 2 ) / ( L 1 L 2 C ) .
According to (17), we further obtain the following resonant frequency property of a fractional-order LCL filter:
Property 2. 
When the magnitude-frequency characteristics of a fractional-order LCL filter have a resonant peak, that is, the condition of α + β = 2 holds, the expression of the resonant frequency is ω r p = ( L 1 + L 2 ) / ( L 1 L 2 C ) . The resonant frequency ω r p is only determined by the values of L 1 , L 2 , and C , and is independent of the order α of the fractional-order inductors and the order β of the fractional-order capacitor.
Property 1 reveals the condition for the existence of the resonance of a fractional-order LCL filter, essentially, and provides a straightforward criterion for determining whether a fractional-order LCL filter exhibits a resonant peak. For example, in traditional integer-order LCL filters, since both α and β are 1, that is, α + β = 2 , their magnitude-frequency characteristics invariably exhibit a resonance peak, as illustrated in Figure 1b. Furthermore, Property 1 offers theoretical ground for effectively circumventing resonance in the design phase of a fractional-order LCL filter. When a fractional-order LCL filter manifests a resonance peak, it can readily lead to system instability [37,38,39]. For an integer-order LCL filter, additional passive or active dampers are commonly incorporated to eliminate resonance and enhance system stability against variations [37,38]. However, the former can incur additional losses, diminishing the attenuation capability of high-frequency harmonics, while the latter can complicate the design of the controller. In contrast, for a fractional-order LCL, we can effectively circumvent the resonance peak by judiciously selecting the order of fractional-order inductors and the fractional-order capacitor, such that α + β 2 , thereby fundamentally enhancing system stability without adding dampers.

3.2. The Corner Frequency and Logarithmic Magnitude-Frequency Characteristics of a Fractional-Order LCL Filter

To analyze the corner frequency of a fractional-order LCL filter, according to the value range of α + β , the following two intermediate variables of angular frequency are defined.
ω t 1 = | A cos ( α + β ) π 2 | 1 α + β ,   α + β ( 0 , 0.5 ] [ 1.5 , 2.5 ] [ 3.5 , 4 ) .
ω t 2 = | A sin ( α + β ) π 2 | 1 α + β ,   α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 ) .
Figure 3 shows the waveforms of | cos ( α + β ) π / 2 | 1 / ( α + β ) in (18) and | sin ( α + β ) π / 2 | 1 / ( α + β ) in (19). From Figure 3, we know that when α + β ( 0 , 0.5 ] [ 1.5 , 2.5 ] [ 3.5 , 4 ) , there is   ω t 1 > ω t 2 ; when α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 ) , there is   ω t 1 < ω t 2 .
Substituting (18) and (19) into (10), the magnitude-frequency characteristics of a fractional-order LCL filter can be expressed as
| G gi ( j ω ) | = 1 / ( L 1 L 2 C ω α ) ( ω α + β + σ ω t 1 α + β ) 2 + ω t 2 2 ( α + β ) ,
where
σ = { 1 , α + β ( 0 , 1 ] [ 3 , 4 ) , 1 , α + β ( 1 , 3 ) .
(1) The analysis when α + β ( 0 , 0.5 ] [ 1.5 , 2.5 ] [ 3.5 , 4 )
According to the value range of ω , (20) can be further expressed as
| G gi ( j ω ) | = { 1 / ( L 1 L 2 C ω t 1 α + β ω α ) [ ( ω / ω t 1 ) α + β + σ ] 2 + ( ω t 2 / ω t 1 ) 2 ( α + β ) ,   ω < < ω t 1 , 1 / [ L 1 L 2 C ω 2 α + β ] [ 1 + σ ( ω t 1 / ω ) α + β ] 2 + ( ω t 2 / ω ) 2 ( α + β ) ,   ω > > ω t 1 .
For the first formula of (22), when ω < < ω t 1 , ( ω / ω t 1 ) α + β 0 , considering that σ 2 = 1 always holds, so it can be simplified as
| G gi ( j ω ) | 1 / ( L 1 L 2 C ω α ) ω t 1 2 ( α + β ) + ω t 2 2 ( α + β ) = 1 L 1 L 2 C A ω α ,   ω < < ω t 1 .
According to (23), the logarithmic magnitude-frequency characteristic and its asymptotic slope of a fractional-order LCL filter when ω < < ω t 1 can be obtained as
L ( ω ) 20 lg ( L 1 L 2 C A ) 20 α lg ω ,   ω < < ω t 1 .
d L ( ω ) d lg ω 20 α dB / dec ,   ω < < ω t 1 .
For the second formula of (22), when ω > > ω t 1 , ( ω t 1 / ω ) α + β 0 . Additionally, considering that there is   ω t 1 > ω t 2 when α + β ( 0 , 0.5 ] [ 1.5 , 2.5 ] [ 3.5 , 4 ) , then ( ω t 2 / ω ) α + β 0 . So, it can be simplified as
| G gi ( j ω ) | 1 L 1 L 2 C ω 2 α + β ,   ω > > ω t 1 .
According to (26), the logarithmic magnitude-frequency characteristic and its asymptotic slope of a fractional-order LCL filter when ω > > ω t 1 can be obtained as
L ( ω ) = 20 lg ( L 1 L 2 C ) 20 ( 2 α + β ) lg ω ,   ω > > ω t 1 .
d L ( ω ) d lg ω 20 ( 2 α + β ) dB / dec ,   ω > > ω t 1 .
(2) The analysis when α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 )
When α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 ) , according to the value range of angular frequency, (20) can be further expressed as
| G gi ( j ω ) | = { 1 / ( L 1 L 2 C ω α ω t 2 α + β ) [ ( ω / ω t 2 ) α + β + σ ( ω t 1 / ω t 2 ) α + β ] 2 + 1 ,   ω < < ω t 2 , 1 / [ L 1 L 2 C ω 2 α + β ] [ 1 + σ ( ω t 1 / ω ) α + β ] 2 + ( ω t 2 / ω ) 2 ( α + β ) ,   ω > > ω t 2 .
For the first formula of (29), when ω < < ω t 2 , ( ω / ω t 2 ) α + β 0 , considering that σ 2 = 1 always holds, it can be simplified as
| G gi ( j ω ) | 1 / ( L 1 L 2 C ω α ) ω t 1 2 ( α + β ) + ω t 2 2 ( α + β ) = 1 L 1 L 2 C A ω α ,   ω < < ω t 2 .
According to (30), we obtain
L ( ω ) 20 lg ( L 1 L 2 C A ) 20 α lg ω ,   ω < < ω t 2 ,
d L ( ω ) d lg ω 20 α dB / dec ,   ω < < ω t 2 ,
For the second formula of (29), when ω > > ω t 2 , ( ω t 2 / ω ) α + β 0 . Additionally, considering that there is ω t 2 > ω t 1 when α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 ) , then ( ω t 1 / ω ) α + β 0 . So, it can be simplified as
| G gi ( j ω ) | 1 L 1 L 2 C ω 2 α + β ,   ω > > ω t 2 .
According to (33), we obtain
L ( ω ) 20 lg ( L 1 L 2 C ) 20 ( 2 α + β ) lg ω ,   ω > > ω t 2 ,
d L ( ω ) d lg ω 20 ( 2 α + β ) dB / dec ,   ω > > ω t 2 .
(3) A summary of the corner frequency and logarithmic magnitude-frequency characteristic
Drawing from the aforementioned theoretical derivation, we can deduce that when α + β ( 0 , 0.5 ] [ 1.5 , 2.5 ] [ 3.5 , 4 ) , the corner frequency of the fractional-order LCL filter is ω t 1 , and when α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 ) , the corner frequency is ω t 2 . ω t 1 and ω t 2 can be condensed into the following expression for the corner frequency.
ω t = { ω t 1 = | A cos ( α + β ) π 2 | 1 α + β , α + β ( 0 , 0.5 ] [ 1.5 , 2.5 ] [ 3.5 , 4 ) , ω t 2 = | A sin ( α + β ) π 2 | 1 α + β , α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 ) .
Furthermore, to facilitate discussions, we define the frequency ranges as the low-frequency band for ω < < ω t and the high-frequency band for ω > > ω t . By comparing (30)–(35) with (23)–(28), we observe that the expressions for the magnitude-frequency characteristics of the fractional-order LCL filter are identical when α + β (0, 0.5] ∪ [1.5, 2.5] ∪ [3.5, 4) and α + β (0.5, 1.5) ∪ (2.5, 3.5). The primary distinction lies in their corner frequency expressions. In summary, the logarithmic magnitude-frequency characteristics of fractional-order LCL filters can be summarized as follows.
Property 3. 
When ω < < ω t , the asymptotic slope of the logarithmic magnitude-frequency characteristic in the low-frequency band of a fractional-order LCL filter is 20   α d B / d e c ; when ω > > ω t , the asymptotic slope of the logarithmic magnitude-frequency characteristic in the high-frequency band is 20 ( 2 α + β ) d B / d e c .
Property 3 is crucial for analyzing the filtering performance of fractional-order LCL filters. From Property 3, it can be observed that the asymptotic slope of the logarithmic magnitude-frequency characteristics of a fractional-order LCL filter in the low-frequency band is solely determined by the order α of the fractional-order inductors and is independent of the order β of the fractional-order capacitor; in the high-frequency band, the asymptotic slope is affected by both α and β . Considering that α , β ( 0 , 2 ) , we know that the asymptotic slope range of the logarithmic magnitude-frequency characteristics of a fractional-order LCL filter in the low-frequency band and high-frequency band is (0 dB/dec, −40 dB/dec) and (0 dB/dec, −120 dB/dec), respectively.
According to (36), we find that the corner frequency of a fractional-order LCL filter is not only affected by the values of L 1 , L 2 and C but also related to the order α of the fractional-order inductors and the order β of the fractional-order capacitor. Additionally, when α + β = 2 , there is ω rp = ω t 1 , which means that the resonant frequency is equal to the corner frequency. In particular, when α = β = 1 , according to (25), (28), and (36), we obtain the corner frequency of an integer-order LCL filter as ω t = A = ( L 1 + L 2 ) / ( L 1 L 2 C ) , and the asymptotic slopes of the low-frequency and high-frequency bands are −20 dB/dec and −60 dB/dec, respectively.

3.3. The Phase-Frequency Characteristics of a Fractional-Order LCL Filter

It is known to all that for a complex number z = x + j y , ( z 0 ) , the value range of π < arg z π is generally called the principal argument of z , and there is
arg z = { arctan ( y / x ) ,   z   is   in   the   first   and   fourth   quadrants , π + arctan ( y / x ) ,   z   is   in   the   second   quadrant ,   π + arctan ( y / x ) ,   z   is   in   the   third   quadrant ,
where π / 2 < arctan ( y / x ) < π / 2 .
From (37), the first term of (11) can be further expressed as
arctan [ tan ( π α / 2 ) ] = { π α / 2 , α ( 0 , 1 ] , ( π α / 2 π ) + π = π α / 2 , α ( 1 , 2 ) ,
which means arctan [ tan ( π α / 2 ) ] = π α / 2 always holds when α ( 0 , 2 ) .
To facilitate the derivation, let the center frequency
ω o = A 1 / ( α + β ) = ( L 1 + L 2 ) / ( L 1 L 2 C ) α + β ,
then, the (11) can be further simplified as
G gi ( j ω ) = π α 2 arctan sin ( α + β ) π 2 cos ( α + β ) π 2 + ( ω o / ω ) α + β .
When ω < < ω o , considering the characteristics of the sin and cos functions, the following relationship always holds.
sin ( α + β ) π 2 < < cos ( α + β ) π 2 + ( ω o / ω ) α + β ,
then, there is
arctan sin ( α + β ) π 2 cos ( α + β ) π 2 + ( ω o / ω ) α + β 0 ,   ω < < ω o .
Therefore, the phase-frequency characteristics of a fractional-order LCL filter when ω < < ω o are
G gi ( j ω ) π α / 2 ,   ω < < ω o .
When ω > > ω o , ω o / ω 0 , so (40) can be further simplified as
G gi ( j ω ) π α 2 arctan [ tan ( α + β ) π 2 ] , ω > > ω o .
According to (37), the second term in (42) can be further expressed as
arctan [ tan ( α + β ) π 2 ] = { ( α + β ) π 2 ,   α + β ( 0 , 2 ] , ( α + β ) π 2 2 π ,   α + β ( 2 , 4 ) .
By synthesizing (42) and (43), the phase-frequency characteristics of a fractional-order LCL filter when ω > > ω o can be obtained as
G gi ( j ω ) { π ( α + β / 2 ) ,   α + β ( 0 , 2 ] , ω > > ω o , π ( α + β / 2 ) + 2 π ,   α + β ( 2 , 4 ) , ω > > ω o .
Based on (42) and (43), the phase-frequency characteristics of a fractional-order LCL filter can be summarized as follows.
Property 4. 
When ω < < ω o , the phase of a fractional-order LCL filter is π α / 2 ; when ω > > ω o , if α + β ( 0 , 2 ] , the phase is ( α + β / 2 ] π . If α + β ( 2 , 4 ) , the phase is ( α + β / 2 ) π + 2 π .
According to Property 4, it can be observed that when ω < < ω o , the phase of the fractional-order LCL filter is solely determined by the order α of the fractional-order inductor, and it is independent of the order β of the fractional-order capacitor. Conversely, when ω > > ω o , the phase is influenced by both α and β .

3.4. The Phase Crossover Frequency and Gain Margin of a Fractional-Order LCL Filter

It is well known that the frequency at which the open-loop phase-frequency characteristic intersects the horizontal line π is commonly referred to as the phase crossover frequency, denoted as ω g . According to (40), we know that the intersection frequency ω g should satisfy the following equation for a fractional-order LCL filter.
G gi ( j ω g ) = π α 2 arctan sin ( α + β ) π 2 cos ( α + β ) π 2 + ( ω o / ω g ) α + β = π .
By solving (45), we obtain
ω g α + β = A sin [ ( 1 α / 2 ) π ] sin [ ( α + β / 2 ) π ] .
Given that 0 < α < 2 , we have 0 < ( 1 α / 2 ) π < π , so the numerator of (46) satisfies A sin [ ( 1 α / 2 ) π ] > 0 . Considering the domains of α and β , the denominator of (46) satisfies sin [ ( α + β / 2 ) π ] > 0 , if and only if α + β / 2 ( 1 , 2 ) . This means the right side of (46) is greater than 0, indicating that (46) has a real number solution, implying the presence of a phase crossover frequency in the phase-frequency characteristic of a fractional-order LCL filter. To sum up, the phase crossover frequency characteristics of a fractional-order LCL filter can be summarized as follows.
Property 5. 
If and only if α + β / 2 ( 1 , 2 ) , the phase-frequency characteristics of a fractional-order LCL filter has a phase crossover frequency ω g , and the expression of ω g is
ω g = A s i n [ ( 1 α / 2 ) π ] s i n [ ( α + β / 2 ) π ] α + β .
In particular, when α = β = 1 , there is ω g = A = ( L 1 + L 2 ) / ( L 1 L 2 C ) = ω rp ; that is, the phase crossover frequency is equal to the resonant frequency.
When α + β / 2 ( 1 , 2 ) , the gain margin of a fractional-order LCL filter and its logarithmic expression can be obtained as
K g = 1 | G gi ( j ω g ) | = L 1 L 2 C ω g α ω g 2 ( α + β ) + 2 A ω g α + β cos ( α + β ) π 2 + A 2 .
20 K g = 20 lg | G gi ( j ω g ) | = 20 lg 1 / L 1 L 2 C ω g α ω g 2 ( α + β ) + 2 A ω g α + β cos ( α + β ) π 2 + A 2 .

3.5. The Gain Crossover Frequency and Phase Margin of a Fractional-Order LCL Filter

As is widely known, the frequency ω c satisfies | G gi ( j ω c ) | = 1 or 20 lg | G gi ( j ω c ) | = 0 is usually referred to as the phase crossover frequency. According to (10), we know that the gain crossover frequency ω c of a fractional-order LCL filter should satisfy the following equation.
1 / ( L 1 L 2 C ω c α ) ( ω c α + β + A cos ( α + β ) π 2 ) 2 + A 2 sin 2 ( α + β ) π 2 = 1 .
By sorting out Equation (50), we obtain
ω c 2 ( 2 α + β ) + 2 A ω c 3 α + β cos ( α + β ) π 2 + A 2 ω c 2 α A 2 ( L 1 + L 2 ) 2 = 0 .
The above formula is a highly complex nonlinear equation, making it challenging to derive an analytical solution for the gain crossover frequency ω c of a fractional-order LCL filter. In order to obtain the trend graph of the change in ω c with the order of the fractional-order inductors and the order of the fractional-order capacitor, we successively used the fsolve function, the particle swarm optimization algorithm [42], and the differential evolution algorithm [43] to solve (49) based on MATLAB 2018a software, and we found that the differential evolution algorithm has the best solution effect and relatively stable results.
Figure 4 shows a change diagram of the digital solution of ω c obtained by the differential evolution algorithm versus α β plane when α , β ( 0.8 , 1.2 ) . The main circuit parameters of a fractional-order LCL filter are as follows: L 1 = 600   μ H / sec ( 1 α ) , L 2 = 150   μ H / sec ( 1 α ) , C = 10   μ F / sec ( 1 β ) .
From Figure 4, we observe that as the combination of α and β approaches α + β = 2 , there is a notable increase in the value of ω c . This is attributed to the presence of a resonance peak in the fractional-order LCL filter when α + β = 2 , which makes the logarithmic magnitude-frequency characteristic curve cross the 0 dB line many times. This phenomenon will be further verified in the following digital simulations.
Additionally, according to Property 3 and Figure 4, we know the following:
(1) Except for the case of α + β = 2 , as the order α of the fractional-order inductors increases, the asymptotic slope of the logarithmic magnitude-frequency characteristic curve in the low-frequency band increases accordingly, and the gain crossover frequency decreases.
(2) Considering that the order β of the fractional-order capacitor does not affect the low-frequency asymptotic slope of the logarithmic magnitude-frequency characteristic curve of a fractional-order LCL filter, altering the order β does not alter the gain crossover frequency, except in the case of α + β = 2 .
The formula for calculating the phase margin of a fractional-order LCL filter is the same as that of an integer-order LCL filter, namely
γ = 180 ° + G gi ( j ω c ) .

4. Simulation Results and Analysis

To validate the theoretical derivation, analyze the frequency characteristics of fractional-order LCL filters, and investigate the impact of the order of the fractional-order inductors and the fractional-order capacitor on these characteristics, a digital model of a fractional-order LCL filter is established and simulated using MATLAB software. The main circuit parameters of the fractional-order LCL filter are consistent with those in Figure 4, and the order of the fractional-order inductors and the fractional-order capacitor will be specified later.

4.1. The Frequency Characteristic Simulation Curves of a Fractional-Order LCL Filter

(1) The frequency characteristic curves when α is constant and β changes
Simulation condition I: Set the order α of the fractional-order inductors to 0.8, 1.0, and 1.2, respectively. For each case, increase the order β of the fractional-order capacitor from 0.6 to 1.4 in 0.2 intervals.
The Bode plots of a fractional-order LCL filter based on Simulation condition I, as shown in Figure 5, and the corresponding frequency characteristic indicators are summarized in Table 1. In this table, NaN indicates the absence of such indicators, and Inf denotes infinite indicators.
(2) The frequency characteristic curves when β is constant and α changes
Simulation condition II: Set the order β of the fractional-order capacitor to 0.8, 1.0, and 1.2, respectively. For each case, increase the order α of the fractional-order inductor from 0.6 to 1.4 in 0.2 intervals.
The Bode plots of a fractional-order LCL filter, as shown in Figure 6, and the corresponding frequency characteristic indicators are summarized in Table 2.

4.2. The Analysis of the Frequency Characteristics of a Fractional-Order LCL Filter

From Figure 5 and Figure 6, as well as Table 1 and Table 2, we can observe that:
(1) When α + β = 2 , for instance, the combinations of α and β are (0.8, 1.2), (1.0, 1.0), and (1.2, 0.8), the magnitude-frequency characteristic curves of these fractional-order LCL filters cross the 0 dB line multiple times and exhibit a resonance peak at the resonant frequency. Moreover, the phase-frequency characteristic curves make a jump of −180° at the resonant frequency. For example, when the combination of α and β is (0.8, 1.2), the phase-frequency characteristic curve makes a jump from −72° to −252°. Furthermore, it can be observed that the resonant frequency ω rp is solely determined by L 1 , L 2 and C and is independent of α and β . In this scenario, the theoretically calculated resonant frequency is
ω rp = ( L 1 + L 2 ) / ( L 1 L 2 C ) 28,867.5 rad / s ,
and the measured resonant frequency from the curves is 28,867 rad/s. Obviously, these two frequencies are essentially identical, with the minor deviation attributable to the discrete points connected by the curve in the Bode plot. In conclusion, the correctness of Properties 1 and 2 is confirmed.
(2) In each subplot of Figure 5, when the order α of the fractional-order inductors remains constant and the order β of the fractional-order capacitor changes, the five logarithmic magnitude-frequency characteristic curves overlap in the low-frequency band and diverge in the high-frequency band. This suggests that the order β of the fractional-order capacitor does not affect the logarithmic magnitude-frequency characteristics of a fractional-order LCL filter in the low-frequency band, but it does have an effect in the high-frequency band. The asymptotic slope of the logarithmic magnitude-frequency characteristic curve of a fractional-order LCL filter in the low-frequency band is 20 α dB / dec , which is solely related to α . For example, when α is 0.8, 1.0, and 1.2, the theoretical low-frequency slopes calculated based on Property 3 are 16 dB / dec , 20 dB / dec , and 24 dB / dec , while the measured values from the curves are 15.96 dB / dec , 19.9 dB / dec , and 23.87 dB / dec , aligning with the theory.
Additionally, the asymptotic slope of the logarithmic magnitude-frequency characteristic curve in the high-frequency band is 20 ( 2 α + β ) dB / dec , which is influenced by α and β . For example, when the combinations of α and β are ( 0.8 , 1.4 ) , ( 1.0 , 1.4 ) , and ( 1.2 , 1.4 ) , the theoretical high-frequency slopes are 60 dB / dec , 68 dB / dec , and 76 dB / dec , with the measured curve values being 60.07 dB / dec , 68.01 dB / dec , and 76.00 dB / dec , which demonstrates that the theoretical calculated values are consistent with the actual measured values, too.
In each subplot of Figure 6, when β remains constant, and α takes on different values, the slopes of the low-frequency and high-frequency bands of the five logarithmic magnitude-frequency characteristic curves differ, further indicating that the asymptotic slope of the logarithmic magnitude-frequency characteristics of a fractional-order LCL filter in the low-frequency band are solely determined by α , while the asymptotic slope in the high-frequency band is influenced by α and β .
Furthermore, when α + β ( 0 , 0.5 ] [ 1.5 , 2.5 ] [ 3.5 , 4 ) , the corner frequency ω t is calculated according to the first formula of (36); for example, when α = β = 1.0 , the theoretically calculated and the measured corner frequencies ω t are 28,868 rad/s and 28,867 rad/s, respectively, which are essentially equivalent. When α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 ) , the corner frequency should be calculated according to the second formula of (36); for example, when the combinations of α and β are ( 0.8 , 0.6 ) and ( 1.2 , 1.4 ) , the theoretically calculated corner frequencies ω t are 2,024,272 rad/s and 2487 rad/s. The corner frequencies derived from the curve are 2,024,102 rad/s and 2487 rad/s, respectively. This suggests that the theoretically calculated values align with the actual measurements.
To sum up, the correctness of Property 3 has been confirmed.
(3) In each subplot of Figure 5, when α is constant and β changes, if ω < < ω o , the five phase-frequency characteristic curves overlap; if ω > > ω o , they diverge. This indicates that β solely influences the phase-frequency characteristics of a fractional-order LCL filter in the high-frequency band but not in the low-frequency band. The low-frequency phase follows the π α / 2 law. For example, when α is 0.8, 1.0, and 1.2, the theoretically calculated and measured low-frequency phases are 72 ° , 90 ° , and 108 ° , respectively. It can be observed that the low-frequency phase is solely determined by α .
Additionally, when α + β ( 0 , 2 ] , the high-frequency phase follows the rule ( α + β / 2 ] π . For example, when the combinations of α and β are ( 0.8 , 0.6 ) , ( 1.0 , 0.6 ) , and ( 1.2 , 0.6 ) , both the theoretically calculated and the measured high-frequency phases equal 198 ° , 234 ° , and 270 ° , respectively. When α + β ( 2 , 4 ) , the high-frequency phase follows the rule ( α + β / 2 ) π + 2 π . For example, when the combinations of α and β are ( 0.8 , 1.4 ) , ( 1.0 , 1.4 ) , and ( 1.2 , 1.4 ) , the theoretically calculated and measured high-frequency phases match 90 ° , 54 ° , and 18 ° , respectively. Overall, the correctness of Property 4 has been verified.
(4) The measured value of the phase crossover frequencies ω g shown in Table 1 and Table 2 equals the theoretical value calculated by (47). For example, when the combinations of α and β are ( 0.8 , 0.8 ) , ( 0.8 , 1.0 ) , ( 1.0 , 1.0 ) , and ( 0.6 , 1.2 ) , the calculated and measured phase-crossing frequencies are equivalent, specifically 508,310 rad/s, 98,864 rad/s, 28,867 rad/s, and 107,916 rad/s, respectively. In conclusion, the correctness of Property 5 is confirmed.
(5) From Table 1 and Figure 5, it is evident that when α remains constant, and β varies, the gain crossover frequency ω c remains relatively constant, except when α + β = 2 . For example, when the combinations of α and β are ( 1.0 , 0.8 ) and ( 1.0 , 1.2 ) , the gain crossover frequencies are nearly identical, specifically 1334 rad/s and 1345 rad/s, respectively. Furthermore, when α + β approaches 2, there is a resonance in the magnitude-frequency characteristics of a fractional-order LCL filter, resulting in a sudden increase in the gain crossover frequency. Similarly, Table 2 and Figure 6 show that when β remains constant and α varies, the gain crossover frequency gradually increases as α decreases, except when α + β = 2 . These simulation results generally align with the observations presented in Figure 5.
Overall, the theoretical derivation and analysis of the frequency characteristics of fractional-order LCL filters have been validated through digital simulations.

5. Conclusions

Building on the mathematical model of fractional-order LCL filters, this article deduces five critical properties of these filters, systematically analyzing and summarizing their working principles and frequency characteristics. The primary differences in frequency characteristics between integer-order LCL filters and fractional-order LCL filters are summarized in Table 3.
Compared with an integer-order LCL filter, since the extra parameters of α and β are included, although a fractional-order LCL filter is more complicated to analyze, it has many novel behaviors and a wider operating range, which brings great flexibility and advantages to the application. The research in this paper provides a theoretical basis for the design of fractional-order LCL filters. By appropriately selecting the values of α and β , it is possible to obtain fractional-order LCL filters that can better coordinate filtering performance and system stability. It is expected that fractional-order LCL filters will be applied to replace traditional integer-order LCL filters in power electronic systems, such as inverters and rectifiers, and to construct better-performing fractional-order power electronic systems.
Furthermore, the frequency characteristic analysis method of fractional-order LCL filters proposed in this paper fills the gap in the research on the frequency characteristics of general fractional systems with 2 α + β -order (where α , β ( 0 , 2 ) ) and has good theoretical significance.

Author Contributions

Conceptualization, J.X.; Methodology, J.X., E.Z. and X.L.; Software, E.Z., G.H. and X.M.; Validation, X.L. and X.M.; Formal analysis, W.L.; Investigation, J.X. and G.H.; Resources, G.H. and W.L.; Data curation, J.X. and E.Z.; Writing—original draft, J.X.; Writing—review & editing, X.L.; Supervision, X.L. 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 authors.

Conflicts of Interest

Author Xuanren Meng was employed by the company Electric Power Research Institute of Guangxi Power Grid Corporation. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Podlubny, I. Fractional Differential Equations; Academic: New York, NY, USA, 1999. [Google Scholar]
  2. Monje, C.A.; Chen, Y.; Vinagre, B.M.; Xue, D.; Feliu, V. Fractional-Order Systems and Controls: Fundamentals and Applications; Springer: London, UK, 2010. [Google Scholar]
  3. Petras, I. Fractional Calculus. In Fractional-Order Nonlinear Systems; Higher Education: Beijing, China, 2011. [Google Scholar]
  4. Radwan, A.G.; Soliman, A.M.; Elwakil, A.S. Design equations for fractional-order sinusoidal oscillators: Four practical circuit examples. Int. J. Circuit Theory Appl. 2008, 36, 473–492. [Google Scholar] [CrossRef]
  5. Elwakil, A.S. Fractional-order circuits and systems: An emerging interdisciplinary research area. IEEE Circuits Syst. Mag. 2010, 10, 40–50. [Google Scholar] [CrossRef]
  6. Radwan, A.G.; Salama, K.N. Fractional-order RC and RL circuits. Circuits Syst. Signal Process. 2012, 31, 1901–1915. [Google Scholar] [CrossRef]
  7. Radwan, A.G.; Fouda, M.E. Optimization of Fractional-Order RLC Filters. Circuits Syst. Signal Process. 2013, 32, 2097–2118. [Google Scholar] [CrossRef]
  8. Xu, J.; Li, X.; Meng, X.; Qin, J.; Liu, H. Modeling and analysis of a single-phase fractional-order voltage source pulse width modulation rectifier. J. Power Sources 2020, 479, 228821. [Google Scholar] [CrossRef]
  9. Xu, J.; Li, X.; Liu, H.; Meng, X. Fractional-order modeling and analysis of a Three-phase Voltage Source PWM Rectifier. IEEE Access 2020, 8, 13507–13515. [Google Scholar] [CrossRef]
  10. Freeborn, T.J.; Maundy, B.; Elwakil, A. Fractional resonance-based RLβCα filters. Math. Probl. Eng. 2013, 2013, 726721. [Google Scholar] [CrossRef]
  11. Diao, L.; Zhang, X.; Chen, D. Fractional-order multiple RL alpha C beta circuit. Acta Phys. Sin. 2014, 63, 038401. [Google Scholar] [CrossRef]
  12. Westerlund, S.; Ekstam, L. Capacitor theory. IEEE Trans. Dielectr. Elect. Insul. 1994, 1, 826–839. [Google Scholar] [CrossRef]
  13. Westerlund, S. Dead Matter Has Memory; Causal Consulting: Kalmar, Sweden, 2002. [Google Scholar]
  14. Valsa, J.; Vlach, J. RC models of a constant phase element. Int. J. Circuit Theory Appl. 2013, 41, 59–67. [Google Scholar] [CrossRef]
  15. Sarafraz, M.S.; Tavazoei, M.S. Passive realization of fractional-order impedances by a fractional element and RLC components: Conditions and procedure. IEEE Trans. Circuits Syst. I Reg. Pap. 2017, 64, 585–595. [Google Scholar] [CrossRef]
  16. Semary, M.S.; Fouda, M.E.; Hassan, H.N.; Radwan, A.G. Realization of fractional-order capacitor based on passive symmetric network. J. Adv. Res. 2019, 18, 147–159. [Google Scholar] [CrossRef]
  17. Tsirimokou, G.; Psychalinos, C.; Elwakil, A.S. Emulation of a constant phase element using operational transconductance amplifiers. Analog Integr. Circuits Signal Process. 2015, 85, 413–423. [Google Scholar] [CrossRef]
  18. Bertsias, P.; Psychalinos, C.; Radwan, A.G.; Elwakil, A.S. High-Frequency Capacitorless Fractional-Order CPE and FI Emulator. Circuits Syst. Signal Process. 2018, 37, 2694–2713. [Google Scholar] [CrossRef]
  19. Jiang, Y.; Zhang, B. High-Power Fractional-Order Capacitor With 1 < alpha < 2 Based on Power Converter. IEEE Trans. Ind. Electron. 2018, 65, 3157–3164. [Google Scholar]
  20. Kapoulea, S.; Psychalinos, C.; Elwakil, A.S.; Radwan, A.G. One-terminal electronically controlled fractional-order capacitor and inductor emulator. AEU-Int. J. Electron. Commun. 2019, 103, 32–45. [Google Scholar] [CrossRef]
  21. Jesus, I.S.; Machado, J.A.T. Development of fractional order capacitors based on electrolyte processes. Nonlinear Dyn. 2009, 56, 45–55. [Google Scholar] [CrossRef]
  22. Krishna, M.S.; Das, S.; Biswas, K.; Goswami, B. Fabrication of a fractional order capacitor with desired specifications: A study on process identification and characterization. IEEE Trans. Electron Devices 2011, 58, 4067–4073. [Google Scholar] [CrossRef]
  23. Haba, T.C.; Ablart, G.; Camps, T.; Olivie, F. Influence of the electrical parameters on the input impedance of a fractal structure realised on silicon. Chaos Solitons Fractals 2005, 24, 479–490. [Google Scholar] [CrossRef]
  24. Mondal, D.; Biswas, K. Packaging of single-component fractional order element. IEEE Trans. Device Mater. Reliab. 2013, 13, 73–80. [Google Scholar] [CrossRef]
  25. Agambayev, A.; Patole, S.P.; Farhat, M.; Elwakil, A.; Bagci, H.; Salama, K.N. Ferroelectric fractional-order capacitors. ChemElectroChem 2017, 4, 2807–2813. [Google Scholar] [CrossRef]
  26. Agambayev, A.; Patole, S.; Bagci, H.; Salama, K.N. Tunable fractional-order capacitor using layered ferroelectric polymers. AIP Adv. 2017, 7, 095202. [Google Scholar] [CrossRef]
  27. Freeborn, T.J.; Maundy, B.; Elwakil, A.S. Measurement of supercapacitor fractional-order model parameters from voltage-excited step response. IEEE J. Emerg. Sel. Top. Circuits Syst. 2013, 3, 367–376. [Google Scholar] [CrossRef]
  28. Machado, J.A.T.; Galhano, A.M.S.F. Fractional order inductive phenomena based on the skin effect. Nonlinear Dyn. 2012, 68, 107–115. [Google Scholar] [CrossRef]
  29. Tripathy, M.C.; Mondal, D.; Biswas, K.; Sen, S. Experimental studies on realization of fractional inductors and fractional-order bandpass filters. Int. J. Circuit Theory Appl. 2015, 43, 1183–1196. [Google Scholar] [CrossRef]
  30. Adhikary, A.; Choudhary, S.; Sen, S. Optimal design for realizing a grounded fractional order inductor using GIC. IEEE Trans. Circuits Syst. I Reg. Pap. 2018, 65, 2411–2421. [Google Scholar] [CrossRef]
  31. Wang, F.; Ma, X. Transfer function modeling and analysis of the open-loop buck converter using the fractional calculus. Chin. Phys. B 2013, 22, 030506. [Google Scholar] [CrossRef]
  32. Wang, F.; Ma, X. 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]
  33. Wei, Z.; Zhang, B.; Jiang, Y. Analysis and Modeling of Fractional-Order Buck Converter Based on Riemann-Liouville Derivative. IEEE Access 2019, 7, 162768–162777. [Google Scholar] [CrossRef]
  34. Tan, C.; Liang, Z. Modeling and simulation analysis of fractional-order Boost converter in pseudo-continuous conduction mode. Acta Phys. Sin. 2014, 63, 070502. [Google Scholar]
  35. Yang, N.; Liu, C.; Wu, C. Modeling and dynamics analysis of the fractional-order Buck-Boost converter in continuous conduction mode. Chin. Phys. B 2012, 21, 080503. [Google Scholar] [CrossRef]
  36. Chen, X.; Chen, Y.; Zhang, B.; Qiu, D. A Modeling and Analysis Method for Fractional-Order DC-DC Converters. IEEE Trans. Power Electron. 2017, 32, 7034–7044. [Google Scholar] [CrossRef]
  37. Twining, E.; Holmes, D.G. Grid current regulation of a three-phase voltage source inverter with an LCL input filter. IEEE Trans. Power Electron. 2003, 18, 888–895. [Google Scholar] [CrossRef]
  38. Wang, X.; Ruan, X.; Liu, S.; Tse, C.K. Full feedforward of grid voltage for grid-connected inverter with LCL filter to suppress current distortion due to grid voltage harmonics. IEEE Trans. Power Electron. 2010, 25, 3119–3127. [Google Scholar] [CrossRef]
  39. Wu, W.; He, Y.; Tang, T.; Blaabjerg, F. A New Design Method for the Passive Damped LCL and LLCL Filter-Based Single-Phase Grid-Tied Inverter. IEEE Trans. Ind. Electron. 2013, 60, 4339–4350. [Google Scholar] [CrossRef]
  40. El-Khazali, R. Fractional-Order LCαL Filter-Based Grid Connected PV Systems. In Proceedings of the 2019 IEEE 62nd International Midwest Symposium on Circuits and Systems (MWSCAS), Dallas, TX, USA, 4–7 August 2019; pp. 533–536. [Google Scholar]
  41. Sun, E. Modeling and Simulation of Fractional Active Power Filter. Master’s Thesis, School of Electrical Engineering, Dalian University of Technology, Dalian, China, 2017. (In Chinese with English abstract). [Google Scholar]
  42. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the IEEE International Conference on Neural Networks, Perth, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  43. Price, K.V.; Storn, R. Differential evolution: A simple evolution strategy for fast optimization. Dr. Dobb’s J. 1997, 22, 18–24. [Google Scholar]
Figure 1. The main circuit topology and Bode plot of an integer-order LCL filter. (a) A grid-connected inverter with an integer-order LCL filter. (b) The Bode plot of an integer-order LCL filter.
Figure 1. The main circuit topology and Bode plot of an integer-order LCL filter. (a) A grid-connected inverter with an integer-order LCL filter. (b) The Bode plot of an integer-order LCL filter.
Fractalfract 08 00038 g001
Figure 2. The main circuit topology and the model structure of a fractional-order LCL filter. (a) The main circuit topology of a fractional-order LCL filter. (b) The model structure of a fractional-order LCL filter.
Figure 2. The main circuit topology and the model structure of a fractional-order LCL filter. (a) The main circuit topology of a fractional-order LCL filter. (b) The model structure of a fractional-order LCL filter.
Fractalfract 08 00038 g002
Figure 3. The waveforms of | cos ( α + β ) π / 2 | 1 / ( α + β ) and | sin ( α + β ) π / 2 | 1 / ( α + β ) versus α + β .
Figure 3. The waveforms of | cos ( α + β ) π / 2 | 1 / ( α + β ) and | sin ( α + β ) π / 2 | 1 / ( α + β ) versus α + β .
Fractalfract 08 00038 g003
Figure 4. The change diagram of the digital solution of ω c obtained by the differential evolution algorithm versus α β plane.
Figure 4. The change diagram of the digital solution of ω c obtained by the differential evolution algorithm versus α β plane.
Fractalfract 08 00038 g004
Figure 5. The frequency characteristic curves of a fractional-order LCL Filter when α is constant and β changes. (a) The frequency characteristics when α = 0.8 . (b) The frequency characteristics when α = 1.0 . (c) The frequency characteristics when α = 1.2 .
Figure 5. The frequency characteristic curves of a fractional-order LCL Filter when α is constant and β changes. (a) The frequency characteristics when α = 0.8 . (b) The frequency characteristics when α = 1.0 . (c) The frequency characteristics when α = 1.2 .
Fractalfract 08 00038 g005
Figure 6. The frequency characteristic curves of a fractional-order LCL Filter when β is constant and α changes. (a) The frequency characteristics when β = 0.8 . (b) The frequency characteristics when β = 1.0 . (c) The frequency characteristics when β = 1.2 .
Figure 6. The frequency characteristic curves of a fractional-order LCL Filter when β is constant and α changes. (a) The frequency characteristics when β = 0.8 . (b) The frequency characteristics when β = 1.0 . (c) The frequency characteristics when β = 1.2 .
Fractalfract 08 00038 g006
Table 1. The frequency characteristic indicators of a fractional-order LCL when α is constant and β changes.
Table 1. The frequency characteristic indicators of a fractional-order LCL when α is constant and β changes.
Serial Number α β ω g (rad/s) 20 lg K g (dB) ω c (rad/s) γ (°) ω t (rad/s)
10.80.65,257,08353.398059107.982,024,102
20.80.8508,31028.808075107.93329,599
30.81.098,8649.068186107.7687,883
40.81.228,867−79.0333,194−72.0028,867
50.81.4NaNInf13,749241.7011,092
61.00.6429,57447.38133490.00329,598
71.00.892,92227.10133489.9987,883
81.01.028,867−75.5829,51090.0028,867
91.01.2NaNInf134590.1611,092
101.01.4NaNInf137991.434772
111.20.687,88345.9640271.9987,882
121.20.828,867−57.7328,950−108.028,867
131.21.0NaNInf40272.0111,092
141.21.2NaNInf40272.074772
151.21.4NaNInf40372.332487
Table 2. The frequency characteristic indicators of a fractional-order LCL when β is constant and α changes.
Table 2. The frequency characteristic indicators of a fractional-order LCL when β is constant and α changes.
Serial Number α β ω g (rad/s) 20 lg K g (dB) ω c (rad/s) γ (°) ω t (rad/s)
10.60.819,437,84650.5116,538124.852,024,630
20.80.8508,31028.808075107.93329,599
31.00.892,92227.10133489.9987,883
41.20.828,867−57.7328,950−108.028,867
51.40.8NaNInf17154.0011,092
60.61.0686,67013.13368,96256.77329,598
70.81.098,8649.068186107.7687,882
81.01.028,867−75.5829,510−90.0028,867
91.21.0NaNInf40272.0111,092
101.41.0NaNInf17154.014772
110.61.2107,916−7.69133,810−19.9887,883
120.81.228,867−93.4233,192−72.0028,867
131.01.2NaNInf134590.1611,092
141.21.2NaNInf40272.074771
151.41.2NaNInf170.754.042487
Table 3. The comparison between an integer-order LCL filter and a fractional-order LCL filter.
Table 3. The comparison between an integer-order LCL filter and a fractional-order LCL filter.
PropertiesInteger-Order LCL FiltersFractional-Order LCL FiltersNotes on Fractional-Order LCL Filters
VariablesThree Variables ( L 1 , L 2 , C ) Five Variables ( L 1 , L 2 , C , α , β ) α is the order of the fractional-order inductors, and β is the order of the fractional-order capacitor.
Range of α and β α = β = 1 α , β ( 0 , 2 ) An integer-order LCL filter is the special case of a fractional-order LCL filter when α = β = 1 .
The transfer function, G gi ( s ) 1 L 1 L 2 C s 1 s 2 + A 1 L 1 L 2 C s α 1 s α + β + A
Resonance peakExists a resonance peakExists a resonance peak when α + β = 2 The necessary and sufficient condition for the existence of a resonance peak is α + β = 2 .
Resonant frequency, ω rp ( L 1 + L 2 ) / ( L 1 L 2 C ) ( L 1 + L 2 ) / ( L 1 L 2 C ) ω rp is determined by the values of L 1 , L 2 , and C , and is independent of α and β .
Corner frequency, ω t ( L 1 + L 2 ) / ( L 1 L 2 C ) { ω t 1 = | A cos [ ( α + β ) π / 2 ] | α + β α + β ( 0 , 0.5 ] [ 1.5 , 2.5 ] [ 3.5 , 4 ) { ω t 2 = | A sin [ ( α + β ) π / 2 ] | α + β α + β ( 0.5 , 1.5 ) ( 2.5 , 3.5 ) According to the different value range of α + β , there are two calculation formulas of ω t . ω t is affected by both α and β .
Slope of the logarithmic magnitude-frequency characteristic, d L ( ω ) / d lg ω { 20 dB / dec ,   ω < < ω t 60 dB / dec ,   ω > > ω t { 20 α dB / dec ,   ω < < L ω t 20 ( 2 α + β ) dB / dec ,   ω > > ω t The slope is only determined by α when ω < < ω t , while is affected by both α and β when ω > > ω t . The range of slope is (0 dB/dec, −40 dB/dec) and (0 dB/dec, −120 dB/dec) when ω < < ω t and ω > > ω t , respectively.
Center frequency, ω o ( L 1 + L 2 ) / ( L 1 L 2 C ) ( L 1 + L 2 ) / ( L 1 L 2 C ) α + β ω o is affected by both α and β .
Phase-frequency characteristic, G gi ( j ω ) { π / 2 ,   ω < ω o 3 π / 2   , ω > ω o { π α / 2 ,   ω < < ω o π ( α + β / 2 ) , α + β ( 0 , 2 ] , ω > > ω o π ( α + β / 2 ) + 2 π , α + β ( 2 , 4 ) , ω > > ω o G gi ( j ω ) is only determined by α when ω < < ω o , while is affected by both α and β when ω > > ω o . The high-frequency phase curve changes from π α / 2 to the more lagging direction when α + β ( 0 , 2 ] , and to the opposite direction when α + β ( 2 , 4 ) .
Phase crossover frequency, ω g ( L 1 + L 2 ) / ( L 1 L 2 C ) A sin [ ( 1 α / 2 ) π ] sin [ ( α + β / 2 ) π ] α + β If and only if α + β / 2 ( 1 , 2 ) , the phase-frequency characteristics of a fractional-order LCL filter has a ω g .
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

Xu, J.; Zeng, E.; Li, X.; He, G.; Liu, W.; Meng, X. Fractional-Order LCL Filters: Principle, Frequency Characteristics, and Their Analysis. Fractal Fract. 2024, 8, 38. https://doi.org/10.3390/fractalfract8010038

AMA Style

Xu J, Zeng E, Li X, He G, Liu W, Meng X. Fractional-Order LCL Filters: Principle, Frequency Characteristics, and Their Analysis. Fractal and Fractional. 2024; 8(1):38. https://doi.org/10.3390/fractalfract8010038

Chicago/Turabian Style

Xu, Junhua, Ermeng Zeng, Xiaocong Li, Guopeng He, Weixun Liu, and Xuanren Meng. 2024. "Fractional-Order LCL Filters: Principle, Frequency Characteristics, and Their Analysis" Fractal and Fractional 8, no. 1: 38. https://doi.org/10.3390/fractalfract8010038

APA Style

Xu, J., Zeng, E., Li, X., He, G., Liu, W., & Meng, X. (2024). Fractional-Order LCL Filters: Principle, Frequency Characteristics, and Their Analysis. Fractal and Fractional, 8(1), 38. https://doi.org/10.3390/fractalfract8010038

Article Metrics

Back to TopTop