About this Journal Submit a Manuscript Table of Contents
Abstract and Applied Analysis
Volume 2013 (2013), Article ID 141058, 16 pages
Research Article

Variable Structure Disturbance Rejection Control for Nonlinear Uncertain Systems with State and Control Delays via Optimal Sliding Mode Surface Approach

1School of Mathematics and Computer Science, Yunnan Nationalities University, Kunming 650500, China
2Key Laboratory in Software Engineering of Yunnan Province, Kunming 650091, China

Received 18 July 2013; Accepted 9 September 2013

Academic Editor: Valery Y. Glizer

Copyright © 2013 Jing Lei et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The paper considers the problem of variable structure control for nonlinear systems with uncertainty and time delays under persistent disturbance by using the optimal sliding mode surface approach. Through functional transformation, the original time-delay system is transformed into a delay-free one. The approximating sequence method is applied to solve the nonlinear optimal sliding mode surface problem which is reduced to a linear two-point boundary value problem of approximating sequences. The optimal sliding mode surface is obtained from the convergent solutions by solving a Riccati equation, a Sylvester equation, and the state and adjoint vector differential equations of approximating sequences. Then, the variable structure disturbance rejection control is presented by adopting an exponential trending law, where the state and control memory terms are designed to compensate the state and control delays, a feedforward control term is designed to reject the disturbance, and an adjoint compensator is designed to compensate the effects generated by the nonlinearity and the uncertainty. Furthermore, an observer is constructed to make the feedforward term physically realizable, and thus the dynamical observer-based dynamical variable structure disturbance rejection control law is produced. Finally, simulations are demonstrated to verify the effectiveness of the presented controller and the simplicity of the proposed approach.

1. Introduction

Various approaches have been proposed to solve disturbance rejection problems, such as control [1], adaptive control [2], internal model control [3], variable structure control (VSC) [4], and optimal control [57]. Optimal disturbance rejection appears in a variety of applications, vehicle engine [5], damped systems [6], and spacecraft attitude control [7], for instance. However, in reality, the factor of nonlinearity, uncertainty, or time delay exhibits much influence on practical system [819], especially on today’s large information communication system. The factor of nonlinearity, uncertainty, or time delay leads the system and computation to be more complex and difficult to calculate. In addition, it brings discrepancy between the real parameter and the precise one. Therefore, it is necessary to take account of nonlinearity, uncertainty, or time delays in system modeling. In recent years, these issues have been given many attention; for example, see [1, 830]. Moreover, there are series techniques appeared to design controllers for such systems. Seeing that it is not an easy work to solve one of these problems, say nothing of concentrating on all of these issues. This study explores a disturbance rejection control design for a nonlinear uncertain time-delay system using the variable structure control approach, for the reason of VSC’s insensitivity to a wide class of uncertainties or disturbances [4, 8, 12, 13, 16, 3034]. In previous VSC investigations, some reports gave the solutions for uncertain or nonlinear systems, for example, [8, 30, 34]; some gave for the time delay systems [12, 13]; but few appeared to focus on systems of both nonlinear uncertain and time-delay. This paper is different from these bodies of literature. It presents a relative simple VSC solution to a nonlinear uncertain and time-delay system.

It is based on the ideas of the approximating sequence method [2527] and the finite spectrum assignment approach [2022]. The contributions of this study are as follows: firstly, the functional transformation method is applied, which transforms the time-delay system to a delay-free one and reduces the original problem from an infinite-dimension space to a finite-dimension one; subsequently, the optimal sliding mode surface (OSMS) is designed by using approximating sequence method, which simplifies the nonlinear OSMS design problem to a linear two-point boundary value (TPBV) one; furthermore, the corresponding compensators are designed in the variable structure disturbance rejection control (VSDC) so that the nonlinearity, uncertainty, and the time-delay effects are entirely compensated consequentially; moreover, the disturbance effect is reduced by the feedforward compensation term; additionally, to realize its physical implementation, a reduced-order observer is constructed to reconstruct the disturbance state vectors, and thus the observer-based dynamical VSDC is produced; finally, the designed VSDC is employed to a quarter-car suspension model which possesses the nonlinear, uncertain, and time-delay properties, and by comparing the system responses with the open-loop system (OLS), the effectiveness and the simplicity of the designed control are demonstrated.

This paper is organized as follows. After an introduction, in Section 2, system description is done. The OSMS and the corresponding VSDC in finite-time horizon are designed in Section 3. In Section 4, the infinite-time horizon ones are given. In Section 5, numerical simulations are illustrated. Concluding remarks are given in Section 6.

