# Advanced Interval Type-2 Fuzzy Sliding Mode Control for Robot Manipulator

Authors: Ji-Hwan Hwang; Young-Chang Kang; Jong-Wook Park; Dong W. Kim
Journal: Computational Intelligence and Neuroscience
 (2017)

---

## Abstract

In this paper, advanced interval type-2 fuzzy sliding mode control (AIT2FSMC) for robot manipulator is proposed. The proposed AIT2FSMC is a combination of interval type-2 fuzzy system and sliding mode control. For resembling a feedback linearization (FL) control law, interval type-2 fuzzy system is designed. For compensating the approximation error between the FL control law and interval type-2 fuzzy system, sliding mode controller is designed, respectively. The tuning algorithms are derived in the sense of Lyapunov stability theorem. Two-link rigid robot manipulator with nonlinearity is used to test and the simulation results are presented to show the effectiveness of the proposed method that can control unknown system well.


---

## Body

## 1. Introduction

In control engineering, the design of robust controller for a class of uncertain nonlinear multiple-input multiple-output (MIMO) systems remains one of the most challenging tasks. When MIMO systems are nonlinear and uncertain, their control problem becomes more challenging.Conventional control theory is well suited to applications, where the control inputs can be generated based on analytical model [1, 2]. Sliding mode control (SMC), which is based on the theory of variable structure systems (VSS), has been widely applied to robust control of nonlinear systems [3–5]. SMC performs well in trajectory tracking of some nonlinear systems. The SMC employs a discontinuous control law to drive the state trajectory toward a specified sliding surface and maintain its motion along the sliding surface in the state space. Hung et al. [3] have made a comprehensive survey of the VSS theory. The dynamic performance of the SMC system has been confirmed as an effective robust control approach with respect to system uncertainties and unknown disturbance when the system trajectories belong to predetermined sliding surface [4].Although the SMC performs well in the nonlinear systems, it suffers from some difficulties. First, due to the highly coupled nonlinear and uncertain dynamics, it is generally difficult or even impossible for many physical systems to obtain accurate mathematical models. Secondly, to operate effectively in the sliding surface, the SMC requires instantaneous change of the control input without sacrificing the robustness against the model uncertainties and external disturbances. The discontinuity in the control action becomes the cause of chattering, which is undesirable in most applications [6]. In the practical implementation, the chattering may cause an unnecessarily large control signal as the system uncertainties are large and may damage system components such as actuators. Thus, the chattering has to be eliminated or alleviated as much as possible. Finally, it is difficult to directly extend the SMC design into a multiple-input multiple-output (MIMO) system, especially when the coupling among the subsystems is unknown.During the last two decades fuzzy logic system (FLS) has been a dominant topic in intelligent systems research or control community. Because the FLS provide a systematic and efficient framework to incorporate linguistic fuzzy information from human expert, it is particularly suitable for those systems with uncertain or complex dynamics. Owing to universal approximation capability [7] of fuzzy system, many FLS schemes have been developed for handling nonlinear systems, especially in the presence of incomplete knowledge of the system [8, 9].Some researchers applied fuzzy system to sliding mode control to improve the performance of SMC. The fuzzy sliding mode control (FSMC) forms the equivalent control of SMC. By employing the FLS, the set of linearized mathematical model can be integrated into a global model that is equivalent to the nonlinear system [10, 11].As an extension of the well-known ordinary fuzzy set (type-1 fuzzy sets), the concept of type-2 fuzzy sets (T2FS) was first introduced by Zadeh [12]. The sets are fuzzy sets whose membership grades themselves are type-1 fuzzy sets. They are very useful in circumstances where it is difficult to determine an exact membership function for a fuzzy set. They are useful for incorporating uncertainties [13].In this paper, we propose a novel advanced interval type-2 fuzzy sliding mode control (AIT2FSMC) for a class of uncertain nonlinear MIMO systems. To inherit the strength of these two methods, we combine IT2FLS and SMC into one methodology. The AIT2FSMC system is comprised of a fuzzy control design and a hitting control design. For resembling a feedback linearization (FL) control law, IT2 fuzzy system is designed. For compensating the approximation error between the FL control law and IT2 fuzzy system, sliding mode controller is designed, respectively.The tuning algorithms are derived in the sense of Lyapunov stability theorem. The two-link robot manipulator is used to test the proposed method and the simulation results show the AIT2FSMC can control the unknown system well.The organization of this paper is as follows. Problem formulation and notation are presented in Section2. In Section 3, IT2FLS is briefly introduced. Section 4 describes the design process and the stability analysis of AIT2FSMC. In Section 5, the simulation results are presented to show the effectiveness of the proposed control for a two-link robot manipulator. Finally, conclusions are given in Section 6.


## 2. Notation and Problem Formulation

