Next Article in Journal
Homotopy Analysis Method for a Fractional Order Equation with Dirichlet and Non-Local Integral Conditions
Next Article in Special Issue
Sufficient Criteria for the Absence of Global Solutions for an Inhomogeneous System of Fractional Differential Equations
Previous Article in Journal
Fractional Cauchy Problems for Infinite Interval Case-II
Previous Article in Special Issue
An Analytical Technique to Solve the System of Nonlinear Fractional Partial Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplified Fractional Order Controller Design Algorithm

by
Eva-Henrietta Dulf
1,2
1
Technical University of Cluj-Napoca, Faculty of Automation and Computer Science, Department of Automation, Memorandumului Str. 28, 400014 Cluj-Napoca, Romania
2
Physiological Controls Research Center, Óbuda University, H-1034 Budapest, Hungary
Mathematics 2019, 7(12), 1166; https://doi.org/10.3390/math7121166
Submission received: 21 August 2019 / Revised: 28 October 2019 / Accepted: 20 November 2019 / Published: 2 December 2019

Abstract

:
Classical fractional order controller tuning techniques usually establish the parameters of the controller by solving a system of nonlinear equations resulted from the frequency domain specifications like phase margin, gain crossover frequency, iso-damping property, robustness to uncertainty, etc. In the present paper a novel fractional order generalized optimum method for controller design using frequency domain is presented. The tuning rules are inspired from the symmetrical optimum principles of Kessler. In the first part of the paper are presented the generalized tuning rules of this method. Introducing the fractional order, one more degree of freedom is obtained in design, offering solution for practically any desired closed-loop performance measures. The proposed method has the advantage that takes into account both robustness aspects and desired closed-loop characteristics, using simple tuning-friendly equations. It can be applied to a wide range of process models, from integer order models to fractional order models. Simulation results are given to highlight these advantages.

1. Introduction

Fractional calculus has become very useful over the last years due to its many applications in almost all applied sciences. There are applications in acoustic wave propagation in inhomogeneous porous material, diffusive transport, fluid flow, dynamical processes in self-similar structures, dynamics of earthquakes, optics, geology, viscoelastic materials, biosciences, bioengineering, medicine, economics, probability and statistics, astrophysics, chemical engineering, physics, splines, tomography, fluid mechanics, electromagnetic waves, nonlinear control, signal processing, control of power electronics, converters, chaotic dynamics, polymer science, proteins, polymer physics, electrochemistry, statistical physics, thermodynamics, neural networks, etc. [1,2,3,4,5,6,7].
Many researchers consider this mathematical tool very useful and provide significant contributions in their field. The work of Podlubny [8] had a major impact in control engineering. He proposed a generalization of the PID controller, namely the PIλDμ controller, involving an integrator of order λ and a differentiator of order μ and of Oustaloup [9], who introduced the CRONE approach for these systems. They also demonstrated that the response of this type of controller is better, in comparison with the classical PID controller, when used for the control of fractional order systems. There are also numerous different forms of fractional order controllers available, proper in some particular cases. For example, in [10,11] is presented a particular fractional-order control scheme, the PDD1/2, which derives from the classical PD scheme with the introduction of the half-derivative term.
The fractional order controller design techniques are in general based on extensions of the classical PID control theory, with an emphasis on the increased flexibility in the tuning strategy resulting better control performances as compared to classical control tuning methods.
Several works approach the tuning of the fractional order PID controller through frequency domain specifications, firstly described by [12]. Tuning the fractional order controller implies solving the system of nonlinear equations composed of the design constraints, usually by optimization method [13] or by approximation methods [14]. There are also available tuning algorithms using time domain cost functions and optimization routines [15], constrained integral optimization methods, the fractional extension of the MIGO algorithm designed by Astrom et al. as an improvement to the Ziegler–Nichols rules [16,17] and auto-tuning methods [18]. Several applications use fractional order techniques, as it is presented in the survey [2].
All of these works use one of the several definitions of fractional order derivative (or integral) described above:
The Riemann–Liouville definition [12]:
D c α f ( t ) = 1 Γ ( α ) c t ( t τ ) α 1 f ( τ ) d τ ,   t > c ,   α R + ,
where Γ ( n ) = 0 t n 1 e t d t is the Euler’s Gamma function which is a generalization of a factorial and n R + is an extension of the fractional integral.
The Caputo definition [12]:
D α C f ( t ) = 1 Γ ( m α ) 0 t f ( m ) ( τ ) ( t τ ) α m + 1 d τ ,   m 1 < α < m ,   m N .
The Grünwald–Letnikov’s definition of the fractional-order derivative [12]:
GL D α f ( t ) = k = 0 m f ( k ) ( 0 + ) t k α Γ ( m + 1 α ) + 1 Γ ( m + 1 α ) 0 t ( t τ ) m α f ( m + 1 ) ( τ ) d τ ,   m > α 1 .
The frequency domain fractional-order controller design methods are generally based on the following design specifications [12]:
  • Phase margin φm and gain crossover frequency ωcg:
  • Iso-damping property;
  • High-frequency noise rejection;
  • Good output disturbance rejection; and
  • Steady-state error cancellation;
