Sliding Mode Control for Singularly Perturbed Systems Using Accurate Reduced Model

In order to deal with unmodeled dynamics in large vehicle systems, which have an ill condition of the state matrix, the use of model order reduction methods is a good approach. This article presents a new construction of the sliding mode controller for singularly perturbed systems. The controller design is based on a linear diagonal transformation of the singularly perturbed model. Furthermore, the use of a single sliding mode controller designed for the slow component of the diagonalized system is investigated. Simulation results indicate the performance improvement of the proposed controllers.


Introduction
The design requirements for high precision and high maneuverable missile, need it to be slender and long. Using light weight material and thin wall structure to design the motor for these long and thin missiles induced elasticity in the body. The rigid body model does not satisfy the needs of precision controller. Singular Perturbation analysis is used by many researchers to overcome such a problem [1] [2] [3] [4] where the fast state variables are the elastic forces and their time derivatives. Singular Perturbation theory first appeared in mathematical literature then in control literature. Many books and research papers are found in this direction, the most famous are the books by Kokotovic [3] and Nadiu [4].
The use of sliding mode control for singularly perturbed systems is explored. A number of articles have been published on this topic. Slotine [5] [6] proposed one of the earliest attempts by applying his sliding mode tracking control algorithm to a robot in singular perturbed form. Heck [7] also introduced sliding mode control for a singular perturbed system. Heck [7] used two sliding mode controllers, one for each of the fast and slow subsystems. Alvarez-Gallegos [1] introduced an approach to find a maximum value of the perturbation parameter ε* when two sliding mode controllers are used for the full order singularly perturbed system. Alvarez-Gallegos [1] uses the Lyapunov approach that was introduced by Kokotovic [3]. These articles propose to design two separate sliding mode controllers for the fast and slow subsystems. The singularly perturbed system is separated into slow and fast modes by using quasi steady state techniques. The result is separate dynamic models for the fast and slow states, the literature proposed designing sliding mode controllers for the slow and fast modes. The difficulty with this method is that the quasi steady state model is only an approximation of the dynamics of the true full order systems.
We propose, in this article, to transform the full order singularly perturbed system into block diagonal form [8]. Based on the transformed systems we construct the slow and fast sliding mode controllers. We compare the performance of the sliding mode controller designed based on the quasi steady state approximation of the dynamics to the sliding mode controllers designed based on the block diagonal transformation. It is shown that sliding mode control design based on the transformed system has superior performance.
We also investigate the design of a single sliding mode controller for the full order system. This single sliding mode controller is designed based on the slow dynamics of the transformed system, the fast dynamics are considered as unmodeled high frequency dynamics.
The outline of this paper is as follows. Section 2 describes the singularly perturbed system and the appropriate transformations. Section 3 illustrates the sliding mode controller design for both reduced subsystems obtained using the diagonalization method. Section 4 describes the design of a single sliding mode controller for the full order system. Section 5 shows the simulation results.

Problem Formulations
The system considered in this paper may be represented in the following form, 11 12 1 and ε is a small positive parameter. This is a form for standard linear time invariant singularly perturbed systems. Due to the high dimensionality of the full order system, singular perturbation theory is useful in decoupling the system into two reduced order subsystems considering the fact that fast modes are important only during a short initial period. After that period they are negligible and the behavior of the system can be described by its slow modes. The model with fast modes neglected is called the quasi steady state model (or zeroth order model) [3] [9]. Neglecting the fast modes is equivalent to assuming that they are infinitely fast, that is setting the perturbation parameter ε in (2) to zero. Without fast modes the system described by equation (1) and (2) are reduced to, where , x z   and u  are tracking vector error. Thus we can get the following reduced system [9] [10],  (1) and (2) can be defined as, The eigenvalues of A 0 are good approximations of the slow modes of the full order system for sufficiently small ε [1] [11]. It is possible to design a controller based on the quasi steady state model to stabilize the full order system for sufficiently small ε. However, for real systems the perturbation parameter ε is not zero. Validation problems for the zero order model arise where for a certain value of ε, the zero order model can be stable while the slow modes of the original system are unstable [1] for this value of perturbation parameter ε. To overcome this problem, we introduce the diagonalization method [11]. This can be achieved by using two stage linear transformations. The first stage is to use the transformation, This transforms the system described by Equations (1) and (2) into upper block diagonal form as shown in Equation (8) ( ) 11 12 12 1 We required that matrix L satisfy ( ) 21 22 Then using the second stage of transformation To obtain, Thus the system has the required block diagonal or decoupled form, where the slow and the fast variables s x and f x can be solved independently of each other. The eigenvalues of A s are exactly the same as the slow poles of the original system. Similarly, the eigenvalues of A f are the same as the fast poles of the full order model. The matrix L is dependent on the perturbation parameter, ε.
Letting 0 ε → gives us Based on the pervious analysis the diagonalization method must be used instead of quasi steady state to decouple the full order system especially when the quasi steady state fails to describe the stability properties of the full order model correctly.
Equation (9) and (13) is an asymmetric Riccati equation which need special procedure to be solved reader can found more detail in the literature [9] [10].

