Robust Stability Analysis of Unstable Second Order Plus Time-Delay (SOPTD) Plant by Fractional-Order Proportional Integral (FOPI) Controllers

: This study investigates the robust stability analysis of an unstable second order plus time-delay (SOPTD) plant by using Fractional-Order Proportional Integral (FOPI) controllers. We assume that there are simultaneous uncertainties in gain, time-constants, and time-delay of the plant. At ﬁrst, a graphical method is provided for a robust stability analysis of the closed-loop system. Then, a robust stability checking function is introduced to facilitate the robust stability analysis. Additionally, new bounds are presented to reduce the computational burden for the robust stability analysis. Finally, two examples are provided to show the correctness of the proposed method.


Introduction
Generally, many real-world systems can be modeled by time-delay transfer functions because the dead time is inevitable and appears in real systems [1]. In practical applications, it was shown that the delay may also contain an interval uncertainty [1][2][3]. Moreover, it has been shown that fractional-order controllers provide more flexibility compared with integerorder controllers [4]. This can be attributed to the extra parameters within fractional-order controllers [5,6]. Therefore, much attention has been paid to design robust fractional-order controllers for fractional-order nonlinear time-delayed systems [7,8] and fractional-order linear time-delayed systems [9,10].
Stability analysis is fundamental to any control system. Hence, some methods have been proposed to analyze the stability of LTI fractional-order systems in [11][12][13][14][15][16]. Lately, the value set-based approach has attracted many researchers to check the robust stability of LTI systems. In [17,18], the value set approach has been provided for robust stability analysis of LTI fractional-order systems having time delay. The works of [19,20] have proven that the value set of fractional-order systems of incommensurate and commensurate orders possessing uncertain coefficients is a parpolygon in the complex plane. Inspired by the results presented in [17,19,20], some value set-based approaches have been presented to inspect the robust stability of fractional-order plants by linear fractional-order controllers in [21,22].
In [23,24], Fractional-Order Proportional Integral Derivative (FOPID) controllers were designed to improve the robustness property of first order in addition to dead time systems. Then, in [25,26], some auxiliary functions have been presented for robust stability analysis of interval fractional order systems having time-delay. Moreover, in [27], a graphical tuning method has been proposed to determine the stabilizing regions of FOPI controllers for fractional systems with time-delays by benefiting from the D-decomposition technique. Moreover, in [28,29], new methods were introduced to design controllers for systems suffering from time-delay.
This study develops the robust stability analysis of an unstable Second Order Plus timedelay (SOPTD) plant by FOPI controllers. In [30][31][32], some methods have been proposed to design various controllers. However, the methods presented in [30][31][32] cannot ensure the robust stability of the closed-loop system in the presence of plant uncertainties in time-delay, gain, and time constants. This issue has prompted the authors of this paper to present a simple method to analyze the robust stability of a closed-loop system consisting of unstable SOPTD plants and FOPI controllers. The contributions and novelties of this paper are mainly summarized as follows: • Robust stability analysis of a closed-loop control system consisting of unstable SOPTD planst and FOPI controllers; • Introducing a frequency range for reducing computational burden; • Presenting a robust stability checking function for robust stability analysis of the closed-loop system.
The remainder of this study is organized as follows: Section 2 reviews preliminaries. Section 3 presents the main results. In Section 4, two examples are provided for demonstration. Finally, Section 5 concludes the paper.

Background and Preliminaries
Among most popular definitions proposed for fractional-order differentiation, Caputo fractional-order derivative is mainly employed and is defined as follows ( [33]): in which γ (n − 1 < γ < n), and n is an integer number. In (1), Γ is the Euler's Gamma function. a and t are the terminals of the integral. Based on the explanations presented in [33], the Laplace transform of the Caputo's fractional-order derivative can be given by the following: for the zero initial conditions can be simplified as follows.
Hence, the transfer function of an FOPI controller can be obtained as follows.
Moreover, the unstable SOPTD plant is considered as follows: in which K is the gain lying in the interval 0 < K − ≤ K ≤ K + ; T 1 and T 2 are time constants within the intervals 0 The characteristic function for the FOPI controller C(s) in (4) and the plant P(s) in (5) can be considered as ∆(s) in (6): (6) in which the following is the case.
Based on the above explanations, the following design problem is formulated. Problem 1: Develop efficient methods to investigate the robust stability of the closedloop system with characteristic function ∆(s) in (6).
Our main results are based on the following principle called the Zero Exclusion Principle.
Zero Exclusion Principle [17]: The closed-loop control system can be robustly stabilized by the FOPI controller C(s) in (4) if and only if ∆(s) contains one stable member and 0 / ∈ ∆(jω) for ω ≥ 0.

