Next Article in Journal
The Soliton Solutions of the Stochastic Shallow Water Wave Equations in the Sense of Beta-Derivative
Previous Article in Journal
Price Prediction of Bitcoin Based on Adaptive Feature Selection and Model Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Sliding Mode Control Algorithm with Elementary Compensation for Input Matrix Uncertainty in Affine Systems

State Key Laboratory for Manufacturing System Engineering, Systems Engineering Institute, Xi’an Jiaotong University, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(6), 1337; https://doi.org/10.3390/math11061337
Submission received: 8 February 2023 / Revised: 27 February 2023 / Accepted: 1 March 2023 / Published: 9 March 2023
(This article belongs to the Section Dynamical Systems)

Abstract

:
This paper aims to develop a sliding mode control (SMC) approach with elementary compensation for input matrix uncertainty in affine systems. As a multiplicative uncertainty regarding the control inputs, input matrix uncertainty adversely modifies the control effort and even further causes the instability of systems. To solve this issue, a sliding mode control algorithm is developed based on a two-step design strategy. The first step is to design a general sliding mode controller for the system without input matrix uncertainty. In the second step, a control term is specially designed to compensate for input matrix uncertainty. In order to realize the elementary compensation for input matrix uncertainty, this term is obtained by solving a nonlinear vector equation which is derived from the Lyapunov function inequality. Theorems and lemmas based on the convex cone theory are proposed to guarantee the existence and uniqueness of the solution to the vector equation. Additionally, an algorithmic process is proposed to solve the vector equation efficiently. In the simulation part, the proposed controller is applied to two systems with different structures and compared with two state-of-the-art SMC algorithms. The comprehensive simulation results demonstrate that the proposed method is able to provide the closed-loop system with a competitive performance in terms of convergence level, overshoot reduction and chattering suppression.

1. Introduction

Sliding mode control (SMC) is an effective approach to affine uncertain systems due to its advantages such as strong robustness, fast response and small overshoot [1,2]. These benefits stem from the property that the system is completely insensitive to matched uncertainties during sliding mode [3,4,5,6,7]. However, such a desirable characteristic depends on a proper control algorithm which is capable of enforcing the sliding mode. Thus, designing a controller to ensure the enforcement of sliding mode is a crucial subject for the SMC method.
Matched uncertainty is a challenging issue that poses a threat to the enforcement of the sliding mode and brings difficulties in the controller design process. In general, there are two types of matched uncertainties, each of which affects the system on different levels: (1) the additive uncertainty influences the system in a parallel way with the control inputs. (2) Input matrix uncertainty, which is multiplicative with respect to the control inputs, refracts the control effort by modifying the gain of the control loop; additionally, the occurrence of the sliding mode will be destroyed if the norm of input matrix uncertainty exceeds a certain value; additionally, for multi-input systems, input matrix uncertainty implies the presence of uncertainty in the control direction [8]. It has been discovered that input matrix uncertainty is more difficult to handle and numerous control strategies have been proposed to suppress its effect [8,9,10,11]. For instance, Bejarano-Rincón et al. [8] design a sliding mode controller for MIMO nonlinear systems with uncertain control direction; nevertheless, the coefficient of each discontinuous control is excessively large as it bounds above the sum of input uncertainties and other system components. One representative method is proposed by Cao and Xu [10], which develops a compensation term for input matrix uncertainty by iteratively relaxing the inequality regarding the derivative of the Lyapunov function. The efficacy of this approach has led to its widespread adoption by many researchers: Chang [12] proposes a control law for the linear MIMO system based on the information of the full-order compensator and the system outputs and addresses input matrix uncertainty with Cao and Xu’s strategy; Shtessel et al. [13] present a unit-vector approach to control the linear systems with matched uncertainties, where input matrix uncertainty is handled by Cao and Xu’s method; Zhang et al. [14] put forward a robust control design based on a new switched robust integral sliding mode (SRISM) for switched systems with uncertain input matrix and compensate for the uncertainties using Cao and Xu’s strategies; more recently, Lin and Hu [15] propose a robust SMC to achieve chaos synchronization for master-slave chaotic systems with matched/mismatched disturbances and handle input uncertainty with Cao and Xu’s scheme.
Nevertheless, the aforementioned methods exhibit a high degree of conservativeness, i.e., the use of triangle inequality results in the continuous amplification of the Lyapunov inequality and further causes an overlarge magnitude of the control gain. Therefore, the excessive control gain will result in amplified amplitude of chattering [16,17]. To address this issue, several continuous SMC strategies have been proposed recently [18,19,20]. For example, Martínez-Fuentes et al. [19] propose a Lipschitz continuous sliding mode controller (LCSMC), which ensures the global finite-time convergence to the sliding manifold despite an unknown control coefficient. Gutiérrez-Oribio et al. [20] propose two continuous second-order sliding mode controllers for a class of underactuated mechanical systems with two degrees of freedom, which are robust against the uncertainty in the control coefficient. However, only the single-input systems are studied in the literature studies above. Recently, Feng et al. [21] proposed a chattering-reduction sliding mode control method for affine uncertain systems, in which a control term is specially designed to compensate for input matrix uncertainty and the term is obtained by solving a nonlinear equation, so that the magnitude of control is effectively decreased and chattering is further suppressed.
Despite the fact that many studies have shown effective compensation for input matrix uncertainty, to the best of our knowledge, few studies have centered on elementarily compensating for input matrix uncertainty. Nasiri et al. [9] propose an adaptive sliding mode control scheme for an affine class of multi-input multi-output (MIMO) nonlinear systems with uncertain control distribution gain; however, the bounds of the input uncertainties are quite strict for a large number of systems. In practical applications, where the input matrix has a physical meaning and specific information for each component, it is critical to be informed about the elements of input matrix uncertainty. In the spacecraft attitude control system, for example, the uncertainty in the control gain matrix represents the faults of actuators and the elements of the matrix represent the execution effectiveness loss of the corresponding actuator (a certain reaction flywheel) [22]. The most common strategy for dealing with uncertainties in the gain matrix is to employ the observer solely for fault estimation and incorporate the gain adaption scheme to realize the controller design. However, the implementation complexity is high and the fault diagnosis introduces an unavoidable time delay, which will inevitably degrade the control performance.
In order to reduce the design conservativeness and simplify the practical implementation, this paper proposes a novel sliding mode algorithm to control the affine systems with input matrix uncertainty. The proposed algorithm is based on a two-step design strategy. More specifically, in the first step, a general SMC is designed for the system without input matrix uncertainty. In the second step, a control term is specifically developed to compensate for input matrix uncertainty by element. By solving a nonlinear vector equation derived from the derivative of the Lyapunov function, the compensation term for input matrix uncertainty is obtained. A theorem and several lemmas based on the convex cone theory are proposed to demonstrate the existence and uniqueness of the solution to the vector equation. In addition, an algorithmic process is proposed to solve the vector equation efficiently. As a result, the main contributions of this work are listed as follows,
  • A two-step design strategy is utilized to develop the controller: the first step aims at the additive matched uncertainties; the second step focuses on addressing input matrix uncertainty.
  • In the proposed method, a new compensation term is specially developed to suppress the input matrix uncertainty. This control term is determined by solving a nonlinear vector equation. By this means, the uncertainty in the input matrix can be compensated for by the element. A theorem and several lemmas based on the principle of the convex cone set are proposed to guarantee the existence and uniqueness of the solution to the vector equation.
  • A practical algorithmic process is proposed to obtain the unique solution to the vector equation efficiently.
  • The proposed method is applied to two systems with different structures. The simulations are conducted by comparing the proposed SMC strategy with the SMC method proposed by Cao and Xu [10] and the SMC method proposed by Feng et al. [21]. The simulation results show that the proposed method can successfully suppress the effect of input matrix uncertainty and is effective in chattering reduction.
  • Compared with SMC studies that integrated with intelligent control approaches such as observers, fuzzy logics and adaptive control, the proposed method has much simpler structure and many fewer control parameters and is useful for practical implementations.
