Sliding mode control

Short description: Method in nonlinear control theory

In control systems, sliding mode control (SMC) is a nonlinear control method that alters the dynamics of a nonlinear system by applying a discontinuous control signal (or more rigorously, a set-valued control signal) that forces the system to "slide" along a cross-section of the system's normal behavior. The state-feedback control law is not a continuous function of time. Instead, it can switch from one continuous structure to another based on the current position in the state space. Hence, sliding mode control is a variable structure control method. The multiple control structures are designed so that trajectories always move toward an adjacent region with a different control structure, and so the ultimate trajectory will not exist entirely within one control structure. Instead, it will slide along the boundaries of the control structures. The motion of the system as it slides along these boundaries is called a sliding mode and the geometrical locus consisting of the boundaries is called the sliding (hyper)surface. In the context of modern control theory, any variable structure system, like a system under SMC, may be viewed as a special case of a hybrid dynamical system as the system both flows through a continuous state space but also moves through different discrete control modes.

Introduction

Figure 1: Phase plane trajectory of a system being stabilized by a sliding mode controller. After the initial reaching phase, the system states "slides" along the line [math]\displaystyle{ s=0 }[/math]. The particular [math]\displaystyle{ s=0 }[/math] surface is chosen because it has desirable reduced-order dynamics when constrained to it. In this case, the [math]\displaystyle{ s=x_1 +\dot{x}_1 = 0 }[/math] surface corresponds to the first-order LTI system [math]\displaystyle{ \dot{x}_1 = -x_1 }[/math], which has an exponentially stable origin.

Figure 1 shows an example trajectory of a system under sliding mode control. The sliding surface is described by [math]\displaystyle{ s=0 }[/math], and the sliding mode along the surface commences after the finite time when system trajectories have reached the surface. In the theoretical description of sliding modes, the system stays confined to the sliding surface and need only be viewed as sliding along the surface. However, real implementations of sliding mode control approximate this theoretical behavior with a high-frequency and generally non-deterministic switching control signal that causes the system to "chatter" in a tight neighborhood of the sliding surface. Chattering can be reduced through the use of deadbands or boundary layers around the sliding surface, or other compensatory methods. Although the system is nonlinear in general, the idealized (i.e., non-chattering) behavior of the system in Figure 1 when confined to the [math]\displaystyle{ s=0 }[/math] surface is an LTI system with an exponentially stable origin. One of the compensatory methods is the adaptive sliding mode control method proposed in which uses estimated uncertainty to construct continuous control law. In this method chattering is eliminated while preserving accuracy (for more details see references [2] and [3]). The three distinguished features of the proposed adaptive sliding mode controller are as follows: (i) The structured (or parametric) uncertainties and unstructured uncertainties (un-modeled dynamics, unknown external disturbances) are synthesized into a single type uncertainty term called lumped uncertainty. Therefore, a linearly parameterized dynamic model of the system is not required, and the simple structure and computationally efficient properties of this approach make it suitable for the real-time control applications. (ii) The adaptive sliding mode control scheme design relies on the online estimated uncertainty vector rather than relying on the worst-case scenario (i.e., bounds of uncertainties). Therefore, a-priory knowledge of the bounds of uncertainties is not required, and at each time instant, the control input compensates for the uncertainty that exists. (iii) The developed continuous control law using fundamentals of the sliding mode control theory eliminates the chattering phenomena without trade-off between performance and robustness, which is prevalent in boundary-layer approach.

Intuitively, sliding mode control uses practically infinite gain to force the trajectories of a dynamic system to slide along the restricted sliding mode subspace. Trajectories from this reduced-order sliding mode have desirable properties (e.g., the system naturally slides along it until it comes to rest at a desired equilibrium). The main strength of sliding mode control is its robustness. Because the control can be as simple as a switching between two states (e.g., "on"/"off" or "forward"/"reverse"), it need not be precise and will not be sensitive to parameter variations that enter into the control channel. Additionally, because the control law is not a continuous function, the sliding mode can be reached in finite time (i.e., better than asymptotic behavior). Under certain common conditions, optimality requires the use of bang–bang control; hence, sliding mode control describes the optimal controller for a broad set of dynamic systems.

One application of sliding mode controller is the control of electric drives operated by switching power converters. Because of the discontinuous operating mode of those converters, a discontinuous sliding mode controller is a natural implementation choice over continuous controllers that may need to be applied by means of pulse-width modulation or a similar technique of applying a continuous signal to an output that can only take discrete states. Sliding mode control has many applications in robotics. In particular, this control algorithm has been used for tracking control of unmanned surface vessels in simulated rough seas with high degree of success.

Sliding mode control must be applied with more care than other forms of nonlinear control that have more moderate control action. In particular, because actuators have delays and other imperfections, the hard sliding-mode-control action can lead to chatter, energy loss, plant damage, and excitation of unmodeled dynamics. Continuous control design methods are not as susceptible to these problems and can be made to mimic sliding-mode controllers.

Control scheme

Consider a nonlinear dynamical system described by

[math]\displaystyle{ \dot{\mathbf{x}}(t)=f(\mathbf{x},t) + B(\mathbf{x},t)\,\mathbf{u}(t) }[/math] ( 1)

where

[math]\displaystyle{ \mathbf{x}(t) \triangleq \begin{bmatrix}x_1(t)\\x_2(t)\\\vdots\\x_{n-1}(t)\\x_n(t)\end{bmatrix} \in \mathbb{R}^n }[/math]

is an n-dimensional state vector and

[math]\displaystyle{ \mathbf{u}(t) \triangleq \begin{bmatrix}u_1(t)\\u_2(t)\\\vdots\\u_{m-1}(t)\\u_m(t)\end{bmatrix} \in \mathbb{R}^m }[/math]

is an m-dimensional input vector that will be used for state feedback. The functions [math]\displaystyle{ f: \mathbb{R}^n \times \mathbb{R} \to \mathbb{R}^n }[/math] and [math]\displaystyle{ B: \mathbb{R}^n \times \mathbb{R} \to \mathbb{R}^{n \times m} }[/math] are assumed to be continuous and sufficiently smooth so that the Picard–Lindelöf theorem can be used to guarantee that solution [math]\displaystyle{ \mathbf{x}(t) }[/math] to Equation (1) exists and is unique.