2. System Descriptions

Consider the uncertain nonlinear system with state and control delays: where , are state and control vector, respectively; , , and are constant matrices; are constant control and state delays, respectively; and is the initial state vector. represents a bounded matching uncertainty produced by the system perturbation or internal disturbance. is the nonlinearity with satisfying Lipschitz condition. In this study, control is assumed unlimited. Meanwhile, the following assumption is needed for the derivation.

Assumption 1. It is assumed that system (1) is spectral controllable.
In system (1), is the external disturbance input whose dynamical characteristic is known and can be described by the following exosystem: where is the disturbance state vector and  and are constant matrices. The pair is observable completely. The initial condition can be unknown. Exosystem (2) can describe most general disturbances, such as step signal with unknown amplitude, sinusoidal signal with known frequency and unknown amplitude and phase [28], or random signal [29]. Due to the dynamical characteristics of the external disturbances, exosystem (2) may be Lyapunov stable or asymptotically stable; that is, or , respectively, where denote the eigenvalues of a matrix.

In what follows, we will reduce the controller design problem of original system (1) to that of a delay-free one referring to the method proposed in [20] which applied this method to a linear time-delay system. Now, we develop it to the nonlinear time-delay system (1).

Define the following functional transformation: where is absolutely continuous and therefore differentiable almost everywhere and is an matrix to be defined. Employing differentiation under the integral, (3) in conjunction with (1) which gives From (4), it is observed that the range of (3) defines an ordinary differential system given byIf the following definitions are adopted: in which is referred to as the characteristic matrix equation and the method of its solution can be found in the researches of Fiagbedzi and Pearson [20, 21] and Zheng, Cheng, and Gao [22], and so forth, system (1) is assumed spectral controllable, which implies that system (5a) and (5b) is completely controllable [20]. Actually, here, ; but we choose the denotation instead of in order to simplify the derivation in the rest of the sections. Hence, time-delay system (1) is transformed into the delay-free system (5a) and (5b).

Notice that system (5a) and (5b) is coupled with and . Indeed, the term is nonlinear. From the previous investigations, it can be seen that by the functional transformation method [2022], is admissible to system (5a) and (5b) if and only if is admissible to system (1). Moreover, the controller stabilizing (1) also stabilizes (5a) and (5b). Based on this idea, in what follows, we will develop this approach to solve the VSDC problem of system (1) through the equivalent delay-free system (5a) and (5b).

First, system (5a) and (5b) will be converted into a regular form. Assume that the uncertainty is matched; that is, there exists , where is of full column rank [34], so that there exists a nonsingular matrix such that where and is nonsingular. Then, denote the following nonsingular transformations: where . Via above transformations, system (5a) and (5b) can be converted into the regular form of

where is the nonlinear and Lipschitz functions with . The pair is controllable and guarantees controllable. And is the uncertain function, where there exists a vector function , subject to in which are scalar functions and the components of are less than or equivalent to the relevant components of , that is, This is supposed a known to be condition to the controller designers [34].

Thereby, we will demonstrate the control law design process of finite-time horizon and infinite-time horizon in Sections 3 and 4, respectively.

3. Design of VSDC in Finite-Time Horizon

This section is to outline the sufficient and necessary condition for the optimality of an optimal sliding mode surface subject to the nonlinear dynamical constraint (9a) and (9b). In what follows, we will give the design procedure in two steps: OSMS design and VSDC design.

3.1. Design of OSMS in Finite-Time Horizon

To design a sliding mode surface for system (9a) and (9b), treat as a virtual control of system (9a). And respecting the finite-time performance index where and are positive semidefinite matrices, is a positive definite matrix, and is assumed observable; the OSMS will be obtained by solving the optimal control problem of system (9a) with performance index (12).

3.1.1. Global OSMS of Nonlinear Systems

Consider the optimal control problem (9a) subject to (12), treating as an excitation term despite its relationship with . Thus, the main result of Theorem 2 is achieved.

Theorem 2. Consider the optimal sliding mode surface design problem described by system (1) under disturbance (2) with quadratic performance index (12). The optimal sliding mode surface is existent and unique which is given by where is the unique positive definite solution of Riccati matrix differential equation: is the unique solution of Sylvester matrix differential equation: and is the unique solution of the adjoint equation: The optimal state is the solution of the closed-loop system:

Proof. In analogy to classical linear quadratic regulator optimal control theory from minimum principle, the Hamiltonian for the linear-quadratic control problem (9a) with respect to (10) and (11) becomes which satisfies the canonical equations: with the transversality condition: and the control equation: giving along an optimal trajectory. Then, the canonical equations (19) and (20) result in the coupled nonlinear TPBV problem: The virtual optimal control law is (23). The optimal sliding mode surface is defined by For . The objective is to design forcing system trajectories reach sliding mode surface (25).
From (21), it shows that the costate has linear relationship with state . Thus, denote the costate vector where unknown matrices , and continuous function are to be determined. Differentiating (26) and substituting the first equation of (24) and (26) into the result follow the following equation: Moreover, putting (26) into the second equation of (24) yields the following equation: Equations (27) and (28) are equivalent directly and give the Riccati differential equation (14), Sylvester differential equation (15), and the adjoint differential equation (16). Then, substituting (26) into the optimal sliding mode surface (25) yields From the denotation in (8), the OSMS (29) becomes Replacing of (30) by that of (3) yields the optimal sliding mode surface (13).
From (23) and (26), it follows the optimal virtual control Substituting (31) into (9a) with (5b) and (8) obtains the closed-loop system (17).
Matrix differential equations (14)–(16) satisfy the conditions of existence and uniqueness which implies that the OSMS (13) is existent and unique. The proof is completed.

3.1.2. Approximations of Sequences of State and Costate Equations

However, noting that (16) and (17) are coupled nonlinear differential equations, they are complex and seldom have analysis solutions. So, the sequence approximation method is adopted to obtain the solutions. Differential equations (16) and (17) can be replaced by sequences of linear time-varying (LTV) approximations [27]: Notice that approximation sequences (32) and (33) are represented as inhomogeneous linear differential equations, which have the following solutions given by the variation of constant formula: in which denotes the transition matrix generated by , that is,

In what follows, with the purpose of proving sequences and , converge to the solution of (16) and (17), and some preliminaries should be carried out.

Firstly, we should arrange the equations in compact forms. Noting that (32) and (33) are equivalent to (34) and (35), respectively, one formula will be chosen for proving. We select the latter giving the proof. In order to simplify the derivation, denote (16) and (17) in a compact form: and combine (32) and (33) in another compact form: where denotes the state transition matrix corresponding to matrix and the nonlinear function with satisfies the uniformly boundedness condition and the Lipschitz conditions, respectively; that is, ,   for all  , ,  for all  , where is some finite constant and .

Secondly, we will prove approximating sequences (38) convergent to the solution of (37). As a result, Lemma 3 is available.

Lemma 3. Let satisfy (A1) and (A2). Then, the limit of the solution of the approximating sequences (38) on converges to the unique solution of (37) on .

Proof. Let Note that . Then, . From (38), it follows that With (A2), (5b), and (8), it gives where . The above inequality gives Subsequently, By analogy, it has According to trigonometry inequality, for any , the following holds: that is, there exists Via the same way, there exists . Hence, . From (5b) and (8), it follows Therefore, and are sequences of Cauchy in Banach space , respectively. Therefore, on and on . The proof of Lemma 4 is completed.

Thirdly, result of (38) and (37) described by Lemma 3 directly carries out the same result of (32)–(35) and (16)-(17) described by Lemma 4.

Lemma 4. Let be bounded and Lipschitz continuous in its arguments. Then, the limits of the solutions of the approximating sequences (32) and (33) (or (34) and (35)) on globally converge to the unique solutions of (16) and (17), respectively, on .

Last but not least, the main result of OSMS by approximation sequences is got as Theorem 5.

Theorem 5. Given the nonlinear system (1) and the cost functional (12), then the optimal sliding mode surface is given by the limit of the sequence where is the unique positive definite solution of Riccati matrix differential equation (14), is the unique solution of the Sylvester matrix differential equation (15), and is given by the converged unique solution of the LTV differential equation sequence (32) or (34). The optimal state is the solution of the closed-loop system (33) or (35).

Remark 6. Actually, in practice, the limit of can be obtained by replacing with a positive integer . Relevantly, the following suboptimal sliding mode surface is achieved: corresponding to the following performance index: where is determined by a small enough error criterion when the following inequality holds:
Consequently, we give the algorithm of OSMS.

Algorithm 7. OSMS of system (1).1. Solve and from (14) and (15). Give some positive real constant . Set and . 2. Obtain the th adjoint vector from (34). 3. Let . Calculate from (35) and from (49).4. Calculate from (50).5. When condition (51) is satisfied, stop and output .6. Calculate from (35).

