# Design of Fuzzy Fractional PD + I Controllers Tuned by a Genetic Algorithm

Authors: Isabel S. Jesus; Ramiro S. Barbosa
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The fractional-order concepts are a useful tool to describe several physical phenomena, and nowadays they are widely used in the field of automatic control. A genetic algorithm (GA) is a search process for finding approximate solutions in optimization problems. The GA provides further flexibility and robustness that are unique for signal process. In this paper we consider the development of an optimal fuzzy fractional PD + I controller in which the parameters are tuned by a GA. The performance of the proposed fuzzy fractional control is illustrated through some application examples.


---

## Body

## 1. Introduction

Fractional calculus (FC) is a generalization of integration and differentiation to a noninteger orderα
∈
C, with the fundamental operator being a
D
t
α, where a and t are the limits of the operation [1, 2]. The FC concepts constitute a useful tool to describe several physical phenomena, such as heat, flow, electricity, magnetism, mechanics, or fluid dynamics. Presently, the FC theory is applied in almost all areas of science and engineering, with its ability being recognized in bettering the modelling and control of many dynamical systems. In fact, during the last years FC has been used increasingly to model the constitutive behavior of materials and physical systems exhibiting hereditary and memory properties. This is the main advantage of fractional-order derivatives in comparison with classical integer-order models, where these effects are simply neglected.In this paper we investigate several control strategies/structures based on fuzzy fractional-order algorithms. The fractional-order PID controller (PI
α
D
β controller) involves an integrator of order α
∈
ℜ
+ and a differentiator of order β
∈
ℜ
+. It was demonstrated the good performance of this type of controller in comparison with the conventional PID algorithms. Recently, there have been a lot of researches in the application of fuzzy PID control [3–9]. The fuzzy method offers a systematic procedure to design controllers for many kinds of systems that often leads to a better performance than that of the conventional PID controller. It is a methodology of intelligent control that mimics human thinking and reacting by using a multivalent fuzzy logic and elements of artificial intelligence.It was proved that the use of the fuzzy fractional controllers improved the results for many kinds of systems, since it gives additional flexibility to the design. In this line of thought many applications of this type of controllers were developed in the last few years. For example, in [3, 4], the authors proved the effectiveness of fuzzy fractional PD and PID controllers in terms of their digital implementation and robustness. In [10], an intelligent robust fractional surface sliding mode control for a nonlinear system is studied. In [6], a fractional-order fuzzy PID controller was proposed and compared to classical PID, fuzzy PID, and even PI
λ
D
μ controllers. Many other applications can be found in [11–16].A genetic algorithm (GA) is a search technique based on the natural selection process. The GA is a particular class of evolutionary algorithms that use techniques inspired by evolutionary biology such as inheritance, mutation, natural selection, and crossover, established by Darwin’s theory of evolution [17–19]. The GA is used in the field of robotics, strategy planning, nonlinear dynamical systems, data analysis, art, evolving pictures, music, and many others in the real world applications [17–20]. The GA provides a unique flexibility and robustness for process optimization. Due to this reason, during the last years many control applications used the GA in order to find better results. For example, in [6], the authors optimized their system with GA while minimizing several integral error indices along with the control signal as the objective function. Hu et al.  [21] developed a methodology for the systematic design of fuzzy PID controllers based on genetic optimization, where they proved that the proposed system always provides the best control performance. In [22], the GA was used for tuning the PI controller for load frequency control. Herrera and Lozano showed in [23] the benefits derived from the synergy between evolutionary algorithms and fuzzy logic systems.Bearing these ideas in mind, the paper is organized as follows. Section2 gives the fundamentals of fractional-order control systems. Section 3 presents the control and optimization strategies. Section 4 gives some simulations results assessing the effectiveness of the proposed methodology. Finally, Section 5 draws the main conclusions.


## 2. Fractional-Order Control Systems

Fractional-order control systems are characterized by differential equations that have, in the dynamical system and/or in the control algorithm, an integral and/or a derivative of fractional order [24, 25]. Due to the fact that these operators are defined by irrational continuous transfer functions, in the Laplace domain, or infinite dimensional discrete transfer functions, in the Z domain, we often encounter evaluation problems in the simulations. Therefore, when analyzing fractional-order systems, we usually adopt continuous or discrete integer-order approximations of fractional-order operators [26–28]. The following two subsections provide a background for the remaining of the paper by giving the fundamental aspects of the FC and the discrete integer-order approximations of fractional-order operators used in this study.


### 2.1. Fundamentals of Fractional Calculus

The mathematical definition of a fractional-order derivative and integral (so-calleddifferintegral) has been the subject of several different approaches [1, 2]. One commonly used definition for the fractional-order derivative is given by the Riemann-Liouville definition (α
>
0):
(1)
a
D
t
α
f
(
t
)
=
1
Γ
(
n
-
α
)
d
n
d
t
n
∫
a
t
f
(
τ
)
(
t
-
τ
)
α
-
n
+
1
d
τ
,
n
-
1
<
α
<
n
,
where n is integer, f
(
t
) is the applied function, and Γ
(
x
) is the Gamma function of x. Another widely used definition is given by the Grünwald-Letnikov approach (α
∈
ℜ
):(2a)
a
D
t
α
f
(
t
)
=
lim
⁡
h
→
0
1
h
α
∑
k
=
0
[
(
t
-
a
)
/
h
]
(
-
1
)
k
(
α
k
)
f
(
t
-
k
h
)
,
(2b)
(
α
k
)
=
Γ
(
α
+
1
)
Γ
(
k
+
1
)
Γ
(
α
-
k
+
1
)
, where h is the time increment and [x] means the integer part of x.Another definition is given by the Caputo approach:(3)
a
C
D
t
α
f
(
t
)
=
1
Γ
(
n
-
α
)
∫
a
t
‍
f
n
(
τ
)
(
t
-
τ
)
α
+
1
-
n
d
τ
,
n
-
1
<
α
<
n
.The “memory” effect of these operators is demonstrated by (1)–(3), where the convolution integral in (1) and (3) and the infinite series in (2a) and (2b) reveal the unlimited memory of these operators, ideal for modelling hereditary and memory properties in physical systems and materials.An alternative definition to (1)–(3), which reveals useful for the analysis of fractional-order control systems, is given by the Laplace transform method. Considering vanishing initial conditions, the fractionaldifferintegration is defined in the Laplace domain, F
(
s
)
=
L
{
f
(
t
)
}, as
(4)
L
{
D
a
t
α
f
(
t
)
}
=
s
α
F
(
s
)
,
α
∈
ℜ
.The fractional operator can be more easily interpreted in the frequency domain. In fact, the open-loop Bode diagrams of amplitude and phase ofs
α have correspondingly a slope of 20α dB/dec and a constant phase positioned at α
π
/
2 rad over the entire frequency domain.