A common task is to design a state-feedback control law [math]\displaystyle{ \mathbf{u}(\mathbf{x}(t)) }[/math] (i.e., a mapping from current state [math]\displaystyle{ \mathbf{x}(t) }[/math] at time t to the input [math]\displaystyle{ \mathbf{u} }[/math]) to stabilize the dynamical system in Equation (1) around the origin [math]\displaystyle{ \mathbf{x} = [0, 0, \ldots, 0]^\intercal }[/math]. That is, under the control law, whenever the system is started away from the origin, it will return to it. For example, the component [math]\displaystyle{ x_1 }[/math] of the state vector [math]\displaystyle{ \mathbf{x} }[/math] may represent the difference some output is away from a known signal (e.g., a desirable sinusoidal signal); if the control [math]\displaystyle{ \mathbf{u} }[/math] can ensure that [math]\displaystyle{ x_1 }[/math] quickly returns to [math]\displaystyle{ x_1 = 0 }[/math], then the output will track the desired sinusoid. In sliding-mode control, the designer knows that the system behaves desirably (e.g., it has a stable equilibrium) provided that it is constrained to a subspace of its configuration space. Sliding mode control forces the system trajectories into this subspace and then holds them there so that they slide along it. This reduced-order subspace is referred to as a sliding (hyper)surface, and when closed-loop feedback forces trajectories to slide along it, it is referred to as a sliding mode of the closed-loop system. Trajectories along this subspace can be likened to trajectories along eigenvectors (i.e., modes) of LTI systems; however, the sliding mode is enforced by creasing the vector field with high-gain feedback. Like a marble rolling along a crack, trajectories are confined to the sliding mode.

The sliding-mode control scheme involves

1. Selection of a hypersurface or a manifold (i.e., the sliding surface) such that the system trajectory exhibits desirable behavior when confined to this manifold.
2. Finding feedback gains so that the system trajectory intersects and stays on the manifold.

Because sliding mode control laws are not continuous, it has the ability to drive trajectories to the sliding mode in finite time (i.e., stability of the sliding surface is better than asymptotic). However, once the trajectories reach the sliding surface, the system takes on the character of the sliding mode (e.g., the origin [math]\displaystyle{ \mathbf{x}=\mathbf{0} }[/math] may only have asymptotic stability on this surface).

The sliding-mode designer picks a switching function [math]\displaystyle{ \sigma: \mathbb{R}^n \to \mathbb{R}^m }[/math] that represents a kind of "distance" that the states [math]\displaystyle{ \mathbf{x} }[/math] are away from a sliding surface.

• A state [math]\displaystyle{ \mathbf{x} }[/math] that is outside of this sliding surface has [math]\displaystyle{ \sigma(\mathbf{x}) \neq 0 }[/math].
• A state that is on this sliding surface has [math]\displaystyle{ \sigma(\mathbf{x}) = 0 }[/math].

The sliding-mode-control law switches from one state to another based on the sign of this distance. So the sliding-mode control acts like a stiff pressure always pushing in the direction of the sliding mode where [math]\displaystyle{ \sigma(\mathbf{x}) = 0 }[/math]. Desirable [math]\displaystyle{ \mathbf{x}(t) }[/math] trajectories will approach the sliding surface, and because the control law is not continuous (i.e., it switches from one state to another as trajectories move across this surface), the surface is reached in finite time. Once a trajectory reaches the surface, it will slide along it and may, for example, move toward the [math]\displaystyle{ \mathbf{x} = \mathbf{0} }[/math] origin. So the switching function is like a topographic map with a contour of constant height along which trajectories are forced to move.

The sliding (hyper)surface/manifold is typically of dimension [math]\displaystyle{ n-m }[/math] where n is the number of states in [math]\displaystyle{ \mathbf{x} }[/math] and m is the number of input signals (i.e., control signals) in [math]\displaystyle{ \mathbf{u} }[/math]. For each control index [math]\displaystyle{ 1 \leq k \leq m }[/math], there is an [math]\displaystyle{ (n-1) }[/math]-dimensional sliding surface given by

[math]\displaystyle{ \left\{ \mathbf{x} \in \mathbb{R}^n : \sigma_k(\mathbf{x}) = 0 \right\} }[/math] ( 2)

The vital part of SMC design is to choose a control law so that the sliding mode (i.e., this surface given by [math]\displaystyle{ \sigma(\mathbf{x})=\mathbf{0} }[/math]) exists and is reachable along system trajectories. The principle of sliding mode control is to forcibly constrain the system, by suitable control strategy, to stay on the sliding surface on which the system will exhibit desirable features. When the system is constrained by the sliding control to stay on the sliding surface, the system dynamics are governed by reduced-order system obtained from Equation (2).

To force the system states [math]\displaystyle{ \mathbf{x} }[/math] to satisfy [math]\displaystyle{ \sigma(\mathbf{x}) = \mathbf{0} }[/math], one must:

1. Ensure that the system is capable of reaching [math]\displaystyle{ \sigma(\mathbf{x}) = \mathbf{0} }[/math] from any initial condition
2. Having reached [math]\displaystyle{ \sigma(\mathbf{x})=\mathbf{0} }[/math], the control action is capable of maintaining the system at [math]\displaystyle{ \sigma(\mathbf{x})=\mathbf{0} }[/math]

Existence of closed-loop solutions

Note that because the control law is not continuous, it is certainly not locally Lipschitz continuous, and so existence and uniqueness of solutions to the closed-loop system is not guaranteed by the Picard–Lindelöf theorem. Thus the solutions are to be understood in the Filippov sense. Roughly speaking, the resulting closed-loop system moving along [math]\displaystyle{ \sigma(\mathbf{x}) = \mathbf{0} }[/math] is approximated by the smooth dynamics [math]\displaystyle{ \dot{\sigma}(\mathbf{x}) = \mathbf{0}; }[/math] however, this smooth behavior may not be truly realizable. Similarly, high-speed pulse-width modulation or delta-sigma modulation produces outputs that only assume two states, but the effective output swings through a continuous range of motion. These complications can be avoided by using a different nonlinear control design method that produces a continuous controller. In some cases, sliding-mode control designs can be approximated by other continuous control designs.

Theoretical foundation

The following theorems form the foundation of variable structure control.

Theorem 1: Existence of sliding mode

Consider a Lyapunov function candidate

[math]\displaystyle{ V(\sigma(\mathbf{x}))=\frac{1}{2}\sigma^\intercal(\mathbf{x})\sigma(\mathbf{x})=\frac{1}{2}\|\sigma(\mathbf{x})\|_2^2 }[/math] ( 3)

where [math]\displaystyle{ \|\mathord{\cdot}\| }[/math] is the Euclidean norm (i.e., [math]\displaystyle{ \|\sigma(\mathbf{x})\|_2 }[/math] is the distance away from the sliding manifold where [math]\displaystyle{ \sigma(\mathbf{x})=\mathbf{0} }[/math]). For the system given by Equation (1) and the sliding surface given by Equation (2), a sufficient condition for the existence of a sliding mode is that

