Next Article in Journal
Percolation on Fractal Networks: A Survey
Next Article in Special Issue
On Reservoir Computing Approach for Digital Image Encryption and Forecasting of Hyperchaotic Finance Model
Previous Article in Journal
On the Existence Results for a Mixed Hybrid Fractional Differential Equations of Sequential Type
Previous Article in Special Issue
Research on Image Encryption Based on Fractional Seed Chaos Generator and Fractal Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison Study of Time-Domain Computation Methods for Piecewise Smooth Fractional-Order Circuit Systems

1
Hubei Provincial Collaborative Innovation Center for New Energy Microgrid, China Three Gorges University, Yichang 443000, China
2
College of Electrical Engineering and New Energy, China Three Gorges University, No. 8 Daxue Road, Yichang 443002, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(3), 230; https://doi.org/10.3390/fractalfract7030230
Submission received: 18 January 2023 / Revised: 2 March 2023 / Accepted: 3 March 2023 / Published: 4 March 2023
(This article belongs to the Special Issue Fractional-Order Circuits, Systems, and Signal Processing)

Abstract

:
The role of fractional calculus in circuit systems has received increased attention in recent years. In order to evaluate the effectiveness of time-domain calculation methods in the analysis of fractional-order piecewise smooth circuit systems, an experimental prototype is developed, and the effects of three typical calculation methods in different test scenarios are compared and studied in this paper. It is proved that Oustaloup’s rational approximation method usually overestimates the peak-to-peak current and brings in the pulse–voltage phenomenon in piecewise smooth test scenarios, while the results of the two iterative recurrence-form numerical methods are in good agreement with the experimental results. The study results are dedicated to provide a reference for efficiently deploying calculation methods in fractional-order piecewise smooth circuit systems. Some quantitative analysis results are concluded in this paper.

1. Introduction

In recent decades, the concepts of fractional calculus and the related techniques have been gaining momentum in circuit system fields [1,2,3]. It has been confirmed that exploring the potential fractional-order characteristics of electronic components is helpful for both the condition monitoring of components and the reliability design of circuit systems [4,5,6,7,8]. An increasing body of evidence suggests that the fractional-order characteristic is widely distributed in electronic components [9,10,11], such as ultracapacitors (UCs), lithium batteries, and non-solid electrolytic capacitors. Moreover, introducing fractional-order elements (or constant phase elements—CPEs) to traditional circuit and control systems can enhance the design flexibility [12,13,14].
In order to better describe the characteristics of circuit systems with fractional-order components, it is of necessity to develop a set of reliable calculation and analysis methods. Basically, the lumped parameter model with fractional-order differential equations is a well-established approach to quantify the characteristics of fractional-order circuit systems, and a number of calculation and analysis works have been proposed for such systems—definition-based methods [15], rational approximation methods [16,17,18], and numerical methods [19,20,21], to name but a few. The methods listed above have been applied to the modeling and analysis of a wide variety of fractional-order circuit systems [22,23,24].
To the best of the authors’ knowledge, most existing calculation and analysis methods are applied only to continuous fractional-order circuit systems, for instance, the discussion on the convergence issue of solutions for continuous cobweb models in Ref. [25] and the study on the existence and uniqueness of the solution of Hadamard fractional Itô–Doob stochastic integral equations in Ref. [26]. However, existing works are rarely applied to piecewise smooth ones [27,28,29]. In practice, the piecewise-smooth characteristic is common in circuit systems, especially for those with semiconductor switching devices, and such a characteristic may bring in effects on the accuracy of numerical methods. Therefore, it is worthwhile to compare and evaluate the effectiveness of existing methods in analyzing piecewise smooth fractional-order systems.
In allusion to the state of the art, this work tries to unfold the applicability of three typical calculation methods for piecewise smooth fractional-order circuit systems. The rest of the work is organized as follows: test scenarios are established in Section 1, in which a non-solid aluminum electrolytic capacitor is employed in the test bench since this kind of component has been confirmed to have frequency-related fractional-order characteristics in a wide frequency band [11]. Section 2 deduces the solutions of the test scenarios by different approaches. Section 3 compares and discusses the results of different approaches, while Section 4 concludes the work.

2. Test Bench Settings and Mathematical Model

2.1. Test Bench Settings

