# Adaptive Control for Nonlinear Systems with Time-Varying Control Gain

Authors: Alejandro Rincon; Fabiola Angulo
Journal: Journal of Control Science and Engineering
 (2012)

---

## Abstract

We propose a scheme for nonlinear plants with time-varying control gain and time-varying plant coefficients, on the basis of a plant model consisting of a Brunovsky-type model with polynomials as approximators. We develop an adaptive robust control scheme for this plant, under the following assumptions: (i) the plant terms involve time-varying but bounded coefficients, being its upper bound unknown; (ii) the control gain is unknown, not necessarily bounded, and only its signum is known. To achieve robustness, we use a combination of robustifying control inputs and dead zone-type update laws. We apply this methodology to the speed control of a permanent magnet synchronous motor (PMSM), and we achieve proper tracking results.


---

## Body

## 1. Introduction

Nonlinear behavior is difficult to model accurately, rendering controller design cumbersome. An approach to handle this is to use linear models, which have local validity in the state space, that is, different values of the parameter set would be required for each region in the state space. A second approach is to use a plant model in either the Brunovsky form defined in [1], or the parametric-pure feedback form defined in [2], using the so-called “function approximation techniques’’, such as locally weighted statistical learning [3], fuzzy sets [4], Fourier series [5], orthonormal functions, [6] or neural networks [7]. In the last case, state and time dependent terms, known as “basis functions,’’ are used to represent the nonlinear and time-varying behaviour. Some identification or learning methods are used to perform the approximation, resulting in a negligible modelling error if there are a sufficient number of basis functions [3]. In some cases, a Brunovsky form model of the plant is used, assuming that there is full-state measurement and that the basis functions are multiplied by unknown constants. Then, a sliding surface model reference adaptive control (SSMRAC) is devised, handling the residual approximation error by means of auxiliary robustifying inputs [7].Common drawbacks of the above-mentioned schemes are(i)
the convergence of the tracking error to some small residual set is not ensured in a strict sense and depends on the value of the approximation error [7];(ii)
upper or lower bounds of the plant coefficients are required to be known [8];(iii)
discontinuous auxiliary inputs are used, which may lead to chattering [5];(iv)
high enough gains are used, which require excessive values [6].Due to environment changes, the coefficients of the plant model may experience time-varying but bounded behavior [9, 10]. Both, the control gain and other coefficients may exhibit this behavior. The time-varying behavior of coefficients different from the control gain may be handled by means of robustness techniques. See for instance the robust MRAC scheme in [10], whose drawback is that projection- or σ-type update laws are used, and hence upper bounds of the plant coefficients are required to be known, in order to achieve the convergence of the tracking error to some residual set of user-defined size. The time-varying control gain may be handled by means of robustness techniques [11] or Nussbaum gain technique [12], a projection-type update law is used, such that a lower bound of the control gain is required to be known. In [11], a σ type update law is used, such that bounds of the plant coefficients are required to be known. In summary, the main drawback of using robustness techniques to handle varying control gain is that upper or lower bounds of some plant coefficients are required to be known. On the other hand, the Nussbaum gain technique relaxes this requirement. The main drawback is that the upper bound of the transient behavior of the state S is significantly altered with respect to that of the disturbance-free case: the value of this bound depends on time integrals of the terms that comprise the Nussbaum function. Another drawback is that the control gain is usually required to be upper-bounded by some constant.The scheme of [13] achieves adequate robustness properties. Mainly, upper bounds of the plant coefficients are not required to be known. The essential element of the technique is to define the quadratic forms related to the sliding surface 𝒮 in terms of a dead zone function of 𝒮 rather than in terms of 𝒮. This scheme has the following benefits.(Ri)
The tracking error converges to a residual set whose size is of the user’s choice;(Rii)
known upper bounds of the plant coefficients are not required, such that high enough gains are not used;(Riii)
all the closed-loop signals are bounded (parameter drifting is avoided);(Riv)
auxiliary control signals are not discontinuous in terms of both the tracking error and the sliding surface, hence input chattering is avoided;(Rv)
upper bounds of time-varying but bounded coefficients are not required to be known.The disadvantages of [13] are (i) the design is only developed for systems with hysteresis nonlinearities in the actuator; (ii) the time-varying control gain is tackled by means of the Nussbaum gain method, which entails higher complexity of the Lyapunov analysis; (iii) the control gain is assumed to be time-varying but bounded.In contrast to the approaches of [3, 4, 14–16], we propose that an adequate regression model for highly nonlinear systems may be obtained from a Brunovsky type model, inserting polynomials to approximate the nonlinear behavior, and considering also: (i) time-varying but bounded behavior of some plant coefficients; (ii) unknown control gain, time-varying, and not necessarily bounded. We develop a robust adaptive scheme for this plant, achieving benefits (Ri) to (Rv). If we compare the robust technique developed here with the technique developed in [13], the main differences, which are also contributions, are(Rvi)
we consider unknown time-varying control gain, not necessarily bounded, not restricted to actuator nonlinearities;(Rvii)
we consider time-varying but bounded behavior of some plant terms;(Rviii)
we tackle the control gain by means of robustness techniques, which gives a simpler design in comparison with the Nussbaum technique.This paper is organized as follows. The outline of the scheme is given in Section2. The function approximation is discussed in Section 3. The plant regression model and the statement of the problem are established in Section 4. The control and update laws are derived in Section 5. The control algorithm and its stability are presented in Section 6. An application of the scheme to a PMSM is given in Section 7. Finally, the conclusions are presented in Section 8.


