• 検索結果がありません。

Robust Output Feedback Disturbance Attenuation of Nonlinear Uncertain Dynamic Systems

N/A
N/A
Protected

Academic year: 2021

シェア "Robust Output Feedback Disturbance Attenuation of Nonlinear Uncertain Dynamic Systems"

Copied!
19
0
0

読み込み中.... (全文を見る)

全文

(1)

ISSN 1344-8803, CSSE-15 June 6, 2001

Robust Output Feedback Disturbance Attenuation of Nonlinear Uncertain Dynamic Systems

via State-Dependent Scaling

¶§

Hiroshi Ito

†‡

Department of Control Engineering and Science, Kyushu Institute of Technology 680-4 Kawazu, Iizuka, Fukuoka 820-8502, Japan

Phone: (+81)948-29-7717, Fax: (+81)948-29-7709 E-mail: [email protected]

Abstract: This paper presents a novel approach to the problem of output feedback stabilization with L

2

disturbance attenuation for nonlinear uncertain systems. A new method of state-dependent scaling is introduced into the output feedback design, which unifies treatment of nonlinear and linear gains. The effect of disturbance on the con- trolled output, which is allowed to be any function of measurement output, can be attenuated to an arbitrarily small level with global asymptotic stability if the plant belongs to a wide class of interconnected systems whose uncertain components unnec- essarily have finite linear-gain. The uncertain dynamics is not limited to input-to-state stable systems either. The approach is not only a natural extension of popular ap- proaches in robust linear control, but also advantageous to numerical computation.

The design procedure proposed in this paper consists of novel recursive calculation of robust observer gain as well as feedback gain.

Keywords: State-dependent scaling, Nonlinear gain, Output feedback, Robust global stabilization, Almost disturbance decoupling, Robust control Lyapunov function

Technical Report in Computer Science and Systems Engineering, Log Number CSSE-15, ISSN 1344-8803.

c

°2001 Kyushu Institute of Technology

§The current version of the paper was completed by May 1, 2001. A primary version was presented at SICE 1st Annual Conference on Control Systems, Kyoto, Japan, May 23, 2001. A brief version of the result has been submitted to IEEE Conference on Decision and Control 2001.

Author for correspondence

(2)

1 Introduction

Recently, an idea of state-dependent scaling has been introduced into dissipativity-based robust control for nonlinear systems in [4, 5]. The primitive state-dependent scaling[4, 5, 6] aimed at nonlinear uncertain systems whose uncertainties were described as L

2

-gain balls, while the known part of the system was allowed to have infinite L

2

-gain. If we have more knowledge on the un- certainty such as nonlinear gain[12, 9, 13, 8] other than linear gain or if the uncertainty does not exhibit finite linear-gain properties, design results based on [4, 5, 6] might be too conservative.

As for output feedback control based on state-dependent scaling, the design procedure presented in [6] did not address the existence of globally stabilizing controllers in the presence of dynamic uncertainties although the existence was proved for static uncertainties. The problem of distur- bance attenuation was not addressed either. Although the recursive observer design proposed in [6] succeeded in extending the output-feedback form defined in [10] to a slightly wide class of systems, it would be possible to cover a more general class of systems by fully exploiting the unique nested structure of robust observer design.

In the last decade, L

2

-gain disturbance attenuation with global internal stabilization using full- state information has been extensively studied for linear and nonlinear systems(See [2] and ref- erences therein). In comparison, when only the output feedback is allowed, the problem is less understood. For essentially nonlinear systems, filtered transformation and backstepping technique were employed in [11] to solve a problem of output feedback tracking with almost disturbance decoupling. Another approach proposed in [1] resorts to solutions of Hamilton-Jacobi partial dif- ferential inequalities and a coupling condition. The recent work [7] has considered a more relaxed class of nonlinear systems than earlier work, and the plant is allowed to involve unmeasured dy- namics which is input-to-state stable. The layout of interconnected uncertain systems for which this paper will give a new characterization of the output feedback disturbance attenuation and the existence of solutions is broader than setups considered in those earlier papers. This paper will also allow systems to have both static and dynamic uncertainties including systems which are not input-to-state stable. Nonlinearities allowed in the plant by this paper is more general than those in [10, 6, 7].

The purpose of this paper is to develop a new method of state-dependent(SD) scaling design in

order to achieve the output feedback disturbance attenuation with global asymptotic stability

for nonlinear systems described by interconnection of nonlinear-gain bounded systems. Thereby,

