Optimal V -Plane Robust Stabilization Method for Interval Uncertain Fractional Order PID Control Systems

: Robust stability is a major concern for real-world control applications. Realization of optimal robust stability requires a stabilization scheme, which ensures that the control system is stable and presents robust performance for a predeﬁned range of system perturbations. This study presented an optimal robust stabilization approach for closed-loop fractional order proportional integral derivative (FOPID) control systems with interval parametric uncertainty and uncertain time delay. This stabilization approach, which is carried out in a v -plane, relies on the placement of the minimum angle system pole to a predeﬁned target angle within the stability region of the ﬁrst Riemann sheet. For this purpose, tuning of FOPID controller coefﬁcients was performed to minimize a root angle error that is deﬁned as the squared difference of minimum angle root of interval characteristic polynomials and the desired target angle within the stability region of the v -plane. To solve this optimization problem, a particle swarm optimization (PSO) algorithm was implemented. Findings of the study reveal that tuning of the target angle can also be used to improve the robust control performance of interval uncertain FOPID control systems. Illustrative examples demonstrated the effectiveness of the proposed v-domain, optimal, robust stabilization of FOPID control systems.


Introduction
Since fractional order system modeling provides more realistic mathematical modeling of real systems, fields of control theory have been extended toward fractional order system modeling and fractional order controller design. Particularly, fractional order systems have paved the way for improved frequency-domain modeling and controller designs. Numerous works have illustrated enhancement of control performance, which is possible by the use of fractional calculus: these research studies have reported robust performance improvements of fractional order controllers in case of parametric perturbations [1][2][3][4][5][6]. Several works have also revealed that fractional order control can enhance the disturbance rejection performance of closed-loop control systems [7][8][9][10]. One of the underlying basic reasons for such performance improvements is the fact that fractional orders of controller functions provide more tuning options that enable one to achieve more restricted frequencydomain design objectives.
System stability is the main concern in control system design, and it is a fundamental stage of control system design efforts. For this reason, stabilization of fractional order control systems has been addressed in several perspectives: stabilization of the systems according to system pole placements [11][12][13][14][15][16][17], closed-loop system stabilization based on stability boundary locus (SBL) analyses [18,19], stabilization based on value set analysis and zero exclusion principle [20,21]. The linear matrix inequalities (LMIs) method has been widely used for the stability checking of fractional order systems [22][23][24]. Graphical stabilization methods have been proposed for robust stabilization of interval plant functions [25] and time-delay system models [26]. Tuning of FOPID controllers for fractional order uncertain systems was shown by using a robust D-stability method [27].
Due to parametric uncertainty of real systems and limitations of mathematical modeling to represent real-world systems, the robust stabilization problem of interval systems becomes one of the essential topics of practical control system design. Ensuring system stability for the interval uncertainty ranges of mathematical models improves the realworld performance of the designed control systems and is an important outcome of robust controller design methodologies. As is well-known, the stability region for integer order systems is the left half-side of the complex s-plane, hence stability analysis based on the checking of system pole placements in the complex s-plane was referred to as the left half-plane (LHP) stability analysis. Previously, stabilization based on the placement of system poles was carried out by residing all roots of interval characteristic polynomials into the stability region of the first Riemann sheet [12][13][14]16,17]. In these works, stability checking for fractional order characteristic polynomials was facilitated by applying a polynomial order expansion approach that is employed by s = v m conformal mapping. Thus, fractional order characteristic polynomials are transformed to the expanded order integer order polynomials in the v-plane. However, this transformation maps the LHP stability region in the s-plane into the region with an angle range of (π/2m, π/m] within the first Riemann sheet of the v-plane [11,12,15,[28][29][30]. This conformal mapping preserves the argument and magnitude relations of system poles, and thus it allows the implications of the edge theorem for the robust stability analysis of interval uncertain fractional order systems and consideration of edge and vertex polynomials in stability checking within the complex v-plane [15]. An application of the edge theorem was shown for the robust stability analysis of fractional order interval systems [14,15]. Later, the edge theorem was discussed for the general interval uncertainty cases, and the necessary and sufficient conditions to check system stability were proposed in [31]. Kharitonov-theorem-based approaches were also utilized for the robust stability checking of fractional order systems [32,33]. Principal characteristic equations and stability relation between the s-plane and v-plane were discussed in [30].
Even though simulation results indicate a satisfactory control performance, the realworld control systems should be designed robust against a certain level of system perturbations and parameter fluctuations. Therefore, robust stability analysis and robust system design have a major significance in the practical fractional order control system design process.
The so-called "fractal robustness" has become a central motivation for studies that have addressed fractional order control system design [34][35][36][37][38][39][40][41]. Fractional order dynamics has been widely harnessed to utilize iso-damping property in robust controller design processes [36][37][38][39][40][41]. The optimal fractional order systems have been designed [42][43][44]. A group of works has proposed graphical approaches for the robust stabilization of interval fractional order systems [25][26][27]. Stabilization of fractional order controllers based on minimum angle characteristic root placement within the first Riemann sheet of the v-plane has been shown [16,17,45]. Optimal FOPID controller design has been demonstrated in the v-plane [46,47]. Effects of pole placements on a fractional order proportional integral (FOPI) controller performance were discussed in [48]. The main advantage of this approach is that the placement of the minimum argument root of an interval characteristic polynomial set can be easily carried out to guarantee the robust stability of interval control systems for a predefined uncertainty range of system parameters, and the uncertainty ranges can be given as a design specification to be achieved. Utilization of these advantages for robust FOPID control system design was our main motivation in the current study. Due to complications in analytically deriving a general rule to govern the placement of system poles, we implemented a meta-heuristic optimization method and developed a computeraided design scheme, which can perform the placement of the minimum angle system pole of a given interval uncertain FOPID control system model to a target angle within the stability region. In addition to the robust stabilization of the interval uncertain system by optimally tuning FOPID controller coefficients, the target angle specification allows the adjustment of the control performance of resulting robust stable FOPID control systems. This property can be useful to improve the control performance besides the robustness of the designed FOPID control system.
In the literature, several heuristic optimization methods, for example, evolutionary methods [42] and particle swarm methods [43,44], have been applied for the tuning of fractional order controllers. The main advantage of heuristic optimization methods comes from the fact that their implications on real-world problems are relatively straightforward, and they can provide satisfactory solutions in the case of complicated optimization problems. For the stabilization of fractional order systems in the v-domain, Alagoz demonstrated that the brute force search could deal with the stabilization of nominal closed-loop control systems according to a minimum angle target root specification in the Riemann sheet [17]. Afterward, Tufenkci et al. implemented a genetic algorithm to stabilize fractional order PID control systems for one-pole fractional order nominal plant models and demonstrated the effectiveness of the v-domain design approach for FOPID controller stabilization [45]. These results became a motivation for extending this approach to the optimal robust stabilization of interval uncertain control systems in the v-domain by using a metaheuristic optimization method. Accordingly, the present study illustrated a design framework that employs computational intelligence for optimal robust stabilization of FOPID control system design in the v-domain. Authors referred to this approach as v-domain design because design efforts were conducted in the v-plane after applying the conformal mapping s = v m to fractional order transfer functions. It can be an alternative to well-known design domains (e.g., time domain, frequency domain and s-domain) when designing fractional order systems.
Essentially, the current study presented a methodology for the optimal robust stabilization of interval uncertain closed-loop FOPID control systems according to the minimum angle system pole placement strategy in the v-plane. For this purpose, the problem of robust stabilization of an interval uncertain fractional order plant function was considered. For a solution to this problem, the purposed approach aimed to stabilize a set of sampled characteristic polynomials of parametric uncertain FOPID control systems that are defined by an interval uncertainty box of the plant parameters. The five design coefficients of the FOPID controller were tuned by the particle swarm optimization (PSO) algorithm in order to ensure robust stability of the system by considering the minimum angle root of all characteristic polynomials that were uniformly sampled from the interval uncertainty box. To guarantee robust stability of the uncertain system, the PSO algorithm places the minimum angle characteristic root of whole interval systems on a target angle line within the stability region. Furthermore, we observed that the robust control performance of the designed FOPID control systems for a given parameter uncertainty range specification can be improved by changing the target angle line specification. This promises an advantage of tuning the control performance of the FOPID control systems by tuning a single design parameter, namely the target angle, and provides optimal robust stabilization of the interval uncertain control systems.

Stability Analysis of Fractional Order Systems According to Minimum Angle Root Placement
Nominal fractional order transfer functions are written in the form of where the denominator polynomial coefficients a i and numerator polynomial coefficients b i are real numbers [28,29]. The fractional orders are denoted by α i ∈ R, α i > α i−1 ≥ 0 (i = 0, 1, 2, 3 . . . , n) and φ i ∈ R, φ i > φ i−1 ≥ 0 (i = 0, 1, 2, 3 . . . , m) [21]. The characteristic polynomial of this function class is expressed as In order to calculate roots of fractional order characteristic polynomials, s = v m mapping is commonly used [12,13,29]. Thus, stability analysis of fractional order models is carried out by examining the root locus of expanded degree integer order characteristic polynomials within the first Riemann sheet [11][12][13][14]16,29]. By applying the s = v m transformation to Equation (2), one obtains expanded degree integer order characteristic polynomials in the v-plane as follows: Roots of the expanded degree integer order characteristic polynomials are calculated, and only the roots that lie in the first Riemann sheet are evaluated for stability analysis [11][12][13][14]16,29]. The rest of the roots, which are not in the first Riemann sheet, are assumed not to be meaningful for system stability. The first Riemann sheet is a portion of the complex plane that is confined within an argument range of (−π/m, π/m) [11,29]. The first Riemann sheet splits into two sections that are stable and unstable root regions by the stability line with the angle of π/2m. Since the first Riemann sheet is symmetrical with respect to the real axis in the complex plane, one can perform stability analysis in the positive angle side of the first Riemann sheet that is in the angle range of [0, π/m]. As a consequence, a fractional order system model is said to be stable when all characteristic roots lie in the range of [π/2m, π/m]. Figure 1 depicts the stability region in the v-plane.
where the denominator polynomial coefficients i a and numerator polynomial co cients i b are real numbers [28,29]. The fractional orders are denoted by i ∈ α ) [21]. The charac istic polynomial of this function class is expressed as  [12,13,29]. Thus, stability analysis of fractional order mode carried out by examining the root locus of expanded degree integer order character polynomials within the first Riemann sheet [11][12][13][14]16,29]. By applying the m v s = tr formation to Equation (2), one obtains expanded degree integer order characteristic p nomials in the v-plane as follows: Roots of the expanded degree integer order characteristic polynomials are calcula and only the roots that lie in the first Riemann sheet are evaluated for stability anal [11][12][13][14]16,29]. The rest of the roots, which are not in the first Riemann sheet, are assum not to be meaningful for system stability. The first Riemann sheet is a portion of the c plex plane that is confined within an argument range of ) / , / ( m m π π − [11,29]. The Riemann sheet splits into two sections that are stable and unstable root regions by stability line with the angle of m 2 / π . Since the first Riemann sheet is symmetrical w respect to the real axis in the complex plane, one can perform stability analysis in the p itive angle side of the first Riemann sheet that is in the angle range of ] . As a c sequence, a fractional order system model is said to be stable when all characteristic r lie in the range of . Figure 1 depicts the stability region in the v-plane. To perform stability analysis, the root set of expanded degree integer order cha teristic polynomials is written as [14,16] The expanded degree integer order characteristic polynomials are said to robust stable polynomials if the minimum argument condition given is satisfied [13,17,29]. This analysis can be extended to interval certain fractional order system models. Interval uncertainty of the characteristic poly mials is expressed by [14][15][16] [ ] where the interval uncertainty of gain coefficients and interval uncertainty of fractio orders are defined by the parameter ranges of ] , [  To perform stability analysis, the root set of expanded degree integer order characteristic polynomials is written as [14,16] The expanded degree integer order characteristic polynomials ∆ m (v) are said to be robust stable polynomials if the minimum argument condition given by min{Arg(R)} > π/2m is satisfied [13,17,29]. This analysis can be extended to interval uncertain fractional order system models. Interval uncertainty of the characteristic polynomials is expressed by [14][15][16] where the interval uncertainty of gain coefficients and interval uncertainty of fractional orders are defined by the parameter ranges of Fractal Fract. 2021, 5, 3

of 21
These uncertainty ranges of parameters define an n-dimensional interval uncertainty box within the parameter space of the system model, and it is represented by a hyperrectangle A that is expressed by the Cartesian product of intervals [14,15], A point of hyperrectangle A is expressed by a design parameter vector, which is defined in the parameter space as Each vector p represents a possible characteristic polynomial of an interval uncertain system and if all polynomials from hyperrectangle A are stable, the interval system is said to be robustly stable. In other words, when min{Arg(R i )} > π/2m is satisfied for all roots of the polynomials of hyperrectangle A that are written by the interval uncertain system is robustly stable [13,15,17].

Brief Introduction of PSO Algorithms
Meta-heuristic optimization methods have an important milestone in computational intelligence practice because of their three prominent advantages: (i) they can be easily implemented to solve complicated optimization problems, (ii) the random search nature of these methods presents a potential of a new solution at each run of algorithms, and (iii) they can be easily applied for the optimization of simulation models or even real systems in offline or online manners. For these reasons, meta-heuristic methods have become popular optimization tools for computer-aided solutions to complex engineering problems.
Due to difficulties in obtaining an analytical solution of fractional calculus, meta-heuristic optimization methods have been widely preferred for optimization problems of fractional order system design [17,[42][43][44]. The PSO algorithm could find satisfactory solutions for the parameter tuning problems of fractional order systems [42][43][44]. The effectiveness of the PSO algorithm in the tuning of control systems mainly originates from the fact that the search space of control systems generally has smoothly changing characteristics that are suitable for the traveling of particles for the search of an optimal point.
The PSO algorithm is an essential swarm-based search algorithm, which can imitate the collective intelligence in animal swarms. The algorithm of PSO benefits from the collective behavior of search agents, namely particles. The particles can move under the effects of interactions with other particles and form a global coherence for the search of a particle swarm [49,50]. Interactions and coherence at particle motions result in a swarm intelligence that can be similar to collective motions of animal swarms in their natural environments. To resemble natural swarm intelligence in the search of an optimal solution, the PSO algorithm presumes that particles move collectively in a multidimensional search space of an optimization problem. Coherence and interactions among the particles of PSO form dynamics in swarm motion, which tends to gather particles to minima of search spaces. Thus, the positions of particles in the search space represent a candidate solution for the PSO algorithm. The motion of a particle is expressed by particle positions x n (t) and the particle velocity v n (t) parameters. The PSO models coherent particle motions by the following formulation [49][50][51].
x n [t where x L,n (t) is the personal best position, and x G,n (t) is the global best position. Parameters c 1 and c 2 are personal learning coefficients of x L,n (t) and global learning coefficient of x G,n (t), respectively. These parameters are used to adjust the strength of local and global interactions among particles. Parameters r 1 and r 2 are pseudo-random numbers. These random parameters enable particles to have a stochastic nature in their motions and allow finding new solutions in the case of subsequent running of the algorithm [49][50][51]. The parameter w is the weight coefficient for particle inertias. It is used to decrease the inertia of particles during the optimization process, and a damping rate ξ is used to update the weight coefficient at each iteration as .ξ [29]. While particles are moving in the search space, the local best position and the global best position are updated at each iteration according to particle scores. These scores are provided by a predefined objective function. The objective function of the optimization evaluates particle positions and decides whether or not to be a valuable point for the solution for an optimization problem [49][50][51].

Methodology for Robust Stabilization of FOPID Control System
Adjusting places of minimum argument roots within the stability region was demonstrated in [17,45]. Optimal FOPI and FOPID controller tuning in the v-domain was illustrated in [46][47][48]. The current study extended this approach to solve the optimal robust stabilization problem of FOPID control systems. Especially, it extended the works in [17] and [45] to stabilization of interval uncertain FOPID control systems, and it demonstrates the effectiveness of the v-domain stabilization approach for optimal stabilization of uncertain FOPID control systems with interval coefficients and interval time-delay uncertainty. Figure 2 illustrates a general block diagram of closed-loop FOPID control systems. The FOPID controller is commonly expressed in the form of where C(s) denotes the transfer function of the FOPID controller. The nominal plant function is assumed to be a time-delay, one-pole fractional order plant, which is commonly written in a general form as where the parameter L stands for the time delay parameter of the plant, and the parameter α is the fractional order of plant function. The parametric interval uncertainty model of this plant function can be expressed in the form of where interval uncertainty of the system is defined with parameter alteration ranges of moving in the search space, the local best position and the global best position ar at each iteration according to particle scores. These scores are provided by a p objective function. The objective function of the optimization evaluates particle and decides whether or not to be a valuable point for the solution for an op problem [49][50][51].

Methodology for Robust Stabilization of FOPID Control System
Adjusting places of minimum argument roots within the stability region w strated in [17,45]. Optimal FOPI and FOPID controller tuning in the v-domain trated in [46][47][48]. The current study extended this approach to solve the optim stabilization problem of FOPID control systems. Especially, it extended the wo and [45] to stabilization of interval uncertain FOPID control systems, and it dem the effectiveness of the v-domain stabilization approach for optimal stabilization tain FOPID control systems with interval coefficients and interval time-delay un Figure 2 illustrates a general block diagram of closed-loop FOPID control sys FOPID controller is commonly expressed in the form of denotes the transfer function of the FOPID controller. The nom function is assumed to be a time-delay, one-pole fractional order plant, whi monly written in a general form as where interval uncertainty of the system is defined with parameter alteration ] , [ . These interval ranges define mensional interval uncertainty box, that is, a hyperrectangle A in the paramet the plant model. Here, each point of hyperrectangle A is expressed by a param ] [  The transfer function of this interval system can be written as To obtain a constant coefficient polynomial, the term e −Ls is substituted with the first order Padé approximation, which is written by where the relation between coefficients and time delay parameter is [θ, θ] = [L/2 L/2]. The characteristic polynomial of the closed-loop control system is the denominator polynomial of the transfer function T(s), and it can be expressed for this interval system model as One applies s = v m mapping and obtains the expanded degree integer order characteristic polynomials as Based on Equation (8), the characteristic root set of the interval system can be written by considering FOPID controller coefficients To achieve placement of the minimum angle root of all polynomials in A near to a target minimum angle root position, the angle error is defined as where φ T is the target angle for minimum angle root placement inside the stability region of the first Riemann sheet [17,45]. The term min|arg(v r )| refers to the minimum angle of roots in the set R i . The squared angle error (SAE) is expressed as The optimization problem that minimizes the SAE is written by Specifying the minimum target angle in the range of π/2m < φ T ≤ π/m ensures the placement of the minimum angle root within the stability region of the v-plane, and thus stabilizes the resulting FOPID control system according to a target angle φ T . The φ T can be configured to obtain improved robust control performance. This optimization problem is solved by using the PSO method according to the flow chart in Figure 3, and the results are discussed in the following section. ract. 2021, 5, x FOR PEER REVIEW 8 of 23 Figure 3. The block diagram describes the implication of the particle swarm optimization (PSO) algorithm for optimal robust stabilization. SAE: squared angle error.

Illustrative Design Examples
This section presents illustrative examples to demonstrate the robust stabilization of FOPID control system models. For the initial PSO configuration, the search ranges of FOPID controller coefficients were set to . The population size of particles was configured to 30 particles. The personal learning coefficient ( 1 c ) was set to 2.0, and the global learning coefficient ( 2 c ) was set to 2.0. The inertia weight ( w ) and damping rate ( ξ ) were 1 and 0.99, respectively [50,51].
The iteration number of PSO is limited to 200. To perform time-domain simulations of control systems, step() function of the FOTF toolbox was used in MATLAB [52]. Figure 4 depicts a design process that includes optimization, simulation and validation stages.
In the case of 10 v s = mapping, the expanded degree integer order characteristic polynomial with interval uncertainty is written by

Illustrative Design Examples
This section presents illustrative examples to demonstrate the robust stabilization of FOPID control system models. For the initial PSO configuration, the search ranges of FOPID controller coefficients were set to k p ∈ [0, 7], k i ∈ [0, 7], k d ∈ [0, 7], λ ∈ [0, 2] and µ ∈ [0, 2]. The population size of particles was configured to 30 particles. The personal learning coefficient (c 1 ) was set to 2.0, and the global learning coefficient (c 2 ) was set to 2.0. The inertia weight (w) and damping rate (ξ) were 1 and 0.99, respectively [50,51]. The iteration number of PSO is limited to 200. To perform time-domain simulations of control systems, step() function of the FOTF toolbox was used in MATLAB [52]. Figure 4 depicts a design process that includes optimization, simulation and validation stages.

Illustrative Design Examples
This section presents illustrative examples to demonstrate the robust stabilization of FOPID control system models. For the initial PSO configuration, the search ranges of FOPID controller coefficients were set to . The population size of particles was configured to 30 particles. The personal learning coefficient ( 1 c ) was set to 2.0, and the global learning coefficient ( 2 c ) was set to 2.0. The inertia weight ( w ) and damping rate ( ξ ) were 1 and 0.99, respectively [50,51].
The iteration number of PSO is limited to 200. To perform time-domain simulations of control systems, step() function of the FOTF toolbox was used in MATLAB [52]. Figure 4 depicts a design process that includes optimization, simulation and validation stages.

Example 1. Suppose that one tries to figure out a stabilizing FOPID controller for an interval uncertain plant function with an interval delay,
In the case of 10 v s = mapping, the expanded degree integer order characteristic polynomial with interval uncertainty is written by We used s = v 10 mapping, and the target angle for the stabilization was set to 3π 40 within the first Riemann sheet. The characteristic polynomial of the interval system was expressed as In the case of s = v 10 mapping, the expanded degree integer order characteristic polynomial with interval uncertainty is written by For 200 iterations of the optimization process, a decrease of the minimum angle error is shown in Figure 5. The figure clearly demonstrates that the PSO algorithm can significantly decrease SAE during the optimization process. This result confirms that the problem can be solved by the PSO algorithm. For 200 iterations of the optimization process, a decrease of the minimum angle error is shown in Figure 5. The figure clearly demonstrates that the PSO algorithm can significantly decrease SAE during the optimization process. This result confirms that the problem can be solved by the PSO algorithm.  Figure 6 shows placements of all roots of sampled characteristic polynomials from the interval uncertainty box, where the characteristic roots are indicated by red asterisks and dots. The figure confirms that the minimum argument root of the interval system was placed on the target angle that is indicated by the green line within the stability region. To illustrate the stabilization of a conventional controller, a PI controller was also optimally stabilized by using the proposed method. The PSO algorithm yielded the design parameters 1.0372 =  After completion of optimization, the PSO algorithm yielded the design parameters k p = 4.4739, k i = 2.6179, k d = 1.3096, λ = 0.8629 and µ = 0.0 for the stabilizing FOPID controller. The SAE is almost zero (at the level of 10 −31 ). Figure 6 shows placements of all roots of sampled characteristic polynomials from the interval uncertainty box, where the characteristic roots are indicated by red asterisks and dots. The figure confirms that the minimum argument root of the interval system was placed on the target angle that is indicated by the green line within the stability region. To illustrate the stabilization of a conventional controller, a PI controller was also optimally stabilized by using the proposed method. The PSO algorithm yielded the design parameters k p = 1.0372 and k i = 5.4914 for the stabilizing PI controller. It is useful to demonstrate step responses of the stabilized control systems to validate robust system stability in transient analysis. Figure 7 shows step responses of the stabilizing FOPID and PI controllers for the nominal plant function that was formed by the parameter values from the mid-points of uncertainty intervals.
To evaluate unit step performances of systems at boundaries of interval parameters, Figure 8 shows step responses of the FOPID and the PI controllers for the stabilization of 16 plant functions that were sampled from vertices of the uncertainty box. These plant functions are listed in Table 1 by considering all possible combinations of interval boundaries. These functions express vertices of the interval uncertainty box, which are the most probable to cause instability of the system. The figure shows that responses of both FOPID and PI control systems for all vertex plant functions are stable, and this indicates that the controllers can control the interval uncertain plants even if plant parameters change to the boundaries of uncertainty intervals. The high overshoot responses in the figure are mainly related to higher magnitude minimum angle roots of the interval system.
The settling of FOPDI control systems is faster than the settling of PI control systems for the control of 16 plant functions because most responses of the PI control system have prolonging ripples. This is an indication of improved robust stabilization of the FOPID controller compared to the PI controller.
The SAE is almost zero (at the level of 31 10 − ). Figure 6 shows placements of all roots of sampled characteristic polynomials from the interval uncertainty box, where the characteristic roots are indicated by red asterisks and dots. The figure confirms that the minimum argument root of the interval system was placed on the target angle that is indicated by the green line within the stability region. To illustrate the stabilization of a conventional controller, a PI controller was also optimally stabilized by using the proposed method. The PSO algorithm yielded the design parameters 1.0372 = p k and for the stabilizing PI controller. It is useful to demonstrate step responses of the stabilized control systems to validate robust system stability in transient analysis. Figure 7 shows step responses of the stabilizing FOPID and PI controllers for the nominal plant function that was formed by the parameter values from the mid-points of uncertainty intervals.        (1 − Q(t)) 2 dt). In the case of s = v 10 mapping, the expanded degree integer order characteristic polynomials can be written by To observe the effect of target angle configuration, PSO optimization was performed for three different target angle specifications within the stability region of the v-plane. These are φ T1 = 4π 60 (for a placement close to the stability boundary by taking d = 1/3), φ T2 = 3π 40 (for a placement in the middle of the stability region by taking d = 1/2) and φ T3 = 7π 80 (for a placement close to the upper bound of the first Riemann sheet by taking d = 3/4). Specification of the target angle can be performed by using the formula

Vertex Plant Forms Vertex Plant Functions
is a partitioning factor of the stability region by a target angle line. For 200 iterations of the optimization process, decreases in SAEs for each target angle specifications are shown in Figure 9. The figures clearly show that the PSO algorithm can significantly decrease SAE during the optimization process. Table 3 shows the stabilizing FOPD controller designs that were obtained after the completion of the PSO algorithm.     Figure 12 shows the step responses of 16 vertex functions of the uncertainty box. For the target angle specification φ T1 , step performances of the control system for the vertex functions are more similar than the results of other target angle specifications. This indicates that the FOPID controller for φ T1 is more robust than other FOPID designs due to less variability among step responses that were obtained for the vertex plant functions. Therefore, we observed in this example that the control performance of φ T1 can be more robust than other target angle specifications against parametric perturbation of the plant function. The FOPID controller design for φ T1 can be more suitable for real control applications.
igure 12 shows the step responses of 16 vertex functions of the uncertainty box. For rget angle specification 1 T  , step performances of the control system for the vertex ions are more similar than the results of other target angle specifications. This indithat the FOPID controller for 1 T  is more robust than other FOPID designs due to ariability among step responses that were obtained for the vertex plant functions. fore, we observed in this example that the control performance of 1 T  can be more t than other target angle specifications against parametric perturbation of the plant ion. The FOPID controller design for 1 T  can be more suitable for real control applis. n Figure 12b,c, the variability of control performance is higher between vertices of terval box, and it indicates more deterioration of the control performance while paters are changing inside the interval box. The step response MSE performances in s 4 and 5 validate these effects.    In Figure 12b,c, the variability of control performance is higher between vertices of the interval box, and it indicates more deterioration of the control performance while parameters are changing inside the interval box. The step response MSE performances in Tables 4 and 5 validate these effects.
In order to evaluate disturbance rejection performances of the robust stabilized FOPID control system designs for the target angle specifications φ T1 ,φ T2 and φ T3 , step disturbance responses of the systems are shown in Figure 13. In these control systems, a set-point prefilter function was used to overcome overshoots, and it improves settling to the step input before applying the additive step disturbance to the input of the plant function. The figure indicates satisfactory disturbance rejection performances of the optimally stabilized robust FOPID controllers in this example.  G 7 ). optimization task. Since classical PSO is a low-complexity search algorithm, it is preferable to reduce the overall computational load of the optimization task.

Conclusions
This study demonstrated an optimal robust stabilization scheme for FOPID control systems. This method is based on the placement of the minimum angle system pole over a stabilization line in the v-domain. To this end, the robust stabilization problem was in- In these example designs, convergence rates of the PSO algorithm were measured by considering the decrease in error (E) per iteration until the convergence of the algorithm and the convergence rates were calculated at levels of 10 −8 error/iteration. However, the computation time can take four to five hours for a standard PC (hardware configuration is I5-3230M 2.60 GHz, 4GB RAM) because it requires the computation of minimum angle roots for all sampled systems (p vectors) from the uncertainty box of the interval plant functions (a 0 ∈ [a 0 , a 0 ], b 0 ∈ [b 0 , b 0 ], b 1 ∈ [b 1 , b 1 ],L ∈ [L, L]) for each particle of PSO as shown in the pseudo code in Appendix A. Searching the minimum angle root for each sample model of the uncertainty box significantly increases the computational load in the optimization task. Since classical PSO is a low-complexity search algorithm, it is preferable to reduce the overall computational load of the optimization task.

Conclusions
This study demonstrated an optimal robust stabilization scheme for FOPID control systems. This method is based on the placement of the minimum angle system pole over a stabilization line in the v-domain. To this end, the robust stabilization problem was introduced as an optimization problem that aims to place the minimum angle root of interval characteristic polynomials on a target angle line that was specified within the stability region of the first Riemann sheet. This study extended stabilization efforts to the optimal robust stabilization of interval uncertain closed-loop FOPID control systems with an interval time delay. It implemented PSO for computational intelligence. Some remarks of the study can be summarized as follows: • This study demonstrated a v-domain design scheme that is straightforward for optimal robust stabilization of fractional order control systems and the development of computer-aided-design tools.

•
The PSO algorithm can find optimal FOPID controller coefficients that lead to minimum angle roots of interval systems placed on the desired angle line within the stability region of the first Riemann sheet. Thus, the proposed approach can ensure the stability of FOPID control systems in uncertainty ranges of plant parameters. The angle of this line is configured by a target angle specification. • Target angle specifications can improve the robust control performance of FOPID control systems. We observed that the target angle specifications could be utilized to reduce the sensitivity of the time response performances of control systems to the variation of plant parameters within the uncertainty ranges. Thus, the method can achieve optimal stabilization of FOPID control system designs.

•
The proposed stabilization scheme may have implications for variable time-delay systems such as the control and robust stabilization problems of the networked control systems.
In summary, this research work addressed optimal robustly stabilizing FOPID controller design for variable time-delay, interval uncertain fractional order plant functions. A straightforward computer-aided design solution for this type of control problems is not frequent in the literature. In future works, the presented v-domain design methodology can be applied to other fractional order model types, and the associated optimization problems can be solved by other search algorithms.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
Pseudo code for the employment of the PSO algorithm in this optimization problem: Receive controller parameters (kd,kp,ki,λµ) as a candidate particle from the PSO algorithm. % This part uniformly samples the hyperrectangle A and calculates minimum angle roots of p sample designs and stores them.
Find minimum argument roots among the roots with positive angle values and store them. End End End End % This part selects the minimum argument root of the whole hyperrectangle A and calculates SAE error accordingly. Find the minimum angle value among the stored minimum angle roots. Calculate squared angle error (SAE) of the minimum angle value according to Equation (20). Return SAE to the PSO algorithm for the next evaluation.