1. Introduction
The term sliding mode has first appeared in the theory of variable structure systems [
1] and referred to a characteristic form of nonlinear, discontinuous, feedback control, as described in [
2,
3]. This comprehensive control scheme may be successfully applied to both linear [
4] and nonlinear systems [
5,
6]. In general, the sliding mode control process consists of two stages [
7]. In the first stage, the so-called reaching phase, the system’s state is driven to a predefined surface in the state space. In the second stage, the so-called sliding phase, the representative point of the system slides along this surface until it reaches the equilibrium point. Various authors discussed the sliding mode design [
8,
9] and its applications [
10] in the literature. Sliding modes are especially widely applied in power electronics, for example, for electronic converters [
11,
12], electrical motor drives [
13], and power systems [
14]. The greatest advantage of the sliding mode control is that once the system’s state becomes bound to the sliding plane, any matched external disturbances and model uncertainties are immediately rejected [
15]. Therefore, in order to improve the system’s robustness, it is highly recommended to make the reaching phase as short as possible [
16]. However, enforcing faster convergence often requires large magnitudes of the control input and therefore, consumes more energy. This observation resulted in the development of time-varying sliding planes.
The time-varying sliding planes have been widely applied for continuous-time systems since the 1990s. Their basic design assumes that the sliding plane must cross the initial position of the system’s representative point at time
t = 0 and for any
t > 0 it shifts in the state space, according to different algorithms, until it crosses the system’s desired position at
t =
t0. Afterward, the sliding plane remains still. As follows, the sliding variable is defined as the distance between the representative point of the system and the sliding plane, and the controller’s task is to keep the system’s state on the sliding plane from the very beginning of the control process. Consequently, the reaching phase is completely eliminated and the benefit of insensitivity to external disturbances and modeling uncertainties is observed in the entire course of motion. The first time-varying sliding plane designs combined rotating [
17] and shifting surfaces [
18], which, for example, were applied in the position control of a DC motor drive [
19].
With the rapid development of digital control systems, most of the sliding mode properties have been transferred to the discrete-time domain resulting in the quasi-sliding mode. As shown by Milosavljevic [
20], due to delays in the control channel, in the discrete-time, the representative point of the system only remains in a certain vicinity of the sliding plane, not on the plane itself, creating the so-called quasi-sliding mode band. The properties [
21] and stability of the quasi-sliding mode [
22,
23] have been studied in numerous further works. Some of the authors concentrated on the equivalent control method [
24], others considered the control saturation [
25], while another group utilized the Lyapunov function-based design [
26]. However, a clear definition of the quasi-sliding mode was first introduced in the seminal study of Gao, Wang, and Homaifa in 1995 [
27]. According to their paper, the quasi-sliding mode exists when the sliding plane is reached and crossed in finite time. After the first crossing of the sliding plane, the sign of the sliding variable changes in each step and its absolute value does not exceed an a priori known value. This definition constituted a solid foundation for future quasi-sliding mode research and will be adopted further in this study. The work of Gao et al. [
27] also established a family of switching type reaching laws, which provide the characteristic zigzagging motion around the sliding plane. Afterward, many new reaching laws have been introduced, such as [
28,
29]. Their generalized form was presented in [
30]. A direct measure of their effectiveness and the system’s robustness has become the width of the achieved band. It was proved in [
31] that for discretized systems, the ultimate bandwidth directly depends on the discretization period and therefore is of
O(
T) order.
Along with the development of digital sliding mode control, the issues connected with the energetical efficiency of control systems have come into the spotlight. The optimal control design has grown in popularity [
32]. As energy saving is one of the main requirements of nowadays control systems, and especially of power electronics systems, it has become an important problem in the quasi-sliding mode design. The obtained switching type of motion around the sliding plane directly causes larger energy consumption of the closed-loop system. Therefore, reducing the width of the quasi-sliding mode band, and as follows reducing the necessary control effort, has become the main area of interest. For that purpose, different authors attempted to design new reaching functions [
33], eliminate the chattering effects [
34,
35,
36], or employ various disturbance compensators [
37,
38]. However, as in the discrete-time, the controller is unable to keep the system’s state exactly on the sliding plane, the idea of time-varying sliding surfaces has not been widely applied in this case.
Sliding mode control has become the center of our efforts in this paper as well. We propose a new time-varying sliding mode controller for sampled-data systems. The originality of this work is manifested in a new, simple but fairly efficient, design of the time-varying sliding surface. It is generated so that the representative point of the system belongs to the surface at the initial step
k = 0. Next, for 0 <
k <
k0, the sliding plane shifts according to a pre-generated reference sliding variable profile
sr(
k), monotonically approaching the system’s demand position. Once, the demand position is crossed, at
k =
k0, the plane shifts by ±
β in each step, ensuring the switching requirement, stated by Gao et al. [
27]. For the generation of such a sliding surface, a reference sliding variable
sr(
k) and a new reaching law are utilized. Next, the control signal for the real perturbed system is designed so that it drives the representative point onto the sliding plane at each time instant from the very beginning of the control process. Our study shows that this control scheme eliminates the reaching phase, provides all the quasi-sliding mode properties, as defined in [
27], including the switching condition, and minimizes the ultimate band-width. In other words, the proposed control method reduces chattering, which results in a significant reduction of energy consumption in the sliding phase. This achievement is particularly important for the control of electric motor drives, where energetical efficiency is one of the main factors in determining the control scheme. Supplemented by a disturbance compensation scheme, our strategy reduces the width of the achieved quasi-sliding mode band from
O(
T) to
O(
T3) order. Moreover, it is proved that the maximum deviations of all state variables from their demand values become of
O(
T3) order as well. In short, our control method improves the robustness of the system in comparison to other control methods, previously proposed in the literature, and increases the energetical efficiency of the closed-loop system.
2. Problem Formulation
This paper presents the design of a sliding mode controller for a linear, continuous-time system subject to matched external disturbances
where appropriate symbols represent:
x(
t)–the system’s
n × 1 state vector,
Φ–the
n ×
n state matrix,
Γ–the
n × 1 input vector,
u(
t)–the control signal and
f(
t)–the disturbance. The aim of the control algorithm is to drive the system from its initial state, denoted with
x(0) =
x0, to the desired position
xd in finite time. The control is feasible, as we assume that the disturbance
f(
t) is lower and upper bounded by ±
fmax, respectively. Moreover, the study is focused on the case when
f(
t) is a smooth function of time and satisfies the matching condition, as in relation (1). As it is most common that all state measurements and the control implementation include sample and hold devices in their channels, the system will be analyzed in the discrete-time domain. Therefore, with sampling period
T, the discretized system becomes
where appropriate matrices are obtained according to the following
Equation (3) indicates that the discrete-time system properties directly depend on the choice of the discretization period
T. Therefore, the obtained input distribution vector
b, as well as the discretized disturbance
d(
k), are of
O(
T) order. Further properties of the discretized disturbance vector have been studied in detail in the noteworthy papers [
31,
37]. With the application of Taylor’s series expansion, the disturbance vector
d(
k) may be expressed as
where
v(
t) =
df(
t)/
dt and
w(
t) =
d2f(
t)/
dt2 are the first and second derivatives of the disturbance and
μ is some moment between the consecutive time steps
kT and (
k + 1)
T. As presented in [
37], from the solution of (4), we get
The first term in (5) proves that the disturbance
d(
k) linearly depends on the discretization period
T, as acknowledged before. However, (5) also shows that in the discrete-time domain, the disturbance contains a mismatched term of
O(
T3) order, represented by the integral in (5). Finally,
d(
k) may be approximated with
Moreover, as the continuous-time disturbance is bounded by ±
fmax, the discretized disturbance range is limited as well. Finally, considering the disturbance rate of change we get
and
The aim of this paper is to present a new unconventional quasi-sliding mode design for the presented system. The key idea behind this work is the implementation of a reference sliding variable profile to determine the time-varying sliding surface evolution, which is described in the following chapter. Next, a discrete-time control law for the original system is proposed, which restricts the system’s state to the quasi-sliding mode band from the beginning of the control process, therefore eliminating the reaching phase. Finally, we introduce a disturbance compensation scheme, which results in the restriction of the sliding mode bandwidth and the state variable error values to O(T3) order. As follows, the proposed control scheme reduces the necessary control effort and energy consumption in the whole regulation process. The paper is concluded with a simulation example.
3. Time-Varying Quasi-SMC
In this chapter, we propose a new discrete-time sliding mode control strategy based on a time-varying sliding surface.
Section 1 presents the generation of the desired sliding surface evolution, employing a reference sliding variable. Next, the control law for the original system (2) is presented and the properties of the resulting sliding mode are analyzed. Finally, in
Section 3 an additional compensation scheme for the matched external disturbances is introduced, which further reduces the quasi-sliding mode bandwidth and therefore, reduces the control effort in the sliding stage. The structure of the proposed closed-loop system is presented in
Figure 1 below.
3.1. Time-Varying Sliding Plane Design
The first step in the sliding mode control design is the proper selection of the sliding hyperplane. In order to limit the impact of external disturbances on the system’s trajectory, a time-varying sliding surface will be utilized. To eliminate the reaching phase, the surface must be designed so that it crosses the initial position of the system’s representative point at
k = 0 and for any
k > 0 it shifts until it crosses the desired state at
k =
k0. Such sliding surface is defined with
where
c is an 1 ×
n control vector, chosen so that
cb ≠ 0 and to ensure stability of the steady state and
e(
k) =
xd −
x(
k) represents the error vector. The first term in (9),
ce(
k) =
s(
k), defines the sliding variable of the system, and
sr(
k) is its pre-generated reference profile. It combines the reference profiles of all state variables of the system as follows
where
xr(
k) represents the reference state vector and
x0 =
xr(0). Therefore, the initial value of the reference sliding variable satisfies
sr(0) =
s(0) =
ce(0), so the system’s representative point is on the sliding plane at the beginning of the control process. The sliding variable defined in (9) physically represents the distance between the current position of the representative point of the system and its desired position in the same step. The task of the control scheme is to keep the system’s state on the sliding plane through the whole regulation process or, in other words, to steer the system’s state vector
x(
k) to the reference
xr(
k). From the initial condition, we conclude that at the beginning of the control process for
k = 0, the representative point of the system belongs to the sliding surface. Therefore, the reachability of the desired position
xd depends on the sliding surface evolution. As the properties of this surface directly depend on the reference sliding variable profile (10), the algorithm for its proper generation will be discussed next.
As the reference sliding variable profile defines the evolution of the sliding surface, it must begin with the initial condition
sr(0) =
ce(0) and monotonically approach zero. Moreover, as the switching type definition of the sliding mode should be satisfied, once the reference crosses zero, its sign must change in each consecutive time instant. To provide the above properties, we propose to generate the reference sliding variable profile according to a new switching type reaching law
where
α,
β > 0 and
α >
β. The introduced reaching law ensures unique favorable properties of the sliding variable evolution, such as restriction of its rate of change, finite time convergence to the vicinity of zero, and satisfaction of the switching condition stated by Gao et al. [
27]. The motion generated with the reaching law (11) may be divided into two modes. When the sliding variable
sr(
k) is outside the interval [−
α,
α], so for |
sr(
k)| >
α, the first formula of the reaching law acts. In that case, from (11) we get
Considering that
α +
β is always positive and
α >
β then
For |
sr(
k)| >
α +
β, the sign of the reference sliding variable remains unchanged, i.e., sgn[
sr(
k + 1)] = sgn[
sr(
k)]. On the other hand, for
α +
β > |
sr(
k)| >
α, the sign of the reference sliding variable changes for the first time, so sgn[
sr(
k + 1)] = −sgn[
sr(
k)]. Moreover, from (12) we conclude that the rate of change of
sr(
k) is constant and expressed with
Therefore, in the whole reaching phase, when |sr(k)| > α, the reference sliding variable is in each step reduced by a constant value. As follows, according to(9), the sliding surface shifts by a constant distance in each step.
We define
k0 as the first time instant when the reference sliding variable enters the interval (−
α,
α). Therefore
And from (11)
Afterward, for any
k ≥
k0 the second formula of the reaching law acts. In that case, (11) takes the form
Therefore, not later than between steps
k0 and
k0 + 1, the sign of
sr(
k) changes for the first time, i.e., sgn[
sr(
k0 + 1)] = −sgn[
s0]. Moreover, the shift of the sliding surface between steps
k0 and
k0 + 1 satisfies
As
β is a positive constant, from (17) we conclude that for any
k ≥
k0 the change of sign occurs, so
and
which fulfills the requirements for the existence of the quasi-sliding mode. Moreover, for
k ≥
k0 + 1, the rate of change becomes
Taking into account (14), (18), and (21), one notices that the rate of change of sr(k) is limited by α + β in the whole process. It may also be concluded that in the reaching phase, the pace of convergence is mainly regulated by the value of α as α > β. In the sliding phase, the reference sliding variable is, in each step, switched between ±β, so β should be as small as possible to minimize the chattering problem and the quasi-sliding mode bandwidth.
3.2. Control Law
The proposed reference sliding variable profile begins with the initial condition of the original system,
sr(0) =
s(0), monotonically converges to zero and crosses it in finite time. Furthermore, after the sign of
sr(
k) changes for the first time, it will change again in each successive time instant and its value will never exceed
β. One may conclude that the reference profile satisfies all the quasi-sliding mode requirements stated by Gao et al. [
27]. Consequently, the task of the sliding mode controller for the original system (2) becomes to steer the system’s state to the reference position in each step, or in other words to drive the system’s representative point onto the sliding plane (9). The required control law may be obtained by assigning
Substituting (2) into (22) and solving for
u(
k) we get
while the unknown disturbance term is neglected. With the control signal (23), the closed-loop system (2) becomes
As follows, from (9) and (24), the sliding variable of the closed-loop system equals
Let the term
cd(
k) be denoted with
D(
k). Considering the boundedness of the disturbance term |
f(
k)| ≤
fmax, we get
From (26) it may be concluded that the designed sliding mode controller drives the system’s representative point to the reference position with an accuracy of the disturbance in each step. Nevertheless, with the appropriate selection of the reference profile, the properties of the quasi-sliding mode may be preserved, which is shown in the following theorem.
Theorem 1. If β > Dmax then, for any k ≥ k0, sgn[s(k + 1)] = −sgn[s(k)], so the change of sign of s(k) = ce(k) occurs and the absolute value of the sliding variable does not exceed β + Dmax.
Proof of Theorem 1. It has already been shown that sr(k) monotonically converges to zero and crosses it in finite time. Moreover, in the reaching phase (14) holds, so the reference sliding variable decreases in each step by α + β. Consequently, if β > Dmax, from (26), we conclude that s(k) is monotonically decreasing for any k < k0.
The reference
sr(
k) changes its sign for the first time in finite time, not later than in step
k0 + 1. Therefore, sgn[
sr(
k0 + 1)] = −sgn(
s0) and for any
k ≥
k0, (19) and (20) hold. Consequently, from (25), for any
k ≥
k0 we get:
The sign of
D(
k) is unknown, however, taking into account (26), we conclude that when
β >
Dmax, then sgn[
s(
k + 1)] = sgn[
sr(
k + 1)] = −sgn[
s(
k)]. Moreover, considering (20), we notice that for any
k ≥
k0 + 1,
which concludes the proof. □
We have shown that, for any
k ≥
k0 + 1, the switching type quasi-sliding mode exists. The sliding variable is driven to the interval of the width 2(
β +
Dmax). As
β determines the quasi-sliding mode bandwidth and is only limited by
β >
Dmax, it shall be chosen as the smallest admissible value, which results in the bandwidth of
O(
T) order. We must also acknowledge that the smaller
β becomes, the smaller the necessary control effort. As the shift of the sliding surface is strictly limited in each step, the proposed control method guarantees a restriction of the rate of change of the sliding variable. Taking into account the rate of change of the reference sliding variable, for any
k ≤
k0, we get
and for any
k ≥
k0 + 1
The restriction of the rate of change of the sliding variable results in a restriction of the necessary control effort and therefore, it increases the energetical efficiency of the closed-loop system.
3.3. O(T3) Accuracy
In this section, in order to improve the robustness of the system, we supplement our control strategy with a disturbance compensation algorithm. The compensation term is constructed from the current and past state measurements and past control signal value in order to calculate the past disturbance impact
D(
k − 1) and
D(
k − 2), as follows
As the disturbance
f(
t) is a continuous function of time, the rate of change of the discretized disturbance vector
d(
k) is limited, as noticed in (7) and (8). Therefore, its impact on the sliding variable is also limited, which, using (8), may be described with
The obtained previous disturbance values may be used as an estimate for the current disturbance
D(
k) in the revised reaching law
According to (8), the error of such estimation remains of
O(
T3) order. With the application of (33) to system (2), we get the control signal
The estimation of the current disturbance value is a dynamic process. Therefore, the implementation of the disturbance compensator introduces an additional state variable to the analyzed system. With the additional dynamics, the problem may be described with the following set of equations
From the implementation of (35) in (36) one gets
Therefore, the closed-loop system is described with the following state equation
where the state matrix is
Let us denote the eigenvalues of the system (38) with
zi, where
i = 1, …,
n + 1. With this notation the characteristic polynomial of the closed-loop system becomes
To determine the eigenvalues
zi of the state matrix (39), the following similarity transformation may be utilized
where
Ac = [
1n×n −
b(
cb)
−1c]
A and
M is an invertible matrix such that
and
From the multiplication of (41) we get
From (44) it follows that
Equation (45) indicates that the placement of eigenvalues z1, …, zn depends on the choice of the sliding plane. The dynamics of the disturbance compensator introduced one extra eigenvalue zn+1 located at the origin of the complex plane. Therefore, the compensator’s dynamics are stable and the whole system’s stability depends on the choice of vector c only. However, this modified control scheme enforces a modification of the chosen control parameters in order to provide minimal width of the quasi-sliding mode band and consequently reduce chattering effects.
Theorem 2. If β > ∆D then, for any k ≥ k0and the sliding variable satisfies Proof of Theorem 2. In
Section 1 of this chapter, we have shown that
the reference sliding variable sr(k) changes its sign for the first time not later than in k0 + 1,
for any k ≥ k0 sgn[sr(k + 1)] = −sgn[sr(k)],
for any k ≥ k0 + 1, the reference satisfies: |sr(k)| = β.
Therefore, the reference sliding variable satisfies all the properties of the quasi-sliding mode as specified in [
27]. In the reaching phase, for any
k <
k0, the reference sliding variable approaches zero by
α +
β in each step. Considering (32), when
α +
β > ∆
D, then (33) leads to the conclusion that
s(
k) is decreasing by minimum
α +
β − ∆
D for any
k <
k0. Considering the properties of
sr(
k), for any
k ≥
k0 from (33) we get
Taking into account (32), we conclude that if
β > ∆
D, then sgn[
s(
k + 1)] = sgn[
sr(
k + 1)] and
As (19) holds for any k ≥ k0, then (46) is satisfied as well, which completes the proof. □
We have shown that the enhanced version of the controller ensures the existence of the quasi-sliding mode with the ultimate bandwidth of 2(
β + ∆
D). We notice that the parameter
β is directly responsible for the width of the quasi-sliding mode band and it is only constrained by
β > ∆
D. Therefore, knowing the bounds of the disturbance, one may always select
β as the smallest admissible value. As follows,
β becomes of
O(
T3) order. Considering that the second-order disturbance difference is of
O(
T3) order, as defined in (8), the width of the whole quasi-sliding mode band becomes of
O(
T3) order and its minimum values is 4∆
D. This clearly shows that with the reduction of the discretization period, the control error resulting from chattering effects is significantly reduced, which leads to limitation of the control effort and energy consumption in the sliding phase. In this control algorithm, the rate of change of the sliding variable, for any
k ≤
k0, is restricted by
and, for any
k ≥
k0 + 1 by
Restraining the rate of change of the sliding variable directly results in the reduction of the control magnitudes in the reaching phase, which leads to better energetical efficiency of the closed-loop system.
Once we have shown that the application of the proposed disturbance compensator ensures a reduction of the quasi-sliding mode bandwidth to order
O(
T3), we may now proceed to estimate the maximum state variable error values. First, we will concentrate on the stage when the sliding surface approaches the desired position. It will be shown that the state variables deviations from their reference values in that stage are maximum of
O(
T3) order. Next, we will derive the maximum error values for all the system’s state variables, while the sliding surface oscillates around the steady-state. As the choice of the sliding plane is arbitrary, in the following calculations it is assumed that
which implies that all the poles are located at the origin of the complex plane.
Theorem 3. For systems controlled according to (34), in the reaching phasewhere ∆xi is the i-th state variable deviation from the reference value, xr(k) represents the reference state vector and pi is a row vector such thatMoreover, if vector c is selected so that (52) holds, then the state error in the quasi-sliding mode for any j ≥ k0 + n, is limited as followswhere ei denotes the i-th state variable error ei(j) = xdi − xi(j). Proof of Theorem 3. To begin the proof, let us introduce the following notation for d(k)
where
g(
k) is the sum of matched disturbance effects from (6) and the mismatched disturbance term is of order
O(
T3). Next, using (32) and (33), we describe the absolute sliding variable deviation from the reference with
Following (9) and (10), we transform (57) to
Considering that vector
c is of
O(1) order, one may conclude that the
i-th state variable deviation from the reference is
which confirms (53). Next, we proceed to calculate
e(
k + 1) from (2) with the control input (34), which yields
Considering that
Ac = [
1n × n −
b(
cb)
−1c]
A, we transform (60) to
Following the same pattern, we write the state error in step
k + 2
In general, for any
j ≥ 0, the state error vector
e(
j) is
One may notice that when (52) holds, then the state matrix
Ac is nilpotent, so
Acj− k = 0 for any
j −
k ≥
n. Therefore, for any
j ≥
k +
n, the sums in (63) may be shortened to
n elements
Moreover, from (56), we may derive that
Taking into account (56), (65), and the reaching law (33), we notice that
With the application of (66), we transform (64) to the following form
Next, we multiply both sides of (67) with the row vector
pi obtaining the
i-th variable’s error
Equation (68) shows that in each step
j, where
j ≥
k +
n and
k ≥ 0, all state errors may be calculated using the
n previous sliding variable values,
s(
j −
n + 1),
s(
j −
n + 2), …,
s(
j). Considering that for any
k ≥
k0 + 1 the value of the sliding variable is bounded by (49), we conclude that the state errors are bounded as well. From (49) and (68), for any
j ≥
k0 +
n, we obtain
which ends the proof. □
Bearing in mind (8), we remember that ∆D is of O(T3) order. Moreover, it has already been shown that the width of the quasi-sliding mode band directly depends on the choice of parameter β. Therefore, β must be selected as a function of T3 to provide the best robustness and lowest energy costs. Consequently, from (69), one may notice that the maximum state error in the sliding phase is of order O(T3) as well.
4. Simulation Example
To conclude this study, a simulation example will be analyzed. We design the sliding mode controller for a simple third-order plant, described with the following state equation
The initial conditions are
x(0) = [10 0 0]
T and the desired state is
xd = [0 0 0]
T. In this example, the disturbance
f(
t) is a sinusoidal function of time with a frequency of 50 Hz and an amplitude of 1. We aim to control the system digitally. Therefore, we discretize the plant with the discretization period
T = 1 ms, obtaining
where
The eigenvalues of matrix
A in (72) are
z1 =
z2 =
z3 = 1, so the system is on the verge of stability. In the first step of the control design, we select the sliding plane according to (52). Therefore, we set
c = [1 0.001 1/3 × 10
−6], which places all the poles at the origin of the complex plane and therefore, ensures the stability of the closed-loop system. The disturbance
f(
t) acts through the control channel and changes between ±1, so the maximum disturbance impact is
Dmax = 1 × 10
−9 and ∆
D = 0.98 × 10
−10. Therefore, according to theorem 2, we choose the control parameter
β = 10
−10. On the other hand, the parameter
α is chosen in order to ensure a control limitation of 3 × 10
6, so
α = 3 × 10
−3. To present the benefits of our control method we compare it with the seminal strategy of Gao et al. [
27], where the control signal was obtained with the reaching law
For the original strategy we set
qT = 3 × 10
−4 and
ε = 6.67 × 10
−6, to ensure fulfillment of the same control constraint. The results of our simulations are presented in
Figure 1,
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9,
Figure 10,
Figure 11,
Figure 12 and
Figure 13.
Figure 2 and
Figure 3 present the evolution of the discrete-time sliding variable
s(
k) =
ce(
k). As visible, both sliding variables converge to zero in finite time and the quasi-sliding motion appears. In Gao’s case, the quasi-steady state is reached for
t ≥ 20.367 s, while with our control strategy the system oscillates around the equilibrium point for any
t ≥ 3.334 s. Therefore, our strategy ensured faster convergence to the sliding manifold without an additional control effort. Gao’s strategy resulted in the ultimate bandwidth of 2 × 6.671 × 10
−6. Applying the trajectory following the strategy with a compensator, we achieved a reduction of this band to 2 × 1.98 × 10
−10. Moreover,
Figure 4 and
Figure 5 clearly show that apart from reducing the ultimate bandwidth and shortening the reaching phase, our new control method satisfies the same control signal limitation as the original strategy and significantly reduces the control effort in the quasi-steady state. Therefore, the closed-loop system with our new control strategy is more energy-saving, which constitutes an important factor in the control design, especially in the case of power systems or electric motor drives.
Finally,
Figure 6,
Figure 7,
Figure 8,
Figure 9,
Figure 10,
Figure 11,
Figure 12,
Figure 13 and
Figure 14 present the evolution of the state variables
x1(
k),
x2(
k), and
x3(
k) for both control methods. It is clear that the errors of all state variables have been reduced as well.
From the presented graphs we may see the maximum state variables error values. For Gao’s method, the errors are:
Our control method reduced the maximum state error values to:
which is consistent with Theorem 3. This shows that the introduced control scheme ensures a significant improvement of the robustness of the system along with a reduction of the control effort. Lastly, to precisely show the improvement we define a performance quality index for the discrete-time system as:
where
i = 1, 2, 3. The obtained index values are presented in
Table 1.
5. Conclusions
This paper proposed a new discrete-time sliding mode control method based on the time-varying sliding surface concept. The control algorithm assumes the application of a time-varying sliding surface, whose shift is obtained according to the reference sliding variable profile. The pre-generated reference profile manifests all the properties required from the actual system, so it converges to zero in finite time and exhibits the zigzagging motion around the equilibrium point. Such a reference profile is obtained with a new-reaching law. Next, the sliding mode control law for the actual system is developed so that it drives the representative point to the desired position with accuracy to the single-step disturbance. It is proved that the proposed control method ensures finite time convergence of the representative point of the system to an a priori known band around the designed sliding plane and the width of this band is of
O(
T) order. Moreover, it provides a restriction of the rate of change of the sliding variable, which results in a reduction of the control effort. As the control magnitudes are limited and chattering effects are reduced, the system becomes more efficient and energy-saving. This achievement may be especially important in practical applications, for example in the control of power electronics systems or electrical motor drives, where reduction of chattering is one of the main issues. Lastly, we supplement our strategy with a disturbance compensator. Its task is to calculate the disturbance value from two previous steps, based on the previous and current state measurement, and use it to compensate for the disturbance in the current step. As the second-order difference of the discretized disturbance is of
O(
T3) order, the application of such a compensator ensures that both the sliding variable and all the state variables belong to the band of
O(
T3) order. Therefore, the robustness of the system is significantly improved and the energy consumption is reduced. In short, the originality of this paper is manifested in a new efficient reaching law combined with a time-varying sliding surface approach and a disturbance compensator, which results in the restriction of the quasi-sliding mode bandwidth as well as the maximum deviations of all state variables from their demand values to
O(
T3) order. As follows, chattering and the control input are minimized, so the system becomes more efficient. Finally, the paper is concluded with a simulation example, which shows that the proposed control strategy not only provides the reduction of the quasi-sliding mode bandwidth but also ensures faster convergence than the conventional method of Gao [
27], without increasing the control input.