Automatic Offline Formulation of Robust Model Predictive Control Based on Linear Matrix Inequalities Method

and Applied Analysis 3 The input and output constraints are satisfied if there exists symmetric matrixX and Y such that [ X ∗ Y T Q ] ≥ 0 with X rr ≤ u 2 r,max, r = 1, 2, . . . , m,


Introduction
Over the last several decades, model predictive control (MPC), for its ability to handle multivariable systems and constraints, has attracted notable attention in chemicals, food, automotive, aerospace applications, and other fields [1,2].MPC consists of a step-by-step optimization technique.It calculates more than one step control signal, but only the first one is implemented, and it repeats these calculations with the new measurements at the next sampling time.Since models used in MPC are only an approximation of the real process, it is extremely important for MPC to be robust to model uncertainty, and for this reason robust MPC is still under investigation [3,4].In this issue, most published research findings are based on the min-max approach, where the performance index to be minimized is computed over the worst case.However, the prohibitive online computation burden limits its application to relatively slow dynamics and small-scale processes.For this reason, the research on reducing the online computation is still a hot topic.
Generally, there are three kinds of strategies to reduce the online computation.One is based on the subdividing of the characteristic regions [5][6][7].The second is to design the invariant ellipsoids offline [8,9], and the third is reducing the free variables in the optimization progress [10][11][12][13].In the first method, the parameter space is systematically subdivided into characteristic regions, piecewise-affine feedback laws are computed offline, and a robust MPC is obtained using multiparametric quadratic programming [5][6][7].In the second strategy, a sequence of explicit control laws corresponding to a sequence of asymptotically stable invariant ellipsoids is constructed offline [8].It only performs the bisection search to find the smallest ellipsoid which contains the current state and applies the designed controller online.Based on the same philosophy, this offline approach is directly extended to the nominal system [9].The closed-loop system's feasibility and optimality are improved.Both of the offline methods in [8,9] are only for uncertain polytopic linear plants in the absence of external disturbance.A family of ellipsoidal sets is precomputed offline by exploiting the viability arguments for uncertain polytopic linear plants subject to persistent disturbance.The ellipsoidal sets are designed from the inner to outer layer, and the proposed procedure is robust in several senses and numerically affordable.The above two methods share the common feature of removing the online computation to offline.In the last method, an efficient robust predictive control (ERPC) is proposed [10] using the degrees of freedom on the control horizon to find large invariant ellipsoids for an augmented system offline, which leaves only a simple optimization problem to be solved online.An improved method for synthesizing the constrained robust model predictive controller is proposed in [12][13][14].They not only dramatically reduce the online computation burden, but also guarantee the control performance by reserving at least one free control move in the whole process.However, the conservative results are usually obtained, due to the fact that almost all the existing methods are based on the known state, and most of system's states are unknown or only partly known in practice.This can hamper the application of robust MPC.
In order to reduce the conservativeness, the output feedback control has been widely investigated.The state feedback robust MPC is combined with a state observer in [15], and the controller gain is obtained by a bisection search in a lookup table online.The system with saturated at transient or steady state is investigated in [16], and an MPC robust stability algorithm is proposed.The synthesis approach for output robust MPC with input/state constraints is proposed followed by the study of the quasi-linear parameter varying system [17,18].The output feedback controllers can reduce the system's conservatism.However, they may increase the online computational burden.So it is important to design a robust MPC method with low computational burden only based on the known initial state.
In this paper, we propose an offline MPC strategy for uncertain polytopic linear plants in the presence of input and output constraints.A nest of geometric proportion asymptotically stable ellipsoids is constructed which can assure the fixed convergence speed.The nested ellipsoids share the same character if the state is in an ellipsoid at the current; it must be shrinking into the inner ellipsoid with the designed controller.Moreover, an automatic offline formulation of robust constrained MPC is proposed, which not only almost reduces the online computation to zero, but is also only based on the known initial states at the expense of a possible reduction of the performance or the memory space requirements for storing the exact regions.Finally, the modified automatic offline robust MPC is proposed to improve the performance.
The rest of this paper is organized as follows.Section 2 gives the problem description and early results.In Section 3, the main results of this contribution are proposed.To illustrate the performance of the proposed controllers, simulation examples are presented in Section 4, and some conclusions are drawn in Section 5.
Notation.The notation used is fairly standard.  is the ndimensional space of real valued vectors.For a matrix  and a vector  ∈   , ‖‖ 2  =   .The matrix inequality  >  ( ≥ ) means that  and  are square symmetric matrices and − is (semi)positive definite.() (or ( | )) denotes the measured or actual value of variable  at real time ,  ∈ {0, 1, 2, . ..}; ( +  | ),  ∈ {0, 1, 2, . ..}, is the predicted value of  at a future prediction time  +  predicted at real time .I is the identity matrix with proper dimension.The symbol * denotes the corresponding transpose of the lower block part of symmetric matrices.For a real number , [] is the biggest integer not more than .{  } is the sequence with nonnegative integer .