### 2.2. Approximations of Fractional-Order Operators

In this study we adopt discrete integer-order approximations to the fundamental elements
α
(
α
∈
ℜ
) of a fractional-order control (FOC) strategy. The usual approach for obtaining discrete equivalents of continuous operators of type s
α adopts the Euler, Tustin, and Al-Alaoui generating functions [26, 29, 30].It is well known that rational-type approximations frequently converge faster than polynomial-type approximations and have a wider domain of convergence in the complex domain [29]. Thus, using the Euler operator w
(
z
-
1
)
=
(
1
-
z
-
1
)
/
T
c and performing a power series expansion of [
w
(
z
-
1
)
]
α
=
[
(
1
-
z
-
1
)
/
T
c
]
α give the discretization formula corresponding to the Grünwald-Letnikov Definition (2a) and (2b):
(5)
D
α
(
z
-
1
)
=
(
1
-
z
-
1
T
c
)
α
=
∑
k
=
0
∞
‍
(
1
T
c
)
α
(
-
1
)
k
(
α
k
)
z
-
k
=
∑
k
=
0
∞
‍
h
α
(
k
)
z
-
k
,
where T
c is the sampling period and h
α
(
k
) is the impulse response sequence.A rational fraction-type approximation can be obtained through a Padé approximation to the impulse response sequenceh
α
(
k
), yielding the discrete transfer function:
(6)
H
(
z
-
1
)
=
b
0
+
b
1
z
-
1
+
⋯
+
b
m
z
-
m
1
+
a
1
z
-
1
+
⋯
+
a
n
z
-
n
=
∑
k
=
0
∞
‍
h
(
k
)
z
-
k
,
where m
≤
n and the coefficients a
k and b
k are determined by fitting the first m
+
n
+
1 values of h
α
(
k
) into the impulse response h
(
k
) of the desired approximation H
(
z
-
1
). Thus, we obtain an approximation that matches the desired impulse response h
α
(
k
) for the first m
+
n
+
1 values of k [26]. Note that the above Padé approximation is obtained by considering the Euler operator but the determination process will be exactly the same for other types of discretization schemes.


## 2.1. Fundamentals of Fractional Calculus

The mathematical definition of a fractional-order derivative and integral (so-calleddifferintegral) has been the subject of several different approaches [1, 2]. One commonly used definition for the fractional-order derivative is given by the Riemann-Liouville definition (α
>
0):
(1)
a
D
t
α
f
(
t
)
=
1
Γ
(
n
-
α
)
d
n
d
t
n
∫
a
t
f
(
τ
)
(
t
-
τ
)
α
-
n
+
1
d
τ
,
n
-
1
<
α
<
n
,
where n is integer, f
(
t
) is the applied function, and Γ
(
x
) is the Gamma function of x. Another widely used definition is given by the Grünwald-Letnikov approach (α
∈
ℜ
):(2a)
a
D
t
α
f
(
t
)
=
lim
⁡
h
→
0
1
h
α
∑
k
=
0
[
(
t
-
a
)
/
h
]
(
-
1
)
k
(
α
k
)
f
(
t
-
k
h
)
,
(2b)
(
α
k
)
=
Γ
(
α
+
1
)
Γ
(
k
+
1
)
Γ
(
α
-
k
+
1
)
, where h is the time increment and [x] means the integer part of x.Another definition is given by the Caputo approach:(3)
a
C
D
t
α
f
(
t
)
=
1
Γ
(
n
-
α
)
∫
a
t
‍
f
n
(
τ
)
(
t
-
τ
)
α
+
1
-
n
d
τ
,
n
-
1
<
α
<
n
.The “memory” effect of these operators is demonstrated by (1)–(3), where the convolution integral in (1) and (3) and the infinite series in (2a) and (2b) reveal the unlimited memory of these operators, ideal for modelling hereditary and memory properties in physical systems and materials.An alternative definition to (1)–(3), which reveals useful for the analysis of fractional-order control systems, is given by the Laplace transform method. Considering vanishing initial conditions, the fractionaldifferintegration is defined in the Laplace domain, F
(
s
)
=
L
{
f
(
t
)
}, as
(4)
L
{
D
a
t
α
f
(
t
)
}
=
s
α
F
(
s
)
,
α
∈
ℜ
.The fractional operator can be more easily interpreted in the frequency domain. In fact, the open-loop Bode diagrams of amplitude and phase ofs
α have correspondingly a slope of 20α dB/dec and a constant phase positioned at α
π
/
2 rad over the entire frequency domain.


## 2.2. Approximations of Fractional-Order Operators