the use of linear gain and nonlinear gain is unified. The design becomes a natural extension of

popular techniques in linear robust control to nonlinear systems. The development is considered as

a global robustification of the previous results[6] against dynamic uncertainties and its extension

to L

2

disturbance attenuation. Dynamic and static uncertainties are treated in a unified way so

that design formulas for the two types of uncertainties are identical. A difference only appears

in classes from which scaling factors are chosen. The state-dependent scaling characterization

does not require systems to fit in some geometric structure. For interconnected uncertain systems

partially in an extended feedback form, the control laws can be systematically generated by

selecting parameters of the observer and the feedback gain recursively. The recursive procedure

proposed in this paper not only allows us to use nonlinear gains, but also brings in a unique way of

(3)

constructing robust observers which enable the output feedback to make the effect of disturbance arbitrarily small with respect to generalized equations of controlled output. The new procedure of observer design proposed in this paper is also unique in that it can be applied to a broader family of nonlinearities in the plant, compared with previous output feedback results in [10, 6, 7]. The design equations are obtained as affine algebraic inequalities with respect to the design parameters, so that the SD scaling approach is advantageous to systematic numerical computation as well as analytical computation.

2 System description

Consider the uncertain nonlinear system Σ

0

:



˙

x = A(y)x + B(y) ¯ w + G(y)u , x(t) ∈ R

n

¯

z = C(y)x , z(t), ¯ w(t) ¯ ∈ R

p+q

y = C

y

x , y(t) ∈ R (1)

The matrices A, B, G, and C are assumed to be C

0

functions of y, and C

y

is a constant row vector. Scalars u(t) and y(t) are control input and measurement output, respectively. The signals

¯

w and ¯ z are partitioned as

w =



w

1

w

2

w

m



, z =



z

1

z

2

z

m



,

w

i

(t), z

i

(t) ∈ R

pi

p

i

0 p =

Pmi=1

p

i

r =



r

1

r

2

r

m



, e=



e

1

e

2

e

m



,

r

i

(t), e

i

(t) ∈ R

qi

q

i

0 q =

Pmi=1

q

i

¯ w

i

=

·

w

i

r

i

¸

, z ¯

i

=

·

z

i

e

i

¸

∈ R

pi+qi

, w ¯ =



¯ w

1

¯ w

2

¯ w

m



, z ¯ =



¯ z

1

¯ z

2

¯ z

m



Suppose that there is a system Σ

described by the following nonlinear mappings between z to w.

i

: z

i

=

·

z

is

z

id

¸

7→ w

i

=

·

w

is

w

id

¸

, w

i

=

·

is

0 0 ∆

id

¸

z

i

. (2)

Here, ∆

is

and ∆

id

represent a time-varying static system and a time-varying dynamic system, respectively. It is unnecessary for ∆

i

to have the both types. These systems are defined by

is

: w

is

= h

is

(z

is

, t) (3)

id

:

½

x ˙

i

= f

id

(x

i

, z

id

, t)

w

id

= h

id

(x

i

, t) (4)

Assume f

id

(0, 0, t) = 0, h

id

(0, t) = 0 and h

is

(0, t) = 0 for all t 0. Functions f

id

, h

i∗

(∗

stands for s or d) are locally Lipschitz in (x

i

, z

i∗

) on R

ni

× R

pi∗

, uniformly in t ∈ R. The

state variable of Σ

is x

= [x

T1

, · · · , x

Tm

]

T

∈ R

n

. The system Σ

represents uncertainty so

that knowledge of f

id

, h

id

and h

is

is unnecessary. We only assume that information about

(4)

nonlinear or linear gain is available in the sense described in Section 3. The interconnected system consisting of Σ

0

and Σ

is denoted by Σ

P

.

Since the state variable x is supposed to be unmeasurable, we employ the following observer to control the uncertain system Σ

P

.

½

x ˙ˆ = A(y)ˆ x + Y (y, x)(y ˆ y) + ˆ G(y)u , x(t) ˆ ∈ R

n

ˆ

y = C

y

x ˆ , y(t) ˆ ∈ R (5)

In this paper, given an arbitrary number τ > 0, we seek the output feedback control consisting of (5) and

u = K(y, x)ˆ ˆ x . (6)

which

globally uniformly asymptotically stabilizes Σ

P

when r 0

makes the mapping between r and e have L

2

-gain less than or equal to τ

The state variables x and x

are not measured for the feedback control. Functions Y and K are C

0

functions which have yet to be determined. The system Σ