In a previous work, it was confirmed that the electrode surface in the capacitor has an infinite self-similar structure, and the particle distribution law in it has a long-tail effect under the electric field, which is suitable to be described by the fractional equivalent impedance model [11]. The internal structure and fractional-order equivalent impedance circuit is presented in Figure 1.
In the above figure, the symbol C is the nominal capacitance of the capacitor, while the symbol α is an estimated fractional order of 0 to 1, and the symbol R Ω is the equivalent series resistance of the capacitor. A test bench is established in this section, which contains such a capacitor. The schematic of this test bench is displayed in Figure 2.
The design of the test bench refers to the previous work [30], in which a type GPS-3303C 3-channel isolated dc power supply is adopted to provide power, and an oscilloscope is employed to record the data. In the platform, a type STP80NF70 power MOSFET S T is adopted, the gate-source voltage of which is controlled by a driving circuit. As a result, the test bench will work in charging and discharging the cycle working mode, so the circuit can be deemed as a piecewise smooth circuit system. On the right side of the red dash line in Figure 2, the resistor R 1 = 5 Ω provides a discharging path for the capacitor, while R 2 = 10 Ω is mainly used as the current sensing and limiting resistor. The circuit in the box with a red background is the fractional-order equivalent circuit of a 10 μ F Rubycon PX series non-solid aluminum electrolytic capacitor. The voltage and the current of the capacitor are v c a p ( t ) and i c a p ( t ) , respectively, and two probes of the oscillascope are used to observe them. In addition, the voltage of the equivalent CPE is assumed to be v c ( t ) . A glimpse of the experimental scene is given in Figure 3.

2.2. Mathematical Model

The test platform can be regarded as a charge–discharge circuit for the non-solid aluminum electrolytic capacitor. One can find that there is only one energy-storage component in the platform, so the interferences and errors are minimized compared with the schemes with rational approximation methods or those with magnetic elements. In the following parts, the test platform will be operated in two test scenarios to verify the applicability of different calculation and analysis methods.
According to the fractional-order equivalent impedance circuit of the capacitor in Figure 1, the current i c a p ( t ) and the voltage v c a p ( t ) of the capacitor satisfy the following relationships:
i c a p ( t ) = C · d α v c ( t ) d t α v c a p ( t ) = v c ( t ) + R Ω · C · d α v c ( t ) d t α ,
It is noteworthy that the nominal capacitance C is 10 μ F in this work, while R Ω and the fractional-order α related to the fractional-order characteristics of the capacitor can be identified by using the method proposed in a previous work [11], wherein their values are 0.8852 , 1.2629 Ω , respectively.
In order to validate different methods, we established two groups of test scenarios, the basic and the advanced. In the basic test group, we consider two operating modes. The first one is sinusoidal mode, where the power supply is a single phase AC power supply with a pre-determined frequency. The second one is step mode, where a DC power supply is employed as the source and the power MOSFET S T is turned on at time t = 0 with a constant forward voltage drop V d s ( o n ) = 0.7 V.
  • The state of the sinusoidal mode can be governed by
    v c ( t ) + ( R 2 + R Ω ) · C · d α v c ( t ) d t α = v i n ( t ) V d s ( o n ) ,
    where v i n ( t ) = K · sin ( ω t ) is a sinusoidal function. K is the predetermined magnitude of the signal, ω = 2 π f s is the angular frequency, and f s is the frequency of the sinusoidal function.
  • The state of the step mode can be governed by
    v c ( t ) + ( R 2 + R Ω ) · C · d α v c ( t ) d t α = V i n V d s ( o n ) ,
    in which V i n is a predetermined constant value.
In the first group of test scenarios, the test bench works under continuous mode or step response mode. This group are used to assess the feasibility of different calculation methods.
In the second group of the test scenario, the power MOSFET S T is controlled by a periodic square waveform with a duty ratio D = 0.5 and a variable frequency f s from 100 Hz to 1 kHz. The performance of the circuit under test in one steady-state cycle T s can be governed by
  • Charge performance in t [ k T s , k T s + D T s ] is
    v c a p ( t ) = V i n V d s ( o n ) R 2 · i c a p ( t ) ,
    where k is the k-th switching cycle, T s is the switching period, and
    V i n V d s ( o n ) = v c a p ( t ) + R 2 · i c a p ( t ) = v c ( t ) + ( R 2 + R Ω ) · C · d α v c ( t ) d t α .
  • Discharge performance in t [ k T s + D T s , ( k + 1 ) T s ] is
    v c a p ( t ) = ( R 1 + R 2 ) · i c a p ( t ) ,
    and
    0 = v c a p ( t ) + ( R 1 + R 2 ) · i c a p ( t ) = v c ( t ) + ( R 1 + R 2 + R Ω ) · C · d α v c ( t ) d t α .
