Tensor Product Model Transformation Based Adaptive Integral-Sliding Mode Controller: Equivalent Control Method

Abstract

This paper proposes new methodologies for the design of adaptive integral-sliding mode control. A tensor product model transformation based adaptive integral-sliding mode control law with respect to uncertainties and perturbations is studied, while upper bounds on the perturbations and uncertainties are assumed to be unknown. The advantage of proposed controllers consists in having a dynamical adaptive control gain to establish a sliding mode right at the beginning of the process. Gain dynamics ensure a reasonable adaptive gain with respect to the uncertainties. Finally, efficacy of the proposed controller is verified by simulations on an uncertain nonlinear system model.

1. Introduction

Sliding mode control (SMC) has advantages like ease of implementation and reduction in the order of the state equation. It also has the ability to withstand external disturbances and model uncertainties satisfying the matching condition, that is, the perturbations that enter the state equation at the same point as the control input. SMC design comprises two steps: a sliding surface is first constructed such that the system trajectories along the sliding surface meet the specified performance and then a discontinuous or continuous control law drives the states towards the sliding surface, and keeps them on the surface thereafter, regardless of disturbances or parasitic uncertainties. The resulting controller, although robust against matched perturbations, still suffers from reaching phase problem, that is, an initial period of time in which the system has not yet reached the sliding surface and it is sensitive to perturbations satisfying the matching condition.

In order to solve the reaching phase problem, an integral-sliding surface was proposed [1, 2]. Successively, many application examples of integral-sliding mode control have been developed in the literature and different systems are studied, such as multi-input multi-output (MIMO) linear system [3], uncertain nonlinear system [4, 5], uncertain stochastic system [6], and uncertain fractional order nonlinear system [7]. Some practical applications are also investigated [8–10]. The basic idea of integral-sliding surface control design is to define the nominal control law, which is responsible for the performance of the nominal system, that is, without perturbations, and a discontinuous or continuous control law, which is responsible for suppression of the parasitic uncertainties and external disturbances or perturbations. In integral-sliding surface control, an integral term is included in the sliding surface, and it allows to define the sliding surface in such a way that the system trajectories start in the sliding motion right at the beginning of the process. Moreover, it is well known that terminal sliding mode control (TSMC) guarantees a finite convergence time, by combining derivative and integral terminal sliding modes in a recursive structure; derivative-integral terminal sliding mode control (DI-TMSC) methods are proposed to achieve exact or approximate finite-time convergence for the output tracking of higher order nonlinear systems [11]. More recently, several works on the TP model transformation based sliding mode control appeared in the literature [12–16].

TP model transformation was first introduced by Baranyi et al. [17, 18]. It allows transforming dynamic system models defined over bounded domains into convex parameter-varying weighted combination of linear time invariant (LTI) systems. This transformation can transform a function into TP function form if it is feasible. If an exact transformation is impossible to achieve, then the method will determine a TP function which can be taken as an approximation. This feature offers the TP model transformation the ability to provide a tradeoff between approximation accuracy and complexity of the resulting TP model. linear parameter varying (LPV) based control has developed based on TP model transformation and linear matrix inequality (LMI) [17–22]. TP model transformation has found many applications in new direction in control theory [16, 23–25]. For further information., please refer to recently published monograph for qLPV control theories [26] and papers [19, 27].

Inspired by the adaptive sliding mode controller [28], the TP model transformation based parallel distributed compensation (TPDC) technique is applied, and a TP model transformation based adaptive integral-sliding mode controller (TPAISMC) is proposed in this paper. By combining with 𝜎-adaptation strategy [29], a real sliding surface is realized; due to uncertainties and external disturbances, the phenomenon that adaptive gain is increasing all the time when ideal sliding mode does not occur is avoided. Chattering caused by perturbations and uncertainties [30] is suppressed by TPAISMC with 𝜎-adaptation. Furthermore, the design process can be numerically realized without analytical derivation due to the characteristics of the TP model transformation based parallel distributed compensation. Thus, this control scheme has potential applications in many practical environments.