Sliding Mode Controller Design
In this section we will design sliding mode controllers based on the reduced subsystems obtained from transforming the systems into block diagonal form. This is different from the methods proposed in the literature [6] [12] [13] [14] [15] that design the sliding mode controllers based on the quasi steady state approximation of the slow and fast modes. The design of the sliding mode control for the full order system is done in two stages. First, design a slow controller based on the slow subsystem. Second, design fast controllers based on fast subsystem, then, use the composite of these controllers to control the full order model.

Slow Sliding Mode Control Design
The slow sliding mode controller is designed using the slow reduced order model in (14). The equivalent control method [3] Using this control law in the slow subsystem obtained from the block diagonalized reduced system (14) we then obtained, To guarantee local stability for this controller we must ensure that the eigenvalues of the equivalent system given by equation (16) have negative real parts. A Lyapunov function is used to determine the discontinuous control law (u N ) that will satisfy the reaching condition [7] [16]. This will ensure the global asymptotic stability of the closed loop system. The control law will have the form, where u N is a discontinuous control action that drives the state to the sliding surface.

Fast Sliding Mode Controller Design
The same argument used to design the slow sliding mode controller can be used to design the fast sliding mode controller. Define the linear fast switching surface The equivalent control method is used to determine the control law for the fast subsystem given by Equation (14). The control law is given by, which ensures global asymptotic stability for fast subsystem. The control law for the full order model will be the composite of the slow and fast controllers as follows, Since the eigenvalues of the reduced order subsystems is the same as the full order model this means that the linear transformation preserves the stability condition for the closed loop system, and taking the stability analysis proposed by Kokotovic [3] into account, we could say that the full order model is globally stable for * 0, ε ε   ∈   almost everywhere given that * 0 ε > .

Slow Sliding Mode Controller for the Full Order Model
In this section, a single sliding mode controller is proposed to control the full order model. The idea is to use the robust properties of the sliding mode controller to counter the effects of the unmodeled high frequency dynamics. We will consider that the unmodeled high frequency dynamics are represented by the fast subsystem. The main assumption is that the matrix A 22 is stable while the full order model given by Equation (1) and Equation (2) where λ is the bandwidth of the system [5] and is taken to be smaller than the lowest frequency of the fast subsystem. The control law is chosen based on an equivalent control method and has the following form, An advantage of using the proposed controller is that the control law in this case will depends only on the slow state. As a result there is no need for measurement of the fast state, which is usually difficult to measure.

Illustrative Example
Sliding mode controller used usually in missiles and spacecraft [15] [19]. However, in analyzing the dynamics of flexible missiles, such as intermediate-range ballistic missiles or intercontinental ballistic missiles, it is convenient to use singular perturbation analysis. Since all detailed information on missiles is classified, we will use a simple example which suffers from flexibility to show the importance of our study. Consider a magnetic tape control system [8] given by This can be put in standard singularly perturbed form as, 11 12 1 The first 2 × 2 transformation matrix L has to satisfy, ( ) The fast subsystem can be described by, We notice that the eigenvalues of the reduced system are the same as the eigenvalues of the full order model. We design two sliding mode controllers based on the slow and fast subsystems. For the slow subsystem the control law will be,  Figure 1 illustrates the closed loop time response of the full order system. It is observed that the absolute value of the error between the desired response and zeroth model response is larger than the absolute error between the desired response and the proposed model response based on diagonalization. The difference in the error can be four times larger.
Since the fast subsystem is stable then based on the decoupled reduced order system we will design a sliding mode controller for the slow subsystem. The fast subsystem will be considered as unmodeled high frequency dynamics. Choosing the surface parameter as,  Two SMC Single SMC a single sliding mode controller. It is shown that it is possible to use one sliding mode controller to control the full order model for singularly perturbed system given that the fast subsystem is stable and we choose the right surface parameters that do not cause excitation of the fast subsystem.

Conclusion
Two sliding mode controller designs for singularly perturbed systems have been proposed. The designs are based on a block diagonal transformation of the system into fast and slow subsystems. The first design method proposes using two separate sliding mode controllers, one for the slow subsystem and a second for the fast subsystem. Simulation results indicate improved performance in comparison to previously published design methods since the errors which are produced during the decoupled process have been avoided. In the second proposed design method a single sliding mode controller is designed only for the slow subsystem and the fast subsystem is considered as high frequency unmodeled dynamics. These allow us to avoid measuring fast state, which is usually difficult to measure. Simulation results indicate good performance where the proposed controllers have less control effort compare to dual controller used before.