In the second group of test scenarios, the capacitor of the test bench will be charged and discharged cycle by cycle, and thus the circuit under test is a typical piecewise smooth circuit system.
One can find that, to reveal the time-domain performance of these two scenarios, the voltage v C should be calculated. However, this voltage is an equivalent quantity of the CPE, and thus one cannot probe it directly and has to solve the fractional-order differential Equations (2)–(7), which are all fractional-order differential equations with constant coefficients. Rewrite these equations to the following generalized form:
a 1 i · d α x d t α + a 0 i · x = u i ,
where i = 1 , 2 , 3 , 4 , x = v C is the only state variable of the circuit, and the coefficients are: a 11 = a 12 = a 13 = ( R 2 + R Ω ) C , a 14 = ( R 1 + R 2 + R Ω ) C , a 01 = a 02 = a 03 = a 04 = 1 , u 1 = K · sin ( ω t ) V d s ( o n ) , u 2 = u 3 = V i n V d s ( o n ) , and u 4 = 0 .
All the coefficients are predetermined in their own time intervals. Meanwhile, the initial conditions of Equation (8) are continuous in their own domains of definition, and hence according to Theorem 3.2 of Ref. [31], fractional-order differential Equation (8) has unique solutions. At this time, the initial values problem of fractional-order differential equations arises.

3. Computation Approaches

In this section, the principles of some related techniques are introduced first. Then three different approaches are applied to calculate Equations (2)–(7), and the results are adopted for validation and comparison. The first calculation method is a numerical calculation method, that is, the fractional Adams–Bashforth–Moulton-type method (F-ABM) [21]. The second calculation method is based on the Gr u ¨ nwald–Letnikov (G-L) definition [15]. To obtain solutions by using the first and the second methods, the stroboscopic map technique should be applied [32]. The third calculation method is Oustaloup’s rational approximation method [16], which will be conducted by using the state-space averaging (SSA) technique [33].

3.1. Preliminaries: Principles of Some Related Techniques

3.1.1. Stroboscopic Map Technique

One can always apply the most calculation approaches to calculate fractional-order differential equations in continuous cases directly. However, in discontinuous cases, or the piecewise smooth case discussed in this work, one may have to face a different situation, for example, the sign function of f ( x ) = s i g n ( x ) , which has the value of 1 for all x > 0 and the value of 1 for all x < 0 , and Figure 4.
The sign function of f ( x ) = s i g n ( x ) appears in a variety of piecewise smooth fractional systems, which creates the initial value problem at the discontinuous point x = 0 , and one needs to physically interpret initial conditions at this point.
In the step mode of the basic test scenario, the proposed test bench also experiences a similar situation like the sign function. In addition, in the advanced test scenario, the capacitor experiences a recurrence of charging and discharging behaviors in each switching cycle, and thus the test bench will be in continuous state switching. Along with the on- and off operations of the power MOSFET S T , there is a set of discontinuous points. Accordingly, a technique called stroboscopic map should be employed in both analytical and numerical calculations. This technique has widely been adopted in the dynamic analysis of piecewise-smooth systems and switching power converters. By this technique, the dynamic behavior of the test circuit at each switching state S n will be collected in one switching cycle, that is, the solution x n of the previous switching state at time t n will be employed as the initial value of the next switching state, and thus a cycle-by-cycle calculation can be carried out. The principle of the stroboscopic map technique is depicted in Figure 5.

3.1.2. State-Space Averaging (SSA) Technique

In practice, some methods are difficult to be directly used in the advanced test scenario of this work. In allusion to this, the SSA technique can be exploited, by which the steady-state characteristic analysis can be achieved. By this technique, the average value of v c ( t ) in one switching cycle can be expressed by
v c T s = 1 T s t 0 t 0 + T s v c ( t ) d t
Accordingly, in the steady state, when the low-frequency hypothesis and the small ripple hypothesis are satisfied, the test platform in the discontinuous case can be governed by the following SSA model:
( R 2 + R Ω ) · C · d α v c T s d t α + v c T s + 1 d ( t ) · R 1 · C · d α v c T s d t α = d ( t ) V i n V d s ( o n ) ,
in which the terms v c T s and d ( t ) can be decomposed into DC terms V c and D, AC terms v ^ c and d ^ , as follows:
v c T s = V c + v ^ c d ( t ) = D + d ^

3.2. Solutions of the Test Bench