[math]\displaystyle{ \underbrace{ \overbrace{\sigma^\intercal}^{\tfrac{\partial V}{\partial \sigma}} \overbrace{\dot{\sigma}}^{\tfrac{\operatorname{d} \sigma}{\operatorname{d} t}} }_{\tfrac{\operatorname{d}V}{\operatorname{d}t}} \lt 0 \qquad \text{(i.e., } \tfrac{\operatorname{d}V}{\operatorname{d}t} \lt 0 \text{)} }[/math]

in a neighborhood of the surface given by [math]\displaystyle{ \sigma(\mathbf{x})=0 }[/math].

Roughly speaking (i.e., for the scalar control case when [math]\displaystyle{ m=1 }[/math]), to achieve [math]\displaystyle{ \sigma^\intercal \dot{\sigma} \lt 0 }[/math], the feedback control law [math]\displaystyle{ u(\mathbf{x}) }[/math] is picked so that [math]\displaystyle{ \sigma }[/math] and [math]\displaystyle{ \dot{\sigma} }[/math] have opposite signs. That is,

• [math]\displaystyle{ u(\mathbf{x}) }[/math] makes [math]\displaystyle{ \dot{\sigma}(\mathbf{x}) }[/math] negative when [math]\displaystyle{ \sigma(\mathbf{x}) }[/math] is positive.
• [math]\displaystyle{ u(\mathbf{x}) }[/math] makes [math]\displaystyle{ \dot{\sigma}(\mathbf{x}) }[/math] positive when [math]\displaystyle{ \sigma(\mathbf{x}) }[/math] is negative.

Note that

[math]\displaystyle{ \dot{\sigma} = \frac{\partial \sigma}{\partial \mathbf{x}} \overbrace{\dot{\mathbf{x}}}^{\tfrac{\operatorname{d} \mathbf{x}}{\operatorname{d} t}} = \frac{\partial \sigma}{\partial \mathbf{x}} \overbrace{\left( f(\mathbf{x},t) + B(\mathbf{x},t) \mathbf{u} \right)}^{\dot{\mathbf{x}}} }[/math]

and so the feedback control law [math]\displaystyle{ \mathbf{u}(\mathbf{x}) }[/math] has a direct impact on [math]\displaystyle{ \dot{\sigma} }[/math].

Reachability: Attaining sliding manifold in finite time

To ensure that the sliding mode [math]\displaystyle{ \sigma(\mathbf{x})=\mathbf{0} }[/math] is attained in finite time, [math]\displaystyle{ \operatorname{d}V/{\operatorname{d}t} }[/math] must be more strongly bounded away from zero. That is, if it vanishes too quickly, the attraction to the sliding mode will only be asymptotic. To ensure that the sliding mode is entered in finite time,

[math]\displaystyle{ \frac{\operatorname{d}V}{\operatorname{d}t} \leq -\mu (\sqrt{V})^{\alpha} }[/math]

where [math]\displaystyle{ \mu \gt 0 }[/math] and [math]\displaystyle{ 0 \lt \alpha \leq 1 }[/math] are constants.

Explanation by comparison lemma

This condition ensures that for the neighborhood of the sliding mode [math]\displaystyle{ V \in [0,1] }[/math],

[math]\displaystyle{ \frac{\operatorname{d}V}{\operatorname{d}t} \leq -\mu (\sqrt{V})^{\alpha} \leq -\mu \sqrt{V}. }[/math]

So, for [math]\displaystyle{ V \in (0,1] }[/math],

[math]\displaystyle{ \frac{ 1 }{ \sqrt{V} } \frac{\operatorname{d}V}{\operatorname{d}t} \leq -\mu, }[/math]

which, by the chain rule (i.e., [math]\displaystyle{ \operatorname{d}W/{\operatorname{d}t} }[/math] with [math]\displaystyle{ W \triangleq 2 \sqrt{V} }[/math]), means

[math]\displaystyle{ \mathord{\underbrace{D^+ \Bigl( \mathord{\underbrace{2 \mathord{\overbrace{\sqrt{V}}^{ {} \propto \|\sigma\|_2}}}_{W}} \Bigr)}_{D^+ W \, \triangleq \, \mathord{\text{Upper right-hand } \dot{W}}}} = \frac{ 1 }{ \sqrt{V} } \frac{\operatorname{d}V}{\operatorname{d}t} \leq -\mu }[/math]

where [math]\displaystyle{ D^+ }[/math] is the upper right-hand derivative of [math]\displaystyle{ 2 \sqrt{V} }[/math] and the symbol [math]\displaystyle{ \propto }[/math] denotes proportionality. So, by comparison to the curve [math]\displaystyle{ z(t) = z_0 - \mu t }[/math] which is represented by differential equation [math]\displaystyle{ \dot{z} = -\mu }[/math] with initial condition [math]\displaystyle{ z(0)=z_0 }[/math], it must be the case that [math]\displaystyle{ 2 \sqrt{V(t)} \leq V_0 - \mu t }[/math] for all t. Moreover, because [math]\displaystyle{ \sqrt{V} \geq 0 }[/math], [math]\displaystyle{ \sqrt{V} }[/math] must reach [math]\displaystyle{ \sqrt{V}=0 }[/math] in finite time, which means that V must reach [math]\displaystyle{ V=0 }[/math] (i.e., the system enters the sliding mode) in finite time. Because [math]\displaystyle{ \sqrt{V} }[/math] is proportional to the Euclidean norm [math]\displaystyle{ \|\mathord{\cdot}\|_2 }[/math] of the switching function [math]\displaystyle{ \sigma }[/math], this result implies that the rate of approach to the sliding mode must be firmly bounded away from zero.

Consequences for sliding mode control

In the context of sliding mode control, this condition means that

[math]\displaystyle{ \underbrace{ \overbrace{\sigma^\intercal}^{\tfrac{\partial V}{\partial \sigma}} \overbrace{\dot{\sigma}}^{\tfrac{\operatorname{d} \sigma}{\operatorname{d} t}} }_{\tfrac{\operatorname{d}V}{\operatorname{d}t}} \leq -\mu ( \mathord{\overbrace{\| \sigma \|_2}^{\sqrt{V}}} )^{\alpha} }[/math]

where [math]\displaystyle{ \|\mathord{\cdot}\| }[/math] is the Euclidean norm. For the case when switching function [math]\displaystyle{ \sigma }[/math] is scalar valued, the sufficient condition becomes

[math]\displaystyle{ \sigma \dot{\sigma} \leq -\mu |\sigma|^{\alpha} }[/math].

