Open access peer-reviewed chapter

Fuzzy Control Systems: LMI-Based Design

By Morteza Seidi, Marzieh Hajiaghamemar and Bruce Segee

Submitted: December 2nd 2011Reviewed: May 14th 2012Published: September 27th 2012

DOI: 10.5772/48529

Downloaded: 3648

1. Introduction

This chapter describes widespread methods of model-based fuzzy control systems. The subject of this chapter is a systematic framework for the stability and design of nonlinear fuzzy control systems. We are trying to build a bridge between conventional fuzzy control and classic control theory. By building this bridge, the strong well developed tools of classic control could be used in model-based fuzzy control systems

Model-based fuzzy control, with the possibility of guaranteeing the closed loop stability, is an attractive method for control of nonlinear systems. In recent years, many studies have been devoted to the stability analysis of continuous time or discrete time model based fuzzy control systems (Takagi & Sugeno, 1985; Rhee & Won, 2006; Chen et al.,1993;Wang et al.,1996; Zhao et al.,1996; Tanaka & Wang, 2001; Tanaka et al.,2001). Among such methods, the method of Takagi-Sugeno (Takagi & Sugeno, 1985) has found many applications for modelling complex nonlinear systems (Tanaka & Sano, 1994;Tanaka & Kosaki, 1997;Li et al., 1998). The concept of sector nonlinearity (Kawamoto et al., 1992) provided means for exact approximation of nonlinear systems by fuzzy blending of a few locally linearized subsystems. One important advantage of using such a method for control design is that the closed-loop stability analysis, using the Lyapunov method, becomes easier to apply. Various stability conditions have been proposed for such systems (Tanaka &Wang, 2001), (Ting, 2006), where the existence of a common solution to a set of Lyapunov equations is shown to be sufficient for guaranteeing the closed-loop stability. Some relaxed conditions are also proposed in (Kim & Lee, 2000; Ding et al, 2006; Fang et al.,2006, Tanaka & Ikeda, 1998). Parallel Distributed Compensator (PDC) is a generalization of the state feedback controller to the case of nonlinear systems, using the Takagi-Sugeno fuzzy model (Wang et al.,1996). This method is based on partitioning nonlinear system dynamics into a number of linear subsystems, for which state feedback gains are designed and blended in a fuzzy sense. Takagi-Sugeno model and parallel distributed compensation have been used in many applications successfully (Sugeno & Kang, 1986, Lee et al., 2006, Hong & Langari, 2000, Bonissone et al., 1995). The Linear Matrix Inequality (LMI) technique offers a numerically tractable way to design a PDC controller with objectives such as stability (Wang et al.,1996; Ding et al, 2006; Fang et al., 2006; Tanaka & Sugeno 1992), H control (Lee et al., 2001), H2 control (Lin & Lo, 2003), pole-placement (Jon et al, 1997; Kang & Lee, 1998), and others ‎( Tanaka & Wang, 2001).

2. Takagi-Sugeno fuzzy model

The main idea of the Takagi-Sugeno fuzzy modeling method is to partition the nonlinear system dynamics into several locally linearized subsystems, so that the overall nonlinear behavior of the system can be captured by fuzzy blending of such subsystems. The fuzzy rule associated with the i-th linear subsystem for the continuous fuzzy system and the discrete fuzzy system, can then be defined as

Continuous fuzzy system

