Next Article in Journal
Study of Flexibility Transformation in Thermal Power Enterprises under Multi-Factor Drivers: Application of Complex-Network Evolutionary Game Theory
Previous Article in Journal
NIGA: A Novel Method for Investigating the Attacker–Defender Model within Critical Infrastructure Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Uniformly Continuous Generalized Sliding Mode Control

by
Aldo Jonathan Muñoz-Vázquez
1,*,†,‡ and
Guillermo Fernández-Anaya
2,†
1
Department of Multidisciplinary Engineering, Texas A&M University, McAllen, TX 78504, USA
2
Department of Physics and Mathematics, Universidad Iberoamericana, Mexico City 01219, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: Higher Education Center at McAllen, 6200 Tres Lagos Blv, McAllen, TX 78504, USA.
Mathematics 2024, 12(16), 2536; https://doi.org/10.3390/math12162536
Submission received: 29 July 2024 / Revised: 11 August 2024 / Accepted: 14 August 2024 / Published: 16 August 2024
(This article belongs to the Section Dynamical Systems)

Abstract

:
This paper explores a general class of singular kernels with the objective of designing new families of uniformly continuous sliding mode controllers. The proposed controller results from filtering a discontinuous switching function by means of a Sonine integral, producing a uniformly continuous control signal, preserving finite-time sliding motion and robustness against continuous but unknown and not necessarily integer-order differentiable disturbances. The principle of dynamic memory resetting is considered to demonstrate finite-time stability. A set of sufficient conditions to design singular kernels, preserving the above characteristics, is presented, and several examples are exposed to propose new families of continuous sliding mode approaches. Simulation results are studied to illustrate the feasibility of some of the proposed schemes.

1. Introduction

Sliding mode techniques are highly regarded by the interested community, as they are reliable, simple and robust control methods [1]. Nevertheless, classical structures induce harmful-chattering (high-frequency components in the control signal) and self-sustained oscillations in the system response. The harmful-chattering problem, due to a discontinuous controller implementation, can be alleviated by using continuous sliding mode schemes that result from integrating the signum function [2,3]. However, in the integer-order case, continuous sliding mode controllers assume the integer-order differentiability of the coupled disturbance, which is very restrictive in several application scenarios. During the last decade, different versions of continuous fractional sliding structures were proposed, which enforce robustness against continuous but not necessarily differentiable disturbances [4], nonetheless requiring an estimation for the least upper bound of the order of differentiability of the disturbance, thus obtaining conservative formulations, and leading to responses that are similar to those obtained by means of a discontinuous controller implementation. From a practical point of view, adaptive fractional nonlinear controllers offer a solution to the estimation of the order of differentiability [5], but from a theoretical point of view, the attractive property of finite-time sliding motion is lost.
On the one hand, integer-order sliding mode control spans a wide range of applications, such as the robust tracking of robotic manipulators [6], the control of chemical processes [7], power control in wind turbines [8,9], anti-lock brake systems [10], glucose regulation [11], etc. On the other hand, fractional-order sliding mode control has demonstrated a superior performance for several different cases: for fractional-order sliding phase-based controllers, we can find successful implementations in robotic manipulators [12], vehicle suspension systems [13], and wind turbines [14], among others; nevertheless, there are fewer examples that consider fractional-order reaching phase-based controllers, since their complexity increases [15,16,17]. Additional successful implementations of integer and fractional sliding mode schemes can be found in [18,19,20,21,22,23,24,25]. Inspired by these aforementioned contributions, a generalized reaching phase control paradigm is considered in this paper, where robustness and finite-time convergence rely on the principle of dynamic memory resetting.
This paper studies the plausibility of a general class of singular kernels to design robust and continuous sliding mode controllers. A particular version of these non-singular kernels are those functions used to define fractional-order operators, leading to the so-called fractional sliding mode control. However, more general structures are still unknown, not to mention their properties such as stability, robustness and convergence in finite time. The motivation is clear, as the proposed generalization produces sliding mode controllers able to reject non-differentiable disturbances with a continuous control signal. The given disturbances are more general than Hölder (fractional-order differentiable) or Lipschitz (integer-order differentiable) functions.
The main contribution of this paper can be itemized as follows:
  • It is demonstrated that a general class of singular kernels can be used to define novel families of uniformly continuous sliding mode controllers.
  • The emerging continuous sliding mode controllers are capable of rejecting non-differentiable disturbances, which include Hölder continuous disturbances as particular cases, being more general than Lipschitz or integer-order differentiable disturbances.
  • Several examples of generalization are given for a broader choice of singular kernels.
The implications of this paper would be of particular interest for designing robust controllers for autonomous systems, which are expected to operate with a high degree of accuracy, and in uncertain environments. Potential applications could be found in robotic manipulators, autonomous vehicles, renewable energy conversion systems, among others.
Integro-differential operators with singular kernels were suggested in [26], and later, additional properties for these operators and some generalization strategies were given in the inspiring works [27,28,29]. The stability and stabilization of a class of dynamic systems described by general operators were studied in [30], although several problems are still unresolved, such as the stability of non-smooth systems that are associated to kernels that do not have a Laplace transform. The use of Sonine operators was considered to design a more general class of robust proportional–integral-like controllers [31], inducing additional advantages with respect to conventional integer- and fractional-order schemes. Other applications of Sonine operators are given in [32] to analyze systems with variable-order operators in a consistent way, such that for every variable-order integral there is a well-suited left-inverse operator, or variable-order derivative. This change in perspective is applied in [33] to study variable-order fractional relaxation processes. It is worth mentioning that the search for new Sonine kernels still goes on [34]. Non-singular kernel based operators are also attractive for several applications [35,36], but these are out of the scope of this document since the methods proposed in this paper rely on the singularity of the kernel to induce a finite-time sliding motion.
The rest of this paper is organized as follows: Section 2 gives the required framework on generalized operators. Section 3 states the required assumptions, presents the control proposal and demonstrates finite-time stability. Section 4 presents some examples of uniformly continuous sliding mode controllers. Section 5 presents two cases of simulation. Finally, Section 6 discusses the obtained results, states some unresolved problems and comments on potential future work, and Section 7 discusses the main conclusions.

2. Sonine Operators

The central tool of this paper is the definition of a kernel pair in the Sonine sense, which is of preponderant importance in the study of generalized calculus [27,28,29]. In this paper, the controller definition relies on a Sonine operator, such that its regularity and robustness can be modulated by means of the design of the kernel function.
The following definition is considered:
Definition 1.
Let κ ( t ) and λ ( t ) L l o c 1 [ 0 , ) be two non-negative and strictly monotonically decreasing functions that satisfy the so-called Sonine condition
0 t λ ( t τ ) κ ( τ ) d τ = 1 f o r   a r b i t r a r y t 0 .
Functions κ ( t ) and λ ( t ) constitute a kernel pair.
In this paper, we also consider that if κ ( t ) belongs to a kernel pair, satisfies the following properties [30]:
  • lim t 0 κ ( t ) = ,
  • lim t κ ( t ) = 0 ,
  • lim t 0 0 t κ ( τ ) d τ = 0 ,
  • lim t 0 t κ ( τ ) d τ = .