In this section, we present the problem formulation for a class of MIMO nonlinear dynamic systems. Consider the following class of MIMO nonlinear dynamic systems:(1)
y
1
r
1
=
f
1
x
+
∑
j
=
1
p
g
1
j
x
u
j
,
⋮
y
p
r
p
=
f
p
x
+
∑
j
=
1
p
g
p
j
x
u
j
,where x
=
y
1
y
˙
1
⋯
y
1
r
1
-
1
⋯
y
p
y
˙
p
⋯
y
p
(
r
p
-
1
)
T
∈
R
r is the fully measurable state vector and r
1
+
⋯
+
r
p
=
r, u
=
u
1
⋯
u
p
T
∈
R
p is the control input vector, y
=
y
1
⋯
y
p
T
∈
R
p is the output vector, and f
i
x, i
=
1
,
…
,
p are continuous nonlinear functions, and g
i
j
(
x
), i
,
j
=
1
,
…
,
p are continuous nonlinear C
1 functions.Let us denote(2)
y
r
=
y
1
r
1
⋯
y
p
r
p
T
,
F
x
=
f
1
x
⋯
f
p
x
T
,
G
x
=
g
11
x
⋯
g
1
p
x
⋮
⋱
⋮
g
p
1
x
⋯
g
p
p
x
.Then, system (1) can be rewritten in the following compact form:(3)
y
r
=
F
x
+
G
x
u
.The control problem is to design a control lawu
(
t
) which assures that the system tracks a p-dimensional desired vector y
d
=
y
d
1
y
d
2
⋯
y
d
p
T
∈
R
p, which belongs to a class of continuous functions on t
0
,
∞. In this paper, we make the following assumption.Assumption 1.
The matrixG
(
x
) is positive definite; then there exists σ
0
>
0, σ
0
∈
R such that G
(
x
)
≥
σ
0
I
p, with I
p being an identity matrix. In the following σ
0 may be known or not.
Although this assumption restricts the considered class of MIMO nonlinear systems, many physical systems, such as robotic systems [5], fulfill the above property.Assumption 2.
The desired trajectoryy
d
i
t, i
=
1
,
…
,
p, is a known bounded function of time with bounded known derivatives, and y
d
i
t is assumed to be r
i-times differentiable.Let us define the tracking error as(4)
e
1
t
=
y
d
1
t
-
y
1
t
,
⋮
e
p
t
=
y
d
p
t
-
y
p
t
,and the sliding surfaces as(5)
s
1
t
=
d
d
t
+
λ
1
r
1
-
1
e
1
t
,
λ
1
>
0
,
⋮
s
p
t
=
d
d
t
+
λ
p
r
p
-
1
e
p
t
,
λ
p
>
0
.The time derivatives of the sliding surfaces can be written as(6)
s
˙
1
=
v
1
-
f
1
x
-
∑
j
=
1
p
g
1
j
x
u
j
,
⋮
s
˙
p
=
v
p
-
f
p
x
-
∑
j
=
1
p
g
p
j
x
u
j
,where v
1
,
…
,
v
p are given as follows:(7)
v
1
=
y
d
1
r
1
+
β
1
,
r
1
-
1
e
1
r
1
-
1
+
⋯
+
β
1,1
e
˙
1
,
⋮
v
p
=
y
d
p
r
p
+
β
p
,
r
p
-
1
e
p
r
p
-
1
+
⋯
+
β
p
,
1
e
˙
p
,where(8)
β
i
,
j
=
r
i
-
1
!
r
i
-
j
!
j
-
1
!
λ
i
r
i
-
j
,
i
=
1
,
…
,
p
,
j
=
1
,
…
,
r
i
-
1
.Denote(9)
s
t
=
s
1
t
⋯
s
p
t
T
,
v
t
=
v
1
t
⋯
v
p
t
T
.Then, (6) can be written in the compact form(10)
s
˙
=
v
-
F
x
-
G
x
u
.If the nonlinear functionsF
(
x
) and G
(
x
) are known, one can use a sliding mode controller. When the closed loop system is in the sliding mode, it satisfies s
˙
=
0, and then the traditional sliding mode control law is obtained by the following equation: (11)
u
=
u
eq
+
u
h
=
G
-
1
x
-
F
x
+
v
+
K
0
sgn
⁡
s
,where u
eq
=
G
-
1
(
x
)
-
F
(
x
)
+
v is an equivalent control law and u
h
=
G
-
1
(
x
)
K
0
sgn
⁡
s is a hitting control law and K
0
=
diag
⁡
k
01
,
…
,
k
0
p with k
0
i
>
0 for i
=
1
,
…
,
p. Using (10) and (11), we can obtain the following equation:(12)
s
˙
=
-
K
0
sgn
⁡
s
.Multiplyings
T to (12) gives(13)
s
T
s
˙
=
-
s
T
K
0
sgn
⁡
s
.Let us consider the following Lyapunov function candidate:(14)
V
=
1
2
s
T
swhose time derivative is given by(15)
V
˙
=
s
T
s
˙
.With (13), (15) can be reexpressed as(16)
V
˙
=
-
s
T
K
0
sgn
⁡
s
=
-
∑
i
=
1
p
k
0
i
s
i
<
0which implies that s
i
t
→
0 as t
→
∞
. Therefore, e
i
t and all its derivatives up to r
i
-
1 converge to zero [5].According to the above analysis, the control law (11) is easily obtained if the nonlinear functions f
i
x and g
i
j
x are known. However, in this paper, these nonlinear functions are assumed to be unknown, so the above design method cannot be applied directly.