3.2.1. Solutions Obtained by G-L Definition

According to the G-L definition, the fractional-order derivative of state variable x ( t ) can be written in the following discrete form [15]:
d α x ( t ) d t α 1 h α j = 0 t / h p j α x t j h = 1 h α x t + j = 1 t / h p j α x t j h ,
where h is the predetermined discrete step size, and p j α can be deduced by
p 0 α = 1 , p j α = 1 1 + α j p j 1 α .
Introducing Equation (13) into Equation (12) leads to the following recurrence form solution:
x t = 1 a 1 i + h α a 0 i h α u i , t j = 1 t / h p j α x t j h .
Then, one can obtain G-L based solutions of Equation (8) by point-by-point iteration of Equation (14). In addition, the numerical solutions of charging and discharging state in the advanced test scenario are collected in each switching cycle by the aforementioned stroboscopic map technique.

3.2.2. Solutions Obtained by F-ABM Method

In the case of 1 α 1 , according to the definition of the F-ABM method [19], the initial value problem of the fractional-order system of Equation (8) can be determined by:
x ( n + 1 ) = x ( 0 ) + h α Γ ( 2 + α ) f t n + 1 , x P ( n + 1 ) + h α Γ ( 2 + α ) j = 0 n A j , n + 1 f t j , x ( j ) ,
where the term n is any integer, and A j , n + 1 is
A j , n + 1 = = n α + 1 ( n α ) ( n + 1 ) α , j = 0 = ( n j + 2 ) α + 1 + ( n j ) α + 1 2 ( n j + 1 ) α + 1 , 1 j n = 1 , j = n + 1
The predictor x P ( n + 1 ) in Equation (15) is
x P ( n + 1 ) = x 0 + 1 Γ ( α ) j = 0 n B j , n + 1 f t j , x ( j ) ,
in which the term B j , n + 1 is
B j , n + 1 = h α α ( n + 1 j ) α ( n j ) α .
Then the numerical solutions of Equation (8) can be obtained by Equation (15). Additionally, the numerical solutions of the charging and discharging state in the advanced test scenario will be collected in each switching cycle by the aforementioned stroboscopic map technique.

3.2.3. Solutions Obtained by Oustaloup’s Rational Approximation Method

By Oustaloup’s rational approximation method, a continuous filter can be designed in the frequency domain to achieve fractional-order calculus operations approximately. More specifically, by the Laplace transform, the fractional-order operation d α d t α can be transformed to an s-domain term s α , then the frequency domain characteristics of s α can be approximated in the pre-defined frequency interval [ ω L , ω H ] by an s-domain rational fraction function:
d α d t α K i = n n s + z i s + p i ,
where the gain K, zeros z i , and poles p i are
K = ( 0.9 · ω H ) α z i = ω L · μ ( ( i + N + 0.5 0.5 · α ) / ( 2 · N + 1 ) ) p i = ω L · μ ( ( k + N + 0.5 + 0.5 · α ) / ( 2 · N + 1 ) ) ,
respectively. The coefficient μ is
μ = ω H ω L .
Then the numerical solutions of Equation (8) can be obtained by encapsulating Equation (19) to a module in MATLAB/Simulink, such as those introduced in literature [16,17,18]. In addition, Equation (10) of the SSA technique should be applied. The block diagrams of the basic test scenario and the advanced test scenario are provided in Figure 6, where the block s α with red background is encapsulated according to Equation (19).

4. Results Comparison and Evaluation

According to the scenario settings and the related derivation, comparative works are carried out in this section, and the numerical results are obtained by the F-ABM method, G-L definition, and Oustaloup’s method, collected in groups. In addition, experiment waveforms are provided as reference.

4.1. Results of Calculation and Simulation