Taking [math]\displaystyle{ \alpha =1 }[/math], the scalar sufficient condition becomes

[math]\displaystyle{ \operatorname{sgn}(\sigma) \dot{\sigma} \leq -\mu }[/math]

which is equivalent to the condition that

[math]\displaystyle{ \operatorname{sgn}(\sigma) \neq \operatorname{sgn}(\dot{\sigma}) \qquad \text{and} \qquad |\dot{\sigma}| \geq \mu \gt 0 }[/math].

That is, the system should always be moving toward the switching surface [math]\displaystyle{ \sigma = 0 }[/math], and its speed [math]\displaystyle{ |\dot{\sigma}| }[/math] toward the switching surface should have a non-zero lower bound. So, even though [math]\displaystyle{ \sigma }[/math] may become vanishingly small as [math]\displaystyle{ \mathbf{x} }[/math] approaches the [math]\displaystyle{ \sigma(\mathbf{x})=\mathbf{0} }[/math] surface, [math]\displaystyle{ \dot{\sigma} }[/math] must always be bounded firmly away from zero. To ensure this condition, sliding mode controllers are discontinuous across the [math]\displaystyle{ \sigma = 0 }[/math] manifold; they switch from one non-zero value to another as trajectories cross the manifold.

Theorem 2: Region of attraction

For the system given by Equation (1) and sliding surface given by Equation (2), the subspace for which the [math]\displaystyle{ \{ \mathbf{x} \in \mathbb{R}^n : \sigma(\mathbf{x})=\mathbf{0} \} }[/math] surface is reachable is given by

[math]\displaystyle{ \{ \mathbf{x} \in \mathbb{R}^n : \sigma^\intercal(\mathbf{x})\dot{\sigma}(\mathbf{x}) \lt 0 \} }[/math]

That is, when initial conditions come entirely from this space, the Lyapunov function candidate [math]\displaystyle{ V(\sigma) }[/math] is a Lyapunov function and [math]\displaystyle{ \mathbf{x} }[/math] trajectories are sure to move toward the sliding mode surface where [math]\displaystyle{ \sigma( \mathbf{x} ) = \mathbf{0} }[/math]. Moreover, if the reachability conditions from Theorem 1 are satisfied, the sliding mode will enter the region where [math]\displaystyle{ \dot{V} }[/math] is more strongly bounded away from zero in finite time. Hence, the sliding mode [math]\displaystyle{ \sigma = 0 }[/math] will be attained in finite time.

Theorem 3: Sliding motion

Let

[math]\displaystyle{ \frac{\partial \sigma}{\partial{\mathbf{x}}} B(\mathbf{x},t) }[/math]

be nonsingular. That is, the system has a kind of controllability that ensures that there is always a control that can move a trajectory to move closer to the sliding mode. Then, once the sliding mode where [math]\displaystyle{ \sigma(\mathbf{x}) = \mathbf{0} }[/math] is achieved, the system will stay on that sliding mode. Along sliding mode trajectories, [math]\displaystyle{ \sigma(\mathbf{x}) }[/math] is constant, and so sliding mode trajectories are described by the differential equation

[math]\displaystyle{ \dot{\sigma} = \mathbf{0} }[/math].

If an [math]\displaystyle{ \mathbf{x} }[/math]-equilibrium is stable with respect to this differential equation, then the system will slide along the sliding mode surface toward the equilibrium.

The equivalent control law on the sliding mode can be found by solving

[math]\displaystyle{ \dot\sigma(\mathbf{x})=0 }[/math]

for the equivalent control law [math]\displaystyle{ \mathbf{u}(\mathbf{x}) }[/math]. That is,

[math]\displaystyle{ \frac{\partial \sigma}{\partial \mathbf{x}} \overbrace{\left( f(\mathbf{x},t) + B(\mathbf{x},t) \mathbf{u} \right)}^{\dot{\mathbf{x}}} = \mathbf{0} }[/math]

and so the equivalent control

[math]\displaystyle{ \mathbf{u} = -\left( \frac{\partial \sigma}{\partial \mathbf{x}} B(\mathbf{x},t) \right)^{-1} \frac{\partial \sigma}{\partial \mathbf{x}} f(\mathbf{x},t) }[/math]

That is, even though the actual control [math]\displaystyle{ \mathbf{u} }[/math] is not continuous, the rapid switching across the sliding mode where [math]\displaystyle{ \sigma(\mathbf{x})=\mathbf{0} }[/math] forces the system to act as if it were driven by this continuous control.

Likewise, the system trajectories on the sliding mode behave as if

[math]\displaystyle{ \dot{\mathbf{x}} = \overbrace{f(\mathbf{x},t) - B(\mathbf{x},t) \left( \frac{\partial \sigma}{\partial \mathbf{x}} B(\mathbf{x},t) \right)^{-1} \frac{\partial \sigma}{\partial \mathbf{x}} f(\mathbf{x},t)}^{f(\mathbf{x},t) + B(\mathbf{x},t) u} = f(\mathbf{x},t)\left( \mathbf{I} - B(\mathbf{x},t) \left( \frac{\partial \sigma}{\partial \mathbf{x}} B(\mathbf{x},t) \right)^{-1} \frac{\partial \sigma}{\partial \mathbf{x}} \right) }[/math]

The resulting system matches the sliding mode differential equation

[math]\displaystyle{ \dot{\sigma}(\mathbf{x}) = \mathbf{0} }[/math]

, the sliding mode surface [math]\displaystyle{ \sigma(\mathbf{x})=\mathbf{0} }[/math], and the trajectory conditions from the reaching phase now reduce to the above derived simpler condition. Hence, the system can be assumed to follow the simpler [math]\displaystyle{ \dot{\sigma} = 0 }[/math] condition after some initial transient during the period while the system finds the sliding mode. The same motion is approximately maintained when the equality [math]\displaystyle{ \sigma(\mathbf{x}) = \mathbf{0} }[/math] only approximately holds.

It follows from these theorems that the sliding motion is invariant (i.e., insensitive) to sufficiently small disturbances entering the system through the control channel. That is, as long as the control is large enough to ensure that [math]\displaystyle{ \sigma^\intercal \dot{\sigma} \lt 0 }[/math] and [math]\displaystyle{ \dot{\sigma} }[/math] is uniformly bounded away from zero, the sliding mode will be maintained as if there was no disturbance. The invariance property of sliding mode control to certain disturbances and model uncertainties is its most attractive feature; it is strongly robust.

As discussed in an example below, a sliding mode control law can keep the constraint

[math]\displaystyle{ \dot{x} + x = 0 }[/math]

in order to asymptotically stabilize any system of the form

