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.

10.2 Constrained Optimal <strong>Control</strong> Problem Formulation 165<br />

where Px,u is a polyhedron in R n+m or mixed input and state constraints over a finite<br />

time<br />

[x(0) ′ , . . . , x(N − 1) ′ , u(0) ′ , . . . , u(N − 1) ′ ] ∈ Px,u,N<br />

(10.27)<br />

where Px,u,N is a polyhedron in R N(n+m) . Note that constraints of the type (10.27)<br />

can arise, for example, from constraints on the input rate ∆u(t) u(t) −u(t −1). In<br />

this chapter, for the sake of simplicity, we will use the less general form (10.25).<br />

Define the cost function<br />

N−1 <br />

J0(x(0), U0) p(xN) + q(xk, uk) (10.28)<br />

where xk denotes the state vector at time k obtained by starting from the state<br />

x0 = x(0) and applying to the system model<br />

k=0<br />

xk+1 = Axk + Buk<br />

(10.29)<br />

the input sequence U0 [u ′ 0, . . . , u ′ N−1 ]′ .<br />

If the 1-norm or ∞-norm is used in the cost function (10.28), then we set<br />

p(xN) = PxNp and q(xk, uk) = Qxkp + Rukp with p = 1 or p = ∞ and P,<br />

Q, R full column rank matrices. Cost (10.28) is rewritten as<br />

N−1 <br />

J0(x(0), U0) PxN p + Qxkp + Rukp<br />

k=0<br />

(10.30)<br />

If the squared euclidian norm is used in the cost function (10.28), then we set<br />

p(xN) = x ′ N PxN and q(xk, uk) = x ′ k Qxk + u ′ k Ruk with P 0, Q 0 and R ≻ 0.<br />

Cost (10.28) is rewritten as<br />

J0(x(0), U0) x ′ N−1 <br />

NPxN +<br />

k=0<br />

x ′ kQxk + u ′ kRuk<br />

Consider the constrained finite time optimal control problem (CFTOC)<br />

J ∗ 0 (x(0)) = minU0 J0(x(0), U0)<br />

subj. to xk+1 = Axk + Buk, k = 0, . . .,N − 1<br />

xk ∈ X, uk ∈ U, k = 0, . . .,N − 1<br />

xN ∈ Xf<br />

x0 = x(0)<br />

(10.31)<br />

(10.32)<br />

where N is the time horizon and Xf ⊆ R n is a terminal polyhedral region. In (10.28)–<br />

(10.32) U0 = [u ′ 0 , . . . , u′ N−1 ]′ ∈ R s , s mN is the optimization vector. We denote<br />

with X0 ⊆ X the set of initial states x(0) for which the optimal control problem<br />

(10.28)–(10.32) is feasible, i.e.,<br />

X0 = {x0 ∈ R n : ∃(u0, . . . , uN−1) such that xk ∈ X, uk ∈ U, k = 0, . . .,N − 1, xN ∈ Xf<br />

where xk+1 = Axk + Buk, k = 0, . . .,N − 1},<br />

(10.33)

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

Saved successfully!

Ooh no, something went wrong!