## 2. Outline of the Scheme

We propose the use of polynomials to approximate the nonlinear behavior, taking into account the fact that polynomials are universal approximators for continuous functions, according to [17]. We consider a Brunovsky-form plant model, into which we introduce the polynomial terms.We devise a robust adaptive controller for this plant, achieving benefits (Ri) to (Rix) mentioned in the introduction. We use the SSMRAC method stated in [18], as the basic framework for designing the control and update laws. To handle the effect of modelling error and time-varying behavior of plant terms, we use a robust continuous control input, whose magnitude is adjusted, and a dead zone-type update law. The resulting controller has the following features: (i) the magnitude of the control input is adjusted to cope with the unknown upper bounds of the time-varying coefficients; (ii) the control input is proportional to some continuous function of the sliding surface 𝒮, so that chattering is avoided; (iii) an inactivation is introduced in all the update laws, which occurs when the sliding surface 𝒮 reaches a target region. If we compare this with projection-type update laws, it has the advantage of not requiring the upper bounds of the plant parameters.For the stability analysis, we use a truncated version of the quadratic form related to the sliding surface, denoted byVs̅. The validity of this technique, including the conditions of the Lyapunov function, is stated in [19]. It is worth noticing that the standard conditions of the Lyapunov theory for time-variable systems, shown in [18] or [20], are not satisfied because of the truncation. We prove the asymptotic convergence of the tracking error in a rigorous manner by means of Barbalat’s Lemma. To that end, we redefine the expression of V̇ as an inequality in terms of Vs̅.


## 3. Function Approximation Based on Polynomials

According to [17, page 116], a continuous real-valued function f̅(x), where x∈D, D⊂Rn, being D a compact set, may be approximated by a polynomial function g̅ in the interval x∈D, with some finite error, being g̅ defined as(1)g̅(x,θo)=∑k1,k2,…,knθo[j]x1k1x2k2⋯xnkn:θo∈Ω,Ω⊂Rp,
where Ω is a convex set and θo[j] means the jth element of the vector θo. Thus, f̅(x) can be expressed in terms of g̅ as follows (cf. [17] page 122):(2)f̅(x)=g̅(x,θa)+ϵo(x),|ϵo(x)|<ϵ,θa={θo∈Ω:θo=argmin⁡θo⁡(sup⁡x∈D⁡|f̅(x)-g̅(x,θo)|)},
where ϵ is a positive constant, and it is known as approximation error, θa is an ideal parameter set and g̅(x,θa) is an ideal representation of f̅(x). The polynomial can be linearly parameterized with respect to its coefficients:(3)f̅(x)=ϕ⊤θa+ϵo(x),
where ϕ contains polynomial terms, whereas θa contains constants. In this work, f̅ represents the nonlinear part of the dynamical equation. We will handle the unknown constant vector θa by means of adjustment rules and ϵo by means of robust techniques.


## 4. Plant Model and Problem Statement

