Output feedback sliding mode control based on adaptive sliding mode disturbance observer

Abstract

In this paper, an adaptive sliding mode disturbance observer is designed to counteract the disturbance actively. By designing the adaptive laws, the assumptions on the disturbance are relaxed in the proposed observer, its first derivative upper bound is considered to be unknown. Based on the proposed disturbance observer, an output feedback sliding mode controller is constructed for the continuous-time linear systems with unknown external disturbance. The proposed controller incorporates only the system output information and has less chattering of the control input. The feasibility of the proposed strategy is shown by numerical simulations.

Keywords

Output feedback control sliding mode disturbance observer adaptive control dynamic sliding surface

Sliding mode control (SMC) has attracted the attention of many researchers because of its robustness to disturbances and design simplicity.^1–4 In the area of SMC systems, many of the theoretical developments assume that the system state vector is available.^5–7 In most practical situations, it is unrealistic or prohibitively expensive to measure the full state information, and only the output information can be physically measured.^8–10 Generally, there are static and dynamic output feedback SMC strategies to work around this limitation.¹¹ The static output feedback SMC design problem is studied for a delay system.¹² Song et al.¹³ address the static output-feedback sliding mode control (SMC) problem for a class of uncertain control systems. In order to improve the control performance of static output feedback SMC, the dynamic output feedback SMC can be designed with a compensator. The additional integral term is introduced in the sliding surface, which can provide one more degree of freedom.¹⁴ An asynchronous output feedback sliding controller is proposed for a class of Markovian jump systems.¹⁵ In the aforementioned results for output feedback SMC, the disturbance rejection problem is ignored, which leads to the chattering phenomenon in the system.

Recently, the observer-based method is utilized to eliminate disturbance and reduce chattering.^16,17 Su et al.¹⁸ present a disturbance observer where the bounds on disturbance derivative are assumed to be known. A new form of the combined observer-controller is designed to provide estimated data of unknown disturbance and unmeasured states in the control law.¹⁹ Lee²⁰ proposed a composite control technique by combining a nonlinear disturbance observer. A fixed-time observer has been put forward, which possessed a better approximation to external disturbances.²¹ However, in the above works the disturbance observer was considered with the assumption of the known maximum upper bound of disturbance. In some practical applications, it is difficult to acquire the disturbance upper bound directly, especially combined with the sliding mode algorithm, too large upper bound will aggravate the chattering phenomenon of the system, and too small upper bound will lead to the instability of the system.

The sliding mode disturbance observer with adaptive control is proposed to relax the restriction on disturbance. The adaptive control is introduced in the sliding mode disturbance observer that use the adaptive method to ensure the control gain is as small as possible whereas sufficient to eliminate the disturbances or uncertainties.²² Negrete-Chávez and Moreno²³ develop an adaptive second-order SMC to interference disturbances. An adaptive scheme based on equivalent control is presented and it requires the disturbances upper-bound is known.²⁴ Adaptive SMC methodology is structured which does not require a priori bounded uncertainty.^25,26 In order to further improve the observation performance of the observer, finite-time control methods were studied.²⁷ Finite-time convergence of the disturbance approximation error is guaranteed using the designed disturbance observer.²⁸ Compared with asymptotic convergence,²⁹ finite-time convergence can track the true value of the disturbance faster. The following three problems are interesting but challenging in this paper, (1) How to reduce the chattering of the output feedback controller, (2) how to design the disturbance observers with adaptive control to relax the constraints on disturbances, and (3) how to achieve finite-time convergence in the proposed observer.

The above-mentioned issue motivated us to devise a new output feedback sliding mode control based on the adaptive sliding mode disturbance observer strategy. Initially, the dynamic sliding surface for parameters analytically is designed, and an additional integral term is introduced in the sliding surface, which can provide one more degree of freedom. Based on the dynamic sliding mode surface, an output feedback controller is organized with power terms. It is worth mentioning that, the controller only contains the output information of the system and improves the convergence speed of the system state. In order to attenuate chattering in the control input, a novel ASMDO is developed. The proposed ASMDO does not require information about the bound on the disturbances and their derivatives. Compared to some existing observers, the conditions on disturbance bounds are relaxed. The effectiveness of the proposed control approaches is illustrated by an inverted pendulum on a cart model.

Preliminaries

Consider the following linear system with matched external disturbance

\begin{matrix} \overset{\cdot}{x} (t) & = Ax (t) + B [u (t) + d (t)] \\ y (t) & = Cx (t) \end{matrix}

(1)

where $x (t) \in R^{n}$ is the state vector, $u (t) \in R^{m}$ is the control input, $y (t) \in R^{q}$ is the output vector, and $m \leq q < n$ . $d (t)$ is the unknown disturbance. $A \in R^{n \times n}$ , $B \in R^{n \times m}$ , $C \in R^{q \times n}$ are given constant matrices.

When matched uncertainty alone is present, it is sufficient to consider the nominal linear system representation when designing the switching function. The sliding motion depends on the choice of the sliding surface, the precise effect is not readily apparent. Therefore, we need to transform the system into a suitable regular form.

Assume $rank (CB) = rank (B) = m$ , thus, there exists a non-singular matrix $T \in R^{n \times n}$ such that⁸: $TB = {[\begin{matrix} 0 & B_{2} \end{matrix}]}^{T}$ , $C T^{- 1} = [\begin{matrix} 0_{q \times (n - q)} & C_{2} \end{matrix}]$ $(C_{2} \in R^{q \times q})$ . Partition $T$ as follows

T = [\begin{matrix} T_{1} \\ T_{2} \end{matrix}] \begin{matrix} ↕ & n - m \\ ↕ & m \end{matrix}, T_{1} = [\begin{matrix} T_{11} \\ T_{12} \\ T_{13} \end{matrix}] \begin{matrix} ↕ & m_{1} \\ ↕ & m_{2} \\ ↕ & m_{3} \end{matrix},

where $m_{1} \geq 0, m_{2} \geq 0, m_{3} \geq 0,$ and $m_{1} + m_{2} + m_{3} = n - m$ .

Under the coordinate transformation $\bar{x} (t) = Tx (t)$ , system (1) is transformed into the following regular form

\overset{\cdot}{\bar{x}} (t) = \bar{A} \bar{x} (t) + \bar{B} [u (t) + d (t)]

(2)

y (t) = \bar{C} \bar{x} (t)

(3)

