30.07.2014 Views

6 Controllability and observability

6 Controllability and observability

6 Controllability and observability

SHOW MORE
SHOW LESS

Transform your PDFs into Flipbooks and boost your revenue!

Leverage SEO-optimized Flipbooks, powerful backlinks, and multimedia content to professionally showcase your products and significantly increase your reach.

Chapter 5 Stability<br />

• An unstable system is useless in practice <strong>and</strong><br />

stability is a basic requirement for all systems.<br />

• 5.2 Input-Output stability of LTI Systems<br />

• Consider a SISO LTI system<br />

y(t)<br />

=<br />

t<br />

∫0 g(t − τ)u(<br />

τ)dτ<br />

= ∫0<br />

g( τ)u(t<br />

− τ)<br />

dτ<br />

• An input u(t) is said to be bounded if<br />

u(t)<br />

≤ u m < ∞<br />

t


• A system is said to be BIBO stable<br />

(bounded-input bounded-output) if every<br />

bounded input excites a bounded output.<br />

• Theorem 5.1 A SISO system is BIBO stable<br />

if <strong>and</strong> only if g(t) is absolutely integrable in<br />

[0, ∞), or<br />

∫0<br />

∞<br />

g(t) dt<br />

≤<br />

M<br />

< ∞<br />

• A function that is absolutely integrable may<br />

not be bounded or may not approach zero as<br />

t → ∞.


• For example:<br />

f (t<br />

−<br />

⎪⎧<br />

n + (t − n)n<br />

n) = ⎨<br />

4<br />

⎪⎩ n − (t − n)n<br />

4<br />

for n -1/n ≤ t ≤ n<br />

for n < t ≤ n + 1/n<br />

The area under each triangle is 1/n.<br />

Thus the absolute integration of the function<br />

equals 2 (1/ 2<br />

∑ ∞ n=<br />

n ) < ∞.<br />

• This function is absolutely integrable but is<br />

not bounded <strong>and</strong> not approach zero as t→∞.<br />

3<br />

3


• Theorem 5.2 If a system with impulse<br />

response g(t) is BIBO stable, then, as t→∞:<br />

1. The output excited by u(t) = a, for t ≥ 0,<br />

approaches ĝ(0)⋅a.<br />

2. The output excited by u(t) = sinω 0 t, for t<br />

≥ 0, approaches<br />

ĝ( jω0<br />

) sin( ω0t<br />

+ ∠ĝ( jω0))<br />

where ĝ(s) is the Laplace transform of g(t).


• Theorem 5.3 A SISO system with proper<br />

rational transfer ĝ(s) is BIBO stable if <strong>and</strong><br />

only if every pole of ĝ(s) has a negative real<br />

part or, equivalently, lies inside the left-half<br />

s-plane.<br />

• Theorem 5.M1 A multivariable system with<br />

impulse response matrix G(t) = [g ij (t)] is<br />

BIBO stable if <strong>and</strong> only if every g ij (t) is<br />

absolutely integrable in [0, ∞).


• Theorem 5.M3 A multivariable system with<br />

proper rational transfer matrix Ĝ (s) = [ĝ ij (s)]<br />

is BIBO stable if <strong>and</strong> only if every pole of<br />

every ĝ ij (s) has a negative real part.<br />

• 5.2.1 Discrete-Time Case<br />

• Consider a discrete-time SISO system<br />

y[k]<br />

=<br />

k<br />

∑ m = 0<br />

g[k − m]u[m] =<br />

k<br />

∑ m = 0<br />

g[m]u[k − m]<br />

• A input sequence u[k] is said to be bounded<br />

if |u[k]| ≤ u m < ∞


• Theorem 5.D1 A discrete-time SISO system<br />

is BIBO if <strong>and</strong> only if g[k] is absolutely<br />

summable in [0, ∞) or<br />

∑<br />

∞ k=0<br />

g[k]<br />

≤<br />

M<br />

<<br />

∞<br />

• Theorem 5.D2 If a discrete-time system<br />

with impulse response sequence g[k] is<br />

BIBO stable, then, as k →∞:<br />

1. The output excited by u[k] = a, for k ≥ 0,<br />

approaches ĝ(1)·a.


2. The output excited by u[k] = sinω 0 k, for k ≥<br />

0, approaches<br />

jω<br />

ĝ(e<br />

) sin( ω k + ∠ĝ(e<br />

0<br />

0 j<br />

0<br />

where ĝ(z) is the z-transform of g[k].<br />

• Theorem 5.D3 A discrete-time SISO system<br />

with proper rational transfer function ĝ(z) is<br />

BIBO stable if <strong>and</strong> only if every pole of ĝ(z)<br />

has a magnitude less than 1, or equivalently,<br />

lies inside the unit circle on the z-plane.<br />

ω<br />

))


• Theorem 5.MD1 A MIMO discrete-time<br />

system with impulse response sequence<br />

matrix G[k] = [g ij [k]] is BIBO stable if <strong>and</strong><br />

only if every g ij [k] is absolutely summable.<br />

• Theorem 5.MD3 A MIMO discrete-time<br />

system with discrete proper rational transfer<br />

matrix Ĝ (z) = [ĝ ij (z)] is BIBO stable if <strong>and</strong><br />

only if every pole of every ĝ ij (z) has a<br />

magnitude less than 1.


5.3 Internal stability<br />

• The BIBO stability is defined for the zerostate<br />

response.<br />

• Internal stability: the stability of zero-input<br />

response or the response of<br />

x &(t)<br />

=<br />

Ax(t)<br />

That is, x(t) = e At x 0 .


• Definition 5.1 The zero-input response or<br />