The following definitions are important for the subsequent analysis and the control proposal: [27,30]:
Definition 2
(Generalized integral). Let f L l o c [ a , ) and κ L l o c 1 [ 0 , ) , with κ ( t ) an absolutely continuous functions on [ t 0 , t ] , for any t 0 > 0 . The function
I a κ ( t ) f ( t ) = a t κ ( t τ ) f ( τ ) d τ ,
for t a , is the generalized integral of f ( t ) with respect to the kernel κ ( t ) .
Definition 3
(Generalized derivative). Let y ( t ) be a sufficiently regular continuous function, such that the expression
D a λ ( t ) y ( t ) = λ ( t a ) [ y ( t ) y ( a ) ] + a t λ ( t τ ) τ [ y ( t ) y ( τ ) ] d τ
is real-valued. Then, D a λ ( t ) y ( t ) is the generalized derivative of y ( t ) with respect to λ ( t ) , and with lower terminal a R .
Definition 4
(Generalized differentiability). Whenever λ ( z ) | y ( t + z ) y ( t ) | is locally bounded, for all t a and z 0 , for some a R , function y ( t ) is called marginally D λ -differentiable (for lower terminal a). In addition, if y ( t ) satisfies
lim z 0 λ ( z ) | y ( t + z ) y ( t ) | = 0 .
Then, y ( t ) is called D λ -differentiable (for lower terminal a).
The following proposition relates the Definitions 3 and 4:
Proposition 1.
Let y ( t ) be a locally bounded and D λ -differentiable function, for some lower terminal a R . Then, D a λ ( t ) y ( t ) is locally bounded.
Proof. 
The term λ ( t a ) [ y ( t ) y ( a ) ] in the right-hand side of (3) is locally bounded since z = t a 0 leads to λ ( z ) [ y ( a + z ) y ( a ) ] , which vanishes at z 0 .
Let us set some t > a . Thus, consider a constant L t > 0 , such that
λ ( t τ ) | y ( t ) y ( τ ) | L t .
Then, for an arbitrarily small ε > 0 , we have
a t ε λ ( t τ ) τ [ y ( t ) y ( τ ) ] d τ a t ε λ ( t τ ) τ | y ( t ) y ( τ ) | d τ L t a t ε λ ( t τ ) τ 1 λ ( t τ ) d τ L t t a ε λ ( z ) λ ( z ) d z = ln λ ( ε ) λ ( t a ) L t ,
which is finite. In addition, consider
t ε t λ ( t τ ) τ [ y ( t ) y ( τ ) ] d τ t ε t λ ( t τ ) τ | y ( t ) y ( τ ) | d τ ,
as well as the sequence of functions
f n ( τ ) = λ ( t τ ) λ ( t [ τ δ / n ] ) δ / n ,
for a sufficiently small δ > 0 . It is clear that f n ( τ ) 0 for all n N , and that f n ( τ ) λ ( t τ ) τ as n . Thus, using the Fatou’s Lemma, we obtain
t ε t λ ( t τ ) τ | y ( t ) y ( τ ) | d τ lim n t ε t λ ( t τ ) λ ( t [ τ δ / n ] ) δ / n | y ( t ) y ( τ ) | d τ = lim n n δ t ε t λ ( t τ ) | y ( t ) y ( τ ) | λ ( t [ τ δ / n ] ) | y ( t ) y ( τ ) | d τ .
Moreover, as λ ( t τ ) > λ ( t [ τ δ / n ] ) for any n N , both λ ( t τ ) and λ ( t [ τ δ / n ] ) are non-negative, and λ ( t τ ) | y ( t ) y ( τ ) | is integrable and goes to zero as τ t , the integral in the right-hand side of the last inequality vanishes as ε 0 . Thus, given η > 0 , we can chose ε > 0 as small as required, such that
t ε t λ ( t τ ) τ | y ( t ) y ( τ ) | d τ < η ,
consequently,
D a λ ( t ) y ( t ) λ ( t a ) | y ( t ) y ( a ) | + ln λ ( ε ) λ ( t a ) L t + η .
The limit case t a is obtained in a similar fashion by choosing t = a + ε , and repeating the above process for t ε t λ ( t τ ) τ | y ( t ) y ( τ ) | d τ , obtaining lim t a D a λ ( t ) y ( t ) = 0 . □
In a more general case when y ( t ) is only marginally D λ -differentiable, it is difficult to show that D a λ ( t ) y ( t ) exists. Nonetheless, in accordance with [30], if y ( t ) can be written as a constant plus the generalized integral with respect to κ ( t ) of some function φ L l o c [ 0 , ) , this is, y ( t ) = y ( a ) + I a κ ( t ) φ ( t ) , it results that φ ( t ) = D a λ ( t ) y ( t ) .
In order to prove that the generalized differentiability notion is weaker than integer-order differentiability, consider the following:
Proposition 2.
Let y ( t ) be a Lipschitz continuous real-valued function. Then y ( t ) is D λ -differentiable.
Proof. 
By virtue of y ( t ) is Lipschitz continuous, there is some constant L > 0 , such that
| y ( t + z ) y ( t ) | L z ,
for any real t and z. Therefore, for t , z > 0 , we have
λ ( z ) | y ( t + z ) y ( t ) | L λ ( z ) z .
Then, as both λ ( z ) and z are continuous functions for z > 0 , the function λ ( z ) z is continuous, and thus locally bounded for any z > 0 .
For the limit case z 0 , consider that λ ( z ) λ ( t ) whenever t z . Whereby, integrating with respect to t from 0 to z, we obtain
λ ( z ) z = λ ( z ) 0 z d t = 0 z λ ( z ) d t 0 z λ ( t ) d t .
Therefore, as the last integral vanishes when z 0 , it follows that the Lipchitz continuous function y ( t ) is D λ -differentiable. □
It is worth remembering that (3) coincides with the generalized Caputo-like derivative   C D a λ ( t ) y ( t ) = I a λ ( t ) y ˙ ( t ) , whenever y ( t ) is Lipschitz continuous.
The following translation properties are also useful to prove that other important properties are preserved independently of the lower terminal [30].
Proposition 3.
Let x ( t ) be a function with well-posed generalized integral I a κ ( t ) x ( t ) and y ( t ) be at least marginally D λ -differentiable function. Then,
I a κ ( t ) x ( t ) | t = t ¯ + a = I 0 κ ( t ¯ ) x ¯ ( t ¯ ) ,
D a λ ( t ) y ( t ) | t = t ¯ + a = D 0 λ ( t ¯ ) y ¯ ( t ¯ ) ,
where x ¯ ( t ¯ ) = x ( t ¯ + a ) and y ¯ ( t ¯ ) = y ( t ¯ + a ) .
An interesting property of the fractional integral in (2) is unveiled as follows:
Proposition 4.
Let f L l o c [ 0 , ) . Then, F ( t ) = I 0 κ ( t ) f ( t ) is continuous over [ 0 , ) .
Proof. 
Let ε > 0 be arbitrarily small. Consider a couple of instants t and t [ 0 , ) , and without loss of generality assume t < t (as the case t = t is trivial and the case t > t is analogous), this is, t = t + δ , for some δ > 0 . Then, we obtain
| F ( t ) F ( t ) | = 0 t κ ( t τ ) f ( τ ) d τ 0 t κ ( t τ ) f ( τ ) d τ = 0 t + δ κ ( t + δ τ ) f ( τ ) d τ 0 t κ ( t τ ) f ( τ ) d τ t t + δ κ ( t + δ τ ) f ( τ ) d τ + 0 t [ κ ( t + δ τ ) κ ( t τ ) ] f ( τ ) d τ ] .
First, if we define γ f = ess sup τ [ 0 , t + δ ) | f ( τ ) | , since κ ( σ ) > 0 for every σ > 0 , we obtain
t t + δ κ ( t + δ τ ) f ( τ ) d τ γ f t t + δ κ ( t + δ τ ) d τ = γ f 0 δ κ ( z ) d z ,
where the change of variable z = t + δ τ and terminals inversion were performed to obtain the last integral. Moreover, as 0 δ κ ( τ ) d τ is absolutely continuous, there is a small enough δ > 0 , such that 0 δ κ ( τ ) d τ < ε 2 γ f .
Secondly, we have
0 t [ κ ( t + δ τ ) κ ( t τ ) ] f ( τ ) d τ ] γ f 0 t [ κ ( t τ ) κ ( t + δ τ ) ] d τ = γ f 0 t κ ( t τ ) d τ 0 t κ ( t + δ τ ) d τ = γ f 0 t κ ( τ ) d τ δ t + δ κ ( τ ) d τ .
In the above expression, we can notice that
δ t + δ κ ( τ ) d τ = 0 t + δ κ ( τ ) d τ 0 δ κ ( τ ) d τ = 0 t κ ( τ ) d τ + t t + δ κ ( τ ) d τ 0 δ κ ( τ ) d τ .
Therefore,
0 t [ κ ( t + δ τ ) κ ( t τ ) ] f ( τ ) d τ ] γ f 0 δ κ ( τ ) d τ γ f t t + δ κ ( τ ) d τ 0 γ f 0 δ κ ( τ ) d τ < ε 2 .
Therefore, | F ( t ) F ( t ) | < ε 2 + ε 2 = ε , and the generalized integral of a locally bounded function is continuous. □
The following two Corollaries expose additional properties:
Corollary 1.
If f L l o c [ 0 , ) is replaced by f L [ 0 , ) , we verify that F ( t ) = I 0 κ ( t ) f ( t ) is uniformly continuous on [ 0 , ) .
Corollary 2.
Let x [ 0 , ) R be a Lebesgue measurable function. Then, the function θ ( t ) = I 0 κ ( t ) sign ( x ( t ) ) is uniformly continuous on [ 0 , ) .

