02.08.2013 Views

Model Predictive Control

Model Predictive Control

Model Predictive Control

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

242 12 Constrained Robust Optimal <strong>Control</strong><br />

Definition 12.1 (Robust Positive Invariant Set) A set O ⊆ X is said to be<br />

a robust positive invariant set for the autonomous system (12.1) subject to the<br />

constraints in (12.3), if<br />

x(0) ∈ O ⇒ x(t) ∈ O, ∀w(t) ∈ W, t ∈ N+<br />

Definition 12.2 (Maximal Robust Positive Invariant Set O∞) The set O∞ ⊆<br />

X is the maximal robust invariant set of the autonomous system (12.1) subject to<br />

the constraints in (12.3) if O∞ is a robust invariant set and O∞ contains all the<br />

robust positive invariant sets contained in X that contain the origin.<br />

Theorem 12.1 (Geometric condition for invariance) A set O ⊆ X is a robust<br />

positive invariant set for the autonomous system (12.1) subject to the constraints<br />

in (12.3), if and only if<br />

O ⊆ Pre(O, W) (12.34)<br />

The proof of Theorem 12.1 follows the same lines of the proof of Theorem 10.1. ✷<br />

It is immediate to prove that condition (12.34) of Theorem 12.1 is equivalent<br />

to the following condition<br />

Pre(O, W) ∩ O = O (12.35)<br />

Based on condition (12.35), the following algorithm provides a procedure for computing<br />

the maximal robust positive invariant subset O∞ for system (12.1)-(12.3)<br />

(for reference to proofs and literature see Chapter 10.3).<br />

Algorithm 12.1 Computation of O∞<br />

input fa , X, W<br />

output O∞<br />

let Ω0 ← X<br />

repeat<br />

k = k + 1<br />

Ωk+1 ← Pre(Ωk, W) ∩ Ωk<br />

until Ωk+1 = Ωk<br />

O∞ ← Ωk<br />

Algorithm 12.1 generates the set sequence {Ωk} satisfying Ωk+1 ⊆ Ωk, ∀k ∈ N and<br />

it terminates if Ωk+1 = Ωk so that Ωk is the maximal robust positive invariant set<br />

O∞ for system (12.1)-(12.3).<br />

Example 12.4 Consider the second order stable system in Example 12.1<br />

x(t + 1) = Ax(t) + w(t) =<br />

0.5 0<br />

1 −0.5<br />

<br />

x(t) + w(t) (12.36)

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

Saved successfully!

Ooh no, something went wrong!