[math]\displaystyle{ \ddot{x}=a(t,x,\dot{x}) + u }[/math]

when [math]\displaystyle{ a(\cdot) }[/math] has a finite upper bound. In this case, the sliding mode is where

[math]\displaystyle{ \dot{x} = -x }[/math]

(i.e., where [math]\displaystyle{ \dot{x}+x=0 }[/math]). That is, when the system is constrained this way, it behaves like a simple stable linear system, and so it has a globally exponentially stable equilibrium at the [math]\displaystyle{ (x,\dot{x})=(0,0) }[/math] origin.

Control design examples

• Consider a plant described by Equation (1) with single input u (i.e., [math]\displaystyle{ m = 1 }[/math]). The switching function is picked to be the linear combination

[math]\displaystyle{ \sigma(\mathbf{x}) \triangleq s_1 x_1 + s_2 x_2 + \cdots + s_{n-1} x_{n-1} + s_n x_n }[/math] ( 4)

where the weight [math]\displaystyle{ s_i \gt 0 }[/math] for all [math]\displaystyle{ 1 \leq i \leq n }[/math]. The sliding surface is the simplex where [math]\displaystyle{ \sigma(\mathbf{x})=0 }[/math]. When trajectories are forced to slide along this surface,

[math]\displaystyle{ \dot{\sigma}(\mathbf{x}) = 0 }[/math]

and so

[math]\displaystyle{ s_1 \dot{x}_1 + s_2 \dot{x}_2 + \cdots + s_{n-1} \dot{x}_{n-1} + s_n \dot{x}_n = 0 }[/math]

which is a reduced-order system (i.e., the new system is of order [math]\displaystyle{ n-1 }[/math] because the system is constrained to this [math]\displaystyle{ (n-1) }[/math]-dimensional sliding mode simplex). This surface may have favorable properties (e.g., when the plant dynamics are forced to slide along this surface, they move toward the origin [math]\displaystyle{ \mathbf{x}=\mathbf{0} }[/math]). Taking the derivative of the Lyapunov function in Equation (3), we have

[math]\displaystyle{ \dot{V}(\sigma(\mathbf{x})) = \overbrace{\sigma(\mathbf{x})^{\text{T}}}^{\tfrac{\partial V}{\partial \mathbf{x}}} \overbrace{\dot{\sigma}(\mathbf{x})}^{\tfrac{\operatorname{d} \sigma}{\operatorname{d} t}} }[/math]

To ensure [math]\displaystyle{ \dot{V} \lt 0 }[/math], the feedback control law [math]\displaystyle{ u(\mathbf{x}) }[/math] must be chosen so that

[math]\displaystyle{ \begin{cases} \dot{\sigma} \lt 0 &\text{if } \sigma \gt 0\\ \dot{\sigma} \gt 0 &\text{if } \sigma \lt 0 \end{cases} }[/math]

Hence, the product [math]\displaystyle{ \sigma \dot{\sigma} \lt 0 }[/math] because it is the product of a negative and a positive number. Note that

[math]\displaystyle{ \dot{\sigma}(\mathbf{x}) = \overbrace{\frac{\partial{\sigma(\mathbf{x})}}{\partial{\mathbf{x}}} \dot{\mathbf{x}}}^{\dot{\sigma}(\mathbf{x})} = \frac{\partial{\sigma(\mathbf{x})}}{\partial{\mathbf{x}}} \overbrace{\left( f(\mathbf{x},t) + B(\mathbf{x},t) u \right)}^{\dot{\mathbf{x}}} = \overbrace{[s_1, s_2, \ldots, s_n]}^{\frac{\partial{\sigma(\mathbf{x})}}{\partial{\mathbf{x}}}} \underbrace{\overbrace{\left( f(\mathbf{x},t) + B(\mathbf{x},t) u \right)}^{\dot{\mathbf{x}}}}_{\text{( i.e., an } n \times 1 \text{ vector )}} }[/math] ( 5)

The control law [math]\displaystyle{ u(\mathbf{x}) }[/math] is chosen so that

[math]\displaystyle{ u(\mathbf{x}) = \begin{cases} u^+(\mathbf{x}) &\text{if } \sigma(\mathbf{x}) \gt 0 \\ u^-(\mathbf{x}) &\text{if } \sigma(\mathbf{x}) \lt 0 \end{cases} }[/math]

where

• [math]\displaystyle{ u^+(\mathbf{x}) }[/math] is some control (e.g., possibly extreme, like "on" or "forward") that ensures Equation (5) (i.e., [math]\displaystyle{ \dot{\sigma} }[/math]) is negative at [math]\displaystyle{ \mathbf{x} }[/math]
• [math]\displaystyle{ u^-(\mathbf{x}) }[/math] is some control (e.g., possibly extreme, like "off" or "reverse") that ensures Equation (5) (i.e., [math]\displaystyle{ \dot{\sigma} }[/math]) is positive at [math]\displaystyle{ \mathbf{x} }[/math]

The resulting trajectory should move toward the sliding surface where [math]\displaystyle{ \sigma(\mathbf{x})=0 }[/math]. Because real systems have delay, sliding mode trajectories often chatter back and forth along this sliding surface (i.e., the true trajectory may not smoothly follow [math]\displaystyle{ \sigma(\mathbf{x})=0 }[/math], but it will always return to the sliding mode after leaving it).

• Consider the dynamic system

[math]\displaystyle{ \ddot{x}=a(t,x,\dot{x})+u }[/math]

which can be expressed in a 2-dimensional state space (with [math]\displaystyle{ x_1 = x }[/math] and [math]\displaystyle{ x_2 = \dot{x} }[/math]) as

[math]\displaystyle{ \begin{cases} \dot{x}_1 = x_2\\ \dot{x}_2 = a(t,x_1,x_2) + u \end{cases} }[/math]

Also assume that [math]\displaystyle{ \sup\{ |a(\cdot)| \} \leq k }[/math] (i.e., [math]\displaystyle{ |a| }[/math] has a finite upper bound k that is known). For this system, choose the switching function

[math]\displaystyle{ \sigma(x_1,x_2)= x_1 + x_2 = x + \dot{x} }[/math]

By the previous example, we must choose the feedback control law [math]\displaystyle{ u(x,\dot{x}) }[/math] so that [math]\displaystyle{ \sigma \dot{\sigma} \lt 0 }[/math]. Here,

[math]\displaystyle{ \dot{\sigma} = \dot{x}_1 + \dot{x}_2 = \dot{x} + \ddot{x} = \dot{x}\,+\,\overbrace{a(t,x,\dot{x})+ u}^{\ddot{x}} }[/math]