The rest of the paper is organized as follows. In Section 2, the concept of TP model transformation is introduced, and the control problem is formulated. In Section 3, TP model transformation based adaptive integral-sliding model control is presented, and TP model transformation based adaptive integral-sliding model control with 𝜎-adaptation method is proposed. In Section 4, the simulations on TP model transformation based adaptive integral-sliding model controllers with different adaptive gain dynamics are carried out, and the control performance is tested. Finally in Section 5 the conclusions are given.

2. Preliminaries

2.1. Tensor Product Model Transformation

It is known that a given bounded continuous function 𝑔(𝑥) (𝑥∈ℝ^𝑁) can be approximated by tensor product over a compact domain Ω. First, divide the compact domain into 𝑁-dimensional hyper rectangular grid based on a given grid size; on each dimension, the locations are defined by the vector 𝑎_𝑛=𝑎_𝑛,1,𝑎_𝑛,2,…,𝑎_{𝑛,𝑀𝑛},𝑛=1,2,…,𝑁, (1) where 𝑎_𝑛,1≤𝑎_𝑛,2≤⋯≤𝑎_{𝑛,𝑀𝑛} and the 𝑁-dimensional hyper rectangular grid belongs to a bounded region Ω=[𝑎_1,1,𝑎_1,𝑀1]×[𝑎_2,1,𝑎_2,𝑀2]×⋯×[𝑎_𝑁,1,𝑎_{𝑁,𝑀𝑁}]. Then function 𝑔(𝑥) can be sampled over the grid points as 𝑔(𝑥)=𝑏_{𝑖1,𝑖2,…,𝑖𝑁}, (2) where 𝑥=[𝑎_1,𝑖1,𝑎_2,𝑖2,…,𝑎_{𝑁,𝑖𝑁}] (for all 𝑛: 𝑖_𝑛=1,2,…,𝑀_𝑛). Sampled matrices are stored into a tensor 𝒮^𝑠 where superscript “𝑠” means “sampled.” Then higher order singular value decomposition (HOSVD) is used to root out the redundant content of the sampled set. It follows that 𝑔(𝑥)≈𝑔(𝑥)=ℬ⊗^𝑁_𝑛=1𝑤_𝑛𝑥_𝑛=^𝐼1_𝑖1^𝐼₌₁²_𝑖2=1⋯^𝐼𝑁_𝑖𝑁^𝑁_=1𝑛=1𝑤_{𝑛,𝑖𝑛}𝑥_𝑛𝑏_{𝑖1,𝑖2,…,𝑖𝑁}, (3) where tensor ℬ∈ℝ^{𝐼1×𝐼2×⋯×𝐼𝑁} and ℝ^{𝐼1×𝐼2×⋯×𝐼𝑁} is the vector space of real-valued 𝐼₁×𝐼₂×⋯×𝐼_𝑁 tensors; 𝐼_𝑛 is the number of univariate basis functions used in the 𝑛th dimension of the parameter vector and 𝑏_{𝑖1,𝑖2,…,𝑖𝑁} is the entry of ℬ. Tensor ℬ satisfies the following properties.(i)All-orthogonality: ⟨ℬ_{𝑖𝑛=𝛼},ℬ_{𝑖𝑛=𝛽}⟩=0 for all possible values of 𝛼, 𝛽, and 𝑛 (𝛽≠𝛼).(ii)Ordering: ‖ℬ_𝑖𝑛=1‖≥‖ℬ_𝑖𝑛=2‖≥⋯≥‖ℬ_{𝑖𝑛=𝐼𝑛}‖≥0 for all possible values of 𝑛, where ‖⋅‖ is the Frobenius norm. Row vector 𝑤_𝑛(𝑥_𝑛)∈ℝ^𝐼𝑛 contains basis functions, and 𝑤_𝑛(𝑥_𝑛)=[𝑤_𝑛,1(𝑥_𝑛),𝑤_𝑛,2(𝑥_𝑛),…,𝑤_{𝑛,𝐼𝑛}(𝑥_𝑛)] is normalized to one, where 𝑤_𝑛,𝑖(𝑥_𝑛)∈[0,1] (𝑛=1,…,𝑁; 𝑖=1,…,𝐼_𝑁). For the LPV state space model, TP model representation can be obtained as 𝑥̇=^𝑟_𝑖=1𝑤_𝑖𝑆_𝑖𝑝𝑥(𝑡)𝑢(𝑡)=^𝑟_𝑖=1𝑤_𝑖𝐴_𝑖𝑥+𝐵_𝑖𝑢, (4)where 𝑆_𝑖 are the convex hull of the LTI vertex systems and 𝑟=∏^𝑁_𝑗=1𝐼_𝑗 is the linear index equivalent of 𝑁-dimensional array’s index 𝑖₁,𝑖₂,…,𝑖_𝑁. Note that if all the discarded singular values are zero, this is the exact form of the original tensor and the approximation error is zero. If nonzero singular values are discarded, TP approximation error results in the following representation: 𝑔(𝑥)−^𝑟_𝑖=1𝑤_𝑖𝐴_𝑖𝑥+𝐵_𝑖𝑢≤𝛿, (5) where 𝛿 is used to symbolize the approximation error; it is a function of discarded singular values and 𝑤_𝑖∈[0,1] is the coefficient function.