where $\bar{x} (t) = [\begin{matrix} {\bar{x}}_{1} (t) \\ {\bar{x}}_{2} (t) \end{matrix}],$ $\bar{A} = TA T^{- 1} = [\begin{matrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{matrix}],$ $\bar{B} = TB = {[\begin{matrix} 0 & B_{2} \end{matrix}]}^{T},$ $\bar{C} = C T^{- 1} = [\begin{matrix} 0_{q \times (n - q)} & C_{2} \end{matrix}],$ ${\bar{x}}_{1} (t) \in R^{n - m},$ ${\bar{x}}_{2} (t) \in R^{m},$ $B_{2} \in R^{m \times m}$ and $C_{2} \in R^{q \times q} .$ Let $C_{1} = [0_{(q - m) \times (n - q)} I_{q - m}],$ substituting $\bar{A}, \bar{B}$ and $\bar{C}$ into (2) and (3), then system (1) is further written as

{\overset{\cdot}{\bar{x}}}_{1} (t) = A_{11} {\bar{x}}_{1} (t) + A_{12} {\bar{x}}_{2} (t)

(4)

{\overset{\cdot}{\bar{x}}}_{2} (t) = A_{21} {\bar{x}}_{1} (t) + A_{22} {\bar{x}}_{2} (t) + B_{2} [u (t) + d (t)]

(5)

y (t) = C_{2} [\begin{matrix} C_{1} {\bar{x}}_{1} (t) \\ {\bar{x}}_{2} (t) \end{matrix}]

(6)

Proposed method

Design of dynamic sliding surface

The dynamic sliding surface is designed as

s (t) = [- w_{1}, I_{m}] C_{2}^{- 1} y (t) + w_{4} ν (t)

(7)

An additional integral term $ν (t)$ is introduced in the sliding surface, which can provide one more degree of freedom. $ν (t)$ is defined as follows

ν (t) = \int_{0}^{t} e^{w_{3} (t - s)} [w_{2}, Φ] C_{2}^{- 1} y (s) ds

(8)

and $w_{1} \in R^{m \times (q - m)},$ $w_{2} \in R^{(n - m) \times (q - m)},$ $w_{3} \in R^{(n - m) \times (n - m)},$ $w_{4} \in R^{m \times (n - m)}$ are the sliding surface parameter matrices, $Φ \in R^{(n - m) \times m}$ is a known matrix. The derivative with respect to (8)

\begin{matrix} \overset{\cdot}{ν} (t) = w_{3} \int_{0}^{t} e^{w_{3} (t - s)} [w_{2}, Φ] C_{2}^{- 1} y (s) ds + [w_{2}, Φ] C_{2}^{- 1} y (t) \\ = w_{3} ν (t) + w_{2} C_{1} {\bar{x}}_{1} (t) + Φ {\bar{x}}_{2} (t) \end{matrix}

(9)

then (7) can be rewritten as

s (t) = {\bar{x}}_{2} (t) - w_{1} C_{1} {\bar{x}}_{1} (t) + w_{4} ν (t) = 0

(10)

According to (10), we have

{\bar{x}}_{2} (t) = w_{1} C_{1} {\bar{x}}_{1} (t) - w_{4} ν (t)

(11)

Substituting (11) into (4) yields the sliding motion equation

{\overset{\cdot}{\bar{x}}}_{1} (t) = (A_{11} + A_{12} w_{1} C_{1}) {\bar{x}}_{1} (t) - A_{12} w_{4} ν (t)

(12)

Define $W_{2} = w_{2} + Φ w_{1}$ and $W_{3} = w_{3} - Φ w_{4}$ , (9) can be rewritten into the following form

\overset{\cdot}{ν} (t) = W_{3} ν (t) + W_{2} C_{1} {\bar{x}}_{1} (t)

(13)

Thus, the new sliding motion is as follows

\overset{\cdot}{X} (t) = A X (t)

(14)

where $A = [\begin{matrix} A_{11} + A_{12} w_{1} C_{1} & - A_{12} w_{4} \\ W_{2} C_{1} & W_{3} \end{matrix}]$ and $X (t) = [\begin{matrix} {\bar{x}}_{1} (t) \\ ν (t) \end{matrix}]$ .

Proposition 1. Define $Q = H_{2}^{T} A H_{1} + (H_{2}^{T} A H_{1})^{T}$ and $F = [\begin{matrix} F_{1} & F_{2}^{T} \\ F_{2} & F_{3} \end{matrix}]$ , $H_{1}$ and $H_{2}$ will be introduced later. The new sliding motion is asymptotically stable on the dynamic sliding surface, if there exist $F_{1} \in R^{(n - m) \times (n - m)}$ , $G_{1} \in R^{(n - m) \times (n - m)}$ are symmetric positive definite matrices and subject to

Q < 0

(15)

[\begin{matrix} F_{1} & I_{n - m} \\ I_{n - m} & G_{1} \end{matrix}] > 0

(16)

Proof. Define $G_{2}$ and $F_{2}$ are nonsingular matrices, and satisfies $G_{2} F_{2} = I - G_{1} F_{1} .$ Coordinate transformation matrices $H_{1}$ and $H_{2}$ are designed as follows

H_{1} = [\begin{matrix} I & G_{1} \\ 0 & G_{2}^{T} \end{matrix}] H_{2} = [\begin{matrix} F_{1} & I \\ F_{2} & 0 \end{matrix}]

where $I$ is identity matrix. Hence, we have $F_{2} = G_{2}^{- 1} (I - G_{1} F_{1}) .$ Define $F_{3} = - F_{2} G_{2}^{- 1} G_{1} = F_{2} (F_{1} - G_{1}^{- 1})^{- 1} F_{2}^{T},$ Applying Schur complement on (16) obtain $F_{1} - G_{1}^{- 1} > 0$ , then

\begin{matrix} F = H_{2} H_{1}^{- 1} = [\begin{matrix} F_{1} & I \\ F_{2} & 0 \end{matrix}] {[\begin{matrix} I & G_{1} \\ 0 & G_{2}^{T} \end{matrix}]}^{- 1} \\ = [\begin{matrix} F_{1} & (I - F_{1} G_{1}) G_{2}^{- T} \\ F_{2} & - F_{2} G_{2}^{- 1} G_{1} \end{matrix}] \\ = [\begin{matrix} F_{1} & F_{2}^{T} \\ F_{2} & F_{3} \end{matrix}] > 0 \end{matrix}

hence exist a positive definite symmetric matrix $F > 0$ such that

\begin{matrix} H_{1}^{- T} {QH}_{1}^{- 1} = A^{T} H_{2} H_{1}^{- 1} + H_{1}^{- T} H_{2}^{T} A \\ = A^{T} F + F A \\ < 0 \end{matrix}

The new sliding motion (14) is asymptotically stable. ■

Design of sliding mode controller

In this subsection, the control method based on output feedback sliding mode with power term is proposed.

Define the generalized inverse matrix of $C_{1}$ is $C_{1}^{-}$ and $ζ_{1} = null (C_{1})^{T}$ satisfying that $ζ_{1} C_{1}^{T} = 0$ , generalized inverse matrix of $ζ_{1}$ is $ζ_{1}^{-}$ , then equation $C_{1}^{-} C_{1} + ζ_{1}^{-} ζ_{1} = I_{n - m}$ holds. Since $A$ is Hurwitz matrix, then exist $ξ > 0$ and $P > 0$ subject to

A^{T} P + P A + 2 ξ P < 0

(17)

Denote $ξ$ as the largest positive scalar such that (17) is feasible. Utilize Cholesky factorization for $P$ to get $P = P_{1} P_{1}^{T} = [\begin{matrix} P_{11} & 0 \\ P_{21} & P_{22} \end{matrix}]$ ${[\begin{matrix} P_{11} & 0 \\ P_{21} & P_{22} \end{matrix}]}^{T}$ .

Remark 1. Inequality (17) holds which can be explained as follows:

Consider the Lyapunov function for the new sliding motion (14): $V = x^{T} Px$

\begin{matrix} \overset{\cdot}{V} = x^{T} PAx + x^{T} A^{T} Px \\ = x^{T} (PA + A^{T} P) x \\ \leq (PA + A^{T} P) x^{2} \end{matrix}

Since $A$ is Hurwitz matrix and $P > 0, ξ > 0$ , hence $\overset{\cdot}{V} < 0$ , $PA + A^{T} P < 0$ . There exists $ξ > 0$ for $(PA + A^{T} P) + 2 ξ P < 0$ holds.

Proposition 2. Define $τ > 0, ϱ > 0, ρ > 0, k_{1} + k_{2} > 1, h_{1} > 1, 1 > h_{2} > 0,$ then the following output feedback sliding mode controller is proposed:

\begin{matrix} u (t) = - B_{2}^{- 1} [τ s (t) + ψ_{1} (t) (k_{1} sg n^{h 1} (s (t)) \\ + k_{2} sg n^{h 2} (s (t))) + w_{4} w_{3} ν (t) \\ - (A_{21} - w_{1} C_{1} A_{11}) ζ_{1}^{-} ζ_{1} P_{11}^{- T} P_{21}^{T} ν (t) \\ + ((A_{21} - w_{1} C_{1} A_{11}) C_{1}^{-} + w_{4} w_{2}) C_{1} {\bar{x}}_{1} (t) \\ + (A_{22} - w_{1} C_{1} A_{12} + w_{4} Φ) {\bar{x}}_{2} (t) - \hat{d}] \end{matrix}

(18)

where $ψ_{1} (t) = ψ_{0} (ϑ_{1} (t) + ϱ)$ and $ψ_{0} = ∥ (A_{21} - w_{1} C_{1} A_{11}) ζ_{1}^{-} ζ_{1} P_{11}^{- T} ∥$ $ϑ_{1} (t)$ can be obtained through ${\overset{\cdot}{ϑ}}_{1} (t) = - ξ ϑ_{1} (t) + ε_{1} ∥ s (t) ∥$ , coefficient $ε_{1} = ∥ P_{1}^{T} \bar{Φ} ∥$ and $\bar{Φ} = [A_{12} Φ]^{T}$ . $\hat{d}$ is the estimation of the disturbance. Then, the system (2)−(3) can reach the sliding surface.

Proof. The first derivative of the sliding mode surface is given as

\begin{matrix} \overset{\cdot}{s} (t) = {\overset{\cdot}{\bar{x}}}_{2} (t) - w_{1} C_{1} {\overset{\cdot}{\bar{x}}}_{1} (t) + w_{4} \overset{\cdot}{ν} (t) \\ = A_{21} {\bar{x}}_{1} (t) + A_{22} {\bar{x}}_{2} (t) + B_{2} [u (t) + d (t)] \\ - w_{1} C_{1} [A_{11} {\bar{x}}_{1} (t) + A_{12} {\bar{x}}_{2} (t)] \\ + w_{4} [W_{3} ν (t) + W_{2} C_{1} {\bar{x}}_{1} (t)] \\ = [(A_{21} - w_{1} C_{1} A_{11}) {\bar{x}}_{1} (t) + w_{4} w_{2} C_{1} {\bar{x}}_{1} (t) \\ + (A_{22} - w_{1} C_{1} A_{12} + w_{4} Φ) {\bar{x}}_{2} (t) + w_{4} w_{3} ν (t) \\ + B_{2} (u (t) + d (t))] \end{matrix}

(19)

In order to prove that $s (t)$ can converge to 0, define the Lyapunov function $V_{s} (t) = \frac{1}{2} s^{T} (t) s (t)$ , its derivative can be determined and substitute (19) into ${\overset{\cdot}{V}}_{s} (t)$ as

\begin{matrix} {\overset{\cdot}{V}}_{s} (t) = s^{T} (t) [(A_{21} - w_{1} C_{1} A_{11}) {\bar{x}}_{1} (t) + w_{4} w_{2} C_{1} {\bar{x}}_{1} (t) \\ + (A_{22} - w_{1} C_{1} A_{12} + w_{4} Φ) {\bar{x}}_{2} (t) + w_{4} w_{3} ν (t) \\ + B_{2} (u (t) + d (t))] \\ = s^{T} (t) [(A_{21} - w_{1} C_{1} A_{11}) ζ_{1}^{-} ζ_{1} P_{11}^{- T} (P_{11}^{T} {\bar{x}}_{1} (t) \\ + P_{21}^{T} ν (t)) - (A_{21} - w_{1} C_{1} A_{11}) ζ_{1}^{-} ζ_{1} P_{11}^{- T} P_{21}^{T} ν (t) \\ + ((A_{21} - w_{1} C_{1} A_{11}) C_{1}^{-} + w_{4} w_{2}) C_{1} {\bar{x}}_{1} (t) \\ + (A_{22} - w_{1} C_{1} A_{12} + w_{4} Φ) {\bar{x}}_{2} (t) \\ + w_{4} w_{3} ν (t) + B_{2} (u (t) + d (t))] \end{matrix}

(20)

Substituting (18) into (20), one obtains

\begin{matrix} {\overset{\cdot}{V}}_{s} (t) \leq - τ ∥ s (t) ∥^{2} - ψ_{1} (t) (k_{1} + k_{2}) ∥ s (t) ∥ \\ + ψ_{0} ∥ s (t) ∥ ∥ P_{11}^{T} {\bar{x}}_{1} (t) + P_{21}^{T} ν (t) ∥ \\ + s^{T} (t) B_{2} (d (t) - k_{d} sgn (B_{2}^{T} s (t))) \\ \leq - τ ∥ s (t) ∥^{2} + χ (t) ∥ s (t) ∥ + (d (t) - \hat{d}) ∥ B_{2}^{T} s (t) ∥ \end{matrix}

(21)

where $χ (t) = ψ_{0} ∥ P_{11}^{T} {\bar{x}}_{1} (t) + P_{21}^{T} ν (t) ∥ - (k_{1} + k_{2}) ψ_{1} (t) .$ Obviously, $- τ ∥ s (t) ∥^{2} < 0$ , and assume $\hat{d} - d (t) = 0$ holds. In order to ensure that ${\overset{\cdot}{V}}_{s} (t) < 0$ holds, we need to discuss the positive and negative of $χ (t)$ .

Substituting (10) into (4) and (9) yields $\overset{\cdot}{X} (t) = A X (t) + \bar{Φ} s (t)$ . Define new Lyapunov function $V_{n} (t) = \sqrt{X^{T} (t) P X (t)}$ , and derivative of $V_{n}$ as follows

{\overset{\cdot}{V}}_{n} (t) = \frac{X^{T} (t) (A^{T} P + P A) X (t) + 2 X^{T} (t) P \bar{Φ} s (t)}{2 \sqrt{X^{T} (t) P X (t)}}

(22)

Substituting $P = P_{1} P_{1}^{T}, V_{n} (t) = \sqrt{X^{T} (t) P X (t)} = ∥ P_{1}^{T} X (t) ∥, ∥ P_{1}^{T} X (t) ∥ = ∥ X^{T} (t) P_{1} ∥,$ and (17) into (22), then

\begin{matrix} {\overset{\cdot}{V}}_{n} (t) \leq \frac{- 2 ξ X^{T} (t) P X (t) + 2 ∥ P_{1}^{T} X (t) ∥ ∥ P_{1}^{T} \bar{Φ} ∥ ∥ s (t) ∥}{2 \sqrt{X^{T} (t) P X (t)}} \\ \leq - ξ V_{n} (t) + ε_{1} ∥ s (t) ∥ \end{matrix}

(23)

Define $e (t) = V_{n} (t) - ϑ_{1} (t)$ and $e (0)$ as the initial value of $e (t)$ . Combining ${\overset{\cdot}{ϑ}}_{1} (t) = - ξ ϑ_{1} (t) + ε_{1} ∥ s (t) ∥$ with (23), we have

\overset{\cdot}{e} (t) \leq - ξ e (t)

(24)

Integrating the formula (24) yields $e (t) \leq e^{- ξ t} e (0) .$ Combining $e (t) = V_{n} (t) - ϑ_{1} (t),$ one obtains $V_{n} (t) \leq ϑ_{1} (t) + e^{- ξ t} e (0) .$ Substitute $P_{1}^{T}$ and $X (t)$ into the equation $V_{n} (t) = ∥ P_{1}^{T} X (t) ∥$ , then

V_{n} (t) = \sqrt{∥ P_{11}^{T} {\bar{x}}_{1} (t) + P_{21}^{T} ν (t) ∥^{2} + ∥ P_{22}^{T} ν (t) ∥^{2}}

(25)

Hence $(k_{1} + k_{2}) ∥ P_{11}^{T} {\hat{\bar{x}}}_{1} (t) + P_{21}^{T} ν (t) ∥ \leq (k_{1} + k_{2}) V_{n} (t) .$ Combining $k_{1} + k_{2} > 1$ with $χ (t) = ψ_{0} ∥ P_{11}^{T} {\bar{x}}_{1} (t) + P_{21}^{T} ν (t) ∥ - (k_{1} + k_{2}) ψ_{1} (t)$ one obtains

\begin{matrix} χ (t) \leq ψ_{0} (k_{1} + k_{2}) V_{n} (t) - (k_{1} + k_{2}) ψ_{1} (t) \\ \leq ψ_{0} (k_{1} + k_{2}) (ψ_{0} V_{n} (t) - ψ_{1} (t)) \end{matrix}

(26)

Since $V_{n} (t) \leq ϑ_{1} (t) + e^{- ξ t} e (0)$ and $ψ_{1} (t) = ψ_{0} (ϑ_{1} (t) + ϱ)$ are known, then

χ (t) \leq ψ_{0} (k_{1} + k_{2}) (e^{- ξ t} e (0) - ϱ)

(27)

In the following, both the cases $e (0) > 0$ and $e (0) \leq 0$ will be discussed.

If $e (0) > 0$ , since $e (t) \leq e^{- ξ t} e (0)$ at $t = 0$ , thus $e (t) \leq e (0)$ and $χ (t) \leq ψ_{0} (e (0) - ϱ)$ . Assume $e (0) > ϱ$ , there exists $t > 0$ such that $e (t) \geq ϱ$ in $t \in [0, T_{1}]$ and $e (t) < ϱ$ in $t > T_{1}$ . When $t \in [0, T_{1}]$ , then

\begin{matrix} {\overset{\cdot}{V}}_{s} (t) \leq - τ ∥ s (t) ∥^{2} + ψ_{0} (e (0) - ϱ) ∥ s (t) ∥ \\ \leq - (τ ∥ s (t) ∥ - ψ_{0} (e (0) - ϱ)) ∥ s (t) ∥ \\ \leq - \sqrt{2} (τ ∥ s (t) ∥ - ψ_{0} (e (0) - ϱ)) V_{s}^{\frac{1}{2}} (t) \end{matrix}

(28)

which indicates that, if $∥ s (t) ∥ \geq \frac{ψ_{0} (e (0) - ϱ)}{τ},$ we have ${\overset{\cdot}{V}}_{s} (t) \leq 0,$ and thus there always exists $T_{2} > 0$ such that $t > T_{2}$ , $∥ s (t) ∥ < \frac{ψ_{0} (e (0) - ϱ)}{τ} .$

Define $T_{m} = \max {T_{1}, T_{2}}$ . If $t > T_{m}$ , then $e (t) < ϱ$ . According to the definition $e (t) = V_{n} (t) - ϑ_{1} (t)$ , thus $V_{n} (t) = ϑ_{1} (t) + e (t) < ϑ_{1} (t) + ϱ$ and substitute the formula into $χ (t) \leq ψ_{0} (k_{1} + k_{2}) V_{n} (t) - (k_{1} + k_{2}) ψ_{1} (t)$ to obtain $χ (t) \leq ψ_{0} (k_{1} + k_{2}) (ϑ_{1} (t) + ϱ) - (k_{1} + k_{2}) ψ_{1} (t) = 0 .$

Since $χ (t) \leq 0$ , hence (21) can be written as

{\overset{\cdot}{V}}_{s} (t) < - τ V_{s} < 0

(29)

which shows that $s (t)$ will converge to origin.

If $e (0) \leq ϱ$ , obviously $χ (t) \leq ψ_{0} (e (0) - ϱ) < 0$ , thus ${\overset{\cdot}{V}}_{s} (t) < 0$ . The inequality (29) can be determined immediately, and the desired result can be obtained.

When $e (0) \leq 0$ , since $χ (t) < 0$ , hence ${\overset{\cdot}{V}}_{s} (t) < 0$ , then (29) is also satisfied.

This implies that if $\hat{d} - d (t) = 0$ is guaranteed, then $s (t)$ will converge to zero and the system can reach the sliding surface. ■

Lemma 1. (Zuo and Lin, 2016)³ Consider the following scalar system

\overset{\cdot}{y} = - c_{1} sgn (y) | y |^{α 1} - c_{2} sgn (y) | y |^{α 2}, y (0) = y_{0}

(30)

where $c_{1}, c_{2}$ are positive constants, $0 < α_{1} < 1$ and $α_{2} > 1$ . The equilibrium point of the above system is finite time stable and its settling time is bounded.

T < T_{\max} \overset{Δ}{=} \frac{1}{c_{1}} \frac{1}{1 - α_{1}} + \frac{1}{c_{2}} \frac{1}{α_{2} - 1}

(31)

Design sliding mode disturbance observer

In this section, the sliding mode disturbance observer is designed to guarantee $s (t)$ converge to zero.

The following auxiliary variable is introduced

e_{o} (t) = z (t) - {\bar{x}}_{2} (t)

(32)

where

\begin{matrix} \overset{\cdot}{z} (t) = A_{21} {\bar{x}}_{1} (t) + A_{22} {\bar{x}}_{2} (t) + B_{2} u (t) + υ_{z} \\ - c_{1 d} sgn (e_{o} (t)) | e_{o} (t) |^{α_{1 d}} - c_{2 d} sgn (e_{o} (t)) | e_{o} (t) |^{α_{2 d}} \end{matrix}

(33)

The coefficients $c_{1 d}, c_{2 d}, α_{1 d}, α_{2 d}$ are positive constants with $α_{1 d} < 1, α_{2 d} > 1$ . $υ_{z}$ will be introduced later.

Proposition 3. Consider the following sliding surface $s_{d}$

\begin{matrix} s_{d} (t) = {\overset{\cdot}{e}}_{o} (t) + c_{1 d} sgn (e_{o} (t)) | e_{o} (t) |^{α_{1 d}} \\ + c_{2 d} sgn (e_{o} (t)) | e_{o} (t) |^{α_{2 d}} \end{matrix}

(34)

If $υ_{z}$ satisfies the following equation

{\overset{\cdot}{υ}}_{z} (t) = - (k_{d} (t) + η_{d}) sgn (s_{d} (t)), υ_{z} (0) = 0

(35)

where $η_{d}$ is a positive constant. $k_{d} (t) > d_{1}$ , $d_{1} > | B_{2} \overset{\cdot}{d} |$ and $d_{0} > | B_{2} d |$ . $d$ and $\overset{\cdot}{d}$ are disturbance and its derivatives, $d_{1}$ and $d_{0}$ are the upper bounds of $| B_{2} \overset{\cdot}{d} |$ and $| B_{2} d |$ , respectively. $e_{o}$ converges to origin in finite time and the estimation of the disturbance as follows

υ_{z} (t) = B_{2} \hat{d}

(36)

Proof. Consider the derivative of $e_{o}$ in equation (32) as follows

\begin{matrix} {\overset{\cdot}{e}}_{o} (t) = υ_{z} (t) - c_{1 d} sgn (e_{o} (t)) | e_{o} (t) |^{α_{1 d}} \\ - c_{2 d} sgn (e_{o} (t)) | e_{o} (t) |^{α_{2 d}} - B_{2} d (t) \end{matrix}

(37)

Substitute (37) into (34) one obtains

s_{d} (t) = υ_{z} (t) - B_{2} d (t)

(38)

Differentiating the above equation (38) and combining with equation (35) yields

\begin{matrix} {\overset{\cdot}{s}}_{d} (t) = {\overset{\cdot}{υ}}_{z} (t) - B_{2} \overset{\cdot}{d} (t) \\ = - (k_{d} (t) + η_{d}) sgn (s_{d} (t)) - B_{2} \overset{\cdot}{d} (t) \end{matrix}

(39)

Consider the Lyapunov function $V_{d} (t) = \frac{1}{2} s_{d}^{2} (t)$ , then $V_{d} (t)$ derivative is as follows

\begin{matrix} {\overset{\cdot}{V}}_{d} (t) = s_{d} (t) {\overset{\cdot}{s}}_{d} (t) \\ = - k_{d} (t) | s_{d} (t) | - η_{d} | s_{d} (t) | - B_{2} \overset{\cdot}{d} (t) s_{d} (t) \\ \leq (| B_{2} \overset{\cdot}{d} (t) | - k_{d} (t)) | s_{d} (t) | - η_{d} | s_{d} (t) | \end{matrix}

(40)

Since $η_{d}$ is a positive constant and $k_{d} (t) > d_{1} > | B_{2} \overset{\cdot}{d} (t) |$ , then (40) can be written as follows

{\overset{\cdot}{V}}_{d} (t) = s_{d} (t) {\overset{\cdot}{s}}_{d} (t) \leq - η_{d} | s_{d} (t) |

(41)

The $s_{d} (t)$ will converges to zero in finite time, by utilizing equation (34) one gets

\begin{matrix} {\overset{\cdot}{e}}_{o} (t) = - c_{1 d} sgn (e_{o} (t)) | e_{o} (t) |^{α_{1 d}} \\ - c_{2 d} sgn (e_{o} (t)) | e_{o} (t) |^{α_{2 d}} \end{matrix}

(42)

By considering Lemma 1, the equilibrium of (42) as $e_{o} (t) = 0$ is finite time stable.

Consider the estimation error of disturbance as follows

e_{d} (t) = B_{2} d (t) - B_{2} \hat{d}

(43)

Both (36) and (38) are satisfied, $e_{d} (t)$ will convergence to origin in finite time.

This implies that if $k_{d} (t) > d_{1}$ is satisfied, the proposed disturbance observer can estimate the disturbance $d (t)$ in finite time.

Design adaptation structure

In this subsection, two adaptation structures are designed, and the assumptions on the disturbance are relaxed in the proposed observer. As is capable of seeing from the proof in the above subsection, the adaptive gain $k_{d} (t)$ must satisfy $k_{d} (t) > d_{1}$ to ensure the system sliding takes place, in which case reachability condition is achieved³⁰ and sliding motion takes place on $s_{d} = {\overset{\cdot}{s}}_{d} = 0$ . In other words, unknown disturbance or uncertainty should be completely eliminated, that is $| u_{eq} (t) | = | d_{1} | .$

In the first adaptive structure, similar to general adaptive control,²⁴ it is assumed that $d_{1}$ is known and $d_{0}$ is unknown. The second adaptive strategy is designed which assumes that both $d_{1}$ and $d_{0}$ are unknown, the assumptions on the disturbance are more relaxed.

Consider utilizing a low pass filter to filter the switching signal to obtain a close approximation. Here, the $sgn (s_{d})$ should take on the average value. Then if ${\hat{u}}_{eq} (t)$ satisfies

{\overset{\cdot}{\hat{u}}}_{eq} (t) = \frac{1}{l} (- (k_{d} (t) + η) sgn (s_{d}) - {\hat{u}}_{eq} (t))

(44)

where $l > 0$ is a time constant, then ${\hat{u}}_{eq} (t)$ almost completely approximates $u_{eq} (t)$ . In order to eliminate the influence of the initial conditions of the filter, assume exist $0 < υ_{1} < 1$ and $υ_{0} > 0$ such that

| | {\hat{u}}_{eq} (t) | - | u_{eq} (t) | | < υ_{1} | u_{eq} (t) | + υ_{0}

(45)

holds. The adaptive algorithm of the control gain is driven by using the equivalent control. Through the following inequality, we introduce the concept of ’safety margin’

k_{d} (t) > \frac{1}{μ} | {\hat{u}}_{eq} (t) | + υ

(46)

where $0 < μ < 1$ and $υ > 0$ are design parameters such that

\frac{1}{μ} | {\hat{u}}_{eq} (t) | + \frac{υ}{2} > | u_{eq} (t) |

(47)

holds. Define error variable $e_{k} (t)$ as follows

e_{k} (t) = k_{d} (t) - \frac{1}{μ} | {\hat{u}}_{eq} (t) | - υ

(48)

which shows that if $e_{k} (t) = 0$ , then $k_{d} (t) > d_{1}$ , that is, the system will maintain sliding motion. Define adaptive scheme

{\overset{\cdot}{k}}_{d} (t) = - g (t) sgn (s_{d})

(49)

where $g (t)$ is auxiliary scalar, it has the following form

g (t) = f_{0} + f (t)

(50)

where $f_{0}$ is a positive scalar and $f (t)$ is elaborated later. In this paper, according to whether $d_{1}$ is known, $f (t)$ will execute different choices. In the next work, we will discuss two situations: $d_{1}$ is known and $d_{1}$ is unknown.

Case 1:

Assuming $d_{1}$ is known but $d_{0}$ is unknown. Define

e_{r} (t) = \frac{α d_{1}}{μ} - f (t)

(51)

\overset{\cdot}{f} (t) = β | e_{k} (t) | + f_{0} \sqrt{β} sgn (e_{r} (t))

(52)

where $β > 0$ , and $α > 1$ is designed to ensure $| {\overset{\cdot}{\hat{u}}}_{eq} (t) | < α d_{1}$ .

Proposition 4. As described above, the problem of maintaining sliding is converted to how to yield $e_{k} (t) \to 0$ in finite time; how to allow $g (t)$ and $k_{d} (t)$ to be bounded. In the following, we will verify the situation when $d_{1}$ is known.

Proof. Consider following Lyapunov function

V (t) = \frac{1}{2} e_{k}^{2} (t) + \frac{1}{2 β} e_{r}^{2} (t)

(53)

It follows from (53) that

\begin{matrix} e_{k} (t) {\overset{\cdot}{e}}_{k} (t) \leq - f_{0} | e_{k} (t) | + (e_{r} (t) - \frac{α d_{1}}{μ}) | e_{k} (t) | \\ + \frac{α d_{1}}{μ} | e_{k} (t) | \\ \leq - f_{0} | e_{k} (t) | + e_{r} (t) | e_{k} (t) | \end{matrix}

(54)

\frac{e_{r} (t) {\overset{\cdot}{e}}_{r} (t)}{β} = - e_{r} (t) | e_{k} (t) | - \frac{f_{0} | e_{k} (t) |}{\sqrt{β}}

(55)

The derivative with respect to (53) and integrate with (54) and (55), one obtains

\begin{matrix} \overset{\cdot}{V} (t) = e_{k} {\overset{\cdot}{e}}_{k} + \frac{e_{r} {\overset{\cdot}{e}}_{r}}{β} \\ \leq - f_{0} | e_{k} (t) | - \frac{f_{0} | e_{r} |}{\sqrt{β}} \\ \leq - \sqrt{2} f_{0} (\frac{| e_{k} (t) |}{\sqrt{2}} + \frac{| e_{r} (t) |}{\sqrt{2 β}}) \\ \leq - f_{0} \sqrt{2 V (t)} \end{matrix}

(56)

Define $V_{0}$ to be the initial value of $V (t)$ and $t_{0}$ to be the time taken for $V (t)$ converge to zero, and hence integrating both sides yields $t_{0} < \frac{\sqrt{2 V_{0}}}{f_{0}}$ , which shows that $e_{k} (t)$ and $e_{r} (t)$ will converge to origin in finite time. Then $k_{d} (t) > d_{1}$ will be guaranteed and the reachability condition is satisfied. Obviously $e_{r} (t)$ is bounded, hence $f (t)$ and $g (t)$ is bounded.

Case 2:

Assuming that both $d_{1}$ and $d_{0}$ are unknown. Define

\overset{\cdot}{f} (t) = {\begin{matrix} β | e_{k} (t) | if | e_{k} (t) | > e_{0} \\ 0 otherwise \end{matrix}

(57)

where $e_{0}$ is a design constant which to eliminate the noise signal in the system.

Proposition 5. Both $d_{1}$ and $d_{0}$ are unknown, if $υ$ is chosen to satisfy

\frac{4 e_{0}^{2}}{υ^{2}} + \frac{4}{υ^{2} β} {(\frac{α d_{1}}{μ})}^{2} < 1

(58)

then $| e_{k} (t) | < \frac{υ}{2}$ is realized in finite time and the sliding motion is guaranteed.

Proof. Consider the Lyapunov function from (53)

e_{k} {\overset{\cdot}{e}}_{k} \leq - f_{0} | e_{k} (t) | + e_{r} (t) | e_{k} (t) |

(59)

Suppose $f (0) = 0$ , from (57) one obtains $\overset{\cdot}{f} (t) \geq 0$ , then $f (t) \geq 0$ holds. Since $e_{r} (t) = \frac{α d_{1}}{μ} - f (t)$ , then $e_{r} (t) \leq \frac{α d_{1}}{μ}$ is satisfied. In (57), if $| e_{k} (t) | > e_{0}$ with (55), ${\overset{\cdot}{e}}_{r} (t) = - \overset{\cdot}{f} (t) = - β | e_{k} (t) |$ is established, (56) can be rewritten as

\begin{matrix} \overset{\cdot}{V} (t) = e_{k} (t) {\overset{\cdot}{e}}_{k} (t) + \frac{e_{r} (t) {\overset{\cdot}{e}}_{r} (t)}{β} \\ \leq - f_{0} | e_{k} (t) | + e_{r} (t) | e_{k} (t) | - \frac{e_{r} (t) \overset{\cdot}{f} (t)}{β} \\ \leq - f_{0} | e_{k} (t) | + e_{r} (t) | e_{k} (t) | - e_{r} (t) | e_{k} (t) | \\ \leq - f_{0} | e_{k} (t) | \end{matrix}

(60)

If $| e_{k} (t) | < e_{0}$ with $e_{r} (t) < 0,$ then $\overset{\cdot}{V} (t) \leq - f_{0} | e_{k} (t) | .$ Since $e_{r} (t) \leq \frac{α d_{1}}{μ},$ then outside of the interval $E = {(e_{k}, e_{r}) : | e_{k} | < e_{0}, 0 \leq e_{r} < \frac{α d_{1}}{μ}}$ and in the solution domain subject to $\overset{\cdot}{V} \leq - f_{0} | e_{k} (t) |$ . Draw the smallest ellipse $F$ with the origin as the center to enclose the $E$ , then $F = {(e_{k}, e_{r}) : V (e_{k}, e_{r}) < r},$ where $r > 0$ represents the ’radius’ of the ellipse, and $r = \frac{e_{0}^{2}}{2} + \frac{1}{2 β} {(\frac{α d_{1}}{μ})}^{2} .$ Choose the appropriate $υ$ to satisfy (58), the sliding motion will is guaranteed.

Since $E \subset F$ and outside of $F$ in the solution domain, $\overset{\cdot}{V} (t) \leq 0$ holds, hence $F$ is an invariant set. If the $V (e_{k}, e_{r})$ enters $F$ , then $V (e_{k}, e_{r})$ will not be able to leave $F$ and from (58), $| e_{k} (t) | < \frac{υ}{2}$ is satisfied. If the $V (e_{k}, e_{r})$ does not enter $F$ , then from the above discussion $\overset{\cdot}{V} \leq - f_{0} | e_{k} (t) |$ and $\int_{0}^{\infty} f_{0} | e_{k} (t) | dt \leq V (0)$ are satisfied. Since $V (e_{k}, e_{r})$ is bounded, then $e_{k} (t), e_{r} (t)$ and ${\overset{\cdot}{e}}_{k} (t), {\overset{\cdot}{e}}_{r} (t)$ are bounded, hence $| e_{k} (t) |$ is uniformly continuous, that is, $e_{k} (t) \to 0$ when $t \to \infty$ . It shows that there exists $t_{0}$ such that $| e_{k} (t) | \leq \frac{υ}{2}$ for $t > t_{0}$ . Hence, $| e_{k} (t) | < \frac{υ}{2}$ has always been established in finite time.

Since the $| e_{k} (t) | < \frac{υ}{2}$ holds, therefore $| e_{k} (t) | = | k_{d} (t) - \frac{1}{μ} | {\hat{u}}_{eq} (t) | - υ | < \frac{υ}{2}$ from (48), that is $k_{d} (t) - \frac{1}{μ} | {\hat{u}}_{eq} (t) | - υ > - \frac{υ}{2}$ , and from (47), one obtains

k_{d} (t) > \frac{1}{μ} | {\hat{u}}_{eq} (t) | + \frac{υ}{2} > | u_{eq} (t) | > d_{1}

(61)

Hence reachability condition is satisfied and the system will maintain sliding motion. Since $e_{k} (t)$ and $e_{r} (t)$ are bounded, consequently $f (t)$ and $g (t)$ remains bounded. From (45)

\begin{matrix} | k_{d} (t) | < | e_{k} (t) | + \frac{1}{μ} | {\hat{u}}_{eq} (t) | + υ \\ < | e_{k} (t) | + (1 + υ_{1}) \frac{d_{0}}{μ} + υ \end{matrix}

(62)

obviously $k_{d} (t)$ is bounded, and the proof is complete.

Remark 2. The method proposed in this paper can be improved by some new fuzzy systems, at the same time, the ASMDO does not require knowledge of the disturbance and its first derivative, hence it can better estimate the disturbances in these systems. ^31,32

Remark 3. The control strategy of this paper can be used to deal with the problem of energy/voltage management in photovoltaic (PV)/battery systems ³³ and interval type-3 fuzzy logic systems. ³⁴ The sliding mode disturbance observer can be used to eliminate the influence of the variation of temperature, radiation, and output load on the system, and effectively guarantee better disturbance rejection performance of the controller.

Simulation example

A practical example is provided in this section to illustrate the efficiency of the obtained result. Consider an inverted pendulum on a cart model given in Edwards and Spurgeon⁸ and the linearization model is as follows

A = [\begin{matrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & - 1.9333 & - 1.9872 & 0.0091 \\ 0 & 36.9771 & 6.2589 & - 0.1738 \end{matrix}]

B = [\begin{matrix} 0 \\ 0 \\ 0.3205 \\ - 1.0095 \end{matrix}] C = [\begin{matrix} 3 & 0 & 0 & 0 \\ 0 & 2 & 0 & 0 \\ 0 & 0 & 2 & 0 \end{matrix}]

The non-singular matrix $T$ can be found as

T = [\begin{matrix} 0 & 0 & 3.1498 & 1.0000 \\ - 3.0000 & 0 & 0 & 0 \\ 0 & 2.0000 & 0 & 0 \\ 0 & 0 & - 2.0000 & 0 \end{matrix}]

then we can obtain the following transformed regular form

\bar{A} = [\begin{matrix} - 0.1454 & 0 & 15.4438 & - 0.2284 \\ 0 & 0 & 0 & 1.5000 \\ 2.0000 & 0 & 0 & 3.1498 \\ - 0.0182 & 0 & 1.9333 & - 2.0159 \end{matrix}]

\bar{B} = [\begin{matrix} 0 \\ 0 \\ 0 \\ - 0.6410 \end{matrix}] \bar{C} = [\begin{matrix} 0 & - 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & - 1 \end{matrix}]

Choose $Φ$ and $G_{2}$ as follows

Φ = [\begin{matrix} - 0.005 \\ - 0.005 \\ 0.2 \end{matrix}] G_{2} = [\begin{matrix} 5 & 0 & 0 \\ 0 & 5 & 0 \\ 0 & 0 & 10 \end{matrix}]

Use the LMI toolbox to obtain these parameter matrices $w_{1}, w_{2}, w_{3}, w_{4}$

w_{1} = [\begin{matrix} 0 & 0.0779 \end{matrix}] w_{2} = [\begin{matrix} 0.0294 & - 1.9900 \\ 0.0913 & - 0.9198 \\ - 0.1462 & 2.3990 \end{matrix}]

w_{3} = [\begin{matrix} - 0.1988 & 1.5439 & 1.3524 \\ 0.3063 & - 0.1739 & 0.8509 \\ 19.5966 & - 2.0623 & - 3.6779 \end{matrix}]

w_{4} = [\begin{matrix} - 16.2153 & 1.6530 & 0.4266 \end{matrix}]

The controller (18) parameters are selected as $τ = 1.8$ , $k_{1} = 2$ , $k_{2} = 2$ , $h_{1} = 3$ , $h_{2} = 0.2$ , $ξ = 0.028$ , $ϱ = 0.48$ . The initial condition is $\bar{x} (0) = [0.1 0.01 0.2 0.15]$ . This system is simulated with disturbance: $d (t) = 2 \sin (2 t) + 2 \cos (2 t)$ . The parameters of the sliding mode disturbance observer are selected as $c_{1 d} = c_{2 d} = 2, α_{1 d} = 0.5, α_{2 d} = 3, η_{d} = 0.01$ . The disturbance is set to $B_{2} \overset{\cdot}{d} (t) = | 2 \sin (π t) + 0.15 \cos (2 π t) |$ . The parameters of adaptive schemes are chosen as $d_{1} = 1$ , $μ = 0.99$ , $f_{0} = 0.6$ , $e_{0} = 0.01$ , $υ = 0.18$ , $l = 0.01$ , and starting with zero initial conditions.

Figures 1 and 2 are the control input of the system obtained from Zhang et al.¹¹ and (18), respectively. Obviously, the proposed controller possesses less chattering. The fluctuation of the control input from 0 to 6 s is caused by the system states ${\bar{x}}_{1} (t)$ and ${\bar{x}}_{2} (t)$ in the controller (18). It can be seen from Figure 4 that when the system state tends to be stable, the input of the controller will also tend to stable. The fluctuation of control input after 6 s is caused by $k_{1} sg n^{h 1} (s (t)) + k_{2} sg n^{h 2} (s (t))$ in the controller (18). Figure 3 shows that after the sliding surface converges to $0$ , there is a slight fluctuation around zero, which causes the $sgn (s)$ in the controller to switch between −1 and 1, thus generating the fluctuations in Figure 2. Figure 4 depicts that system state by Zhang et al.¹¹ and (18), respectively. Obviously, the system states obtained with the proposed controller have less overshoot and faster convergence performance. The coefficient of additional term $k_{1} + k_{2} > 1$ located in denominator, although the item is omitted in $t_{s}$ , actually, it must enable the system to converge faster.

Figure 1.

Control input by applying controller.¹¹

Figure 2.

Control input by applying proposed controller.

Figure 3.

Trajectory of sliding surface.

Figure 4.

Comparisons of system states: (a and b) are “position” signals, (c and d) are “velocity” signals.

Figure 5 is evolutions of disturbance and its estimation in controller. Compared with Hwang et al.,²⁹ the estimation of disturbance can faster track its true value in finite time under the proposed observer (33). Figures 6 and 7 shown that the adaptive gain $k_{d} (t)$ still closely follows $B_{2} \overset{\cdot}{d} (t)$ . As in (48), $k_{d} (t)$ converges to a safety margin which depends on the parameters $μ$ and $υ$ . $k_{d} (t)$ is always above $B_{2} \overset{\cdot}{d} (t)$ , the sliding motion will be maintained and the conditions on disturbance bounds are relaxed, the proposed ASMDO does not require information about the bound on the disturbances and their derivatives.

Figure 5.

Disturbance $d (t)$ and its estimation $\hat{d} (t)$ in controller.

Figure 6.

Adaptive gain $k_{d} (t)$ and disturbance $B_{2} \overset{\cdot}{d} (t)$ in the proposed observer ( $d_{1}$ is known but $d_{0}$ is unknown).

Figure 7.

Adaptive gain $k_{d} (t)$ and disturbance $B_{2} \overset{\cdot}{d} (t)$ in the proposed observer ( $d_{1}$ and $d_{0}$ are unknown).

Conclusion

This paper has addressed a disturbance observer-based control method for continuous-time linear systems with unknown external disturbance. Based on a novel adaptive sliding mode disturbance observer, the output feedback sliding mode controller has been designed, which guarantees that reachability condition holds strictly. In the absence of upper bound information of the disturbance and its first derivative, the restrictive restraints on disturbance have been relaxed by designing the adaptive laws. It is worth mentioning that the estimation of disturbance can track its true value in finite time under the proposed observer. An inverted pendulum on a cart model has been exploited to illustrate the effectiveness of the proposed controller and observer. Simulation results show that the convergence performance of the controller can be further improved, and the finite time control observer-based method will be studied in future.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported in part by the Natural Science Foundation of Tianjin under Grant Number: 19JCTPJC48600.

ORCID iD

Chao Jiang

References

. Digital integral terminal sliding mode predictive control of piezoelectric-driven motion system. IEEE Trans Ind Electron 2016; 63(6): 3976–3984.

Sun

. Discrete-time fractional order terminal sliding mode tracking control for linear motor. IEEE Trans Ind Electron 2018; 65(4): 3386–3394.

Zuo

Lin

. Distributed robust finite-time nonlinear consensus protocols for multi-agent systems. Int J of Syst Sci 2016; 47(6): 1366–1375.

Balootaki

Rahmani

Moeinkhah

, et al. On the synchronization and stabilization of fractional-order chaotic systems: Recent advances and future perspectives. Phys A Stat Mech Appl 2020; 551: 124203.

Man

. Model reference adaptive control systems with terminal sliding modes. Int J Control 1996; 64(6): 1165–1176.

Feng

Man

. Terminal sliding mode control – an overview. IEEE Open J Ind Electron Soc 2021; 2: 36–52.

Lee

Chang

Jin

. Adaptive integral sliding mode control with time-delay estimation for robot manipulators. IEEE Trans Ind Electron 2017; 64(8): 6796–6804.

Edwards

Spurgeon

. Sliding mode control: theory and applications. New York, NY: CRC Press, 1998.

Zhang

Xia

. Design of static output feedback sliding mode control for uncertain linear systems. IEEE Trans Ind Electron 2010; 57(6): 2161–2170.

10.

Chang

. Static output feedback sliding mode controller design for minimum phase uncertain systems. IMA J Math Control Inf 2015; 32(3): 485–495.

11.

Zhang

Shi

Xia

. Design of sliding mode output feedback controllers via dynamic sliding surface. Automatica 2021; 124: 109310.

12.

Seuret

Edwards

Spurgeon

, et al. Static output feedback sliding mode control design via an artificial stabilizing delay. IEEE Trans Automat Contr 2009; 54(2): 256–265.

13.

Song

Wang

Niu

. Static output-feedback sliding mode control under round-robin protocol. Int J Robust Nonlinear Control 2018; 28(18): 5841–5857.

14.

Lai

Edwards

Spurgeon

. On output tracking using dynamic output feedback discrete-time sliding-mode controllers. IEEE Trans Automat Contr 2007; 52(10): 1975–1981.

15.

Yang

, et al. A novel asynchronous control for artificial delayed Markovian jump systems via output feedback sliding mode approach. IEEE Trans Syst Man Cybern Syst 2019; 49(2): 364–374.

16.

Chen

Yang

Guo

, et al. Disturbance-observer-based control and related methods—an overview. IEEE Trans Ind Electron 2016; 63(2): 1083–1095.

17.

Wang

Luo

Wang

, et al. Integral sliding mode control using a disturbance observer for vehicle platoons. IEEE Trans Ind Electron 2020; 67(8): 6639–6648.

18.

Yang

. Continuous finite-time anti-disturbance control for a class of uncertain nonlinear systems. Trans Inst Meas Contr 2014; 36(3): 300–311.

19.

Alinaghi Hosseinabadi

Ordys

Soltani Sharif Abadi

, et al. State and disturbance observers-based chattering-free fixed-time sliding mode control for a class of high-order nonlinear systems. Adv Control Appl 2021; 3(3): e81.

20.

Lee

. Nonlinear disturbance observer-based robust control of attitude tracking of rigid spacecraft. Nonlinear Dyn 2017; 88(2): 1317–1328.

21.

Gao

. Fixed-time formation control of AUVs based on a disturbance observer. Acta Autom Sin 2019; 45(6): 1094–1102.

22.

Plestan

Shtessel

Brégeault

, et al. New methodologies for adaptive sliding mode control. Int J Control 2010; 83(9): 1907–1919.

23.

Negrete-Chávez

Moreno

. Second-order sliding mode output feedback controller with adaptation. Int J Adapt Control Signal Process 2016; 30(8–10): 1523–1543.

24.

Utkin

Poznyak

. Adaptive sliding mode control with application to super-twist algorithm: Equivalent control method. Automatica 2013; 49(1): 39–47.

25.

Roy

Baldi

Fridman

. On adaptive sliding mode control without a priori bounded uncertainty. Automatica 2020; 111: 108650.

26.

Abadi

ASS

Hosseinabadi

Mekhilef

. Fuzzy adaptive fixed-time sliding mode control with state observer for a class of high-order mismatched uncertain systems. Int J Control Autom Syst 2020; 18(10): 2492–2508.

27.

Tan

Man

. Terminal sliding mode observers for a class of nonlinear systems. Automatica 2010; 46(8): 1401–1404.

28.

Mobayen

Javadi

. Disturbance observer and finite-time tracker design of disturbed third-order nonholonomic systems using terminal sliding mode. J Vib Control 2017; 23(2): 181–189.

29.

Hwang

Park

Joo

. Disturbance observer-based integral fuzzy sliding-mode control and its application to wind turbine system. IET Control Theory Appl 2019; 13(12): 1891–1900.

30.

Utkin

. Sliding modes in control and optimization. Berlin: Springer-Verlag, 1992.

31.

Mohammadzadeh

Kumbasar

. A new fractional-order general type-2 fuzzy predictive control system and its application for glucose level regulation. Appl Soft Comput 2020; 91: 106241.

32.

Mohammadzadeh

Rathinasamy

. Energy management in photovoltaic battery hybrid systems: a novel type-2 fuzzy control. Int J Hydrogen Energy 2020; 45(41): 20970–20982.

33.

Mosavi

Qasem

Shokri

, et al. Fractional-order fuzzy control approach for photovoltaic/battery systems under unknown dynamics, variable irradiation and temperature. Electronics 2020; 9(9): 1455.

34.

Liu

Mohammadzadeh

Turabieh

, et al. A new online learned interval type-3 fuzzy control system for solar energy management systems. IEEE Access 2021; 9: 10498–10508.