3. Generalized Sliding Mode Control

Consider a dynamical system described by the first-order differential equation
x ˙ ( t ) = u ( t ) + Δ ( t , x ( t ) ) ,
where x : [ 0 , ) R is the system output, u : [ 0 , ) R is the control input and Δ : [ 0 , ) × R R is an unknown but continuous disturbance. Additionally, the following assumption is taken on Δ ( t , x ( t ) ) .
Assumption 1.
Disturbance Δ ( t , x ( t ) ) is globally bounded and satisfies the condition
[ Δ ( t , x ( t ) ) Δ ( a , x ( a ) ) ] I a κ ( t ) L [ a , ) ,
for arbitrary lower terminal a R .
The above implies that there is a function d L [ a , ) , with
Δ ( t , x ( t ) ) = Δ ( a , x ( a ) ) + I a κ ( t ) d ( t , x ( t ) ) .
In accordance with [30], the above condition on Δ , can be guaranteed (strengthened) whenever both the disturbance Δ ( t , x ( t ) ) and the kernel κ ( t ) have well-defined Laplace transforms and
lim z 0 λ ( z ) | Δ ( t + z , x ( t + z ) ) Δ ( t , x ( t ) ) | = 0 ,
for t , z > 0 , and λ ( t ) the kernel function that is paired with κ ( t ) . The condition on Δ ( t , x ( t ) ) in (25) is a weaker condition (for functions that have a Laplace transform) than integer-order, or even fractional-order, differentiability. In such a case, the function d ( t , x ( t ) ) acts as the generalized derivative of Δ ( t , x ( t ) ) with respect to λ ( t ) .
In addition to the properties stated in the above section, the following condition on the kernel κ ( t ) is required to develop the control scheme of this paper.
Assumption 2.
The solution ξ ( t ) of the non-autonomous and singular integro-differential equation
t ξ ( t ) + ξ ( t ) = 0 t κ ( τ ) d τ ,
with initial condition ξ ( 0 ) = 0 , is locally Lipschitz continuous, and satisfies the following differential inequality for any t 0 :
ξ ( t ) 1 1 + α κ ( t ) ,
and some α ( 0 , 1 ) .
The above assumption induces some interesting characteristics on κ ( t ) , which are presented in the following:
Proposition 5.
If κ ( t ) and ξ ( t ) satisfy (27), then ξ ( t ) is strictly monotonically increasing and concave for any t 0 . In addition, the range of ξ ( t ) is [ 0 , ) .
Proof. 
Since κ ( t ) is a real positive function for t > 0 , function ξ ( t ) is strictly monotonically increasing. Furthermore, taking the derivative in both sides of (26), we obtain
t ξ ( t ) + 2 ξ ( t ) = κ ( t ) .
Then, considering (27), we arrive at
t ξ ( t ) = κ ( t ) 2 ξ ( t ) κ ( t ) 2 1 + α κ ( t ) = 1 α 1 + α κ ( t ) < 0 ,
concluding that ξ ( t ) is a concave function for t 0 .
Finally, integrating (27) produces
lim t ξ ( t ) lim t 1 1 + α 0 t κ ( τ ) d τ = ,
and ξ ( t ) being continuous means that it assumes all the points in [ 0 , ) . □
The uniformly continuous controller and the stability of the closed-loop system are presented in the next main result:
Theorem 1.
Consider system (22) closed by the uniformly continuous controller
u ( t ) = u ( t n ) γ t n t κ ( t τ ) sign ( x ( τ ) ) d τ ,
with the kernel function κ ( t ) designed in accordance with Assumption 2, for some α > 0 , and γ the feedback gain that satisfies
γ > 3 + α 1 α ,
with t n , such that x ( t n ) = 0 . Then, there is a finite time t s R , such that x ( t ) 0 as t t s , and remains invariant thereafter.
Proof. 
The proof of Theorem 1 is broken into Lemmas 1–5. □
The closed-loop system can be rewritten as
x ˙ ( t ) = x ˙ ( t n ) I t n κ ( t ) [ γ sign ( x ( t ) ) d ( t , x ( t ) ) ] .
It is noticeable that the existence of solutions, x ( t ) , for the above nonlinear and singular integro-differential equation is supported, as sign ( x ( t ) ) and d ( t , x ( t ) ) are locally integrable, and these functions are uniformly continuous in every time interval [ t n , t n + 1 ] . In the case of this paper, we are not particularly interested in the uniqueness of the solution ( x ( t ) , x ˙ ( t ) ) , which cannot be taken for granted, but rather in the property that all solutions of (33) converge to the origin before a finite time, and remain invariant thereafter.
Lemma 1.
For the case x ˙ ( t 0 ) = x ( t 0 ) = 0 , we have x ( t ) = 0 for all t t 0 .
Proof. 
This case is straightforward as it can be appreciated that
x ( t ) x ˙ ( t ) = x ( t ) I t 0 κ ( t ) [ γ sign ( x ( t ) ) d ( t , x ( t ) ) ]
Then, as γ > sup | d ( t , x ( t ) ) | , we verify that the sliding mode condition
d d t 1 2 x 2 ( t ) = x ( t ) x ˙ ( t ) 0
sustains for all t t 0 . Therefore, x ˙ ( t 0 ) = x ( t 0 ) = 0 is invariant for any t t 0 . □
From now on, we consider the case x ˙ ( t 0 ) 0 . As x ( t ) is a continuous function and [ 0 , ) is a Lindelöf space, the set { t n } n N 0 of points where x ( t ) crosses, or hits, the zero value is at most a countable set. Thus, the problem is in demonstrating that t s = lim n t n < and lim n x ˙ ( t n ) = 0 , producing lim t t s x ( t ) = 0 , with x ( t ) = 0 an invariant sliding mode for any t t s .
It is clear that t 0 , the first instant with x ( t 0 ) = 0 , is finite. Then, let us start the analysis of t [ t 0 , t 1 ] , considering without loss of generality that x ˙ ( t 0 ) > 0 . Consider the functions ψ 1 ( t ) and ψ 2 ( t ) , with
ψ ˙ 1 ( t ) = x ˙ ( t 0 ) ( γ + γ d ) I κ ( t t 0 ) , ψ ˙ 2 ( t ) = x ˙ ( t 0 ) ( γ γ d ) I κ ( t t 0 ) , ψ 1 ( t ) = x ˙ ( t 0 ) ( t t 0 ) ( γ + γ d ) I 2 κ ( t t 0 ) , ψ 2 ( t ) = x ˙ ( t 0 ) ( t t 0 ) ( γ γ d ) I 2 κ ( t t 0 ) .
for γ d = ess sup t [ 0 , ) | d ( t , x ( t ) ) | , I κ ( t ) = 0 t κ ( τ ) d τ and I 2 κ ( t ) = 0 t I κ ( τ ) d τ .
Lemma 2.
For x ( t 0 ) > 0 , the time t 1 , when x ( t 1 ) = 0 , is bounded by
ξ ( t 1 t 0 ) ξ ( t ψ 2 t 0 ) = x ˙ ( t 0 ) γ γ d ,
and
I κ ( t ψ 2 t 0 ) 1 + α γ γ d x ˙ ( t 0 ) .
with t ψ 2 > t 0 , such that ψ 2 ( t ψ 2 ) = 0 .
Proof. 
It is possible to determine for t [ t 0 , t 1 ] , that
ψ ˙ 1 ( t ) x ˙ ( t ) ψ ˙ 2 ( t ) , ψ 1 ( t ) x ( t ) ψ 2 ( t ) .
In addition, it can be deducted that t 1 t ψ 2 , where t ψ 2 can be resolved from ψ 2 ( t ψ 2 ) = 0 , producing
ξ ( t ψ 2 t 0 ) = x ˙ ( t 0 ) γ γ d
since the function ξ ( t ) in (26) satisfies ξ ( t ) = I 2 κ ( t ) t , and its range is [ 0 , ) . Thus, combining the integral of (27) from 0 to ( t ψ 2 t 0 ) we obtain
ξ ( t ψ 2 t 0 ) 1 1 + α I κ ( t ψ 2 t 0 ) .
Thus, the first part of the proof follows from the fact ξ ( t ) is non-decreasing, while the second part follows by substituting (40) into (41). □
The next step is to analyze the behavior of x ˙ ( t ) around t 1 . which is bounded from below by ψ ˙ 1 ( t ψ 2 ) .
Lemma 3.
Under the conditions of Theorem 1 and Lemma 2, we have
| x ˙ ( t 1 ) | μ | x ˙ ( t 0 ) | ,
for
μ = 1 + ( 1 + α ) γ + γ d γ γ d ,
with μ ( 0 , 1 ) .
Proof. 
Signal x ˙ ( t ) is bounded from below by ψ ˙ 1 ( t ψ 2 ) in [ t 0 , t ψ 2 ] , where
ψ ˙ 1 ( t ψ 2 ) = x ˙ ( t 0 ) ( γ + γ d ) I κ ( t ψ 2 t 0 ) x ˙ ( t 0 ) ( 1 + α ) γ + γ d γ γ d x ˙ ( t 0 ) μ x ˙ ( t 0 ) ,
for
μ = 1 + ( 1 + α ) γ + γ d γ γ d .
Besides, as γ > γ d from (32), it is clear that μ > 0 . In addition, the tuning of γ in (32) guarantees that μ < 1 . □
The convergence time t s is analyzed as follows:
Lemma 4.
Under the conditions of Theorem 1 and Lemma 2, we have x ˙ ( t ) and x ( t ) 0 as t t s , for
t s t 0 + 1 1 μ ξ 1 | x ˙ ( t 0 ) | γ γ d .
Proof. 
This proof is divided in three parts:
Part I. Convergence x ˙ ( t ) 0 :
Considering the dynamic memory resetting and repeating the process for any time interval [ t n 1 , t n ] , we have | x ˙ ( t n ) | μ | x ( t n 1 ) | , which implies that | x ˙ ( t n ) | μ n | x ( t 0 ) | , producing x ˙ ( t n ) 0 as n .
Part II. Convergence of x ( t ) 0 as t t s :
We try finding max [ t 0 , t 1 ] x ( t ) , which can be estimated from max [ t 0 , t ψ 2 ] ψ 2 ( t ) . Additionally, as ψ 2 ( t ) is a concave function (its second derivative is negative almost everywhere), we can rely on
max [ t 0 , t 1 ] x ( t ) ψ 2 ( t ψ 2 ) ,
for t ψ 2 resolved from ψ ˙ 2 ( t ψ 2 ) = 0 (note that the continuity of ψ ˙ 2 ( t ) and the fact that ψ 2 ( t ) becomes negative after t ψ 2 imply that ψ 2 M < , and assure that t ψ 2 [ t 0 , t ψ 2 ] ), where
( t ψ 2 t 0 ) = I κ 1 x ˙ ( t 0 ) γ γ d ,
which leads to
ψ 2 ( t ψ 2 ) x ˙ ( t 0 ) I κ 1 x ˙ ( t 0 ) γ γ d .
Remembering that I κ 1 ( t ) is monotonically increasing, and its greatest lower bound is zero, the value of max [ t 0 , t 1 ] | x ( t ) | increases or decreases as | x ˙ ( t 0 ) | does. Repeating this process for every interval [ t n + 1 , t n ] , we conclude that
max [ t n , t n + 1 ] | x ( t ) | μ n | x ˙ ( t 0 ) | I κ 1 μ n | x ˙ ( t 0 ) | γ γ d .
Therefore, max [ t n , t n + 1 ] | x ( t ) | 0 as n .
Part III. Time of convergence:
Since ξ ( z ) is non-decreasing and concave for z 0 , function ξ 1 ( z ) is convex for z 0 . Now, as ξ ( 0 ) = 0 and ξ ( t ) is bijective for t 0 , we have
ξ 1 μ z μ ξ 1 z ,
for any μ ( 0 , 1 ) and z 0 . Thus, the bound for the time of convergence can be estimated from
t s = t 0 + lim N T N ,
where the partial sum T N = n = 0 N 1 ( t n + 1 t n ) = t N t 0 satisfies
T N = n = 0 N 1 ( t n + 1 t n ) n = 0 N 1 ξ 1 | x ˙ ( t n ) | γ γ d n = 0 N 1 ξ 1 μ n | x ˙ ( t 0 ) | γ γ d ξ 1 | x ˙ ( t 0 ) | γ γ d n = 0 N 1 μ n ξ 1 | x ˙ ( t 0 ) | γ γ d 1 μ N 1 μ .
Therefore, the results follow by taking the limit N .
The following result completes the proof of Theorem 1:
Lemma 5.
Under the conditions of Theorem 1 and Lemma 2, we have x ˙ ( t ) = x ( t ) = 0 for all t t s .
Proof. 
Consider the time instant t > t s , such that x ( t ) 0 . This leads to a contradiction, since that would imply that x ( t ) will cross the zero value at some time t ω 0 + 1 > t that is not in the set { t n } n N 0 . This is also evident from considering the candidate Lyapunov function V ( x ( t ) ) = 1 2 x 2 ( t ) , which has the time derivative V ˙ ( x ( t ) ) = x ( t ) I t s κ ( t ) [ γ sign ( x ( t ) ) d ( t , x ( t ) ) ] . Therefore, the result is evident after realizing that V ( x ( t s ) ) = 0 and V ˙ ( x ( t ) ) 0 after t = t s . □
Remark 1.
Considering the limit case α 1 in the requirement of Theorem 1 would imply adjusting γ to an infinite value, which contradicts the fact that γ R . For this reason, we should refrain from proposing kernel functions that satisfy (27), for positive α which are not strictly lesser than some positive number, in turns lesser than 1. Therefore, appealing kernels of the form κ ( t ) = 0 1 t α 1 d α are not supported by Theorem 1, and its study still remains as an unresolved problem.