2.2. Problem Statement

Consider a class of 𝑛th-order uncertain nonlinear systems. The control system is represented in a state-space form as 𝑥̇=𝑓(𝑥)+𝑏(𝑥)𝑢+𝑤, (6) where 𝑥=[𝑥₁,𝑥₂,…,𝑥_𝑛]^𝑇∈ℝ^𝑛 is the state vector, 𝑢∈ℝ is the control input, 𝑓(𝑥)∈ℝ^𝑛 and 𝑏(𝑥)∈ℝ^𝑛 are the smooth vector fields, and 𝑤∈ℝ^𝑛 is the external disturbance or parasitic uncertainties. Note that (6) is in a controllable form. The system and input matrices 𝑓(𝑥) and 𝑏(𝑥) are not exactly known but assumed to be in the form 𝑓(𝑥)=𝑓₀(𝑥)+Δ𝑓(𝑥) and 𝑏(𝑥)=𝑏₀(𝑥)+Δ𝑏(𝑥), where 𝑓₀(𝑥) and 𝑏₀(𝑥) are the nominal parts of 𝑓(𝑥) and 𝑏(𝑥), respectively. Assume that the nominal part 𝑓₀(𝑥)+𝑏₀(𝑥)𝑢 of (6) can be transformed into TP models; that is, 𝑆𝑝𝑥(𝑡)𝑢(𝑡)≈^𝑟_𝑖=1𝑤_𝑖𝐴_𝑖𝑥+𝐵_𝑖𝑢≜𝐴𝑝(𝑡)𝑥+𝐵𝑝(𝑡)𝑢. (7) Let 𝐸(𝑥,𝑢) be the lumped uncertainty and denoted as 𝐸(𝑥,𝑢)=Δ𝑓(𝑥)+Δ𝑏(𝑥)𝑢+𝛿+𝑤. (8) Then the uncertain nonlinear system (6) can be represented as 𝑥̇=𝐴𝑝(𝑡)𝑥+𝐵𝑝(𝑡)𝑢+𝐸(𝑥,𝑢). (9) Assume that the following hypothesis holds for the uncertain nonlinear system (6).

Assumption 1. The nominal parts 𝑓₀(𝑥) and 𝑏₀(𝑥) of uncertain nonlinear system (6) can be transformed into TP model.

Assumption 2. The matrix product 𝐶𝐵(𝑝) is invertible.

Assumption 3. The actual values of Δ𝑓(𝑥), Δ𝑏(𝑥), and 𝑤 are unknown, but their bound should exist for all 𝑥 and 𝑡.

