Next Article in Journal
Multi-Level Dependent-Chance Model for Hydropower Reservoir Operations
Next Article in Special Issue
Identification of Lightning Overvoltage in Unmanned Aerial Vehicles
Previous Article in Journal
Electric Vehicle Charging Model in the Urban Residential Sector
Previous Article in Special Issue
Methods for Enhancing the Electrical Properties of Epoxy Matrix Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximation of Fractional Order PIλDμ-Controller Transfer Function Using Chain Fractions

1
Institute of Power Engineering and Control Systems, Lviv Polytechnic National University, 12 Bandera Street, 79013 Lviv, Ukraine
2
Department of Electrical and Computer Engineering Fundamentals, Faculty of Electrical and Computer Engineering, Rzeszow University of Technology, 35-959 Rzeszow, Poland
3
Institute of Technical Engineering, State University of Technology and Economics in Jaroslaw, Czarnieckiego Str. 16, 37-500 Jaroslaw, Poland
4
Asseco Poland S.A., ul. Olchowa 14, 35-322 Rzeszow, Poland
*
Author to whom correspondence should be addressed.
Energies 2022, 15(13), 4902; https://doi.org/10.3390/en15134902
Submission received: 30 March 2022 / Revised: 22 May 2022 / Accepted: 30 June 2022 / Published: 4 July 2022

Abstract

:
The approximation of a fractional order PIλDμ-controller transfer function using a chain fraction theory is considered. Analytical expressions for the approximation of s ± α components of the transfer functions of PIλDμ-controllers were obtained through the application of the chain fraction theory. Graphs of transition functions and frequency characteristics of Dμ (α = μ = 0.5) and Iλ (α = λ = −0.5) parts for five different decomposition orders were obtained and analyzed. The results showed the possibility of applying the approximation of the PIλDμ-controller transfer function by the method of chain fractions with different valuesof λ and μ. For comparison, the transfer functions with the same order polynomials, obtained by the methods of Oustaloup transformation and chain fractions, were approximated for α = ±0.5. The analysis proved the advantages of using the chain fraction method to approximate the transfer function of the PIλDμ-controller. The performed approximation opens up the possibility of developing engineering methods for the technical implementation of PIλDμ-controllers. The accuracy of the same order transfer function approximation is higher when the method of chain fractions is used. It has been established that the adequacy of the frequency characteristics of the transfer functions obtained by the chain fraction method also depends on the approximation order.

1. Introduction

The use of proportional–integral–differential (PIλDμ) fractional order controllers in electromechanical systems (EMSs) improves the quality of transition processes [1] and increases their marginal stability [2] compared to similar systems using classical (integer order) controllers. With the increasing demand for adjustable average current (AC) drives, researchers are working on the wider use of fractional controllers PIλDμ with transfer functions (TFs), shownin (1).
The implementation of arbitrary controllers, including integer order PID-controllers, has been thoroughly worked out and is well known. At the same time, the implementation of fractional order controllers requires further studies and research.
One way of solving the problem of the technical implementation of fractional order controllers is by the use of the Oustaloup transformation (approximation) [3]. Due to this approximation, fractional order TFs are represented by integer order expressions, whose frequency characteristics correspond with a certain accuracy to similar characteristics of fractional order controllers. In this case, the problem of the technical implementation of the fractional controllers comes down to the construction of the controllers of orders of integers, which are usually regarded as controllers with proportional, differentiating, and integrating properties expressed as the following transmittance:
W c ( s ) = K p + 1 T i s λ + T d s μ .
The problem of implementing PIλDμ fractional order controllers in microcontrollers was partly solved in [4], where an original approach for the implementation of a digital differential–integral fractional order controller for electromechanical systems using first- and second-order Oustaloup transformation was proposed. Unfortunately, increasing the accuracy of the frequency characteristics of the approximated TFs (approximation accuracy) leads to a significant increase in the orders of the numerator and denominator polynomials of the obtained integer TFs. Studies on the controller model based on the microcontrollers Atmel ATMega2560 and Atmel ARM ATSAM3 × 8E showed that a satisfactory approximation accuracy can be obtained by applying the second-order Oustaloup transformation. Such approximation resulted in obtaining integer order TFs with fifth-order polynomials in the numerator and denominator. The application of the fifth-order Oustaloup transformation significantly increases the approximation accuracy, but in this case, the order of the numerator and denominator polynomials equals 11. Technical implementations of such equivalent integer order controllers is much more complicated.
Thus, the implementation of fractional order PIλDμ-controllers by approximating their fractional TFs and the study of their accuracy and speed in automated EMSs presents a topical problem, which requires a solution.
The aim of this workis to approximate the transfer functions of a fractional order PIλDμ-controller using the theory of chain fractions and comparative analysis of different ways of approximation of their fractional TFs, thus, paving the way for the development of engineering methods for PIλDμ-controller technical implementations.
The achievement of this goal requires solving the following tasks:
-
to obtain analytical expressions for the approximation of the s ± α components of PIλDμ-controller transfer functions through the application of the chain fraction theory;
-
to analyze the accuracy of the approximation of fractional order PIλDμ-controller TFs by the method of chain fractions by analyzing the transition functions and frequency characteristics of the units Dμ(α = μ = 0.5) and Iλ (α = λ = −0.5) for five different orders of decomposition by the method of chain fractions;
-
to compare the approximated TFs with polynomials of the same order obtained by the methods of Oustaloup and chain fractions for α = ±0.5 and to establish the advantages of using the method of chain fractions to approximate the TFs of the PIλDμ-controllers.

