• No results found

S-581 83 Linkoping, Sweden tel: +46 13 281622 fax: +46 13 282622 email:

N/A
N/A
Protected

Academic year: 2021

Share "S-581 83 Linkoping, Sweden tel: +46 13 281622 fax: +46 13 282622 email:"

Copied!
13
0
0

Loading.... (view fulltext now)

Full text

(1)

Mixed Uncertainties

A. Helmersson

Department of Electrical Engineering Linkoping University

S-581 83 Linkoping, Sweden tel: +46 13 281622 fax: +46 13 282622 email:

andersh@isy.liu.se

September 26, 1994 1995 ECC

Abstract

This paper presents the solution to the mixed synthesis problem, and how to design gain scheduling controllers with linear fractional transfor- mations (LFTs). The system is assumed to have a parameter dependences described by LFTs The goal of the controller is, with real-time knowledge of the parameters, to provide disturbance and error attenuation. The paper treats both complex and real parametric dependencies.

Keywords:

Gain sceduling, linear fractional transforms, structured singular values, parametric dependent systems, synthesis.

1 Introduction

During the last couples of year synthesis methods for gain scheduling using linear fractional transformations (LFTs) and structured singular values have been developed, see e.g 10, 9, 7]. The idea behind this approach is to let the controller have access to some of the uncertainties or perturbations of the system to be controlled. The design leads to two linear matrix inequalities (LMIs), which are coupled by a third LMI and non-convex rank conditions.

This is analogous with the two-Riccati algorithm for solving

H1

synthesis.

In the referred papers 10, 9, 7], the shared uncertainties have been assumed

to be complex, possibly repeated scalar blocks. In this paper the mixed



design

problem (that is with both real and complex uncertainties) is elaborated. The

synthesis algorithm has the same structure as in the pure complex case with two

LMIs coupled by a third LMI and rank conditions. The di erence to the pure

(2)

M



-

Figure 1: System with uncertainties.

complex case is that blocks corresponding to real uncertainties in the connecting LMI disappear only the rank conditions remain.

The structure of the problem is treated in a unied and general way cover- ing both discrete and continuous systems, which are equivalent using bilinear transformation between the

z

-domain and the

s

-domain 2].

1.1 Notations

X



denotes the complex conjugate transpose of

X



X >

( ) 0 a hermitian (

X

=

X



) positive denite (semidenite) matrix

X;

= (

X

)

;1



Xy

is the Moore- Penrose pseudo inverse of

X

 diag

X1X2

] a block-diagonal matrix composed of

X1

and

X2

 rank

X

denotes the rank of the matrix

X

 herm

X

=

12

(

X

+

X

)

S

(

::

) denotes the Redhe er star product 



(

X

) the maximal singular value of

X

.

2 -analysis and LFTs

This section gives a short review on structured singular values and linear frac- tional transformations (LFTs), see also e.g. 2].

2.1 Denitions

We will here use a similar notations and denitions as is used in 12, 13]. The denition of



depends upon the underlying block structure

N

of the uncer- tainties , which could be either real or complex, see Figure 1. For notational convenience we assume that all uncertainty blocks are square. This can be done without loss of generality by adding dummy inputs or outputs.

Given a matrix

M 2C

n



n and three non-negative integers

f

r ,

f

c and

f

C , with

f

=

f

r +

f

c +

f

C

n

the block structure is an

f

-tuple of positive integers

N

= 

n1:::n

f

rn

f

r+1:::n

f

r+

f

cn

f

r+

f

c+1::: n

f

r+

f

c+

f

C

] (1)

where

P

f i

=1n

i =

n

for dimensional compatibility. The set of allowable per-

turbations is dened by a set of block diagonal matrices

X 2 C

n



n dened

(3)

by

X

=

f

 = diag



r

1I

n

1:::

rn

frI

n

fr



c

1

I

n

fr+1:::

cf

cI

n

fr+fc

 C

1:::

 Cn

C

] :



ri

2R

ci

2C

 Ci

2C

n

fr+fc+i

n

fr+fc+ig:

(2) Assuming the uncertainty structure

N

, the structured singular value



of a matrix

M2C

n



n is dened by



=



min

2X

f





() : det(

I;



M

) = 0

g



;1