In this study we adopt discrete integer-order approximations to the fundamental elements
α
(
α
∈
ℜ
) of a fractional-order control (FOC) strategy. The usual approach for obtaining discrete equivalents of continuous operators of type s
α adopts the Euler, Tustin, and Al-Alaoui generating functions [26, 29, 30].It is well known that rational-type approximations frequently converge faster than polynomial-type approximations and have a wider domain of convergence in the complex domain [29]. Thus, using the Euler operator w
(
z
-
1
)
=
(
1
-
z
-
1
)
/
T
c and performing a power series expansion of [
w
(
z
-
1
)
]
α
=
[
(
1
-
z
-
1
)
/
T
c
]
α give the discretization formula corresponding to the Grünwald-Letnikov Definition (2a) and (2b):
(5)
D
α
(
z
-
1
)
=
(
1
-
z
-
1
T
c
)
α
=
∑
k
=
0
∞
‍
(
1
T
c
)
α
(
-
1
)
k
(
α
k
)
z
-
k
=
∑
k
=
0
∞
‍
h
α
(
k
)
z
-
k
,
where T
c is the sampling period and h
α
(
k
) is the impulse response sequence.A rational fraction-type approximation can be obtained through a Padé approximation to the impulse response sequenceh
α
(
k
), yielding the discrete transfer function:
(6)
H
(
z
-
1
)
=
b
0
+
b
1
z
-
1
+
⋯
+
b
m
z
-
m
1
+
a
1
z
-
1
+
⋯
+
a
n
z
-
n
=
∑
k
=
0
∞
‍
h
(
k
)
z
-
k
,
where m
≤
n and the coefficients a
k and b
k are determined by fitting the first m
+
n
+
1 values of h
α
(
k
) into the impulse response h
(
k
) of the desired approximation H
(
z
-
1
). Thus, we obtain an approximation that matches the desired impulse response h
α
(
k
) for the first m
+
n
+
1 values of k [26]. Note that the above Padé approximation is obtained by considering the Euler operator but the determination process will be exactly the same for other types of discretization schemes.


## 3. Control and Optimization Strategies

### 3.1. Fractional PID Control

The generalized PID controller,G
c
(
s
), has a transfer function of the form [28]
(7)
G
c
(
s
)
=
U
(
s
)
E
(
s
)
=
K
p
+
K
i
s
α
+
K
d
s
β
,
α
,
β
>
0
,
where α and β are the orders of the fractional integrator and differentiator, respectively. The parameters K
p, K
i, and K
d are correspondingly the proportional, integral, and derivative gains of the controller. Clearly, taking (
α
,
β
)
=
{
(
1,1
)
,
(
1,0
)
,
(
0,1
)
,
(
0,0
)
} we get the classical {
PID
,
PI
,
PD
,
P
} controllers, respectively [24, 31, 32]. Other PID controllers are possible, namely, PD
β controller, PI
α controller, PID
β controller, and so on. The fractional-order controller is more flexible and gives the possibility of adjusting more carefully the closed-loop system characteristics [2, 33].In the time domain thePI
α
D
β is represented by
(8)
u
(
t
)
=
K
p
e
(
t
)
+
K
i
D
0
t
-
α
e
(
t
)
+
K
d
D
0
t
β
e
(
t
)
.The fractional-order differential operators in (8) are implemented using the approximations (5) and (6), yielding the discrete transfer function:
(9)
G
c
(
z
)
=
U
(
z
)
E
(
z
)
=
K
p
+
K
i
H
i
(
z
-
1
)
+
K
d
H
d
(
z
-
1
)
,
where H
i
(
z
-
1
) and H
d
(
z
-
1
) are the fraction-type approximations to fractional-order integral and derivative, respectively.


### 3.2. Fuzzy Fractional PD + I Control