Noting with Hp(s) the transfer function of the process and with Hc(s) the transfer function of the controller, these design specifications can be mathematically described as:
  • | H C ( j ω g c ) · H P ( j ω g c ) | = 0 d B ; arg ( H C ( j ω g c ) · H P ( j ω g c ) ) = π + ϕ m ;
  • d arg ( H C ( j ω ) · H P ( j ω ) ) d ω | ω = ω g c = 0 ;
  • | T ( j ω ) = H C ( j ω ) · H P ( j ω ) 1 + H C ( j ω ) · H P ( j ω ) | A d B , with A the desired noise attenuation for frequencies ω ω T rad/s;
  • | S ( j ω ) = 1 1 + H C ( j ω ) · H P ( j ω ) | B d B , with B the desired value of the sensitivity function for frequencies ω ω S rad/s.
Using the frequency definition of fractional order [12]:
( j ω ) α = e j π α 2 = ω α cos π α 2 + j ω α sin π α 2 ,
it is easy to imagine the complexity of the resulted inequality system.
In contrast, the proposed method offers new, simple and tuning-friendly rules for fractional order PID controllers, with guaranteed phase margin and gain crossover frequency, while the fractional order offers an excellent tradeoff between dynamic performances and stability robustness.
The paper is structured as follows. After this first, brief introductory part, the second section describes the proposed controller design method, followed by case studies for different process models, from integer order model to fractional order model. The work ends with concluding remarks.

2. The Proposed Controller Design Method

The design method is inspired by the ‘symmetrical optimum’ introduced by Kessler [19]. The plant to be controlled is assumed to be of the form:
H P ( s ) = K 0 i = 1 m ( 1 + T p i s ) j = 1 n ( 1 + T p j s ) e T m s ,
where Tpi correspond to large (compensable) time constants with respect to the sum of the ‘parasitic’ time constants Tpj and time delay Tm, i.e.,:
T p i > > j = 1 n T j + T m = T Σ .
Therefore, for the frequencies below 1 / T Σ the plant transfer function can be approximated by:
H P ( s ) = K 0 ( 1 + T Σ s ) i = 1 m ( 1 + T p i s )
and in the region of the crossover frequency is furthermore approximated by a cascade of pure integrators:
H P ( s ) K 0 ( 1 + T Σ s ) i = 1 m ( T p i s ) = K 0 i = 1 m ( T p i ) s m ( 1 + T Σ s ) = K 0 s m ( 1 + T Σ s ) .
The Kessler’s ‘symmetrical optimum’ can be expressed as follows: ‘The crossover frequency of the compensated system is ωcr = 1/(2T) and the PI(D) is adjusted such that a region with a slope of −20 dB/s is secured for one octave on the right and m octaves on the left of the crossover frequency” [20].
The resulting tuning rules for a PI controller, in the case when m = 1, are the well-known Kessler’s equations [20]:
H C ( s ) = K C · T c ( 1 + 1 T c s )   with   K C = 1 8 T Σ 2 K 0 ;   T C = 4 T Σ
obtained from the open loop: H o l ( s ) = 1 + 4 T Σ s 8 T Σ 2 s 2 ( 1 + T Σ s ) , with gain crossover frequency and open loop gain:
ω g c = 1 2 T Σ ;   k = 1 8 T Σ 2 .
In the same research Voda and Landau conclude that in the case of a pure integrator plant (l/sT1), a damped response with 43% overshoot is obtained–due to the zero (1 + 4Ts) in the closed loop–having a rise time of 3.1T and a settling time of 16.3T. The gain margin is GM >2.7 and the phase margin (PM) is 36.8°. These performances, excepting phase margin which cannot be modified, can be corrected by using a reference filter or by using a PI with the proportional part acting only on the output. The above-mentioned performance becomes unacceptable due to a large sensitivity with respect to the modification of the plant gain accompanied by an alleviation of the phase margin. This shortcoming can be much stronger if T corresponds to the sum of parasitic time constants [18]. A way for control system performance enhancement, including the value of the phase margin, is obtained by the generalization of the tuning rules in terms of [21]:
K C = 1 β β T Σ 2 K 0 ;   T C = β T Σ
with the recommended values for β constrained to the range [4,18]. The exact value of β is chosen as a result of a compromise between the imposed closed loop performances (overshoot, settling time, etc.) and the desired phase margin. The case of β = 4 is the solution presented in [21].
The gain crossover frequency and the gain of the open loop in this case are:
ω g c = 1 β T Σ ;   k = 1 β β T Σ 2 .
A more generalized form of this method and a new approach, including one more degree of freedom using fractional order derivatives are proposed in the present work.