(3)

and if no 

2X

satises det(

I;



M

) = 0 then



(

M

) = 0.

2.2 Upper and Lower Bounds

Generally the structured singular value cannot be exactly computed, and instead we have to resort to upper and lower bounds, which are usually sucient for most practical applications. A tutorial review of the complex structured singular value is given in 8].

An upper bound can be determined using convex methods, either involving minimization of singular values with respect to a scaling matrix or by solving a linear matrix inequality (LMI) problem. The upper bound is conservative in the general case, but can be improved by branch and bound schemes.

A lower bound can be found by maximizing the real eigenvalue of a scaled matrix. This bound is nonconservative in the sense that if the true global maximum is found it is equal to



. However, since the problem is not convex, we cannot guarantee that we nd the global maximum.

We will here focus on the computation of the upper bound, which we here denote

 

, in order to distinguish it from the true



function. The upper bound



can be computed as a convex optimization problem. For complex uncertainties it is determined by



(

M

) = inf D

2D





(

DMD;1

) (4)

where

D

is the set of block diagonal Hermitian matrices that commute with

X

, that is

D

=

f

0

<D

=

D2C

n



n :

D

 = 

D8



2Xg:

(5) This problem is equivalent to an LMI problem



(

M

) = inf >

P

2D0

f

:

MPM<2Pg:

(6) Real uncertainties can be included in the LMI problem for computing the upper bound (see e.g. 3, 12, 13]):



(

M

) = inf >

P G

2D2G0

f

:

MPM

+

j

(

GM;MG

)

<2Pg

(7)

(4)

where

G

=

fG

:

G

=

G2C

n



n

G

 = 

G8



2Xg:

(8) Every

G2G

is block diagonal with zero blocks for complex uncertainties. If we let

G

=

f

0

g

in (7) we recover the complex upper bound (6).

We can reformulated (7) as a positive real property



(

M

) = inf >

W

2W0





: herm

;

(

I

+

1



M

)

W

(

I;1



M

)

>

0



(9) where herm

X

=

12

(

X

+

X

) and

W

=

fW

=

P

+

jG

:

P 2DG2Gg

. Note that herm

W >

0 always and that

W

=

W

for complex uncertainties. Another equivalent reformulation of (7) is



(

M

) = inf >

P G

2D2G0

f

: 

;1



DMD;1;jG

(

I

+

G2

)

;12 <

1

g:

(10)

2.3 Linear Fractional Transformations (LFTs)

Suppose

M

is a complex matrix partitioned as

M

=

M

11 M

12

M

21 M

22

2C

(

p

1+

p

2)(

m

1+

m

2)

(11) and let  u

2C

m

1

p

1

and  l

2C

m

2

p

2

. The upper and lower linear fractional transformations (LFTs) are dened by

F

u (

M

 u ) =

M22

+

M21

 u (

I;M11

 u )

;1M12

(12) and

F

l (

M

 l ) =

M11

+

M12

 l (

I;M22

 l )

;1M21

(13) respectively. Clearly, the existence of the LFTs depends on the invertibility of

I;M

11

 u and

I;M22

 l respectively.

The Redhe er star product 11] is a generalization of the LFTs. Assume that

Q

is partitioned similarly to

M

. Then the star product is dened by

S

(

QM

) =

F

l (

QM11

)

Q12

(

I;M11Q22

)

;1M12

M

21

(

I;Q22M11

)

;1Q21 F

u (

MQ22

)

:

(14) A block diagram illustrating the star product is given in Figure 2.

Note that the denition above is dependent on the partitioning of the ma- trices

Q

and

M

. The LFTs can be dened by the

S

notation, as

F

u (

M

 u ) =

S

( u

M

) and

F

l (

M

 l ) =

S

(

M

 l )

:

The star product is associative, that is

S

(

AS

(

BC

)) =

S

(

S

(

AB

)

C

)

:

(5)

M

!

!

!

!

!

a a a a a

,

S

(

QM

)

Figure 2: The Redhe er star product.

3 Synthesis

In this section we will treat the



synthesis problem when uncertainties in the controller may be shared by the controller.

3.1 Structure of shared uncertainties

Formally we treat this by letting the controller

K