P

is said to be globally uniformly asymptotically stabilized if the equilibrium x

cl

= [x

T

, x

T

, x ˆ

T

]

T

= 0 is globally uniformly asymp- totically stable. In this paper, the system Σ

P

is said to have L

2

-gain less than or equal to τ if there exists a storage function V (x

cl

) which is positive definite and radially unbounded such that for all initial states x

cl

(0) ∈ R

2n+n

, and all r ∈ L

2

[0, T ], the inequality

V (x

cl

(T)) V (x

cl

(0)) +

Z T

0

2

krk

2

− kek

2

)dt holds for all T 0.

3 Nonlinearly bounded uncertainty

In this paper, the uncertainty Σ

is supposed to belong to the following class of nonlinearly bounded systems.

Assumption 1 For each i = 1, 2, . . . , m, the uncertain system Σ

satisfies the following.

(a) There exists a C

0

function ψ

is

: [0, ∞) [0, ∞) such that

kw

is

k

2

ψ

is

(kz

is

k)kz

is

k

2

(7) holds for all t [0, ∞).

(b) There exists a C

0

function ψ

id

: [0, ∞) [0, ∞) and a C

1

function W

∆i

: [0, ∞) × R

ni

→ R such that

β

i

(kx

i

k) W

∆i

(t, x

i

) β ¯

i

(kx

i

k) (8)

∂W

∆i

∂t + ∂W

∆i

∂x

i

f

id

≤− β

i

(x

i

)−kw

id

k

2

id

(kz

id

k)kz

id

k

2

(9)

(5)

hold for all (t, x

i

, z

id

) [0, ∞) × R

ni

× R

pid

, where β

i

and β ¯

i

are class K

functions, and β

i

is a positive definite C

0

function of x

i

.

A system Σ

is said to be admissible if Assumption 1 is true. The assumption does not require uncertain systems to have finite L

2

-gain. Instead, they are supposed to have finite nonlinear-gain.

When ∆

is

(∆

id

) exhibits finite L

2

-gain, the parameter ψ

is

id

, respectively) reduces to a positive constant. In such a case, we obtain ψ

is

= γ

is2

and ψ

id

= γ

id2

, where γ

is

and γ

id

are L

2

-gain. The new class of uncertain systems is broad and it includes input-to-state stable(ISS) systems[12] in the following sense.

Lemma 1 (a) Suppose that a static system

is

admits class K

functions α

i

and σ

i

such that

α

i

(kw

is

k) σ

i

(kz

is

k) (10)

holds for all t [0, ∞) and

s→0

lim

+

σ

i

(s)

α

i

(s) < +∞ (11)

holds. Then, there exists a C

0

function ψ

is

such that (7) holds for all t [0, ∞).

(b) Suppose that a dynamic system

id

admits a C

1

function V

∆i

: [0, ∞) × R

ni

→ R such that α

i

(kx

i

k) V

∆i

(t, x

i

) α ¯

i

(kx

i

k) (12)

∂V

∆i

∂t + ∂V

∆i

∂x

i

f

id

≤ −α

i

(kx

i

k) + σ

i

(kz

id

k) (13) are satisfied for all (t, x

i

, z

id

) [0, ∞) × R

ni

× R

pid

where α

i

, α ¯

i

and α

i

are class K

functions and σ

i

is a class K function and they satisfy

kx

lim

ik→0

kw

id

k

2

α

i

(kx

i

k) < +∞, lim

kzidk→0

σ

i

(kz

id

k)

kz

id

k

2

<+∞ (14)

uniformly in t. Then, there exists a C

0

function ψ

id

, a C

1

function W

∆i

, class K

functions β

i

and β ¯

i

and a positive definite C

0

function β

i

such that (8) and (9) hold for all (t, x

i

, z

id

) [0, ∞) × R

ni

× R

pid

.

The functions ψ

is

and ψ

id

are obtained easily from α

i

, σ

i

α

i

and ¯ α

i

[3]. Note that β

i

(x

i

) can be always chosen as a class K

function of kx

i

k for ISS systems defined in Lemma 1(b). It is emphasized that Assumption 1 admits systems which are not ISS. An example of nonlinearly bounded static mappings which violate (11) is w

is

= h

is

(z

is

) =

p

kz

is

k which is not Lipschitz at z

is

= 0. Indeed, if h

is

(z

is

, t) is Lipschitz at z

is

= 0 uniformly in t as assumed in Section 2, there always exists a class K

