Robust Output Feedback Controller Scheme For A Class of Uncertain Nonlinear Systems
Robust Output Feedback Controller Scheme For A Class of Uncertain Nonlinear Systems
J. Kuvulmaz* E. Zergeroglu**
Abstract: In this study, we present a new continuous output feedback type controller mechanism for the
tracking problem of a class of uncertain nonlinear systems. The proposed strategy requires the
uncertainties of the dynamical system to be first order differentiable and achieves semi-global asymptotic
tracking when only the system outputs are measurable. The Controller design is based on a Lyapunov-type
stability argument. Simulation studies on a two link planar robotic system are presented to illustrate the
feasibility of the proposed strategy.
the stability result was semi-global2 asymptotic. On the other differentiable nonlinear, uncertain function and u (t ) is the
hand, a limitation that exists in almost all of the proposed control input.
OFB link position tracking controllers is the semi-global
nature of the stability results. To give a few examples, a Property 1: The inertia matrix can be upper and lower
model-based observer was used to construct a semi-global bounded by the following inequalities (Lewis et al., 1993)
exponential link position tracking controller (Lim et al.,
1996). Variable structure OFB controllers were designed to m1 ≤ m ( x) ≤ m 2 (2)
compensate for parametric uncertainty (Canudas de Wit and where m1 and m2 are positive constants.
Fixot, 1991; Canudas de Wit and Slotine 1991). Filter-based
robust control schemes which produce semi-global, Our control objective is to ensure that the state signal x(t )
uniformly ultimately bounded (UUB) were designed in would track the given smooth reference trajectory x d (t ) . To
(Berghuis and Nijmeijer, 1994; Qu et al., 1995; Yuan and
Stepanenko, 1991). Then adaptive OFB controllers were quantify this objective we defined the tracking error signal
presented in (Burg et al.,1994; Burg et al., 1996; Kaneko and e(t ) in the following form
Horowitz, 1997; Zergeroglu et al., 1999) that yield semi-
e = xd − x (3)
global asymptotic link position tracking in the presence of
parametric uncertainty. In 1996, Loria developed a model- In our analysis we will utilize the common assumption that
based controller which produces global uniform asymptotic the reference trajectory signal x d (t ) and its first three time
tracking but the proposed method is only valid for a one derivatives are always bounded (i.e.
degree-of-freedom (DOF) system. Then, a global OFB x d (t ), x& d (t ), &x&d (t ), &x&&d (t ) ∈ L∞ ).
adaptive controller was designed for n-DOF robot
manipulator (Zhang et al.,2000) Finally, a robust OFB Velocity variables can not be measured because of that the
tracking controller was proposed with a global, uniformly related variables are obtained by filtering technique of
ultimately bounded link position tracking (Dixon et al.,2004). position error. While position tracking error is the input, the
In this work, we extend our previous result (Kuvulmaz and velocity tracking error is the output of the system. The
Zergeroglu, 2007) to output feedback case. When compared proposed filter dynamics is as follows
to (Ortega et al., 2002) our approach does not possess any r f = p − (k1 + 1) e (4)
singularities and when compared to (Xian et al., 2004) can
compensate the uncertainties for a larger class of nonlinear p& = − p + e − e f − ( k1 + 1) r f (5)
systems. However, our approach requires a high gain e& f = − e f + r f (6)
condition on the feedforward compensation gain which might
be considered as a theoretical weakness. where the filter output r f (t ) ∈ℜ will be used for the link
velocity variable. p(t ), e f (t ) ∈ ℜ are the auxiliary variables
The remaining of the paper is organized as follows: the
model under consideration and the control problem are stated used to establish the velocity variable. k1 is the positive
in Section 2. Section 3 contains the error system control gain which is defined as
development, while controller design with the stability
analysis to ensure asymptotic tracking and boundedness of 1
k1 = (k n + 1) (7)
the closed loop system are given in section 4. Simulations m1
performed on a two-link, planar robotic mechanism are where k n is the positive nonlinear damping gain. In order to
presented in Section 5 and lastly some concluding remarks
with possible future research are presented in Section 6. form the open loop position error system, dynamics of the
filter output would be obtained. Taking the time derivative of
2. PROBLEM STATEMENT (4) and inserting for p& from (5) we get the dynamics of the
filter output
For the ease of presentation3, we consider a second order
single-input single-output dynamical system having the r& f = − (k1 + 1)e& − p + e − e f − (k1 + 1) r f (8)
following form Inserting the value of p from (4) in equation (8),
m( x) &x& + f ( x, x& , θ ) = u (t ) (1)
r&f = − (k1 + 1)η + e − e f − r f (9)
where x is the output and the notations x& , &x& are used to
identify the first and second order derivatives of the output where the filtered error signal η (t ) ∈ ℜ is defined as
respectively. m( x) 〉 0 and f ( x, x& , θ ) is first order η = e& + e + rf (10)
If (10) is rearranged, we get the error dynamics as
e& = − e − r f + η (11)
2
In a semi-global stability result, a control gain often has to
In addition, we define an auxiliary term named integral effect
be adjusted according to the “size” of the initial conditions.
3 injection term ξ as
Extension to multi-input multi output and higher order
versions are also possible with a considerably small effort.
14126
17th IFAC World Congress (IFAC'08)
Seoul, Korea, July 6-11, 2008
At this point we define the desired version of the auxiliary wI :=η ( N d − β tanh(e + e f )) (25)
signal N , N d such that where the auxiliary constant term δ b , explicitly given in the
Nd := N (17) following form
x& = x& d , x = x d
δ b := β (ln cosh(e(0) + e f (0)) +1)
Note that due to assumption that the reference trajectory term (26)
x d is third order differentiable, the newly defined “desired” − (e(0) + e f (0)) N d (0).
version of the auxiliary term can be proven to be at least first When the constant scalar control gain β of (26) is selected to
order differentiable (i. e. N d (t ), N& d (t ) ∈ L∞ ) . Adding and satisfy
subtracting N d to the right hand side of (17), we obtain
β 〉 k Nd [ ∞
+ N& d
∞
] (27)
~ 1
m η& = − k1 mη + N + N d − u − m& η (18) with the high but bounded design gain k defined as
2
~ ⎧⎪ ⎫⎪
where the function N is defined as 1
k = max ⎨(1 + ε ), ⎬ (28)
⎪⎩ tanh( w)
~ 1 w≠ 0 ⎪⎭
N := N − N d − m (2 rf + e f ) + m& η (19)
2 where ε 〉 0 . Then β I (t ) will always be lower bounded by
Remark 1: Since the auxiliary function N defined in (16) is
~ zero (i. e. β I ≥ 0 ) or equivalently
continuously differentiable, we can show that N can be upper
t
bounded in the following manner:
~ ∫ w (τ ) dτ ≤ δ
I b
(29)
N ≤ ρ ( j) j (20) 0
and ρ (⋅) is a positive defined non-decreasing bounding To prove the Lemma, we take only the part inside the integral in
(24) and substitute (30) to obtain
function.
14127
17th IFAC World Congress (IFAC'08)
Seoul, Korea, July 6-11, 2008
t t
and s (t ) is the vector function defined as
∫ w (τ ) dτ = ∫ [ w (τ ) ( N
I d (τ ) − β tanh( w (τ )))] dτ
0
t
0
d w (τ ) (32)
[
s = jT z ξ ]
T
(39)
+ ∫
0
dτ
N d (τ ) dτ
Proof: To prove the Theorem we define a non negative
t
function of the form
d w (τ )
− β ∫
0
dτ
tanh ( w (τ ))dτ . 1 1
V ( t ) = mη 2 + Z + k i ξ 2 (40)
2 2
evaluating the third term and integrating the second term on
the right hand side by parts, we have with Z (t ) is defined to have the form
t t
⎡ dN d (τ ) ⎤ Z := β P + β I (41)
∫ wI (τ ) dτ =
0
∫ w (τ ) ⎢⎣ N
0
d (τ ) −
dτ
− β tanh( w(τ ) ⎥ dτ
⎦ where β P is selected as
+ [ w (t ) N d (t ) − β (ln (cosh ( w(t ))) + 1 ] (33)
1
1 1
+ [ β (ln (cosh ( w ( 0 ))) + 1) − w ( 0 ) N d ( 0) ] , β P = e 2 + e 2f + r f2 (42)
22 2
where the β term has been added and subtracted from the
and β I term is the lower bounded function defined in (24).
right hand side. Now we can upper bound (33) in the
following way We can upper and lower bound (40) by using λ1 and λ 2 terms
in the following way
t t
⎡ dN (τ ) ⎤
∫0 I
w (τ ) dτ ≤ ∫0 w(τ ) ⎢⎢⎣ Nd (τ ) + d (dτ ) − β tanh( w(t)) ⎥⎥⎦ dτ 2
λ1 s ≤V ≤ λ2 s
2
(43)
[
+ w (t ) Nd (t ) − β ( ln(cosh ( w (t ))) +1) ] (34) Taking the time derivative of (40) substituting (7) and (23)
we obtain
+ [β ( ln(cosh ( w (0))) +1) − w (0) Nd (0)]
~
V& = − e2f − rf2 − e2 −η 2 − kn η 2 +η N (44)
Where the w (t ) tanh( w (t )) ≥ 0 has been utilized. Notice that And from (20), we can state the following upper bound for
in (34), the integral term (the first line) on the right hand side the time derivative of function V given in (40)
of the inequality will exactly be zero when w(t ) = 0 and will
have negative values for all other values of w(t ) when the
2
[
V& ≤ − j + η ρ ( j ) j − kn η 2 ] (45)
2
controller gain β is selected to satisfy (27). Similarly the ρ 2 ( j) j
Adding and subtracting term to the right hand
same selection of β will also ensure the negative semi 4 kn
definiteness of second term on the right hand side of (34) (as side of (45) yields
k (ln (cosh ( w)) + 1) − w ≥ 0 for k ≥ (1 + ε ) ). And from the 2
2 ρ 2 ( j) j
definition of δ b given in (26), it is straight forward to show V& ≤ − j +
4 kn
that (29) holds. (46)
⎡ ρ 2 ( j) j ⎤
2
14128
17th IFAC World Congress (IFAC'08)
Seoul, Korea, July 6-11, 2008
14129
17th IFAC World Congress (IFAC'08)
Seoul, Korea, July 6-11, 2008
14130