We assume that the dynamical nonlinear system can be represented by a Brunovsky type model, as defined in [1], with polynomial functions:(4)y(n)=-cn-1y(n-1)-⋯-c1ẏ-coy+f̅(x)+bu,(5)x=[y,ẏ,…,y(n-1)]⊤,
where y(t)∈ℝ1 is the system output u(t)∈ℝ1 is the input and the coefficients cn-1,…, and co are unknown and time-varying but bounded. The relative degree n may be established by means of previous step response analysis. Inserting the function (3) into (4) gives(6)y(n)=bu+ϕ⊤θa+d,(7)d=ϵo(x)-cn-1y(n-1)-⋯-c1ẏ-coy,|cn-1|≤μn-1,…,|c1|≤μ1,|co|≤μn,|ϵo(x)|≤μn+1=ϵ.
We make the following assumptions:(Ai)
the entries ofθa and the terms μ1,…, and μn+1 are constant and unknown;(Aii)
the control gainb varies with respect to the state x, or time, so that it satisfies the following:
(8)(i)0<bm≤|b|≤fb,∀t≥to,(9)(ii)sign⁡(b)constant,∀t≥to,
where the value of sign(b) is known, bm is an unknown positive constant, and fb is an unknown positive function of time or x. Notice that condition (8) implies that b is always different from zero;(Aiii)
the entries of the vectorx are available for measurement;(Aiv)
the entries of the vectorϕ are known functions of x;(Av)
the value of the desired trajectoryyd(t) and its derivatives yd(n-1),…,ẏd is bounded.The desired outputyd is specified in terms of a bounded external command r(t) as follows:(10)yd(n)+am,n-1yd(n-1)+⋯+am,oyd=am,or,
where am,n-1,…,am,o are constant coefficients prespecified by the user, such that the polynomial K(p) is Hurwitz with at least one real root, being K(p) defined as K(p)=p(n)+am,n-1p(n-1)+⋯+am,o. The external reference signal r(t) must be bounded. The objective of the MRAC design is to formulate a control law u(t) such that the tracking error e(t)=y(t)-yd(t) asymptotically converges to the residual set De, defined as follows:(11)De={e:|e|≤Cbe},
where Cbe is a user-defined bound.


## 5. The Control and Update Laws