4. Examples

Example 1
(Fractional sliding mode control). The simple kernel κ ( t ) = t α 1 , for α ( 0 , 1 ) , stands for a suitable design, satisfying the constraints of Theorem 1. It is possible to verify that ξ ( t ) = t α α ( 1 + α ) , resulting in
ξ ( t ) 1 1 + α κ ( t ) = t α 1 1 + α t α 1 1 + α = 0 .
The control formulation that arises in Equation (31) is referred to as continuous fractional sliding mode [4].
In order to design a new family of uniformly continuous sliding mode controllers, consider the following definition and proposition.
Definition 5.
The Mittag–Leffler function of two parameters is defined as
E α , β ( σ ) = k = 0 σ k Γ ( k α + β ) ,
for parameters α , β C , and variable argument σ C . The function Γ : C C is the Euler’s gamma function that satisfies Γ ( ν + 1 ) = ν Γ ( ν ) whenever { ν } > 0 .
Proposition 6.
Let α, β ( 0 , 1 ) and ϱ > 0 . Then,
( 1 + β ) E α , β + 2 ( ϱ t α ) E α , β + 1 ( ϱ t α ) 0 ,
for all t 0 .
The proof of the above proposition is given in the Appendix A.
It is worth mentioning that there are other values of α and β that would work in the above proposition, but those are not of interest for the controller design suggested below.
Example 2
(Mittag–Leffler sliding mode control). An interesting generalization arises by considering 0 < α β < 1 and ϱ > 0 . Then, the kernel
κ ( t ) = t β 1 E α , β + 1 ( ϱ t α )
produces
ξ ( t ) = t β 1 E α , β + 2 ( ϱ t α ) .
It can be noticed that the last proposition guarantees ξ ( t ) 1 1 + β κ ( t ) .
It is also clear that κ ( t ) satisfy the first two properties given in Section 2. For the third property, we have
0 t κ ( τ ) d τ = k = 0 ( ϱ ) k t k α + β ( k α + β ) Γ ( k α + β + 1 ) ,
indicating that lim t 0 0 t κ ( τ ) d τ = 0 . Finally, for the fourth property, we consider that
0 f ( t ) t d τ = 0 F ( s ) d s ,
where F ( s ) is the Laplace transform of f ( t ) . Thus, noticing that
κ ( t ) = t β E α , β + 1 ( ϱ t α ) t ,
and that the Laplace transform of t β E α , β + 1 ( ϱ t α ) is
L t β E α , β + 1 ( ϱ t α ) = s α ( β + 1 ) s α + ϱ ,
we have
0 κ ( τ ) d τ = 0 s α ( β + 1 ) s α + ϱ d s .
For the case of α = β , we have for any ε ( 0 , 1 )
0 κ ( τ ) d τ ε 1 1 s ( s α + ϱ ) d s 1 1 + ϱ ε 1 d s s = ln ( ε 1 ) 1 + ϱ ,
concluding that 0 κ ( τ ) d τ is greater than any real number.
Now, for the case of α < β , and an arbitrarily small positive number ε, we obtain
0 κ ( τ ) d τ ε 1 s α β 1 s α + ϱ d s 1 1 + ϱ ε 1 s α β 1 d s = ε ( β α ) 1 ( β α ) ( 1 + ϱ ) ,
leading again to 0 κ ( τ ) d τ = . The condition β α induces a non-integrable singularity around the origin, as was shown above.
The form of the Mittag–Leffler function-based kernel in the above example is given to facilitate the demonstration of the properties in Assumption 2. Other choices would, possibly, satisfy Assumption 2, but induce cumbersome derivations.
Remark 2.
In order to see that α > β induces ill-behaved kernel pairs, consider
0 κ ( τ ) d τ = 0 s α β 1 s α + ϱ d s = 0 1 s α β 1 s α + ϱ d s + 1 s α β 1 s α + ϱ d s 1 ϱ 0 1 s α β 1 d s + 1 s α β 1 s α d s = 1 ϱ ( α β ) + 1 β .
Thus, for α > β function κ ( τ ) fails to be a Sonine kernel, which is catastrophic as we cannot guarantee that the range of ξ ( t ) is [ 0 , ) , and the derivative of the disturbance with respect to λ ( t ) possesses intricate characteristics, beyond the scope of this paper.
Example 3
(Multi-order sliding mode control). It can be seen that if κ ( t ) satisfies (27), then c 0 κ ( t ) also satisfies (27) for arbitrary constant c 0 > 0 . Additionally, if the kernels κ 1 ( t ) and κ 2 ( t ) satisfy (27) for α 1 , α 2 ( 0 , 1 ) . Then, the kernel function κ ( t ) = c 1 κ 1 ( t ) + c 2 κ 2 ( t ) also satisfies (27) for α = max ( α 1 , α 2 ) . Following this reasoning, for the finite set of kernels { κ i ( t ) } i = 1 N that satisfy (27) for corresponding { α i } i = 1 n in ( 0 , 1 ) , we have κ ( t ) = i = 1 N c i κ i ( t ) satisfying (27) for α = max { α i } i = 1 N and a set of non-negative and not all zero constants { c i } i = 1 n .
Example 4
(Distributed-order sliding mode control). To further generalize the above designs, consider the kernel κ α ( t ) , satisfying (27) for some α ( α m , α M ) ( 0 , 1 ) , that is ξ α ( t ) 1 1 + α κ α ( t ) , where ξ α ( t ) is computed from κ α ( t ) . Multiplying the above inequality by a distribution function c ( α ) : ( 0 , 1 ) [ 0 , ) , with c ( α ) > 0 for some α ( α m , α M ) , and integrating with respect to α over ( α m , α M ) , we obtain
c ( α ) ξ α ( t ) 1 1 + α c ( α ) κ α ( t ) α m α M c ( α ) ξ α ( t ) d α α m α M c ( α ) 1 + α κ α ( t ) d α α m α M c ( α ) ξ α ( t ) d α 1 1 + α M α m α M c ( α ) κ α ( t ) d α ξ ( t ) 1 1 + α M κ ( t )
where κ ( t ) = α m α M c ( α ) κ α ( t ) d α , and ξ ( t ) is computed accordingly.