2. Application of the Oustaloup Transformation for Approximations of Fractional Order TFs

The Oustaloup transformation provides an approximation of the fractional order operator s ± α using transfer functions of the integer order in a given frequency range [ ω l , ω h ], where there is a lower and upper limit of the frequency range, respectively.
s α = ( ω u ω h ) α k = N k = N 1 + s / ω k 1 + s / ω k ,
where ω u = ω l ω h , N is the approximation order and ω k , ω k are the zeros and poles of the equivalent integer order TFs, respectively.
The calculation of the zeros and poles of the integer order approximating the TFs is made in accordance with the following expressions:
ω k = ω l ( ω h ω l ) ( k + N + 0.5 0.5 α ) / ( 2 N + 1 )
ω k = ω l ( ω h ω l ) ( k + N + 0.5 + 0.5 α ) / ( 2 N + 1 )
The coefficient of the integer 2N + 1 order approximating the TFs is as follows:
k = ( ω u ω h ) α ,
and its orderin approximating the TFs is n = 2N + 1.
Approximating TFs can be written as the ratio of n-order polynomials as follows:
W ( s ) = k b n s n + b n 1 s n 1 + + b 1 s + b 0 a n s n + a n 1 s n 1 + + a 1 s + a 0 = P ( s ) Q ( s ) .  
Below, there are exemplified appropriate integer order transfer function expressions obtained by applying the Oustaloup transformation (2) with n = 1, 2 with regard to fractional order differential and integral units, which are adjusted to the following form (5):
(1)
integral fractional unit W ( s ) = s 0.5
N = 1
s 0.5 = 0.1 s 3 + 4.867 s 2 + 10.49 s + 1 s 3 + 10.49 s 2 + 4.867 s + 0.1 ;
N = 2
s 0.5 = 0.1 s 5 + 7.497 s 4 + 76.85 s 3 + 121.8 s 2 + 29.85 s + 1 s 5 + 29.85 s 4 + 121.8 s 3 + 76.85 s 2 + 7.497 s + 0.1 ;  
(2)
differential fractional unit W ( s ) = s 0.5
N = 1
s 0.5 = 10 s 3 + 104.9 s 2 + 48.67 s + 1 s 3 + 48.67 s 2 + 104.9 s + 10 ;
N = 2
s 0.5 = 10 s 5 + 298.5 s 4 + 1218 s 3 + 768.5 s 2 + 74.97 s + 1 s 5 + 74.97 s 4 + 768.5 s 3 + 1218 s 2 + 298.5 s + 10

3. Application of Chain Fractions for Approximation ofthe Fractional Order Operator s α