The first group of results are from the basic test scenario. By using the aforementioned three approaches, one can obtain numerical simulation results as depicted in Figure 7. The second group of results is from the advanced test scenario, with the simulation results given in Figure 8. In both cases, the blue solid line represents the voltage of power supply, while the red dash-and-dot line corresponds to the results of Oustaloup’s rational approximation method, the black dash line corresponds to the results of G-L definition, and the green solid line corresponds to the results of the F-ABM method.
In the above two figures, the sub-graphs with red and blue backgrounds refer to the enlarged parts.
In the sinusoidal case, it can be seen that the results of the three methods mesh well with each other in the steady state. There is a lag of around 0.5 ms between the voltage v c ( t ) and the voltage of the power supply, which is determined by the capacitance C = 10 μ F and the fractional order α = 0.8852 of the equivalent impedance model of the capacitor.
In step case, it can be seen that at the discontinuous point, the peak currents obtained by the three methods are almost the same, the value is around 1.0395 A. However, note that, the results obtained by three methods are slightly different at a steady state. In detail, the steady-state voltage calculated by Oustaloup’s method is about 0.1 V lower than the results obtained by the G-L definition and F-ABM method. In addition, the steady-state current calculated by Oustaloup’s method is about 0.01 A higher than the results obtained by the G-L definition and F-ABM method.
In advance test scenarios, it can be seen that the results of three methods mesh well with each other. However, note that there are spikes of the voltages v c ( t ) and v c a p ( t ) at discontinuous points if one uses Oustaloup’s method. With the acceleration frequency f s of charge–discharge state transition, this spike phenomenon becomes obvious. In addition, the steady-state current i c a p ( t ) and the voltage v c a p ( t ) calculated by the three methods are slightly different. With the acceleration frequency f s of the charge–discharge state transition, this difference becomes obvious.

4.2. Evaluation and Discussion

In order to further validate the calculation results, experiment waveforms are provided in Figure 9 and Figure 10.
In Figure 9, the blue solid lines are the voltage v 2 ( t ) of the resistor R 2 . It can be regarded as the sampling for the current i c a p ( t ) flowing through the capacitor. The red solid lines correspond to the voltage of the power supply, the green solid line corresponds to the voltage v c a p ( t ) of the capacitor.
It can be found that experimental waveforms are very similar to the theoretical calculation and simulation results. In order to assess the results of the calculation, statistical analysis was used to qualify different methods.
First of all, we compare the spike current i s p i k e , the averaged steady-state capacitor voltage v ¯ c a p , and the averaged steady-state capacitor current i ¯ c a p in the step case in Table 1.
From Table 1, it can be seen that most of the calculation results of the three methods have little deviation from the experimental results in the step case. However, compared with the other results, the averaged current i ¯ c a p obtained by Oustaloup’s method is too large. This current value is unreasonable because in both theory and practice, if the electric field applied to the capacitor is constant, the current of the capacitor should be in a trickle state.
Then, we compare the calculations results in sinusoidal case and advanced scenarios and list the peak-to-peak values obtained by different methods in Table 2.
From Table 2, it can be seen that the calculation results of the three methods have little deviation from the experimental results in the sinusoidal case, which means that the three methods are applicable for continuous situations.
However, compared with other results, the peak-to-peak current i p p obtained by Oustaloup’s method is larger in advanced scenarios. In addition, one can observe pulse voltage signals at some edge time points in the simulation results obtained by Oustaloup’s method, which cannot be observed in both experiments and the results of the other two methods. In principle, this pulse voltage phenomenon mainly occurs at the on and off operation edges of the power MOSFET S T . At these time points, the test circuit switches between charging and discharging states, just like the high-level and low-level signals of the sign function in Figure 6 undergo periodic switching. As a result, the state equation of the test circuit switches between Equation (5) and (7), and the test bench is a typical piecewise smooth circuit system.
Basically, the expressions of most frequency-domain rational approximation methods are the approximation of the amplitude–frequency or phase–frequency characteristics of ideal fractional calculus operators in a predetermined frequency band, and usually considering the situation of being continuously differentiable. However, at the boundary points of piecewise smooth scenarios, the state equation of the test bench is discontinuous and non-derivable. Therefore, the approximation results of the rational approximation method are not that effective. This pulse–voltage phenomena can also be explained by the Fourier transform of a step function.

5. Conclusions

This paper compares the effectiveness of the G-L definition based method, the F-ABM method, and Oustaloup’s rational approximation method for fractional-order piecewise smooth circuit systems. A test platform for verification is developed, in which a fractional-order component, a non-solid electrolytic capacitor, is adopted. The test platform can be operated under basic and advanced test scenarios. The basic scenario contains a continuous sinusoidal case and a step case, while the advanced scenario is the periodic charging and discharging operation of the capacitor.
Three computational methods work well in basic test scenarios, especially in the sinusoidal case, which indicates that these methods are effective in continuous situations. However, in the step case, the averaged current obtained by Oustaloup’s method is too large. In advance scenarios, the waveforms obtained by three methods are similar to the experimental results. Two iterative numerical calculation methods, G-L definition method and F-ABM method, perform well in f s = 100 Hz and f s = 200 Hz advanced scenarios, the deviation between their calculation results and the experimental results is small. However, applying Oustaloup’s method in advanced scenarios leads to large calculation deviations. Moreover, if one employs Oustaloup’s method in advanced scenarios, there is a pulse–voltage phenomenon when the circuit changes from the charging state to discharging state, which cannot be observed in both experiments and the results of the two iterative numerical calculation methods. The comparison and experimental verification results show that the G-L definition based method and F-ABM method are effective for fractional-order piecewise-smooth circuit systems.
The results provided by this research will provide more confidence for understanding the dynamics of real-world systems governed by fractional calculus. For those discrete- or continuous-time piecewise-smooth fractional-order dynamical systems, such as power electronic converters, which are typical piecewise smooth fractional-order circuit systems, if the fractional-order characteristics of passive components used in topology are considered, the discussion in the work provides a reference for effective deployment calculation methods.