Assumption 4. The known nominal nonlinear plant of (6), 𝑥̇(𝑡)=𝑓₀(𝑥,𝑡)+𝑏₀(𝑥)𝑢(𝑡), is globally asymptotically stabilizable via a nominal control 𝑢_eq; that is, there is a nonempty set of Lyapunov function 𝑉 such that for any choice of ℂ¹ function 𝑉(𝑥,𝑡)∈𝑉∶ℝ^𝑛×ℝ⁺→ℝ⁺ satisfies 𝛾₁(‖𝑥‖)≤𝑉(𝑥,𝑡)≤𝛾₂(‖𝑥‖),𝜕𝑉+𝜕𝑡𝜕𝑉𝜕𝑥^𝑇𝑓(𝑥)+𝑏(𝑥)𝑢_eq≤−𝛾(‖𝑥‖), (10) where 𝛾₁,𝛾₂∶ℝ⁺→ℝ⁺ are class 𝒦_∞ functions and 𝛾∶ℝ⁺→ℝ⁺ is a class 𝒦 function, where 𝒦 is the set of all continuous functions 𝛼∶[0,∞)→[0,∞) that are zero at zero, strictly increasing, and continuous and 𝒦_∞ is the subset of 𝒦 consisting of those functions that are unbounded.
The following nonlinear integral-type sliding surface [31, 32] is used in the adaptive sliding mode control design: 𝜎(𝑥,𝑡)=𝐶𝑥(𝑡)−𝐶𝑥𝑡₀−𝐶^𝑡_𝑡0𝐴𝑝(𝜏)𝑥+𝐵𝑝(𝜏)𝑢_eq𝑑𝜏=0. (11)

Remark 1. The reaching phase is eliminated for nonlinear integral-type sliding surface defined by (11). It is obvious that the sliding mode exists from the very beginning. The additional integral ∫^𝑡_𝑡0[𝐴(𝑝)𝑥+𝐵(𝑝)𝑢_eq]𝑑𝜏+𝑥(𝑡₀) provides another degree of freedom than the linear sliding surface, and this term can be thought of as a trajectory of the nominal system, that is, the system in the absence of perturbations and in the presence of the nominal control 𝑢_eq. In this work, equivalent control 𝑢_eq is given by TP model transformation based parallel distributed compensation (PDC) technique [17]. Moreover, 𝑓₀(𝑥) is assumed to be in the linear form; that is, 𝑓₀(𝑥)=𝐴𝑥 in [31, 32] is removed. The requirement of 𝑓₀(𝑥) is just to satisfy Assumption 1 in this paper. Thus, our method can be applied to a variety of nonlinear systems in which TP model transformation can be applied, for example, inverted pendulum system, underactuated robot manipulators [33], translational oscillator with a rotational actuator (TORA) [34], vertical takeoff and landing (VTOL) aircraft [35], and quadrotor helicopter [36].

3. TP Model Transformation Based Adaptive Integral-Sliding Mode Control Design

The adaptive sliding mode control law using the integral-sliding surface (11) can be designed as 𝑢(𝑡)=𝑢_eq+𝑢_as,𝑢_eq=−^𝑟_𝑖=1𝑤_𝑖𝑝(𝑡)𝐾_𝑖𝑢𝑥,_as=−𝑘(𝑡)𝐶𝐵𝑝⁻¹sgn(𝜎(𝑥,𝑡)), (12) where the gain 𝑘(𝑡) is the adaptive value preserving the sliding mode and sgn(⋅) is the signum function. Updating law is designed as 𝑘̇(𝑡)=𝑘|𝜎(𝑥,𝑡)|, (13) where 𝑘>0 and 𝑘(0)>0. For the nominal feedback, TP model transformation based parallel distributed compensation (TPDC) is used; gains 𝐾_𝑖 are computed by the following LMI-based stability.

Theorem 2 (see [21]). The polytopic model (4) with the equivalent control 𝑢_𝑒𝑞 (12) is asymptotically stable if there exist 𝑋>0 and 𝑀_𝑖 that satisfy −𝑋𝐴^𝑇_𝑖−𝐴_𝑖𝑋+𝑀^𝑇_𝑖𝐵^𝑇_𝑖+𝐵_𝑖𝑀_𝑖>0, (14) for 𝑖=1,2,…,𝑟, and −𝑋𝐴^𝑇_𝑖−𝐴_𝑖𝑋−𝑋𝐴^𝑇_𝑠−𝐴_𝑠𝑋+𝑀^𝑇_𝑠𝐵^𝑇_𝑖+𝐵_𝑖𝑀_𝑠+𝑀^𝑇_𝑖𝐵^𝑇_𝑠+𝐵_𝑠𝑀_𝑖≥0, (15) for 𝑖<𝑠≤𝑟, except the pairs (𝑖,𝑠) such that 𝑤_𝑖(𝑝(𝑡))𝑤_𝑠(𝑝(𝑡))=0, for all 𝑝(𝑡), and the feedback gains are determined by the solutions 𝑋 and 𝑀_𝑖 as 𝐾_𝑖=𝑀_𝑖𝑋⁻¹. (16)