The rest of this paper is organized as follows. Section 2 is devoted to system, sliding surfaces descriptions. In Section 3, the proposed control algorithm with theoretical analyses and the calculation process are given. Section 4 presents two numerical examples with simulation results and comparisons with other state-of-the-art methods. Section 5 draws the conclusion.
Notation 1. 
The notation used in this paper is fairly standard. Specially, R + denotes the set of real numbers and R n denotes the set of n × 1 real column vectors. Furthermore, · T is written for transpose, · denotes the absolute value for a real number and the absolute value of each element in a matrix. · denotes the Euclidean norm and its induced matrix norm, · 1 denotes the Taxicab norm. n denotes the set 1 , , n . If x R n is a vector, [ x ] i denotes the i-th element of x , i n . If A R n × m is a matrix, [ A ] i , j denotes the element in the i-th row and j-th column of A and i n , j m . n denotes a matrix set n = S S = diag ζ 1 , , ζ n , ζ i { 1 , 1 } , which contains 2 n members of n × n matrix. S ( x ) : R n n denotes matrix function in which each element fills the sign of the corresponding component in x , i.e., S ( x ) = diag sign [ x ] i and sign [ x ] i = 1 , [ x ] i 0 1 , [ x ] i < 0 .

2. Preliminaries of Sliding Mode Controller Design

Consider the following affine uncertain systems [10,23]:
x ˙ ( t ) = f ( x , t ) + B ( x , t ) I + Δ B ( x , t ) u ( t ) + Δ f m ( x , t ) ,
where t [ t 0 , ) represents time, x R n is a vector of measurable states, u R m is a vector of control inputs, f ( x , t ) and B ( x , t ) are known vector and matrix with appropriate size and their elements are piecewise continuous functions with respect to x and t [13]. Δ f m ( x , t ) is an unknown vector representing the additive matched uncertainties and Δ B ( x , t ) represents the input matrix uncertainty.
Assumption 1 ([24]). 
Assume that B ( x , t ) has full rank, i.e., rank ( B ( x , t ) ) = m .
Assumption 2 ([10,23]). 
For all x , t R n × R + , Δ f m ( x , t ) is bounded by a known nonlinear function as Δ f m x , t ρ m x , t .
Assumption 3 ([10]). 
For all x , t R n × R + , Δ B ( x , t ) R m × m satisfies
[ Δ B ( x , t ) ] i , j [ Δ B ¯ ] i , j , i m , j m Δ B ¯ ρ b < 1 .
where [ Δ B ( x , t ) ] i , j is the element in the i-th row and the j-th column of Δ B ( x , t ) , [ Δ B ¯ ] i , j is the upper bound of [ Δ B ( x , t ) ] i , j and is the element in the i-th row and the j-th column of Δ B ¯ . ρ b is a known positive constant which bounds the upper bound matrix Δ B ¯ .
In this paper, two kinds of sliding surfaces are investigated as follows.
  • If system (1) can be put in the regular canonical form [24,25,26]:
    x ˙ 1 ( t ) = f 1 ( x 1 , x 2 ) , x ˙ 2 ( t ) = f 2 ( x 1 , x 2 ) + B 2 ( x 1 , x 2 ) I + Δ B ( x 1 , x 2 ) u ( t ) + Δ f m ( x 1 , x 2 ) ,
    where x 1 R n m , x 2 R m , x T = x 1 T , x 2 T . B 2 ( x 1 , x 2 ) R m and det B 2 ( x 1 , x 2 ) 0 . f T ( x , t ) = f 1 T ( x 1 , x 2 ) , f 2 T ( x 1 , x 2 ) . B T ( x , t ) = 0 m × ( n m ) , B 2 T ( x 1 , x 2 ) , where 0 m × ( n m ) is the zero matrix of size m × ( n m ) [27]. Then, a linear sliding surface with the following form is adopted [24,28]:
    σ ( t ) = x 2 ( t ) + G 1 x 1 ( t ) = Gx ( t ) = 0 ,
    where G 1 R m × ( n m ) is a diagonal positive definite matrix. G = G 1 , I m × m R m × n , where I m × m is an m × m identity matrix.
  • Otherwise, the integral sliding surface is adopted as [29]
    σ ( t ) = Gx ( t ) Gx ( 0 ) 0 t G f ( x , τ ) + B ( x , τ ) κ ( x , τ ) d τ = 0 ,
    where G R m × n is a constant matrix; κ ( x , t ) is a nominal control under which the system x ˙ ( t ) = f ( x , t ) + B ( x , t ) κ ( x , t ) is globally asymptotically stable [10,23].
The derivative of sliding surface (4) and (5) can be integrated into one equation, Equation (6), as follows:
σ ˙ ( t ) = Gf σ ( x , t ) + GB Δ f m ( x , t ) + G B ( x , t ) ( I + Δ B ( x , t ) ) u ( t ) ,
where GB ( x , t ) is uniformly invertible and
f σ ( x , t ) = f ( x , t ) , if   σ ( t )   is   linear B ( x , t ) κ ( x , t ) , if   σ ( t )   is   integral .
Item (7) covers the different results corresponding to the linear surface and integral surface, respectively. Therefore, it can be directly used in the following derivation of the controller design. Namely, The controller design strategy proposed in this paper is applicable to both linear and integral sliding surfaces.

3. Sliding Mode Controller Design

In this section, we present the sliding mode controller design for system (1). The design proceeds in two step: In the first step, a general sliding mode controller is developed for the system without considering input matrix uncertainty; in the second step, a control term is designed specially to compensate for input matrix uncertainty. Therefore, the final controller for system (1) comprises two parts and is formulated as
u ( t ) = u d ( t ) + u c ( t ) ,
where u d ( t ) is the general sliding mode controller for a system only with uncertainty Δ f m ( x , t ) and u c ( t ) is the compensation control term for Δ B ( x , t ) .

3.1. General Control Law

The first step aims at the following system only with additive uncertainty Δ f m ( x , t ) :
x ˙ ( t ) = f ( x , t ) + B ( x , t ) u ( t ) + Δ f m ( x , t ) .
Since GB ( x , t ) is uniformly invertible, u d ( t ) can be designed as [10]
u d ( t ) = ( G B ( x , t ) ) 1 G f σ ( x , t ) ρ ( x , t ) v 0
where v = ( GB ( x , t ) ) T σ , v 0 = v / v , v > 0 0 , v = 0 and ρ ( x , t ) satisfies
ρ ( x , t ) > ρ m ( x , t ) .
To verify that u d ( t ) in (10) is able to enforce the sliding mode with regard to system (9), the following Lyapunov function is selected [10]:
V σ = 1 / 2 σ T σ .
Differentiating V σ with respect to time along with (9)–(11), one obtains
V ˙ σ = σ T Gf σ ( x , t ) + GB ( x , t ) ( G B ( x , t ) ) 1 G f σ ( x , t ) ρ ( x , t ) v 0 + Δ f m ( x , t ) = v T ρ ( x , t ) v 0 + Δ f m ( x , t ) ρ ( x , t ) v + Δ f m ( x , t ) v ρ m ( x , t ) v ρ ( x , t ) v < 0 .
Hence, the enforcement of the sliding mode is guaranteed by the control law u d ( t ) . Alternatively, u d ( t ) can be designed with another form [3]:
u d ( t ) = ( GB ( x , t ) ) 1 Gf σ ( x , t ) + ρ ( x , t ) sign ( σ ) ,
where ρ ( x , t ) > GB ( x , t ) ρ m ( x , t ) .
Choosing Lyapunov function as (12) and differentiating the function with respect to time along with (9), (11) and (14), one obtains:
V ˙ σ = σ T Gf σ ( x , t ) + GB ( x , t ) ( G B ( x , t ) ) 1 Gf σ ( x , t ) + ρ ( x , t ) sign ( σ ) + σ T GB ( x , t ) Δ f m ( x , t ) = σ T ρ ( x , t ) sign ( σ ) + GB ( x , t ) Δ f m ( x , t ) ρ ( x , t ) σ 1 + σ · GB ( x , t ) Δ f m ( x , t ) ρ ( x , t ) σ + σ · GB ( x , t ) ρ m ( x , t ) < 0 .
Therefore, u d ( t ) in (14) can also guarantee the enforcement of the sliding mode with respect to system (9).