Rule i: IF Z1(t)  is Mi1. . . and Zl(t) is MilTHEN{x˙(t)=Aix(t)+Biu(t)y(t)=Cix(t) i=1,2,...,r E1

Discrete Fuzzy System

Rule i: IF Z1(t)  is Mi1. . . and Zl(t) is Mil THEN{x(t+1)=Aix(t)+Biu(t)y(t)=Cix(t) i=1,2,...,rE2

where, x(t)Rnis the state vector, u(t)Rmis the input vector,AiRn×n, BiRn×m,CiRq×n; {z1(t),z2(t),...,zp(t)}are nonlinear functions of the state variables obtained from the original nonlinear equation, and Mij(zi)are the degree of membership of zi(t)in a fuzzy setMij. Whenever there is no ambiguity, the time argument in z(t) is dropped. The overall output, using the fuzzy blend of the linear subsystems, will then be as follows:

Continuous fuzzy system


Discrete Fuzzy System




It is also true, for all t, that


2.1. Building a fuzzy model

There are generally three approaches to build the fuzzy model: "sector nonlinearity," "local approximation," or a combination of the two.

2.1.1. Sector nonlinearity

Figure 1 illustrates the concept of global and local sector nonlinearity. Suppose the original nonlinear system satisfies the sector non-linearity condition (Kawamoto et al., 1992, as cited in Tanaka & Wang, 2001), i.e., the values of nonlinear terms in the state-space equation remain within a sector of hyper-planes passing through the origin. This model guarantees the stability of the original nonlinear system under the control law. A function Φ: R→R is said to be sector [a,c] if for all xϵR, y= Φ(x) lies between b1xandb2x.

Figure 1.

a) Global sector nonlinearity, b) Local sector nonlinearity

Example 1

The well-known nonlinear control benchmark, the ball-and-beam system is commonly used as an illustrative application of various control methods (Wang & Mendel, 1992) depicted in figure 2. Let x1(t) and x2(t) denote the position and the velocity of the ball and let x3(t) and x4(t) denote the angular position and the angular velocity of the beam Then, the system dynamics can be described by the following state-space equation

Figure 2.

The ball and beam system


Where x(t)=[x1(t)x2(t)x3(t)x4(t)]Tand u(t) is torque.

sin(x3)and x1x42are nonlinear terms in the state-space equation. We define z1=sin(x3)andz2=x1x42. Assume x3[π2π2]and x1x4[dd]as the region within which the system will operate. Figure 3 shows that z1(t)=sin(x3(t))and its local sector operating region.The sector [b1, b2] consists of two lines blxl and b2xl, where the slopes are bl = 1 and b2=2π. It follows that

Figure 3.

Formula: Eqn029.wmf>and its local sector

We present sin(x3(t))is represented as follows:


From the property of membership functionsz1=sin(x3(t))=(i=12Mi(z1(t))bi)x3(t), we can obtain the membership functions