Theorem 3 (see [21]). Assume that ‖𝑥(0)‖≤𝜙, where 𝑥(0) is unknown, but the upper bound 𝜙 is known. The constraint ‖𝑢(𝑡)‖₂≤𝛾 is enforced at all times 𝑡≥0 if the following LMIs hold: 𝜙²𝐼≤𝑋,𝑋𝑀^𝑇_𝑖𝑀_𝑖𝛾²𝐼≥0. (17)

The aforementioned control algorithm can be summarized as follows.

Theorem 4. Consider the uncertain nonlinear system (6). Given the integral-type sliding surface (11) and the adaptive control law (12), the sliding mode is guaranteed to be reached in finite time.

Proof. Assume that there exists a final gain 𝑘^∗ such that 𝑢_as=−𝑘^∗(𝐶𝐵(𝑝))⁻¹sign(𝜎(𝑥,𝑡)) is the optimal solution for 𝑢_as, and 𝑘^∗>|𝐶||𝐸(𝑥,𝑢)|. It is always true since adaptive gain 𝑘(𝑡) is increasing. Differentiating the sliding surface (11) with respect to time, it yields 𝜎̇=𝐶𝑥̇(𝑡)−𝐶𝐴𝑝(𝑡)𝑥−𝐶𝐵𝑝(𝑡)𝑢_eq+𝐶𝐸(𝑥,𝑢). (18) A Lyapunov function is chosen as 𝑉=𝜎²/2+𝑘²/(2𝑘), where 𝑘=𝑘−𝑘^∗. Differentiating 𝑉 with respect to time using (12) and (18), it yields ̇̇𝑘̇𝑘𝑉=𝜎𝜎+𝑘=𝜎𝐶𝑥̇(𝑡)−𝐶𝐴𝑝𝑥−𝐶𝐵𝑝𝑢_eq+𝑘−𝑘^∗̇𝑘𝑘=𝜎𝐶𝐸(𝑥,𝑢)+𝐶𝐵𝑝𝑢_eq+𝑢_as−𝐶𝐵𝑝𝑢_eq+𝑘−𝑘^∗|𝜎(𝑥,𝑡)|. (19) By applying 𝑢_as=−𝑘(𝑡)(𝐶𝐵(𝑝))⁻¹sgn(𝜎(𝑥,𝑡)), (19) becomes 𝑉̇=𝜎|𝐶𝐸(𝑥,𝑢)|−𝑘^∗|𝜎|<0. (20) It can be inferred that 𝜎 and 𝑘 reach zero in finite time; that is, 𝜎→0 and 𝑘−𝑘^∗=𝑘→0⇒𝑘→𝑘^∗. The convergence of 𝜎 and 𝑘 is proven by Lyapunov stability criterion. In practice, the discontinuous switching function sgn(𝜎) in (12) may cause chattering in the sliding mode. In order to alleviate the input chattering, switching function sgn(𝜎) is replaced with the saturation function sat(𝜎/𝜓), where Φ is the boundary layer width which should be carefully chosen according to each particular control system. This completes the proof.

