Abstract

For a class of linear MIMO uncertain systems, a dynamic sliding mode control algorithm that avoids the chattering problem is proposed in this paper. Without using any differentiator, we develop a modified asymptotically stable second-order sliding mode control law in which the proposed controller can guarantee the finite time convergence to the sliding mode and can show that the system states asymptotically approach to zero. Finally, a numerical example is explained for demonstrating the applicability of the proposed scheme.

1. Introduction

Sliding mode control (SMC) has been successfully used in controlling many uncertain systems [1, 2]. For a system with the matched disturbance, SMC can obtain the perfect disturbance rejection during the sliding mode. The control objectives are attained by constraining the system dynamics on a properly chosen sliding variable by means of discontinuous control laws. In theory, SMC offers robust stability to systems through high-gain control with the infinite fast switching action. However, high-gain control designs suffer from the drawback of peaking phenomenon, in which the control input peaks to an extremely large value during the transient stage. The peaking phenomenon can easily violate the control saturation constraint. The SMC scheme is often discontinuous and the feedback control input needs to switch with infinite switching frequency. The discontinuous high speed switching action results in the chattering problem due to the inherent delay and other problems [1, 2]. The chattering action may excite the unmodeled high order dynamics, which probably leads to unforeseen instability [15]. There are two major approaches reported to cope with the chattering problem. The first approach is to insert a fixed or variable boundary layer near the sliding variable [6], so that a continuous control replaces the discontinuous one when the system is inside the boundary layer. Another approach to eliminate the chattering is carried out by adding an auxiliary control input into the system such as using of fuzzy control [7, 8] and adaptive fuzzy control [9]. These mentioned methods [15] can give a chattering-free system, but a finite steady-state error may exist. Hence, the property of perfect disturbance rejection cannot be guaranteed.

Dynamic sliding mode control [1014] where an integrator is used in front of the system is a special approach to eliminate the control chattering. The time derivative of the control input is treated as the new control variable for the augmented system in which the augmented system is includes the original system and the integrator. Since no boundary layer is used in the controller, the advantages of a dynamic sliding mode controller are that chattering reduction is obtained by using an integrator and the property of perfect disturbance rejection is guaranteed. Hence, dynamic sliding mode control not only removes some of the fundamental limitations of the traditional approach but also provides improved tracking accuracy under sliding mode. The main problem in implementation of dynamic sliding mode controller is the increasing information demand in which the knowledge of the derivative of the sliding variable is required. Bartolini et al. [11, 12] presented a suboptimal version of the twisting algorithm [15, 16] to cope with the chattering problem. However, this method requires at least the knowledge of the sign of the derivative of the sliding variable. Chen et al. [14] applied an LTR observer to overcome the problem of sliding variable estimation. Recently, Levant [1517] presented the twisting algorithm to stabilize second-order nonlinear systems but required the knowledge of the derivative of the sliding variable. The super twisting algorithm [15, 16] does not require the output derivative to be measured but it has been originally developed and analyzed for system with relative degree one. Levant [17] proposed an exact finite time convergence differentiator which can successfully estimate the derivative of the sliding variable but required a priori bounded trajectory.

For a linear MIMO system with the matched disturbance, an alternative dynamic sliding mode controller for avoiding the chattering problem is established in this paper. We first present a modified second-order sliding mode control method to stabilize the perturbed system. Introducing a proportional-integral term of the sliding variable and an integral sign function term into the control law, the resulting control forces are chattering-free. Based on the developed second-order sliding mode technique, it is shown that the finite time convergence to the sliding mode is provided theoretically and the zero steady-state error can be guaranteed by applying the proposed control law. As a result, the control accuracy is better than those performed by the conventional boundary layer control [35]. Moreover, a sufficient condition for the closed-loop stability is given and the implementation of the proposed control algorithm is simple. Comparing with the conventional dynamic sliding mode controllers [1014], the proposed method does not need any observer structure to estimate the derivative of the sliding variable. Finally, the feasibility of the proposed method is illustrated by a numerical example.

The work of this paper is organized as follows. Section 2 describes a class of uncertain MIMO linear systems and gives the problem formulation. Section 3 presents the modified second-order sliding mode method and develops the dynamic sliding mode controller design. The simulation result is included in Section 4. Section 5 offers a brief conclusion.

2. Problem Formulation

