multivariable control systems ecse 6460 fall 2009 lecture 26: 4 december 2009

39
Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Upload: matthew-lamb

Post on 18-Dec-2015

229 views

Category:

Documents


1 download

TRANSCRIPT

Page 1: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Multivariable Control SystemsECSE 6460

Fall 2009

Lecture 26: 4 December 2009

Page 2: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Uncertain MIMO Plant

Gp(s)

Page 3: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Frequency Response

-80

-60

-40

-20

0

From: In(1)

To:

Out

(1)

10-2

100

102

-80

-60

-40

-20

0

To:

Out

(2)

From: In(2)

10-2

100

102

Bode Diagram

Frequency (rad/sec)

Mag

nitu

de (

dB)

Page 4: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Feedback Control Loop

Gp(s)K(s)+

-

• practically zero steady state error for step reference• rise time less than 1 second• overshoot less than 5%• robustly for any value of c

Page 5: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Menu for the day

Representing the uncertaintya la unstructured uncertaintya la structured uncertainty

Assessment of robust stabilitya la unstructured uncertaintya la structured uncertainty

Assessment of robust stabilitya la unstructured uncertaintya la structured uncertainty

Page 6: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009
Page 7: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Unstructured Uncertainty

G (s)

Wo (s) (s)

G (s)

Wi (s) (s)

++

++

Output uncertainty

Input uncertainty

Page 8: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Output Uncertainty

G (s)

Wo (s) (s)

++

Output uncertainty

We need to select a weight such that:

Page 9: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Input Uncertainty

We need to select a weight such that:

G (s)

Wi (s) (s)

++

Input uncertainty

Page 10: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Weight Wo(s)

10-5

100

105

10-3

10-2

10-1

100

101

Output Uncertainty

Refer to uncert1.m

Page 11: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Structured Uncertainty

Page 12: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Structured Uncertainty

G (s)

Ws (s)

++

relax to

Page 13: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Relaxation from real to complex uncertainty

Page 14: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009
Page 15: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Assessment of Robustness

In this topic, we first design a controller for the nominal system, then we assess the robustness properties of the system.

Two control designs are explored, decentralized inverse-based, and H-infinity synthesis.

Gp(s)K(s)+

-

Page 16: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Decentralized Inverse Based

Based on nominal plant model

Refer to decent.m

Page 17: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Nominal Performance

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 20

0.5

1

1.5

System: TnomOutput: Out(2)Time (sec): 0.451Amplitude: 0.9

To:

Out

(2)

Nominal Close Loop Plant Step Response

Time (sec)

Am

plitu

de 0

0.5

1

1.5 System: TnomOutput: Out(1)Time (sec): 0.556Amplitude: 0.901

To:

Out

(1)

Page 18: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

H synthesis for nominal plant

Use S/T synthesis with:

Refer to hinf_ex.m

Page 19: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Nominal Performance

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 20

0.2

0.4

0.6

0.8

1

System: TnomOutput: Out(2)Time (sec): 0.465Amplitude: 0.901

To:

Out

(2)

Nominal Close Loop Plant Step Response

Time (sec)

Am

plitu

de 0

0.2

0.4

0.6

0.8

1

System: TnomOutput: Out(1)Time (sec): 0.465Amplitude: 0.901

To:

Out

(1)

Page 20: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Stability Assessment

K(s)+

-

G (s)

Wo (s) (s)

++

Unstructured Uncertainty

Page 21: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Stability Assessment

K(s)

-

G (s)

Wo (s)

++

Page 22: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Stability Assessment

10-5

100

105

10-4

10-3

10-2

10-1

100

101

102

max

(T) (blue) vs 1/|Wo| (red)

Refer to decent.m

Page 23: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Stability Assessment

10-5

100

105

10-6

10-5

10-4

10-3

10-2

10-1

100

101

102

max

(T) (blue) vs 1/|Wo| (red)

Refer to hinf_ex.m

Page 24: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Stability Assessment

K(s)+

-

G (s)

Ws (s)

++

Structured Uncertainty

Page 25: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Stability Assessment

K(s)+

-

G (s)

Ws (s)

++

Structured Uncertainty

Page 26: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Structured Singular Value

10-5

100

105

-90

-80

-70

-60

-50

-40

-30

-20

-10M

agni

tude

(dB

)

(M(j))

(rad/sec)

Peak at 0.130

For decentralized controller

Page 27: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Perturbed system responses

-0.5

0

0.5

1

1.5T

o: O

ut(1

)

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 20

0.5

1

1.5

To:

Out

(2)

Perturbed Close Loop Plant Step Responses

Time (sec)

Am

plitu

de

Page 28: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

10-5

100

105

-110

-100

-90

-80

-70

-60

-50

-40

-30

-20

-10M

agni

tude

(dB

)

(M(j))

(rad/sec)

Structured Singular Value

Peak at 0.150

Page 29: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Perturbed system responses

-0.5

0

0.5

1

1.5

To:

Out

(1)

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 20

0.2

0.4

0.6

0.8

1

To:

Out

(2)

Perturbed Close Loop Plant Step Responses

Time (sec)

Am

plitu

de

Page 30: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009
Page 31: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Nominal Performance

10-5

100

105

10-6

10-5

10-4

10-3

10-2

10-1

100

101

max

(S) (blue) vs 1/|Wp| (red)

Refer to hinf_ex.m

Page 32: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Performance Assessment

K(s)+

-

G (s)

Wo (s)

++

Unstructured Uncertainty

+WP (s)

Page 33: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Performance Assessment

K(s)+

-

G (s)

Wo (s)

++

Unstructured Uncertainty

+WP (s)

Page 34: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Performance Assessment

N

P

Page 35: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Performance Assessment

10-5

100

105

-30

-25

-20

-15

-10

-5

0

5

10

Mag

nitu

de (

dB)

(N(j))

(rad/sec)

Page 36: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Performance Assessment

K(s)+

-

G (s)

Ws (s)

+

+

Structured Uncertainty

+

WP (s)

Page 37: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Performance Assessment

K(s)+

-

G (s)

Ws (s)

+

+

Structured Uncertainty

+

WP (s)

Page 38: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Performance Assessment

N

P

Page 39: Multivariable Control Systems ECSE 6460 Fall 2009 Lecture 26: 4 December 2009

Robust Performance Assessment

10-5

100

105

-60

-50

-40

-30

-20

-10

0

10

Mag

nitu

de (

dB)

(N(j)) Structured vs Unstructured

(rad/sec)