• When [math]\displaystyle{ x + \dot{x} \lt 0 }[/math] (i.e., when [math]\displaystyle{ \sigma \lt 0 }[/math]), to make [math]\displaystyle{ \dot{\sigma} \gt 0 }[/math], the control law should be picked so that [math]\displaystyle{ u \gt |\dot{x} + a(t,x,\dot{x})| }[/math]
• When [math]\displaystyle{ x + \dot{x} \gt 0 }[/math] (i.e., when [math]\displaystyle{ \sigma \gt 0 }[/math]), to make [math]\displaystyle{ \dot{\sigma} \lt 0 }[/math], the control law should be picked so that [math]\displaystyle{ u \lt -|\dot{x} + a(t,x,\dot{x})| }[/math]

However, by the triangle inequality,

[math]\displaystyle{ |\dot{x}| + |a(t,x,\dot{x})| \geq |\dot{x} + a(t,x,\dot{x})| }[/math]

and by the assumption about [math]\displaystyle{ |a| }[/math],

[math]\displaystyle{ |\dot{x}| + k + 1 \gt |\dot{x}| + |a(t,x,\dot{x})| }[/math]

So the system can be feedback stabilized (to return to the sliding mode) by means of the control law

[math]\displaystyle{ u(x,\dot{x}) = \begin{cases} |\dot{x}| + k + 1 &\text{if } \underbrace{x + \dot{x}} \lt 0,\\ -\left(|\dot{x}| + k + 1\right) &\text{if } \overbrace{x + \dot{x}}^{\sigma} \gt 0 \end{cases} }[/math]

which can be expressed in closed form as

[math]\displaystyle{ u(x,\dot{x}) = -(|\dot{x}|+k+1) \underbrace{\operatorname{sgn}(\overbrace{\dot{x}+x}^{\sigma})}_{\text{(i.e., tests } \sigma \gt 0 \text{)}} }[/math]

Assuming that the system trajectories are forced to move so that [math]\displaystyle{ \sigma(\mathbf{x})=0 }[/math], then

[math]\displaystyle{ \dot{x} = -x \qquad \text{(i.e., } \sigma(x,\dot{x}) = x + \dot{x} = 0 \text{)} }[/math]

So once the system reaches the sliding mode, the system's 2-dimensional dynamics behave like this 1-dimensional system, which has a globally exponentially stable equilibrium at [math]\displaystyle{ (x,\dot{x})=(0,0) }[/math].

Automated design solutions

Although various theories exist for sliding mode control system design, there is a lack of a highly effective design methodology due to practical difficulties encountered in analytical and numerical methods. A reusable computing paradigm such as a genetic algorithm can, however, be utilized to transform a 'unsolvable problem' of optimal design into a practically solvable 'non-deterministic polynomial problem'. This results in computer-automated designs for sliding model control.

Sliding mode observer

Sliding mode control can be used in the design of state observers. These non-linear high-gain observers have the ability to bring coordinates of the estimator error dynamics to zero in finite time. Additionally, switched-mode observers have attractive measurement noise resilience that is similar to a Kalman filter. For simplicity, the example here uses a traditional sliding mode modification of a Luenberger observer for an LTI system. In these sliding mode observers, the order of the observer dynamics are reduced by one when the system enters the sliding mode. In this particular example, the estimator error for a single estimated state is brought to zero in finite time, and after that time the other estimator errors decay exponentially to zero. However, as first described by Drakunov, a sliding mode observer for non-linear systems can be built that brings the estimation error for all estimated states to zero in a finite (and arbitrarily small) time.

Here, consider the LTI system

[math]\displaystyle{ \begin{cases} \dot{\mathbf{x}} = A \mathbf{x} + B \mathbf{u}\\ y = \begin{bmatrix}1 & 0 & 0 & \cdots & \end{bmatrix} \mathbf{x} = x_1 \end{cases} }[/math]

where state vector [math]\displaystyle{ \mathbf{x} \triangleq (x_1, x_2, \dots, x_n) \in \mathbb{R}^n }[/math], [math]\displaystyle{ \mathbf{u} \triangleq (u_1, u_2, \dots, u_r) \in \mathbb{R}^r }[/math] is a vector of inputs, and output y is a scalar equal to the first state of the [math]\displaystyle{ \mathbf{x} }[/math] state vector. Let

[math]\displaystyle{ A \triangleq \begin{bmatrix} a_{11} & A_{12} \\ A_{21} & A_{22} \end{bmatrix} }[/math]

where

• [math]\displaystyle{ a_{11} }[/math] is a scalar representing the influence of the first state [math]\displaystyle{ x_1 }[/math] on itself,
• [math]\displaystyle{ A_{21} \in \mathbb{R}^{(n-1)} }[/math] is a row vector corresponding to the influence of the first state on the other states,
• [math]\displaystyle{ A_{22} \in \mathbb{R}^{(n-1) \times (n-1)} }[/math] is a matrix representing the influence of the other states on themselves, and
• [math]\displaystyle{ A_{12} \in \mathbb{R}^{1\times(n-1)} }[/math] is a column vector representing the influence of the other states on the first state.

The goal is to design a high-gain state observer that estimates the state vector [math]\displaystyle{ \mathbf{x} }[/math] using only information from the measurement [math]\displaystyle{ y=x_1 }[/math]. Hence, let the vector [math]\displaystyle{ \hat{\mathbf{x}} = (\hat{x}_1,\hat{x}_2,\dots,\hat{x}_n) \in \mathbb{R}^n }[/math] be the estimates of the n states. The observer takes the form

[math]\displaystyle{ \dot{\hat{\mathbf{x}}} = A \hat{\mathbf{x}} + B \mathbf{u} + L v(\hat{x}_1 - x_1) }[/math]

where [math]\displaystyle{ v: \R \to \R }[/math] is a nonlinear function of the error between estimated state [math]\displaystyle{ \hat{x}_1 }[/math] and the output [math]\displaystyle{ y=x_1 }[/math], and [math]\displaystyle{ L \in \mathbb{R}^n }[/math] is an observer gain vector that serves a similar purpose as in the typical linear Luenberger observer. Likewise, let

[math]\displaystyle{ L = \begin{bmatrix} -1 \\ L_{2} \end{bmatrix} }[/math]

where [math]\displaystyle{ L_2 \in \mathbb{R}^{(n-1)} }[/math] is a column vector. Additionally, let [math]\displaystyle{ \mathbf{e} = (e_1, e_2, \dots, e_n) \in \mathbb{R}^n }[/math] be the state estimator error. That is, [math]\displaystyle{ \mathbf{e} = \hat{\mathbf{x}} - \mathbf{x} }[/math]. The error dynamics are then