## 3. Interval Type-2 Fuzzy Logic System

The theory and design of interval type-2 fuzzy logic systems (FLS) are presented well in [13–15]. The brief description of the interval type-2 FLS is depicted here. Detailed descriptions can be found in [13–15]. In particular, refer to [13, 15] for more notations and calculations of type-2 fuzzy logic equations.A T2FS in the universal setX is denoted as A
~ which is characterized by a type-2 membership function u
A
~
(
x
) in (17). u
A
~
(
x
) can be referred to as a secondary membership function (MF) or also referred to as secondary set, which is a type-1 set in [
0,1
]. In (17) f
x
(
u
) is a secondary grade, which is the amplitude of a secondary MF; that is, 0
≤
f
x
u
≤
1. The domain of a secondary MF is called the primary membership of x. In (17), J
x is the primary membership of x, where u
∈
J
x
⊆
0,1 for ∀
x
∈
X; u is a fuzzy set in [
0,1
], rather than a crisp point in [
0,1
]. (17)
A
~
=
∫
x
∈
X
u
A
~
x
x
=
∫
x
∈
X
∫
u
∈
J
x
f
x
u
/
u
x
J
x
⊆
0,1
.Whenf
x
u
=
1, ∀
u
∈
J
x
⊆
0,1, then the secondary MFs are interval sets such that u
A
~
x in (17) can be called an interval type-2 MF [13]. Therefore, T2FS A
~ can be rewritten as(18)
A
~
=
∫
x
∈
X
u
A
~
x
x
=
∫
x
∈
X
∫
u
∈
J
x
1
/
u
x
J
x
⊆
0,1
.Also, a Gaussian primary MF with uncertain mean and fixed standard deviation having an interval type-2 secondary MF can be called an interval type-2 Gaussian MF. A 2D interval type-2 Gaussian MF with an uncertain mean in[
m
1
,
m
2
] and a fixed standard deviation σ is shown in Figure 1. It can be expressed as (19)
u
A
~
x
=
exp
⁡
-
1
2
x
-
m
σ
2
,
m
∈
m
1
,
m
2
.Figure 1
Interval type-2 Gaussian fuzzy set with uncertain mean.It is obvious that the T2FS in a region is called a footprint of uncertainty (FOU) and bounded by an upper MF and a lower MF [13], which are denoted as u
-
A
~
x and u
_
A
~
x, respectively. Both of them are type-1 MFs. Hence, (18) can be reexpressed as (20)
A
~
=
∫
x
∈
X
∫
μ
∈
u
_
A
~
x
,
u
-
A
~
x
1
/
u
x
.A T2FLS is very similar to a T1FLS as shown in Figure2 [13], the major structure difference being that the defuzzifier block of a T1FLS is replaced by the output processing block in a T2FLS, which consists of type-reduction followed by defuzzification.Figure 2
The structure of T2FLS.There are five main parts in a T2FLS: fuzzifier, rule base, inference engine, type-reducer, and defuzzifier. A T2FLS is a mappingf
:
R
p
→
R
1. After fuzzification, fuzzy inference, type-reduction, and defuzzification, a crisp output can be obtained.Consider a T2FLS havingp inputs x
1
∈
X
1
,
…
,
x
p
∈
X
p and one output y
∈
Y. The type-2 fuzzy rule base consists of a collection of IF-THEN rules. We assume there are M rules and the rule of a type-2 relation between the input space X
1
×
X
2
×
⋯
×
X
p and the output space Y can be expressed as(21)
R
u
l
e
i
:
I
F
x
1
i
s
F
~
1
i
a
n
d
…
a
n
d
x
p
i
s
F
~
p
i
,
T
H
E
N
y
i
s
G
~
i
,
i
=
1,2
,
…
,
M
,where F
~
j
is are antecedent T2FSs (j
=
1,2
,
…
,
p) and G
~
is are consequent T2FSs.The inference engine combines rules and gives a mapping from input T2FSs to output T2FSs. To achieve this process, we have to compute unions and intersections of type-2 set, as well as compositions of type-2 relations. The output of inference engine block is a type-2 set. By using the extension principle of type-1 defuzzification method, type-reduction takes us from type-2 output sets of the FLS to a type-1 set called the “type-reduced set.” This set may then be defuzzified to obtain a single crisp value.In Figure2, we only consider singleton input fuzzification throughout this paper. Similar to T1FLS, the firing strength F
i in (22) can be obtained by following inference process:(22)
F
i
=
∐
x
∈
X
∏
k
=
1
p
u
F
~
k
i
x
k
,where ∏ is the meet operation and ∐ is the join operation [13].For Gaussian IT2FS as shown in Figure1, the upper MF is a subset that has the maximum membership grade and the lower MF is a subset that has the minimum membership grade. The join operation in (22) leads to joining the result from meet operations, which is using maximum value. The result of join operation can be an interval type-1 set [13] as(23)
F
i
=
f
_
i
f
-
i
T
,where(24)
f
_
i
=
u
_
F
~
1
i
x
1
∗
⋯
∗
u
_
F
~
p
i
x
p
,
f
-
i
=
u
-
F
~
1
i
x
1
∗
⋯
∗
u
-
F
~
p
i
x
p
.There are many kinds of type-reduction, such as centroid, height, modified weight, and center-of-sets [13]. The center-of-sets type-reduction will be used in this paper and can be expressed as (25)
Y
c
o
s
x
=
y
l
,
y
r
=
∫
y
1
∈
y
l
1
,
y
r
1
⋯
∫
y
M
∈
y
l
M
,
y
r
M
∫
f
1
∈
f
_
1
,
f
-
1
⋯
∫
f
M
∈
f
_
M
,
f
-
M
1
∑
i
=
1
M
f
i
y
i
/
∑
i
=
1
M
f
i
,where Y
c
o
s is the interval set determined by two end points y
l and y
r, and firing strengths f
i
∈
F
i
=
f
_
i
,
f
-
i. The interval set y
l
y
r should be computed or set first before the computation of Y
c
o
s
x. For any value y
∈
Y
c
o
s, y can be expressed as (26)
y
=
∑
i
=
1
M
f
i
y
i
∑
i
=
1
M
f
i
,where y is a monotonic increasing function with respect to y
i. Also, y
l in (25) is the minimum associated only with y
l
i, and y
r in (25) is the maximum associated only with y
r
i. Note that y
l and y
r depend only on mixture of f
_
i or f
-
i values. Hence, left-most point y
l and right-most point y
r can be expressed as [13](27)
y
l
=
∑
i
=
1
M
f
l
i
y
l
i
∑
i
=
1
M
f
l
i
,
y
r
=
∑
i
=
1
M
f
r
i
y
r
i
∑
i
=
1
M
f
r
i
.For illustrative purpose, we briefly provide the computation procedure fory
r. Without loss of generality, assume y
r
is are arranged in ascending order; that is, y
r
1
≤
y
r
2
⋯
≤
y
1
M.Step 1.
Computey
r in (27) by initially using f
r
i
=
f
_
i
+
f
-
i
/
2 for i
=
1
,
…
,
M, where f
_
i and f
-
i are precomputed by (24); and let y
r
′
=
y
r.Step 2.
FindR
(
1
≤
R
≤
M
-
1
) such that y
r
R
≤
y
r
′
≤
y
r
R
+
1.Step 3.
Computey
r in (27) with f
r
i
=
f
_
i for i
≤
R and f
r
i
=
f
-
i for i
>
R, and let y
r
′
′
≡
y
r.Step 4.
Ify
r
′
′
≠
y
r
′, then go to Step 5. If y
r
′
′
=
y
r
′, then stop and set y
r
′
′
=
y
r.Step 5.
Sety
r
′ equal to y
r
′
′, and return to Step 2.This algorithm decides the point to separate two sides by the numberR, one side using lower firing strengths f
_
r’s and another side using upper firing strengths f
-
r’s. Hence, y
r in (27) can be reexpressed as(28)
y
r
=
y
r
f
_
1
,
…
,
f
_
R
,
f
-
R
+
1
,
…
,
f
-
M
,
y
r
1
,
…
,
y
r
M
=
∑
i
=
1
R
f
_
i
y
r
i
+
∑
i
=
R
+
1
M
f
-
i
y
r
i
∑
i
=
1
R
f
_
i
+
∑
i
=
R
+
1
M
f
-
i
.The procedure to computey
l is similar to computing y
r. In Step 2, it only needs to find L
1
≤
L
≤
M
-
1, such that y
l
L
≤
y
l
′
≤
y
l
L
+
1. In Step 3, let f
l
i
=
f
-
i for i
≤
L, and f
r
i
=
f
_
i for i
>
L. The y
l in (27) can be also rewritten as(29)
y
l
=
y
l
f
-
1
,
…
,
f
-
L
,
f
_
L
+
1
,
…
,
f
_
M
,
y
l
1
,
…
,
y
l
M
=
∑
i
=
1
L
f
-
i
y
l
i
+
∑
i
=
L
+
1
M
f
_
i
y
l
i
∑
i
=
1
L
f
-
i
+
∑
i
=
L
+
1
M
f
_
i
.The defuzzified crisp output from an IT2FLS is the average of(30)
y
x
=
y
l
+
y
r
2
.


## 4. Interval Type-2 Fuzzy Sliding Mode Control

In this section, we propose an adaptive interval type-2 fuzzy sliding mode controller (AIT2FSMC) for nonlinear unknown MIMO systems. Due to unknown functionsf
i
x and g
i
j
x in our problem, it is impossible to obtain the control law (11). We use the interval type-2 fuzzy system to approximate unknown functions f
i
x and g
i
j
x. First, let the nonlinear functions f
i
x and g
i
j
x be approximated, over a compact set D
X, by interval type-2 fuzzy systems as follows:(31)
f
^
i
x
,
α
~
f
i
=
ξ
f
i
T
x
α
~
f
i
,
i
=
1
,
…
,
p
,
g
^
i
j
x
,
α
~
g
i
j
=
ξ
g
i
j
T
x
α
~
g
i
j
,
i
,
j
=
1
,
…
,
p
,where ξ
f
i
(
x
) and ξ
g
i
j
(
x
) are fuzzy basis vectors fixed by the designer and α
~
f
i and α
~
g
i
j are the corresponding adjustable parameter vectors of each interval type-2 fuzzy system.Let us define(32)
α
~
f
i
∗
=
a
r
g
m
i
n
α
~
f
i
s
u
p
x
∈
D
X
f
i
x
-
f
^
i
x
,
α
~
f
i
,
α
~
g
i
j
∗
=
a
r
g
m
i
n
α
~
g
i
j
s
u
p
x
∈
D
X
g
i
j
x
-
g
^
i
j
x
,
α
~
g
i
jas the optimal parameters of α
~
f
i and α
~
g
i
j, respectively. Notice that optimal parameters α
~
f
i
∗ and α
~
g
i
j
∗ are artificial constant quantities introduced only for analytical purpose, and their values are not needed for the implementation. Define(33)
α
~
-
f
i
=
α
~
f
i
∗
-
α
~
f
i
,
α
~
-
g
i
j
=
α
~
g
i
j
∗
-
α
~
g
i
jas the parameter estimation errors, and(34)
ε
f
i
x
=
f
i
x
-
f
^
i
x
,
α
~
f
i
∗
,
ε
g
i
j
x
=
g
i
j
x
-
g
^
i
j
x
,
α
~
g
i
j
∗as the minimum fuzzy approximation errors, which correspond to approximation errors obtained when optimal parameters are used.In this paper, we assume that the used interval type-2 fuzzy systems do not infringe the universal approximation property on the compact setD
X, which is assumed large enough so that state variables remain within D
X under closed loop control. Therefore, it is reasonable to assume that the minimum approximation errors are bounded for all x
∈
D
X; that is,(35)
ε
f
i
x
≤
ε
-
f
i
,
ε
g
i
j
x
≤
ε
-
g
i
j
,
∀
x
∈
D
X
,where ε
-
f
i and ε
-
g
i
j are given constants.Denote(36)
F
^
x
,
α
~
f
=
f
^
1
x
,
α
~
f
1
⋯
f
^
p
x
,
α
~
f
p
T
,
G
^
x
,
α
~
g
=
g
^
11
x
,
α
~
g
11
⋯
g
^
1
p
x
,
α
~
g
1
p
⋮
⋱
⋮
g
^
p
1
x
,
α
~
g
p
1
⋯
g
^
p
p
x
,
α
~
g
p
p
,
ε
f
x
=
ε
f
1
x
⋯
ε
f
p
x
T
,
ε
g
x
=
ε
g
11
x
⋯
ε
g
1
p
x
⋮
⋱
⋮
ε
g
p
1
x
⋯
ε
g
p
p
x
,
ε
-
f
=
ε
-
f
1
⋯
ε
-
f
p
T
,
ε
-
g
=
ε
-
g
11
⋯
ε
-
g
1
p
⋮
⋱
⋮
ε
-
g
p
1
⋯
ε
-
g
p
p
.From the above analysis, we have(37)
F
x
-
F
^
x
,
α
~
f
=
F
^
x
,
α
~
f
∗
-
F
^
x
,
α
~
f
+
ε
f
x
,
G
x
-
G
^
x
,
α
~
g
=
G
^
x
,
α
~
g
∗
-
G
^
x
,
α
~
g
+
ε
g
x
.Now, let us consider the control law,u
=
u
s, where u
s is a sliding mode control term [4] defined as (38)
u
s
=
G
^
-
1
x
,
α
~
g
-
F
^
x
,
α
~
f
+
v
+
K
0
sgn
⁡
s
.The above control term results from (11) by using the adaptive interval type-2 fuzzy approximation F
^
x
,
α
~
f and G
^
x
,
α
~
g instead of actual functions F
(
x
) and G
(
x
), respectively.The sliding mode control law (38) is not well-defined when the estimated matrix G
^
x
,
α
~
g is singular. The matrix G
^
x
,
α
~
g is generated online via the estimation of the parameters α
~
g. In order to implement this controller, additional precautions have to be made to guarantee that α
~
g remains in a feasible region in which G
^
x
,
α
~
g is regular. Therefore, we modify the sliding mode control term (38) as follows [9]:(39)
u
s
=
G
^
x
,
α
~
g
ε
0
I
p
+
G
^
x
,
α
~
g
G
^
T
x
,
α
~
g
-
1
-
F
^
x
,
α
~
f
+
v
+
K
0
sgn
⁡
s
,where ε
0 is a small positive constant.Within the sliding mode control term (39), we have used the regularized inverse of G
^
-
1
x
,
α
~
g defined as (40)
G
^
x
,
α
~
g
ε
0
I
p
+
G
^
x
,
α
~
g
G
^
T
x
,
α
~
g
-
1
.In fact, the regularized inverse (40) is well-defined even when G
^
x
,
α
~
g is singular, and the sliding mode control term (38) is always well-defined.Even though the control law (39) is always well-defined, it cannot guarantee alone the stability of the closed loop system. It is due, partly, to the approximation of G
^
-
1
x
,
α
~
g by the regularized inverse and, partly, to the unavoidable reconstruction errors of the unknown functions F
(
x
) and G
x. For these reasons, and hoping for the cancellation of these approximations errors, we append to the controller (39) a robustifying control term u
r [8](41)
u
=
u
s
+
u
r
.The controller (41) is the sum of two control terms: a modified sliding mode control term, u
s
(42)
u
s
=
G
^
x
,
α
~
g
ε
0
I
p
+
G
^
x
,
α
~
g
G
^
T
x
,
α
~
g
-
1
-
F
^
x
,
α
~
f
+
v
+
K
0
sgn
⁡
s
,and a robustifying control term, u
r
(43)
u
r
=
s
s
T
ε
-
f
+
ε
-
g
u
s
+
u
0
σ
0
s
2
+
δ
,where u
0 is(44)
u
0
=
ε
0
ε
0
I
p
+
G
^
x
,
α
~
g
G
^
T
x
,
α
~
g
-
1
-
F
^
x
,
α
~
f
+
v
+
K
0
sgn
⁡
sand δ is a design time-varying parameter defined below.In order to meet the control objectives, the adaptive parametersα
~
f
i, α
~
g
i
j, and the design parameter δ are updated by the following adaptive laws:(45)
α
~
˙
f
i
=
-
η
f
i
ξ
f
i
x
s
i
,
(46)
α
~
˙
g
i
j
=
-
η
g
i
j
ξ
g
i
j
x
s
i
u
s
i
,
(47)
δ
˙
=
-
η
0
s
T
ε
-
f
+
ε
-
g
u
s
+
u
0
σ
0
s
2
+
δ
,where η
f
i
,
η
g
i
j
,
η
0
,
δ
0
>
0.Then, we can prove the following theorem.Theorem 3 (consider system (1)).
Suppose that Assumptions1 and 2 are satisfied. Then the control law defined by (41) and (42), with adaptation laws given by (45)–(47), guarantees the following properties:(3.1)
All signals in the closed loop system are bounded.
(3.2)
The tracking errors and its derivatives decrease asymptotically to zero; that is,e
i
j
(
t
)
→
0 as t
→
∞ for i
=
1
,
…
,
p and j
=
0,1
,
…
,
r
i
-
1.Proof.
Using the control law (51), (10) can be rewritten as(48)
s
˙
=
v
-
F
x
-
G
x
-
G
^
x
,
α
~
g
u
c
-
G
^
x
,
α
~
g
u
s
-
G
x
u
r
.
By introducing the control term (42)–(48), we obtain(49)
s
˙
=
-
K
0
sgn
⁡
s
-
F
x
-
F
^
x
,
α
~
f
-
G
x
-
G
^
x
,
α
~
g
u
s
+
u
0
-
G
x
u
r
.
Here, we have used the fact that(50)
G
^
x
,
α
~
g
G
^
T
x
,
α
~
g
ε
0
I
p
+
G
^
x
,
α
~
g
G
^
T
x
,
α
~
g
-
1
=
I
p
-
ε
0
ε
0
I
p
+
G
^
x
,
α
~
g
G
^
T
x
,
α
~
g
-
1
.
From (37), one can write (49) as(51)
s
˙
=
-
K
0
sgn
⁡
s
-
F
^
x
,
α
~
f
∗
-
F
^
x
,
α
~
f
-
G
^
x
,
α
~
g
∗
-
G
^
x
,
α
~
g
u
s
+
u
0
-
G
x
u
r
+
u
0
-
ε
f
x
-
ε
g
x
u
s
.
Multiplyings
T to (51) gives(52)
s
T
s
˙
=
-
s
T
K
0
sgn
⁡
s
-
∑
i
=
1
p
ξ
f
i
T
x
α
~
-
f
i
s
i
-
∑
i
=
1
p
∑
j
=
1
p
ξ
g
i
j
T
x
α
~
-
g
i
j
s
i
u
s
j
-
s
T
G
x
u
r
+
s
T
u
0
-
s
T
ε
f
x
-
s
T
ε
g
x
u
s
.
Let us now consider the following Lyapunov function candidate:(53)
V
=
1
2
s
T
s
+
1
2
∑
i
=
1
p
1
η
f
i
α
~
-
f
i
T
α
~
-
f
i
+
1
2
∑
i
=
1
p
∑
j
=
1
p
1
η
g
i
j
α
~
-
g
i
j
T
α
~
-
g
i
j
+
1
2
η
0
δ
2
,whose time derivative is given by(54)
V
˙
=
s
T
s
˙
-
∑
i
=
1
p
1
η
f
i
α
~
-
f
i
T
α
~
-
˙
f
i
-
∑
i
=
1
p
∑
j
=
1
p
1
η
g
i
j
α
~
-
g
i
j
T
α
~
-
˙
g
i
j
+
1
η
0
δ
δ
˙
.
With (52), (53) can be expressed as (55)
V
˙
=
s
T
K
0
sgn
⁡
s
+
V
˙
1
+
V
˙
2
,where(56)
V
˙
1
=
-
∑
i
=
1
p
α
~
-
f
i
T
ξ
f
i
x
s
i
+
1
η
f
i
α
~
-
˙
f
i
-
∑
i
=
1
p
∑
j
=
1
p
α
~
-
g
i
j
T
ξ
g
i
j
x
s
i
u
s
j
+
1
η
g
i
j
α
~
-
˙
g
i
j
,
(57)
V
˙
2
=
-
s
T
G
x
u
r
+
s
T
u
0
-
s
T
ε
f
x
-
s
T
ε
g
x
u
s
+
1
η
0
δ
δ
˙
.
Substituting the parameter adaptive laws (45) and (46) into (56) gives (58)
V
˙
1
=
0
.
Using (43), we can write(59)
s
T
G
x
u
r
≥
s
T
ε
-
f
+
ε
-
g
u
s
+
u
0
-
s
T
ε
-
f
+
ε
-
g
u
s
+
u
0
σ
0
s
2
+
δ
.
Here, we have used the inequality(60)
s
T
G
x
s
≥
σ
0
s
2
,which is true because G
(
x
) is assumed positive definite and satisfies G
(
x
)
≥
σ
0
I
p.
Equation (57) can be bounded as follows:(61)
V
˙
2
≤
-
s
T
G
x
u
r
+
s
T
ε
-
f
+
ε
-
g
u
s
+
u
0
+
1
η
0
δ
δ
˙
.
With (60), (61) becomes(62)
V
˙
2
≤
-
δ
s
T
ε
-
f
+
ε
-
g
u
s
+
u
0
σ
0
s
2
+
δ
+
1
η
0
δ
δ
˙
.
Using (47) in (62) yields(63)
V
˙
2
≤
0
.
From (58) and (63), it follows that(64)
V
˙
≤
-
s
T
K
0
sgn
⁡
s
=
-
∑
i
=
1
p
k
o
i
s
i
.
By Barbalat’s lemma [5], it can conclude that s
→
0 as t
→
∞. In spite of the demonstrated properties of the controller, the hitting control law leads to the well-known chattering phenomenon. In order to overcome the undesirable chattering effects, the sign function is replaced with the saturation function [5].