pair of {α(s), σ(s)} satisfying (11) and (10). As for a dynamic nonlinear mapping ∆

id

: z

id

7→ w

id

, the condition (14) together with (13) is common in asymptotic analysis based on the nonlinear small-gain technique[7, 8, 9]. It is known that (14) is always satisfied for appropriate functions α

i

∈ K

and σ

i

∈ K if the Jacobian linearization of ∆

id

at x

i

= 0 is uniformly asymptotically stable.

Example 1 An example of admissible uncertain dynamics z

id

7→ w

id

is

½

x ˙

i

= −x

i

(1 z

idρ

)

w

id

= sat(x

i

)

(6)

with ρ 2. Indeed, it satisfies (9) for ψ

id

= 2.2|z

id

|

ρ−2

and W

i

=

Z x2

i

0

1.1

s ds, β

i

= x

2i

5(x

2i

+ 1)

Clearly, the system is not input-to-state stable although it is globally asymptotically stable when z

id

0.

Example 2 The following input-to-state stable system

(

x ˙

i

= −x

3i

+ z

2id

w

id

= x

ρi

which is used in [7] with ρ 2 is also an admissible dynamic system. The functions ψ

id

and β

i

are obtained as

ψ

id

= τ π

ρ/2−1

3

ρ/2

2

−2ρ/3

|z

id

|

4ρ/3−2

β

i

= (1 τ

−1

π

−1

|x

i

|

for any τ > 1 and π > τ /(τ 1).

4 SD scaling characterization

This section derives a characterization of global robustness properties of the output feedback system described in Section 2 via a new concept of state-dependent (SD) scaling which incorporates the nonlinear gain.

First, a set of scaling factors associated with static uncertain components ∆

is

is defined by Φ

is

=

½

Φ

is

(y, x ˆ ) = φ

is

(y, x)I ˆ : φ

is

(·) ∈ C

0

φ

is

(y, x) ˆ >0, ∀(y, x) ˆ ∈ R

n+1

¾

(15) The identity matrix I is compatible in size with z

is

. The scaling factors are functions of the output and the state estimate. For dynamic uncertain components ∆

id

, a set of scaling factors is defined by

Φ

id

=

½

Φ

id

= φ

d

·

φ ˇ

id

I 0 0 I

¸

: φ

d

, φ ˇ

id

>0

¾

(16) The block partition of Φ

id

is compatible in size with that of [z

Tid

, e

Ti

]

T

. All sets Φ

id

, i = 1, 2, . . . , m are defined with a common constant φ

d

. For i = 1, 2, . . . , m, define Φ

i

(x) as

Φ

i

=

½

Φ

i

(y, x) = ˆ

·

Φ

is

(y, x) 0 ˆ 0 Φ

id

¸

: Φ

is

Φ

is

Φ

id

Φ

id

¾

Using C

0

functions ψ

id

, ψ

is

: [0, ∞) [0, ∞) in Assumption 1, define ¯ Ψ(x) as Ψ(x) = ¯ block-diag

m

i=1

Ψ ¯

i

(x) (17)

Ψ ¯

i

(x) =

ψ

is

(kz

is

k)

1/2

I 0 0 0 ψ

id

(kz

id

k)

1/2

I 0

0 0 τ

−1

I

(18)

(7)

The block diagonal structure of ¯ Ψ

i

is conformable in size to the partition ¯ z

i

= [z

isT

, z

idT

, e

Ti

]

T

. The scalar τ is a positive number to describe the level of disturbance attenuation. We are now ready to define three sets of SD scaling matrices by

Φ =

½

Φ(y, x) = ˆ block-diag

m

i=1

Φ

i

(y, x), ˆ Φ

i

Φ

i

¾

(19) Θ =

½

Θ(y, x) : ˆ R

n+1

→ R

(p+q)×(p+q)

, Θ(·) ∈ C

0

Θ(y, x) ˆ > 0 ∀(y, x) ˆ ∈ R

n+1

¾

(20) Ψ =

½

Ψ(x, x) : ˆ R

2n

→ R

(p+q)×(p+q)

, Ψ(·) ∈ C

0

Ψ(x, x) ˆ 0 ∀(x, x) ˆ ∈ R

2n

¾

(21) All scaling matrices Φ, Θ and Ψ are ‘state-dependent’.

Based on the triplet of these scaling matrices, we shall characterize stability and L

2

distur- bance attenuation of Σ

P

. Consider a global diffeomorphism between [ˆ x

T

, x ˆ

T

x

T

]

T

∈ R

2n

and [ ˆ χ

T

, η

T

]

T

∈ R

2n

as follows:

·

χ ˆ η

¸

=

·

S(y, x) 0 ˆ

0 W

¸ ·

x ˆ ˆ x x

¸

(22) where W is a constant matrix. The time-derivative of ˆ χ is obtained as

˙ˆ χ = ∂S

∂y xC ˆ

y

x ˙ +

·

∂S

x ˆ

1

x, ˆ ∂S

x ˆ

2

x, , ˆ ∂S

x ˆ

n

x ˆ

¸

˙ˆ

x + S(y, x) ˙ˆ ˆ x

= X(y, x) ˙ ˆ x + T (y, x) ˙ˆ ˆ x Define

A ¯ =

h

C

yT

A

Ti

, W ˆ =

·

−Y

T

W

T

W

T

¸

, S ˆ =

·

S

−1

KS

−1

¸

, A ˆ = [A G]

Then, we obtain the following theorem.

Theorem 1 If there exist P > 0, P > ˜ 0 and scaling matrices Φ Φ, Θ Θ, Ψ Ψ such that

M (y, x) = ˆ



S ˆ

T

A ˆ

T

(X + T )

T

P + P (X + T ) ˆ A S P XB ˆ

B

T

X

T

P −Θ

ΦΨCS

−1

0

−W

−T

(XA + T Y C

y

)

T

P P W B ˜ S

−T

C

T

ΨΦ −P (XA + T Y C

y

)W

−1

0 −B

T

W

T

P ˜

−Φ −ΦΨCW

−1

−W

−T

C

T

ΨΦ W

−T

A ¯ W ˆ P ˜ + ˜ P W ˆ

T

A ¯

T

W

−1



< 0 (23)

Θ Φ (24)

Ψ ¯ Ψ (25)

hold for all (x, x) ˆ ∈ R

2n

, the output-feedback law (5-6) globally uniformly asymptotically stabilizes Σ

P

for all admissible uncertainties and Σ

P

has L

2

-gain less than or equal to τ .

In the case of {q = 0, φ

d

= 1, Θ = Φ, Ψ = I}, Theorem 1 reduces to the primitive result[6].

(8)

5 Recursive design

This section defines a class of systems Σ

P

and presents a recursive procedure of the output feedback control design for this class. Suppose that Σ

0

is in the following triangular structure.

y = x

1

, C

y

= [ 1 0 0 ] (26)

A(y) =







a

11

a

12

0 0

a

21

a

22

a

23

0 0 a

n−1,1

a

n−1,2

a

n−1,n

0

a

n1

a

n2

a

nn







, G (y) =



0 a

n,n+1

0



(27)

a

i,i+1

(y) 6= 0, 1 i n, ∀y ∈ R (28)

B (y) =

"

B

11

B

n1

#

, C(y) = [ C

11

0 0 ] (29)

where B

i1

(y) ∈ R

1×(p1+q1)

, C

11

(y) ∈ R

(p1+q1)×1

, m = 2n, p

1

= p, q

1

= q. In addition, we assume that

sup

y∈R

¯¯

¯¯

¯

a

2ij

(y) a

k−1,k

(y)a

l,l+1

(y)

¯¯

¯¯

¯

< +∞, 2 k n 1 k l n 1

k i n, k j n (30) sup

y∈R

a

nn

(y)

|a

n−1,n

(y)| < +∞ (31)

The conditions (30) and (31) will be used for ensuring the existence of global solutions to the observer design problem described later. In this paper, a system Σ

P

consisting of Σ

0

and Σ

which fulfill these structural assumptions (26-31) and Assumption 1, respectively, is said to be in the generalized robust output-feedback form. Compared with a standard output-feedback form defined in [10], the generalized robust output-feedback form not only allows for disturbance signals and dynamic uncertain components which unnecessarily have finite linear-gain, but also nonlinearity is not restricted to A(y)x = A

0

x + A

1

(y) where A

0

is a constant matrix. The class of generalized robust output-feedback form is also broader than an extended class considered in [6]. When the nonlinearity is limited to A(y)x = A

0

x + A

1

(y) + A

2

(y)x

2

with a constant matrix A

0

, the conditions (30) and (31) reduces to

sup

y∈R

¯¯

¯

a

2i2

(y)/a

12

(y))

¯¯¯

< +∞, 2 i n

y∈R

inf |a

12

(y)| 6= 0 which are assumptions employed in [6].

In order to solve the disturbance attenuation problem with robust stability for the above class of systems, we first pick any constant matrices P and ˜ P of the form

P = diag

n

i=1

P

i

> 0, P ˜ = diag

n

i=1

P ˜

i

> 0 (32)

Define

S

−1

(x

1

, x ˆ

[n−2]

) =





1 0 0 0

s

1

1 0 0

0 s

2

1 0

0 0 s

n−1

1





(33)

(9)

u = s

n

(x

1

, x ˆ

[n−1]

) ˆ χ

n

(34) Here, s

1

(x

1

), s

2

(x

1

, x ˆ

1

), ..., s

n

(x

1

, x ˆ

[n−1]

) are smooth functions to be determined. The notation

ˆ

x

[k]

= [ ˆ x

1

x ˆ

2

· · · x ˆ

k

]

T

is used. Let W be

W =





1 0 0 0

w

2

1 0 0 0 w

3

1 0 0 0 w

n

1





(35)

whose components w

i

for 2 i n are constant. Let the observer gain Y be of the form Y (x

1

) = −W

−1

·

w

1

(x

1

) 0

¸

(36) where w

1

is a smooth scalar function of x

1

. The parameters w

1

, · · · , w

n

have yet to be determined.

Candidates of state-dependent scaling matrices are parameterized as follows:

Φ = {Φ = φ

1

(x

1

)I

p+q

: φ

1

(x

1

) > 0, ∀x

1

∈ R} (37) Θ = {Θ = φ

1

(x

1

)I

p+q

: φ

1

(x

1

) > 0, ∀x

1

∈ R} (38) The scalar function φ

1

has yet to be determined. Choose a matrix Ψ from Ψ so that (25) holds and Ψ depends only on x

1

. Such a SD scaling matrix Ψ exists due to the definition of C. A simple choice is Ψ = ¯ Ψ. Extract M

[k]

from M as

M

[k]

=

·

Q

Tk

0 0 I

n

¸

M

·

Q

k

0 0 I

n

¸

=

·

M

[k]11

M

[k]12

M

[k]21

H

¸

H = W

−T

A ¯ W ˆ P ˜ + ˜ P W ˆ

T

A ¯

T

W

−1

Q

k

=



I

k

0 0

0 0 0

0 I

p+q

0 0 0 I

p+q



, Q

n

= I

n+2(p+q)

where I

k

is a k × k identity matrix. This matrix M

[k]

has the following properties.

(a-i) M

[k]

is independent of x

k

, x ˆ

k+1

, · · · , x ˆ

n

}.