Remark 5. Briefly, the proposed TPAISMC scheme provides the following three main advantages. First, the knowledge of the upper bound of the system uncertainties is not required. Second, TP model transformation has the ability to provide a tradeoff between approximation accuracy and complexity of the resulting TP function. If an exact transformation is impossible, approximation error is bounded by discarded singular values, then this error can be incorporated into lumped uncertainty 𝐸(𝑥,𝑢). Third, the tracking performance of the control system is guaranteed.
The above adaptive method does not require knowledge of uncertainties or perturbations bound. Adaptive gain 𝑘(𝑡) may increase all the time when 𝜎(𝑥,𝑡)≠0. However, this happens all the time in real applications due to sampled computation, noisy measurements, or other nonidealities, sliding surface can hardly achieve zero exactly. The method referred to as “𝜎-adaptation” proposed in [29] is adopted to counteract the uncertainties and perturbations effects; based on integral-sliding surface and TP model transformation, the adaptive sliding mode controller with 𝜎-adaptation is given as follows: 𝑢(𝑡)=𝑢_eq+𝑢_as,𝑢_eq=−^𝑟_𝑖=1𝑤_𝑖𝑝(𝑡)𝐾_𝑖𝑢𝑥,_as=−𝑘(𝑡)𝐶𝐵𝑝⁻¹sgn(𝜎(𝑥,𝑡)), (21) where the gain 𝑘(𝑡) is defined as 𝑘̇(𝑡)𝑘|𝜎(𝜎()(𝑡)|𝜎((𝑡)𝑥,𝑡)|sgn(|𝑥,𝑡)|−𝜖if𝑘>𝜇,𝑥,𝑡)|if𝑘≤𝜇, (22) where 𝜇 is a prescribed positive parameter. By Lemma 2 [29], there exists a positive constant 𝑘^∗ such that 𝑘(𝑡)≤𝑘^∗ for all 𝑡>0 and 𝑘^∗>|𝐶𝐸(𝑥,𝑢)|. Furthermore, (18) can be simplified as 𝜎̇=𝐶𝐸(𝑥,𝑢)+𝐶𝐵𝑝𝑢_as (23) by uncertain nonlinear system (6) and 𝑢(𝑡) is given by (21). Choose Lyapunov function as 𝑉=𝜎²/2+𝑘²/(2𝑘). Differentiating 𝑉 with respect to time, using (21) and (23), it yields ̇̇𝑘̇𝑘𝑉=𝜎𝜎+𝑘|𝐶𝐸(|𝜎|+≤𝜎||𝑥,𝑢)|−𝑘𝑘−𝑘^∗|𝜎|sgn(|𝜎|−𝜖), (24) when 𝑘>𝜇. Suppose that |𝜎|>𝜖, we have 𝑉̇≤𝜎|𝐶||𝐸(𝑥,𝑢)|−𝑘^∗|𝜎|<0. (25) Suppose that |𝜎|<𝜖, 𝑉̇ could be sign indefinite, |𝜎| can increase over 𝜖, then adaptive gain will increase to guarantee the sliding motion, and then the control strategy ensures that |𝜎|<𝜖 again. When 𝑘≤𝜇, differentiating 𝑉 with respect to time, using (22), it yields 𝑉̇≤𝜎|𝐶||𝐸(𝑥,𝑢)|−𝑘|𝜎|+𝑘−𝑘^∗|𝜎(𝑥,𝑡)|<0. (26) In the adaptive sliding mode controller design, 𝜎-adaptation is used to regulate the adaptive gain, combining with TPPDC, the following theorem holds.

Theorem 6. Consider the uncertain nonlinear system (6). Given the integral-type sliding surface (11) and the control law (21), the sliding mode is guaranteed to be reached in finite time.

4. Simulations

Consider the following nonlinear system: 𝑥̇₁=𝑥₂,𝑥̇₂=𝑥₃,𝑥̇₃=−(1+0.3sin𝑡)𝑥²₁−(1.5+0.2cos𝑡)𝑥₂−(1+0.4sin𝑡)𝑥₃+3+cos𝑥₁𝑢+𝑤. (27) The nominal parts of system are 𝑓₀(𝑥)=[𝑥₂,𝑥₃,𝑥²₁−1.5𝑥₂−𝑥₃]^𝑇 and 𝑏₀(𝑥)=[0,0,3]^𝑇. TP model transformation is used to decompose [𝑓₀(𝑥),𝑏₀(𝑥)], and the LTI systems are given by 𝐴₁𝐵₁=010000103−1.5−13,𝐴₂𝐵₂=01000010−3−1.5−13. (28) Since the equivalent control method used in ASMC [28] is replaced with TPDC, the newly proposed method is referred to as TP based adaptive sliding mode control (TPASMC). States responses with TPASMC are shown in Figure 1(a), the sliding surface function is shown in Figure 1(b), and control input and adaptive gain parameter are given in Figures 2(a) and 2(b), respectively. We can conclude that the equivalent control based on TPDC used in TPASMC has the same performance as ASMC which is proposed by Huang et al. [28]. However, since the TPDC is a numerical method, we can design the adaptive sliding mode controller without knowing the analytical system model. There is no modeling error for the nominal system when all the nonzero singular values are used in TP model transformation. The boundary layer technique is applied to avoid high-frequency control activity; the signum function in 𝑢_as is replaced with sat(𝜎/𝜓), and we have 𝜎sat𝜓=𝜎𝜓𝜎𝜓𝜎≤1,sgn𝜓𝜎𝜓>1, (29) where 𝜓 is the boundary layer width. In the simulation, 𝛾=0.5, 𝜙=0.1, 𝐶=[12,7,1],Γ₀=5, 𝜓=0.05, the initial state vector for all the simulations is 𝑥(0)=[0,0,−0.35]^𝑇, and 𝛼=0.1 is the decay rate for TPDC. The desired trajectory is supposed to be 𝑥_𝑑(𝑡)=[sin𝑡,cos𝑡,−sin𝑡]^𝑇. The disturbance 𝑤 is supposed to be a random noise possessing a mean value of 0.5 and |𝑤|≤0.1.