have access to a copy of the system's uncertainties. Not all blocks in the uncertainty matrix  are accessible and we denote the accessible uncertainties by

L

R

() = diag

1I

r

12I

r

2:::

f

1I

r

f1

] (15) where

f1

=

f

r +

f

c . We only allow blocks that are scalar repeated blocks of the type



i

I

r

i

where



could be either real or complex. The structure of the controller accessible uncertainties are given by

R

= 

r1:::r

f

rr

f

r+1:::r

f

r+

f

c

0

:::

0] (16) Note that all blocks corresponding to full complex blocks are not allowed and, consequently,

r

i = 0 for those.

3.2 The General Synthesis Problem

The general synthesis problem can be depicted as in Figure 3. Find a controller

K

, possibly with a given order, to the system

M

(or the augmented system ~

M

), such that



(

S

(

MK

)) =



(

S

( ~

MK

)) is minimized or below a specied value.

The diagrams below show three equivalent representations of the



-synthesis problem. In the left diagram

M

and

K

have separated uncertainty blocks,  and

LR

() respectively. In the middle diagram the uncertainties are explicitly shared. Here ~ = diag

LR

()]. In the right diagram through connections between

K

and  have been included in ~

M

.

Here we use a somewhat more general structure of ~ than used so far. With this formulation it is possible to catch most of the relevant synthesis problems in the



setting. The generalized structure of ~ is nothing more than a reordering of rows and columns, so that both

M

and

K

may share the same uncertainties.

For instance, if  includes frequency, either

s

or

z

, we can handle continuous-

time and discrete-time problems respectively.

(6)

L

R

()

K

- -

M



-

,

K M

!

!

!

!

!

a a a a a

~

- -

,

K -

~

M

~

-

Figure 3: Equivalent representations of a system

M

controlled by

K

that have access to a reduced copy

LR

() of the system's uncertainties . The augmented uncertainty block ~ = diag

LR

()].

3.3 An Ane Problem

We assume that

M

is partitioned as:

M

=

M

11 M

12

M

21 M

22

2C

(

n

+

m

)(

n

+

p

):

(17) Without restriction we may assume that

M22

= 0 (

M

is strictly proper). If

M

22

6

= 0, we replace

M

with



M

=

S

(

MN

M ) =

M

11 M

12

M

21

0

(18) with

N

M =

0

I

I ;M

22

(19) Thus, if we can nd a controller 

K

to the modied system 

M

, then

K

=

S

(

N

M

K

 ) will solve the original problem. Since the star product inverse of

N

M

exists for all

M

, the modied problem is equivalent to the original one as long as

I;M22K

 is not singular.

Due to the simpler structure of

M

, the star product can be rewritten as the following matrix expression, which is ane in 

K

:

S

( 

MK

 ) =

M11

+

M12KM



21

=

Q

+

UKV:

 (20)

3.4 Mathematical Preliminaries

In this section we give two important and related lemmas to be used for solv- ing



synthesis problems. The lemmas, theorems and proofs are made on the assumption that the matrices are complex. The real case can be obtained by replacing

C

with

R

without a ecting the validity of the theory.

We rst consider an ane optimization problem 1].

(7)

Lemma 3.1 Assume that

Q2C

n



n and that both

U 2C

n



p and

V 2C

m



n have full column rank and row rank respectively. Suppose

U? 2C

n

(

n

;

p

)

and

V

? 2 C

(

n

;

m

)

n are chosen such that

U U?

,

V

V

?

are both invertible, and that

UU?

= 0 and

VV?

= 0. Then

K

2C

inf

p m





(

Q

+

UKV

)

<

(21) if and only if

V

?

;

Q



Q; 2

I



V



?

<

0 (22)

and

U



?

;

QQ



; 2

I



U

?

<

0

:

(23)

Proof: See 1].

2

A related lemma 4, 5], goes as follows.

Lemma 3.2 Given a hermitian matrix

Q 2 C

n



n and

U

,

V

,

U?

and

V?

as above. Then

Q

+

UKV

+ (

UKV

)

 <

0 (24) is solvable for

K2C

p



m if and only if



U



? QU

?

<

0

V

? QV



?

<

0

:

(25)

Proof: Necessity of (25) is clear: for instance,

U?U

= 0 implies

U?QU?<