(a-ii) M

[k]

does not include {s

k+1

, · · ·, s

n−1

, s

n

}.

(a-iii) M

[k]

< 0 implies M

[k−1]

< 0.

(a-iv) M

[n]

= M

(a-v) M

[1]

is jointly affine in {s

1

, φ

1

}. For k 2, M

[k]

is affine in s

k

. (a-vi) M

[k]

< 0 implies H < 0.

For achieving M < 0, the properties suggests a recursive procedure in which

M

[k]

(x

1

, x ˆ

[k−1]

) < 0, ∀(x

1

, x ˆ

[k−1]

) ∈ R × R

k−1

(39)

is solved for s

k

, and φ

1

(when k = 1), recursively from k = 1 through k = n. The property (a-vi)

claims that H < 0 should be secured beforehand. In order to obtain H < 0, employing an idea

(10)

which is analogus to recursive design of robust observers [6], the parameters {w

1

(x

1

), w

2

, . . ., w

n

} are selected by solving

H

hki

(x

1

) < −Γ

k

(x

1

)

−1

, ∀x

1

∈ R (40) for w

k

recursively from k = n down to k = 1. The notation H

hki

is defined as

H

hki

=

·

H

k,k

H

k,∗

H

∗,k

H

hk+1i

¸

, H

hni