the equation is marginally stable or stable in<br />

the sense of Lyapunov if every finite initial<br />

state x 0 excites a bounded response.<br />

• It is asymptotically stable if every finite<br />

initial state excites a bounded response,<br />

approaches 0 as t →∞.<br />

• Theorem 5.4<br />

1. The equation x & = Ax is marginally stable if<br />

<strong>and</strong> only if all eigenvalues of A have zero or<br />

negative real parts.


2. The equation x & = Ax is asymptotically stable<br />

if <strong>and</strong> only if all eigenvalues of A have<br />

negative real parts.<br />

• 5.3.3 Discrete-Time Case<br />

• The internal stability of discrete-time<br />

system: the stability of x[k+1]=Ax[k].<br />

• Its solution is x[k] = A k x 0 .


• Theorem 5.D4<br />

1. The equation x[k+1] = Ax[k] is marginally<br />

stable if <strong>and</strong> only if all eigenvalues of A<br />

have magnitude less than or equal to 1.<br />

2. The equation x[k+1] = Ax[k] is<br />

asymptotically stable if <strong>and</strong> only if all<br />

eigenvalues of A have magnitude less than<br />

1.


5.4 Lyapunov Theorem<br />

• Theorem 5.5 all eigenvalues of A have<br />

negative real parts if <strong>and</strong> only if for given<br />

positive definite symmetric matrix N, the<br />

Lyapunov equation A’M + MA = -N has a<br />

unique symmetric solution M <strong>and</strong> M is<br />

positive definite.


• Corollary 5.5 All eigenvalues of an n×n<br />

matrix A have negative real parts if <strong>and</strong><br />

only if for any given m×n matrix N with m <<br />

n <strong>and</strong> with the property<br />

⎡ N<br />

⎢<br />

NA<br />

rank O : = rank⎢<br />

⎢ .<br />

⎢ n<br />

⎣NA<br />

−1<br />

⎤<br />

⎥<br />

⎥ = n (full column rank)<br />

⎥<br />

⎥<br />

⎦<br />

where O is an nm×n matrix, the Lyapunov<br />

equation<br />

A′<br />

M +<br />

MA<br />

= −N′<br />

N<br />

= : −N<br />

has unique symmetric solution M <strong>and</strong> N is<br />

positive


• Theorem 5.6 If all eigenvalues of A have<br />

negative real parts, then the Lyapunov<br />

equation A’M + MA = -N has a unique<br />

solution for every N, <strong>and</strong> the solution can be<br />

expressed as<br />

M<br />

= ∫0<br />

∞<br />

e<br />

A't<br />

Ne<br />

At<br />

dt


5.4.1 Discrete-Time case<br />

• Consider the Lyapunov equation<br />

M - AMB = C<br />

where A: n×n, B: m×m M <strong>and</strong> C: n×m<br />

• As (3.60), the above equation can be expressed<br />

as Ym = c.<br />

• Let η k be an eigenvalue of Y. Then we have<br />

η k = 1 - λ i µ j<br />

where λ i <strong>and</strong> µ j are the eigenvalues of A <strong>and</strong> B


• Let A(M) := M - AMB. Let Au = λ i u <strong>and</strong><br />

vB = vµ j . A(uv) = uv - AuvB = (1 - λ i µ j )uv.<br />

• The eigenvalues of the above equation are<br />

1 - λ i µ j<br />

• If there are no i <strong>and</strong> j such that λ i µ j = 1, then<br />

(5.25) is nonsigular <strong>and</strong>, for any C, a unique<br />

solution M exists.


• Theorem 5.D5 All eigenvalues of an n×n<br />

matrix A have magnitudes less than 1 if <strong>and</strong><br />

only if the discrete Lyapunov equation<br />

M - A’MA = N<br />

has a unique symmetric solution M <strong>and</strong> N is<br />

positive definite.<br />

• Theorem 5.D6 The solution of the discrete<br />

Lyapunov equation can be expressed as<br />

M<br />

= ∞ m=0<br />

∑ (A ′ )<br />

m NA<br />

m


5.5 Stability of LTV Systems<br />

• Consider a SISO LTV system<br />

y(t) = ∫t t g(t, τ)u(<br />

τ)<br />

dτ<br />

0<br />

• The system is said to be BIBO if every<br />

bounded input excites a bounded output.<br />

∫<br />

t<br />

t 0<br />

g(t, τ) dτ<br />

≤ M < ∞<br />

• For the multivariable case,<br />

y(t) = ∫t t G(t, τ)<br />

dτ<br />

≤ M < ∞<br />

0


• The zero-state response<br />

G(t, τ) = C(t)Φ(t, τ)B(τ) + D(t)δ(t-τ)<br />

is BIBO stable if <strong>and</strong> only if there exist<br />

constants M 1 <strong>and</strong> M 2 such that<br />

D(t)<br />

≤ M1 < ∞,<br />

t<br />

t<br />

∫ G(t, τ)<br />

dτ<br />

≤ M2<br />

0<br />

<strong>and</strong><br />

< ∞<br />

• The zero-input response is marginally stable<br />

if <strong>and</strong> only if<br />

Φ( t, t0)<br />

≤<br />

M<br />

< ∞


• The equation x & = A(t)xis asymptotically stable<br />

if the response excited by every finite initial<br />

state is bounded <strong>and</strong> approaches zero as t →<br />

∞.<br />

• Theorem 5.7 Marginal <strong>and</strong> asymptotic<br />

stabilities of x & = A(t)x are invariant under any<br />

Lyapunov transformation.

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

Saved successfully!

Ooh no, something went wrong!