0 when pre- and post-multiplying (24) by

U?

and

U?

respectively. For details on

the suciency part, see 4, 5].

2

In 5] the set of

K

solving (24) is parametrized.

3.5 Complex -synthesis

If we restrict the problem to only have complex uncertainties (real uncertainties can be conservatively treated us complex ones), the following lemma from 10]

applies.

Lemma 3.3 Let

Q

,

U

,

V

,

U?

and

V?

be given as above. Then

K

2C

D inf

2Dp m 



;D

(

Q

+

UKV

)

D;1<

(26) if and only if there exists a

P 2D

such that

V

?

;

Q



PQ; 2

P



V



?

<

0 (27)

and

U



?

;

QP

;1

Q



; 2

P

;1



U

?

<

0

:

(28)

Proof: The necessity of (27) follows by rewriting the inequality in (26) as (

Q

+

UKV

)

P

(

Q

+

UKV

)

;2P<

0

and pre and post-multiplying by

V?

and

V?

respectively. The necessity of (28)

is obtained analogously. For the suciency part refer to 10].

2

(8)

3.6 Mixed -synthesis

If both real and complex uncertainties are involved, we have a mixed



problem.

The following lemma then applies.

Lemma 3.4 Let

Q

,

U

,

V

,

U?

and

V?

be given as above. Then

K

2C

D G inf

2D2Gp m







;

1



D

(

Q

+

UKV

)

D;1;jG

(

I

+

G2

)

;12 <

1



(29) if and only if there exists an

W 2W

where

W

=

fW

=

P

+

jG

:

P 2DG2Gg

, such that

herm

;V?

(

I

+ 

1Q

)

W

(

I ;



1Q

)

V?>

0 (30) and

herm

;U?

(

I

+

1



Q

)

W;1

(

I;



1Q

)

U?>

0 (31) where herm(

X

) =

12

(

X

+

X

).

Proof: Let ^

Q

=

1



Q

, ~

Q

= (

DQD

^

;1 ;jG

)(

I

+

G2

)

;12

, ~

U

=

DU

, ~

V

=

VD

;1

(

I

+

G2

)

;12

and ^

K

= 

1K

. Then





Q

~ + ~

UK

^

V

~

<

1 (32) if and only if

~

U



?

( ~

QQ

~

;I

)~

U?

=

U?D;1

(

DQD

^

;1;jG

)(

I

+

G2

)

;1

(

D;Q

^

D

+

jG

)

;I D;U?

=

U?QP

^ U

Q

^



+

j

( ^

QG

U

;G

U

Q

^



)

;P

U

U

?

= herm

U?

( ^

Q

+

I

)(

P

U

;jG

U )( ^

Q;I

)

U? <

0

(33) with

P

U =

D;1

(

I

+

G2

)

;1D;

and

G

U =

D;1

(

I

+

G2

)

;1GD;

. Analogously we obtain

~

V

?

( ~

QQ

~

;I

)~

V?

=

V?D;1

(

D;Q

^

D

+

jG

)(

DQD

^

;1;jG

)

;

(

I

+

G2

)

D;V?

=

V?Q

^

P

V

Q

^ +

j

(

G

V

Q

^

;Q

^

G

V )

;P

V

V



?

= herm

V?

( ^

Q

+

I

)(

P

V +

jG

V )( ^

Q;I

)

V? <

0

(34) with

P

V =

DD

and

G

U =

DGD

. Next, by introducing

X

=

P

V +

jG

V

2W

and observing that

W;1

=

P

U

;jG

U , we can conclude the proof.

2

Remark 3.1 Lemma 3.4 follows by letting

W

=

W

=

P

.

(9)

3.7 Shared uncertainties

We will now consider the problem with shared uncertainties. We augment the system ~

M

with through connections between the uncertainty block ~ and

K

. We assume that

M22

= 0.

~

M

=

2

4

~

M

11 M

~

12

~

M

21

0

3

5

=

2

6

6

6

4 M

11

0

M12

0

0 0 0

I

r

M

21

0 0 0

0

I

r 0 0

3

7

7

7

5

=

2

6

6

6

4

Q

0

U

0 0 0 0

I

r

V

0 0 0

0

I

r 0 0

3

7

7

7

5