3.2. Design of VSDC in Finite-Time Horizon

In this section, we will present the VSDC design of finite-time horizon based on the proposed OSMS in the last section.

3.2.1. VSDC Design

To begin with, a trending law should be chosen. As known to all, different variable structure control laws can be designed corresponding to different sliding mode reachable conditions. In this paper, we will adopt the exponential trending law introduced by Gao [30, 34] as follows: where and are parameters to be selected which determine system instant response and convergence speed of the sliding mode reaching, respectively. The related reachable condition should be satisfied: for , that is, . For the reason of the reachable condition (54), variable structure control satisfying condition (54) will drive the state trajectories to reach the sliding mode surface in exponential velocity and then remain on it [30, 34]. Consequentially, in the basis of trending law (52), the VSDC is presented as Theorem 8.

Theorem 8. Consider system (1) under disturbance (2) with quadratic performance index (12). The sliding mode surface is given by (13). Then, the state trajectories can reach the sliding mode surface in finite time and remain on it by the variable structure control law:

Proof. Differentiating sliding mode surface (29) with respect to time and substituting system (9a) and (9b) into it give Substituting (26) into (25) and putting the result into (52) yield Comparing (56) and (57), the variable structure control law is designed: Replacing in (58) by of (8) and (3) yields the VSDC (55).
Besides, the reachable condition (54) is verified in what follows. Substituting the VSDC (58) into (56) gives If , the following holds: for all  . Since and . Thus, (59) becomes Alternatively, if , employing the same results, Therefore, combining (62) and (63), the VSDC (55) which satisfies reachable condition , can make the state trajectories reach the sliding mode surface (13) in finite time and remain on it [30, 34]. The proof is completed.

Remark 9. Actually, for the reason in Remark 6, the suboptimal sliding mode surface (SOSMS) (49) can be applied in practice. Thus, the relevant VSDC for SOSMS (49) can be gotten:

Remark 10. Summarily, on the ideal sliding mode surface, there exists . And then, putting virtual control (31) into system (9a) yields the ideal sliding mode equations of system (1): At this time, taking into VSDC (55) leads to the equivalent control of system (1):

3.2.2. Closed-Loop Stability Analysis

Select the OSMS function as a Lyapunov function candidate, that is, . Then, the increment of it is . As known, the OSMS satisfies the reachable condition (54), that is, . Thus, is negative definite, which implies . Moreover, from (3) and (8), it can be proved that as indicating the closed-loop system states asymptotically stable.

3.2.3. Physical Realization of VSDC

The VSDC in (55) contains disturbance state which is physically unmeasurable. To solve this problem, one can construct a reduced-order observer for disturbance vector so as to reconstruct the disturbance variables [28, 29, 35].

Since , there exists an arbitrary matrix such that is nonsingular. Note that , where and . It is obvious that where denotes the -order identity matrix. Consequently, the following result can be obtained.

Theorem 11. Consider system (1) under disturbance (2) with quadratic performance index (12). Respecting the sliding mode surface (13), the dynamical sliding mode control law is givenby which the state trajectories can reach the sliding mode surface in finite time and remain on it.

Proof. Defining the nonsingular transformation it gets Introduce a new variable where is the gain matrix to be selected. Then, differentiating in (71) with respect to time and substituting (70) into it yield the observation state equation: where is the component state to be reconstructed. Note that the other component state can be measured by . On the other hand, from (69), it follows that Thus, substituting the second equation in (72) into (73) yields Relevant to (74), the reduced-order observer can be written as where , represent the observed values of , , respectively. Denote the related errors as , . Subtracting (74) from (75) results in the error equation
Now, the problem is converted into finding the observer gain such that error system (76) is asymptotically stable. That the pair is observable ensures the pair observable; thus there exists the gain such that all eigenvalues of matrix can be assigned to the desired position in left-half plane, which guarantees error system (76) exponentially stable. It implies , that is, . Hence, rewriting VSDC (57) by replacing the disturbance state by the estimated state of (75):yields the dynamical VSDC (68). The proof is completed.

At this time, the VSDC design in finite-time horizon is completed. In the next section, we will discuss the corresponding results of infinite-time horizon.

4. VSDC Design in Infinite-Time Horizon

To design an OSMC in the infinite-time horizon, regarding the stability situations of exosystem (2), two different performance indexes in infinite-time horizon can be chosen. When the disturbance is asymptotically stable, that is, , the following general one should be selected: where is a positive semidefinite matrix, is a positive definite matrix, and is observable. On the other hand, when the disturbance is Lyapunov stable, that is, , since in (78) concluding disturbance may make the general infinite-horizon performance index (78) not convergent, in this case, the following average one can be selected:

Then, our object is to design OSMS and VSDC with respect to performance index (78) and (79) subject to dynamic constraint (9a) and (9b).

In reality, when , the time-variable matrices and of finite-time horizon approach to some constant matrices denoted by and , respectively, that is, and [35]. Therefore, the results of the infinite-time horizon can be directly obtained from those of the finite-time horizon. Yet, for the sake of briefness, the corresponding results are omitted here.

5. Simulation Example

In this section, a 2DOF quarter-car model is applied to simulate the system responses under the designed VSDC comparing with those of OLS. The governing dynamic equations are where is the sprung mass representing car chassis; is the unsprung mass representing wheel assembly; is the linear stiffness coefficients; is the damping of the uncontrolled suspension; stands for the stiffness of tire respectively; are displacements of sprung and unsprung masses, respectively; is the road displacement input. The actuator force acts between sprung and unsprung masses. Let represent the spring nonlinearity and , and let   indicate the uncertainty parameters of spring mass, stiffness of spring, and tire with known bounds. Summarizing equations (80), one can get the nonlinearity of the system, which was employed in some references, for example, [14, 15] and the uncertainties of the system which are referred to [1618]: More characteristic details about these kinds of nonlinearity and uncertainties can be found in abovementioned [1418]. By defining the set of state variable where is the suspension deflection, the tire deflection, the sprung mass velocity, and the unsprung mass velocity, the state vector is denoted by . Then, the vehicle suspension system is rewritten in state-space representation (1) with

Adopting the parameter values listed in Table 1 [19], the associated matrices and vectors of (84) are as follows:

Table 1: Parameters and values of quarter-car suspension.

Set the initial value and time-delays , . Take , for ;  , where indicates the percentage of variation allowed around its nominal value and determine the actual parameter derivation changing in interval [18]. Referring to [16, 17], given the possible variations of uncertain parameters: , , , and , note they represent 28.5% uncertainty in , 15% in and , and 25% in , which yields the uncertainty bound parameters:  kg,  N/m,  N/m,  N/m. Consequentially, the bound vector function is produced:

Then, transform the system in the form of (9a) and (9b). To simulate the road profiles, two cases of sinusoidal signal and attenuated signal are employed.

Case 1. Consider sinusoidal disturbance described by exosystem (2) with which produces the eigenvalues of are ,  . Take the performance index (79) with . First of all, with (14) and (15) the matrices , are computed. Moreover, from (32) and (33), and are solved. Consequentially, setting and in trending law (52), the OSMS (49) is gotten as well as VSDC (77) is obtained. Finally, the suspension responses by VSDC and OLS are demonstrated in Figures 1 and 2.

Figure 1: States trajectories of VSDC and OLS under sinusoidal disturbance.
Figure 2: Control inputs and sliding mode surface under sinusoidal disturbance.

From Figure 1, it can be seen that the states are stabilized and become stable by using VSDC and the magnitudes of the suspension response are obviously reduced comparing with those of the OLS. The reason can be concluded from Figure 2, that by using VSDC, the states were controlled towards the OSMS; while being forced onto the OSMS, the equivalent control drove them stabilized and hold good dynamical characteristics.

Case 2. Consider attenuated disturbance described by exosystem (2) with which yields that the eigenvalues of are ,  . Take performance index (12) with . Through the same way, from (14) and (15), it yields the values of matrix , . And from (32) and (33), and are solvable. Then, the OSMS (49) and VSDC (77) are obtained taking and in trending law (52). Hence, as displayed in Figures 3 and 4, the suspension responses by VSDC and OLS are illustrated.

Figure 3: States trajectories of VSDC and OLS under attenuated disturbance.
Figure 4: Control inputs and sliding mode surface under attenuated disturbance.

It is shown that the state variables are reduced by VSDC comparing with OLS, which are controlled asymptotically stable at about 30 sec. Clearly, the VSDC rejected the disturbance which is faced by the system and makes the state variable of OLS vibrate.

Summarily, comparing the simulation results between Figures 1 and 2 and Figures 3 and 4, we can see that, although the states are reduced by VSDC and lower than those of OLS when facing sinusoidal disturbance, it still vibrates a little; whereas when facing attenuated disturbance, they converge to zero gradually. The reason lies in the main different feature of these two kinds of signals. Notice that in VSDC (55), there exists the feedforward compensation term described as