Let us consider the possibility of applying the theory of chain fractions [5,6,7,8,9] to the approximate s ± α TF components of PIλDμ-controllers and their subsequent technical implementation. According to this theory [5,10,11,12,13,14], any power function can be extended to a chain fraction. In this case, for a binomial of arbitrary degree α, where the Laplace operator does not yet act as a variable, we can write the following expression:
( 1 + x ) α = 1 1 x α 1 + ( 1 + α ) x 2 + ( 1 α ) x 3 + ( 2 + α ) x 2 + ( 2 α ) x 5 + .
For simplicity of the analysis, we limit ourselves to the number of elements of the chain fraction in order to obtain the highest order of the approximated TF polynomials equal to one. Then, we obtain the following:
( 1 + x ) α = 1 1 α x 1 + ( 1 + α ) x 2 = x + α x + 2 x α x + 2 .
Let us replace the above-mentioned Expression (11) by using the Laplace operator as a variable x = s 1 as follows:
( s ) α = s α + α s + 1 s + α α s + 1
Thus, if α = 0.5, we obtain the following:
( s ) 0.5 = 3 s + 1 s + 3 .
Accordingly, if α = −0.5, we obtain the following:
( s ) 0.5 = s + 3 3 s + 1 .
The obtained results are much simpler than the results obtained on the basis of the Oustaloup transformation.
For higher orders of TFs from n = 2 to n = 5, under the condition that 0 < α < 1, we obtain the following:
If n = 2, we obtain the following:
( 1 + x ) α = 1 1 α x 1 + ( 1 + α ) x 2 + ( 1 α ) x 3 + ( 2 + α ) x 2 = = α 2 x 2 + 3 α x 2 + 6 α x + 2 x 2 + 12 x + 12 α 2 x 2 3 α x 2 6 α x + 2 x 2 + 12 x + 12 .
Having replaced x = s 1 , Expression (14) looks as follows:
( s ) α = α 2 s 2 + 3 α s 2 + 2 s 2 2 α 2 s + 8 s + α 2 3 α + 2 α 2 s 2 3 α s 2 + 2 s 2 2 α 2 s + 8 s + α 2 + 3 α + 2 .
Thus, if α = 0.5, we obtain the following: it is corect
( s ) 0.5 = 5 s 2 + 10 s + 1 s 2 + 10 s + 5 .
Accordingly, if α = −0.5, we obtain the following:
( s ) 0.5 = s 2 + 10 s + 5 5 s 2 + 10 s + 1 .
The obtained results are simpler than the results obtained on the basis of the Oustaloup transformation (n = 1).
In the same way, TFs from n = 3 to n = 5 are calculated as follows. If n = 3, we obtain the following:
( 1 + x ) α = 1 1 α x 1 + ( 1 + α ) x 2 + ( 1 α ) x 3 + ( 2 + α ) x 2 + ( 2 α ) x 5 + ( 3 + α ) x 2 = α 3 x 3 + 6 α 2 x 3 + 12 α 2 x 2 + 11 α x 3 + 60 α x 2 + 60 α x + α 3 x 3 6 α 2 x 3 12 α 2 x 2 + 11 α x 3 + 60 α x 2 + 60 α x 6 x 3 + 72 x 2 + 180 x + 120 6 x 3 72 x 2 180 x 120 .
Having replaced x = s 1 , Expression (17) looks as follows:
( s ) α = α 3 s 3 + 6 α 2 s 3 + 11 a s 3 + 6 s 3 3 α 3 s 2 6 α 2 s 2 + 27 α s 2 + 6 α 2 s 3 α 3 s 3 11 a s 3 + 6 s 3 + 3 α 3 s 2 6 α 2 s 2 27 α s 2 + + 54 s 2 + 3 α 3 s 6 α 2 s 27 α s + 54 s α 3 + 6 α 2 11 α + 6 + 54 s 2 3 α 3 s 6 α 2 s + 27 α s + 54 s + α 3 + 6 α 2 + 11 α + 6 ,
If n = 4, we obtain the following:
( 1 + x ) α = = 1 1 α x 1 + ( 1 + α ) x 2 + ( 1 α ) x 3 + ( 2 + α ) x 2 + ( 2 α ) x 5 + ( 3 + α ) x 2 + ( 3 α ) x 7 + ( 4 + α ) x 2 = = α 4 x 4 + 10 α 3 x 4 + 20 α 3 x 3 + 35 α 2 x 4 + 180 α 2 x 3 + α 4 x 4 10 α 3 x 4 20 α 3 x 3 + 35 α 2 x 4 + 180 α 2 x 3 + + 180 α 2 x 2 + 50 α x 4 + 520 α x 3 + 1260 α x 2 + 840 α x + 24 x 4 + + 180 α 2 x 2 50 α x 4 520 α x 3 1260 α x 2 840 α x + 24 x 4 + + 480 x 3 + 2160 x 2 + 3360 x + 1680 + 480 x 3 + 2160 x 2 + 3360 x + 1680
Having replaced x = s 1 , Expression (18) looks as follows:
( s ) α = α 4 s 4 + 10 α 3 s 4 + 35 α 2 s 4 + 50 α s 4 + 24 s 4 4 α 4 s 3 20 α 3 s 3 + 40 α 2 s 3 + α 4 s 4 10 α 3 s 4 + 35 α 2 s 4 50 α s 4 + 24 s 4 4 α 4 s 3 + 20 α 3 s 3 + 40 α 2 s 3 + 320 α s 3 + 384 s 3 + 6 α 4 s 2 150 α 2 s 2 + 864 s 2 4 α 4 s + 20 α 3 s + 40 α 2 s 320 α s 3 + 384 s 3 + 6 α 4 s 2 150 α 2 s 2 + 864 s 2 4 α 4 s 20 α 3 s + 40 α 2 s + 320 α s + 384 s + α 4 10 α 3 + 35 α 2 50 α + 24 + 320 α s + 384 s + α 4 + 10 α 3 + 35 α 2 + 50 α + 24
If n = 5, we obtain the following:
( 1 + x ) α = = 1 1 α x 1 + ( 1 + α ) x 2 + ( 1 α ) x 3 + ( 2 + α ) x 2 + ( 2 α ) x 5 + ( 3 + α ) x 2 + ( 3 α ) x 7 + ( 4 + α ) x 2 + ( 4 α ) x 9 + ( 5 + α ) x 2 = α 5 x 5 + 15 α 4 x 5 + 30 α 4 x 4 + 85 α 3 x 5 + 420 α 3 x 4 + α 5 x 5 15 α 4 x 5 30 α 4 x 4 + 85 α 3 x 5 + 420 α 3 x 4 + + 420 α 3 x 3 + 225 α 2 x 5 + 2130 α 2 x 4 + 5040 α 2 x 3 + 3360 α 2 x 2 + + 420 α 3 x 3 225 α 2 x 5 2130 α 2 x 4 5040 α 2 x 3 3360 α 2 x 2 + + 274 α x 5 + 4620 α 1 x 4 + 19740 α 1 x 3 + 30240 α 1 x 2 + 15120 α 1 x 1 + + 274 α x 5 + 4620 α 1 x 4 + 19740 α 1 x 3 + 30240 α 1 x 2 + 15120 α 1 x 1 + 120 x 5 + 3600 x 4 + 25200 x 3 + 67200 x 2 + 75600 x 1 + 30240 120 x 5 3600 x 4 25200 x 3 67200 x 2 75600 x 1 30240
Having replaced x = s 1 , Expression (19) looks as follows:
( s ) α = α 5 s 5 + 15 α 4 s 5 + 85 α 3 s 5 + 225 α 2 s 5 + 274 α s 5 + 120 s 5 5 α 5 s 4 15 α 4 s 5 α 5 s 5 85 α 3 s 5 + 225 α 2 s 5 274 α s 5 + 120 s 5 + 5 α 5 s 4 45 α 4 s 4 5 α 3 s 4 + 1005 α 2 s 4 + 3250 α s 4 + 3000 s 4 + 10 α 5 s 3 + 30 α 4 s 3 45 α 4 s 4 + 5 α 3 s 4 + 1005 α 2 s 4 3250 α s 4 + 3000 s 4 10 α 5 s 3 + 30 α 4 s 3 + 410 α 3 s 3 1230 α 2 s 3 + 4000 α s 3 + 12000 s 3 10 α 5 s 2 + 30 α 4 s 2 + 410 α 3 s 2 + 410 α 3 s 3 1230 α 2 s 3 4000 α s 3 + 12000 s 3 + 10 α 5 s 2 + 30 α 4 s 2 410 α 3 s 2 1230 α 2 s 2 4000 α 1 s 2 + 12000 s 2 + 5 α 5 s 1 45 α 4 s 1 + 5 α 3 s 1 + 1005 α 2 s 1 1230 α 2 s 2 + 4000 α 1 s 2 + 12000 s 2 5 α 5 s 1 45 α 4 s 1 5 α 3 s 1 + 1005 α 2 s 1 + 3250 α 1 s 1 + 3000 s 1 s 5 + 15 s 4 85 s 3 + 225 s 2 274 s + 120 + 3250 α 1 s 1 + 3000 s 1 + s 5 + 15 s 4 + 85 s 3 + 225 s 2 + 274 s + 120
So, for example, if α = 0.5 we obtain the following:
n = 3:
( s ) 0.5 = 7 s 3 + 35 s 2 + 21 s + 1 s 3 + 21 s 2 + 35 s + 7 ;
n = 4:
( s ) 0.5 = 9 s 4 + 84 s 3 + 126 s 2 + 36 s + 1 s 4 + 36 s 3 + 126 s 2 + 84 s + 9 ;
n = 5:
( s ) 0.5 = 11 s 5 + 165 s 4 + 462 s 3 + 330 s 2 + 55 s + 1 s 5 + 55 s 4 + 330 s 3 + 462 s 2 + 165 s + 11 .
Accordingly, if α = 0.5 , we obtain the following:
n = 3:
( s ) 0.5 = s 3 + 21 s 2 + 35 s + 7 7 s 3 + 35 s 2 + 21 s + 1 ;
n = 4:
( s ) 0.5 = s 4 + 36 s 3 + 126 s 2 + 84 s + 9 9 s 4 + 84 s 3 + 126 s 2 + 36 s + 1 ;
n = 5:
( s ) 0.5 = s 5 + 55 s 4 + 330 s 3 + 462 s 2 + 165 s + 11 11 s 5 + 165 s 4 + 462 s 3 + 330 s 2 + 55 s + 1 .
Figure 1 shows the graphs of the transition functions of theDμ parts (α = μ = 0.5), which are described by the TFs (12), (15), (20)–(22). Figure 2 demonstrates the frequency characteristics (Bode diagrams) of these TFs.
Similarly, the graphs of the transition functions of theIλ parts (α = λ = −0.5), which are described by the TFs (13), (16), (23)–(25) (Figure 3), and the frequency characteristics of these TFs (Figure 4) were obtained.
Using the approximations of theIλ and Dμ parts, as components of the fractional order controllers, we will approximate the PIλDμ-controllers with different λ and μ values of their TFs. For example, if α = ±0.5, then (1) turns into the following:
W ( ω ) = 1 + 1 ( j ω ) 0.5 + ( j ω ) 0.5
For five different decomposition orders of the PIλDμ-controller (26), the frequency characteristics of such controller were obtained by the method of chain fractions (Figure 5). Here are the corresponding frequency characteristics (red graphs), which are calculated for the fractional TFs using the Mathcad application package. We consider these graphs as reference.
It follows from the obtained graphs that the position of the frequency characteristics in relation to the reference characteristics for frequencies at most 1rad/s significantly depends on the order of approximations. The difference between the frequency and the reference characteristics starts to considerably decrease from the fifth order of approximation. Thus, in this case, the standard deviation σ for the amplitude–frequency characteristic is σ1 = 0.0596 and for the phase–frequency characteristic is σ2 = 6.316.
If the frequency exceeds 1 rad/s, then for each order of approximation there is a frequency range, for which the lowest value of σ is provided. Thus, if n = 1, then the smallest value σ1 = 1087 for the amplitude–frequency characteristic remains in the frequency range between 1rad/s and 10rad/s. With respect to the phase–frequency characteristics, the lowest value σ2 = 3.85 remains in the frequency range between 1 rad/s and 4 rad/s. If we apply the approximation order n = 5, then the same values of σ can be provided in the frequency range between 1 rad/s and 100 rad/s for the amplitude–frequency characteristic and between 1 and 40 rad/s for the phase–frequency characteristic.
Therefore, it is recommended to choose the order of approximation based on the analysis of the development of the expected frequency effects of the system.
For clarity of the controller operation, Figure 6 shows its transition processes for the part described by Expression (26).
Due to the greater informativeness of the frequency characteristics, we will further analyze these characteristics.
Similarly to the case of α = ±0.5, the frequency characteristics for the TF were obtained when α = ± 0.1, as follows:
W ( ω ) = 1 + 1 ( j ω ) 0.1 + ( j ω ) 0.1 .
These graphs are shown in Figure 7. Similarly, the studies were performed to change α in the range from ±0.1 to ±0.9 in various combinations.
For all investigated systems, the same conclusions are made for obtaining the specified accuracy of the frequency characteristics in the appropriate frequency range, as those made for the case of α = ±0.5.
The obtained results showed the possibility of applying the TF approximation of the PIλDμ-controller by the method of chain fractions with different values for λ and μ.
In order to compare the methods of Oustaloup and chain fractions, we analyzed the frequency characteristics approximated from the condition n = 5 for the method of Oustaloup and n = 4 for the method of chain fractions if α = ±0.5 is inherent in the TF.
The chain fraction method is as follows:
W 4 c f ( ω ) = 1 + ( j ω ) 4 + 36 ( j ω ) 3 + 126 ( j ω ) 2 + 84 ( j ω ) + 9 9 ( j ω ) 4 + 84 ( j ω ) 3 + 126 ( j ω ) 2 + 36 ( j ω ) + 1 + + 9 ( j ω ) 4 + 84 ( j ω ) 3 + 126 ( j ω ) 2 + 36 ( j ω ) + 1 ( j ω ) 4 + 36 ( j ω ) 3 + 126 ( j ω ) 2 + 84 ( j ω ) + 9 .
The Oustaloup method is as follows:
W 5 o ( ω ) = 1 + 0.1 ( j ω ) 5 + 7.497 ( j ω ) 4 + 76.85 ( j ω ) 3 + 121.8 ( j ω ) 2 + 29.85 ( j ω ) + 1 1 ( j ω ) 5 + 29.85 ( j ω ) 4 + 121.8 ( j ω ) 3 + 76.85 ( j ω ) 2 + 7.497 ( j ω ) + 0.1 +
+ 10 ( j ω ) 5 + 298.5 ( j ω ) 4 + 1218 ( j ω ) 3 + 768.5 ( j ω ) 2 + 74.97 ( j ω ) + 1 ( j ω ) 5 + 74.97 ( j ω ) 4 + 768.5 ( j ω ) 3 + 1218 ( j ω ) 2 + 298.5 ( j ω ) + 10 .
The frequency characteristics for these two cases are shown in Figure 8.
It follows from the obtained graphs that σ1 ≈ 9.09, σ2 = 14.13 in almost the whole investigated frequency range for the method of chain fractions and that σ1 ≈ 9.42, σ2 = 14.31 in this frequency range for the Oustaloup method.
The analysis of the obtained results showed almost the same accuracy of approximation of the TF PIλDμ-controllers by both methods. It should be noted that this result is achieved by the method of chain fractions of the lower order than the order of the Oustaloup method. This simplifies the implementation of the PIλDμ-controller.