2.1. The Generalized Optimum Method

For the plant transfer function as in (8), with an integral behavior, the ideal form of the open loop transfer function which can reject a step disturbance is:
H o l ( s ) = k s 2 .
The closed loop in this case will ensure perfect disturbance rejection in steady state, but the phase margin is 0°, the system is at stability limit, a highly oscillatory system. To correct this problem, a positive phase element is added to the open loop:
H o l ( s ) = k s 2 · T 1 s + 1 T 2 s + 1 ,   with   T 1 > T 2 .
To ensure maximum stability (maximum value of phase margin), the gain crossover frequency is imposed to be at the maximum value of the open loop phase characteristic. Analytically this can be expressed by the equations:
{ | H o l ( j ω g c ) | = 0 d B = 1 d H o l ( j ω ) d ω | ω = ω g c = 0 .
Solving this equation system yields the gain crossover frequency and the open loop gain:
ω g c = 1 T 1 · T 2   and   k = 1 T 1 · T 2 T 2 T 1
being a generalization of (12), with one more degree of freedom, ensuring better performances than the classical method.
Choosing the time constants:
T 2 = T Σ   and   T 1 = β · T Σ ,
the particular tuning rules are obtained in terms of Preitl and Precup [21] and for β = 4 the Kessler’s optimum method [20].

2.2. Fractional Order Optimum Method

The above presented method has the disadvantage of compromise between the desired closed loop performances and the desired phase margin. To eliminate this disadvantage, one more degree of freedom can be added using a fractional order correction element in Equation (13):
H o l ( s ) = k s 2 · β 2 T s α + 1 T s α + 1 , α ,
where α is the fractional order, the generalization of the classical operation of derivation and integration to orders other than integer [12]. Theoretically, this parameter can take any real, positive value. However, for the controller to have physical meaning, the interval of the fractional orders of integration and differentiation is usually limited to (0, 2) [12].
The multiplication term of the time constant is chosen β2 instead of β to avoid the square root from the controller’s tuning equations.
Considering this fractional order form of the open loop, the system Equation (14) becomes:
{ | k ( j ω g c ) 2 · β 2 T ( j ω g c ) α + 1 T ( j ω g c ) s α + 1 | = 0 d B = 1 d d ω | k ( j ω g c ) 2 · β 2 T ( j ω g c ) α + 1 T ( j ω g c ) s α + 1 | ω = ω g c = 0 .
The explicit form of the equations for the gain crossover frequency and phase conditions, using the frequency definition of fractional order, Equation (4), are:
{ K ω 2 · 1 + 2 β 2 T ω α cos α π 2 + ( β 2 T ω α ) 2 1 + 2 T ω α cos α π 2 + ( T ω α ) 2 | ω = ω g c = 1 d d ω ( π + arctan β 2 T ω α sin α π 2 1 + β 2 T ω α cos α π 2 arctan T ω α sin α π 2 1 + T ω α cos α π 2 ) | ω = ω g c = 0 .
The solution of this system, the fractional order generalization of Equation (14) is:
ω g c = ( 1 β T ) 1 α   and   k = 1 β ( 1 β T ) 2 α .
Using these equations, for any chosen value of the fractional order, at the crossover frequency the phase always reaches its maximum value:
arctan β sin α π 2 1 + β cos α π 2 arctan 1 β sin α π 2 1 + 1 β cos α π 2
or, expressed in terms of phase margin:
PM = arctan ( β 2 1 ) tan α π 2 ( β + 1 ) 2 + β tan 2 α π 2 .
With the particular values T = T Σ and α = 1 the results from (8), while for β = 2 , the Kessler’s form (Equation (9)) is obtained.
The obtained phase margin and gain crossover frequency as function of the fractional order is plotted in Figure 1 and Figure 2. For simplicity reasons is considered T = T = 1, without affecting the conclusions.
It can be observed that both gain crossover frequency and phase margin increases with the fractional order α, meaning increased stability and smaller settling time as higher the α is. For α = 1 the “classical” Kessler’s phase margin value of 36.8° and gain crossover frequency is obtained, as in Equation (10).
From the step response of the closed control loop with different values of α the corresponding overshoots can be determined, resulting the plot from Figure 3. This plot reveals that the overshoot decreases with the increasing fractional order.
The main advantage of the Kessler’s optimum method-the steady state speed error cancellation-is maintained with the fractional order system as well. Moreover, the higher the fractional order, the better the transient response of the closed loop system, as it is presented in Figure 4.
Having two more degrees of freedom-due to the parameter β and fractional order α–the controller design problem becomes an optimization problem: select the proper β and α to ensure the desired closed loop performance measures.
In Figure 5 the gain crossover frequency evolution (which is inversely proportional with the settling time) with respect to parameters β and α is presented. In a similar manner the maximum phase margin, Figure 6, or any other desired performance measure of the closed loop system can be represented. With this optimization technique the desired performances of the system can be ensured, no matter how rigorous they are.
The controller parameter tuning algorithm can be described as follows:
  • Having the process mathematical model of the form of Equation (5), the open loop transfer function form is imposed as in Equation (16) to provide zero steady-state position and velocity error.
  • Using Equations (18) and (19) the tuning parameters K, α and β for the desired values of gain crossover frequency and phase margin are computed.
  • Having the open loop in Equation (17) and the process model in Equation (5), the transfer function of the fractional order controller in one of the forms presented in [22] is obtained.
The controller obtained with the proposed method being a fractional order one, engineers are faced with the problem of implementation. Actually, the fractional-order controller itself is an infinite-dimensional linear filter due to the fractional-order differentiator. A band-limit implementation is important in practice. Finite dimensional approximation of the fractional order controller should be used in a proper range of frequency of practical interest. A possible approximation method, the most widely applicable, is the Oustaloup recursive algorithm [23].

3. Case Studies

The results obtained with the controller tuning method presented in the previous section are illustrated.

3.1. Integer Order Plant, without Zero

It is considered a plant described by the transfer function:
H P ( s ) = 1 s ( T s + 1 ) .   with   T =   1 ,
a typical process model from mechatronics. This model could be the transfer function from position to armature voltage in a DC motor. Applying the general tuning rules presented in the introductory section, results the following nonlinear inequalities system:
K P 1 + 2 K i ω g c λ cos π · λ 2 + ω g c 2 λ · K i 2 = Re H P 2 + Im H P 2 ,
K i sin ( π · λ 2 ) ω c g λ + K i cos ( π · λ 2 ) = tan ( π ϕ m + arctan ( Im H P Re H P ) ) ,
K i · λ · ω g c λ 1 · sin ( π · λ 2 ) K i 2 + 2 ω g c λ · K i · cos ( π · λ 2 ) + ω g c 2 λ = Re P | ω = ω g c · d d ω Im P | ω = ω g c Im P | ω = ω g c · d d ω Re P | ω = ω g c Re P 2 | ω = ω g c + Im P 2 | ω = ω g c ,
K P 2 · ( ω 2 λ + 2 K i · ω λ · cos ( π · λ 2 ) + K i 2 ) [ ω λ ( Re P + K P ) + K i K P · cos ( π · λ 2 ) ] 2 + [ ω λ · Im P K i K P · sin ( π · λ 2 ) ] 2 A 2 , ω ω T ,
ω 2 λ · ( Re P 2 + Im P 2 ) [ ω λ ( Re P + K P ) + K i K P · cos ( π · λ 2 ) ] 2 + [ ω λ · Im P K i K P · sin ( π · λ 2 ) ] 2 B 2 , ω ω S ,
where Kp, Ki, and λ are the tuning parameters of the fractional PI controller H C ( s ) = K P ( 1 + K i s λ ) . The resulted system can be solved by optimization routines or approximation methods. A possible solution for the present case study, obtained with the fmincon command in Matlab® (20192-academic use, MathWorks, Inc., Natick, MA, USA) is Kp = 0.97, Ki = 0.1 and λ = 2.245. With these parameters the gain crossover frequency is 0.6 (rad/s) and the phase margin is 60°.
As opposed to this complex design, following the proposed method at least the same performances can be achieved using simple, user-friendly equations, more suitable for industrial applications. The open loop structure in this case is:
H o l ( s ) = k s 2 · β 2 s α + 1 s α + 1 ,   α .
The parameters α and β can be established using the equations:
ω g c = ( 1 β ) 1 α   and   k = ( 1 β ) 2 + α α .
The maximum achievable phase margin results from Equation (19):
arctan ( β 2 1 ) tan α π 2 ( β + 1 ) 2 + β tan 2 α π 2 .
The design problem consists now in an optimization between the desired performances: finding the proper values of β and α to have maximum value of the phase margin, maximum value of gain crossover frequency. With the chosen values of β and α, the controller structure results from the open loop and the process transfer function.
In Table 1 the frequency domain performance measures are presented (phase margin, gain crossover frequency) and the control system main performance measures for a unit step input (overshoot, rise time) for different values of the parameters α and β.
Figure 7 highlights the closed loop performance improvement for different fractional order in comparison with the Kessler’s optimum method values (α = 1, β = 2), while Figure 8 deals with the frequency domain measures. These simulation results are obtained using a second order “crone” approximation of the fractional order in the frequency domain (10−2, 102) rad/s [23,24].
Comparing the results with the Kessler’s performance measures—indicated in the first row of the table, obtained for α = 1 , β = 2 , and with the blue line on Figure 7 and Figure 8—the advantages are obvious. For example for β = 2 and α = 1.5 the achieved phase margin value is 77.65°, instead of 36.87° with the classical method, the overshoot is 22.4% instead of 43.2%, while the rise time is 1.59 s instead of 2.15 s. The integer order approximation of the designed controller is:
H C ( s ) = 0.79311 ( s + 9.961 ) ( s + 1.461 ) s ( s + 9.844 ) ( s + 0.1016 ) ( s 2 + 1.113 s + 1 ) + ( s + 0.1071 ) ( s 2 + 0.3966 s + 0.161 ) s ( s + 9.844 ) ( s + 0.1016 ) ( s 2 + 1.113 s + 1 ) .
For β = 4 and α = 1.5 a phase margin of 100.78°, overshoot 12.9% is obtained, but the rise time increases to 3.35 s. The best solution for the controller parameters can be achieved by an optimization procedure based on the imposed performance values.

3.2. Integer Order Plant with Zero

Another advantage of the method consists in the possibility to apply for a large variety of process models. For a process model having a zero, which is typical for parallel connected systems, the following transfer function is considered:
H P ( s ) = 0.5 s + 1 s ( s + 1 ) .
The obtained performances are: for α = 1.3 and β = 2 the overshoot is 26.9%, rise time 1.74 s, phase margin 55.5°, gain crossover frequency 0.587 rad/s, while for α = 1.3 and β = 4 the overshoot decreases to 15.9%, phase margin increases to 83°, but the gain crossover frequency became 0.345 rad/s yielding a rise time of 4.23 s. These results are depicted in Figure 9, Figure 10, Figure 11 and Figure 12.
The integer order approximation of the designed controller for α = 1.3 and β = 2 , the considered optimum in this case, is:
H C ( s ) = 0.94585 ( s + 12.1 ) ( s + 1.881 ) s ( s + 11.83 ) ( s + 2 ) ( s + 0.08454 ) ( s 2 + 1.729 s + 1 ) + ( s + 0.2315 ) ( s 2 + 0.1845 s + 0.01243 ) s ( s + 11.83 ) ( s + 2 ) ( s + 0.08454 ) ( s 2 + 1.729 s + 1 ) .

3.3. Fractional Order Plant

The more general case considered is a fractional order plant. Such a process model is typical for effective modeling of high order control plants, for example of modeling experimental heat plant.
For:
H P ( s ) = 1 s 1.5 ( s + 1 )
is considered β = 2 and α = 1.2 in the proposed algorithm. The obtained performances are: overshoot 33%, rise time 1.74 s, phase margin 48°, gain crossover frequency 0.561 rad/s, while for α = 1.2 and β = 4 the overshoot decreases to 15%, phase margin increases to 67.9°, but the gain crossover frequency became 0.315 rad/s yielding a rise time of 4.05 s. The steady-state velocity error is zero, as it is expected. All these results are highlighted in Figure 13, Figure 14, Figure 15 and Figure 16. The integer order approximation of the designed controller in this case is:
H C ( s ) = B ( s ) A ( s ) ,
where:
B ( s ) = 793.1058 ( s + 10 ) ( s + 9.961 ) ( s + 1.585 ) ( s + 1.461 ) · ( s + 0.2512 ) ( s + 0.1071 ) ( s + 0.03981 ) ( s + 0.00631 ) ( s 2 + 0.3966 s + 0.161 ) ,
A ( s ) = s ( s + 158.5 ) ( s + 25.12 ) ( s + 9.844 ) ( s + 3.981 ) · ( s + 1 ) ( s + 0..631 ) ( s + 0.1016 ) ( s + 0.1 ) ( s 2 + 1.113 s + 1 ) .
The cost of such a good results is the implementation of a fractional order controller, instead of a simple PI or PID controller, but with the present hardware possibilities this is not a real problem.

3.4. Experimental Case Study

In order to prove the efficiency of the proposed method, an experimental case study is included. The DC motor is a versatile execution element which requires a certain degree of robustness due to varying operation conditions, load changes and other varying variables linked to it, making it a linear parameter varying system. It is explicitly chosen this example due to its simplicity in dynamics and operation. The experimental unit consists in the modular servo system designed by Inteco [25] used in the particular configuration indicated in Figure 17. The plant is composed of a tachogenerator (used to measure the rotational speed), inertia load, backlash, incremental encoder, and gearbox with output disk.
The mathematical model of the modular servo system without backlash and loads has been determined experimentally for the operating point of 100 rad/s as:
H P ( s ) = Ω ( s ) u ( s ) = k s ( T s + 1 ) = 194 s ( 0.6 s + 1 ) ,
where Ω is the angular position of the rotor and u is the input voltage.
Applying the proposed method, it is imposed:
H o l ( s ) = k s 2 · β 2 T s α + 1 T s α + 1 .
For a gain crossover frequency ωgc = 0.75 and phase margin PM = 72°, using Equations (19) and (20), the parameters α = 1.5, β = 2 and k = 0.198 are obtained. Computing the controller transfer function from the well-known equation H C ( s ) = H o l ( s ) H P ( s ) results H C ( s ) = 0.001 · ( 2.4 s 1.5 + 1 ) ( 0.6 s + 1 ) s ( 0.6 s 1.5 + 1 ) . Approximating this fractional order transfer function with the Oustaloup recursive approximation, results:
H C ( s ) = 0.793 s 5 + 9.459 s 4 + 16.270 s 3 + 7.672 s 2 + 2.505 s + 0.199 s 5 + 11.060 s 4 + 13.070 s 3 + 11.060 s 4 + s .
This controller was implemented using a specialized RT-DAC/PCI-D I/O board and the Real Time Windows Target toolbox from Matlab®. The obtained results are in accordance with the imposed values, as it is presented in Figure 18.

4. Conclusions

A novel tuning method for fractional order controllers is presented, inspired by the Kessler’s optimum method. The major advantages of the method are:
  • Ensures practically any closed loop performance measures, given the possibility to choose the most convenient solution by optimization for the tuning parameters α and β. Each solution will ensure the maximum possible value of the phase margin.
  • It is a simple method of the same complexity as the Kessler’s optimum method.
  • It can be applied for practically any type of process model, from integer order models to fractional order models, which can be approximated as in Equation (9).
All the main points of this work were verified both by numerical simulation and experimental results. The simulated case studies include process transfer functions with integer order, with zeros and fractional order transfer functions as well. As an experimental case study it was chosen as an example of great simplicity in mechatronic applications and basic loop control in the manifold of production systems: the DC motor.

Funding

This research was funded by Hungarian Academy of Science, Janos Bolyai Grant (BO/ 00313/17) and Bolyai+ Grant (ÚNKP-19-4-OE-64).

Acknowledgments

This research was supported by Hungaryan Academy of Science, Janos Bolyai Grant and the ÚNKP-19-4-OE-64 NEW NATIONAL EXCELLENCE PROGRAM OF THE MINISTRY FOR INNOVATION AND TECHNOLOGY.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Chen, Y.Q.; Petras, I.; Xue, D. Fractional order control-A tutorial. In Proceedings of the 2009 American Control Conference, St. Louis, MO, USA, 10–12 June 2009; pp. 1397–1411. [Google Scholar]
  2. Tejado, I.; Vinagre, B.M.; Traver, J.E.; Prieto-Arranz, J.; Nuevo-Gallardo, C. Back to Basics: Meaning of the Parameters of Fractional Order PID Controllers. Mathematics 2019, 7, 530. [Google Scholar] [CrossRef]
  3. Petráš, I.; Terpák, J. Fractional Calculus as a Simple Tool for Modeling and Analysis of Long Memory Process in Industry. Mathematics 2019, 7, 511. [Google Scholar] [CrossRef]
  4. Gutiérrez, R.E.; Rosário, J.M.; Machado, J.T. Fractional Order Calculus: Basic Concepts and Engineering Applications. Math. Probl. Eng. 2010, 2010, 375858. [Google Scholar] [CrossRef]
  5. Fathalla, A.R. Numerical Modeling of Fractional-Order Biological Systems. Abstr. Appl. Anal. 2013, 2013, 816803. [Google Scholar] [CrossRef]
  6. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, Singapore, 2000. [Google Scholar] [CrossRef]
  7. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; John Wiley & Sons: New York, NY, USA, 1993. [Google Scholar]
  8. Podlubny, I. Fractional-order systems and PIλDμ-controllers. IEEE Trans. Autom. Control. 1999, 44, 208–213. [Google Scholar] [CrossRef]
  9. Oustaloup, A. La Commande CRONE; Editions HERMES: Paris, France, 1991. [Google Scholar]
  10. Bruzzone, L.; Fanghella, P. Fractional-Order Control of a Micrometric Linear Axis. J. Control. Sci. Eng. 2013, 2013, 947428. [Google Scholar] [CrossRef]
  11. Bruzzone, L.; Fanghella, P. Comparison of PDD1/2 and PDu Position Controls of a Second Order Linear System. In Proceedings of the 33rd IASTED International Conference on Modelling, Identification and Control MIC 2014, Innsbruck, Austria, 17–19 February 2014; pp. 182–188. [Google Scholar]
  12. Monje, C.A.; Chen, Y.Q.; Vinagre, B.M.; Xue, D.; Feliu-Batlle, V. Fundamentals of Fractional-Order Systems; Springer: London, UK, 2010. [Google Scholar]
  13. Dulf, E.H.; Timis, D.; Muresan, C.I. Robust Fractional Order Controllers for Distributed Systems. Acta Polytech. Hung. 2017, 14, 163–176. [Google Scholar]
  14. Muresan, C.I.; Dulf, E.H.; Both, R. Vector-based tuning and experimental validation of fractional-order PI/PD controllers. Nonlinear Dyn. 2016, 84, 179–188. [Google Scholar] [CrossRef]
  15. Padula, F.; Visioli, A. Tuning rules for optimal PID and fractional-order PID controllers. J. Process Control 2011, 21, 69–81. [Google Scholar] [CrossRef]
  16. Valerio, D.; da Costa, J.S. Tuning of fractional PID controllers with Ziegler Nichols-type rules. Signal Process. 2006, 86, 2771–2784. [Google Scholar] [CrossRef]
  17. Hmed, A.B.; Amairi, M.; Aoun, M.; Hamdi, S.E. Comparative study of some fractional PI controllers for first order plus time delay systems. In Proceedings of the 2017 18th International Conference on Sciences and Techniques of Automatic Control and Computer Engineering (STA), Monastir, Tunisia, 21–23 December 2017; pp. 278–283. [Google Scholar]
  18. Garrappa, R.; Kaslik, E.; Popolizio, M. Evaluation of Fractional Integrals and Derivatives of Elementary Functions: Overview and Tutorial. Mathematics 2019, 7, 407. [Google Scholar] [CrossRef]
  19. Kessler, C. Das symmetrische Optimum. Regelungstechnik 1958, 6, 432–436. [Google Scholar]
  20. Voda, A.A.; Landau, I.D. A Method for the Auto-calibration of PID Controllers. Automatica 1995, 31, 41–53. [Google Scholar] [CrossRef]
  21. Preitl, S.; Precup, R.E. An extension of tuning relations after symmetrical optimum method for PI and PID controllers. Automatica 1999, 35, 1731–1736. [Google Scholar] [CrossRef]
  22. Ahmadi Dastjerdi, A.; Vinagre, B.M.; Chen, Y.Q.; HosseinNia, S.H. Linear fractional order controllers; A survey in the frequency domain. Annu. Rev. Control 2019, 47, 51–70. [Google Scholar] [CrossRef]
  23. Oustaloup, A.; Levron, F.; Mathieu, B.; Nanot, F.M. Frequency-band complex noninteger differentiator: Characterization and synthesis. IEEE Trans. Circuits Syst. I 2000, 47, 25–39. [Google Scholar] [CrossRef]
  24. Oustaloup, A.; Sabatier, J.; Lanusse, P.; Malti, R.; Melchior, P.; Moreau, X.; Moze, M. An overview of the crone approach in system analysis, modeling and identification, observation and control. IFAC Proc. 2008, 41, 14254–14265. [Google Scholar] [CrossRef] [Green Version]
  25. Inteco, Poland. Modular Servo System-User’s Manual. Available online: www.inteco.com.pl (accessed on 26 August 2019).
Figure 1. Variation of phase margin with respect to the fractional order α.
Figure 1. Variation of phase margin with respect to the fractional order α.
Mathematics 07 01166 g001
Figure 2. Variation of gain crossover frequency with respect to the fractional order α.
Figure 2. Variation of gain crossover frequency with respect to the fractional order α.
Mathematics 07 01166 g002
Figure 3. Overshoot variation with respect to the fractional order α.
Figure 3. Overshoot variation with respect to the fractional order α.
Mathematics 07 01166 g003
Figure 4. Transient for ramp input for different fractional order.
Figure 4. Transient for ramp input for different fractional order.
Mathematics 07 01166 g004
Figure 5. Gain crossover frequency evolution with respect to parameters β and α.
Figure 5. Gain crossover frequency evolution with respect to parameters β and α.
Mathematics 07 01166 g005
Figure 6. Phase margin evolution with respect to parameters β and α.
Figure 6. Phase margin evolution with respect to parameters β and α.
Mathematics 07 01166 g006
Figure 7. Step response of the closed loop with integer order plant for different fractional order α and β = 2.
Figure 7. Step response of the closed loop with integer order plant for different fractional order α and β = 2.
Mathematics 07 01166 g007
Figure 8. Bode plot of the control loop with integer order plant for different fractional order α and β = 2.
Figure 8. Bode plot of the control loop with integer order plant for different fractional order α and β = 2.
Mathematics 07 01166 g008
Figure 9. Step response of the closed loop with integer order plant with zero for different fractional order α and β = 2.
Figure 9. Step response of the closed loop with integer order plant with zero for different fractional order α and β = 2.
Mathematics 07 01166 g009
Figure 10. Bode plot of the control loop with integer order plant with zero for different fractional order α and β = 2.
Figure 10. Bode plot of the control loop with integer order plant with zero for different fractional order α and β = 2.
Mathematics 07 01166 g010
Figure 11. Step response of the closed loop with integer order plant with zero for different fractional order α and β = 4.
Figure 11. Step response of the closed loop with integer order plant with zero for different fractional order α and β = 4.
Mathematics 07 01166 g011
Figure 12. Bode plot of the control loop with integer order plant with zero for different fractional order α and β = 4.
Figure 12. Bode plot of the control loop with integer order plant with zero for different fractional order α and β = 4.
Mathematics 07 01166 g012
Figure 13. Step response of the closed loop with fractional order plant for different fractional order α and β = 2.
Figure 13. Step response of the closed loop with fractional order plant for different fractional order α and β = 2.
Mathematics 07 01166 g013
Figure 14. Bode plot of the control loop with fractional order plant for different fractional order α and β = 2.
Figure 14. Bode plot of the control loop with fractional order plant for different fractional order α and β = 2.
Mathematics 07 01166 g014
Figure 15. Step response of the closed loop with fractional order plant for different fractional order α and β = 4.
Figure 15. Step response of the closed loop with fractional order plant for different fractional order α and β = 4.
Mathematics 07 01166 g015
Figure 16. Bode plot of the control loop with fractional order plant for different fractional order α and β = 4.
Figure 16. Bode plot of the control loop with fractional order plant for different fractional order α and β = 4.
Mathematics 07 01166 g016
Figure 17. The experimental unit: the modular servo system.
Figure 17. The experimental unit: the modular servo system.
Mathematics 07 01166 g017
Figure 18. Response of the DC motor for a square input signal.
Figure 18. Response of the DC motor for a square input signal.
Mathematics 07 01166 g018
Table 1. Performance measures for different β and α values.
Table 1. Performance measures for different β and α values.
αGain Crossover Frequency (rad/s)Phase Margin
(°)
Overshoot
(%)
Rise Time
(s)
β = 210.5036.8743.22.15
1.10.5342.6338.22.01
1.20.5649.2933.01.95
1.30.5857.0828.01.79
1.40.6066.3823.41.69
1.50.6377.6522.41.59
β = 310.3350.9024.93.39
1.10.3658.0018.53.26
1.20.4065.8515.43.19
1.30.4274.6215.02.62
1.40.4584.5014.92.43
1.50.4895.7614.82.22
β = 410.2556.3017.34.91
1.10.2863.6414.54.70
1.20.3171.5913.34.58
1.30.3480.3113.14.08
1.40.3789.9613.03.74
1.50.39100.7812.93.35

Share and Cite

MDPI and ACS Style

Dulf, E.-H. Simplified Fractional Order Controller Design Algorithm. Mathematics 2019, 7, 1166. https://doi.org/10.3390/math7121166

AMA Style

Dulf E-H. Simplified Fractional Order Controller Design Algorithm. Mathematics. 2019; 7(12):1166. https://doi.org/10.3390/math7121166

Chicago/Turabian Style

Dulf, Eva-Henrietta. 2019. "Simplified Fractional Order Controller Design Algorithm" Mathematics 7, no. 12: 1166. https://doi.org/10.3390/math7121166

APA Style

Dulf, E. -H. (2019). Simplified Fractional Order Controller Design Algorithm. Mathematics, 7(12), 1166. https://doi.org/10.3390/math7121166

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