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 and be two non-negative and strictly monotonically decreasing functions that satisfy the so-called Sonine condition Functions and constitute a kernel pair.
In this paper, we also consider that if
belongs to a kernel pair, satisfies the following properties [
30]:
The following definitions are important for the subsequent analysis and the control proposal: [
27,
30]:
Definition 2 (Generalized integral)
. Let and , with an absolutely continuous functions on , for any . The functionfor , is the generalized integral of with respect to the kernel . Definition 3 (Generalized derivative)
. Let be a sufficiently regular continuous function, such that the expressionis real-valued. Then, is the generalized derivative of with respect to , and with lower terminal . Definition 4 (Generalized differentiability)
. Whenever is locally bounded, for all and , for some , function is called marginally -differentiable (for lower terminal a). In addition, if satisfiesThen, is called -differentiable (for lower terminal a).
The following proposition relates the Definitions 3 and 4:
Proposition 1. Let be a locally bounded and -differentiable function, for some lower terminal . Then, is locally bounded.
Proof. The term
in the right-hand side of (
3) is locally bounded since
leads to
, which vanishes at
.
Let us set some
. Thus, consider a constant
, such that
Then, for an arbitrarily small
, we have
which is finite. In addition, consider
as well as the sequence of functions
for a sufficiently small
. It is clear that
for all
, and that
as
. Thus, using the Fatou’s Lemma, we obtain
Moreover, as
for any
, both
and
are non-negative, and
is integrable and goes to zero as
, the integral in the right-hand side of the last inequality vanishes as
. Thus, given
, we can chose
as small as required, such that
consequently,
The limit case is obtained in a similar fashion by choosing , and repeating the above process for , obtaining . □
In a more general case when
is only marginally
-differentiable, it is difficult to show that
exists. Nonetheless, in accordance with [
30], if
can be written as a constant plus the generalized integral with respect to
of some function
, this is,
, it results that
.
In order to prove that the generalized differentiability notion is weaker than integer-order differentiability, consider the following:
Proposition 2. Let be a Lipschitz continuous real-valued function. Then is -differentiable.
Proof. By virtue of
is Lipschitz continuous, there is some constant
, such that
for any real
t and
z. Therefore, for
, we have
Then, as both and z are continuous functions for , the function is continuous, and thus locally bounded for any .
For the limit case
, consider that
whenever
. Whereby, integrating with respect to
t from 0 to
z, we obtain
Therefore, as the last integral vanishes when , it follows that the Lipchitz continuous function is -differentiable. □
It is worth remembering that (
3) coincides with the generalized Caputo-like derivative
, whenever
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 be a function with well-posed generalized integral and be at least marginally -differentiable function. Then,where and . An interesting property of the fractional integral in (
2) is unveiled as follows:
Proposition 4. Let . Then, is continuous over .
Proof. Let
be arbitrarily small. Consider a couple of instants
t and
, and without loss of generality assume
(as the case
is trivial and the case
is analogous), this is,
, for some
. Then, we obtain
First, if we define
, since
for every
, we obtain
where the change of variable
and terminals inversion were performed to obtain the last integral. Moreover, as
is absolutely continuous, there is a small enough
, such that
.
In the above expression, we can notice that
Therefore, and the generalized integral of a locally bounded function is continuous. □
The following two Corollaries expose additional properties:
Corollary 1. If is replaced by , we verify that is uniformly continuous on .
Corollary 2. Let be a Lebesgue measurable function. Then, the function is uniformly continuous on .
3. Generalized Sliding Mode Control
Consider a dynamical system described by the first-order differential equation
where
is the system output,
is the control input and
is an unknown but continuous disturbance. Additionally, the following assumption is taken on
.
Assumption 1. Disturbance is globally bounded and satisfies the conditionfor arbitrary lower terminal . The above implies that there is a function
, with
In accordance with [
30], the above condition on
, can be guaranteed (strengthened) whenever both the disturbance
and the kernel
have well-defined Laplace transforms and
for
, and
the kernel function that is paired with
. The condition on
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
acts as the generalized derivative of
with respect to
.
In addition to the properties stated in the above section, the following condition on the kernel is required to develop the control scheme of this paper.
Assumption 2. The solution of the non-autonomous and singular integro-differential equationwith initial condition , is locally Lipschitz continuous, and satisfies the following differential inequality for any :and some . The above assumption induces some interesting characteristics on , which are presented in the following:
Proposition 5. If and satisfy (27), then is strictly monotonically increasing and concave for any . In addition, the range of is . Proof. Since
is a real positive function for
, function
is strictly monotonically increasing. Furthermore, taking the derivative in both sides of (
26), we obtain
Then, considering (
27), we arrive at
concluding that
is a concave function for
.
Finally, integrating (
27) produces
and
being continuous means that it assumes all the points in
. □
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 controllerwith the kernel function designed in accordance with Assumption 2, for some , and γ the feedback gain that satisfieswith , such that . Then, there is a finite time , such that as , and remains invariant thereafter. Proof. The proof of Theorem 1 is broken into Lemmas 1–5. □
The closed-loop system can be rewritten as
It is noticeable that the existence of solutions,
, for the above nonlinear and singular integro-differential equation is supported, as
and
are locally integrable, and these functions are uniformly continuous in every time interval
. In the case of this paper, we are not particularly interested in the uniqueness of the solution
, 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 , we have for all .
Proof. This case is straightforward as it can be appreciated that
Then, as
, we verify that the sliding mode condition
sustains for all
. Therefore,
is invariant for any
. □
From now on, we consider the case . As is a continuous function and is a Lindelöf space, the set of points where crosses, or hits, the zero value is at most a countable set. Thus, the problem is in demonstrating that and , producing , with an invariant sliding mode for any .
It is clear that
, the first instant with
, is finite. Then, let us start the analysis of
, considering without loss of generality that
. Consider the functions
and
, with
for
,
and
.
Lemma 2. For , the time , when , is bounded byandwith , such that . Proof. It is possible to determine for
, that
In addition, it can be deducted that
, where
can be resolved from
, producing
since the function
in (
26) satisfies
, and its range is
. Thus, combining the integral of (
27) from 0 to
we obtain
Thus, the first part of the proof follows from the fact
is non-decreasing, while the second part follows by substituting (
40) into (
41). □
The next step is to analyze the behavior of around . which is bounded from below by .
Lemma 3. Under the conditions of Theorem 1 and Lemma 2, we haveforwith . Proof. Signal
is bounded from below by
in
, where
for
Besides, as
from (
32), it is clear that
. In addition, the tuning of
in (
32) guarantees that
. □
The convergence time is analyzed as follows:
Lemma 4. Under the conditions of Theorem 1 and Lemma 2, we have and as , for Proof. This proof is divided in three parts:
Part I. Convergence :
Considering the dynamic memory resetting and repeating the process for any time interval , we have , which implies that , producing as .
Part II. Convergence of as :
We try finding
, which can be estimated from
. Additionally, as
is a concave function (its second derivative is negative almost everywhere), we can rely on
for
resolved from
(note that the continuity of
and the fact that
becomes negative after
imply that
, and assure that
), where
which leads to
Remembering that
is monotonically increasing, and its greatest lower bound is zero, the value of
increases or decreases as
does. Repeating this process for every interval
, we conclude that
Therefore, as .
Part III. Time of convergence:
Since
is non-decreasing and concave for
, function
is convex for
. Now, as
and
is bijective for
, we have
for any
and
. Thus, the bound for the time of convergence can be estimated from
where the partial sum
satisfies
Therefore, the results follow by taking the limit .
□
The following result completes the proof of Theorem 1:
Lemma 5. Under the conditions of Theorem 1 and Lemma 2, we have for all .
Proof. Consider the time instant , such that . This leads to a contradiction, since that would imply that will cross the zero value at some time that is not in the set . This is also evident from considering the candidate Lyapunov function , which has the time derivative . Therefore, the result is evident after realizing that and after . □
Remark 1. Considering the limit case in the requirement of Theorem 1 would imply adjusting γ to an infinite value, which contradicts the fact that . 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 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 , for , stands for a suitable design, satisfying the constraints of Theorem 1. It is possible to verify that , resulting inThe 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 asfor parameters , and variable argument . The function is the Euler’s gamma function that satisfies whenever . Proposition 6. Let α, and . Then,for all . 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 and . Then, the kernelproducesIt can be noticed that the last proposition guarantees .
It is also clear that satisfy the first two properties given in Section 2. For the third property, we haveindicating that . Finally, for the fourth property, we consider thatwhere is the Laplace transform of . Thus, noticing thatand that the Laplace transform of iswe have For the case of , we have for any concluding that is greater than any real number. Now, for the case of , and an arbitrarily small positive number ε, we obtainleading again to . 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 Thus, for function fails to be a Sonine kernel, which is catastrophic as we cannot guarantee that the range of is , and the derivative of the disturbance with respect to possesses intricate characteristics, beyond the scope of this paper.
Example 3 (Multi-order sliding mode control)
. It can be seen that if satisfies (27), then also satisfies (27) for arbitrary constant . Additionally, if the kernels and satisfy (27) for . Then, the kernel function also satisfies (27) for . Following this reasoning, for the finite set of kernels that satisfy (27) for corresponding in , we have satisfying (27) for and a set of non-negative and not all zero constants . Example 4 (Distributed-order sliding mode control)
. To further generalize the above designs, consider the kernel , satisfying (27) for some , that is where is computed from . Multiplying the above inequality by a distribution function , with for some , and integrating with respect to α over , we obtainwhere , and 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
and
. 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
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 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 , with , and the kernel for the Mittag–Leffler sliding mode controller relies on the same (for a fair comparison), and and .
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 is fulfilled.
5.1. System without Disturbances
In this case, the feedback gain
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
to face system uncertainties. The disturbance is given by
where
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 and ISC .
Mittag–Leffler sliding mode control:
ISE and ISC .
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.