4. Conclusions

  • Analytical expressions based on the application of the chain fraction theory to approximate the components of the transfer functions of PIλDμ-controllers s ± α have been obtained.
  • The accuracy of the approximation of the fractional transfer function by the method of chain fractions has been researched by analyzing the transition functions and frequency characteristics of the units Dμ and Iλ for five different orders of decomposition by the method of chain fractions.
  • The comparison of the transfer functions, approximated by different order polynomials, has been conducted. As a result, the possibility of obtaining the same approximation accuracy has been established under the condition of using the method of chain fractions of the lower order. Using this approximation method simplifies the implementation of PIλDμ.
  • The analysis proved the possibility of using the method of chain fractions to approximate the transfer function of PIλDμ-controllers. The accuracy of the same order transfer function approximation is higher when the method of chain fractions is used. Based on the analysis of the system’s performance of the expected frequency effects, it is necessary to choose the appropriate approximation order to ensure the specified accuracy of the approximation.

Author Contributions

Conceptualization, Y.M., B.K. (Bohdan Kopchak) and D.M.; methodology, Y.M., T.K. and B.K. (Bogdan Kwiatkowski); software B.K. (Bohdan Kopchak), Y.M., T.K. and M.K.; validation, Y.M., B.K. (Bogdan Kwiatkowski) and D.M.; formal analysis, B.K. (Bogdan Kwiatkowski), Y.M., D.M. and M.K.; investigation, D.M. and B.K. (Bogdan Kwiatkowski); resources and data curation, Y.M., B.K. (Bogdan Kwiatkowski) and D.M.; writing—original draft preparation, T.K., Y.M. and D.M.; writing—review and editing, Y.M., B.K. (Bohdan Kopchak) and D.M.; visualization, B.K. (Bogdan Kwiatkowski), T.K. and D.M.; supervision, B.K. (Bogdan Kwiatkowski), Y.M. and D.M.; project administration Y.M., B.K. (Bohdan Kopchak) and D.M.; funding acquisition, D.M. All authors have read and agreed to the published version of the manuscript.