Consider an uncertain system satisfying the matched condition of the form ??(??)=????(??)+??(??(??)+??(??,??)),(1) where ???R?? is the state vector, ???R?? is the control input vector, ???R?? is the output vector, and ???R?? is the unknown matched disturbance vector with the known upper bounds ???(??,??)?=??1 and ???(??,??)?=??2. Suppose that the system states are accessible for measurement and the pair (??,??) is stabilizable. Let the sliding variable be chosen as ??(??)=????(??),(2) where the matrix ???R??×?? is designed to stabilize the reduced-order system. To satisfy the reaching and sliding condition, the control input for the conventional sliding mode controller is designed as ??(??)=-(????)-1(??????(??)+??sign(??(??))),(3) where ??>0 is a high gain to design such that the system reaches and slides on the sliding variable in finite time. However, a phenomenon called chattering is generated due to the discontinuous function sign(??(??)). It can be considered as the undesired chattering effect produced by the high switching action of the control input. As a result, the chattering becomes the main implementation problem of SMC. Numerous techniques have been proposed to eliminate this phenomenon in SMC [314]. The one of the most common solution to reduce the chattering is the boundary layer techniques [35]. However, the boundary layer thickness has the trade-off relation between the control performance and the chattering migration. Another drawback with applying the boundary layer methods [35] is the reduction of the control accuracy.

For a linear MIMO system with the matched disturbance, in this paper we propose a dynamic sliding mode control algorithm in which the proposed procedure can effectively reduce the chattering effect. We first develop a modified second-order sliding mode control method to stabilize the perturbed systems. Introducing a proportional-integral term of the sliding variable and an integral sign function term into the controller, the finite time convergence to the sliding mode is guaranteed by applying the developed second-order sliding mode technique. Moreover, all the states of the system asymptotically approach to zero once the system is in the sliding mode.

3. Dynamic Sliding Mode Controller Design

In this section, we propose a dynamic sliding mode control algorithm which can successfully avoid the chattering in the linear MIMO system with the matched disturbance. A modified second-order sliding mode control algorithm that does not require the derivative of the sliding variable is presented. The proposed control law can guarantee the finite time convergence to the sliding mode and stabilize the reduced-order system in which the system states asymptotically approach to zero.

For the sliding variable (2), we design the dynamic sliding mode controller as ??(??)=-(????)-1???????(??)+??1??(??)+??2???0???(??)????+????0?,sign(??(??))????(4) where ??1?R??×??, ??2?R??×??, and ???R??×?? are the positive definite diagonal matrix given by ??1???=diag11,,??1???,??2???=diag21,,??2???,?????=diag1,,?????.(5) Taking the time derivative of (2) and substituting the control input (4) into it can obtain ??(??)=????(??)=??????(??)+????(??(??)+??(??,??))=-??1??(??)-??2???0???(??)????-????0sign(??(??))????+??????(??,??).(6) Further differentiating (6) with time yields the following dynamics: ??(??)+??1??(??)+??2??(??)=-??sign(??(??))+??(??),(7) where ??(??)=??????(??,??)=[??1?????]???R??. According to the linear algebraic theory, we have ???(??)?=???????(??,??)?=??2?????? and |????(??)|=???(??)? for 1=??=??.

Lemma 1. Consider the unperturbed system as ??(??)+??1??(??)+??2??(??)=-??sign(??(??)),(8) where the gains ??1, ??2, and ?? are positive constants. If the parameters ??1>0 and ??2>0 are chosen to satisfy the condition ??2<??21/4, then ??(??) and ??(??) asymptotically converge to zero in finite time for a sufficiently large value of ??.