Author Contributions

Conceptualization, methodology, writing—original draft preparation, and supervision, X.C.; validation, formal analysis, and data curation, F.Z.; writing—review and editing, and project administration, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Hubei Provincial Natural Science Foundation, China, grant number 2020CFB248.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Elwakil, A. Fractional-Order Circuits and Systems: An Emerging Interdisciplinary Research Area. IEEE Circuits Syst. Mag. 2010, 10, 40–50. [Google Scholar] [CrossRef]
  2. Zhang, B.; Shu, X. Fractional-Order Electrical Circuit Theory (CPSS Power Electronics Series), 1st ed.; Springer: Singapore, 2022; pp. 39–54. [Google Scholar]
  3. Petráš, I. Fractional-Order Nonlinear Systems; Springer: Berlin/Heidelberg, Germany, 2011; pp. 43–54. [Google Scholar]
  4. Zou, C.; Zhang, L.; Hu, X.; Wang, Z.; Wik, T.; Pecht, M. A review of fractional-order techniques applied to lithium-ion batteries, lead-acid batteries, and supercapacitors. J. Power Sources 2018, 390, 286–296. [Google Scholar] [CrossRef] [Green Version]
  5. Zhang, L.; Kartci, A.; Elwakil, A.; Bagci, H.; Salama, K.N. Fractional-Order Inductor: Design, Simulation, and Implementation. IEEE Access 2021, 9, 73695–73702. [Google Scholar] [CrossRef]
  6. Allagui, A.; Elwakil, A.S.; Fouda, M.E. Revisiting the Time-Domain and Frequency-Domain Definitions of Capacitance. IEEE Trans. Electron Devices 2021, 68, 2912–2916. [Google Scholar] [CrossRef]
  7. Huang, Y.; Chen, X. A fractional-order equivalent model for characterizing the interelectrode capacitance of MOSFETs. IEEE COMPEL—Int. J. Comput. Math. Electr. Electron. Eng. 2022, 41, 1660–1676. [Google Scholar] [CrossRef]
  8. Chen, X.; Pei, M. Enhancing parameter identification of electrochemical double layer capacitors by fractional-order equivalent impedance models and Levy flight strategy. Int. J. Circuit Theory Appl. 2022. [Google Scholar] [CrossRef]
  9. Mondal, D.; Biswas, K. Packaging of Single-Component Fractional Order Element. IEEE Trans. Device Mater. Reliab. 2013, 13, 73–80. [Google Scholar] [CrossRef]
  10. Allagui, A.; Freeborn, T.J.; Elwakil, A.S.; Elwakil, A.S.; Fouda, M.E.; Maundy, B.J.; Radwan, A.G.; Said, Z.; Abdelkareem, M.A. Review of fractional-order electrical characterization of supercapacitors. J. Power Sources 2018, 400, 457–467. [Google Scholar] [CrossRef]
  11. Chen, X.; Xi, L.; Zhang, Y.; Ma, H.; Huang, Y.; Chen, Y. Fractional techniques to characterize non-solid aluminum electrolytic capacitors for power electronic applications. Nonlinear Dyn. 2019, 98, 3125–3141. [Google Scholar] [CrossRef]
  12. Freeborn, T.J. A Survey of Fractional-Order Circuit Models for Biology and Biomedicine. IEEE J. Emerg. Sel. Top. Circuits Syst. 2013, 3, 416–424. [Google Scholar] [CrossRef]
  13. 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]
  14. Jiang, Y.; Zhang, B.; Shu, X.; Wei, Z. Fractional-order autonomous circuits with order larger than one. J. Adv. Res. 2020, 25, 217–225. [Google Scholar] [CrossRef]
  15. Scherer, R.; Kalla, S.L.; Tang, Y.; Huang, J. The Grünwald–Letnikov method for fractional differential equations. Comput. Math. Appl. 2011, 62, 902–917. [Google Scholar] [CrossRef] [Green Version]
  16. Oustaloup, A.; Levron, F.; Mathieu, B.; Nanot, F.M. Frequency-band complex noninteger differentiator: Characterization and synthesis. IEEE Trans. Circuits Syst. I-Fundam. Theory Appl. 2000, 47, 25–39. [Google Scholar] [CrossRef]
  17. Malti, R.; Victor, S. CRONE Toolbox for system identification using fractional differentiation models. IFAC-Pap. OnLine 2015, 48, 769–774. [Google Scholar] [CrossRef]
  18. Xue, D.; Chen, Y. System Simulation Techniques with MATLAB® and Simulink®, 1st ed.; John Wiley & Sons: West Sussex, UK, 2013; pp. 85–90. [Google Scholar]
  19. Diethelm, K. The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type; Springer: Berlin/Heidelberg, Germany, 2010; pp. 195–225. [Google Scholar]
  20. Momani, S.; Odibat, Z. Numerical approach to differential equations of fractional order. J. Comput. Appl. Math. 2007, 1, 96–110. [Google Scholar] [CrossRef] [Green Version]
  21. Owolabi, K.M.; Atangana, A. Numerical Methods for Fractional Differentiation; Springer Nature: Singapore, 2019; pp. 203–208. [Google Scholar]
  22. Chen, X.; Chen, Y. A Modeling and Analysis Method for Fractional-order DC-DC Converters. IEEE Trans. Power Electron. 2017, 32, 7034–7044. [Google Scholar] [CrossRef]
  23. 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]
  24. Li, G.; Lehman, B. Averaging Theory for Fractional Differential Equations. Fract. Calc. Appl. Anal. 2021, 24, 621–640. [Google Scholar] [CrossRef]
  25. Nagy, A.M.; Assidi, S.; Makhlouf, A.B. Convergence of solutions for perturbed and unperturbed cobweb models with generalized Caputo derivative. Bound. Value Probl. 2022, 89, 1–13. [Google Scholar] [CrossRef]
  26. Kahouli, O.; Makhlouf, A.B.; Mchiri, L.; Mchiri, L.; Rguigui, H. Hyers—Ulam stability for a class of Hadamard fractional Itô—Doob stochastic integral equations. Chaos. Soliton. Fract. 2021, 103, 2855–2866. [Google Scholar] [CrossRef]
  27. Danca, M.F. Synchronization of piecewise continuous systems of fractional order. Nonlinear Dyn. 2014, 78, 2065–2084. [Google Scholar] [CrossRef] [Green Version]
  28. Liu, X.; Hong, L.; Tang, D.; Yang, L. Crises in a fractional-order piecewise system. Nonlinear Dyn. 2021, 103, 2855–2866. [Google Scholar] [CrossRef]
  29. Alazman, I.; Alkahtani, B.S.T. Investigation of Novel Piecewise Fractional Mathematical Model for COVID-19. Fractal Fract. 2022, 6, 661. [Google Scholar] [CrossRef]
  30. Yang, Z.; Xi, L.; Zhang, Y.; Chen, X. An Online Parameter Identification Method for Non-solid Aluminum Electrolytic Capacitors. IEEE Trans. Circuits Syst. II Express Briefs 2022, 69, 3475–3479. [Google Scholar] [CrossRef]
  31. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Academic Press: London, UK, 1998; pp. 124–125. [Google Scholar]
  32. Tse, C.K.; Xi, L.; Bernardo, M.D. Complex Behavior in Switching Power Converters. Proc. IEEE 2002, 90, 768–781. [Google Scholar] [CrossRef] [Green Version]
  33. Erickson, R.W.; Maksimovic, D. Fundamentals of Power Electronics, 2nd ed.; Kluwer Academic Publishers: Norwell, MA, USA, 2001; pp. 213–246. [Google Scholar]