[math]\displaystyle{ \begin{align} \dot{\mathbf{e}} &= \dot{\hat{\mathbf{x}}} - \dot{\mathbf{x}}\\ &= A \hat{\mathbf{x}} + B \mathbf{u} + L v(\hat{x}_1 - x_1) - A \mathbf{x} - B \mathbf{u}\\ &= A (\hat{\mathbf{x}} - \mathbf{x}) + L v(\hat{x}_1 - x_1)\\ &= A \mathbf{e} + L v(e_1) \end{align} }[/math]

where [math]\displaystyle{ e_1 = \hat{x}_1 - x_1 }[/math] is the estimator error for the first state estimate. The nonlinear control law v can be designed to enforce the sliding manifold

[math]\displaystyle{ 0 = \hat{x}_1 - x_1 }[/math]

so that estimate [math]\displaystyle{ \hat{x}_1 }[/math] tracks the real state [math]\displaystyle{ x_1 }[/math] after some finite time (i.e., [math]\displaystyle{ \hat{x}_1 = x_1 }[/math]). Hence, the sliding mode control switching function

[math]\displaystyle{ \sigma(\hat{x}_1,\hat{x}) \triangleq e_1 = \hat{x}_1 - x_1. }[/math]

To attain the sliding manifold, [math]\displaystyle{ \dot{\sigma} }[/math] and [math]\displaystyle{ \sigma }[/math] must always have opposite signs (i.e., [math]\displaystyle{ \sigma \dot{\sigma} \lt 0 }[/math] for essentially all [math]\displaystyle{ \mathbf{x} }[/math]). However,

[math]\displaystyle{ \dot{\sigma} = \dot{e}_1 = a_{11} e_1 + A_{12} \mathbf{e}_2 - v( e_1 ) = a_{11} e_1 + A_{12} \mathbf{e}_2 - v( \sigma ) }[/math]

where [math]\displaystyle{ \mathbf{e}_2 \triangleq (e_2, e_3, \ldots, e_n) \in \mathbb{R}^{(n-1)} }[/math] is the collection of the estimator errors for all of the unmeasured states. To ensure that [math]\displaystyle{ \sigma \dot{\sigma} \lt 0 }[/math], let

[math]\displaystyle{ v( \sigma ) = M \operatorname{sgn}(\sigma) }[/math]

where

[math]\displaystyle{ M \gt \max\{ |a_{11} e_1 + A_{12} \mathbf{e}_2| \}. }[/math]

That is, positive constant M must be greater than a scaled version of the maximum possible estimator errors for the system (i.e., the initial errors, which are assumed to be bounded so that M can be picked large enough; al). If M is sufficiently large, it can be assumed that the system achieves [math]\displaystyle{ e_1 = 0 }[/math] (i.e., [math]\displaystyle{ \hat{x}_1 = x_1 }[/math]). Because [math]\displaystyle{ e_1 }[/math] is constant (i.e., 0) along this manifold, [math]\displaystyle{ \dot{e}_1 = 0 }[/math] as well. Hence, the discontinuous control [math]\displaystyle{ v(\sigma) }[/math] may be replaced with the equivalent continuous control [math]\displaystyle{ v_{\text{eq}} }[/math] where

[math]\displaystyle{ 0 = \dot{\sigma} = a_{11} \mathord{\overbrace{e_1}^{ {} = 0 }} + A_{12} \mathbf{e}_2 - \mathord{\overbrace{v_{\text{eq}}}^{v(\sigma)}} = A_{12} \mathbf{e}_2 - v_{\text{eq}}. }[/math]

So

[math]\displaystyle{ \mathord{\underbrace{v_{\text{eq}}}_{\text{scalar}}} = \mathord{\underbrace{A_{12}}_{1 \times (n-1) \atop \text{ vector}}} \mathord{\underbrace{\mathbf{e}_2}_{(n-1) \times 1 \atop \text{ vector}}}. }[/math]

This equivalent control [math]\displaystyle{ v_{\text{eq}} }[/math] represents the contribution from the other [math]\displaystyle{ (n-1) }[/math] states to the trajectory of the output state [math]\displaystyle{ x_1 }[/math]. In particular, the row [math]\displaystyle{ A_{12} }[/math] acts like an output vector for the error subsystem

[math]\displaystyle{ \mathord{\overbrace{ \begin{bmatrix} \dot{e}_2\\ \dot{e}_3\\ \vdots\\ \dot{e}_n \end{bmatrix} }^{\dot{\mathbf{e}}_2}} = A_2 \mathord{\overbrace{ \begin{bmatrix} e_2\\ e_3\\ \vdots\\ e_n \end{bmatrix} }^{\mathbf{e}_2}} + L_2 v(e_1) = A_2 \mathbf{e}_2 + L_2 v_{\text{eq}} = A_2 \mathbf{e}_2 + L_2 A_{12} \mathbf{e}_2 = ( A_2 + L_2 A_{12} ) \mathbf{e}_2. }[/math]

So, to ensure the estimator error [math]\displaystyle{ \mathbf{e}_2 }[/math] for the unmeasured states converges to zero, the [math]\displaystyle{ (n-1)\times 1 }[/math] vector [math]\displaystyle{ L_2 }[/math] must be chosen so that the [math]\displaystyle{ (n-1)\times (n-1) }[/math] matrix [math]\displaystyle{ ( A_2 + L_2 A_{12} ) }[/math] is Hurwitz (i.e., the real part of each of its eigenvalues must be negative). Hence, provided that it is observable, this [math]\displaystyle{ \mathbf{e}_2 }[/math] system can be stabilized in exactly the same way as a typical linear state observer when [math]\displaystyle{ A_{12} }[/math] is viewed as the output matrix (i.e., "C"). That is, the [math]\displaystyle{ v_{\text{eq}} }[/math] equivalent control provides measurement information about the unmeasured states that can continually move their estimates asymptotically closer to them. Meanwhile, the discontinuous control [math]\displaystyle{ v = M \operatorname{sgn}( \hat{x}_1 - x ) }[/math] forces the estimate of the measured state to have zero error in finite time. Additionally, white zero-mean symmetric measurement noise (e.g., Gaussian noise) only affects the switching frequency of the control v, and hence the noise will have little effect on the equivalent sliding mode control [math]\displaystyle{ v_{\text{eq}} }[/math]. Hence, the sliding mode observer has Kalman filter–like features.

The final version of the observer is thus