5. Simulations

Simulation and numerical implementations of continuous-time controllers face several challenges, such as discretization, fast dynamics and non-smooth effects, as well as other issues that are inherent to each control scheme.
The plant is the disturbed first-order integrator Equation (22), with initial conditions x ( 0 ) = 0 and x ˙ ( 0 ) = 1 . It is worth mentioning that the objective of this paper is to show the validity of the proposed scheme in a particular application scenario, not to show superiority with respect to previously reported results. However, the single integrator system (22) represents a large class of sliding dynamics, where x ( t ) stands for the sliding variable that must be driven to zero to assure that the whole system evolves ideally, without the effect of disturbances and uncertainties.
The simulations were programmed in Matlab, where the controller was obtained by considering the Euler method to compute the convolution. In the same sense, the Euler method with a sampling time of d t = 0.1 ms is considered.
Two different case are considered: (i) a disturbance free case and (ii) a disturbed case, to allow for a better comparison of the controllers that are presented. For these cases, the controllers of Examples 1 and 2 are studied. The fractional sliding mode controller considers the kernel κ ( t ) = t β 1 , with β = 0.5 , and the kernel for the Mittag–Leffler sliding mode controller relies on the same β (for a fair comparison), and α = 0.4 and ρ = 0.1 .
The same feedback gain value γ is chosen for both schemes, depending on each case, whose value can be determined based on a heuristic way (trial and error or educated guess), by gradually increasing its value to obtain an acceptable performance. The principle of dynamic memory resetting is applied every time when the condition x ( t ) x ( k + d t ) < 0 is fulfilled.