## 5. Simulation Results

In this section, we test the AIT2FSMC design on the tracking control of a two-link robot. Consider a two-link rigid robot manipulator moving a horizontal plant in Figure3. The first link is mounted on a rigid base by means of frictionless hinges and the second is mounted at the end of first link by means of a frictionless ball bearing. The dynamic equations of this MIMO system are given by [5] (65)
q
¨
1
q
¨
2
=
M
11
M
12
M
21
M
22
-
1
u
1
u
2
-
-
h
q
˙
2
-
h
q
˙
1
+
q
˙
2
h
q
˙
1
0
q
˙
1
q
˙
2
,where(66)
M
11
=
a
1
+
2
a
3
cos
⁡
q
2
+
2
a
4
sin
⁡
q
2
,
M
22
=
a
2
,
M
12
=
M
21
=
a
2
+
a
3
cos
⁡
q
2
+
a
4
sin
⁡
q
2
h
=
a
3
s
i
n
q
2
-
a
4
c
o
s
q
2
,with(67)
a
1
=
I
1
+
m
1
l
c
1
2
+
I
e
+
m
e
l
c
e
2
+
m
e
l
1
2
,
a
2
=
I
e
+
m
e
l
c
e
2
,
a
3
=
m
e
l
1
l
c
e
cos
⁡
δ
e
,
a
4
=
m
e
l
1
l
c
e
sin
⁡
δ
e
.Figure 3
Planar model of the two-link manipulator.In the simulation, the following parameter values are used:(68)
m
1
=
1
,
m
e
=
2
,
l
1
=
1
,
l
c
l
=
0.5
,
l
c
e
=
0.6
,
I
1
=
0.12
,
I
e
=
0.25
,
δ
e
=
3
0
°
.Lety
=
y
1
y
2
T
=
q
1
q
2
T, u
=
u
1
u
2
T, x
=
q
1
q
˙
1
q
2
q
˙
2
T, and(69)
F
x
=
f
1
x
f
2
x
=
-
M
-
h
q
˙
2
-
h
q
˙
1
+
q
˙
2
h
q
˙
1
0
q
˙
1
q
˙
2
,
G
x
=
g
11
x
g
12
x
g
21
x
g
22
x
=
M
-
1
=
M
11
M
12
M
21
M
22
,and then, the robot system (65) can be described as follows:(70)
y
¨
=
F
x
+
G
x
u
,which is the input-output form given by (3). Since the matrix M is positive definite [5], then it is always regular and G
(
x
)
=
M
-
1 is positive definite.The control objective is to force the system outputsq
1 and q
2 to track the sinusoidal desired trajectories y
d
1
=
sin
⁡
t and y
d
2
=
sin
⁡
t. In order to analyze the performance of the AIT2FSMC, we compared the AIT2FSMC with the A-Fuzzy Sliding Mode Controller (AFSMC) which used the type-1 FLS to approximate the nonlinear F
x and G
x. The external disturbances cos
⁡
t
sin
⁡
t
Tare added to system (65). Since the components of F
x and G
x are assumed unknown, two fuzzy systems in the form of (30) are used to approximate the elements of F
x, and four are used to approximate the elements of G
x. In the AIT2FSMC and the AFSMC, the sliding surface is selected as with λ
1
,
λ
2
=
5 and the design parameters used in this simulation are chosen as follows: K
0
=
0.3
I
2, ε
0
=
0.1, η
f
i
=
0.5, η
g
i
j
=
0.5 for i
,
j
=
2 and the initial conditions of robot are selected as x
(
0
)
=
0.5
0
0.25
0. The fuzzy systems used to describe F
x have q
1
(
t
), q
˙
1
(
t
), q
2
(
t
), and q
˙
2
(
t
) as inputs. The input membership functions and parameters for the AIT2FSMC and AFSMC are shown in Table 1.Table 1
The parameters of input membership functions for the AIT2FSMC and the AFSMC.
Negative
Zero
Positive
m
N
1
m
N
2
σ
N
m
Z
1
m
Z
2
σ
Z
m
P
1
m
P
2
σ
P
AIT2FSMC
-
1.15
-
0.35
0.6
-
0.01
0.01
0.6
1.15
1.35
0.6
AFSMC
-
1.25
-
1.25
0.6
0
0
0.6
1.25
1.25
0.6As shown in Figures4 and 5, the AIT2FSMC shows the better performance than the AFSMC. In the AFSMC, a type-1 FLS, which is not able to handle rule uncertainties, is used for control of unknown nonlinear MIMO system. Therefore, the system performance is deteriorated by the disturbance. Meanwhile, the proposed AIT2FSMC utilizes the interval type-2 FLS. The simulation results show that the interval type-2 FLS is able to handle rule uncertainties, and thus the system performance is compensated by the interval type-2 FLS [13].Figure 4
Tracking results of link 1.
(a)
Tracking curve ofq
1
(b)
Tracking curve ofq
˙
1
(c)
Tracking error ofq
1
(d)
Tracking error ofq
˙
1Figure 5
Tracking results of link 2.
(a)
Tracking curve ofq
2
(b)
Tracking curve ofq
˙
2
(c)
Tracking error ofq
2
(d)
Tracking error ofq
˙
2


## 6. Discussion

In this paper, we propose a novel advanced interval type-2 fuzzy sliding mode control (AIT2FSMC) for a class of uncertain nonlinear MIMO systems with external disturbances. The parameters of the proposed AIT2FSMC system, as well as the approximation error bound, are tuned online. The control laws are obtained in the Lyapunov sense to ensure the stability of the control system.Unlike the conventional SMCs, the design of the proposed AIT2FSMC is independent of the mathematical model of the system and can be applied to both unknown and uncertain nonlinear MIMO systems. Furthermore, the uncertainty bound is not needed to be available beforehand. Simulation results performed on a two-link robot manipulator demonstrate the feasibility of the proposed control system.



---

*Source: 9640849-2017-02-08.xml*