Figure 1. Internal structure and fractional-order equivalent impedance circuit model of non-solid aluminum electrolytic capacitors.
Figure 1. Internal structure and fractional-order equivalent impedance circuit model of non-solid aluminum electrolytic capacitors.
Fractalfract 07 00230 g001
Figure 2. Schematic of the test bench.
Figure 2. Schematic of the test bench.
Fractalfract 07 00230 g002
Figure 3. A glimpse of the experimental scene.
Figure 3. A glimpse of the experimental scene.
Fractalfract 07 00230 g003
Figure 4. A graph of typical s i g n ( x ) function.
Figure 4. A graph of typical s i g n ( x ) function.
Fractalfract 07 00230 g004
Figure 5. Principle of stroboscopic map technique.
Figure 5. Principle of stroboscopic map technique.
Fractalfract 07 00230 g005
Figure 6. Principle block diagrams: (a) basic test scenario, (b) advanced test scenario.
Figure 6. Principle block diagrams: (a) basic test scenario, (b) advanced test scenario.
Fractalfract 07 00230 g006
Figure 7. Comparison of basic test scenarios: (a) sinusoidal case, (b) step case.
Figure 7. Comparison of basic test scenarios: (a) sinusoidal case, (b) step case.
Fractalfract 07 00230 g007
Figure 8. Comparison of advanced test scenarios: (a) advanced scenario with f s = 100 Hz, (b) advanced scenario with f s = 200 Hz.
Figure 8. Comparison of advanced test scenarios: (a) advanced scenario with f s = 100 Hz, (b) advanced scenario with f s = 200 Hz.
Fractalfract 07 00230 g008
Figure 9. Experiment results of basic test scenarios: (a) sinusoidal case, (b) step case.
Figure 9. Experiment results of basic test scenarios: (a) sinusoidal case, (b) step case.
Fractalfract 07 00230 g009
Figure 10. Experiment results of advanced test scenarios: (a) advanced scenario with f s = 100 Hz, (b) advanced scenario with f s = 200 Hz.
Figure 10. Experiment results of advanced test scenarios: (a) advanced scenario with f s = 100 Hz, (b) advanced scenario with f s = 200 Hz.
Fractalfract 07 00230 g010
Table 1. Comparison of the spike current i s p i k e , the averaged steady-state capacitor voltage v ¯ c a p , and the averaged steady-state capacitor current i ¯ c a p in step case.
Table 1. Comparison of the spike current i s p i k e , the averaged steady-state capacitor voltage v ¯ c a p , and the averaged steady-state capacitor current i ¯ c a p in step case.
v ¯ cap i ¯ cap i spike
G-L definition11.2894 V≈1.29 mA1.0012 A
F-ABM11.3025 V≈−0.2 mA1.0308 A
Oustaloup’s11.2816 V0.0110 A1.0379 A
Experiments11.5570 V≈1 mA0.9869 A
Table 2. Comparison of peak-to-peak values of capacitor voltage v c p p and current i c p p .
Table 2. Comparison of peak-to-peak values of capacitor voltage v c p p and current i c p p .
In Sinusoidal CaseIn Advanced Scenario with
f s = 100 Hz
In Advanced Scenario with
f s = 200 Hz
G-L definition(0.5851 A, 21.9082 V)(1.6848 A, 11.0712 V)(1.6656 A, 10.6078 V)
F-ABM(0.5908 A, 21.9414 V)(1.6853 A, 10.7952 V)(1.6686 A, 10.0280 V)
Oustaloup’s(0.5514 A, 22.3235 V)(1.9266 A, 11.6760 V)(1.9088 A, 12.2848 V)
Experiments(0.6211 A, 21.1210 V)(1.6991 A, 11.6547 V)(1.6624 A, 11.3748 V)
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

Chen, X.; Zheng, F.; Wei, Y. A Comparison Study of Time-Domain Computation Methods for Piecewise Smooth Fractional-Order Circuit Systems. Fractal Fract. 2023, 7, 230. https://doi.org/10.3390/fractalfract7030230

AMA Style

Chen X, Zheng F, Wei Y. A Comparison Study of Time-Domain Computation Methods for Piecewise Smooth Fractional-Order Circuit Systems. Fractal and Fractional. 2023; 7(3):230. https://doi.org/10.3390/fractalfract7030230

Chicago/Turabian Style

Chen, Xi, Feng Zheng, and Yewen Wei. 2023. "A Comparison Study of Time-Domain Computation Methods for Piecewise Smooth Fractional-Order Circuit Systems" Fractal and Fractional 7, no. 3: 230. https://doi.org/10.3390/fractalfract7030230

APA Style

Chen, X., Zheng, F., & Wei, Y. (2023). A Comparison Study of Time-Domain Computation Methods for Piecewise Smooth Fractional-Order Circuit Systems. Fractal and Fractional, 7(3), 230. https://doi.org/10.3390/fractalfract7030230

Article Metrics

Back to TopTop