5.1. System without Disturbances

In this case, the feedback gain γ = 1 is considered, as there is no disturbance present in the system. The comparison results of Fractional sliding mode control vs. Mittag–Leffler sliding mode control are shown in Figure 1. It is possible to appreciate that the overshot is somehow lesser in the Mittag–Leffler case than in the fractional one. The control signals are similar and the phase diagram also shows a slightly faster convergence in the case of the Mittag–Leffler controller.

5.2. System with Disturbances

In this case, the feedback gain is set to γ = 2.0 to face system uncertainties. The disturbance is given by
d ( t ) = 0.5 sin ( 2 t ) + 0.25 sin ( 4 t ) + r ( t )
where r : [ 0 , ) [ 0.1 , 0.1 ] is a random valued function that updates every 10 ms, see Figure 2. The random pattern is the same for both controllers.
The comparison results are shown in Figure 3. As in the previous case, the output value function is slightly better for the Mittag–Leffler sliding mode controller. Both control signals are similar and reject the disturbance in a good extent.
This follows from the known robustness properties of fractional sliding mode controllers.
A more accurate comparison is obtained by relying on the integral square error (ISE) and integral square control (ISC) norms, obtaining the following results for both schemes:
  • Fractional sliding mode control:
    ISE = 7.4037 × 10 5 and ISC = 6.7932 .
  • Mittag–Leffler sliding mode control:
    ISE = 3.3341 × 10 5 and ISC = 6.7834 .
This shows that the Mittag–Leffler sliding mode control provides better regulation with less control energy, thus improving the closed-loop performance. It can be mentioned that considering kernels with lower values of β , the accuracy of the regulation task improves significantly. However, for real-world applications, the order should be kept as high as possible to improve the regularity of the control signal, such that the controller implementation is closer to the theoretical framework, even under the action of moderately fast actuators and conservative computational resources.

6. Discussions

The simulation results show comparable performances as those obtained by implementing well-established control techniques; nonetheless, the contribution of this paper is not limited to a particular class of sliding mode control methodology, rather, it stands for a generalization, and it opens the door to new families of robust control methods. Some limitations of the proposed work are the numerical implementation of the Sonine integral, as it depends on a convolution operation. Nonetheless, the kernel function can be evaluated beforehand to alleviate computational cost, such that the convolution can be approximated by a finite-impulse-response filter. Future research considers some open problems, such as the design of sliding mode controllers with the following kernels: variable-order kernels; Prabhakar functions; Non-Laplace transformable kernels; non-singular kernels; kernels with bounded integrals, etc. Furthermore, the applications of the present methodology in different systems, justifying its implementation in the presence of a large class of non-differentiable disturbances, constitute one of the most interesting avenues of this work for engineering professionals. The potential application cases consider, but are not limited to, the control of unmanned aerial robots, which are subject to gust winds and turbulence effects; robotic manipulators in free-motion with backlash; robots in constrained motion in contact with rough surfaces and diverse tribological phenomena; vehicles immersed in non-Newtonian flows; super-capacitive effects in electric networks; physical and engineering process with noisy measurement; renewable systems (wind and solar energy) subject to non-smooth wind speed and light intensity patterns, etc.