Let𝒮 the dynamics imposes over the tracking error given by(12)S=(p+λ)n-1e=pn-1e+λn-2pn-2e+⋯+λ1ė+λoe,
where λ is a positive constant chosen by the user. Having defined 𝒮, we establish the dynamic equation for 𝒮 by differentiating with respect to time:(13)Ṡ=pne+λn-2pn-1e+⋯+λ1ë+λoė,(14)Ṡ=y(n)-yd(n)+λn-2pn-1e+⋯+λ1ë+λoė,(15)Ṡ=y(n)+φa,(16)φa=-yd(n)+λn-2pn-1e+⋯+λ1ë+λoė.
We insert the expression for y(n) of (6) into (15):(17)Ṡ=bu+θa⊤ϕ+d+φa.
Define(18)u*=-θa⊤ϕ-φa-amS=φ⊤θ*,θ*=[θa⊤,1]⊤,(19)φ=[-ϕ⊤,-φa-amS]⊤,
where φa is defined in (16), am is a positive constant. We express (17) in terms of u* as follows:(20)Ṡ=-amS+bu+d-u*=-amS+bu+d-φ⊤θ*.
Let(21)l=n+1,(22)c*=2+l2Cbvs(12bmn)2(23)Cbvs≜(12)(λn-1Cbe)2,(24)f1=|ẏ|,…,fl-2=|y(n-1)|,fl-1=|y|,fl=1,μ1=max⁡{|c1|},…,μl-2=max⁡{|cn-1|},μl-1=max⁡{|co|},μl=ϵ.
Now, we multiply (20) by 𝒮 and add and subtract the term c*S2:(25)SṠ=-amS2+bSu+Sd-Sφ⊤θ*+c*S2-c*S2,
where the terms c*𝒮2-c*𝒮2, -𝒮φ⊤θ*, and 𝒮d can be expressed in terms of adjustment errors:(26)(i)c*S2-c*S2=-c*S2-c̃S2+ĉS2,c̃=ĉ-c*;(27)(ii)-Sφ⊤θ*=Sφ⊤θ̃-Sφ⊤θ̂,θ̃=θ̂-θ*;(28)(iii)Sd≤(-1)c̃1|S|f1+ĉ1|S|f1+⋯+(-1)c̃l|S|fl+ĉl|S|fl,(29)c̃1=ĉ1-μ1,…,c̃l=ĉl-μl,
where ĉ,θ̂,ĉ1,…,ĉl are adjusted parameters whose update laws will be defined later. Inserting the above properties into (25) gives(30)SṠ≤-amS2-c̃S2+Sφ⊤θ̃+(-1)c̃1|S|f1+⋯+(-1)c̃l|S|fl+ĉS2-Sφ⊤θ̂+ĉ1|S|f1+⋯+ĉl|S|fl+bSu-c*S2.
If b were constant, we would choose the control u so as to cancel the terms involving updated parameters ĉ,θ̂,ĉ1,…,ĉl. Since b is varying, we chose u to attenuate the effect of adjusted parameters, being the remaining error rejected by -c*𝒮2:(31)u=(-1)sgn⁡(b)ĉ2S3+(-1)sgn⁡(b)(φ⊤θ̂)2S+(-1)sgn⁡(b)ĉ12f12S+⋯+(-1)sgn⁡(b)ĉl2fl2S,
where ĉ,θ̂,ĉ1,…,ĉl are adjusted parameters whose update laws will be defined later. Replacing the above control law into (20) and (30) gives(32)Ṡ=-amS-|b|ĉ2S3-|b|(φ⊤θ̂)2S-|b|ĉ12f12S+⋯-|b|ĉl2fl2S+d-φ⊤θ*,(33)SṠ≤-amS2-c̃S2+Sφ⊤θ̃+(-1)c̃1|S|f1+⋯+(-1)c̃l|S|fl+ĉS2-Sφ⊤θ̂+ĉ1|S|f1+⋯+ĉl|S|fl-|b|ĉ2S4-|b|(φ⊤θ̂)2S2-|b|ĉ12f12S2-⋯-|b|ĉl2fl2S2-c*S2,
where the last terms of the above equation satisfy the following inequality (see the proof in Appendix A):(34)ĉS2-Sφ⊤θ̂+ĉ1|S|f1+⋯+ĉl|S|fl-|b|ĉ2S4-|b|(φ⊤θ̂)2S2-|b|ĉ12f12S2+⋯-|b|ĉl2fl2S2-c*S2≤0ifS2≥2Cbvs.
Equation (34) expresses the attenuation of the effect of the adjusted parameters. Substituting (34) into (33) gives(35)SṠ≤-amS2-c̃S2+Sφ⊤θ̃-c̃1|S|f1-⋯-c̃l|S|flifS2≥2Cbvs,
or equivalently,(36)SṠ≤-amS2-c̃S2+Sφ⊤θ̃-|S|C̃d⊤fifS2≥2Cbvs,(37)C̃d=[c̃1,…,c̃l]⊤,f=[f1,…,fl]⊤.
We choose the following update laws:(38)ĉ̇={γcS2ifS2≥2Cbvs,0otherwise,Ĉ̇d={Γdf|S|ifS2≥2Cbvs,0otherwise,θ̂̇={-ΓφSifS2≥2Cbvs,0otherwise,
where 𝒮 is defined in (12), Cbvs in (23), whereas γc, and the diagonal entries of Γ, Γd are positive constants chosen by the user, being Γ and Γd diagonal matrices.


## 6. The Control Algorithm

Now we recall the equations corresponding to the controller and establish the tracking convergence theorem. The control law is given by (31); the update laws are given by (38); the terms φ and 𝒮 are defined in (19) and (12); φa is defined in (16); Cbvs is defined in (23) and f is defined in ((37)). The adjusted parameters ĉ1,…,ĉl, required to compute u, are the entries of the vector Ĉd, that is, Ĉd=[ĉ1,…,ĉl]⊤.


### 6.1. Theorem: Boundedness and Tracking Convergence

If the controller designed in Section5 is applied to the plant (6), then the tracking error e(t) converges to De asymptotically, and the closed-loop signals remain bounded.Proof.
Now we proceed to analyze the stability of the controlled system using the direct Lyapunov method and the Barbalat’s Lemma. First, we establish the boundedness of the closed-loop signals on the basis of the time derivative of the Lyapunov function. Then, we establish the convergence of the tracking error to the target regionDe defined in (11), on the basis of the Barbalat’s Lemma. The validity of the procedure can be stated using the theory developed and applied in [18, 19, 21, 22].
The closed-loop dynamics is given by (32), and (38). We define the following truncated quadratic form, on the basis of the truncation presented in [19]:
(39)V̅s={VsifVs≥Cbvs,CbvsifVs<Cbvs,(40)Vs≜(12)S2,
where Cbvs is defined in (23). The form Vs̅ has the following properties:
(41)-Vs<-Vs̅+Cbvs<0ifVs>Cbvs,-Vs<-Vs̅+Cbvs=0ifVs=Cbvs,-Vs≤0=-Vs̅+Cbvs=0ifVs<Cbvs.
The states of the closed-loop system may be grouped in the following vector:
(42)x̅=[S,c̃,C̃d⊤,θ̃⊤]⊤,
where c̃, C̃d, and θ̃ are defined in (26), (37), (27). We define the following function:
(43)V(x̅(t))=Vs̅+Vθ+Vc+Vd,Vθ=(12)θ̃⊤Γ-1θ̃,Vc=(12)γc-1c̃2,Vd=(12)C̃d⊤Γd-1C̃d.
Notice that V(x̅) does not satisfy standard conditions of Lyapunov theory for time-variable systems, shown in [18]. It satisfies the conditions of [19] instead, which we use to prove the boundedness of the closed-loop signals and the convergence of the tracking error to the residual set De (11). The time derivative of Vs along the trajectory (32) is
(44)V̇s=SṠ=-amS2-|b|ĉ2S4-|b|(φ⊤θ̂)2S2+(-1)|b|ĉ12f12S2+⋯+(-1)|b|ĉl2fl2S2+S(d-φ⊤θ*)≤-amS2-c̃S2+Sφ⊤θ̃+(-1)|S|C̃d⊤fifVs≥Cbvs,
where the last inequality is obtained from (36). The time derivative of Vs̅ can be derived from (39) and expressed on the basis of the above equation:
(45)Vs̅̇={V̇s≤-amS2-c̃S2+Sφ⊤θ̃+(-1)|S|C̃d⊤fifVs≥Cbvs,0ifVs<Cbvs.
The time derivative V̇θ+V̇c+V̇d along trajectories (38) is
(46)V̇θ+V̇c+V̇d={(-1)Sθ̃⊤φ+c̃S2+C̃df|S|ifVs≥Cbvs,0otherwise.
Therefore, V̇ is given by
(47)V̇=Vs̅̇+V̇θ+V̇c+V̇d={V̇s+V̇θ+V̇c+V̇d≤-amS2ifVs≥Cbvs,0ifVs<Cbvs.
The above equation can be rewritten in terms of Vs defined in (40):
(48)V̇≤-2amVsifVs≥Cbvs,V̇=0ifVs<Cbvs.
This implies that V̇≤0 for all t≥to. Thus, according to [19], all the closed-loop signals are bounded, that is, (𝒮, θ̃, c̃, C̃d)∈L∞, or equivalently, x̅∈L∞, where x̅ is defined in (42). As a consequence, (V,  V̅s, Vs, Vθ, Vc, Vd)∈L∞, according to the definitions in (39) and (43). The boundedness of x̅ implies (e, …, e(n-1))∈L∞ according to (12); hence, 𝒮̇∈L∞ [22], and V̅̇s∈L∞. To establish the convergence of the tracking error, we begin by expressing the equation ahead in terms of Vs̅, according to the definition of Vs̅ in (39) and the properties of (41):
(49)V̇≤-2amVs<-2amVs̅+2amCbvs≤0ifVs≥Cbvs,V̇=0=2am(-Vs̅+Cbvs)=0ifVs<Cbvs,⟹V̇≤-2amVs̅+2amCbvs≤0.
We reorganize the above equation and integrate as follows:
(50)2am∫tot(Vs̅-Cbvs)dτ+V(t)≤V(to).
Thus, (Vs̅-Cbvs)∈L1. Recall that (Vs̅,Vs̅̇)∈L∞. Thus, (Vs̅-Cbvs)∈L∞∩L1, (d/dt)(Vs̅-Cbvs)∈L∞. By invoking the Barbalat Lemma [23], we obtain (Vs̅-Cbvs)→0 as t→∞. In turn, this implies that Vs→Dvs, being Dvs defined as
(51)Dvs={Vs:Vs≤Cbvs},
where Cbvs is defined in (23). 𝒮 can be expressed in terms of Vs, on the basis of the definition of Vs in (40):
(52)S=2Vs.
In the following, we analyze the convergence of 𝒮. Taking into account the definition of Cbvs in (23), the fact that Vs converges to Dvs defined in (51), and (52) the convergence of 𝒮 is:
(53)lim⁡t→∞S=Ds,Ds={S:|S|≤Cbs},Cbs=2Cbvs=λn-1Cbe.
To analyze the convergence of e, we begin by expressing e in terms of 𝒮, according to (12):
(54)e=1(p+λ)n-1S.
The convergence of the tracking error may be established on the basis of (53) and (54) [18, 21]:
(55)lim⁡t→∞⁡e(t)=De,De={e:|e|≤Cbe},
which means that the tracking error e converges to a residual set whose size is of the user’s choice.


## 6.1. Theorem: Boundedness and Tracking Convergence

If the controller designed in Section5 is applied to the plant (6), then the tracking error e(t) converges to De asymptotically, and the closed-loop signals remain bounded.Proof.
Now we proceed to analyze the stability of the controlled system using the direct Lyapunov method and the Barbalat’s Lemma. First, we establish the boundedness of the closed-loop signals on the basis of the time derivative of the Lyapunov function. Then, we establish the convergence of the tracking error to the target regionDe defined in (11), on the basis of the Barbalat’s Lemma. The validity of the procedure can be stated using the theory developed and applied in [18, 19, 21, 22].
The closed-loop dynamics is given by (32), and (38). We define the following truncated quadratic form, on the basis of the truncation presented in [19]:
(39)V̅s={VsifVs≥Cbvs,CbvsifVs<Cbvs,(40)Vs≜(12)S2,
where Cbvs is defined in (23). The form Vs̅ has the following properties:
(41)-Vs<-Vs̅+Cbvs<0ifVs>Cbvs,-Vs<-Vs̅+Cbvs=0ifVs=Cbvs,-Vs≤0=-Vs̅+Cbvs=0ifVs<Cbvs.
The states of the closed-loop system may be grouped in the following vector:
(42)x̅=[S,c̃,C̃d⊤,θ̃⊤]⊤,
where c̃, C̃d, and θ̃ are defined in (26), (37), (27). We define the following function:
(43)V(x̅(t))=Vs̅+Vθ+Vc+Vd,Vθ=(12)θ̃⊤Γ-1θ̃,Vc=(12)γc-1c̃2,Vd=(12)C̃d⊤Γd-1C̃d.
Notice that V(x̅) does not satisfy standard conditions of Lyapunov theory for time-variable systems, shown in [18]. It satisfies the conditions of [19] instead, which we use to prove the boundedness of the closed-loop signals and the convergence of the tracking error to the residual set De (11). The time derivative of Vs along the trajectory (32) is
(44)V̇s=SṠ=-amS2-|b|ĉ2S4-|b|(φ⊤θ̂)2S2+(-1)|b|ĉ12f12S2+⋯+(-1)|b|ĉl2fl2S2+S(d-φ⊤θ*)≤-amS2-c̃S2+Sφ⊤θ̃+(-1)|S|C̃d⊤fifVs≥Cbvs,
where the last inequality is obtained from (36). The time derivative of Vs̅ can be derived from (39) and expressed on the basis of the above equation:
(45)Vs̅̇={V̇s≤-amS2-c̃S2+Sφ⊤θ̃+(-1)|S|C̃d⊤fifVs≥Cbvs,0ifVs<Cbvs.
The time derivative V̇θ+V̇c+V̇d along trajectories (38) is
(46)V̇θ+V̇c+V̇d={(-1)Sθ̃⊤φ+c̃S2+C̃df|S|ifVs≥Cbvs,0otherwise.
Therefore, V̇ is given by
(47)V̇=Vs̅̇+V̇θ+V̇c+V̇d={V̇s+V̇θ+V̇c+V̇d≤-amS2ifVs≥Cbvs,0ifVs<Cbvs.
The above equation can be rewritten in terms of Vs defined in (40):
(48)V̇≤-2amVsifVs≥Cbvs,V̇=0ifVs<Cbvs.
This implies that V̇≤0 for all t≥to. Thus, according to [19], all the closed-loop signals are bounded, that is, (𝒮, θ̃, c̃, C̃d)∈L∞, or equivalently, x̅∈L∞, where x̅ is defined in (42). As a consequence, (V,  V̅s, Vs, Vθ, Vc, Vd)∈L∞, according to the definitions in (39) and (43). The boundedness of x̅ implies (e, …, e(n-1))∈L∞ according to (12); hence, 𝒮̇∈L∞ [22], and V̅̇s∈L∞. To establish the convergence of the tracking error, we begin by expressing the equation ahead in terms of Vs̅, according to the definition of Vs̅ in (39) and the properties of (41):
(49)V̇≤-2amVs<-2amVs̅+2amCbvs≤0ifVs≥Cbvs,V̇=0=2am(-Vs̅+Cbvs)=0ifVs<Cbvs,⟹V̇≤-2amVs̅+2amCbvs≤0.
We reorganize the above equation and integrate as follows:
(50)2am∫tot(Vs̅-Cbvs)dτ+V(t)≤V(to).
Thus, (Vs̅-Cbvs)∈L1. Recall that (Vs̅,Vs̅̇)∈L∞. Thus, (Vs̅-Cbvs)∈L∞∩L1, (d/dt)(Vs̅-Cbvs)∈L∞. By invoking the Barbalat Lemma [23], we obtain (Vs̅-Cbvs)→0 as t→∞. In turn, this implies that Vs→Dvs, being Dvs defined as
(51)Dvs={Vs:Vs≤Cbvs},
where Cbvs is defined in (23). 𝒮 can be expressed in terms of Vs, on the basis of the definition of Vs in (40):
(52)S=2Vs.
In the following, we analyze the convergence of 𝒮. Taking into account the definition of Cbvs in (23), the fact that Vs converges to Dvs defined in (51), and (52) the convergence of 𝒮 is:
(53)lim⁡t→∞S=Ds,Ds={S:|S|≤Cbs},Cbs=2Cbvs=λn-1Cbe.
To analyze the convergence of e, we begin by expressing e in terms of 𝒮, according to (12):
(54)e=1(p+λ)n-1S.
The convergence of the tracking error may be established on the basis of (53) and (54) [18, 21]:
(55)lim⁡t→∞⁡e(t)=De,De={e:|e|≤Cbe},
which means that the tracking error e converges to a residual set whose size is of the user’s choice.


## 7. Application to a Permanent Magnet Synchronous Motor

A PMSM is a kind of highly efficient and high-powered motor. The benefits of the PMSM are discussed in [24, 25]. The PMSM has the following difficulties [26]: (i) it is highly nonlinear; (ii) the parameters of the physical model experience unknown time-varying behavior, for example, the stator resistance R and the friction coefficient B; (iii) unknown external disturbances appear, that is, the load torque disturbance TL. Moreover, by varying the permanent magnet flux λaf the state ω exhibits a pitchfork bifurcation [27]. All the mentioned facts lead to response deterioration of many controllers, specially for high-speed and high-precision tasks in real applications [28].In view of the complex behavior, with parameters varying with respect to time and state plane, the scheme developed in this work is suitable. We apply the developed scheme by simulation, to a PMSM whose model and parameters are presented in [29]. Therein, the possible manipulated inputs are ud and uq; and id, iq, and ω are possible outputs. We choose the motor angular frequency ω as the output to be controlled and the d-axis voltage ud as the control input, that is, y=ω, u=ud, υd=0. A similar choice is made in [24], where an adaptive controller is derived for a PMSM. After performing a step response analysis, the variable ω exhibits a second-order behavior when a step change is introduced in ud. This second order behavior was already noticed in [29]. Thus, we use the regression model (6) with n=2. We can summarize the basic features of the plant as(56)n=2,sgn⁡(b)=+1,l=3,f1=|ẏ|,f2=|y|,f3=1,⟹f=[|ẏ|,|y|,1]⊤,ϕ=[y2,y3]⊤,φa=-ÿd+λoė,⟹φ=[-ϕ⊤,-φa-amS]⊤.
We used a factorization with terms y2 and y3 to take into account the presence of the pitchfork bifurcation. According to [30] the normal form of the pitchfork bifurcation gives a description of the system behavior in a tight neighborhood of the bifurcation point. For the case of the pitchfork bifurcation, the normal form includes the terms y2 and y3. Although the system usually works in different regions far from the bifurcation point, we wanted to include the behavior in the neighborhood of the bifurcation point. Thus, we included the terms of the normal form for the pitchfork bifurcation [30], that is, y2 and y3. The terms 1 and y are not included in ϕ because they are already present in d. The remaining parameters of the controller are defined on the basis of the parameters in (56):(57)e=y-yd,S=(p+λ)e=ė+λe⟹λo=λ,Cbvs=(12)(λCbe)2,u=(-1)sgn⁡(b)ĉ2S3+(-1)sgn⁡(b)(φθ̂)2S+(-1)sgn⁡(b)S(ĉ12f12+ĉ22f22+ĉ32f32),=(-1)sgn⁡(b)S(ĉ2S2+(φθ̂)2+ĉ12f12+ĉ22f22+ĉ32f32).
In addition, we choose:(58)λ=3,am=10,am,1=70,am,o=1225(for the reference model),Cbe=0.1,γc=20,Γd=diag⁡[20,20,20],Γ=diag⁡[0.4,0.4,0.003].The results are shown in Figure1. The simulation shows expected results: all the closed-loop signals remain bounded, and transient values of the tracking error remain in a small interval. It is worthy of note that the transient values of e depend on its initial value and the initial values of the adjustment errors. To show the effect of the time-varying behavior of the model coefficients, we change the stator resistance R from 1.4 to 1.7 Ω at the time instant 0.6 sec. Results are shown in Figure 2.Transient behavior of the output and control input.
(a)(b)Performance of the tracking error and control input varying the resistance value.
(a)(b)In addition, we consider the variation of the damping constantB from 0.00038818 to 0.00046582 Nm/(rad/s) at the instant 0.6 sec. Results are shown in Figure 3. Notice that in the three cases the final value of |e| is less than Cbe=0.1. Moreover, the control input u belongs to the interval [-200200] V. Notice in Figures 2 and 3 that the effect of the disturbance on the tracking error is almost negligible. Nevertheless, the control input experiences a large variation, as can be seen in Figure 3.Performance of the tracking error and control input varying the damping constant value.
(a)(b)


## 8. Conclusions

In this work, we have proposed a control scheme for highly nonlinear plants, based on a simple plant model with polynomial approximators, which provides an adequate description of transient behavior. It is worth noticing the fact that benefits Ri to Rviii (Section2) are achieved at the same time, with minimal requirements on the knowledge of the plant. Indeed, the relative degree and the signum of the control gain can be established by a previous step response analysis. Many techniques are combined at the same time: sliding surface MRAC, dead zone-type update law, robustifying auxiliary control, approximation techniques, and truncation of the quadratic forms.The disadvantage of polynomials as approximators is that they may be less accurate than other techniques, for example, neural networks or fuzzy sets, leading to higher approximation error. Since the approximation error is bounded, it can be handled by means of robustness techniques without requiring the upper bound to be known. Moreover, we considered the coefficients of the termsy,…,y(n-1) as time-varying but bounded, with constant and unknown bounds. Then, we handled this by means of robust control, without requiring the upper bounds to be known.We handled the time-varying behavior of the control gain by means of robustness techniques, without using the Nussbaum gain method. The redefinition of the plant terms in terms of adjustment errors and adjustment parameters is a fundamental step. The resulting expression allows a straightforward definition of the control law. The variation of the control gain implies that the terms involving adjusted parameters cannot be cancelled. Rather, we attenuate its effect by means of squared terms and handle the residual error by means of an additional control term that is only a function of the sliding surface. The resulting expression forV, V̇, and the design is simpler in comparison with the Nussbaum technique.



---

*Source: 269346-2012-06-06.xml*