= H

n,n

The next section shows how to determine appropriate matrices Γ

k

(x

1

) > 0 of C

0

functions which guarantee the existence of solutions to (39) for all k = 1, 2, . . . , n. The matrix H

hki

also satisfies the following.

(b-i) H

hki

does not include {w

k−1

, · · ·, w

2

, w

1

}.

(b-ii) H

hki

< −Γ

−1k

implies H

hk+1i

< −[Γ

−1k

]

hk+1i

(b-iii) H

h1i

= H

(b-iv) H

hki

is affine in w

k

.

The properties (a-v) and (b-iv) are advantageous to numerical computation of (39) and (40).

6 Existence of solution

Let [H

−1

]

11

denote the (1, 1)-component of the matrix H

−1

. The following can be obtained by modifying a result in [6] properly.

Lemma 2 Suppose that H(x

1

) < 0 is satisfied for all x

1

∈ R.

(i) Case k = 1 : There exist smooth functions {s

1

(x

1

), φ

1

(x

1

)} such that (39) is satisfied if

−[H

−1

]

11

λ

max³

−B

T

W

T

P H ˜

−1

P W B ˜

´

λ

max³

ΨC

11

C

11T

Ψ

´

< 1

4 (41)

holds for all x

1

∈ R.

(ii) Case k 2 : Assume that M