Proof. Since the condition ??2<??21/4 holds, the roots of the characteristic equation, ??2+??1??+??2=0, are real, distinct, and negative. Assume now for simplicity that the initial conditions are ??(??0)=0 and ??(??0)>0. Thus the trajectory enters the half-plane ??(??)>0 (quadrant I), as shown in Figure 1. When ??(??)>0, we from (8) obtain ??(??)+??1??(??)+??2??(??)=-?? and know its equivalent point as (??,??)=(-??/??2,0). Let the function ??(??) be generated by ????(??)=??(??)+????(??)+??sign(??(??)),(9) where ??+??=??1 and ????=??2. The parameters ??>0 and ??>0 are real constants. Since ??(??)>0, it follows from ??(??)+??1??(??)+??2??(??)+??=0 that ????(??)??????+????(??)=?????????(??)+????(??)+??????+????(??)+????(??)+???=??(??)+(??+??)??(??)+??????(??)+??=??(??)+??1??(??)+??2??(??)+??=0.(10) Hence, we can obtain ??(??)=??(0)??-????. It follows that the upper bound of ??(??) is |??(??)|=????-???? where ??>0 is a constant. Choose a Lyapunov function as ??(??)=(1/2)??2(??) and then obtain its time derivative as ?????(??)=??(??)??(??)=??(??)??(??)-???||||??sign(??(??))-????(??)=????(??)-????-????||||??(??)-????2(??).(11) Since ??>0, there exists a finite time, ??1>0, such that ????-????<(??/??)-?? for a sufficiently large gain ?? and ??>??1, where ??>0 is a constant. Hence, ||||??(??)=-????(??)-????2||||(??)=-????(??)for??>??1.(12) The above equation implies that the function ??(??) converges to zero in finite time. Let the trajectory of (8) intersect next time with the axis ??(??)=0 at the point ??(??1). Since the roots of the characteristic equation, ??2+??1??+??2=0, are stable, we know that the spiral trajectories converge to the equivalent point and the behavior of ??(??) changes monotonously. Hence, ||?????1?||||?????0?||=??<1.(13) Extending the trajectory into the half plane ??(??)<0 after a similar reasoning achieves that the successive crossing the axis ??(??)=0 satisfies the inequality |??(????+1)|/|??(????)|=??<1. Therefore, its solutions cross the axis ??(??)=0 from quadrant II to quadrant I, and from quadrant IV to quadrant III. Every trajectory of the system crosses the axis ??(??)=0 in finite time. After gluing these paths along the line ??(??)=0, we obtain the phase portrait of the system, as shown in Figure 1. This algorithm features a twisting of the phase portrait around the origin and an infinite number encircling the origin occurs. According to Levant’s papers [1416], the total convergence time is estimated as ?||?????=?????||=||?????0?||?1+??+??2?=||???+???0?||.1-??(14) As a result, we can show that the trajectories perform rotations around the origin while converging in finite time to the origin of the phase plane. The finite time convergence to the origin is due to switching between two different control amplitudes as the trajectory comes nearer to the origin [18]. The proof of the theorem is finished.

Lemma 2. Consider the following system: ??(??)+??1??(??)+??2??(??)=-??sign(??(??))+??(??),(15) where the function ??(??) has the upper bound |??(??)|=?? and ??>0 is a known constant. If the parameters ??1>0 and ??2>0, and the gain ?? are chosen to satisfy the following conditions: ??2<??214,??>??,(16) then ??(??) and ??(??) converge to zero in finite time.

Proof. We take the parameters ??1 and ??2 to satisfy the condition ??2<??21/4 and then obtain the characteristic polynomial of system (15) having two distinct real roots ??1,2=-??,-??, where ??>0, ??>0, ??>??, ??1=??+??, and ??2=????. When the condition ??(??)>0 holds, (15) becomes ??(??)+??1??(??)+??2??(??)=-??+??(??).(17) Let ??1=??+(??/??2) and ??2=??1=??. It follows that ??????(??)=01-??2-??1??????????1??(??)2(????+????0??????)??(??)=F??(??)+????(??),(18) where ?????=??1????2????R2????,F=01-??2-??1????????01????,??=.(19) Write the above dynamic equation as its solution in an explicit form ??(??)=??F????(0)+??F?????0??F??????(??)????=??F?????(0)+??0??F??????(??-??)????,(20) where ??F??=?????1???-??????-????-????-?????1?????-??-????-??-??????????????-??-????-??-?????1???-??????-????-????-??????????.(21) The upper bounds of ??1(??) and ??2(??) can be constructed as ||??1(||??)=??1??-????+1???-????0||??-????-??-????||||||??(??-??)????=??1??-????+?????-????0||??-????-??-????||????=??1??-????+?????-??80???-????-??-?????????=??1??-????+?????-????-???????=??1??-????+????2,||??2||(??)=??2??-????+??????-????0||????-????-????-????||?????=??2??-????+??×????-??80||????-????-????-????+(??-??)??-????||?????=??2??-????+??×????-??80?????-????-????-??????????+(??-??)80??-?????????=??2??-????+???????-????-??+??????-?????=??2??-????+2??????2,(22) where ??1>0 and ??2>0 are constants. It follows that ||??1||(??)=??1??-????+????2,||??2||(??)=??2??-????+2??????2.(23) Since ??1=??+(??/??2), we can from (23) obtain ||||????(??)+??2||||=??1??-????+????2.(24) This part shows that the ball of radius ??=??/??2 with center located at (??/??2,0) is an attractor ????1. Similar to the work, we have, when ??(??)<0, the ball of radius r, with center located at (-??/??2,0) as another attractor ????2. Choose the gain ??to satisfy the inequality ??>?? and then we have ?????2?=??-????2-????2?-??>0,??2?=+??-????2+????2<0.(25) It follows from the above two inequalities that the two attractors ????1 and ????2 do not intersect each other, and the behavior of the perturbed system (15) will be qualitatively similar to the behavior of the nominal system. Therefore, the perturbed system converges to the origin in the same way of the nominal system and the condition that ??(??) and ??(??) converge to zero in finite time can be guaranteed. We finish the proof of the lemma.