(35)

:

Thus, ~

U

= diag

M12I

r ] and ~

V

= diag

M21I

r ]. Then, ~

U?

= diag

U?

0] and

~

V

?

= diag

V?

0]. Applying Lemma 3.4 we obtain

herm

;V?

(

I

+

1



Q

)

S

(

I;1



Q

)

V?>

0 (36) and

herm

;U?

(

I

+

1



Q

)

R

(

I;



1Q

)

U?>

0 (37) where

S

,

R2C

n



n are the upper left block of

W

and

W;1

respectively.

W

=

S 

 

(38) and

W

;1

=

R 

 

:

(39)

The matrices

R

and

S

are mutually constrained by these two equations. Next, we will show that this constraint can be reformulated as a rank condition on

I;R S

and, in the case of complex uncertainties, an additional LMI.

Before going much further we need to specify the structure of ~,

W

and

W;1

more explicitly. The uncertainty matrix ~ is assumed to be a block diagonal matrix such that

 = diag

LR

()]

= diag

1I

n

1:::

f

1I

n

f1

 f

1+1:::

 f

1I

r

1:::

f

1I

r

f1

] (40) where

f1

=

f

r +

f

c is the number of scalar blocks and

f

is the total number of blocks. Shared blocks are restricted to the type



i

I

where



i is either a real or a complex scalar.

We assume that

W

have the following block structure.

W

=

2

6

6

6

6

6

6

6

6

6

6

6

6

6

6

6

4 W

111

0

0

W121

0

0 0

W112

0 0

W122

0 ... ... ... ... ... ... ... ...

0 0

W11

f 0 0

W12

f

W

211

0

0

W221

0

0 0

W212

0 0

W222

0 ... ... ... ... ... ... ... ...

0 0 f 0 0 f

3

7

7

7

7

7

7

7

7

7

7

7

7

7

7

7

5

(41)

(10)

For notational convenience, all

W12

i ,

W21

i and

W22

i are included even if they are not shared and are in such case empty matrices. The structure of

W;1

will be identical to the structure of

W

.

We now denote

W

i =

W

11

i

W12

i

W

21

i

W22

i

:

Specically, for uncertainties that are not shared we have

W

i =

W11

i .

The structure of

W

i is given by the type of uncertainty it refers to. We will now look specically at two kinds of uncertainties that may be shared: complex and real. They are both of the type  i =



i

I

, where



i is either complex or real. For complex uncertainties we have

W

i =

W

i

 >

0 for real uncertainties

W

i could be any square matrix such that herm

W

i

>

0.

We will now assume that

R

and

S

are diagonal block matrices with a struc- ture identical to

W11

. Depending on the structure of , and consequently of

W

, the matrices

R

and

S

are related to each other. The next two sections treats the two cases of shared uncertainties: complex and real. For notational convenience we drop the sux

i

and consider each block at a time.

3.7.1 Complex uncertainties

If we are dealing with complex uncertainties we restrict

W

to be hermitian and positive denite

W

=

W

=

P >

0. Then we have the following lemma from

7]. Lemma 3.5 Suppose that

R

=

R 2C

n



n ,

S

=

S 2C

n



n , with

R >

0 and

S >

0. Let

r

be a positive integer. Then there exists matrices

N 2C

n



r and

L2C

r



r such that

L

=

L

and

P

=

S N

N



L

>

0 and

P;1

=

R 

 



if and only if

R I

n

I

n

S

0 and rank

;S;R;1r:

Proof: (

)

) Using Schur complements,

R

=

S;1

+

S;1N

(

L;NS;1N

)

;1NS;1

. Inverting, using the matrix inversion lemma, gives that

R;1

=

S;NL;1N

. Hence

S;R;1

=

NL;1N

0, and indeed, rank(

S;R;1

) = rank(

NL;1N

)



r

