Next Article in Journal
Pore-Scale Modeling of Air–Water Two Phase Flow and Oxygen Transport in Gas Diffusion Layer of Proton Exchange Membrane Fuel Cell
Next Article in Special Issue
Research on Compound Sliding Mode Control of a Permanent Magnet Synchronous Motor in Electromechanical Actuators
Previous Article in Journal
Degradation-Conscious Equivalent Consumption Minimization Strategy for a Fuel Cell Hybrid System
Previous Article in Special Issue
Zero-Width Quasi-Sliding Mode Band in the Presence of Non-Matched Uncertainties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Time-Varying Sliding Surface for Switching Type Quasi-Sliding Mode Control

by
Katarzyna Adamiak
* and
Andrzej Bartoszewicz
Institute of Automatic Control, Łódź University of Technology, 18/22 Bohdana Stefanowskiego St., 90-924 Łódź, Poland
*
Author to whom correspondence should be addressed.
Energies 2021, 14(13), 3811; https://doi.org/10.3390/en14133811
Submission received: 27 May 2021 / Revised: 12 June 2021 / Accepted: 17 June 2021 / Published: 24 June 2021
(This article belongs to the Special Issue Sliding Mode Control in Electromechanical Systems)

Abstract

:
This study considers the problem of energetical efficiency in switching type sliding mode control of discrete-time systems. The aim of this work is to reduce the quasi-sliding mode band-width and, as follows, the necessary control input, through an application of a new type of time-varying sliding hyperplane in quasi-sliding mode control of sampled time systems. Although time-varying sliding hyperplanes are well known to provide insensitivity to matched external disturbances and uncertainties of the model in the whole range of motion for continuous-time systems, their application in the discrete-time case has never been studied in detail. Therefore, this paper proposes a sliding surface, which crosses the system’s representative point at the initial step and then shifts in the state space according to the pre-generated demand profile of the sliding variable. Next, a controller for a real perturbed plant is designed so that it drives the system’s representative point to its reference position on the sliding plane in each step. Therefore, the impact of external disturbances on the system’s trajectory is minimized, which leads to a reduction of the necessary control effort. Moreover, thanks to a new reaching law applied in the reference profile generator, the sliding surface shift in each step is strictly limited and a switching type of motion occurs. Finally, under the assumption of boundedness and smoothness of continuous-time disturbance, a compensation scheme is added. It is proved that this control strategy reduces the quasi-sliding mode band-width from O(T) to O(T3) order from the very beginning of the regulation process. Moreover, it is shown that the maximum state variable errors become of O(T3) order as well. These achievements directly reduce the energy consumption in the closed-loop system, which is nowadays one of the crucial factors in control engineering.

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
x ˙ ( t ) = Φ x ( t ) + Γ u ( t ) + Γ f ( t ) ,
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
x ( k + 1 ) = A x ( k ) + b u ( k ) + d ( k ) ,
where appropriate matrices are obtained according to the following
A = e Φ T ,   b = 0 T e Φ τ d τ Γ ,   d ( k ) = 0 T e Φ τ Γ f [ ( k + 1 ) T τ ] d τ .
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
d ( k ) = 0 T e Φ τ Γ f ( k T ) d τ + 0 T e Φ τ Γ v ( k T ) ( T τ ) d τ + 0 T e Φ τ Γ 1 2 ! w ( μ ) ( T τ ) 2 d τ + ,
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
d ( k ) = b f ( k T ) + 1 2 b v ( k T ) T + 0 T e Φ τ Γ 1 2 ! w ( μ ) ( T τ ) 2 d τ + .
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
d ( k ) = b f ( k T ) + 1 2 b v ( k T ) T + O ( T 3 ) .
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
d ( k ) d ( k 1 ) = O ( T 2 )
and
d ( k ) 2 d ( k 1 ) + d ( k 2 ) = O ( T 3 ) .
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
c e ( k ) s r ( k ) = 0 ,
where c is an 1 × n control vector, chosen so that cb ≠ 0 and to ensure stability of the steady state and e(k) = xdx(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
s r ( k ) = c [ x d x r ( k ) ] ,
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
s r ( k + 1 ) = { s r ( k ) ( α + β ) sgn [ s r ( k ) ] for   | s r ( k ) | > α β sgn [ s r ( k ) ] for   | s r ( k ) | α ,
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
s r ( k + 1 ) = [ | s r ( k ) | ( α + β ) ] sgn [ s r ( k ) ] .
Considering that α + β is always positive and α > β then
| s r ( k + 1 ) | < | s r ( k ) | .
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
| s r ( k ) s r ( k + 1 ) | = α + β .
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
| s r ( k 0 ) | < α
And from (11)
k 0 = floor ( | s 0 | α α + β ) + 1 .
Afterward, for any kk0 the second formula of the reaching law acts. In that case, (11) takes the form
sgn [ s r ( k + 1 ) ] | s r ( k + 1 ) | = β sgn [ s r ( k ) ] .
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
| s r ( k 0 ) s r ( k 0 + 1 ) | α + β .
As β is a positive constant, from (17) we conclude that for any kk0 the change of sign occurs, so
sgn [ s r ( k + 1 ) ] = sgn [ s r ( k ) ]
and
| s r ( k + 1 ) | = β ,
which fulfills the requirements for the existence of the quasi-sliding mode. Moreover, for kk0 + 1, the rate of change becomes
| s r ( k ) s r ( k + 1 ) | = 2 β .
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
s ( k + 1 ) s r ( k + 1 ) = 0 .
Substituting (2) into (22) and solving for u(k) we get
u ( k ) = ( c b ) 1 [ c x d c A x ( k ) s r ( k + 1 ) ] ,
while the unknown disturbance term is neglected. With the control signal (23), the closed-loop system (2) becomes
x ( k + 1 ) = [ 1 b ( c b ) 1 c ] A x ( k ) + b ( c b ) 1 c x d b ( c b ) 1 s r ( k + 1 ) + d ( k ) .
As follows, from (9) and (24), the sliding variable of the closed-loop system equals
s ( k + 1 ) = s r ( k + 1 ) c d ( k ) .
Let the term cd(k) be denoted with D(k). Considering the boundedness of the disturbance term |f(k)| ≤ fmax, we get
| s ( k + 1 ) s r ( k + 1 ) | D max = c b f max .
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 kk0, (19) and (20) hold. Consequently, from (25), for any kk0 we get:
sgn [ s ( k + 1 ) ] | s ( k + 1 ) | = β sgn [ s r ( k + 1 ) ] D ( k ) .
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 kk0 + 1,
| s ( k ) | β + D max .
which concludes the proof. □
We have shown that, for any kk0 + 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 kk0, we get
| s ( k ) s ( k + 1 ) | α + β + 2 D max
and for any kk0 + 1
| s ( k ) s ( k + 1 ) | 2 β + 2 D max .
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
D ( k 1 ) = c [ x ( k ) A x ( k 1 ) b u ( k 1 ) ] , D ( k 2 ) = c [ x ( k 1 ) A x ( k 2 ) b u ( k 2 ) ] .
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
| D ( k ) 2 D ( k 1 ) + D ( k 2 ) | Δ D .
The obtained previous disturbance values may be used as an estimate for the current disturbance D(k) in the revised reaching law
s ( k + 1 ) s r ( k + 1 ) = D ( k ) + 2 D ( k 1 ) D ( k 2 ) .
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
u ( k ) = ( c b ) 1 [ c x d c A x ( k ) s r ( k + 1 ) 2 D ( k 1 ) + D ( k + 2 ) ] .
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
x ( k + 1 ) = A x ( k ) + b u ( k ) + 2 d ( k 1 ) d ( k 2 )
u ( k + 1 ) = ( c b ) 1 [ c x d c A x ( k + 1 ) 2 c d ( k 1 ) + c d ( k 2 ) ] .
From the implementation of (35) in (36) one gets
u ( k + 1 ) = ( c b ) 1 c x d + [ ( c b ) 1 c A 2 ( c b ) 1 c A b ] [ x ( k ) u ( k ) ] ( c b ) 1 c ( A + 1 n × n ) [ 2 d ( k 1 ) d ( k 2 ) ] .
Therefore, the closed-loop system is described with the following state equation
[ x ( k + 1 ) u ( k + 1 ) ] = A [ x ( k ) u ( k ) ] + [ 1 n × n ( c b ) 1 c ( A + 1 n × n ) ] [ 2 d ( k 1 ) d ( k 2 ) ] + [ 0 n × 1 ( c b ) 1 c x d ] ,
where the state matrix is
A = [ A b ( c b ) 1 c A 2 ( c b ) 1 c A b ] .
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
0 = det ( 1 ( n + 1 ) × ( n + 1 ) z i A ) = | 1 n × n z i A b ( c b ) 1 c A 2 z i + ( c b ) 1 c A b | .
To determine the eigenvalues zi of the state matrix (39), the following similarity transformation may be utilized
A = [ A b ( c b ) 1 c A 2 ( c b ) 1 c A b ] = M 1 [ A c b 0 1 × n 0 ] M ,
where Ac = [1n×nb(cb)−1c]A and M is an invertible matrix such that
M = [ 1 n × n 0 n × 1 ( c b ) 1 c A 1 ]
and
M 1 = [ 1 n × n 0 n × 1 ( c b ) 1 c A 1 ] .
From the multiplication of (41) we get
M 1 [ A c b 0 1 × n 0 ] M = [ 1 n × n 0 n × 1 ( c b ) 1 c A 1 ] [ ( 1 n × n b ( c b ) 1 c ) A b 0 1 × n 0 ] [ 1 n × n 0 n × 1 ( c b ) 1 c A 1 ] = [ ( 1 n × n b ( c b ) 1 c ) A b ( c b ) 1 c A ( 1 n × n b ( c b ) 1 c ) A ( c b ) 1 c A b ] [ 1 n × n 0 n × 1 ( c b ) 1 c A 1 ] = [ A b ( c b ) 1 c A 2 ( c b ) 1 c A b ] = A .
From (44) it follows that
z 1 ,   ... ,   n = eig { [ 1 n × n b ( c b ) 1 c ] A }     z n + 1 = 0 .
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 ≥ k0
sgn [ s ( k + 1 ) ] = sgn [ s ( k ) ]  
and the sliding variable satisfies
| s ( k ) | β + Δ D .
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 kk0 sgn[sr(k + 1)] = −sgn[sr(k)],
  • for any kk0 + 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 kk0 from (33) we get
sgn [ s ( k + 1 ) ] | s ( k + 1 ) | = β sgn [ s r ( k + 1 ) ] D ( k ) + 2 D ( k 1 ) D ( k 2 ) .
Taking into account (32), we conclude that if β >D, then sgn[s(k + 1)] = sgn[sr(k + 1)] and
| s ( k + 1 ) | β + Δ D .
As (19) holds for any kk0, 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 kk0, is restricted by
| s ( k ) s ( k + 1 ) | α + β + 2 Δ D
and, for any kk0 + 1 by
| s ( k ) s ( k + 1 ) | 2 β + 2 Δ D .
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
z 1 ,   ,   n = eig { [ 1 n × n b ( c b ) 1 c ] A } = 0 ,
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 phase
| Δ x i ( k ) | = | p i [ x ( k ) x r ( k ) ] | = O ( T 3 ) ,
where ∆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 that
p i = [ 0 0 i 1 1 0 0 n i ] .
Moreover, 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 follows
| e i ( j ) | p i | x d l = 0 n 1 A c l b ( c b ) 1 c x d | + ( β + Δ D ) l = 0 n 1 | p i A c l b ( c b ) 1 | ,
where 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)
d ( k ) = b g ( k ) + O ( T 3 ) ,
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
| s ( k + 1 ) s r ( k + 1 ) | = | D ( k ) 2 D ( k 1 ) + D ( k 2 ) | Δ D .
Following (9) and (10), we transform (57) to
| c [ x r ( k + 1 ) x ( k + 1 ) ] | = | c [ d ( k ) 2 d ( k 1 ) + d ( k 2 ) ] | Δ D = O ( T 3 ) .
Considering that vector c is of O(1) order, one may conclude that the i-th state variable deviation from the reference is
| p i [ x ( k ) x r ( k ) ] | = O ( T 3 ) ,
which confirms (53). Next, we proceed to calculate e(k + 1) from (2) with the control input (34), which yields
e ( k + 1 ) = x d A x ( k ) b ( c b ) 1 c x d + b ( c b ) 1 c A x ( k ) + b ( c b ) 1 s r ( k + 1 ) b ( c b ) 1 D ( k 2 ) + 2 b ( c b ) 1 D ( k 1 ) d ( k ) .
Considering that Ac = [1n × nb(cb)−1c]A, we transform (60) to
e ( k + 1 ) = x d A c 0 b ( c b ) 1 x d A c 1 x ( k ) + A c 0 b ( c b ) 1 s r ( k + 1 ) + A c 0 b ( c b ) 1 [ 2 D ( k 1 ) D ( k 2 ) ] A c 0 d ( k ) .
Following the same pattern, we write the state error in step k + 2
e ( k + 2 ) = x d A c 0 b ( c b ) 1 c x d A c 1 b ( c b ) 1 c x d A c 2 x ( k ) + A c 1 b ( c b ) 1 s r ( k + 1 ) + + A c 1 b ( c b ) 1 [ 2 D ( k 1 ) D ( k 2 ) ] A c 1 d ( k ) + A c 0 b ( c b ) 1 s r ( k + 2 ) + A c 0 b ( c b ) 1 [ 2 D ( k ) D ( k 1 ) ] A c 0 d ( k + 1 ) .
In general, for any j ≥ 0, the state error vector e(j) is
e ( j ) = x d l = 0 j k 1 A c l b ( c b ) 1 c x d A c j k x ( k ) + + l = 0 j k 1 A c l { b ( c b ) 1 s r ( j l ) + b ( c b ) 1 [ 2 D ( j l 2 ) D ( j l 3 ) ] d ( j l 1 ) } .
One may notice that when (52) holds, then the state matrix Ac is nilpotent, so Acj− k = 0 for any jkn. Therefore, for any jk + n, the sums in (63) may be shortened to n elements
e ( j ) = x d l = 0 n 1 A c l b ( c b ) 1 c x d + l = 0 n 1 A c l { b ( c b ) 1 s r ( j l ) + b ( c b ) 1 [ 2 D ( j l 2 ) D ( j l 3 ) ] d ( j l 1 ) } .
Moreover, from (56), we may derive that
b ( c b ) 1 D ( k ) = b ( c b ) 1 c [ b g ( k ) + O ( T 3 ) ] = b g ( k ) + O ( T 3 ) .
Taking into account (56), (65), and the reaching law (33), we notice that
b ( c b ) 1 s r ( j l ) + b ( c b ) 1 [ 2 D ( j l 2 ) D ( j l 3 ) ] d ( j l 1 ) = b ( c b ) 1 s ( j l ) + O ( T 3 ) .
With the application of (66), we transform (64) to the following form
e ( j ) = x d l = 0 n 1 A c l b ( c b ) 1 c x d + l = 0 n 1 A c l b ( c b ) 1 s ( j l ) + O ( T 3 ) .
Next, we multiply both sides of (67) with the row vector pi obtaining the i-th variable’s error
e i ( j ) = x d i l = 0 n 1 p i A c l b ( c b ) 1 c x d l = 0 n 1 p i A c l b ( c b ) 1 s ( j l ) + O ( T 3 ) .
Equation (68) shows that in each step j, where jk + n and k ≥ 0, all state errors may be calculated using the n previous sliding variable values, s(jn + 1), s(jn + 2), …, s(j). Considering that for any kk0 + 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 jk0 + n, we obtain
| e i ( j ) | | x d i l = 0 n 1 p i A c l b ( c b ) 1 c x d | + ( β + Δ D ) l = 0 n 1 | p i A c l b ( c b ) 1 | + O ( T 3 ) ,
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
x ˙ ( t ) = [ 0 1 0 0 0 1 0 0 0 ] x ( t ) + [ 0 0 1 ] [ u ( t ) + f ( t ) ] ,
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
x ( k + 1 ) = A x ( k ) + b u ( k ) + d ( k ) ,
where
A = [ 1 0.001 5 10 7 0 1 0.001 0 0 1 ] ,   b = [ 1 6 10 9 5 10 7 0.001 ] .
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 × 106, 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
s ( k + 1 ) = ( 1 q T ) s ( k ) ε sgn [ s ( k ) ] D ( k ) .
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:
| e 1 ( k ) | 1.9 × 10 6 ,   | e 2 ( k ) | 1.3 × 10 6 ,   | e 3 ( k ) | 13.5 .
Our control method reduced the maximum state error values to:
| e 1 ( k ) | 1.4 × 10 10 ,   | e 2 ( k ) | 3.2 × 10 8 ,   | e 3 ( k ) | 4.1 × 10 4 ,
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:
P Q I = k = 0 50,000 | e i ( k ) | ,
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.

Author Contributions

Conceptualization, A.B. and K.A.; methodology, A.B. and K.A.; validation, A.B. and K.A.; formal analysis, K.A.; investigation, K.A; writing—original draft preparation, K.A.; writing—review and editing, A.B.; visualization, K.A.; supervision, A.B.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Emelyanov, S.V. Variable Structure Control Systems; Nauka: Moscow, Russia, 1967. [Google Scholar]
  2. Itkis, U. Control Systems of Variable Structures; Wiley: New York, NY, USA, 1976. [Google Scholar]
  3. Utkin, V.I. Variable structure systems: Present and future. Autom. Remote Control 1984, 44, 1105–1220. [Google Scholar]
  4. Edwards, C. A practical method for the design of sliding mode controllers using linear matrix inequalities. Automatica 2004, 40, 1761–1769. [Google Scholar] [CrossRef]
  5. DeCarlo, R.A.; Żak, S.; Matthews, G. Variable structure control of nonlinear multivariable systems: A tutorial. Proc. IEEE 1988, 76, 212–232. [Google Scholar] [CrossRef] [Green Version]
  6. Slotine, J.E. Sliding controller design for nonlinear system. Int. J. Control 1989, 40, 421–434. [Google Scholar] [CrossRef]
  7. Hung, J.Y.; Gao, W.B.; Hung, J.C. Variable structure control: A survey. IEEE Trans. Ind. Electron. 1993, 40, 2–22. [Google Scholar] [CrossRef] [Green Version]
  8. Zinober, A.S.I. Variable Structure and Lyapunov Control; Springer: London, UK, 1994. [Google Scholar]
  9. Edwards, C.; Colet, E.F.; Fridman, L. Advances in variable structure and sliding mode control. In Lecture Notes in Control and Information Sciences; Springer: Berlin, Germany, 2006; Volume 334. [Google Scholar]
  10. Edwards, C.; Spurgeon, S. Sliding Mode Control: Theory and Applications; Taylor & Francis: London, UK, 1998. [Google Scholar]
  11. Mobayen, S.; Bayat, F.; Lai, C.-C.; Taheri, A.; Fekih, A. Adaptive global sliding mode controller design for perturbed DC-DC buck converters. Energies 2021, 14, 1249. [Google Scholar] [CrossRef]
  12. Oucheriah, S.; Guo, L. PWM-based adaptive sliding-mode control for boost DC–DC converters. IEEE Trans. Ind. Electron. 2013, 60, 3291–3294. [Google Scholar] [CrossRef]
  13. Gao, P.; Zhang, G.; Lv, X. Model-free control using improved smoothing extended state observer and super-twisting nonlinear sliding mode control for PMSM drives. Energies 2021, 14, 922. [Google Scholar] [CrossRef]
  14. Tran, A.-T.; Minh, B.L.N.; Huynh, V.V.; Tran, P.T.; Amaefule, E.N.; Phan, V.-D.; Nguyen, T.M. Load frequency regulator in interconnected power system using second-order sliding mode control combined with state estimator. Energies 2021, 14, 863. [Google Scholar] [CrossRef]
  15. Draženović, B. The invariance conditions in variable structure systems. Automatica 1969, 5, 287–295. [Google Scholar] [CrossRef]
  16. Levant, A. Sliding order and sliding accuracy in sliding mode control. Int. J. Control 1993, 58, 1247–1263. [Google Scholar] [CrossRef]
  17. Choi, S.B.; Cheong, C.C.; Park, D.W. Moving switching surfaces for robust control of second-order variable structure systems. Int. J. Control 1993, 58, 229–245. [Google Scholar] [CrossRef]
  18. Choi, S.B.; Park, D.W.; Jayasuriya, S. A time-varying sliding surface for fast and robust tracking control of second-order uncertain systems. Automatica 1994, 30, 899–904. [Google Scholar] [CrossRef]
  19. Betin, F.; Pinchon, D.; Capolino, G. A time-varying sliding surface for robust position control of a DC motor drive. IEEE Trans. Ind. Electron. 2002, 49, 462–473. [Google Scholar] [CrossRef]
  20. Milosavljević, C. General conditions for the existence of a quasi-sliding mode on the switching hyperplane in discrete variable structure systems. Autom. Remote Control 1985, 46, 307–314. [Google Scholar]
  21. Sira-Ramirez, H. Non-linear discrete variable structure systems in quasi-sliding mode. Int. J. Control 1991, 54, 1171–1187. [Google Scholar] [CrossRef]
  22. Sarpturk, S.Z.; Istefanopulos, Y.; Kaynak, O. On the stability of discrete-time sliding mode control systems. IEEE Trans. Autom. Control 1987, 22, 930–932. [Google Scholar] [CrossRef]
  23. Kotta, U.; Sarpturk, S.Z.; Istefanopulos, Y. Comments on ‘On the stability of discrete-time sliding mode control systems’ [with reply]. IEEE Trans. Autom. Control 1989, 34, 1021–1022. [Google Scholar] [CrossRef]
  24. Utkin, V.I.; Drakunov, S.V. On discrete-time sliding modes. IFAC Proc. Vol. 1989, 22, 273–278. [Google Scholar]
  25. Bartolini, G.; Ferrara, A.; Utkin, V.I. Adaptive sliding mode control in discrete-time systems. Automatica 1995, 31, 769–773. [Google Scholar] [CrossRef]
  26. Furuta, K. Sliding mode control of a discrete system. Syst. Control Lett. 1990, 14, 145–152. [Google Scholar] [CrossRef]
  27. Gao, W.B.; Wang, Y.; Homaifa, A. Discrete-time variable structure control systems. IEEE Trans. Ind. Electron. 1995, 42, 117–122. [Google Scholar]
  28. Veselić, B.; Perunicic-Draženović, B.; Milosavljević, C. Improved discrete-time sliding-mode position control using Euler velocity estimation. IEEE Trans. Ind. Electron. 2010, 57, 3840–3847. [Google Scholar] [CrossRef]
  29. Ma, H.; Wu, J.; Xiong, Z. Discrete-time sliding-mode control with improved quasi-sliding-mode domain. IEEE Trans. Ind. Electron. 2016, 63, 6292–6304. [Google Scholar] [CrossRef]
  30. Chakrabarty, S.; Bandyopadhyay, B. A generalized reaching law for discrete-time sliding mode control. Automatica 2015, 52, 83–86. [Google Scholar] [CrossRef]
  31. Wu-Chung, S.; Drakunov, S.V.; Ozguner, U. An O(T2) boundary layer in sliding mode for sampled-data systems. IEEE Trans. Autom. Control 2000, 45, 482–485. [Google Scholar] [CrossRef]
  32. Rezig, S.; Ezzeddine, W.; Turki, S.; Rezg, N. Mathematical Model for Production Plan Optimization—A Case Study of Discrete Event Systems. Mathematics 2020, 8, 955. [Google Scholar] [CrossRef]
  33. Golo, G.; Milosavljević, C. Robust discrete-time chattering free sliding mode control. Syst. Control Lett. 2000, 41, 19–28. [Google Scholar] [CrossRef]
  34. Du, H.; Yu, X.; Chen, M.Z.Q.; Li, S. Chattering-free discrete-time sliding mode control. Automatica 2016, 68, 87–91. [Google Scholar] [CrossRef]
  35. Shah, D.; Mehta, A. Discrete-Time Sliding Mode Control for Networked Control System. In Studies in Systems, Decision and Control; Springer: Singapore, 2018. [Google Scholar]
  36. Adamiak, K. Reference Sliding Variable Based Chattering-Free Quasi-Sliding Mode Control. IEEE Access 2020, 8, 133086–133094. [Google Scholar] [CrossRef]
  37. Abidi, K.; Xu, J.; Yu, X. On the discrete-time integral sliding-mode control. IEEE Trans. Autom. Control 2007, 52, 709–715. [Google Scholar] [CrossRef] [Green Version]
  38. Qu, S.; Xia, X.; Zhang, J. Dynamics of discrete-time sliding-mode-control uncertain systems with a disturbance compensator. IEEE. Trans. Ind. Electron. 2013, 61, 3502–3510. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Structure of the closed-loop system.
Figure 1. Structure of the closed-loop system.
Energies 14 03811 g001
Figure 2. Comparison of the discrete-time sliding variables.
Figure 2. Comparison of the discrete-time sliding variables.
Energies 14 03811 g002
Figure 3. Comparison of the discrete-time sliding variables in the sliding phase.
Figure 3. Comparison of the discrete-time sliding variables in the sliding phase.
Energies 14 03811 g003
Figure 4. Evolution of the control signal u(k).
Figure 4. Evolution of the control signal u(k).
Energies 14 03811 g004
Figure 5. Evolution of the control signal in the sliding phase.
Figure 5. Evolution of the control signal in the sliding phase.
Energies 14 03811 g005
Figure 6. x1(k) for both control strategies.
Figure 6. x1(k) for both control strategies.
Energies 14 03811 g006
Figure 7. x1(k) for both control strategies in the sliding phase.
Figure 7. x1(k) for both control strategies in the sliding phase.
Energies 14 03811 g007
Figure 8. x1(k) for both control strategies in the sliding phase (zoomed).
Figure 8. x1(k) for both control strategies in the sliding phase (zoomed).
Energies 14 03811 g008
Figure 9. x2(k) for both control strategies.
Figure 9. x2(k) for both control strategies.
Energies 14 03811 g009
Figure 10. x2(k) for both control strategies in the sliding phase.
Figure 10. x2(k) for both control strategies in the sliding phase.
Energies 14 03811 g010
Figure 11. x2(k) for both control strategies in the sliding phase (zoomed).
Figure 11. x2(k) for both control strategies in the sliding phase (zoomed).
Energies 14 03811 g011
Figure 12. x3(k) for both control strategies.
Figure 12. x3(k) for both control strategies.
Energies 14 03811 g012
Figure 13. x3(k) for both control strategies in the sliding phase.
Figure 13. x3(k) for both control strategies in the sliding phase.
Energies 14 03811 g013
Figure 14. x3(k) for both control strategies in the sliding phase (zoomed).
Figure 14. x3(k) for both control strategies in the sliding phase (zoomed).
Energies 14 03811 g014
Table 1. Performance Quality Index.
Table 1. Performance Quality Index.
Control StrategyPQI
x1x2x3
Gao’s control 32.8910.17401.36
New control method16.68106.02
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adamiak, K.; Bartoszewicz, A. New Time-Varying Sliding Surface for Switching Type Quasi-Sliding Mode Control. Energies 2021, 14, 3811. https://doi.org/10.3390/en14133811

AMA Style

Adamiak K, Bartoszewicz A. New Time-Varying Sliding Surface for Switching Type Quasi-Sliding Mode Control. Energies. 2021; 14(13):3811. https://doi.org/10.3390/en14133811

Chicago/Turabian Style

Adamiak, Katarzyna, and Andrzej Bartoszewicz. 2021. "New Time-Varying Sliding Surface for Switching Type Quasi-Sliding Mode Control" Energies 14, no. 13: 3811. https://doi.org/10.3390/en14133811

APA Style

Adamiak, K., & Bartoszewicz, A. (2021). New Time-Varying Sliding Surface for Switching Type Quasi-Sliding Mode Control. Energies, 14(13), 3811. https://doi.org/10.3390/en14133811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop