Simpliﬁed Fractional Order Controller Design Algorithm

: 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 speciﬁcations 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 ﬁrst 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, o ﬀ ering 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.

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 PDD 1/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]: where Γ(n) = ∞ 0 t n−1 e −t dt 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]: The Grünwald-Letnikov's definition of the fractional-order derivative [12]: The frequency domain fractional-order controller design methods are generally based on the following design specifications [12]:
Good output disturbance rejection; and 5.
Steady-state error cancellation; Noting with H p (s) the transfer function of the process and with H c (s) the transfer function of the controller, these design specifications can be mathematically described as: 1.
S( jω) = 1 1+H C ( jω)·H P ( jω) ≤ BdB, with B the desired value of the sensitivity function for frequencies ω ≤ ω S rad/s. Using the frequency definition of fractional order [12]: 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.

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: where T pi correspond to large (compensable) time constants with respect to the sum of the 'parasitic' time constants T pj and time delay T m , i.e.,: Therefore, for the frequencies below 1/T Σ the plant transfer function can be approximated by: and in the region of the crossover frequency is furthermore approximated by a cascade of pure integrators: 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]: obtained from the open loop: , with gain crossover frequency and open loop gain: In the same research Voda and Landau conclude that in the case of a pure integrator plant (l/sT 1 ), a damped response with 43% overshoot is obtained-due to the zero (1 + 4T s) 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]: 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: 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.

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: 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: 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: Solving this equation system yields the gain crossover frequency and the open loop gain: being a generalization of (12), with one more degree of freedom, ensuring better performances than the classical method. Choosing the time constants: the particular tuning rules are obtained in terms of Preitl and Precup [21] and for β = 4 the Kessler's optimum method [20].

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): 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: The explicit form of the equations for the gain crossover frequency and phase conditions, using the frequency definition of fractional order, Equation (4), are: The solution of this system, the fractional order generalization of Equation (14) is: Using these equations, for any chosen value of the fractional order, at the crossover frequency the phase always reaches its maximum value: or, expressed in terms of phase margin: Mathematics 2019, 7, 1166 6 of 21 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 Figures 1 and 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. 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. 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 d The controller parameter tuning algorithm can be described as follows:

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

2.
Using Equations (18) and (19) the tuning parameters K, α and β for the desired values of gain crossover frequency and phase margin are computed.

3.
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].

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

Integer Order Plant, without Zero
It is considered a plant described by the transfer function: . 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: where K p , K i , 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 K p = 0.97, K i = 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: The parameters α and β can be established using the equations: The maximum achievable phase margin results from Equation (19): 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 β.

Integer Order Plant with Zero
Another advantage of the method consists in process models. For a process model having a zero, w the following transfer function is considered: 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.

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: The          The integer order approximation of the designed controller for α = 1.3 and β = 2, the considered optimum in this case, is:

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: 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 Figures 13-16. The integer order approximation of the designed controller in this case is:          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.

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

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.