Similarly we obtain membership functions associated withM1(z1(t))={z1(t)(2π)Ssin(z1(t))(12π)sin1(z1(t))z1(t)01,otherwise.M2(z1(t))={sin1(z1(t))z1(t)(12π)sin1(z1(t))z1(t)00,otherwise.. Assume z2(t)=x1(t)x4(t)and max(z2(t))=d=α1we have:


The exact TS-fuzzy model-based dynamic system of the ball and beam system can be obtained as following:


The fuzzy model has the following 4 rules:



Rule 1: ifz1(t)isM1andz2(t)isN1Thenx˙(t)=A1x(t)+B1u(t),Rule 2: if z1(t)isM1andz2(t)isN2Then x˙(t)=A2x(t)+B2u(t),Rule 3: if z1(t)isM2andz2(t)isN1Then x˙(t)=A3x(t)+B3u(t),Rule 4: if z1(t)isM2andz2(t)isN2Then x˙(t)=A4x(t)+B4u(t)E15

2.1.2. Local approximation

The original system can be partitioned into subsystems by approximation of nonlinear terms about equilibrium points. This approach can have fewer rules and of course less complexity but it cannot guarantee the stability of the original system under the controller. Usually in this approach, construction of a fuzzy membership function requires knowledge of the behavior of the original system and of course different types of membership functions can be selected.

3. Parallel distributed compensation

Parallel distributed compensation (PDC) is a model-based design procedure introduced in (Wang et al,. 1995). Using the Takagi-Sugeno fuzzy model, a fuzzy combination of the stabilizing state feedback gains, A1=[010000GblDα100010000],A1=[010000GblDα200010000],A1=[010000Gb2Dα100010000],A1=[010000Gb2Dα200010000]B1=B2=B3=B4=B=[0001],z1=sin(x3)andz2=x1x4associated with every linear subsystem is used as the overall state feedback controller. The general structure of the controller is then as


The output of the controller is represented by


The Takagi-Sugeno model and the Parallel Distributed Compensation have the same number of fuzzy rules and use the same membership functions.

4. Stability conditions and control design

4.1. LMI

A variety of problems arising in system and control theory can be reduced to a few standard convex or quasi-convex optimization problems involving linear matrix inequalities (LMIs). Lyapunov published his theory in 1890 and showed that u=i=1rωi(z)Fix(t)i=1rωi=i=1rhi(z)Fix(t).is stable if and only if there exists a positive-definite matrix P such thatddtx(t)=Ax(t). The Lypanov inequality, ATP+PA<0and P>0is a form of an LMI.

An LMI has the form


Where F(x)F0+i=1mxiFi>0,are the given symmetric matrices and FiRn×n,i=0,...,mis the variable and the inequality symbol shows that xRmis positive definite (Boyd, 1994).

4.2. Stability conditions

There are a large number of works on stability conditions and control design of fuzzy systems in the literature. A sufficient stability condition for ensuring stability of PDC was derived by Tanaka and Sugeno (Tanaka & Sugeno, 1990; 1992 ).

By substituting the controller output (15) into the TS model for the continuous fuzzy control (4), we have:


wherex˙(t)=j=1rhi(z(t))hi(z(t))Giix(t)          +2i=1ri<jhi(z(t))hj(z(t)){Gij+Gji2}x(t), Similarly for the discrete fuzzy system we have


Theorem 1: The equilibrium of the continuous fuzzy system (3) with u(t) = 0 is globally asymptotically stable if there exists a common positive definite matrix P such that

x˙(t+1)=j=1rhi(z(t))hi(z(t))Giix(t)          +2i=1ri<jhi(z(t))hj(z(t)){Gij+Gji2}x(t)E23

that is, a common P has to exist for all subsystems.

Theorem 2: The equilibrium of the discrete fuzzy system (4) with u(t) = 0 is globally asymptotically stable i f there exists a common positive definite matrix P such that


that is, a common P has to exist for all subsystems.

The stability of the closed loop system can be derived by using theorem 1 and 2.

Theorem 3: The equilibrium of the continuous fuzzy control system described by (18) is globally asymptotically stable if there exists a common positive definite matrix P such that


Theorem 4: The equilibrium of the discrete fuzzy control system described by (20) is globally asymptotically stable if there exists a common positive definite matrix P such that

i<j s.t. hihjϕE27

4.3. Stable controller design

By using the following conditions, the solution of the LMI problem for continuous and discrete fuzzy systems gives us the state feedback gains Fi and the matrix P (if the problem is solvable).

Consider a new variable i<j s.t. hihjϕthen the stable fuzzy controller design problem is:

Continuous fuzzy system

Find X=P1andX>0, Mi

XAiTAiX+MiTBiT+BiMi>0,XAiTAiXXAjTAjX           +MjTBiT+BiMj+MiTBjT+BjMi0.E30

The conditions (27) and (28) gives us a positive definite matrix X=P1i<j s.t. hihjϕand X(or that there is no solution). From the solution MiandX, a common P and the feedback gains can be found as:


Similarly for a discrete fuzzy system the design problem is

Find P=X1,Fi=MiX1andX>0, Mi


4.4. Decay rate

Decay rate is associated with the speed of response. The decay rate fuzzy controller design helps to find feedback gains that provide better setteling time (Tanaka et al,. 1996; 1998a; 1998b).

Continuous fuzzy system: The condition that X(AiXBiMi)TX1(AiXBiMi)>0,X14X(AiXBiMi+AjXBjMi)TX1×(AiXBiMj+AjXBjMi)X0.(Ichikawa et al, 1993, as cited in Tanaka & Wang, 2001) for all V˙(x(t))2αV(x(t))can be written as




Therefore, by solving the following generalized eigenvalue minimization problem in X, the largest lower bound on the decay rate that can be found by using a quadratic Lyapunov function:

maximize Gij=AiBiFi,α> 0 andi<j s.t. hihjϕsubject to

X>0,XAiT+AiX+MiTBiT+BiMi2αX>0,XAiTAiXXAjTAjX+MjTBjT+BiMj          +MiTBjT+BjMi4αX>0,E36

Similarly for a discrete fuzzy system:

The condition that i<j s.t. hihjφ,whereX=P1,    Mi=FiX.(Ichikawa et al, 1993, as cited in Tanaka & Wang, 2001) for all ΔV(x(t))(α21)V(x(t))can be written as


The generalized eigenvalue minimization can be found in (Tanaka & Wang, 2001).

4.5. Constraint on control

Theorem 5: Assume that the initial condition x(0) is known. The constraint i<j s.t. hihjϕandα<1is satisfied at all times u(t)2μif the LMIs


Hold, where [1x(0)Tx(0)X]0[XMiTMiμ2I]0andX=P1.

The above LMI design conditions depend on the initial states. Thus, if the initial states  Mi=FiXchange, this means that the feedback gains Fi must be again determined. To overcome this disadvantage, modified LMI constraints on the control input have been developed, where x(0)is unknown but the upper bound x(0)of ϕis known, i.e.,x(t).

Theorem 6: Assume thatx(t)ϕ, where x(0) is unknown but the upper bound x(t)ϕis known. Then,




Proofs of theorem 1 and 2 are given in (Tanaka & Wang, 2001)

4.6. Performance-oriented parallel distributed compensation

In the modified PDC proposed in (Seidi & Markazi, 2011), unlike the conventional PDC, state feedback gains associated with every linear subsystem, are not assumed fixed. Instead, based on some pre-specified performance criteria, several feedback gains are designed and used for every subsystem. The overall gain associated with each of the subsystems, is then determined by a fuzzy blending of such gains, so that a better closed-loop performance can be achieved. The required membership functions are chosen based on some pre-specified performance indices, for example, a faster response or a smaller control input. In general, the rest of the method for calculating the overall state feedback gain remains similar to the conventional PDC method, as in (14) and (15). Figure 4, depicts the general framework for the proposed method, through which and depending on various performance criteria, different characteristics for the controller can be specified. For example, two different feedback gains could be designed for a typical subsystem; one providing a lower control input with a longer settling time response, and the other a faster response but with a larger control input. The idea is then to select the overall feedback gain for this subsystem as a weighted sum of such gains, where the weights are appropriately adjusted, in a fuzzy sense, during the time evolution of the system response, so that as a whole, a faster response with a lower control input can be achieved. For this purpose, when the magnitude of the control input becomes large, the relative weight of the first feedback gain is increased, so that the magnitude of the control input is kept within the permissible limits. On the other hand, when the control input is well below the permissible limit, the weight of the second feedback gain is increased, for a faster response. The dynamics of the resulting closed-loop control system can be analyzed as follows:

Consider the following Takagi–Sugeno model of the plant


The following structure is proposed for the fuzzy controller rules


Wherei th rule: IfZ1(t)is Mi1andZ2(t)isMi2,.......,Zp(t)isMip,J(t)isHi1,....andJ(t)isHiqthenui(t)={n=1qmin(J(t))Kin}x(t), i=1,2,...,ris the number of gain coefficients in the ith subsystem, qiis the relevant membership degree for J(t), minis the nth state feedback gain associated with the ith subsystem, Kinis the n th membership function for J(t), defined in the ith rule. Here Hiqis a term depicting a selected performance index, for instance, if one wants to limit the magnitude of the control signalJ(t), thenu(t). Where the control input generated by the PDC controller is in the form of


Figure 4.

General methodology in the proposed PDC method

Lemma: The fuzzy control system (39), with the control strategy (41) is globally, asymptotically stable, if there exists a common positive definite matrix P such that


where GiinTP+PGiin<0(Gijn+Gjin2)TP+P(Gijn+Gjin2)0

i<j,    hihjϕ,E46

Example 2

Consider a single link robot with flexible joint as in Figure 5. This benchmark problem is introduced in (Spong et al., 1987).

Figure 5.

A single link robot with a flexible joint

The state space equations for the system of Figure 4 are


In order to apply the PDC methodology, the fuzzy Takagi-Sugeno Model is developed first (Seidi & Markazi, 2008). The nonlinear expression{x˙1=x3(t)x˙2=x4(t)x˙3=1I(k(x2(t)x1(t))mgLsin(x1(t)))x˙4=1J(u(t)k(x2(t)x1(t))), forZ=sin(x1(t)), can be expressed as


Where, z=sin(x1(t))=(i=12Mi(z)bi)x1(t)and, hence, the membership functions for b1=1,b2=0are obtained as


The resulting fuzzy model would then have the following fuzzy rules:

M1(z)={zSin1z,  z(t)0        1,   OtherwiseM2(z)={Sin1zzSin1z,  z(t)01,   OtherwiseE50


Rule 1:Ifz(t)isM1(z),thenx˙(t)=A1x(t)+B1u(t)Rule 1:Ifz(t)isM2(z),thenx˙(t)=A2x(t)+B2u(t)

AssumeB1=B2=B=[0,0,0,1]T., k=100Nm/radand other parameters are assumed unity then we have


The final output of the controller is

Control Rule1:Ifz(t)isM1(z),thenu(t)=F1x(t)Control Rule2:Ifz(t)isM2(z),thenu(t)=F2x(t)E55

Case 1: Stable controller design

Using conditions (27) and (28) the stable controller can be obtained by solving below conditions


Using the MATLAB LMI Control Toolbox we obtain

X>0 [X A1A1X+M1TBT+B M1]>0,[X A2A2X+M2TBT+B M2]>0,[X A1TA1XX A2TA2X+M2TBT+B M2+M1TBT+B M1]>0E57

Figures 6 and 7 show the response of the system and control effort, respectively.

Case 2: The decay rate

Using conditions (31) and (32) the stable controller can be obtained by solving the conditions:

Figure 6.

Response of flexible joint robots x1(t), case 1.

Figure 7.

Control input for flexible joint robots, case 1.


Considering [X A1TA1X+M1T BT+B M12αX]>0[X A2TA2X+M2TBT+B M22αX]>0[X A1TA1XX A2TA2X+M2TBT+B M2+M1TBT+B M14αX]>0and by using the MATLAB LMI Control Toolbox we obtain:


Figures 8 and 9 show the response of the system and control effort, respectively.

Figure 8.

Response of flexible joint robots x1(t), case 2.

Figure 9.

Control input for flexible joint robots, case 2.

Case 3: The decay rate with the constraint on the input

We design a stable fuzzy controller by considering the decay rate and the constraint on the control input. The design problem of the FJR is defined as follows:


F1=[4108.86545.21271.3127.77]F2=[4066.96502.61261.7127.1]P=[36.5087 24.01406.21350.335224.014030.13416.32230.50136.21356.32231.42600.09950.33520.50130.0995 0.0099]E60


WhereX>0[XA1TA1X+M1TBT+BM12αX]>0[XA2TA2X+M2TBT+BM22αX]>0[XA1TA1XXA2TA2X+M2TBT           +BM2+M1TBT+BM14αX]>0[XM1TM1μ2I]>0[XM2TM2μ2I]>0[Xϕ2I]>0, X=P1, Mi=FiX,


Using the MATLAB LMI toolbox to solve the LMI conditions (50), we can get the positive definite matrix and a set of gains (51), that make the system stable.


Figures 10 and 11 show the response of the system and control effort, respectively.

Figure 10.

System responses of the single-link flexible joint, case 3.

Figure 11.

Control input for flexible joint robots, case 3.

Case 4: Performance-oriented parallel distributed compensation

The following stabilizing feedback gains are chosen using the pole placement method, so that P=[0.73010.324860.0967940.00345520.324860.554830.106160.0102090.0967940.106160.0230490.00171390.00345520.0102090.00171390.00023565]F1=[327.571745261.8657.475]F2=[356.051739.2259.7757.5]and K11produce large magnitude inputs for subsystems 1 and 2, respectively, and K21and K22induce low magnitude inputs for those subsystems. In particular,


The required simple membership functions are selected as in Figure 12, so that, with a decrease in the corresponding plant input, in subsystems 1 and 2 respectively, the overall feedback gains come closer to K11=[6667.24411.91052.492.6]K12=[-33.3211413.7191.6351.2]K21=[6658.74332.41025.491.1]K22=[72.31389.8189.650.6]andK11, and with an increase in the corresponding control input respectively, the overall feedback gains come closer to K21andK21. Now, the fuzzy rules for the controller are constructed as follows:

Rule 1: If K22is z(t)and M1(z)is "small" then |u(t)|

Rule 2: If u(t)=K11x(t)is z(t)and M1(z)is "large" then |u(t)|

Rule 3: If u(t)=K12x(t)is z(t)and M2(z)is "small" then |u(t)|

Rule 4: If u(t)=K21x(t)is z(t)and M2(z)is "large" then |u(t)|

Figure 12.

Membership functions for the control effort in the flexible joint robots.

A common positive definite matrix, P, satisfying the stability conditions (42) is obtained by solving the LMI problems:


Applying a unit step reference signal forP=104×[121710158582558.563.525158588624.41458.4105.362558.51458.4702.2442.52963.525105.3642.5295.0962], the response history and the corresponding control input are shown in Figures (13) and (14), respectively. Simulation results are investigated for the following three controllers:

Figure 13.

Response of flexible joint robot x1(t), case 4.

Figure 14.

Control input for flexible joint robot, case 4.

  1. A PDC controller with feedback gains x1(t)and K11providing a high speed response, and with possible high control inputs (HPDC controller).

  2. A PDC controller with feedback gains K21and K22providing a low speed response, and with a lower control input, as compared with the HPDC case (LPDC controller).

  3. Proposed modified PDC controller, providing a fast response, yet with an acceptable level of control input (NPDC controller).

It is observed that the new controller provides a settling time similar to the HPDC case, with a much lower magnitude for the control input.

3. Conclusion

This chapter deals with approximation of the nonlinear system using Takagi-Sugeno (T-S) models with linear models as rule consequences and a construction procedure of T-S models. Also, the stability conditions and stabilizing control design of parallel distributed compensation (PDC) are discussed. It is seen that PDC a linear control method can be used to control the nonlinear system. Moreover, the stability analysis and control design problems for both continuous and discrete fuzz control systems can be transformed to linear matrix inequality (LMI) problems and they can be solved efficiently by convex programming techniques for LMIs. Design examples demonstrate the effectiveness of the LMI-based designs.

© 2012 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Morteza Seidi, Marzieh Hajiaghamemar and Bruce Segee (September 27th 2012). Fuzzy Control Systems: LMI-Based Design, Fuzzy Controllers - Recent Advances in Theory and Applications, Sohail Iqbal, Nora Boumella and Juan Carlos Figueroa Garcia, IntechOpen, DOI: 10.5772/48529. Available from:

chapter statistics

3648total chapter downloads

5Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

New Results on Robust H∞ Filter for Uncertain Fuzzy Descriptor Systems

By Wudhichai Assawinchaichote

Related Book

First chapter

New Strategy to Approach the Inverse Kinematics Model for Manipulators with Rotational Joints

By José-Emilio Vargas-Soto, Efren Gorrostieta-Hurtado, Saúl Tovar- Arriaga, Jesús-Carlos Pedraza-Ortega and Juan-Manuel Ramos- Arreguín

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us