. (

(

) By assumption, there is a matrix

N 2C

n



r such that

S;R;1

=

NN

.

Dening

L

=

I

r completes the construction.

2

3.7.2 Real uncertainties

For real uncertainties we have a more general set of matrices such that herm

W >

0. Note that herm

W >

0 implies that

W;1

exists and that herm(

W;1

)

>

0.

The following lemma then applies.

Lemma 3.6 Suppose that

R 2 C

n



n ,

S 2 C

n



n , with herm

R >

0 and herm

S>

0. Let

r

be a positive integer. Then there exists matrices

N 2C

n



r ,

M 2C

r



n and

L2C

r



r such that

L

is nonsingular and

W

=

S N

M L



herm

W >

0 and

W;1

=

R 

 



(11)

if and only if

rank

;S;R;1r:

Proof: (

)

) Using Schur complements,

R

=

S;1

+

S;1N

(

L;MS;1N

)

;1MS;1

. Inverting, using the matrix inversion lemma, gives that

R;1

=

S;NL;1M

. Hence,

S;R;1

=

NL;1M

, and indeed, rank(

S;R;1

) = rank(

NLM

)

r

. (

(

) Dene

t

= rank(

S;R;1

). Find any

M2C

t



n of full row rank, such that

S;R

;1

=

KM

t , with

K2C

n



t having full column rank. Choose

N

t as a linear function of

L

t by letting

N

t = (

S;R;1

)

M

t

yL

t where

M

t

y

=

M

t



(

M

t

M

t



)

;1 2

C

n



t is the Moore-Penrose inverse of

M

t . By dening

M

=

M

t

0

,

N

=

N

t 0

and

L

= diag

L

t

I

r

;

t ], we obtain

W

=

S N

M L

=

2

4

S

0 0

M

t 0 0 0 0

I

r

;

t

3

5

+

2

4

(

S;R;1

)

M

t

y

I

t

0

3

5

L

t

0

I

t 0

:

Then, herm

W

i

>

0 is equivalent to

S

+

S M

t



M

t 0

+

(

S;R;1

)

M

t

y

I

t

L

t

0

I

t

+ ]

>

0



(42) where ]



denotes the complex conjugate transpose of the previous term. Using lemma 3.2 with

U?

=

I

n (

R;1;S

)

M

t

y

and

V?

=

I

n 0

, this is in turn equivalent to

U



?

S

+

S M

t



M

t 0

U

?

=

S

+

S

+ (

M

t

yM

t )



(

R;1;S

)



+ (

R;1;S

)

M

t

yM

t

=

R;1

+

KM

t +

R;

+ (

KM

t )

;

(

KM

t

M

t

yM

t )

;KM

t

M

t

yM

t

=

R;1

+

R;

= 2herm(

R;1

)

>

0



and

V

?

S

+

S M

t



M

t 0

V



?

=

S

+

S

= 2herm

S>

0

:

Thus (42) is satised implying that

L

is nonsingular, which concludes the proof.

2

3.8 Summary

It has been shown that the matrix inequality

S R;1

or equivalently

S I

I R

0 (43)

disappears in the case of real uncertainties. When used for gain scheduling the parameters are most naturally described by real uncertainites. Real uncertain- ties are parameters while complex uncertainties can be interpreted as single- input-single-output dynamics linear systems (rational transfer functions) with time-varying state-space matrices. Assuming real uncertainties LMI condition (43) used in for instance 9, 7] can be dispensed with.

Thus we have arrived at the following mixed- synthesis theorem

References

Related documents

Mitt uppdrag som chef för Program Management var att säkerställa en stabil/växande volym och lönsam affär med våra befintliga kunder med god kunddialog och en nöjd kund.. I

Multiconfigurational quantum chemistry methods, and especially the multiconfigurational self- consistent field (MCSCF) and multireference perturbation theory (MRPT2), are

To explore the usefulness of symbolic and algebraic methods, we use polynomials over finite fields (see section 2) applied to DEDS with industrial sized complexity: The landing

However, employing the inherent structure of some important problems, such as H 1 and gain-scheduling synthesis, convexity can be recovered and the existence of a controller K can

Dynamic uncertainties are treated by using constant scaling matrices in the nonlinear case and dynamic scaling systems in the linear (complex  ) case.. The nonlinear problem

Time-varying uncertainties with bounds on rate of variation can also be included in the same framework by left and right non-square multipliers that correspond to an

We have shown that it is possible to transfer the state between two such points in nite time if they can be joined by a continuous curve of equilibria, provided the linearization

This along with the ensured monotone behavior suggest that the fuzzy approach might be a good alternative to neural nets when applied to predictive