International Journal of Systems Science
International Journal of Systems Science
To cite this article: Anouar Bouazza , Anis Sakly & Mohamed Benrejeb (2013) Order reduction of complex systems described
by TSK fuzzy models based on singular perturbations method, International Journal of Systems Science, 44:3, 442-449, DOI:
10.1080/00207721.2011.602480
Taylor & Francis makes every effort to ensure the accuracy of all the information (the “Content”) contained
in the publications on our platform. However, Taylor & Francis, our agents, and our licensors make no
representations or warranties whatsoever as to the accuracy, completeness, or suitability for any purpose of the
Content. Any opinions and views expressed in this publication are the opinions and views of the authors, and
are not the views of or endorsed by Taylor & Francis. The accuracy of the Content should not be relied upon and
should be independently verified with primary sources of information. Taylor and Francis shall not be liable for
any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever
or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of
the Content.
This article may be used for research, teaching, and private study purposes. Any substantial or systematic
reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any
form to anyone is expressly forbidden. Terms & Conditions of access and use can be found at http://
www.tandfonline.com/page/terms-and-conditions
                                                                     International Journal of Systems Science
                                                                     Vol. 44, No. 3, March 2013, 442–449
                                                                             Order reduction of complex systems described by TSK fuzzy models based on singular
                                                                                                            perturbations method
                                                                                                             Anouar Bouazza*, Anis Sakly and Mohamed Benrejeb
                                                                                      Unite´ de recherche LARA Automatique, Ecole Nationale d’Inge´nieurs de Tunis, BP 37 Le Belve´de`re,
                                                                                                                            Tunis 1002, Tunisia
                                                                                                           (Received 13 November 2008; final version received 22 June 2011)
                                                                              This article deals with the concept of order reduction of linear complex systems described by TSK fuzzy models.
                                                                              The use of singular perturbations technique for process modelling and the choice of Benrejeb arrow form
                                                                              characteristic matrix provide the decoupling of dynamics of linear systems in the continuous and discrete case.
                                                                              An original contribution is to apply nonconventional TSK fuzzy approach on a direct current motor in order, on
Downloaded by [Florida State University] at 00:29 29 December 2014
                                                                              one hand, to highlight the order reduction of this system presented locally in the form of linear models and, on
                                                                              the other hand, to show the efficiency of the proposed approaches.
                                                                              Keywords: order reduction technique; method of singular perturbations; two-time scale system; discretisation of
                                                                              continuous systems; TSK fuzzy models; Benrejeb arrow form matrix
                                                                     Borne 2006) and can be written in the following form:                              Then, the matrix P is considered
                                                                                                                                                              2                                                 3
                                                                            x_ 1 ðtÞ     A11       A12        x1 ðtÞ     B1                                        0                    1                1       0
                                                                                      ¼                                 þ    uðtÞ              ð1Þ               6 1
                                                                            x_ 2 ðtÞ     A21       A22        x2 ðtÞ    B2                                     6                       2              n1     07 7
                                                                                                                                                                 6 2                                                 .. 7
                                                                                                                                                             P¼6
                                                                                                                                                                 6 1                    22             2n1     .77           ð5Þ
                                                                                                              x1 ðtÞ                                           6 ..                     ..          ..                7
                                                                                          yðtÞ ¼ C1         
                                                                                                           C2                                  ð2Þ               4 .                       .              .       15
                                                                                                                 x2 ðtÞ
                                                                                                                                                                   n1
                                                                                                                                                                    1                   n1
                                                                                                                                                                                         2               n1
                                                                                                                                                                                                              n1    1
                                                                     with
                                                                                                                                                     The transformation
                                                                              x1 ðtÞ 2 Rn1     slow vectors of the system
                                                                              x1 ðtÞ 2 Rn2     fast vectors of the system                                                            z ¼ P1 x                                      ð6Þ
                                                                                                                                                     deals up to the new description:
                                                                                 Aij              C2                  B2
                                                                        Aij ¼       ;    C2 ¼      ;      B2 ¼        ;       n1 þ n2 ¼ n                                                 ~
                                                                                                                                                                                     z_ ¼ AðÞz                                  ð7Þ
Downloaded by [Florida State University] at 00:29 29 December 2014
                                                                               2                                                           3
                                                                                     0            1            0                 0
                                                                             6       0            0            1                 0      7
                                                                             6       ..           ..           ..        ..         ..     7         3. Reduced TSK fuzzy models
                                                                      AðÞ ¼ 6
                                                                             6        .            .            .           .        .     7 ð3Þ
                                                                                                                                           7
                                                                             4                                                             5         The TSK fuzzy approach is a multi-model technique
                                                                                     0          0            0                 1
                                                                                   a0 ðÞ    a1 ðÞ      a2 ðÞ          an1 ðÞ             used for modelling and control of complex systems.
                                                                                                                                                     The idea is to decompose the process in linear
                                                                     i ðÞ is a coefficient of the instantaneous characteristic                     submodels in space state correlated by membership
                                                                     polynomial of the matrix AðÞ such that                                         functions. The TSK fuzzy approach is a very interest-
                                                                                                                                                     ing approach for the description of nonlinear systems.
                                                                                                               X
                                                                                                               n1
                                                                                                                                                         The TSK fuzzy system under duress the rule, Ri, for
                                                                                          Pð, Þ ¼ n þ              i ðÞi                 ð4Þ
                                                                                                                                                     i ¼ 1, 2, . . . , r, is given by Benrejeb et al. (2006) and
                                                                                                                i¼0
                                                                                                                                                     Benrejeb, Sakly, Ben Othman, and Borne (2008)
                                                                     with 1 , 2 , . . . , n1 : ðn  1Þ being distinct real or
                                                                     complex numbers.                                                                         If x1 ðtÞ is Mi1 and . . . and xn ðtÞ is Min
                                                                                                                                                                         (
                                                                         The        suitable         choice   of    the    vectors                                         _ ¼ Ai xðtÞ þ Bi uðtÞ
                                                                                                                                                                           xðtÞ                                                    ð12Þ
                                                                      ¼ ½1 , 2 , . . . , n1  has high consequences for mod-                                Then
                                                                     elling and analysis of the reduced order TSK model.                                                   yðtÞ ¼ Ci xðtÞ
                                                                     444                                                      A. Bouazza et al.
                                                                     with Mij , j ¼ 1, 2, . . . , n, the ith fuzzy subset relating to    The inferred model of the reduced TSK fuzzy models
                                                                     the state Vector xi, x 2 Rn , u 2 R.                                can be described by Benrejeb and Abdelkrim (2003)
                                                                         The final outputs of the fuzzy model are inferred as                                     Xr
                                                                     follows:                                                                                A~ ¼    hi A~ i            ð22Þ
                                                                                                                                                                            i¼1
                                                                                         X  r
                                                                                  _ ¼
                                                                                  xðtÞ          hi ðxðtÞÞðAi xðtÞ þ Bi uðtÞÞ     ð13Þ
                                                                                                                                             A~ i ¼ A22,i    A21,i A1
                                                                                                                                                                     11,i   A12,i   i ¼ 1, 2, . . . , r   ð23Þ
                                                                                         i¼1
                                                                                                X
                                                                                                r
                                                                                       yðtÞ ¼         hi ðxðtÞÞCi xðtÞ          ð14Þ     4. Discretisation of the continuous systems
                                                                                                i¼1
                                                                                                                                         The discretisation concerns the process of transferring
                                                                                                     wi ðxðtÞÞ                           continuous models and equations into discrete coun-
                                                                                       hi ðxðtÞÞ ¼ Pr                           ð15Þ
                                                                                                    i¼1 wi ðxðtÞÞ                        terparts. This process is usually carried out as a first
                                                                     The hi verifies the property of a sum convex, which is              step towards making them suitable for numerical
                                                                     to say that                                                         evaluation and implementation on digital computers.
                                                                                                Xr                                           Among the methods of discretisation, the homo-
                                                                                       hj  0;        h ¼1             ð16Þ
Downloaded by [Florida State University] at 00:29 29 December 2014
                                                                                                                                                              x_~ ¼ A~ i x~ þ Bu
                                                                                                                                                                              ~    i ¼ ð1, 2Þ             ð35Þ
                                                                     5. Exploitation of the TKS fuzzy approach for order                 with
                                                                        reduction of DC motor
                                                                                                                                                       A~ i ¼ A22,i  A21,i A1
                                                                                                                                                                             11,i A12,i      i ¼ ð1, 2Þ   ð36Þ
                                                                     The singular perturbation method presupposes the
                                                                     knowledge of both slow and fast state vectors.                      such as
                                                                     Generally, the partitioning is based on some physical                                                        
                                                                     properties of the process.                                                              ~      0:1321 0:0770
                                                                                                                                                             A1 ¼                                         ð37Þ
                                                                         An example of a very important application is that                                         0:0103 0:0066
                                                                     of electrical machinery rotating.
                                                                                                                                                                                  
                                                                         Indeed, they are modelled by two types of                                                  0:1321 0:0770
                                                                                                                                                             ~
                                                                                                                                                             A2 ¼                                         ð38Þ
Downloaded by [Florida State University] at 00:29 29 December 2014
                                                                     defined and, now, a fuzzy controller is designed for                  If ðx1 is M2 Þ and ðx2 is M2 Þ and ðx3 is M2 Þ
                                                                     them.                                                                                                                  ð50Þ
                                                                                                                                                Then u ¼ K2 x
                                                                         For the global third order system, the following
                                                                     fuzzy rules are defined.                                          with Ki ði ¼ 1, 2Þ: Vectors of gains, defined from the
                                                                                                                                       quadratic optimal control relating to the global system
                                                                          . For the open loop’s system considered in
                                                                                                                                       described locally in the form of two linear submodels.
                                                                     autonomous regime:
                                                                                                                                            . For the global system:
                                                                                                                                                                         
                                                                           If ðx1 is M1 Þ and ðx2 is M1 Þ and ðx3 is M1 Þ                      K1 ¼ 10:9898 0:3419 1:5019                   ð51Þ
                                                                                                                             ð45Þ
                                                                                     _ ¼ F1 xðtÞ
                                                                              Then xðtÞ                                                                                      
                                                                                                                                                   K2 ¼ 11:4055 0:2593 1:5260               ð52Þ
                                                                           If ðx1 is M2 Þ and ðx2 is M2 Þ and ðx3 is M2 Þ
                                                                                                                             ð46Þ
                                                                                     _ ¼ F2 xðtÞ
                                                                              Then xðtÞ                                                     . For the reduced system:
                                                                                                                                                                                
                                                                                                                                                       K~ 1 ¼ 0:3402    1:4392              ð53Þ
Downloaded by [Florida State University] at 00:29 29 December 2014
                                                                     Figure 3. Responses of the continuous and discrete global system before introducing the TSK fuzzy controller.
                                                                                                              International Journal of Systems Science                                     447
Downloaded by [Florida State University] at 00:29 29 December 2014
Figure 4. Responses of the continuous and discrete global system after introducing the TSK fuzzy controller.
Figure 5. Responses of the continuous and discrete reduced system before introducing the TSK fuzzy controller.
                                                                     In the same way, the fuzzy rules for the second order                According to Figures 4 and 6, the robustness of the
                                                                     reduced system (37), (38) are defined.                           proposed models is examined; the responses of
                                                                         Similarly, vectors of gains and fuzzy rules for the          the continuous and discrete global system and the
                                                                     global (40), (41) and reduced (43) systems are defined           responses of the continuous and discrete reduced
                                                                     in the discrete case (Bouazza et al. 2008).                      system are similar. The TSK fuzzy approach conserves
                                                                         The membership functions relating to global and              the state variables of the initial system and, therefore,
                                                                     reduced systems are represented as shown in Figure 2.            the real physical meaning.
                                                                         The simulation of the global system gives the
                                                                     following responses (Figures 3 and     4). The initial
                                                                                                       0:5
                                                                     conditions are such that xð0Þ ¼ 0:5    .                         6. Conclusion
                                                                                                        0:5
                                                                         The simulation of the reduced system gives the               This work is interested in TSK fuzzy approach for the
                                                                     following responses (Figures 5 and 6). The initial             modelling and control of nonlinear systems presented
                                                                     conditions are such that xð0Þ ¼ 0:5
                                                                                                     0:5 .                            locally in the form of linear submodels.
                                                                     448                                                 A. Bouazza et al.
Downloaded by [Florida State University] at 00:29 29 December 2014
Figure 6. Responses of the continuous and discrete reduced system after introducing the TSK fuzzy controller.
                                                                       Analysis and Generalized Dimensional Analysis’,               in 15th Mediterranean Conference on Control and
                                                                       International Journal for Computer-Aided Engineering and      Automation (MED’07), Athens, Greece, pp. 1–6.
                                                                       Software, 25, 443–463.                                      Lubuma, J.M.–S., and Patidar, K.C. (2009), ‘Towards the
                                                                     Benrejeb, M., and Abdelkrim, M.N. (2003), ‘On Order             Implementation of the Singular Function Method for
                                                                       Reduction and Stabilization of TSK Nonlinear Fuzzy            Singular Perturbation Problems’, Applied Mathematics and
                                                                       Models by Using Arrow Form Matrix’, Systems Analysis          Computation, 209, 68–74.
                                                                       Modeling Simulation, 43, 977–991.                           Sampaio, R., and Soize, C. (2006), ‘Remarks on the
                                                                     Benrejeb, M., Sakly, A., Ben Othman, K., and Borne, P.          Efficiency of POD and KL Methods for Model
                                                                       (2008), ‘Choice of Conjunctive Operator of TSK Fuzzy          Reduction in Nonlinear Dynamics of Continuous
                                                                       Systems and Stability Domain Study’, Mathematics and          Systems’, International Journal for Numerical Methods in
                                                                       Computers in Simulation, 76, 410–421.                         Engineering, 72, 22–45.
                                                                     Benrejeb, M., Soudani, D., Sakly, A., and Borne, P. (2006),   Satakshi, Mukherjee, S., and Mittal, R.C. (2005), ‘Order
                                                                       ‘New Discrete Tanaka Sugeno Kang Fuzzy Systems                Reduction of Linear Discrete Systems Using a
                                                                       Characterization and Stability Domain’, International         Genetic Algorithm’, Applied Mathematical Modelling, 29,
                                                                       Journal of Computers, Communications and Control, 1,          565–578.
                                                                       9–19.                                                       Shafrir, I. (2004), ‘On a Class of Singular Perturbation
                                                                     Bouazza, A., Sakly A., and Benrejeb, M. (2008), ‘On Order       Problems’, in Handbook of Differential Equations:
Downloaded by [Florida State University] at 00:29 29 December 2014
                                                                       Reduction of Discrete Complex Systems Described by            Stationary Partial Differential Equations (Vol. 1,
                                                                       TSK Fuzzy Models’, in 9th International Conference on         Chap. 5), eds. M. Chipot, and P. Quittner, Amsterdam:
                                                                       Sciences and Techniques of Automatic Control and              Elsevier, pp. 297–383.
                                                                       Computer Engineering (STA’08), Sousse, Tunisia.             Taniguchi, T., and Tanaka, K. (2000), ‘Rule Reduction and
                                                                     Coluccio, L., Eisinberg, A., and Fedele, G. (2007),             Robust Control of Generalized Takagi–Sugeno Fuzzy
                                                                       ‘A Prony-like Polynomial-based Approach to Model              Systems’, Journal of Advanced Computational Intelligence
                                                                       Order Reduction’, in 15th Mediterranean Conference on         and Intelligent Informatics, 4, 373–379.
                                                                       Control and Automation (MED’07), Athens, Greece,            Wang, H.O., Tanaka, K., and Griffin, M. (1996),
                                                                       pp. 1–6.                                                      ‘An Approach to Fuzzy Control of Non-linear Systems,
                                                                     Frangos, M., and Jaimoukha, I.M. (2007), ‘Rational              Stability and Design Issues’, IEEE Transactions of Fuzzy
                                                                       Interpolation: Modified Rational Arnoldi Algorithm and        Systems, 4, 14–23.
                                                                       Arnoldi-like Equations’, in 46th IEEE Conference on         Zhang, B., and Li, Y.F. (2005), ‘An Efficient Method for
                                                                       Decision and Control, New Orleans, LA, USA.                   Dynamic Calibration and 3D Reconstruction Using
                                                                     Kherraz, K. (2007), ‘Application of the Genetic Algorithms      Homographic Transformation’, Sensors and Actuators A:
                                                                       to the Reduction of Order of an Engine with DC Current’,      Physical, 119, 349–357.