[k−1]

< 0 holds for all (x

1

, x ˆ

[k−2]

) ∈ R × R

k−2

. Then, there exists a smooth function s

k

(x

1

, x ˆ

[k−1]

) such that (39) is satisfied.

Here, λ

max

(·) denotes the maximum eigenvalue of a matrix. The inequality (39) is solvable recursively from k = 1 through k = n if H(x

1

) < 0 and (41) are satisfied for all x

1

. Let

Γ

k

=

·

γ

k

0 0 ν

k

Γ

k+1

¸

, k = 1, 2, . . . , n 1, Γ

n

= γ

n

where γ

i

, i = 1, 2, . . . , n and ν

i

, i = 1, 2, . . . , n −1 are real scalars. In order to achieve (40) and (41)

simultaneously, the following lemma is useful, which successfully extends the previous result of

the robust observer design [6] to much more general systems defined with the relaxed assumptions

(30-31).

(11)

Lemma 3 Let

1

(x

1

), γ

2

(x

1

), . . ., γ

n

(x

1

)} be any C

0

functions satisfying

0 < γ

i

(x

1

), ∀x

1

∈ R, i = 1, 2, . . . , n (42) 0 < inf

x1∈R

i

(x

1

)|a

i−1,i

(x

1

)|} , i = 2, 3, . . . , n (43) sup

x1∈R

i

(x

1

)|a

i−1,i

(x

1

)|} < +∞, i = 3, 4, . . . , n (44) Let

1

, ν

2

, . . ., ν

n−1

} be any constants satisfying

ν

1

1, ν

i

> 1, i = 2, 3, . . . , n 1 (45) Then, there exist a smooth function w

1

(x

1

) and constants w

2

, w

3

, . . ., w

n

which solve (40) se- quentially in descending order of k, where the existence of w

k

is independent of

1

, γ

2

, . . ., γ

k−1

} and {w

1

, w

2

, . . ., w

k−1

}.

Pick a C

0

function γ

1

(x

1

) such that

γ

1

λ

max³

B

T

W

T

P ˜ Γ

1

P W B ˜

´

λ

max³

ΨC

11

C

11T

Ψ

´

< 1

4 (46)

If −H

−1

< Γ

1

is achieved, this inequality implies (41). Thus, if we select γ

1

, . . ., γ

n

as (46), (42-44), Lemma 2 proves that s

k

and φ

1

solving M < 0 can be constructed recursively from k = 1 up to k = n.

According to the proof of Lemma 2, any C

0

function satisfying

¯

e

< φ

1

(x

1

) < e ¯

+

, x

1

∈ R (47) is a solution of M

[1]

< 0. The real numbers ¯ e

and ¯ e

+

are given by

¯

e

±

= 1 + ¯ c ¯ b ±

q

(1 + ¯ c ¯ b)

2

c

c (48)

where

¯

a = λ

max³

−B

T

W

T

P H ˜

−1

P W B ˜

´

¯ b = λ

max

(Z

bT

Z

b

), Z

b

= −B

T

W

T

P H ˜

−1

C

T

Ψ

¯

c = λ

max

(−ΨC

11

[H

−1

]

11

C

11T

Ψ)

Set φ

1s

= φ

1

. Then, according to Theorem 1, we achieve the condition (23) for global asymptotic stability when Σ

P

has neither dynamic uncertain components nor exogenous disturbances. How- ever, Theorem 1 requires a constant φ

1

when either dynamic uncertain components or exogenous disturbances is involved. Lemma 2 does not guarantee that the set (¯ e

, e ¯

+

) admits a constant solution φ

1

globally in x

1

. The following new result is the key to the existence of constant φ

1

. Lemma 4 Let

1

, ν

2

, . . ., ν

n−1

} be any positive real constants. Suppose that