Fuzzy control emerged on the foundations of Zadeh’s fuzzy set theory [3, 4, 9]. This kind of control is based on the ability of a human being to find solutions for particular problematic situations. It is well known from our experience that humans have the ability to simultaneously process a large amount of information and make effective decisions, although neither input information nor consequent actions is precisely defined. Through multivalent fuzzy logic, linguistic expressions in antecedent and consequent parts of IF-THEN rules describing the operator’s actions can be efficaciously converted into a fully structured control algorithm.The fuzzy logic controllers are not dependent on accurate mathematical models, which are one of the most important advantages in its use, particularly in applications where systems are difficult to model or contain significant nonlinearities.In the system of Figure1, we apply a fuzzy logic control (FLC) for the PD
β actions and the integral of the error is added to the output in order to find a fuzzy PD
β
+
I controller [3, 34]. This kind of configuration eliminates the steady state error due to the integer integrative action. The block diagram of Figure 2 illustrates the configuration of the proposed fuzzy controller.Figure 1
Block diagram of the fuzzy control system.Figure 2
FuzzyPD
β
+
I controller.In this controller, the control actions are the errore, the fractional derivative of e, and the integral of e. The U represents the controller output. Also, the controller has four gains to be tuned, K
e, K
ie, K
ce corresponding to the inputs and K
u to the output.The control actionU is generally a nonlinear function of error E, fractional change of error CE, and integral of error IE:
(10)
U
(
k
)
=
[
f
(
E
,
CE
)
+
IE
]
K
u
=
[
f
(
K
e
e
(
k
)
,
K
ce
D
β
e
(
k
)
)
+
K
ie
I
e
(
k
)
]
K
u
,
where D
β is the discrete fractional derivative implemented as rational fraction approximation (6) using the Euler scheme (5); the integral of error is calculated by rectangular integration:
(11)
I
(
z
-
1
)
=
T
c
1
-
z
-
1
.In fact, we can adopt an integral action of fractional order,I
α, yielding a fuzzy fractional PD
β
+
I
α controller [3]. However, in this work we consider only values of α
=
1.To further illustrate the performance of the fuzzyPD
β
+
I a saturation nonlinearity is included in the closed-loop system of Figure 1 and inserted in series with the output of the fuzzy controller. The saturation element is defined as
(12)
n
(
u
)
=
{
u
,
|
u
|
<
δ
,
δ
sign
⁡
(
u
)
,
|
u
|
≥
δ
,
where u and n are, respectively, the input and the output of the saturation block and sign
⁡
(
u
) is the signum function.Here we give an emphasis of the proposed FLC presented in Figure2. The basic structure for FLC is illustrated in Figure 3 [35].Figure 3
Structure for fuzzy logic controller.The fuzzy rule base, which reflects the collected knowledge about how a particular control problem must be treated, is one of the main components of a fuzzy controller. The other parts of the controller perform make up the tasks necessary for the controller to be efficient.For the fuzzyPD
β
+
I controller illustrated in Figure 2, the rule base can be constructed in the following form (see Table 1):ifE is NM and CE is NS, then v is NL, where NL, NM, NS, ZR, PS, PM, and PL are linguistic values representing “negative large,” “negative medium,” “negative small,” “zero,” “positive small,” “positive medium,” and “positive large,” respectively. E is the error, CE is the fractional derivative of error, and v is the output of the fuzzy PD
β controller. The membership functions for the premises and consequents of the rules are shown in Figure 4.Table 1
Fuzzy control rules.
CE
E
NL
NM
NS
ZR
PS
PM
PL
NL
NL
NL
NL
NL
NM
NS
ZR
NM
NL
NL
NL
NM
NS
ZR
PS
NS
NL
NL
NM
NS
ZR
PS
PM
ZR
NL
NM
NS
ZR
PS
PM
PL
PS
NM
NS
ZR
PS
PM
PL
PL
PM
NS
ZR
PS
PM
PL
PL
PL
PL
ZR
PS
PM
PL
PL
PL
PLFigure 4
Membership functions forE, CE, and v.With two inputs and one output the input-output mapping of the fuzzy logic controller is described by a nonlinear surface, as presented in Figure5.Figure 5
Control surface.The fuzzy controller will be adjusted by changing the parameter values ofK
e, K
ce, K
ie, and K
u. The fuzzy inference mechanism operates by using the product to combine the conjunctions in the premise of the rules and in the representation of the fuzzy implication. For the defuzzification process we use the centroid method.


### 3.3. Genetic Optimization

The evolutionary computing was introduced in the 60s by I. Rechenberg, and the GA was invented by John Holland who published a book in 1975 about this subject.These algorithms begin with a set of solutions, represented by chromosomes, called population (P). Initially, a population is generated randomly. Solutions from one population are taken (parents) and used to form a new P. This is motivated by the hope that the new P will be better than the old one. Individuals are then selected to form new solutions according to their fitness; therefore, the more suitable they are the more chances they have to reproduce. This is repeated until some condition is satisfied. Figure 6 presents the block diagram representative of the methodology used in the GA.Figure 6
Block diagram of a GA.In these algorithms the crossover (C) and the mutation (M) operators are the most important parts. The C is a recombination operator that combines subparts of two parent chromosomes to produce offspring that contain some parts of both parents genetic material. The simplest way to do it is to choose some random C point, copy everything before this point from the first parent, and then to copy everything after the C point from the other parent. There are other ways to make C; namely, we can choose more C points. The most used way of encoding is a binary string; however, there are many other ways of encoding, such as to encode directly through real numbers.The selection of a better encoding technique depends on the problem we have to solve. TheM operation randomly changes the offspring resulting from C. This procedure intended to prevent falling of all solutions in the P into a local optimum. In case of binary encoding, we can switch a few randomly chosen bits from 1 to 0 or from 0 to 1.Another important concept in GA is the Elitism. The Elitism strategy (ES) was introduced by Kenneth De Jong in 1975 and is an addition to many selection methods that forces the GA to retain some number of the best individuals at each generation (G). With this tool, such individuals can be lost if they are not selected to reproduce or if they are destroyed by C or M. Many researchers have found that ES improves significantly the GA’s performance [17–20].The advantage of GA is in its parallelism. GA is travelling in a search space using more individuals than other methods. However, GA also has disadvantages, namely, the computational cost, because many times these algorithms are slower than other methodologies.In this work we propose a fuzzy fractionalPD
β
+
I controller, where the gains will be tuned through the application of a GA, in order to achieve a superior control performance of the control system of Figure 1. The optimization fitness function corresponds to the minimization of the integral time absolute error (ITAE) criterion that measures the response error as defined as [20]
(13)
J
(
K
e
,
K
ce
,
K
ie
,
K
u
)
=
∫
0
∞
‍
t
|
r
(
t
)
-
c
(
t
)
|
d
t
,
where (
K
e, K
ce, K
ie, K
u) are the PD
β
+
I controller parameters to be optimized.During the minimization of the fitness function in GA, the chromosomes (potential solutions) that lead to unstable responses make the error take very high values penalizing the fitness function. Therefore, these chromosomes are less likely to be selected for reproduction or elitism, this way being destroyed or lost in the process. This mechanism ensures that only the best chromosomes are chosen to produce new population and that the GA converges to an optimal solution with a stable closed-loop system.We can use other integral performance criteria such as the integral absolute error (IAE), the integral square error (ISE), or the integral time square error (ITSE). In the present study the ITAE criterion produced good results and is adopted in the sequel. Furthermore, the ITAE criterion enables us to study the influence of time in the error generated by the system.


## 3.1. Fractional PID Control

The generalized PID controller,G
c
(
s
), has a transfer function of the form [28]
(7)
G
c
(
s
)
=
U
(
s
)
E
(
s
)
=
K
p
+
K
i
s
α
+
K
d
s
β
,
α
,
β
>
0
,
where α and β are the orders of the fractional integrator and differentiator, respectively. The parameters K
p, K
i, and K
d are correspondingly the proportional, integral, and derivative gains of the controller. Clearly, taking (
α
,
β
)
=
{
(
1,1
)
,
(
1,0
)
,
(
0,1
)
,
(
0,0
)
} we get the classical {
PID
,
PI
,
PD
,
P
} controllers, respectively [24, 31, 32]. Other PID controllers are possible, namely, PD
β controller, PI
α controller, PID
β controller, and so on. The fractional-order controller is more flexible and gives the possibility of adjusting more carefully the closed-loop system characteristics [2, 33].In the time domain thePI
α
D
β is represented by
(8)
u
(
t
)
=
K
p
e
(
t
)
+
K
i
D
0
t
-
α
e
(
t
)
+
K
d
D
0
t
β
e
(
t
)
.The fractional-order differential operators in (8) are implemented using the approximations (5) and (6), yielding the discrete transfer function:
(9)
G
c
(
z
)
=
U
(
z
)
E
(
z
)
=
K
p
+
K
i
H
i
(
z
-
1
)
+
K
d
H
d
(
z
-
1
)
,
where H
i
(
z
-
1
) and H
d
(
z
-
1
) are the fraction-type approximations to fractional-order integral and derivative, respectively.


## 3.2. Fuzzy Fractional PD + I Control

Fuzzy control emerged on the foundations of Zadeh’s fuzzy set theory [3, 4, 9]. This kind of control is based on the ability of a human being to find solutions for particular problematic situations. It is well known from our experience that humans have the ability to simultaneously process a large amount of information and make effective decisions, although neither input information nor consequent actions is precisely defined. Through multivalent fuzzy logic, linguistic expressions in antecedent and consequent parts of IF-THEN rules describing the operator’s actions can be efficaciously converted into a fully structured control algorithm.The fuzzy logic controllers are not dependent on accurate mathematical models, which are one of the most important advantages in its use, particularly in applications where systems are difficult to model or contain significant nonlinearities.In the system of Figure1, we apply a fuzzy logic control (FLC) for the PD
β actions and the integral of the error is added to the output in order to find a fuzzy PD
β
+
I controller [3, 34]. This kind of configuration eliminates the steady state error due to the integer integrative action. The block diagram of Figure 2 illustrates the configuration of the proposed fuzzy controller.Figure 1
Block diagram of the fuzzy control system.Figure 2
FuzzyPD
β
+
I controller.In this controller, the control actions are the errore, the fractional derivative of e, and the integral of e. The U represents the controller output. Also, the controller has four gains to be tuned, K
e, K
ie, K
ce corresponding to the inputs and K
u to the output.The control actionU is generally a nonlinear function of error E, fractional change of error CE, and integral of error IE:
(10)
U
(
k
)
=
[
f
(
E
,
CE
)
+
IE
]
K
u
=
[
f
(
K
e
e
(
k
)
,
K
ce
D
β
e
(
k
)
)
+
K
ie
I
e
(
k
)
]
K
u
,
where D
β is the discrete fractional derivative implemented as rational fraction approximation (6) using the Euler scheme (5); the integral of error is calculated by rectangular integration:
(11)
I
(
z
-
1
)
=
T
c
1
-
z
-
1
.In fact, we can adopt an integral action of fractional order,I
α, yielding a fuzzy fractional PD
β
+
I
α controller [3]. However, in this work we consider only values of α
=
1.To further illustrate the performance of the fuzzyPD
β
+
I a saturation nonlinearity is included in the closed-loop system of Figure 1 and inserted in series with the output of the fuzzy controller. The saturation element is defined as
(12)
n
(
u
)
=
{
u
,
|
u
|
<
δ
,
δ
sign
⁡
(
u
)
,
|
u
|
≥
δ
,
where u and n are, respectively, the input and the output of the saturation block and sign
⁡
(
u
) is the signum function.Here we give an emphasis of the proposed FLC presented in Figure2. The basic structure for FLC is illustrated in Figure 3 [35].Figure 3
Structure for fuzzy logic controller.The fuzzy rule base, which reflects the collected knowledge about how a particular control problem must be treated, is one of the main components of a fuzzy controller. The other parts of the controller perform make up the tasks necessary for the controller to be efficient.For the fuzzyPD
β
+
I controller illustrated in Figure 2, the rule base can be constructed in the following form (see Table 1):ifE is NM and CE is NS, then v is NL, where NL, NM, NS, ZR, PS, PM, and PL are linguistic values representing “negative large,” “negative medium,” “negative small,” “zero,” “positive small,” “positive medium,” and “positive large,” respectively. E is the error, CE is the fractional derivative of error, and v is the output of the fuzzy PD
β controller. The membership functions for the premises and consequents of the rules are shown in Figure 4.Table 1
Fuzzy control rules.
CE
E
NL
NM
NS
ZR
PS
PM
PL
NL
NL
NL
NL
NL
NM
NS
ZR
NM
NL
NL
NL
NM
NS
ZR
PS
NS
NL
NL
NM
NS
ZR
PS
PM
ZR
NL
NM
NS
ZR
PS
PM
PL
PS
NM
NS
ZR
PS
PM
PL
PL
PM
NS
ZR
PS
PM
PL
PL
PL
PL
ZR
PS
PM
PL
PL
PL
PLFigure 4
Membership functions forE, CE, and v.With two inputs and one output the input-output mapping of the fuzzy logic controller is described by a nonlinear surface, as presented in Figure5.Figure 5
Control surface.The fuzzy controller will be adjusted by changing the parameter values ofK
e, K
ce, K
ie, and K
u. The fuzzy inference mechanism operates by using the product to combine the conjunctions in the premise of the rules and in the representation of the fuzzy implication. For the defuzzification process we use the centroid method.


## 3.3. Genetic Optimization

The evolutionary computing was introduced in the 60s by I. Rechenberg, and the GA was invented by John Holland who published a book in 1975 about this subject.These algorithms begin with a set of solutions, represented by chromosomes, called population (P). Initially, a population is generated randomly. Solutions from one population are taken (parents) and used to form a new P. This is motivated by the hope that the new P will be better than the old one. Individuals are then selected to form new solutions according to their fitness; therefore, the more suitable they are the more chances they have to reproduce. This is repeated until some condition is satisfied. Figure 6 presents the block diagram representative of the methodology used in the GA.Figure 6
Block diagram of a GA.In these algorithms the crossover (C) and the mutation (M) operators are the most important parts. The C is a recombination operator that combines subparts of two parent chromosomes to produce offspring that contain some parts of both parents genetic material. The simplest way to do it is to choose some random C point, copy everything before this point from the first parent, and then to copy everything after the C point from the other parent. There are other ways to make C; namely, we can choose more C points. The most used way of encoding is a binary string; however, there are many other ways of encoding, such as to encode directly through real numbers.The selection of a better encoding technique depends on the problem we have to solve. TheM operation randomly changes the offspring resulting from C. This procedure intended to prevent falling of all solutions in the P into a local optimum. In case of binary encoding, we can switch a few randomly chosen bits from 1 to 0 or from 0 to 1.Another important concept in GA is the Elitism. The Elitism strategy (ES) was introduced by Kenneth De Jong in 1975 and is an addition to many selection methods that forces the GA to retain some number of the best individuals at each generation (G). With this tool, such individuals can be lost if they are not selected to reproduce or if they are destroyed by C or M. Many researchers have found that ES improves significantly the GA’s performance [17–20].The advantage of GA is in its parallelism. GA is travelling in a search space using more individuals than other methods. However, GA also has disadvantages, namely, the computational cost, because many times these algorithms are slower than other methodologies.In this work we propose a fuzzy fractionalPD
β
+
I controller, where the gains will be tuned through the application of a GA, in order to achieve a superior control performance of the control system of Figure 1. The optimization fitness function corresponds to the minimization of the integral time absolute error (ITAE) criterion that measures the response error as defined as [20]
(13)
J
(
K
e
,
K
ce
,
K
ie
,
K
u
)
=
∫
0
∞
‍
t
|
r
(
t
)
-
c
(
t
)
|
d
t
,
where (
K
e, K
ce, K
ie, K
u) are the PD
β
+
I controller parameters to be optimized.During the minimization of the fitness function in GA, the chromosomes (potential solutions) that lead to unstable responses make the error take very high values penalizing the fitness function. Therefore, these chromosomes are less likely to be selected for reproduction or elitism, this way being destroyed or lost in the process. This mechanism ensures that only the best chromosomes are chosen to produce new population and that the GA converges to an optimal solution with a stable closed-loop system.We can use other integral performance criteria such as the integral absolute error (IAE), the integral square error (ISE), or the integral time square error (ITSE). In the present study the ITAE criterion produced good results and is adopted in the sequel. Furthermore, the ITAE criterion enables us to study the influence of time in the error generated by the system.


## 4. Simulation Results

In this section we analyze the closed-loop system of Figure1 with the fuzzy fractional PD
β
+
I controller of Figure 2. The systems used correspond to typical plants [36, 37], namely, a high order process, a high order process with a zero, and a system with a time delay. In all the experiments, the fractional-order derivative D
β is implemented using a 4th-order Padé discrete rational transfer function (
m
=
n
=
4
) of type (5). It is used as sampling period of T
c
=
0.01 s. The PD
β
+
I controller is tuned through the optimization of fitness function corresponding to the minimization of ITAE (13) using a GA. We use δ
=
15. We establish the following values for the GA parameters: population size P
=
20, crossover probability C
=
0.8, mutation probability M
=
0.05, and number of generations N
G
=
100. The gene codification adopts a decimal code. It is important to refer that a reliable execution and analysis of a GA usually require a large number of simulations to guarantee that stochastic effects are properly considered [17, 19]. Therefore, the experiments consist in executing the GA several times, for generating a good combination of controller parameters. In this study the GA is repeated 10 times and we get the best result, that is, the simulation that leads to the smaller J. We set the search space of K
e, K
ce, K
ie, and K
u
∈
[
0,5
].In the first case, we compare a fuzzy fractionalPD
β
+
I controller (
β
=
0.9
), with a fuzzy integer PD
+
I controller (
β
=
1
). Figure 7 shows the unit step responses of both controllers. The plant system G
1
(
s
) used is represented by the transfer function:
(14)
G
1
(
s
)
=
1
(
s
+
1
)
(
1
+
α
s
)
(
1
+
α
2
s
)
(
1
+
α
3
s
)
with
α
=
0.5
.Figure 7
Step responses of closed-loop system with fuzzyPD
+
I and PD
β
+
I
(
β
=
0.9
) controllers.The controller parameters, corresponding to the minimization of the ITAE index, lead to the values for the fuzzy integerPD
+
I controller: {
K
e
,
K
ce
,
K
ie
,
K
u
}
≡
{
1.0808
,
0.3408
,
0.3442
,
4.2095
}, with J
=
0.8235, and for the fuzzy fractional PD
β
+
I controller to the following values: {
K
e
,
K
ce
,
K
ie
,
K
u
}
≡
{
0.7581
,
0.3510
,
0.3038
,
4.3276
}, with J
=
0.7693. These values lead us to conclude that the fuzzy fractional-order controller produced similar results to integer one; however, the error J is smaller, as can be seen in Figure 8, where it shows the ITAE error as function of β. The graph reveals that for this process we have a lower error for a fractional value of β
=
0.9. We verify that the fractional controller is better (in terms of error J) than the integer version, only in a narrow region of 0.8
≤
β
<
1. Figure 9 illustrates the variation of FLC parameters (K
e, K
ce, K
ie, K
u) as function of the order’s derivative β, while Figure 10 shows the variation of the transient response parameters, namely, the settling time t
s, rise time t
r, peak time t
p, and overshoot o
v(%) versus β, for the closed-loop step response. The variation of FLC parameters and the transient response parameters reveals a smooth variation with β.Figure 8
ErrorJversusβ for G
1
(
s
).Figure 9
ThePID
β parameters (K
e, K
ce, K
ie, K
u) versus β for G
1
(
s
).Figure 10
Parameterst
s, t
r, t
p, o
v(%) versus β of the step responses of the closed-loop system with G
1
(
s
) and with a fuzzy PD
β
+
I controller.In a second experiment, we consider a fuzzyPD
β
+
I (β
=
0.2) controller, for process G
2
(
s
) with a right-half plane zero, represented by the transfer function:
(15)
G
2
(
s
)
=
1
-
α
s
(
s
+
1
)
3
with
α
=
5.0
.Once again, we consider for comparison the corresponding integer version (β
=
1). Figure 11 shows the unit step responses of both controllers.Figure 11
Step responses of closed-loop system with fuzzyPD
+
I and PD
β
+
I (β
=
0.2) controllers.The controller parameters, corresponding to the minimization of the ITAE index, lead to the values for the fuzzy integer controller:{
K
e
,
K
ce
,
K
ie
,
K
u
}
≡
{
0.5393
,
0.4647
,
0.3439
,
0.2486
}, with J
=
75.4509, and for the fuzzy fractional controller: {
K
e
,
K
ce
,
K
ie
,
K
u
}
≡
{
0.6061
,
0.0326
,
0.3175
,
0.2909
}, with J
=
55.9414. These values lead us to conclude that the fuzzy fractional-order controller produced better results than the integer ones, since the transient response (in particular the settling time and rise time) and the error J are smaller. Figure 12 shows the ITAE error as function of β. The graph reveals that for this process we have a lower error for β
=
0.2. Also note that the fractional controller is better (in terms of error J) than the integer one for 0.1
≤
β
<
1. Figure 13 illustrates the variation of FLC parameters (K
e, K
ce, K
ie, K
u) as function of the order’s derivative β, while Figure 14 shows the variation of t
s, t
r, t
p and o
v(%) versus β, for the closed-loop step response. In this case, the K
ce and the K
ie have the smaller value near the best case (i.e., β
=
0.2) and t
s, t
r, and t
p reveal similar values with the variation of β.Figure 12
ErrorJversusβ for G
2
(
s
).Figure 13
ThePD
β
+
I parameters (K
e, K
ce, K
ie, K
u) versusβ for G
2
(
s
).Figure 14
Parameterst
s, t
r, t
p, o
v(%) versus β of the step responses of the closed-loop system with G
2
(
s
) and with a fuzzy PD
β
+
I (0
<
β
≤
1) controller.In a third study we consider a fuzzyPD
β
+
I (β
=
0.9) controller, for process G
3
(
s
), represented by the transfer function (16), where the time delay is T
=
1 [s]:
(16)
G
3
(
s
)
=
e
-
s
T
(
1
+
α
s
)
2
with
α
=
2.0
.For comparison purposes, we consider the corresponding integer version (β
=
1). Figure 15 shows the unit step responses of both controllers.Figure 15
Step responses of closed-loop system with fuzzyPD
+
I and PD
β
+
I(β
=
0.9) controllers.The controller parameters, corresponding to the minimization of the ITAE index, lead to the values for the fuzzy integer controller:{
K
e
,
K
ce
,
K
ie
,
K
u
}
≡
{
1.0496
,
0.9300
,
0.1757
,
2.2949
}, with J
=
7.2953, and for the fuzzy fractional controller: {
K
e
,
K
ce
,
K
ie
,
K
u
}
≡
{
0.6112
,
1.9142
,
0.1731
,
1.7278
}, with J
=
6.8251. These values lead us again to remain the previous conclusions drawn for G
1
(
s
) and G
2
(
s
), namely, that the fuzzy fractional-order controller produced better results than the integer ones, since the transient response (viz., the overshoot and settling time) and the error J are smaller. Figure 16 shows the ITAE error as function of β. The graph reveals that for this process we have a lower error for β
=
0.9. The fractional controller presents smaller values of J for 0.8
<
β
<
1. Figure 17 illustrates the variation of FLC parameters (K
e, K
ce, K
ie, K
u) as function of the order’s derivative β, while Figure 18 shows the variation of t
s, t
r, t
p, and o
v(%) versus β for the closed-loop step response. The FLC parameters reveal a slight variation with β, as well the transient response parameters.Figure 16
ErrorJversus β for G
3
(
s
).Figure 17
ThePD
β
+
I parameters (K
e, K
ce, K
ie, K
u) versus β for G
3
(
s
).Figure 18
Parameterst
s, t
r, t
p, o
v(%) versus β of the step responses of the closed-loop system with G
3
(
s
) and with a fuzzy PD
β
+
I (0
<
β
≤
1) controller.In conclusion, with the fuzzy fractionalPD
β
+
I controller we get the best controller tuning, superior to the performance revealed by the integer-order scheme. Moreover, we prove the effectiveness of this control structure when used in systems with time delay. In fact, systems with time delay are more difficult to be controlled with the classical methodologies; however, the proposed algorithm reveals to be very effective in the control of this type of systems.


### 4.1. Fuzzy Fractional PID Structures

In this subsection we study the impact of different fuzzy fractional PID structures on system performance. For that, we compare the structure of Figure2  (
S
1
) used in the previous section with the structures of Figure 19  (
S
2
) and Figure 20  (
S
3
) applied to process G
3
(
s
). The GA was executed 10 times, and we get the result that leads to the smaller J. We consider the same fuzzy control rules base of Table 1 in all the structures (S
1, S
2, and S
3).Figure 19
Fuzzy fractional PID controller(
S
2
).Figure 20
Fuzzy fractionalPD
+
PI controller (
S
3
).The control action for these new structures is, respectively, for(
S
2
) and (
S
3
):
(17)
U
(
k
)
=
K
u
v
+
K
i
e
I
v
=
K
u
[
f
(
E
,
CE
)
]
+
K
ie
I
[
f
(
E
,
CE
)
]
=
K
u
[
f
(
K
e
e
(
k
)
,
K
ce
D
β
e
(
k
)
)
]
+
K
i
e
I
[
f
(
K
e
e
(
k
)
,
K
ce
D
β
e
(
k
)
)
]
(18)
U
(
k
)
=
K
PD
v
1
+
K
PI
I
v
2
,
wherev
1
=
f
(
E
1
,
C
E
1
)
=
f
(
K
e
1
e
(
k
)
,
K
ce
1
D
β
e
(
k
)
),
v
2
=
f
(
E
2
,
C
E
2
)
=
f
(
K
e
2
e
(
k
)
,
K
ce
2
D
β
e
(
k
)
).Figure21 shows the unit step responses of fuzzy fractional controller for S
1, S
2, and S
3, for the best J. Figure 22 depicts the corresponding error as function of β. We verify that the lower error occurs for β
=
0.9 in S
1, S
2, and S
3, leading to errors J
=
6.8251, 5.8000, and 6.2900, respectively.Figure 21
Step responses of closed-loop system with fuzzy fractional PID structuresS
1, S
2, and S
3 for G
3
(
s
).Figure 22
ErrorJ versusβ with the structures S
1, S
2, and S
3 for G
3
(
s
).Analyzing Figure21 we observe that the step response is less oscillatory for S
1, with smaller values of t
r and t
s than other structures. Figure 22 shows that all three structures present similar results in terms of error J as function of order β. Structures S
1 and S
2 give the most similar results while S
3 differs mostly on the regions around β
=
0.1, 0.5, and 0.6.The different structures analyzed lead to similar results both in the step responses and errorJ. In fact, it is interesting to notice that the best β
=
0.9 for all the structures when minimizing the ITAE index.Based on this analysis we verify that the proposedPD
β
+
I controller (
S
1
) gives comparable and, in some situations, even better results than other two structures (
S
2 and S
3
). Moreover, the adopted structure is simpler and easier to analyze and implement. These facts justify the choice of the structure S
1 used in this study.


## 4.1. Fuzzy Fractional PID Structures

In this subsection we study the impact of different fuzzy fractional PID structures on system performance. For that, we compare the structure of Figure2  (
S
1
) used in the previous section with the structures of Figure 19  (
S
2
) and Figure 20  (
S
3
) applied to process G
3
(
s
). The GA was executed 10 times, and we get the result that leads to the smaller J. We consider the same fuzzy control rules base of Table 1 in all the structures (S
1, S
2, and S
3).Figure 19
Fuzzy fractional PID controller(
S
2
).Figure 20
Fuzzy fractionalPD
+
PI controller (
S
3
).The control action for these new structures is, respectively, for(
S
2
) and (
S
3
):
(17)
U
(
k
)
=
K
u
v
+
K
i
e
I
v
=
K
u
[
f
(
E
,
CE
)
]
+
K
ie
I
[
f
(
E
,
CE
)
]
=
K
u
[
f
(
K
e
e
(
k
)
,
K
ce
D
β
e
(
k
)
)
]
+
K
i
e
I
[
f
(
K
e
e
(
k
)
,
K
ce
D
β
e
(
k
)
)
]
(18)
U
(
k
)
=
K
PD
v
1
+
K
PI
I
v
2
,
wherev
1
=
f
(
E
1
,
C
E
1
)
=
f
(
K
e
1
e
(
k
)
,
K
ce
1
D
β
e
(
k
)
),
v
2
=
f
(
E
2
,
C
E
2
)
=
f
(
K
e
2
e
(
k
)
,
K
ce
2
D
β
e
(
k
)
).Figure21 shows the unit step responses of fuzzy fractional controller for S
1, S
2, and S
3, for the best J. Figure 22 depicts the corresponding error as function of β. We verify that the lower error occurs for β
=
0.9 in S
1, S
2, and S
3, leading to errors J
=
6.8251, 5.8000, and 6.2900, respectively.Figure 21
Step responses of closed-loop system with fuzzy fractional PID structuresS
1, S
2, and S
3 for G
3
(
s
).Figure 22
ErrorJ versusβ with the structures S
1, S
2, and S
3 for G
3
(
s
).Analyzing Figure21 we observe that the step response is less oscillatory for S
1, with smaller values of t
r and t
s than other structures. Figure 22 shows that all three structures present similar results in terms of error J as function of order β. Structures S
1 and S
2 give the most similar results while S
3 differs mostly on the regions around β
=
0.1, 0.5, and 0.6.The different structures analyzed lead to similar results both in the step responses and errorJ. In fact, it is interesting to notice that the best β
=
0.9 for all the structures when minimizing the ITAE index.Based on this analysis we verify that the proposedPD
β
+
I controller (
S
1
) gives comparable and, in some situations, even better results than other two structures (
S
2 and S
3
). Moreover, the adopted structure is simpler and easier to analyze and implement. These facts justify the choice of the structure S
1 used in this study.


## 5. Conclusion

This paper presented the fundamental aspects of application of the FC theory in fuzzy control systems. In this line of thought, several typical plants were studied. The dynamics of the closed-loop systems were analyzed in the perspective of FC, with the use of a fuzzyPD
β
+
I controller in which the parameters were tuned through a GA algorithm.In general, the control strategies presented give better results than those obtained with conventional integer control structures, showing their effectiveness in the control of nonlinear systems.However, much research is still needed to understand their effective application in systems with larger time delays, and in the presence of noise and different actuator saturation levels.



---

*Source: 676121-2014-05-07.xml*