3.2. Compensation Term

However, u d ( t ) in (10) and (14) cannot guarantee the enforcement of the sliding mode with respect to system (1), since its control effect will be refracted by Δ B ( x , t ) . To be specific, considering Lyapunov function (12), differentiating it with respect to time along with the system (1) and general control law u d ( t ) (10) (or (14)), the following inequality is obtained:
V ˙ σ = σ T Gf σ ( x , t ) + GB ( x , t ) I + Δ B u d ( t ) + Δ f m ( x , t ) = σ T Gf σ ( x , t ) + GB ( x , t ) u d ( t ) + GB ( x , t ) Δ B u d ( t ) + Δ f m ( x , t ) < σ T GB ( x , t ) Δ B u d ( t ) = v T Δ B ( x , t ) u d ( t ) .
It is evident that a residual term v T Δ B ( x , t ) u d ( t ) appears in the last line of inequality V ˙ σ and destroys the negative definiteness of V ˙ σ . To address this issue, the compensation term u c ( t ) is needed.
To design the compensation term u c ( t ) , differentiating the Lyapunov function V ˙ σ in (12) with respect to time based on (1) and (8), one obtains
V ˙ σ = σ T Gf σ ( x , t ) + GB ( x , t ) I + Δ B u d ( t ) + u c ( t ) + Δ f m ( x , t ) < v T Δ B ( x , t ) u d ( t ) + v T ( I + Δ B ( x , t ) ) u c ( t ) = v T u c ( t ) + v T Δ B ( x , t ) u d ( t ) + u c ( t ) .
Let
ψ = v T u c ( t ) + v T Δ B ( x , t ) u d ( t ) + u c ( t ) .
It is apparent that if ψ is nonpositive, the negative definiteness of V ˙ σ is guaranteed. According to Assumption 3, ψ can be further derived as follows:
ψ = v T u c + v T Δ B u v T u c + i = 1 m [ v ] i j = 1 m [ Δ B ] i , j [ u ] j v T u c + i = 1 m [ v ] i sign [ v ] i j = 1 m [ Δ B ¯ ] i , j sign [ u ] j [ u ] j = v T u c + v T S ( v ) Δ B ¯ S ( u ) u = v T u c + S ( v ) Δ B ¯ S ( u ) u .
Let
u c + S ( v ) Δ B ¯ S ( u ) u = 0 ,
then ψ is nonpositive. To be clear, change (20) into the following form
S ( u ) + S ( v ) Δ B ¯ u = u d .
Since u d has been designed in the first step, the only unknown in (21) is u (if u is determined, u c is determined). Then, the key problems are the following: (1) Whether there exists solution u to Equation (21)? (2) If the solution exists, is it unique?
Before handling the above two problems, change (21) into a the following compact form:
S + H p = u d ,
where S = S ( u ) , p = u 0 . According to Assumption 3, H = S ( v ) Δ B ¯ and satisfies
H = S ( v ) Δ B ¯ Δ B ¯ < 1 .
Based on the definition of S, there is a total of 2 m different selections of S. Let the i-th selection of S be S i , i [ 2 m ] and let D = S + H ; then, there exists 2 m corresponding to D and the i-th selection of D can be denoted by D i = S i + H , i [ 2 m ] . Thus, the following definitions are given
Definition 1. 
Regarding D i , define the following closed convex cone of real vector space R m as
C i = D i p p 0 .
Correspondingly, define the interior open region of C i as
C i o = D i p p > 0 ,
and define the conical surface of C i as
C ^ i = C i C i o .
The objective of the following text is to solve S p from (22) which is equivalent to finding the unique solution u to Equation (21).
Based on the definition of the convex cone, the following Theorem 1 analyzes the existence and uniqueness of the solution to Equation (22). In order to support the proof of Theorem 1, three lemmas describing the properties of the convex cone C i , the interior open region C i o and the conical surface C ^ i are first proposed.
Lemma 1. 
If H < 1 , C i o C j o = Ø for all i j , i , j 2 m .
Proof. 
Proof of contradiction is adopted here. Assume that an m-dimensional vector y belongs to both C i o and C j o , i j , there exist two positive p s, denoted by p a and p b , such that y = D i p a = D j p b . Then
S i p a S j p b = S ( v ) Δ B ¯ p b p a .
If p a = p b , S i = S j can be derived from (26), then
D i = S i + S ( v ) Δ B ¯ = S j + S ( v ) Δ B ¯ = D j .
C i o and C j o are the same cone.
If p a p b , the following inequality holds:
S i p a S j p b S ( v ) Δ B ¯ p b p a < p b p a .
However, the following inequality also holds:
S i p a S j p b p b p a ,
since p a and p b are positive. Then, it is easy to obtain that
S i p a S j p b p b p a ,
which derives the contradiction with the previous inequality (27). This contradiction indicates that y cannot possibly simultaneously belong to C i o and C j o if i j . Consequently, Lemma 1 is proved. □
Lemma 2. 
Any closed convex cone C i does not have an independent conical surface; namely, if an m-dimensional vector y C ^ i , there exists another conical surface C ^ j , i j that y C ^ j .
Proof. 
Suppose an m-dimensional vector y C ^ i . Then, y = D i p and there is at least one element of p being zero. Therefore, there exists S i and S j , i j so that S i p = S j p . It follows that y C ^ j . Consequently, the conical surface C ^ i is not independent. Lemma 2 is proved. □
Lemma 3. 
If H < 1 , C = i = 1 2 m C i is equal to R m .
Proof. 
C = i = 1 2 m C i is the union of finite closed convex cones; thus, C is a closed cone. Correspondingly, the conical surface of C should be the union of the independent surfaces of C i , i 2 m ; however, there is no independent conical surface for any one C i according to Lemma 2; therefore, the surface of C is an empty set. Since R m itself is a closed cone with no conical surface, uniquely, C = R m . □
Based on the above three lemmas, the following theorem is proposed in order to ensure the existence and uniqueness of the solution to Equation (22).
Theorem 1. 
For arbitrarily given vector u d R m , the solution S p of (22) exists and is unique if H < 1 .
Proof. 
If u d = 0 , since H < 1 , D = S + H in (22) is nonsingular, the unique solution to the equation D p = 0 exists, which is p = 0 .
If u d 0 , the following proof will be divided into two parts: proof of existence and proof of uniqueness. □
Proof of Existence. 
According to Lemma 3, any u d R m belongs to at least one closed convex cone. Suppose a u d belongs to the convex cone C i , then u d = D i p holds with p 0 . Hence, the solution S i p exists. □
Proof of Uniqueness. 
If u d C i o , u d C j o , j i , according to Lemma 1. Then, u d can be uniquely expressed as u d = D i p = S i + H p with p > 0 and u = S i p is unique.
If u d C ^ i , at least one element in p is zero according to the definition of conical surface. Based on Assumption 2, there exist at least two selections of S, denoted by S i and S j , i j , such that S i p = S j p . It follows that S p is unique. □
According to Theorem 1, there exists a unique solution to Equation (22) if H < 1 . Namely, if S ( v ) Δ B ¯ < 1 , there exists a solution u to Equation (21) and this solution is unique. Subsequently, the key point is to solve the solution from (21). Algorithm 1 presents an algorithmic process, in which u can be obtained by the traversal of the selections of S. Usually, the low dimension of the control input u speeds up the traversal process. Additionally, the compensation control term u c can be obtained by u u d .
Algorithm 1 Calculation of the control law u
Input: m, Δ B ¯
Output:  u ;
  1:
Design u d and calculate H;
  2:
i = 1 ;
  3:
f l a g = T r u e ;
  4:
while  f l a g   do
  5:
    Calculate p = S i + H 1 u d ;
  6:
    if sign( p > 0 and i 2 m  then
  7:
         f l a g = F a l s e ;
  8:
         u = S i p ;
  9:
    else
10:
         i = i + 1 ;
11:
    end if
12:
end while
13:
return  u ;
Remark 1. 
In Assumption 3, Δ B ¯ is bounded in Euclidean norm. In fact, Δ B ¯ can also be bounded in L 1 norm and L norm, which are all suitable for Theorem 1.
Remark 2. 
In the proposed method, u d and u c are designed in separate two cascade steps. Such a design strategy reduces the difficulties of the control design for uncertain system (1). Furthermore, u d and u c are designed based on the Lyapunov function and Lyapunov stability theorem, so u d can also be obtained by any other control method following the Lyapunov stability theorem, such as the Backstepping method, the LMI robust method, etc, not limited to SMC.

4. Stability Analysis

The stability of the system in sliding mode is discussed under two scenarios.
  • If system (1) can be transformed into the regular canonical form in (3), linear sliding surface is defined as (4). During the sliding mode,
    x 2 = G 1 x 1 .
    Taking (4) into the subsystem of (3), i.e., x ˙ 1 = f 1 ( x 1 , x 2 ) , one obtains
    x ˙ 1 = f 1 ( x 1 , G 1 x 1 ) .
    Since G 1 is a diagonal positive definite matrix, x 2 = G 1 · 0 = 0 . Supposing that the origin of system (31), i.e., x 1 = 0 is globally asymptotically stable, if not considering Δ B ( x , t ) , there exists a gain K ( · ) providing a global asymptotic stabilization of the origin with respect to (3) by means of the control [24,25,30]:
    u r ( t ) = GB x , t 1 Gf x , t K sign ( σ )
    One can observe that u d ( t ) in (14) has the same form as u r ( t ) in (32) and K in (32) is equal to GB x , t 1 ρ ( x , t ) in (14). Therefore, u d ( t ) in (14), which ensures the enforcement of sliding mode, is able to provide a global asymptotic stabilization of the origin with respect to (3) without Δ B ( x , t ) [25].
    Since Δ B ( x , t ) also satisfies the matched condition and the compensation control term proposed in this paper can suppress its effect and guarantee the system reachability, in conclusion, the global asymptotic stabilization of system (3) can be achieved by combining the general control u d in (14) and the compensation term u c proposed in this paper.
  • If integral sliding surface (5) is adopted here, considering Assumption 3, as well as the non-singularity of GB ( x , t ) , the equivalent control can be solved from equation σ ˙ = 0 [10]:
    u e q ( t ) = I + Δ B ( x , t ) 1 κ ( x , t ) Δ f m ( x , t ) .
    Then, the expression for the sliding motion can be obtained by substituting the value of u e q ( t ) from (1) into (33), which yields
    x ˙ ( t ) = f ( x , t ) + B ( x , t ) κ ( x , t ) ,
    One can find the origin of system (1) in the sliding mode is globally asymptotically stable [10], since κ ( x , t ) can globally asymptotically stabilize system (34).

5. Numerical Simulation Results

To evaluate the effectiveness of the proposed method, we consider two systems with different structures: a nonlinear plant and a two-link robot manipulator. An SMC proposed by Cao and Xu [10] and another SMC proposed by Feng et al. [21] are used to compare the proposed method. In the following text, we denote Cao and Xu’s method cSMC, denote the method proposed by Feng et al. fSMC and denote the proposed method pSMC.

5.1. Nonlinear Plant with Matched Uncertainties

Consider a nonlinear plant with the following parameters:
f = 1 2 3 4 5 6 7 8 9 x + 0 0 0.2 x 3 sin 2 x 1 ,
B = 1 2 3 + 0.1 cos 3 x 2 4 5 6 + 0.1 cos x 3 2 ,
Δ f m = 0.1 x 1 2 + 0.1 0.1 x 1 2 + 0.2 x 2 3 + 0.1 x 3 4 + 0.1 ,
Δ B = ρ b sin x 2 0.15 cos x 3 0.01 sin 3 x 1 cos 2 x 1 .
where x = x 1 x 2 x 3 T R 3 is the state vector. The initial value x 0 = 1 1 1 T . The control variable u = u 1 u 2 T R 2 . ρ m is the upper bound of the Euclidean norm of Δ f m .
The simulation is performed based on the actual control system. The control signal is calculated by a digital controller based on the current states at each sampling time and is kept unchanged until the next sampling time. Since the parameters of the controlled system are all continuous functions of time, the Runge–Kutta (RK4) method is applied in the simulation. The sampling interval is T s ; the trial duration is T.
Since it is difficult to find a diffeomorphism to transfer the system into the regular form, integral sliding surface is utilized here. Select the integral sliding surface as
σ = G x G x ( 0 ) 0 t ( G f + G B κ ) d τ = 0 ,
where σ = σ 1 σ 2 T R 2 is the sliding variables. G is constant matrix; GB is non-singular. κ = Kx and K is constant matrix. The system parameter values are listed in Table 1. These settings are cited from [10] for a fair comparison. Based on the stability analysis in [10], the value of K is able to ensure the global asymptotical stability of system x ˙ = f + B κ .
The first simulation is implemented to demonstrate the influence of input matrix uncertainty on the system. Figure 1 shows the state variable performances in the cases of ρ b = 0.1 and ρ b = 0.9 by using cSMC [10]. The control law of cSMC is formulated as [10]
u ( t ) = κ ρ v 0 ,
where
ρ > 1 1 ρ b ρ b κ + ρ m .
Based on the second case (integral sliding surface) in Section 4, the control law of cSMC can ensure the uncertain system is globally asymptotically stable.
It can be found that three state variables x 1 , x 2 and x 3 show obvious oscillations with large amplitude when ρ b = 0.9 , versus much slighter oscillations when ρ b = 0.1 . It is more clear from the insets over 2 s to 2.09 s that the amplitude and frequency of the oscillations in state variables are higher when ρ b = 0.9 . Similarly, in Figure 2, the magnitude of control variables u 1 and u 2 are much larger when ρ b = 0.9 than those when ρ b = 0.1 . The reason is that a larger input matrix uncertainty brings a more severe disturbance to the control. Therefore, a more powerful control with larger magnitude of control is needed to suppress the disturbance, which will necessarily increase the chattering amplitude.
The second simulation is conducted to compare the effects of cSMC, fSMC and pSMC under ρ b = 0.9 . The control law of cSMC is formulated as (36). The control law of fSMC is formulated as [21]
u d = κ ρ v 0 u c = β 2 μ 2 + β β μ u d v 0 ,
where ρ > ρ m , β = ρ b 2 1 ρ b 2 , μ = v 0 T u d ( t ) / u d ( t ) , u d ( t ) > 0 0 , u d ( t ) = 0 , .
For pSMC, u d is the same as u d in (37). Then, the control law of pSMC can be obtained according to the algorithmic process in Algorithm 1, where the upper bound matrix Δ B ¯ = ρ b 0.15 ρ b 0.01 ρ b ρ b and Δ B ¯ = 0.9742 < 1 , which satisfies Assumption 3. Based on the second case (integral sliding surface) in Section 4, all three methods can ensure that the uncertain system is globally asymptotically stable.
Figure 3 shows the performances of control variables u 1 and u 2 under cSMC, fSMC and pSMC. One can find that the oscillations of both u 1 and u 2 under pSMC are obviously the slightest among the three methods, versus fSMC and cSMC with much more severe oscillations existing in the trajectories. From the insets over 14 s to 14.09 s, it is evident that pSMC outperforms fSMC and cSMC in reducing both the magnitude and frequency of the control impulse. Such an impulse can be viewed as an adverse impact on the system, which leads to chattering. To measure the impulse and statistically analyze the impact, an index termed Mean Jump Value (MJV) is introduced. Assume that the value of u i , i = 1 , 2 at the j-th sampling point is u i ( j ) , j = 1 , 2 , , n = T / T s , the MJV across the interval T a to T b is defined as ϕ i ( T a , T b ) and computed according to
ϕ i ( T a , T b ) = 1 ( T b T a ) / T s 1 j = T a / T s T b / T s 1 u i ( j + 1 ) u i ( j ) ,
where 0 T a < T b T , T a and T b can be divisible by T s . A small MJV indicates small impulse of the control, which causes lighter chattering.
The MJVs of the control variables in the first five seconds ( T a = 0 s, T b = 5 s) and for the entire time ( T a = 0 s, T b = 15 s) with two methods are listed in Table 2.
Table 2 presents the MJV comparisons among cSMC, fSMC and pSMC in the first five seconds, the last five seconds and the entire trial time. A clear advantage of pSMC is shown: during the first five seconds, the MJVs under cSMC are more than 8 times and 20 times larger than those under pSMC for u 1 and u 2 , respectively. The MJVs under fSMC are more than 2 times and 3 times larger than those under pSMC, respectively. During the last five seconds, the MJVs under cSMC and fSMC are approximately three times larger than those under pSMC. Throughout the entire trial time, the MJVs for u 1 and u 2 under cSMC are greater than one-fourth and one-sixth proportional to those under pSMC, respectively. The MJVs under fSMC are approximately three times larger than those under pSMC. Figure 4a,b depict the comparisons on sliding variables σ 1 and σ 2 , respectively, which illustrates the smaller magnitude of trajectories under pSMC in comparison to those under cSMC and fSMC. Figure 5 presents the performances of the state variables under all three methods. One can observe that chattering under cSMC and fSMC are obviously more severe, particularly for state variables x 2 in Figure 5b and x 3 in Figure 5c. A clearer result can be demonstrated in the insets, which cover the time period from 14 s to 14.09 s.
Finally, the boundary layer solution is applied to v 0 in the discontinuous control of cSMC, fSMC and pSMC as follows:
v 0 = v 0 , v > ε v ε , v ε ,
where ε = 0.011 .
Figure 6 depicts the comparisons of the control variables among the three methods. It is evident that under cSMC, there is no oscillation present in the lines of variables, whereas under cSMC and fSMC, oscillations are present in the entire trajectories. Figure 7 presents the trajectories of the state variables and confirms the superiority of pSMC: there is no chattering in the state trajectories under pSMC, whereas chattering is present in every state variable under cSMC and fSMC. Based on these results, In conclusion, pSMC outperforms cSMC and fSMC in reducing chattering.

5.2. A Two-Link Robot Manipulator

The dynamic equation of the two-link robot manipulator shown in Figure 8 is given by [9,24] and formulated as
M ( q ) q ¨ + C ( q , q ˙ ) q ˙ + G ( q ) = τ ,
where M ( q ) is the moment matrix of inertia, C ( q , q ˙ ) q ˙ is the Coriolis centripetal torque vector and G ( q ) is the the gravitational torque vector. q = q 1 q 2 T represents the angular positions; τ is the vector of applied torques.
For convenience, s i = sin ( q i ) , c i = cos ( q i ) for i = 1 , 2 . q d = q 1 d q 2 d T is the target angular positions. The control objective is to control the robot angular positions to track the given target trajectories. Define the vectors of state, target and input as
x = q 1 q 2 q ˙ 1 q ˙ 2 T ,
x d = q 1 d q 2 d q ˙ 1 d q ˙ 2 d T ,
u = τ .
Then, the dynamic equation of the two-link manipulator in state-space form can be written as
x ˙ = f ( x ) + B ( x ) u ,
where
1 = [ x ] 3 , [ f ( x ) ] 2 = [ x ] 4 , [ f ( x ) ] 3 = s 1 c 3 c 1 s 3 × m 2 l 1 l 2 s 1 s 3 + c 1 c 3 x 2 2 m 2 l 2 2 x 4 2 l 1 l 2 m 1 + m 2 m 2 s 1 s 3 + c 1 c 3 2 + m 1 + m 2 l 2 g s 1 m 2 l 2 g s 3 s 1 s 3 + c 1 c 3 l 1 l 2 m 1 + m 2 m 2 s 1 s 3 + c 1 c 3 2 ,
4 = s 1 c 3 c 1 s 3 × m 1 + m 2 l 1 x 2 2 + m 2 l 1 l 2 s 1 s 3 + c 1 c 3 x 4 2 l 1 l 2 m 1 + m 2 m 2 s 1 s 3 + c 1 c 3 2 + m 1 + m 2 l 1 g s 1 s 1 s 3 + c 1 c 3 + m 1 + m 2 l 1 g s 3 l 1 l 2 m 1 + m 2 m 2 s 1 s 3 + c 1 c 3 2 ,
[ B ( x ) ] 1 , 1 = [ B ( x ) ] 1 , 2 = [ B ( x ) ] 2 , 1 = [ B ( x ) ] 2 , 2 = 0 ,
[ B ( x ) ] 3 , 1 = m 2 l 2 2 m 2 l 1 2 l 2 2 m 1 + m 2 m 2 s 1 s 3 + c 1 c 3 2 ,
[ B ( x ) ] 3 , 2 = m 2 l 1 l 2 s 1 s 3 + c 1 c 3 m 2 l 1 2 l 2 2 m 1 + m 2 m 2 s 1 s 3 + c 1 c 3 2 ,
[ B ( x ) ] 4 , 1 = m 2 l 1 l 2 s 1 s 3 + c 1 c 3 m 2 l 1 2 l 2 2 m 1 + m 2 m 2 s 1 s 3 + c 1 c 3 2 ,
[ B ( x ) ] 4 , 2 = m 1 + m 2 l 1 2 m 2 l 1 2 l 2 2 m 1 + m 2 m 2 s 1 s 3 + c 1 c 3 2 .
where m 1 and m 2 are masses of the links and l 1 and l 2 are lengths of the links.
In fact, uncertainties exist in the system. The additive matched uncertainty and input matrix uncertainty are
Δ f m = 3 [ f ] 3 [ f ] 4 T + 0.015 0.015 T ,
Δ B = ρ b s 1 c 1 2 0 s 3 2 .
The upper bound of Δ f m is ρ m = 3 [ f ] 3 + [ f ] 4 + 0.03 , the upper bound matrix Δ B ¯ = ρ b 1 1 0 1 . With consideration of uncertainties, system (43) can be written as
x ˙ = f + B I + Δ B u + Δ f m .
Since system (44) can be transformed into the following regular form [25,26,30]:
x ˙ 1 = x 2 , x ˙ 2 = f 2 ( x 1 , x 2 ) + B 2 ( x 1 , x 2 ) I + Δ B ( x 1 , x 2 ) u + Δ f m ( x 1 , x 2 )
where x 1 = q 1 q 2 T , x 2 = q ˙ 1 q ˙ 2 T , f 2 ( x 1 , x 2 ) = [ f ] 2 [ f ] 4 T and B 2 ( x 1 , x 2 ) = [ B ] 2 , 1 [ B ] 2 , 2 [ B ] 4 , 1 [ B ] 4 , 2 T . Then, define the linear sliding surface as
σ = x 2 + G 1 x 1 = Gx .
The initial state is x 0 = 0.01 0.01 0 0 T . The target state is x d = 0 0 0 0 T .
The simulation is performed based on the actual control system. The control signal is calculated by a digital controller based on the current states at each sampling time and is kept unchanged until the next sampling time. Since the parameters of the controlled system are all continuous functions of time, the Runge–Kutta (RK4) method is applied in the simulation. The sampling interval is T s ; the trial duration is T. The parameters’ values are listed in Table 3.
In the first simulation, we attempt to investigate the influence of input matrix uncertainty on the system. Under cSMC [10], the cases of ρ b = 0 and ρ b = 0.607 are compared. For ρ b = 0.607 , Δ B ¯ ρ b 1 1 0 1 = 0.9821 < 1 , which satisfies Assumption 3. The control law of cSMC is
u d = GB 1 Gf + ρ sign ( σ ) u c = ρ b 1 ρ b u d ( t ) v 0
where
ρ > GB ρ m .
Based on the first case in Section 4, cSMC is able to ensure the global asymptotical stability of the system with uncertainties.
Figure 9 depicts the performances of state variables q 1 and q 2 in the case of ρ b = 0.607 and ρ b = 0 . One can observe that obvious chattering and severe overshoots exist before 2 s in the state trajectories when ρ b = 0.607 . The states converge to zero without chattering and overshoot when ρ b = 0 . Figure 10 plots the control performances when ρ b = 0.607 and ρ b = 0 . It is evident that before 2 s, the magnitude and oscillations of the control variables when ρ b = 0.607 are larger than those when ρ b = 0 . In conclusion, input matrix uncertainty has adverse effects on systems and leads to the large magnitude of control gain, which aggravates the chattering effect.
The second simulation is conducted to investigate the performance of pSMC by comparing with cSMC and fSMC. ρ b = 0.6 and Δ B ¯ ρ b 1 1 0 1 = 0.9708 < 1 , which also satisfies Assumption 3.
The control law of fSMC is formulated as
u d = GB 1 Gf + ρ sign ( σ ) u c = β 2 μ 2 + β β μ u d v 0 ,
where ρ > GB ρ m , β = ρ b 2 1 ρ b 2 , μ = v 0 T u d / u d , u d > 0 0 , u d = 0 , .
For pSMC, u d is the same as u d in cSMC and fSMC; u can be calculated according to the process Algorithm 1. Based on the first case in Section 4, all three methods can ensure the global asymptotical stability of the system.
Figure 11 shows the comparisons of the system states q 1 and q 2 among cSMC, fSMC and pSMC. One can observe that compared to cSMC and fSMC, pSMC demonstrates superior performance in terms of convergence level, overshoot reduction and chattering suppression. To be specific, for cSMC, both q 1 and q 2 exhibit severe chattering and large overshoots before 2 s, whereas no chattering and overshoots exist in the state trajectories under fSMC and pSMC. The insets over the period from 5 s to 5.09 s visualize that pSMC outperforms fSMC and cSMC in reducing the steady-state errors.
Figure 12 shows the comparisons of the control variables u 1 and u 2 under cSMC, fSMC and pSMC. As illustrated from the original images and the insets over 5 s to 5.09 s, the oscillations in the variable trajectories under pSMC are featured with smaller magnitude and notably lower frequency in comparison with the oscillations under cSMC and fSMC. The consistent results are reported in Table 4. It is shown that the MJVs of u 1 and u 2 under cSMC are more than 12 times and 55 times greater than those under pSMC, respectively, and the MJVs under fSMC are more than 6 times and 26 times larger than those under pSMC, respectively. Similarly, in Figure 13, larger magnitude and higher frequency of the oscillations are shown in the performances of sliding variables σ 1 and σ 2 under cSMC and fSMC than those under pSMC. Consequently, pSMC is more effective than cSMC and fSMC at mitigating the steady-state errors, while simultaneously reducing the magnitude and frequency of oscillations in the discontinuous control, which aids in the suppression of chattering.
In the third simulation, the saturation function sat ( σ / ε ) is utilized to replace the discontinuous controls in cSMC, fSMC and pSMC for continuous approximation, where ε = 0.0005 .
Figure 14 and Figure 15 depict the performance comparisons of state variables and control variables among three control methods. As shown in Figure 15, only under pSMC, trajectories do not exhibit oscillations, whereas obvious oscillations appear in the entire lines under cSMC and fSMC. Figure 14 demonstrates that the continuous approximation in control reduces chattering in all three methods. However, cSMC still exhibits some minor chattering in the state trajectories and the steady-state error under pSMC remains the smallest of the three approaches. Based on the above analysis, it can be concluded that pSMC is able to eliminate chattering with smaller ε .
Finally, simulation is implemented in the case where the initial point x 0 is selected far from the target state x d = 0 . Choose x 0 = 0.1 0.1 0 0 T . The saturation function sat ( σ / ε ) is applied to three control methods for continuous approximation, where ε = 0.0005 .
Figure 16 presents the performances of state variables q 1 and q 2 under cSMC, fSMC and pSMC. The comparisons indicates that pSMC maintains the same high quality performances, including good convergence to the desired states and the absence of chattering. However, larger steady-state errors are present in the trajectories under fSMC and chattering exists in the trajectories before 2 s under cSMC. Figure 17 illustrates the comparisons of control variables among three methods. Evidently, with the exception of a slight oscillation at the very beginning of the line in u 1 , there is no oscillation present in the control variables under pSMC. Both cSMC and fSMC display continuous oscillations throughout the entire trial time. In conclusion, pSMC is effective in improving convergence, reducing overshoot and suppressing chattering.

6. Conclusions

In this paper, a novel SMC method is proposed for affine systems with input matrix uncertainty. To facilitate the controller design, the proposed method builds on a two-step strategy. In the first step, a general control law is adopted for the system without considering input matrix uncertainty; in the second step, a novel control term is designed to specially address input matrix uncertainty. This term is able to compensate for every element of the uncertain input matrix and can be obtained by solving a nonlinear vector equation. A theorem and three lemmas are proposed to guarantee the existence and the uniqueness of the equation solution and an algorithmic process is proposed to solve the vector equation. Numerical simulations on two affine uncertain systems illustrate that the proposed method can provide desired closed-loop system performances in the presence of input uncertainties. Furthermore, the simulation results revealed that the proposed method (pSMC) outperforms the SMC approach proposed by Cao and Xu (cSMC) and by Feng et al. (fSMC) in terms of improving convergence, reducing overshoot and alleviating the chattering problem.

Author Contributions

Methodology, Z.F.; Formal analysis, R.S.; Writing—original draft, R.S.; Writing—review & editing, R.S.; Supervision, Z.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the Project of Precision Swing Table with Project Number 32199000000154.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ji, W.; Qiu, J.; Karimi, H.R.; Fu, Y. New Results on Fuzzy Integral Sliding Mode Control of Nonlinear Singularly Perturbed Systems. IEEE Trans. Fuzzy Syst. 2021, 29, 2062–2067. [Google Scholar] [CrossRef]
  2. Ji, W.; Ji, W.; Chadli, M. A Novel Fuzzy Output Feedback Dynamic Sliding Mode Controller Design for Two-Dimensional Nonlinear Systems. IEEE Trans. Fuzzy Syst. 2021, 29, 2869–2877. [Google Scholar] [CrossRef]
  3. Chen, W.-H.; Deng, X.; Zheng, W.X. Sliding-Mode Control for Linear Uncertain Systems with Impulse Effects via Switching Gains. IEEE Trans. Autom. Control 2022, 67, 2044–2051. [Google Scholar] [CrossRef]
  4. Kao, Y.; Ma, S.; Xia, H.; Wang, C.; Liu, Y. Integral Sliding Mode Control for a Kind of Impulsive Uncertain Reaction-Diffusion Systems. IEEE Trans. Autom. Control 2022, 68, 1154–1160. [Google Scholar] [CrossRef]
  5. Utkin, V.I. Variable Structure Systems with Sliding Modes. IEEE Trans. Autom. Control 1977, 22, 212–222. [Google Scholar] [CrossRef]
  6. Utkin, V.I. Sliding Mode Control Design Principles and Applications to Electric Drives. IEEE Trans. Ind. Electron. 1993, 40, 23–36. [Google Scholar] [CrossRef] [Green Version]
  7. Sharma Kumar, N.; Janardhanan, S. Discrete-Time Higher-Order Sliding Mode Control of Systems with Unmatched Uncertainty. Int. J. Robust Nonlinear Control 2019, 29, 135–152. [Google Scholar] [CrossRef] [Green Version]
  8. Bejarano-Rincón, A.; Estrada, A.; Herrera-Hernández, E.C.; Orozco, J.M.A. Control Design for a Class of Multivariable Nonlinear System with Uncertain Control Direction: A Laser Cladding Case Study. Eur. J. Control 2021, 60, 114–124. [Google Scholar] [CrossRef]
  9. Nasiri, A.; Nguang, S.K.; Swain, A. Adaptive Sliding Mode Control for a Class of MIMO Nonlinear Systems with Uncertainties. J. Frankl. Inst. 2014, 351, 2048–2061. [Google Scholar] [CrossRef]
  10. Cao, W.-J.; Xu, J. Nonlinear Integral-Type Sliding Surface for Both Matched and Unmatched Uncertain Systems. IEEE Trans. Autom. Control 2004, 49, 1355–1360. [Google Scholar] [CrossRef]
  11. Chen, L.; Edwards, C.; Alwi, H. Integral Sliding Mode Fault-Tolerant Control Allocation for a Class of Affine Nonlinear System. Int. J. Robust Nonlinear Control 2019, 29, 565–582. [Google Scholar] [CrossRef]
  12. Chang, J.-L. Dynamic Output Integral Sliding-Mode Control with Disturbance Attenuation. IEEE Trans. Autom. Control 2009, 54, 2653–2658. [Google Scholar] [CrossRef]
  13. Shtessel, Y.; Edwards, C.; Fridman, L.; Levant, A. Sliding Mode Control and Observation. In Control Engineering; Springer: New York, NY, USA, 2014. [Google Scholar] [CrossRef]
  14. Zhang, X.; Xiao, L.; Li, H. Robust Control for Switched Systems with Unmatched Uncertainties Based on Switched Robust Integral Sliding Mode. IEEE Access 2020, 8, 138396–138405. [Google Scholar] [CrossRef]
  15. Lin, C.-H.; Hu, G.-H.; Yan, J.-J. Estimation of Synchronization Errors between Master and Slave Chaotic Systems with Matched/Mismatched Disturbances and Input Uncertainty. Mathematics 2021, 9, 176. [Google Scholar] [CrossRef]
  16. Lee, H.; Utkin, V.I. Chattering Suppression Methods in Sliding Mode Control Systems. Annu. Rev. Control 2007, 31, 179–188. [Google Scholar] [CrossRef]
  17. Utkin, V.I.; Lee, H. Chattering Problem in Sliding Mode Control Systems. In Proceedings of the International Workshop on Variable Structure Systems, Alghero, Sardinia, 5–7 June 2006; pp. 346–350. [Google Scholar] [CrossRef]
  18. Castillo, I.; Fridman, L.; Moreno, J.A. Super-Twisting Algorithm in Presence of Time and State Dependent Perturbations. Int. J. Control 2018, 91, 2535–2548. [Google Scholar] [CrossRef]
  19. Martínez-Fuentes, A.C.; Seeber, R.; Fridman, L.; Moreno, J.A. Saturated Lipschitz Continuous Sliding Mode Controller for Perturbed Systems with Uncertain Control Coefficient. IEEE Trans. Autom. Control 2021, 66, 3885–3891. [Google Scholar] [CrossRef]
  20. Gutiérrez-Oribio, D.; Mercado-Uribe, A.; Moreno, J.A.; Fridman, L. Robust Global Stabilization of a Class of Underactuated Mechanical Systems of Two Degrees of Freedom. Int. J. Robust Nonlinear Control 2021, 31, 3908–3928. [Google Scholar] [CrossRef]
  21. Feng, Z.-R.; Sha, R.-Z.; Ren, Z.-G. A Chattering-Reduction Sliding Mode Control Algorithm for Affine Systems with Input Matrix Uncertainty. IEEE Access 2022, 10, 58982–58996. [Google Scholar] [CrossRef]
  22. Hui, L.; Liu, L.; Wang, Y.; Cheng, Z.; Luo, Q. Active Fault-Tolerant Attitude Tracking Control with Adaptive Gain for Spacecrafts. Aerosp. Sci. Technol. 2020, 98, 105706. [Google Scholar] [CrossRef]
  23. Utkin, V.I.; Shi, J. Integral Sliding Mode in Systems Operating under Uncertainty Conditions. In Proceedings of the 35th IEEE Conference on Decision and Control, Kobe, Japan, 13 December 1996; Volume 4, pp. 4591–4596. [Google Scholar] [CrossRef]
  24. Sira-Ramírez, H. Differential Geometric Methods in Variable-Structure Control. Int. J. Control 1988, 48, 1359–1390. [Google Scholar] [CrossRef]
  25. Perruquetti, W.; Barbot, J.-P. Sliding Mode Control in Engineering; CRC Press: New York, NY, USA, 2002. [Google Scholar] [CrossRef]
  26. Utkin, V.I.; Guldner, J.; Shi, J. Sliding Mode Control in Electro-Mechanical Systems, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar] [CrossRef]
  27. Loukianov, A.; Utkin, V.I. Methods of Reducing Equations of Dynamic Systems to Regular Form. Autom. Remote. Control 1981, 42, 413–420. [Google Scholar]
  28. Rubagotti, M.; Estrada, A.; Castanos, F.; Ferrara, A.; Fridman, L. Integral Sliding Mode Control for Nonlinear Systems with Matched and Unmatched Perturbations. IEEE Trans. Autom. Control 2011, 56, 2699–2704. [Google Scholar] [CrossRef]
  29. Hamayun, T.M.; Edwards, C.; Alwi, H. Integral Sliding Mode Control In Fault Tolerant Control Schemes Using Integral Sliding Modes; Hamayun, M.T., Edwards, C., Alwi, H., Eds.; Studies in Systems, Decision and Control; Springer International Publishing: Cham, Switzerland, 2016; pp. 17–37. [Google Scholar] [CrossRef]
  30. Perruquetti, W.; Borne, P.; Richard, J. A Generalized Regular Form for Sliding Mode Stabilization of MIMO Systems. In Proceedings of the 36th IEEE Conference on Decision and Control, San Diego, CA, USA, 13 December 1997; Volume 2, pp. 957–961. [Google Scholar] [CrossRef]