Problem Description and Early Results
Consider the following time-varying uncertain system [4]  ( + 1) =  ()  () +  ()  () ,  ≥ 0 where  ∈   is the state of the plant,  ∈   is the control input, and The purpose is to design a controller that drives the system (1) with constraints (2) to the equilibrant state and minimize the worst case infinite horizon quadratic objective function: where () is the state feedback gain.Condition (6) should guarantee cost monotone decreasing and the system's robust stability.Summing ( 6) from  = 0 to  = ∞ the result of max [(+),(+)]∈Ω,≥0  ∞ () ≤ ‖()‖ 2 () ≤  can be obtained.With the definitions  = () −1 and () =  −1 , (6) and ‖()‖ 2 () ≤  are equivalent to the following LMIs: Abstract and Applied Analysis The aforementioned method is computationally prohibitive for fast and/or high-dimensional dynamic systems.Wan and Kothare [8] proposed the offline robust MPC strategy to reduce the online computation.They defined the asymptotically stable invariant ellipsoid and designed a sequence of invariant ellipsoids one inside another.A set of stabilizing state feedback laws corresponding to the ellipsoids is computed offline, so no optimization is involved online.The online computation is only the selecting of the proper feedback law using bisection search.This method supposes the state to be known, but in practice the state is usually unknown or partly known.So we propose our automatic offline MPC which only supposes the initial state is known.

Automatic Offline Robust MPC
In this section, we define the geometric proportion asymptotically stable ellipsoid and present an automatic offline approach based on it.The proposed method not only reduces the online computation dramatically, but also decreases the conservativeness.
Definition 1.Given a discrete dynamical system ( + 1) = (()), a subset E = { ∈   |    −1  ≤ 1} of the space   is said to be a geometric proportion asymptotically stable ellipsoid (GPASE), if it has the property that, whenever () ∈ E, then ( + 1) ∈ E  , where Remark 2. If the state of a system is in GPASE, the slowest convergence speed is more than .
Remark 3. The condition of system's state in GPASE increases the conservatism, but it can be decreased through selecting a bigger  near to 1.
Remark 6. Algorithm 5 is independent of the state.The control law is not optimal as it can be applied to all the states within the ellipsoid.The information of the exact smallest ellipsoid containing the current state is not needed before the application of the controller compared with the method in paper [8].
Lemma 7. Suppose the positive matrices satisfy the following condition: From above, it is easy to get E , ⊂ E +1 ⊂ E  .Theorem 8. Given a dynamical system (1) and a feasible initial state  0 satisfying ,  ∈ {1, . . ., }, the automatic offline robust MPC Algorithm 5 robustly asymptotically stabilizes the closed-loop system.
≤  2 }.From Lemma 7, the control law corresponding to the ellipsoid E  is guaranteed to keep the state within E  and converge into the ellipsoid E +1 , and so on.Lastly,   =   () is guaranteed to keep the state within E  and converge to the origin.
Compared with the strategy in paper [8], the proposed algorithm reduces the conservatism but may sacrifice some optimality.For example, according to Algorithm 5, the state is in a smaller ellipsoid, but the feedback control law may be computed based on a larger ellipsoid, which is showed in Figure 1.So we propose a modified algorithm.
Algorithm 9 (modified automatic offline robust MPC).Consider the same system as Algorithm 5. Given a sufficiently large feasible state  1 , generate a sequence of minimizers   ,   ,   ,   , and   as follows offline.Let  := 1.
(1) Compute the minimizer   ,   ,   ,   ,   at   by solving the problem (9) with additional LMIs (11) and The modified control strategy is selected online.Suppose the initial state ,  ∈ {1, . . ., }.Let the state be () at time .(3)  :=  + 1, select Remark 10.Algorithm 9 enforces the optimality an as it sets observation point.It ought to set more points to get more points to get the better performance.But the added points may increase the online computation and need more information of the system.The number of the observation points depends on the real system.How to select the proper number of observation points depends on the control engineering practice.
Remark 11.The online execution time for Algorithm 5 does not depend on the size of the system, so it has constant time (also written as O (1) time).For the method in paper [8], the total number of flops required to calculate an input move is ( 2  log 2 ) for a sequence of , with   being the number of state variables.So we can conclude that this offline approach can dramatically reduce the online computational burden.

Numerical Example
In this section, two examples are presented to illustrate the effectiveness of the automatic offline approach.The simulations were performed on a Sony PC with an Intel Core i5 3317U processor (1.7 GHz, memory 4 GB) and using the software LMI control Toolbox in the MATLAB 7.0 environment to compute the solution of the linear minimization problem.
We choose the  1 -axis as the one-dimensional subspace and discretize it according to the fixed geometric proportion, for example, with  = 0.96.  ,  = 1, 2, 3, . ... Figure 2 shows the ellipsoids defined by  −1  for all 150 discrete points.As it is too dense, an ellipsoid every other ten ellipsoids is selected and it is illustrated in Figure 3.
Figure 4 shows the computed offline control law along the  1 -axis.Both plots are obtained using the offline part of Algorithm 5.As the step  increases, the feedback almost keeps constant ( > 90).This is corresponding to the system's unconstrained control.
Given the initial state (0) = [0.5 0]  , Figure 5 shows the closed-loop responses of the system under the traditional model predictive control [4] (TMPC), offline model predictive control (OLMPC) [8], and automatic offline model predictive control (AOMPC) respectively.It shows that the AOLMPC, as it is based on the geometric proportion asymptotically stable ellipsoid, has the less overshoot.The average online computation time of the proposed algorithm compared with the other two methods is listed in Table 1.
The system's convergence speed at step  = 1, 10, 20, 30 is shown in Table 2.As 0.96 10 = 0.6648, it is concluded that the average convergence rate is more than 0.96 from Table 2.

Conclusions
An automatic stable closed-loop robust MPC strategy for uncertainty system with input and output constraints has been developed.This method not only reduces the conservativeness dramatically, but also reduces the online computation.The stability of this technique is proved based on GPASE.Two examples indicate the effectiveness of this method.

Figure 4 :Figure 5 :
Figure 4: The offline control law F for Example 12.

Table 1 :
The corresponding average computation time.