Funding

This project is financed by the Minister of Education and Science of the Republic of Poland within the “Regional Initiative of Excellence” program for years 2019–2022. Project number 027/RID/2018/19, amount granted 11 999 900 PLN.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kopchak, B.; Marushchak, Y.; Kushnir, A. Devising a procedure for the synthesis of electromechanical systems with cascade-enabled fractional-order controllers and their study. East.-Eur. J. Enterp. Technol. Inf. Technol. Ind. Control. Syst. 2019, 5, 65–71. [Google Scholar] [CrossRef] [Green Version]
  2. Lozynskyy, O.; Mazur, D.; Marushchak, Y.; Kwiatkowski, B.; Lozynskyy, A.; Kwater, T.; Kopchak, B.; Hawro, P.; Kasha, L.; Pękala, R.; et al. Formation of characteristic polynomials on the basis of fractional powers j of dynamic systems and stability problems of such systems. Energies 2021, 14, 7374. [Google Scholar] [CrossRef]
  3. Oustaloup, A.; Levron, F.; Nanot, F.; Mathieu, B. Frequency band complex non-integer differentiator: Characterization and synthesis. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 2000, 47, 25–40. [Google Scholar] [CrossRef]
  4. Kopchak, B. Development of fractional order differential-integral controller by using Oustaloup transformation. In Proceedings of the 2016 XII International Conference on Perspective Technologies and Methods in MEMS Design (MEMSTECH), Lviv, Ukraine, 20–24 April 2016; pp. 62–65. [Google Scholar] [CrossRef]
  5. Khovanskii, A.N. The Application of Continued Fractions and Their Generalizations to Problems in Approximation Theory; P. Noordhoff: Groningen, The Netherlands, 1963. [Google Scholar]
  6. Krishna, B.T.; Reddy, K.V.V.S. Active and passive realization of fractance device of order ½. Act. Passiv. Electron. Compon. 2008, 2008, 369421. [Google Scholar] [CrossRef] [Green Version]
  7. Kumar, R.; Perumalla, S.; Vista, J.; Ranjan, A. Realization of single and double cole tissue models using higher order approximation. In Proceedings of the 1st International Conference on Electronics, Materials Engineering and Nano-Technology (IEMENTech), Science City, Kolkata, India, 28–29 April 2017; pp. 1–4. [Google Scholar]
  8. Pu, Y.; Yuan, X.; Liao, K. Structuring analog fractance circuit for 1/2 order fractional calculus. In Proceedings of the 6th International Conference on ASIC (ASICON’05), Shanghai, China, 24–27 October 2005; Volume 2, pp. 1136–1139. [Google Scholar]
  9. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Section 5.2. Evaluation of Continued Fractions. In Numerical Recipes: The Art of Scientific Computing, 3rd ed.; Cambridge University Press: New York, NY, USA, 2007. [Google Scholar]
  10. Cuyt, A.; Brevik Petersen, V.; Verdonk, B.; Waadeland, H.; Jones, W.B. Handbook of Continued Fractions for Special Functions; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  11. Rockett, A.M.; Szüsz, P. Continued Fractions; World Scientific Press: Singapore, 1992. [Google Scholar]
  12. Scheinerman, E.; Pickett, T.; Coleman, A. Another Continued Fraction for π. Am. Math. Mon. 2008, 115, 930–933. [Google Scholar]
  13. Collins, D.C. Continued Fractions. MIT Undergrad. J. Math. 1999, 1, 11–20. [Google Scholar]
  14. Martin, R.M. Electronic Structure: Basic Theory and Practical Methods; Cambridge University Press: Cambridge, UK, 2004; p. 557. [Google Scholar]