Main Results
By the zero exclusion principle, the position relationship between the origin and the value set of ∆(jω) has to be investigated for ω ≥ 0. Therefore, in Theorem 1, a frequency range is obtained instead of [0, ∞) for reducing the computational cost. Theorem 1. The origin is not contained in the value set of ∆(jω) for the frequency interval ω ∈ (ω max , ∞), where ω max is obtained as follows.
It is apparent that Based on the zero exclusion principle, it is important to check whether or not 0 ∈ ∆(jω) in (6) for ω ≥ 0. Therefore, the next theorem presents a graphical method to check the distance between the origin and the value set of ∆(jω) in (6) for 0 ≤ ω ≤ ω max . (6) for 0 ≤ ω ≤ ω max , if the value sets of −∆ 1 (jω) and ∆ 2 (jω) do not have any overlap in the complex plane.
Now, by using Theorem 2 and the zero exclusion principle, the following lemma presents a graphical procedure to check the robust stability of ∆(s) in (6).
Generally, the graphical method presented in Lemma 1 for robust stability analysis may be difficult, since we have to plot value sets ∆ b 1 (ω) in (12) and ∆ b 2 (ω) in (13) at each frequency. Hence, in the next theorem, a robust stability checking function is offered in order to easily inspect the overlap between ∆ b 1 (ω) in (12) and ∆ b 2 (ω) in (13).

Theorem 3.
Assume that a nominal characteristic function of ∆(s) in (6) as∆(s) is stable. Then, the characteristic function ∆(s) in (6) is robust and stable, if the inequality D(ω) > 0 in (14) holds for ω ∈ [0, ω max ]: 14) in which the following is the case: where N is an arbitrary number and N ∈ N ≥3 .
Proof. Figure 1 depicts the value set of K(K p s λ + K i )e −Ls and the vertices P h (s) and P h (s) for s = jω, ω ∈ [0, 2π L + −L − ). As it is seen from Figure 1, the value set of ∆ 2 (jω) is located inside of H g (jω). If it is proven that in the complex plane, the value sets of −∆ 1 (jω) and H g (jω) do not have any overlap, then based on Lemma 1, it can be concluded that this system is robust stable.
It is obvious that the triangle inequality holds for two consecutive vertices of −∆ 1 (jω) and any vertices of H g (jω), then no vertices of H g (jω) intersect −∆ 1 (jω). Moreover, if the inequality D N h,g (jω) > 0 holds, then no vertices of −∆ 1 (jω) intersect H g (jω). Therefore, in the complex plane, the value sets of −∆ 1 (jω) and H g (jω) do not have any overlap for , it is apparent that if inequality D 2 (ω) is held, then 0 / ∈ ∆(jω). Accordingly, based on Lemma 1 and the zero exclusion principle, the proof can be simply completed. Figure 1. The value set of ∆ 2 (jω) and the vertices P h (jω) and P h (jω).

Remark 1.
The robust stability of ∆(jω) can be checked by performing the following steps:

1.
Check the stability of one arbitrary member of ∆(s) as∆(s). If∆(s) is not robust stable, then clearly ∆(s) is not robust stable. Otherwise, continue the following step; 2.
Calculate the value of ω max presented in Theorem 1; 3.

Illustrative Examples
In this Section, two examples show the use of the obtained results.

Example 1.
Consider of the following unstable plant discussed in [32].
Let us consider the plant P(s) in (16) as an uncertain plant in (17).
In the following, it is shown that FOPI controller C(s) = 2 + 0.02/s 0.1 can robustly stabilize the closed-loop system. The approach proposed in [12] is employed to check the stability of∆(s) in (18). Based on [12] and Figure 2, it can be concluded that∆(s) is stable, since, in the complex plane, ψ(jω) in (19) does not encircle the origin.
Moreover from Theorem 1, one can obtain ω max = 4.7361 rad/s. Figure 3 shows the values of the auxiliary function D(ω) in (14). It is clear that the inequality D(ω) > 0 holds for ω ∈ [0, ω max ]. From Figure 3 and Theorem 3, the closed-loop system is robust stable. As this example shows, the work in [32] cannot guarantee the robust stability of such systems. On the other hand, the results of the present paper can be used for robust stability analysis of unstable processes having uncertain time delay by FOPI controllers. Moreover, Theorem 1 helps to reduce the computational burden by presenting the upper frequency bound ω max .
and an FOPI controller as follows.
Again, the method presented in [12] is used to check the stability of∆(s) in (22). Figure 4 shows the values of ψ(jω) in (23). Based on [12] and Figure 4, it can be deduced that∆(s) is stable since, in the complex plane, ψ(jω) in (23) does not encircle the origin.

Conclusions
This study was dedicated to check the robust stability of systems consisting of an unstable SOPTD plant and FOPI controllers. Theorem 1 was presented to reduce the computational cost for the robust stability analysis, since, instead of frequency [0, ∞), a robust stability analysis can be performed in ω ∈ [0, ω max ]. Moreover, Lemma 1 was offered to present a graphical method for the robust stability analysis. In Theorem 3, an auxiliary function was provided to investigate robust stability. Moreover, auxiliary function D(ω) facilitates the process of robust stability analysis. this is because the sign of this function can determine the robust stability of the system. Moreover, two examples were presented to show the effectiveness and simplicity of the results. Future works may include presenting the stabilizing region of FOPI controllers for such systems having uncertain time delays.