Fuzzy Logic Control System Stability Analysis Based On Lyapunov's Direct Method
Fuzzy Logic Control System Stability Analysis Based On Lyapunov's Direct Method
Fuzzy Logic Control System Stability Analysis Based On Lyapunov's Direct Method
, (2)
where r is the total number of rules, X
i,1
, X
i,:
, .., X
i,n
are fuzzy sets that describe the linguistics terms (LTs)
of the input variables x
k
, k = 1, n, u = u
i
(x) is the control signal of rule i, similar to the case of parallel
distributed compensation, and the function AND is a t-norm. u
i
can be a single value or a function of the
state vector, x.
The structure presented in Figure 1 cam be viewed as a nonlinear state-feedback control system.
However other input variables (to the FLC) can be considered as well instead of the state variables x
k
, k =
1, n. One simple design of the fuzzy logic control system can be done in terms of parallel distributed
compensation.
Each fuzzy rule generates the ring strength dened in (3):
i
(x)=AND(
i,1
(x
1
),
i,:
(x
:
). . .
i,n
(x
n
)) [o, 1] , x X, i = 1, r. (3)
Fuzzy Logic Control System Stability Analysis Based on Lyapunovs Direct Method 417
It is assumed that for any x X there exists among all rules at least one
i
(o, 1] , i = 1, r. The control
signal u is a function of
i
and u
i
. Applying the weighted sum defuzzication method the output of the
FLC is given by
u =
r
i=1
i
u
i
r
i=1
i
. (4)
Denition 1. For any input x
o
X if the ring strength
i
(x
o
) corresponding to the fuzzy rule i is zero,
that fuzzy rule i, i = 1, r, is called an inactive fuzzy rule for the input x
o
; otherwise, it is called an active
fuzzy rule.
It should be noted that with x = x
o
an inactive fuzzy rule will not affect the controller output u(x
o
).
Hence (4) can be rewritten as follows aiming the consideration of all active fuzzy rules only:
u(x
o
) =
r
i=1,
i
=o
i
(x
o
)u
i
(x
o
)
r
i=1,
i
=o
i
(x
o
)
. (5)
Denition 2. An active region of the fuzzy rule i is dened as a set
X
A
i
= {x X |
i
(x) = o} , i = 1, r. (6)
3 Stability Analysis
The stability analysis presented in this paper is based on LaSalles invariance principle cited and
analyzed in [11]. This Section is concentrated on the formulation and proof of Theorem 1 that ensures
sufcient conditions for the stability of nonlinear processes controlled by T-S FLCs.
The Lyapunov function candidate V : R
n
R, V (x) = x
T
Px is considered. It is positive and un-
bounded, where P R
nn
is a positive denite matrix. Considering the state trajectories fullling (1) in
order to obtain the closed-loop system dynamics, it results that V has continuous partial derivatives and
the derivatives of V with respect to time expressed in terms of (7):
V (x) = x
T
Px +x
T
P x = ( f (x) +b(x)u(x))
T
Px+
+x
T
P( f (x)+b(x)u(x)) = F (x) +B(x)u(x),
(7)
where:
F (x) = f (x)
T
Px +x
T
Pf (x), B(x) = b(x)
T
Px +x
T
Pb(x). (8)
The following sets are dened to be used in the stability analysis:
B
o
= {x X |B(x) = o} , B
+
= {x X |B(x) o} , B
, i = 1, r,
418 R.-E. Precup, M.-L. Tomescu, St. Preitl
3. the set
x X
V (x) = o
i=1,
i
=o
i
(x
o
)u
i
(x
o
)
r
i=1,
i
=o
i
(x
o
)
F(x
o
)
B(x
o
)
r
i=1,
i
=o
i
(x
o
)
r
i=1,
i
=o
i
(x
o
)
=
F (x
o
)
B(x
o
)
V (x
o
)=F (x
o
) +B(x
o
)u(x
o
) F (x
o
) +B(x
o
)
F (x
o
)
B(x
o
)
= o. (10)
Therefore,
u
i
(x
o
)
F (x
o
)
B(x
o
)
V (x
o
) o. (11)
Case 2: B(x
o
) is strictly negative. Once more, from the condition 2 of Theorem 1 it results that
u
i
(x
o
)
F (x
o
)
B(x
o
)
u(x
o
) =
r
i=1,
i
=o
i
(x
o
)u
i
(x
o
)
r
i=1,
i
=o
i
(x
o
)
F(x
o
)
B(x
o
)
r
i=1,
i
=o
i
(x
o
)
r
i=1,
i
=o
i
(x
o
)
=
F (x
o
)
B(x
o
)
V (x
o
) = F (x
o
) +B(x
o
)u(x
o
)F (x
o
) +B(x
o
)
F (x
o
)
B(x
o
)
= o. (12)
Therefore,
u
i
(x
o
)
F (x
o
)
B(x
o
)
V (x
o
) o. (13)
Case 3: x
o
B
o
. In this case using the condition 1 in Theorem 1 the result will be F (x
o
) o. Hence,
V (x
o
) = F (x
o
) +B(x
o
)u(x
o
) = F (x
o
) o. (14)
From the above three cases it is obtained that
V (x) o, x X. (15)
In conclusion, the derivative with respect to time of the Lyapunov function candidate,
V, is negative
semi-denite.
The condition 3 ensures the fullment of LaSalles invariance principle. This justies the fact that
the equilibrium point at the origin is globally asymptotically stable.
The proof is now complete.
The stability theorem presented here ensures sufcient conditions for the stability of the fuzzy logic
control system described in Section 2. So it has been proved that if the Lyapunov function candidate is
negative semi-denite in the active region of each fuzzy rule then, the closed-loop system will be globally
asymptotically stable in the sense of Lyapunov.
Fuzzy Logic Control System Stability Analysis Based on Lyapunovs Direct Method 419
The conditions 1 and 2 in Theorem 1 guarantee that the function
V is negative semi-denite in the
active region of each fuzzy rule. The condition 3 proves that the set {o} is the largest invariance set in
x X
V (x) = o
. By LaSalles invariance principle it has been guaranteed that the fuzzy logic control
system, comprising the nonlinear process described by (1) and the T-S FLC, is globally asymptotically
stable in the sense of Lyapunov at the origin.
The stability analysis algorithm ensuring the stability of the class of fuzzy logic control systems
considered in Section 2 is based on Theorem 1. It consists of the following steps:
1. Set the Lyapunov function candidate V (i.e. set P).
2. Check that the set
x X
V (x) = o
, B
+
.
4. If F (x) o, x B
o
then go to step 5. Else go to step 1.
5. For each fuzzy control rule i determine u
i
such that u
i
(x)
F(x)
B(x)
for x X
A
i
B
+
and u
i
F(x)
B(x)
for x X
A
i
B
, i = 1, r.
The application of this algorithm will be illustrated in the next Section.
4 Illustrative example
This Section is dedicated to the validation of the theoretical results derived in Section 3 by the de-
sign of a stable fuzzy logic control system with T-S FLC controlling a nonlinear process, the inverted
pendulum on a cart system. This simple mechanical system is representative to model a class of attitude
control problems whose goal is to maintain permanently the desired vertically oriented position. Since
the inverted pendulum is a nonlinear system, the basic balance equations for the system are derived rstly
and put into the standard state-space form. Given an inverted pendulum mounted on a cart as shown in
Figure 2, the rst principle nonlinear equations are applied in the sequel. Assuming that the rod is mass-
less and that the cart mass and the point mass at the upper end of the inverted pendulum are denoted as M
and m, respectively, there is an externally x-directed force on the cart, F(t), and the gravity force acts on
the point mass at all times. The coordinate system is dened according to Figure 2, where x(t) represents
the cart position and (t) is the tilt angle referenced to the vertical upward direction.
Figure 2: Variables related to the inverted pendulum on a cart system.
The differential equation that describes the behavior of the simplied system, playing the role of
controlled process, is usually written as
(m+M) l
:
(m+M) l g sin() = u, (16)
420 R.-E. Precup, M.-L. Tomescu, St. Preitl
where:
M - the mass of the cart,
m - the mass of the pendulum,
l - the length of pendulum (distance to the center of mass),
x - the cart position coordinate,
- the pendulum angle with respect to the vertical position,
u - the control signal, equal to the externally x-directed force, u = F.
The state vector consists of the angle, , and the angular velocity of the pendulum,
. Therefore, the
two state variables are dened as z
1
and z
:
, where z
1
[8o, 8o], z
:
[o, o], z
1
(t) = (t) and
z
:
(t) =
(t). In order to write equation (19) in terms of state variables, they are substituted resulting in
z = f (z) +b(z)u, (17)
where: z =
z
1
z
:
z
:
g
l
sin(z
1
)
, b(z) =
1
(m+M)l
:
.
The goal of fuzzy logic control system design, to be presented as follows, is to ensure the upright
stabilization of the pendulum aiming the setpoint value of z, z = o. The design starts with setting the
fuzzication module of the T-S FLC. Figures 3 and 4 illustrates the membership functions corresponding
to the LTs of the linguistic variables z
1
and z
:
. The three LTs representing Positive, Zero and Negative
values are noted by P, Z and N, respectively.
Figure 3: Membership functions of z
1
.
Figure 4: Membership functions of z
:
.
Fuzzy Logic Control System Stability Analysis Based on Lyapunovs Direct Method 421
The inference engine of the FLC employs the MIN and MAX operators and it is assisted by the
complete rule base illustrated in Table 1. The weighted sum defuzzication method is used in the T-S
FLC structure. Summarizing, the only parameters to be calculated are the consequents u
i
in the 9 fuzzy
control rules.
Table 1
Fuzzy Control Rule Base
Rule Antecedent Consequent
z
1
z
:
u
1 P P u
1
2 N N u
:
3 P N u
4 N P u
5 P Z u
6 N Z u
6
7 Z P u
8 Z N u
8
9 Z Z u
The algorithm presented in Section 3 will be applied as follows in order to nd the values of u
i
for which
the system (19) can be stabilized by the above described T-S FLC.
Step 1: The Lyapunov function candidate V (z) = z
T
Pz = z
:
1
+z
:
:
is considered, where P =
1 o
o 1
.
Therefore V is positive. If z then V (z) . V (z) o, z = o.
Step 2: The derivative is
V(z)= :z
:
z
1
+
g
l
sin(z
1
)
1
(m+M)l
:
u
and
V (o) = o. Assume that there is a
trajectory with z
:
(t) =o and z
1
(t) =o. Then
d
dt
z
:
(t) =
g
l
sin(z
1
(t))
1
(m+M)l
:
u(t) =o, which means that
z
:
(t) can not stay constant. Hence, z(t) = o is the only possible state trajectory for which
V (z) = o. So
the set
z X
V (z) = o
contains no trajectory of the system except the trivial trajectory z(t) = o for
t o.
Step 3: The expressions of F and B are:
F (z) = :z
:
z
1
+
g
l
sin(z
1
)
, B(z) =
:z
:
(m+M)l
:
. (18)
The following elements necessary in Theorem 1 obtain the particular values expressed in (23) to (26):
B
o
= {(z
1
, o) X |z
1
[1, 1] } , B
+
= {(z
1
, z
:
) X |z
:
< o} , B
= {(z
1
, z
:
) X |z
:
o} , (19)
F (z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)). (20)
Step 4: If z B
o
then z
:
= o and F (z) = o.
Step 5: Each rule will be analyzed further on. This is not a complex task since only 9 rules are involved.
For rule 1: z
1
IS P, z
:
IS P. So X
A
1
= (o, 8o] (o, o], X
A
1
B
+
= / 0 and X
A
1
B
= / 0. Thus u
:
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)). It is taken u
:
(z) = l (m+M)(z
1
l g), and this
function fullls the condition 2 in Theorem 1.
422 R.-E. Precup, M.-L. Tomescu, St. Preitl
For rule 3: z
1
IS P, z
:
IS N. So X
A
= / 0 and X
A
B
+
= (o, 8o] (o, o].
Thus u
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)). It is taken u
(z) = z
1
, and this function fullls again
the given condition.
For rule 4: z
1
IS N, z
:
IS P. So X
A
B
+
= / 0 and X
A
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)). It is taken u
(z) = z
1
, and this function fullls the
given condition.
For rule 5: z
1
IS P, z
:
IS Z. So X
A
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)) and
b) for z X
A
B
+
= (o, 8o] (, o] u
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)).
In order to satisfy both conditions it is chosen u
(z) = l (m+M)(z
1
l +gsin(z
1
)).
For rule 6: z
1
IS N, z
:
IS Z. So X
A
6
= [8o, o) (, ). The result will be:
a) for z X
A
6
B
= [8o, o) (o, ) u
6
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)) and
b) for z X
A
6
B
+
= [8o, o) (, o) u
6
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)).
In order to satisfy both conditions it is chosen u
6
(z) = l (m+M)(z
1
l +gsin(z
1
)).
For rule 7: z
1
IS Z, z
:
IS P. So X
A
B
+
= / 0 and X
A
= (1o, 1o)
(o, o] . Thus u
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)). It is set u
(z) = l (m+M)(z
1
l +g) to fulll
the condition 2 in Theorem 1.
For rule 8: z
1
IS Z, z
:
IS N. So X
A
8
= (1o, 1o) [o, o), X
A
8
B
+
= (1o, 1o) [o, o) and
X
A
8
B
= / 0. Thus u
8
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)). It is set u
8
(z) = l (m+M)(z
1
l g) to
fulll the condition 2 in Theorem 1.
For rule 9: z
1
IS Z, z
:
IS Z. So X
A
B
+
= (1o, 1o) (, o) u
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)) and
b) for z X
A
(z)
F(z)
B(z)
= l (m+M)(z
1
l +gsin(z
1
)).
In order to satisfy both conditions it is set u
(z) = l (m+M)(z
1
l +gsin(z
1
)).
Concluding, from Theorem 1 it results that the closed-loop system composed by the nonlinear pro-
cess modeled in (19) and the T-S FLC designed here is globally asymptotically stable in the sense of
Lyapunov at the origin. Considering the values of process parameters m = o., M = o., l = 1, g = .8,
the responses of z
1
and z
:
versus time in the closed-loop system are presented in Figures 5 to 8 for
different initial conditions.
5 Summary and Conclusions
A new approach to the global asymptotic stability analysis of fuzzy logic control systems employing
T-S FLCs dedicated to a class of nonlinear processes has been introduced. The example proves how the
stability analysis algorithm suggested here can be applied to the design of a stable fuzzy logic control
system for a nonlinear process. The new stability approach can be applied also in situations when the
system has an equilibrium point different to the origin and / or the setpoint is nonzero by an appropriately
dened state transformation [15].
Fuzzy Logic Control System Stability Analysis Based on Lyapunovs Direct Method 423
>V@
Figure 5: State variables versus time for fuzzy logic control system with T-S FLC in the condition
z
1
(o) = o and z
:
(o) = o.
>V@
Figure 6: State variables versus time for fuzzy logic control system with T-S FLC in the condition
z
1
(o) = o and z
:
(o) = o.
>V@
Figure 7: State variables versus time for fuzzy logic control system with T-S FLC in the condition
z
1
(o) = o and z
:
(o) = 6o.
>V@
Figure 8: State variables versus time for fuzzy logic control system with T-S FLC in the condition
z
1
(o) = :o and z
:
(o) = 1o.
The stability analysis algorithm suggested in this paper canbe applied also when the rule base (2)
of the T-S FLC is not complete. However interpolation techniques [10, 19] are needed in the imple-
424 R.-E. Precup, M.-L. Tomescu, St. Preitl
mentation of the T-S FLC. They require the re-assessment of the stability conditoins derived prior to the
implementation.
Further research will be concentrated on newapplications of the proposed algorithmto several classes
of processes [4, 5, 7, 8, 9, 17, 20]. The complex applications require the computer-aided design of the
Takagi-Sugeno fuzzy logic controllers employing the stability analysis algorithm proposed in this paper
to strive for increased generality.
Acknowledgement
The support from the CNCSIS and CNMP of Romania is acknowledged.
Bibliography
[1] M. M. B ala s, V. E. B ala s, World Knowledge for Control Applications by Fuzzy-Interpolative Sys-
tems, International Journal of Computers, Communications & Control, Vol. III, Suppl. Issue: Pro-
ceedings of ICCCC 2008, pp. 28-32, 2008.
[2] A. Bazoula, M.S. Djouadi, H. Maaref, Formation Control of Multi-Robots via Fuzzy Logic Tech-
nique, International Journal of Computers, Communications & Control, Vol. III, Suppl. Issue: Pro-
ceedings of ICCCC 2008, pp. 179-184. 2008.
[3] D. Bellomo, D. Naso, and R. Babuka, Adaptive fuzzy control of a non-linear servo-drive: theory
and experimental results, Engineering Applications of Articial Intelligence , Vol. 21, Issue 6, pp.
846-857, 2008.
[4] S. Blai c and I. krjanc, Design and stability analysis of fuzzy model-based predictive control - a
case study, Journal of Intelligent and Robotic Systems, Vol. 49, Issue 3, pp. 279-292, 2007.
[5] S. Ertugrul, Predictive modeling of human operators using parametric and neuro-fuzzy models by
means of computer-based identication experiment, Engineering Applications of Articial Intelli-
gence , Vol. 21, Issue 2, pp. 259-268, 2008.
[6] G. Feng, A survey on analysis and design of model-based fuzzy control systems, IEEE Transactions
on Fuzzy Systems, Vol. 14, Issue 5, pp. 676-697, 2006.
[7] D. Hladek, J. Vascak, and P. Sincak, Multi-robot control system for pursuit-evasion problem, Journal
of Electrical Engineering, Vol. 60, Issue 3, pp. 143-148, 2009.
[8] L. Horvth and I. J. Rudas, Towards content oriented integration of product and robot system models,
International Symposium on Logistics and Industrial Informatics, Wildau, pp. 7-12, 2007.
[9] G.-H. Hwang, D.-W. Kim, J.-H. Lee, and Y.-J. An, Design of fuzzy power system stabilizer using
adaptive evolutionary algorithm, Engineering Applications of Articial Intelligence, Vol. 21, Issue
1, pp. 86-96, 2008.
[10] Z. C. Johanyk and S. Kovcs: Sparse fuzzy system generation by rule base extension, 11th IEEE
International Conference of Intelligent Engineering Systems, Budapest, pp. 99-104, 2007.
[11] H. K. Khalil, Nonlinear Systems, Third Edition, Prentice Hall, Englewood Cliffs, NJ, 2002.
[12] K. Michels, F. Klawonn, R. Kruse, and A. Nrnberger, Fuzzy Control: Fundamentals, Stability and
Design of Fuzzy Controllers, Springer Verlag, Berlin, Heidelberg, New York, 2006.
Fuzzy Logic Control System Stability Analysis Based on Lyapunovs Direct Method 425
[13] B. M. Mohan and A. Sinha, Analytical structure and stability analysis of a fuzzy PID controller,
Applied Soft Computing, Vol. 8, Issue 1, pp. 749-758, 2008.
[14] A. G. Perry, G. Feng, Y.-F. Liu, and P. C. Sen, A design method for PI-like fuzzy logic controllers
for DC-DC converter, IEEE Transactions on Industrial Electronics, Vol. 54, Issue 5, pp. 2688-2696,
2007.
[15] R.-E. Precup, S. Preitl, and G. Faur, PI predictive fuzzy controllers for electrical drive speed control:
methods and software for stable development, Computers in Industry, Vol. 52, Issue 3, pp. 253-270,
2003.
[16] R.-E. Precup, S. Preitl, I. J. Rudas, M. L. Tomescu, and J. K. Tar, Design and experiments for a
class of fuzzy controlled servo systems, IEEE/ASME Transactions on Mechatronics, Vol. 13, Issue
1, pp. 22-35, 2008.
[17] E. Tian and C. Peng, Delay-dependent stability analysis and synthesis of uncertain T-S fuzzy sys-
tems with time-varying delay, Fuzzy Sets and Systems, Vol. 157, Issue 4, pp. 544-559, 2006.
[18] T. Vesselenyi, S. Dzitac, I. Dzitac and M.-J. Manolescu, Fuzzy and Neural Controllers for a Pneu-
matic Actuator, International Journal of Computers, Communications & Control, Vol. II, No. 4, pp.
375-387, 2007.
[19] K. W. Wong, D. Tikk, T. D. Gedeon, and L. T. Kczy, Fuzzy rule interpolation for multidimensional
input spaces with applications: a case study, IEEE Transactions on Fuzzy Systems, Vol. 13, Issue 6,
pp. 809-819, 2005.
[20] J. Yoneyama, Robust stability and stabilization for uncertain Takagi-Sugeno fuzzy time-delay sys-
tems, Fuzzy Sets and Systems, Vol. 158, Issue 2, pp. 115-134, 2007.
Radu-Emil Precup was born in Lugoj, Romania, in 1963. He received the Dipl.Ing. degree
(Hons.) in Automation and Computers from the "Traian Vuia" Polytechnic Institute of Timi soara,
Romania, in 1987, the diploma in Mathematics from the West University of Timi soara in 1993,
and the PhD degree in Automatic Systems from the "Politehnica" University of Timisoara (PUT),
Romania, in 1996. From 1987 to 1991, he was with Infoservice S.A., Timi soara. He is currently
with the PUT, where he became a Professor in the Department of Automation and Applied
Informatics, in 2000, and he is currently a Doctoral Supervisor of Automation and Systems
Engineering. He is also an Honorary Professor with the Budapest Tech Polytechnical Institution,
Budapest, Hungary. He is the author or coauthor of more than 120 papers published in various
scientic journals, refereed conference proceedings, and contributions to books.
Marius L. Tomescu was born in Arad, Romania, in 1967. He received his MSc degree in
Computer Science from the West University of Timi soara, Romania, in 1993, and the PhD degree
in Automation from the "Politehnica" University of Timi soara, Romania, in 2008. From 1993 to
2008, he was with Grup Service S.R.L., Arad. He is with the "Aurel Vlaicu" University of Arad,
where he has been a Lecturer in the Department of Computer Science since 2002. His research
interests include fuzzy control applications and stability analysis of fuzzy control systems. He
is the author or coauthor of more than 30 papers published in various scientic journals and
refereed conference proceedings. Dr. Tomescu is a member of the Romanian Society of Control
Engineering and Technical Informatics.
426 R.-E. Precup, M.-L. Tomescu, St. Preitl
Stefan Preitl was born in Timi soara, Romania, in 1943. He received the Dipl.Ing. degree in
Electrical Engineering, and the PhD degree in Measurement Techniques from the "Traian Vuia"
Polytechnic Institute of Timi soara in 1966 and 1983, respectively. From 1967 to 1972, he was
with Electromotor S.A., Timi soara. He is currently with the "Politehnica" University of Tim-
i soara, Timi soara, where he became a Professor in the Department of Automation and Applied
Informatics, in 1992, and he is currently a Doctoral Supervisor of Automation and Systems En-
gineering. He is also an Honorary Professor with the Budapest Tech Polytechnical Institution,
Budapest, Hungary. He is the author or coauthor of more than 170 papers published in various
scientic journals, refereed conference proceedings, and books in the eld of automatic control.