(a)

(b)

Figure 1 

(a) States responses with TPASMC. (b) Sliding surface function with TPASMC.

(a)

(b)

Figure 2 

(a) Control input with TPASMC. (b) Adaptive gain parameter.

Figure 3(a) shows the states responses of TPASMC with integral-sliding surface (TPAISMC), sliding surface function is shown in Figure 3(b), and Figures 4(a) and 4(b) show the control input and adaptive gain parameter, respectively. To implement TPAISMC, the following parameters are used: 𝐶=[5,8,3],𝜓=0.005,Γ₀=2,𝛾=5,𝜙=0.5,𝛼=1. (30) For integral-type sliding surface, reaching phase is eliminated; this can be verified by Figure 3(b). However, in real application, the ideal sliding surface does not exist, then the adaptive gain may increase forever (as shown in Figure 4(b)) when the sliding surface reaches nonzero; the adaptive strategy (13) may cause chattering when the adaptive gain grows large.

(a)

(b)

Figure 3 

(a) States responses with TPAISMC. (b) Sliding surface function with TPAISMC.

(a)

(b)

Figure 4 

(a) Control input with TPASMC. (b) Adaptive gain parameter.

To avoid large adaptive gain and chattering, we apply the adaptation law (22) to reduce the chattering. For TPAISMC with 𝜎 adaptation, the following parameters are used: 𝐶=[5,8,3],𝑘𝑡₀=0.5,𝑘=0.05,𝜖=0.002,𝜇=0.6,𝜓=0.005,𝜙=0.5,𝛾=5,𝛼=1. (31) Figure 5(a) shows the states responses with TPAISMC when 𝜎-adaptation is used. Sliding surface function is shown in Figure 5(b) and control input and adaptive gain parameter are shown in Figures 6(a) and 6(b), respectively. Here it is clearly seen from Figure 5(b) that TPAISMC with 𝜎-adaptation is effective in chattering suppression. Adaptive gain 𝑘(𝑡), defining the chattering amplitude in the sliding mode, increases to a reasonable level to guarantee a real sliding mode.

(a)

(b)

Figure 5 

(a) States responses with TPAISMC (𝜎-adaptation). (b) Control input with TPAISMC (𝜎-adaptation).

(a)

(b)

Figure 6 

(a) Sliding surface function with TPAISMC (𝜎-adaptation). (b) Adaptive gain parameter (𝜎-adaptation).

5. Conclusion

In this work, an adaptive siding mode controller scheme for a class of uncertain nonlinear system is studied, and an integral-sliding surface was adopted in the designing of TPAISMC. Instead of using the traditional feedback control, TP model transformation based parallel distributed compensation controller was applied to stabilize the states. It is noted that the modeling error caused by TP model transformation can be zero if nonzero singular value is discarded. Nonzero modeling error can be modeled as lumped uncertainties, combining modeling error with perturbations and external disturbances, and can be suppressed effectively by adaptive gain controller with 𝜎-adaptation strategy.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 11101066; no. 61074044; no. 61374118; no. 61104038; no. 60834004) and the Fundamental Research Funds for the Central Universities (no. DUT13LK32), and partially supported by the Chinese National Basic Research 973 Program (2009CB320602).