Figure 1. Transition processes of parts (12)—1, (15)—2, (20)—3, (21)—4, (22)—5.
Figure 1. Transition processes of parts (12)—1, (15)—2, (20)—3, (21)—4, (22)—5.
Energies 15 04902 g001
Figure 2. Frequency characteristics of parts (12)—1, (15)—2, (20)—3, (21)—4, (22)—5.
Figure 2. Frequency characteristics of parts (12)—1, (15)—2, (20)—3, (21)—4, (22)—5.
Energies 15 04902 g002
Figure 3. Transition processes of parts (13)—1, (16)—2, (23)—3, (24)—4, (25)—5.
Figure 3. Transition processes of parts (13)—1, (16)—2, (23)—3, (24)—4, (25)—5.
Energies 15 04902 g003
Figure 4. Frequency characteristics of parts (13)—1, (16)—2, (23)—3, (24)—4, (25)—5.
Figure 4. Frequency characteristics of parts (13)—1, (16)—2, (23)—3, (24)—4, (25)—5.
Energies 15 04902 g004
Figure 5. Reference frequency characteristics and frequency characteristics for five approximation orders, if α = ±0.5.
Figure 5. Reference frequency characteristics and frequency characteristics for five approximation orders, if α = ±0.5.
Energies 15 04902 g005
Figure 6. Transition processes of the PIλDμ-controllers (26) for n = 1, 2, …, 5.
Figure 6. Transition processes of the PIλDμ-controllers (26) for n = 1, 2, …, 5.
Energies 15 04902 g006
Figure 7. Reference frequency characteristics 1 and frequency characteristics for the fifth approximation order if α = ±0.1; 2 for n = 5.
Figure 7. Reference frequency characteristics 1 and frequency characteristics for the fifth approximation order if α = ±0.1; 2 for n = 5.
Energies 15 04902 g007
Figure 8. Reference frequency characteristics 1 and frequency characteristics for the TF, approximated by the methods of Oustaloup 2 and chain fractions 3.
Figure 8. Reference frequency characteristics 1 and frequency characteristics for the TF, approximated by the methods of Oustaloup 2 and chain fractions 3.
Energies 15 04902 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marushchak, Y.; Mazur, D.; Kwiatkowski, B.; Kopchak, B.; Kwater, T.; Koryl, M. Approximation of Fractional Order PIλDμ-Controller Transfer Function Using Chain Fractions. Energies 2022, 15, 4902. https://doi.org/10.3390/en15134902

AMA Style

Marushchak Y, Mazur D, Kwiatkowski B, Kopchak B, Kwater T, Koryl M. Approximation of Fractional Order PIλDμ-Controller Transfer Function Using Chain Fractions. Energies. 2022; 15(13):4902. https://doi.org/10.3390/en15134902

Chicago/Turabian Style

Marushchak, Yaroslav, Damian Mazur, Bogdan Kwiatkowski, Bohdan Kopchak, Tadeusz Kwater, and Maciej Koryl. 2022. "Approximation of Fractional Order PIλDμ-Controller Transfer Function Using Chain Fractions" Energies 15, no. 13: 4902. https://doi.org/10.3390/en15134902

APA Style

Marushchak, Y., Mazur, D., Kwiatkowski, B., Kopchak, B., Kwater, T., & Koryl, M. (2022). Approximation of Fractional Order PIλDμ-Controller Transfer Function Using Chain Fractions. Energies, 15(13), 4902. https://doi.org/10.3390/en15134902

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

Article Metrics

Back to TopTop