7. Conclusions

The contribution of this paper was proving that a class of singular kernels, even more general than fractional-order, multi-fractional and distributed-order kernels, can be used to propose several different families of uniformly uniformly continuous sliding mode controllers. As mentioned before, the continuous fractional sliding mode control, a powerful and robust methodology, constitutes a very particular case of the larger class of schemes studied in this paper. The present results are of potential interest for a wide range of control applications, where the plant is subject to a general class of disturbances and uncertainties, such as those with multi-fractal indices, or even non-differentiable disturbances that change their regularity over time.

Author Contributions

Conceptualization, A.J.M.-V.; methodology, A.J.M.-V. and G.F.-A.; software, A.J.M.-V.; validation, A.J.M.-V.; formal analysis, A.J.M.-V. and G.F.-A.; investigation, A.J.M.-V. and G.F.-A.; writing—original draft preparation, A.J.M.-V.; writing—review and editing, A.J.M.-V. and G.F.-A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Any data related to this work will be provided upon explicit request to the corresponding author of this article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Proof of Proposition 6

Proof. 
For the case of t = 0 , we have
E α , β + 2 ( 0 ) = 1 Γ ( β + 2 ) and E α , β + 1 ( 0 ) = 1 Γ ( β + 1 ) .
Therefore ( 1 + β ) E α , β + 2 ( 0 ) = E α , β + 1 ( 0 ) .
For the case of t > 0 , we obtain
( 1 + β ) E α , β + 2 ( ϱ t α ) E α , β + 1 ( ϱ t α ) = ( 1 + β ) k = 0 ( ϱ ) k t k α Γ ( k α + β + 2 ) k = 0 ( ϱ ) k t k α Γ ( k α + β + 1 ) = k = 0 ( ϱ ) k t k α 1 + β Γ ( k α + β + 2 ) 1 Γ ( k α + β + 1 ) = k = 0 ( ϱ ) k t k α 1 + β Γ ( k α + β + 2 ) k α + β + 1 Γ ( k α + β + 2 ) = k = 0 ( ϱ ) k ( k α ) t k α Γ ( k α + β + 2 ) .
It can be appreciated that
( k α ) t k α = t d t k α d t for   any t > 0 .
Then,
( 1 + β ) E α , β + 2 ( ϱ t α ) E α , β + 1 ( ϱ t α ) = t d d t k = 0 ( ϱ ) k t k α Γ ( k α + β + 2 ) = t E α , β + 2 ( ϱ t α ) .
Finally, since the function E α , β + 2 ( σ ) , with σ 0 , is completely monotonic for α ( 0 , 1 ) and β 0 [37], we have E α , β + 2 ( ϱ t α ) 0 for every t > 0 . Thereby, ( 1 + β ) E α , β + 2 ( ϱ t α ) E α , β + 1 ( ϱ t α ) , completing the proof. □

