15.07.2014 Views

solutions - UBC Mechanical Engineering

solutions - UBC Mechanical Engineering

solutions - UBC Mechanical Engineering

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Foundations in Control <strong>Engineering</strong> (MECH550P-201)<br />

Final exam (50 points). April 17, 3:30pm–6pm<br />

Sketch of <strong>solutions</strong><br />

1. Realize the following transfer functions in both controllable and observable<br />

canonical forms. Select a minimal realization for each transfer function. Explain<br />

the reason why your selected realization is minimal.<br />

(10pt)<br />

[ ]<br />

s<br />

(a) G(s) =<br />

2 + 3s + 3<br />

s 2 + 3s + 2 1 s2 + 3s + 4<br />

. (3-inputs, 1-output)<br />

s 2 + 3s + 2<br />

⎡ [ ] [ ] ⎤<br />

[ ]<br />

03 I A B<br />

3<br />

03<br />

CCF:<br />

= ⎣ −2I<br />

C D [ [ 3 −3I 3 I 3<br />

⎦<br />

] [ ] ] [ ]<br />

1 0 2 0 0 0 1 1 1<br />

OCF:<br />

[ A B<br />

C D<br />

⎡<br />

]<br />

= ⎣<br />

[ ] [ [ ]<br />

0 −2<br />

[ 1 0 2<br />

]<br />

[<br />

1 −3<br />

] [ 0 0 0<br />

]<br />

0 1 1 1 1<br />

]<br />

⎤<br />

⎦<br />

Observable canonical form is minimal, since it is controllable.<br />

⎡<br />

2s + 3<br />

⎤<br />

⎢ s + 1 ⎥<br />

(b) G(s) = ⎣ ⎦ . (1-input, 2-outputs)<br />

s + 3<br />

s + 2<br />

⎡<br />

] ⎤<br />

[ ] A B<br />

CCF:<br />

=<br />

C D<br />

⎢<br />

⎣<br />

[ ] [ 0 1 0<br />

[ [<br />

−2<br />

] [<br />

−3<br />

] ] [<br />

1<br />

]<br />

2 1 2<br />

1 1 1<br />

⎥<br />

⎦<br />

OCF:<br />

[ A B<br />

C D<br />

⎡<br />

]<br />

=<br />

⎢<br />

⎣<br />

[ ]<br />

02 −2I 2<br />

I 2 −3I 2<br />

[<br />

02 I 2<br />

]<br />

⎡<br />

⎢<br />

⎣<br />

[ 2<br />

1<br />

[ 1<br />

1<br />

[ 2<br />

1<br />

]<br />

]<br />

]<br />

⎤<br />

⎥<br />

⎦<br />

⎤<br />

⎥<br />

⎦<br />

Controllable canonical form is minimal, since it is observable.<br />

1


2. Consider a system: ⎧<br />

⎪⎨<br />

⎪ ⎩<br />

ẋ(t) =<br />

[ 0 2<br />

1 −1<br />

]<br />

x(t),<br />

y(t) = [ 1 1 ] x(t).<br />

For each of the following specified locations of the observer poles:<br />

• eig(A − LC) = −2, −3.<br />

• eig(A − LC) = −3, −3.<br />

• eig(A − LC) = −1 ± j.<br />

is it possible to design an observer? If the answer is “yes”, design an observer.<br />

For the system given above, what is the condition for given observer pole<br />

locations to be achievable?<br />

(10pt)<br />

By the direct method,<br />

det(λI − (A − LC)) = · · · = (λ + 2)(λ + l 1 + l 2 − 1).<br />

Therefore, for whatever choice of L, the eigenvalue of s = −2 is impossible<br />

to move. Thus, the first observer poles are possible to achieve with L such<br />

that l 1 + l 2 = 4, but other two observer poles are not possible to achieve.<br />

The condition for given observer pole locations to be achievable is that they<br />

include s = −2.<br />

3. Consider a system {<br />

ẋ(t) = αx(t) + u(t),<br />

y(t) = x(t).<br />

(a) For α = −1, design a state feedback with an integrator such that the<br />

closed-loop poles are s = −1, −1.<br />

(10pt)<br />

det(λI−(A aug −B aug<br />

[<br />

K<br />

Ka<br />

]<br />

) = · · · = λ 2 +(K−α)λ−K a = λ 2 +2λ+1<br />

When α = −1, K = 1 and K a = −1.<br />

(b) Suppose that our modeling is inaccurate and that the actual α is not −1.<br />

For the designed controller in (a), what is the range of the parameter<br />

α that results in zero steady-state tracking error for any step reference<br />

input?<br />

(10pt)<br />

For the designed controller parameters K = 1 and K a = −1, α must be<br />

less than one to keep the closed-loop stability.<br />

2


4. Solve the following continuous-time infinite-horizon LQR problem: (10pt)<br />

min<br />

u(·)<br />

∫ ∞<br />

0<br />

subject to<br />

(<br />

x<br />

2<br />

1 (t) + u 2 (t) ) dt,<br />

[<br />

ẋ1 (t)<br />

ẋ 2 (t)<br />

]<br />

=<br />

[ 0 1<br />

0 −1<br />

] [<br />

x1 (t)<br />

x 2 (t)<br />

] [ 0<br />

+<br />

1<br />

]<br />

u(t).<br />

Verify the stability of the closed-loop system.<br />

Algebraic Riccati equation:<br />

[ √ ] 3 1<br />

A T P + P A − P BR −1 B T P + Q = 0 ⇒ P = √ > 0<br />

1 3 − 1<br />

u(t) = −R −1 B T P x(t) = − [ 1<br />

√ 3 − 1 ] x(t)<br />

The closed-loop system is stable since<br />

A cl := A − BR −1 B T P =<br />

[ 0 1<br />

−1 − √ 3<br />

]<br />

3

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!