Figure 1. The state variable performance of the system controlled by cSMC in the case of ρ b = 0.9 (gray solid line) and ρ b = 0.1 (red solid line).
Figure 1. The state variable performance of the system controlled by cSMC in the case of ρ b = 0.9 (gray solid line) and ρ b = 0.1 (red solid line).
Mathematics 11 01337 g001
Figure 2. The control variable performance of the system controlled by cSMC in the case of ρ b = 0.9 (gray solid line) and ρ b = 0.1 (red solid line).
Figure 2. The control variable performance of the system controlled by cSMC in the case of ρ b = 0.9 (gray solid line) and ρ b = 0.1 (red solid line).
Mathematics 11 01337 g002
Figure 3. The trajectories of control variables u 1 and u 2 when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 3. The trajectories of control variables u 1 and u 2 when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g003
Figure 4. The trajectories of sliding variables σ 1 and σ 2 when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 4. The trajectories of sliding variables σ 1 and σ 2 when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g004
Figure 5. The trajectories of state variables x 1 , x 2 and x 3 when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 5. The trajectories of state variables x 1 , x 2 and x 3 when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g005
Figure 6. The trajectories of control variables u 1 and u 2 with continuous approximation when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 6. The trajectories of control variables u 1 and u 2 with continuous approximation when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g006
Figure 7. Thetrajectories of state variables x 1 , x 2 and x 3 with continuous approximation when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 7. Thetrajectories of state variables x 1 , x 2 and x 3 with continuous approximation when ρ b = 0.9 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g007
Figure 8. Two-link robot manipulator dynamics.
Figure 8. Two-link robot manipulator dynamics.
Mathematics 11 01337 g008
Figure 9. The state variable performance of the system controlled by cSMC in the cases of ρ b = 0.607 (Gray solid line) and ρ b = 0 (Red solid line).
Figure 9. The state variable performance of the system controlled by cSMC in the cases of ρ b = 0.607 (Gray solid line) and ρ b = 0 (Red solid line).
Mathematics 11 01337 g009
Figure 10. The control variable performance of the system controlled by cSMC in the cases of ρ b = 0.607 (Gray solid line) and ρ b = 0 (Red solid line).
Figure 10. The control variable performance of the system controlled by cSMC in the cases of ρ b = 0.607 (Gray solid line) and ρ b = 0 (Red solid line).
Mathematics 11 01337 g010
Figure 11. The trajectories of state variables q 1 and q 2 when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 11. The trajectories of state variables q 1 and q 2 when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g011
Figure 12. The trajectories of control variables u 1 and u 2 when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 12. The trajectories of control variables u 1 and u 2 when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g012
Figure 13. The trajectories of sliding variables σ 1 and σ 2 when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 13. The trajectories of sliding variables σ 1 and σ 2 when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g013
Figure 14. Thetrajectories of state variables q 1 and q 2 with continuous approximation in control when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 14. Thetrajectories of state variables q 1 and q 2 with continuous approximation in control when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g014
Figure 15. The trajectories of control variables u 1 and u 2 with continuous approximation in control when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 15. The trajectories of control variables u 1 and u 2 with continuous approximation in control when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g015
Figure 16. The trajectories of state variables q 1 and q 2 with a far initial state when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 16. The trajectories of state variables q 1 and q 2 with a far initial state when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g016
Figure 17. The trajectories of control variables u 1 and u 2 with a far initial state in control when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Figure 17. The trajectories of control variables u 1 and u 2 with a far initial state in control when ρ b = 0.6 : cSMC (Blue solid line), fSMC (Yellow solid line) and pSMC (Red solid line).
Mathematics 11 01337 g017
Table 1. Parameter Values.
Table 1. Parameter Values.
ParameterValue
ρ m 0.2 x 1 2 + 0.2 x 2 3 + 0.1 x 3 4 + 0.2
ρ b ρ b < 1
T s 1 (ms)
T15 (s)
G 0 0.1 0.2 0.4 0.1 0.2
K 1.7902 1.8765 2.4429 0.0444 0.2676 0.0344
Table 2. MJVs of cSMC, fSMC and pSMC at time intervals [0 s, 5 s], [10 s, 15 s] and [0 s, 15 s].
Table 2. MJVs of cSMC, fSMC and pSMC at time intervals [0 s, 5 s], [10 s, 15 s] and [0 s, 15 s].
ϕ 1 ( 0 , 5 ) ϕ 2 ( 0 , 5 ) ϕ 1 ( 10 , 15 ) ϕ 2 ( 10 , 15 ) ϕ 1 ( 0 , 15 ) ϕ 2 ( 0 , 15 )
cSMC37.150123.433721.50196.904926.631712.4571
fSMC8.98724.091420.58886.664116.63335.8035
pSMC4.20521.12946.44532.24335.68501.8903
cSMC/pSMC8.834320.74943.33613.07804.68466.5898
fSMC/pSMC2.13723.62273.19442.97072.92583.0701
Table 3. Parameter Values.
Table 3. Parameter Values.
ParameterValue
m 1 1 (kg)
m 2 1 (kg)
l 1 1 (m)
l 2 1 (m)
g9.8 (m/s 2 )
T6 (s)
T s 1 (ms)
G 1 5 0 0 5
G 5 0 1 0 0 5 0 1
Table 4. MJVs of cSMC, fSMC and pSMC.
Table 4. MJVs of cSMC, fSMC and pSMC.
ϕ 1 ( 0 , 6 ) ϕ 2 ( 0 , 6 )
cSMC7.551712.0323
fSMC3.86515.7159
pSMC0.605126.2142
cSMC/pSMC12.481055.1825
fSMC/pSMC6.388126.2142
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

Sha, R.; Feng, Z. A Sliding Mode Control Algorithm with Elementary Compensation for Input Matrix Uncertainty in Affine Systems. Mathematics 2023, 11, 1337. https://doi.org/10.3390/math11061337

AMA Style

Sha R, Feng Z. A Sliding Mode Control Algorithm with Elementary Compensation for Input Matrix Uncertainty in Affine Systems. Mathematics. 2023; 11(6):1337. https://doi.org/10.3390/math11061337

Chicago/Turabian Style

Sha, Ruizhi, and Zuren Feng. 2023. "A Sliding Mode Control Algorithm with Elementary Compensation for Input Matrix Uncertainty in Affine Systems" Mathematics 11, no. 6: 1337. https://doi.org/10.3390/math11061337

APA Style

Sha, R., & Feng, Z. (2023). A Sliding Mode Control Algorithm with Elementary Compensation for Input Matrix Uncertainty in Affine Systems. Mathematics, 11(6), 1337. https://doi.org/10.3390/math11061337

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