Theorem 3. Consider system (1) with the sliding variable (2) and the control input (4). Let ??=??2??????. If the elements of these matrices ??1, ??2, and ?? satisfy the following conditions: ??2??<??21??4,????>??,for??=1,,??,(26) then one can show that the system states asymptotically approach to zero.

Proof. We first express system (7) as a set of second-order systems with the form ????(??)+??1????(??)+??2????(??)=-???????sign???(??)+????(??)for??=1,,??,(27) where |????(??)|=?? for ??=1,,?? where ??=??2??????>0 is a known constant. Applying the result of Lemma 2 into (27), if the parameters ??1??, ????2, and ???? can satisfy ??2??<??21??4,????>??,for??=1,,??,(28) then ????(??) and ????(??) asymptotically converge to zero in finite time according to Lemma 2. When the condition ??(??)=?? is guaranteed, it follows from the concept of equivalent control [1, 2] that the system dynamics in the sliding mode is ???(??)=??-??(????)-1???????(??).(29) Since the nonzero eigenvalues of the matrix ??-??(????)-1???? are placed in the left-half plane, we can conclude that the system states asymptotically approach to zero and finally finish the proof of the theorem.

4. Numerical Example

To demonstrate the effectiveness of the proposed method, we consider an unstable batch reactor where the matched disturbance is introduced into the system and its state space form can be given by ????????????????+???????????????????????????(??)=1.38-0.20776.715-5.676-0.5814-4.2900.6751.0674.273-6.6545.8930.0484.2731.343-2.104??(??)005.67901.136-3.1461.1360??(??)+1.5sin(0.5??)1cos????2?????????.(30) The sliding variable is chosen as ??????????(??)=-0.01060.168400.0383-0.19310.0155-0.31790.2404??(??),(31) where the nonzero eigenvalues of the system in the sliding mode are assigned as {-2,-3}. We design the conventional sliding mode controller and the boundary layer controller, respectively, as ??????????????????(??)=--0.1108-0.5565-0.020.0934-0.6031-0.35761.1415-1.2726??(??)-3sign(??(??)),??(??)=--0.1108-0.5565-0.020.0934-0.6031-0.35761.1415-1.2726??(??)-3sat(??(??),0.05).(32) The control input with applying the proposed algorithm is given by ???????????(??)=--0.1108-0.5565-0.020.0934-0.6031-0.35761.1415-1.2726??(??)-10??(??)-2??0???(??)????-3??0sign(??(??))????.(33) Three cases are simultaneously simulated under the initial condition ??(0)=[12-21]?? and the simulation is carried out at a fixed step size of 1 milliseconds. The time responses of the system states in the three cases are shown in Figures 2, 4, and 6, respectively. As can be seen from these figures, the proposed method obtains the zero steady-state error. Figures 3, 5, and 7 show the control inputs in the three cases, respectively. As illustrated, our control force does not exhibit any chattering effects. The responses of the sliding variable with applying the proposed method are given in Figure 8. It is clear from the figure that our control law can guarantee the finite time convergence to the sliding mode. As a result, the simulation results verify the good performance of the proposed strategy.

5. Conclusions

In this paper we have proposed a modified second-order sliding mode control algorithm to avoid the chattering problem for a MIMO uncertain system. The developed approach does not require the derivative of the sliding variable, thus eliminating the requirement of designing a differentiator. Under the proposed dynamic sliding mode controller, we show that the finite time convergence to the sliding mode is guaranteed and the system states can asymptotically approach to zero. Simulation results demonstrate that the proposed control scheme exhibits reasonably good system performance.