1

, γ

2

, . . ., γ

n

} satisfy (42) and

sup

x1∈R

i

(x

1

)|a

i−1,i

(x

1

)|} < +∞, i = 2, 3, . . . , n (49)

(12)

(a) If B(x) satisfies

sup

x1∈R

λ

max

(B

i1

( x

1

)B

i1T

( x

1

))

|a

i,i+1

( x

1

)| < +∞, i = 1, 2, ..., n−1 (50) sup

x1∈R

λ

max

(B

n1

( x

1

)B

n1T

( x

1

))

|a

n−1,n

( x

1

)| < +∞ (51)

then, there exists a C

0

function γ

1

(x

1

) such that (42) and

λ

max³

B

T

W

T

˜

1

P W B ˜

´

α, ∀x

1

∈ R (52) γ

1

λ

max³

ΨC

11

C

11T

Ψ

´

< 1

, ∀x

1

∈ R (53)

hold for a finite constant α > 0.

(b) If γ

1

(x

1

) satisfies (42), (52) and (53) and H < −Γ

−11

holds, there exists a positive constant φ

1

such that (47) holds.

According to the proof, a constant solution φ

1

fulfilling (47) for ‘all’ x

1

∈ R is any real number belonging to (ˆ e

, ˆ e

+

), where constants ˆ e

and ˆ e

+

are given by

ˆ

e

±

= 1 + ˆ c ˆ b ±

q

(1 + ˆ c ˆ b)

2

c

c (54)

ˆ

a = sup

x1∈R

¯ a(x

1

), ˆ b = sup

x1∈R

¯ b(x

1

), ˆ c = sup

x1∈R

¯ c(x

1

)

Here, ˆ a, ˆ b and ˆ c are guaranteed to exist. Thanks to Lemma 4, the recursive design of output feedback controllers which accommodate nonlinear dynamic uncertainties proceeds as follows.

1) solve (40) for w

k

recursively in descending order of k with γ

n

, . . ., γ

2

, γ

1

given in (49) (43), (42), (52), (53) and (45).

2) solve (39) for s

k

and a constant φ

1

(when k = 1) recursively in ascending order of k.

Using Schur complements formula, design equations (39) and (40) in each step k reduce to scalar inequalities which are affine in the decision variables. We thereby arrive at the main result by setting φ

1s

= φ

d

= φ

1

and ˇ φ

1d

= 1.

Theorem 2 Suppose that the system Σ

P

is in the generalized robust output-feedback form and satisfies (50-51). Then, the system Σ

P

can be globally uniformly asymptotically stabilized, and the L

2

-gain from r to e can be rendered less than or equal to τ for all admissible uncertainties Σ

by the output-feedback law (5-6).

It is emphasized that the disturbance attenuation level τ can be made arbitrarily small. The scaling factor φ

1

can be chosen as either a constant or a function of x

1

. The uncertainty Σ

is not allowed to be dynamic unless φ

1

is a constant. Non-constant φ

1

may often lead us to a less complicated feedback gain whose growth order and local gain are not very large.

It is worth noting that Theorem 2 does not restrict Σ

to ISS dynamic systems. When we restrict

Σ

to ISS dynamics considered in [7], the outcome of Theorem 2 is similar to Corollary 1 of [7].

Figure 1: State transition of output feedback control.

参照

関連したドキュメント

It can be shown that cubic graphs with arbitrarily large girth exist (see Theorem 3.2) and so there is a well-defined integer µ 0 (g), the smallest number of vertices for which a

The purpose of this paper is to guarantee a complete structure theorem of bered Calabi- Yau threefolds of type II 0 to nish the classication of these two peculiar classes.. In

The main task of this paper is to relax regularity assumptions on a shape of elastic curved rods in a general asymptotic dynamic model and to derive this asymptotic model from a

Keywords: continuous time random walk, Brownian motion, collision time, skew Young tableaux, tandem queue.. AMS 2000 Subject Classification: Primary:

This paper is devoted to the investigation of the global asymptotic stability properties of switched systems subject to internal constant point delays, while the matrices defining

36 investigated the problem of delay-dependent robust stability and H∞ filtering design for a class of uncertain continuous-time nonlinear systems with time-varying state

The purpose of this paper is to apply a new method, based on the envelope theory of the family of planes, to derive necessary and sufficient conditions for the partial

This paper presents an investigation into the mechanics of this specific problem and develops an analytical approach that accounts for the effects of geometrical and material data on