[math]\displaystyle{ \begin{align} \dot{\hat{\mathbf{x}}} &= A \hat{\mathbf{x}} + B \mathbf{u} + L M \operatorname{sgn}(\hat{x}_1 - x_1)\\ &= A \hat{\mathbf{x}} + B \mathbf{u} + \begin{bmatrix} -1\\L_2 \end{bmatrix} M \operatorname{sgn}(\hat{x}_1 - x_1)\\ &= A \hat{\mathbf{x}} + B \mathbf{u} + \begin{bmatrix} -M\\L_2 M\end{bmatrix} \operatorname{sgn}(\hat{x}_1 - x_1)\\ &= A \hat{\mathbf{x}} + \begin{bmatrix} B & \begin{bmatrix} -M\\L_2 M\end{bmatrix} \end{bmatrix} \begin{bmatrix} \mathbf{u} \\ \operatorname{sgn}(\hat{x}_1 - x_1) \end{bmatrix}\\ &= A_{\text{obs}} \hat{\mathbf{x}} + B_{\text{obs}} \mathbf{u}_{\text{obs}} \end{align} }[/math]

where

• [math]\displaystyle{ A_{\text{obs}} \triangleq A, }[/math]
• [math]\displaystyle{ B_{\text{obs}} \triangleq \begin{bmatrix} B & \begin{bmatrix} -M\\L_2 M\end{bmatrix} \end{bmatrix}, }[/math] and
• [math]\displaystyle{ u_{\text{obs}} \triangleq \begin{bmatrix} \mathbf{u} \\ \operatorname{sgn}(\hat{x}_1 - x_1) \end{bmatrix}. }[/math]

That is, by augmenting the control vector [math]\displaystyle{ \mathbf{u} }[/math] with the switching function [math]\displaystyle{ \operatorname{sgn}(\hat{x}_1-x_1) }[/math], the sliding mode observer can be implemented as an LTI system. That is, the discontinuous signal [math]\displaystyle{ \operatorname{sgn}(\hat{x}_1-x_1) }[/math] is viewed as a control input to the 2-input LTI system.

For simplicity, this example assumes that the sliding mode observer has access to a measurement of a single state (i.e., output [math]\displaystyle{ y=x_1 }[/math]). However, a similar procedure can be used to design a sliding mode observer for a vector of weighted combinations of states (i.e., when output [math]\displaystyle{ \mathbf{y} = C \mathbf{x} }[/math] uses a generic matrix C). In each case, the sliding mode will be the manifold where the estimated output [math]\displaystyle{ \hat{\mathbf{y}} }[/math] follows the measured output [math]\displaystyle{ \mathbf{y} }[/math] with zero error (i.e., the manifold where [math]\displaystyle{ \sigma(\mathbf{x}) \triangleq \hat{\mathbf{y}} - \mathbf{y} = \mathbf{0} }[/math]).

See also

• Variable structure control
• Variable structure system
• Hybrid system
• Nonlinear control
• Robust control
• Optimal control
• Bang–bang control – Sliding mode control is often implemented as a bang–bang control. In some cases, such control is necessary for optimality.
• H-bridge – A topology that combines four switches forming the four legs of an "H". Can be used to drive a motor (or other electrical device) forward or backward when only a single supply is available. Often used in actuator in sliding-mode controlled systems.
• Switching amplifier – Uses switching-mode control to drive continuous outputs
• Delta-sigma modulation – Another (feedback) method of encoding a continuous range of values in a signal that rapidly switches between two states (i.e., a kind of specialized sliding-mode control)
• Pulse-density modulation – A generalized form of delta-sigma modulation.
• Pulse-width modulation – Another modulation scheme that produces continuous motion through discontinuous switching.

References

• Zinober, Alan S.I., ed (1994). Variable Structure and Lyapunov Control. Lecture Notes in Control and Information Sciences. 193. London: Springer-Verlag. doi:10.1007/BFb0033675. ISBN 978-3-540-19869-7. 

• Edwards, C.; Spurgeon, S. (1998). Sliding Mode Control. Theory and Applications. London: Taylor and Francis. ISBN 978-0-7484-0601-2. 
• Sabanovic, Asif; Fridman, Leonid; Spurgeon, Sarah (2004). Variable Structure Systems: From Principles to Implementation. London: The Insitite of Electrical Engineers. ISBN 0-86341-350-1. 
• Edwards, C.; Fossas Colet, E.; Fridman, L. (2006). Advances in Variable Structure and Sliding Mode Control. Berlin: Springer Verlag. ISBN 978-3-540-32800-1. 
• Bartolini, G.; Fridman, L.; Pisano, A.; Usai, E. (2008). Modern Sliding Mode Control Theory - New perspectives and applications.. Berlin: Springer Verlag. ISBN 978-3-540-79016-7. 
• Fridman, L.; Barbot, J.-P.; Plestan, F. (2016). Recent Trends in Sliding Mode Control. London: The Institution of Engineering and Technology. ISBN 978-1-78561-076-9. 
• Shtessel, Y.; Edwards, C.; Fridman, L.; Levant, A. (2014). Sliding Mode Control and Observation. Basel: Birkhauser. ISBN 978-0-81764-8923. 
• Sira- Ramirez, Hebertt (2015). Sliding Mode Control. The Delta-Sigma Modulation Approach. Basel: Birkhauser. ISBN 978-3-319-17256-9. 
• Fridman, L.; Moreno, J.; Bandyopadhyay, B.; Kamal Asif Chalanga, S.; Chalanga, S. (2015). Continuous Nested Algorithms: The Fifth Generation of Sliding Mode Controllers. In: Recent Advances in Sliding Modes: From Control to Intelligent Mechatronics, X. Yu, O. Efe (Eds.), pp. 5–35. Studies in Systems, Decision and Control. 24. London: Springer-Verlag. doi:10.1007/978-3-319-18290-2_2. ISBN 978-3-319-18289-6. 
• Li, S.; Yu, X.; Fridman, L.; Man, Z.; Wang, X. (2017). Advances in Variable Structure Systems and Sliding Mode Control—Theory and Applications.. Studies in Systems, Decision and Control. 24. London: Springer-Verlag. ISBN 978-3-319-62895-0. 
• Ferrara, A.; Incremona, G.P.; Cucuzzella, M. (2019). Advanced and Optimization Based Sliding Mode Control. SIAM. doi:10.1137/1.9781611975840. ISBN 978-1611975833. 
• Steinberger, M.; Horn, M.; Fridman, L., eds (2020). Variable-Structure Systems and Sliding-Mode Control. Studies in Systems, Decision and Control. 271. London: Springer-Verlag. doi:10.1007/BFb0033675. ISBN 978-3-030-36620-9. 

Further reading

• Drakunov S.V., Utkin V.I.. (1992). "Sliding mode control in dynamic systems". International Journal of Control 55 (4): 1029–1037. doi:10.1080/00207179208934270.