References

  1. Utkin, V.I.; Vadim, I. Sliding mode control. Var. Struct. Syst. Princ. Implement. 2004, 66, 1. [Google Scholar]
  2. Parra-Vega, V.; Hirzinger, G. Chattering-free sliding mode control for a class of nonlinear mechanical systems. Int. J. Robust Nonlinear Control IFAC Affil. J. 2001, 11, 1161–1178. [Google Scholar] [CrossRef]
  3. Utkin, V.; Poznyak, A.; Orlov, Y.; Polyakov, A. Conventional and high order sliding mode control. J. Frankl. Inst. 2020, 357, 10244–10261. [Google Scholar] [CrossRef]
  4. Muñoz-Vázquez, A.J.; Parra-Vega, V.; Sánchez-Orta, A. A novel continuous fractional sliding mode control. Int. J. Syst. Sci. 2017, 48, 2901–2908. [Google Scholar] [CrossRef]
  5. Muñoz-Vázquez, A.J.; Fernández-Anaya, G.; Sánchez-Torres, J.D. Adaptive parallel fractional sliding mode control. Int. J. Adapt. Control Signal Process. 2022, 36, 751–759. [Google Scholar] [CrossRef]
  6. Parra-Vega, V.; Arimoto, S.; Liu, Y.H.; Hirzinger, G.; Akella, P. Dynamic sliding PID control for tracking of robot manipulators: Theory and experiments. IEEE Trans. Robot. Autom. 2003, 19, 967–976. [Google Scholar] [CrossRef]
  7. Chen, C.T.; Peng, S.T. Design of a sliding mode control system for chemical processes. J. Process Control 2005, 15, 515–530. [Google Scholar] [CrossRef]
  8. Barambones, O. Sliding mode control strategy for wind turbine power maximization. Energies 2012, 5, 2310–2330. [Google Scholar] [CrossRef]
  9. Evangelista, C.; Valenciaga, F.; Puleston, P. Active and reactive power control for wind turbine based on a MIMO 2-sliding mode algorithm with variable gains. IEEE Trans. Energy Convers. 2013, 28, 682–689. [Google Scholar] [CrossRef]
  10. Sánchez-Torres, J.D.; Loukianov, A.G.; Galicia, M.I.; Domínguez, J.R. Robust nested sliding mode integral control for anti-lock brake system. Int. J. Veh. Des. 2013, 62, 188–205. [Google Scholar] [CrossRef]
  11. Hernández, A.G.G.; Fridman, L.; Levant, A.; Shtessel, Y.; Leder, R.; Monsalve, C.R.; Andrade, S.I. High-order sliding-mode control for blood glucose: Practical relative degree approach. Control Eng. Pract. 2013, 21, 747–758. [Google Scholar] [CrossRef]
  12. Nojavanzadeh, D.; Badamchizadeh, M. Adaptive fractional-order non-singular fast terminal sliding mode control for robot manipulators. IET Control Theory Appl. 2016, 10, 1565–1572. [Google Scholar] [CrossRef]
  13. Wang, H.; Mustafa, G.I.; Tian, Y. Model-free fractional-order sliding mode control for an active vehicle suspension system. Adv. Eng. Softw. 2018, 115, 452–461. [Google Scholar] [CrossRef]
  14. Ardjal, A.; Mansouri, R.; Bettayeb, M. Fractional sliding mode control of wind turbine for maximum power point tracking. Trans. Inst. Meas. Control 2019, 41, 447–457. [Google Scholar] [CrossRef]
  15. Efe, M.Ö. Fractional order sliding mode control with reaching law approach. Turk. J. Electr. Eng. Comput. Sci. 2010, 18, 731–748. [Google Scholar] [CrossRef]
  16. Efe, M.Ö. Integral sliding mode control of a quadrotor with fractional order reaching dynamics. Trans. Inst. Meas. Control 2011, 33, 985–1003. [Google Scholar] [CrossRef]
  17. Yin, C.; Chen, Y.; Zhong, S.m. Fractional-order sliding mode based extremum seeking control of a class of nonlinear systems. Automatica 2014, 50, 3173–3181. [Google Scholar] [CrossRef]
  18. Yu, L.; Huang, J.; Luo, W.; Chang, S.; Sun, H.; Tian, H. Sliding-Mode Control for PMLSM Position Control—A Review. Actuators 2023, 12, 31. [Google Scholar] [CrossRef]
  19. Ye, Z.; Zhang, D.; Deng, C.; Yan, H.; Feng, G. Finite-time resilient sliding mode control of nonlinear UMV systems subject to DoS attacks. Automatica 2023, 156, 111170. [Google Scholar] [CrossRef]
  20. Rezaei, A.; Tabatabaei, M. Ship roll stabilization using an adaptive fractional-order sliding mode controller. Ocean. Eng. 2023, 287, 115883. [Google Scholar] [CrossRef]
  21. Giap, V.N. Text message secure communication based on fractional-order chaotic systems with Takagi–Sugeno fuzzy disturbance observer and sliding mode control. Int. J. Dyn. Control 2023, 11, 3109–3123. [Google Scholar] [CrossRef]
  22. Giap, V.N.; Pham, D.H.; Lin, C.M.; Trinh, M.H.; Le, M.C.; Nguyen, Q.D. Fractional Order T–S Fuzzy Chaotic Models for Secure EEG Signal via a Wireless Communication Protocol Using a Disturbance Observer and Sliding Mode Control Technique. Int. J. Fuzzy Syst. 2024, 1–19. [Google Scholar] [CrossRef]
  23. Do, T.A.; Pham, Q.H.; Giap, V.N.; Nguyen, Q.D.; Vu, P. Homogeneous Disturbance Observer Based on Sliding Mode Observer and Controller for T-Type Inverter. J. Electr. Eng. Technol. 2024, 19, 1671–1683. [Google Scholar] [CrossRef]
  24. Shi, H.; Zeng, J.; Guo, J. Disturbance observer-based sliding mode control of active vertical suspension for high-speed rail vehicles. Veh. Syst. Dyn. 2024, 1–24. [Google Scholar] [CrossRef]
  25. Dong, H.Q.; Gam, N.T.; Cuong, H.M.; Tuan, L.A. Fractional-order fast terminal back-stepping sliding mode control of autonomous robotic excavators. J. Frankl. Inst. 2024, 361, 106686. [Google Scholar] [CrossRef]
  26. Sonine, N. Sur la généralisation d’une formule d’Abel. Acta Math. 1884, 4, 171–176. [Google Scholar] [CrossRef]
  27. Luchko, Y. General fractional integrals and derivatives with the Sonine kernels. Mathematics 2021, 9, 594. [Google Scholar] [CrossRef]
  28. Luchko, Y. General fractional integrals and derivatives of arbitrary order. Symmetry 2021, 13, 755. [Google Scholar] [CrossRef]
  29. Luchko, Y. Operational calculus for the general fractional derivative and its applications. Fract. Calc. Appl. Anal. 2021, 24, 338–375. [Google Scholar] [CrossRef]
  30. Muñoz-Vázquez, A.J.; Fernández-Anaya, G.; Martínez-Fuentes, O. Stability analysis of a class of integral equations with not necessarily differentiable solutions. J. Comput. Appl. Math. 2021, 398, 113702. [Google Scholar] [CrossRef]
  31. Muñoz-Vázquez, A.J.; Martínez-Fuentes, O.; Fernández-Anaya, G. Generalized PI control for robust stabilization of dynamical systems. Math. Comput. Simul. 2022, 202, 22–35. [Google Scholar] [CrossRef]
  32. Garrappa, R.; Giusti, A.; Mainardi, F. Variable-order fractional calculus: A change of perspective. Commun. Nonlinear Sci. Numer. Simul. 2021, 102, 105904. [Google Scholar] [CrossRef]
  33. Beghin, L.; Cristofaro, L.; Garrappa, R. Renewal processes linked to fractional relaxation equations with variable order. J. Math. Anal. Appl. 2024, 531, 127795. [Google Scholar] [CrossRef]
  34. Ortigueira, M.D. Searching for Sonin kernels. Fract. Calc. Appl. Anal. 2024, 1–29. [Google Scholar] [CrossRef]
  35. Baleanu, D.; Shiri, B. Collocation methods for fractional differential equations involving non-singular kernel. Chaos Solitons Fractals 2018, 116, 136–145. [Google Scholar] [CrossRef]
  36. Fernandez, A.; Al-Refai, M. A Rigorous analysis of integro-differential operators with non-singular kernels. Fractal Fract. 2023, 7, 213. [Google Scholar] [CrossRef]
  37. Miller, K.S. A note on the complete monotonicity of the generalized Mittag-Leffler function. Real Anal. Exch. 1999, 22, 753–755. [Google Scholar] [CrossRef]
Figure 1. Disturbance free case comparison: Fractional sliding mode control (left column) vs. Mittag–Leffler sliding mode control (right column).
Figure 1. Disturbance free case comparison: Fractional sliding mode control (left column) vs. Mittag–Leffler sliding mode control (right column).
Mathematics 12 02536 g001
Figure 2. Unknown disturbance.
Figure 2. Unknown disturbance.
Mathematics 12 02536 g002
Figure 3. Disturbed case comparison: Fractional sliding mode control (left column) vs. Mittag–Leffler sliding mode control (right column).
Figure 3. Disturbed case comparison: Fractional sliding mode control (left column) vs. Mittag–Leffler sliding mode control (right column).
Mathematics 12 02536 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Muñoz-Vázquez, A.J.; Fernández-Anaya, G. Uniformly Continuous Generalized Sliding Mode Control. Mathematics 2024, 12, 2536. https://doi.org/10.3390/math12162536

AMA Style

Muñoz-Vázquez AJ, Fernández-Anaya G. Uniformly Continuous Generalized Sliding Mode Control. Mathematics. 2024; 12(16):2536. https://doi.org/10.3390/math12162536

Chicago/Turabian Style

Muñoz-Vázquez, Aldo Jonathan, and Guillermo Fernández-Anaya. 2024. "Uniformly Continuous Generalized Sliding Mode Control" Mathematics 12, no. 16: 2536. https://doi.org/10.3390/math12162536

APA Style

Muñoz-Vázquez, A. J., & Fernández-Anaya, G. (2024). Uniformly Continuous Generalized Sliding Mode Control. Mathematics, 12(16), 2536. https://doi.org/10.3390/math12162536

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