Abstract

For quasi-linear parameter varying (quasi-LPV) systems with bounded disturbance, a synthesis approach of dynamic output feedback robust model predictive control (OFRMPC) with the consideration of input saturation is investigated. The saturated dynamic output feedback controller is represented by a convex hull involving the actual dynamic output controller and an introduced auxiliary controller. By taking both the actual output feedback controller and the auxiliary controller with a parameter-dependent form, the main optimization problem can be formulated as convex optimization. The consideration of input saturation in the main optimization problem reduces the conservatism of dynamic output feedback controller design. The estimation error set and bounded disturbance are represented by zonotopes and refreshed by zonotopic set-membership estimation. Compared with the previous results, the proposed algorithm can not only guarantee the recursive feasibility of the optimization problem, but also improve the control performance at the cost of higher computational burden. A nonlinear continuous stirred tank reactor (CSTR) example is given to illustrate the effectiveness of the approach.

1. Introduction

In control system design, input saturation nonlinearity can severely degrade closed-loop system control performance and may lead to the instability of a closed-loop system [1, 2]. Hence, stability analysis of controller design with the consideration of input saturation has received extensive attention. An approach to deal with input saturation is to penalize the control input such that input constraint is never violated. This approach is common in the synthesis approach of MPC; for example, see [39]. In the synthesis approach of MPC, at each sampling time, the on-line optimization problem is solved to obtain an optimal controller, which considers physical constraints (e.g., input constraint and output constraint) and stability condition. Here, we refer to MPC with guaranteed stability as the synthesis approach, which is usually based on recursive feasibility; that is, the optimization problem is feasible for all time if it is feasible at the initial time [8, 10].

Linear parameter varying (LPV) systems offer a promising framework for modeling and control of a large class of nonlinear systems [11]. By applying plenty of methods and techniques developed for linear systems, LPV models provide efficient ways to deal with some complex nonlinear systems [1214]. When the scheduling parameters of LPV systems are exactly known at the current time but unknown in future, it is quasi-LPV system [79]. In robust MPC (RMPC) studies, a real nonlinear system is often approximated or included by polytopic uncertainty and then represented by LPV description. The synthesis approach of RMPC is often solved by on-line min–max optimization, which considers all the possible realizations of polytopic model parametric uncertainty, constraints, and/or not bounded disturbance; for example, see [3, 5, 79]. For the control of nonlinear systems subject to uncertainty, constraints, and faults in the control actuators, the authors in [15, 16] investigate a RMPC design for fault-tolerant control. Furthermore, when true states are unmeasurable, an output feedback MPC controller based on state observer is also considered in [16]. The feedback controller design in the synthesis of RMPC is often enforced by stringent constraints, which is usually conservative for the optimized controller. To reduce the conservatism of the controller design, the authors in [17] take a parameter-dependent controller, which improves the control performance in [3]. Free control moves (the control moves are the decision variables in on-line optimization problem) are introduced in [46] for state feedback controller design, which can also improve the control performance in [3]. Furthermore, the time-varying parameters of LPV systems having bounds on their rate of variation are considered in [5, 6], which further reduce the conservatism of the controller design. By introducing one free control move, the authors in [18] generalize the procedure in [4] to the case when the state is unmeasurable.

In [79], the dynamic OFRMPC approach considers quasi-LPV systems with bounded disturbance. By taking controller parameters with a parameter-dependent form, the main optimization problem is solved as convex optimization. The estimation error sets in [79] are represented by polyhedral sets (in [9], the estimation error set is represented by a zonotope, which is a special case of polyhedral set). In [7, 9], at each sampling time, the auxiliary optimization problem can not only refresh the bounds of estimation error set, but also determine whether to solve the main optimization at the next sampling time or not. However, the iterative auxiliary optimization problem in [7] increases the on-line computational burden, and the recursive feasibility of the optimization problem in [7, 9] is not guaranteed. The algorithm in [8] can guarantee the recursive feasibility of the optimization problem, while the control performance is worse than that in [9]. The estimation error set in [9] is refreshed by zonotopic set-membership estimation, which improves the control performance in [7, 8]. Furthermore, in [79], an ellipsoidal bounded disturbance set is considered in the main optimization problem. However, an outer approximation polyhedral set of the ellipsoidal bounded disturbance is utilized in the process of refreshing estimation error set, which may introduce conservatism.

To reduce the conservatism of OFRMPC controller design in [79] and then improve the control performance, the present paper considers a synthesis approach of dynamic OFRMPC with input saturation for the quasi-LPV system with bounded disturbance. Different from [15, 16], we assume that the system is under nominal operation and control actuator faults are not considered. The saturated dynamic output feedback controller is expressed by a convex hull involving the parameter-dependent form of both the actual dynamic output feedback controller and the auxiliary feedback controller, which is different from the saturated output feedback controller in [1, 19, 20]. The parameter-dependent form of controller can help to formulate the main optimization problem as convex optimization. Then the main optimization problem can be solved by a linear matrix inequality (LMI) tool. In the main optimization problem, only the constraints on the auxiliary feedback controller are considered, which introduce more optimization freedom for the the actual dynamic output feedback controller. By this way, the conservatism of output feedback controller design is reduced such that the control performance can be improved. In the process of refreshing estimation error set, based on the estimation error equation, the polyhedral estimation error set is refreshed by zonotopic set-membership estimation. Then the obtained polyhedral estimation error set compares with the ellipsoidal estimation error set obtained from the invariance condition of the augmented closed-loop system. By properly refreshing the bounds of estimation error set, it not only obtains the bounds of polyhedral estimation error set, but also guarantees the recursive feasibility of the optimization problem. Furthermore, different from [79], both the main optimization problem and the refreshment of estimation error set consider componentwise bounded disturbance, which is confined by a polyhedral set and also can be represented by a zonotope. Hence, the outer approximation of bounded disturbance in the process of refreshing estimation error set is not involved. Compared with [8, 9], the consideration of input saturation in the main optimization problem improves the control performance at the cost of higher on-line computational burden.

Notations. For any vector and positive-definite matrix , . is the value of at time , predicted at time . is the identity matrix with appropriate dimension. denotes the ellipsoid associated with the symmetric positive-definite matrix . All vector inequalities are interpreted in a componentwise sense. An element belonging to means that it is a convex combination of the elements in , with the scalar combing coefficients being nonnegative and their sum is equal to 1. The symbol “” induces a symmetric structure in the matrix inequalities. A value with superscript “” means that it is the optimal solution of the optimization problem. The time dependence of the MPC decision variables is often omitted for brevity.

Let an interval be defined as the set . The unitary interval is . A box is an interval vector. A unitary box in , denoted by , is a box composed by unitary intervals. The Minkowski sum of two sets and is defined by .

Definition 1. An -zonotope in can be defined as the linear image of an -dimensional hypercube in , where is the order of the zonotope and satisfies . Given a vector and a matrix , an -zonotope is the following set:

This is the Minkowski sum of the -segments defined by columns of matrix in . The center of zonotope is vector ; , with , are the generators of zonotope . Zonotopes are the special case of polyhedral sets and can be represented by the convex combination of their vertices.

2. Problem Statement

Consider the following discrete-time uncertain LPV system:where ,  ,  , and are the input, state, output, and disturbance, respectively. The disturbance is componentwise bounded satisfies ; that is, is bounded in a polyhedral set, which can be represented by the following zonotope:The output constraint is required to satisfywhere , , , and . The function is a vector valued standard saturation function; that is, , and , . Here, we have slightly abused the notation by using to denote both the scalar valued and vector valued saturation functions. The time-varying system matrices are known to vary within a polytope ; that is, , and there exist nonnegative coefficients ’s with and , where scheduling parameters ’s are exactly known at the current time but unknown in future.

For above system (2), the dynamic output feedback controller is of the following form:where is the controller state; are the controller gain matrices; are the feedback gain matrices.

To deal with the saturated feedback control, the method in [1, 20] is utilized. At each time , for matrices and , we definewhere and represent the th row of and , respectively. represents the region in , where the auxiliary feedback controller does not saturate.

Furthermore, denote as the set of diagonal matrices whose diagonal elements are either 1 or 0. There are elements in . Suppose that the elements of are labeled as . Denote ,  , and . Clearly, if . The following lemma adopted from [1, 20] is utilized to represent the saturated output feedback controller by a convex hull.

Lemma 2. For any ,  , the saturated input can be represented by the convex hull of the actual controller and the auxiliary controller as

By applying Lemma 2, the saturated output feedback controller, with , can be expressed aswhere ’s are exactly known at the current time and satisfy .

The above controller parameters take the parameter-dependent form asHence, the augmented closed-loop system, based on (5), (6), (9), (10), and the predictions made by system (2), is where

At each time , the true state is unmeasurable. Hence it is necessary to use its bounds to represent the true state. Suppose at time that is confined by -order zonotope represented by Then the true state can be represented by the following convex set:where , are the vertices of ; is the total number of vertices of zonotope at time .

Remark 3. The method in [9] can be applied to obtain the vertices of ; that is, calculate vertices of zonotope , and apply Quickhull Algorithm in [21] to remove the redundant vertices; then vertices of can be obtained.

Definition 4 (robust positively invariant set). For the following dynamical system,where , and bounded disturbance . If for all , , the condition holds for all , then the convex set is said to be robust positively invariant set.

3. Main Optimization Problem and Constraints Handling

3.1. Main Optimization Problem

At each time , for , the following main optimization problem with the consideration of input saturation to be solved iswhere is the performance cost; (17) is the stability/optimality condition; (18) is the augmented state constraint; (19) guarantees the auxiliary feedback controller does not saturate; (20) is the output constraint; and are the positive-definite weighting matrices; satisfies the following conditions:Choose ; then , , and .

Remark 5. The saturated output feedback controller with the parameter-dependent form is considered in problem (16)–(20), which is different from the main optimization problem in [79]. Furthermore, the present paper considers componentwise bounded disturbance. Hence the methods to deal with the constraints in the main optimization problem are different from those in [79].
Similar to the constraint on the auxiliary feedback controller in [2, 22, 23], the satisfaction of (19) guarantees that the auxiliary feedback controller does not saturate such that the saturated output feedback controller can be represented by the convex hull as (9). The constraint on the actual dynamic output feedback controller is not considered in the main optimization problem. This enables us to design with fewer restrictions. By this way, the conservatism of actual dynamic output feedback controller design can be reduced.

3.2. Constraints Handling

In the section, the following methods are utilized to deal with (17)–(20), which consider saturated output feedback controller with the parameter-dependent form and componentwise bounded disturbance.

Lemma 6. Equation (17) is satisfied if there exist scalars , matrices , positive dialog matrix , and symmetric matrices such thatFurthermore, if (22) is feasible, then (10) can be parameterized as

Proof. By adopting (24), it is equivalent toIf and , , such that , then, based on Definition 4, the ellipsoidal set is a robust positively invariant set for system (10). This condition is satisfied if there exist and a positive dialog matrix such thatClearly, (26) is satisfied for all and disturbance if the following conditions are satisfied:By applying the Schur complement, it is shown that (27) is guaranteed by Let and let . By pre- and postmultiplying the left-hand side of (29) with and , respectively, considering the convexity of the polytopic description of system (10) and then applying (25), one can obtain (22).

Denote the estimation error signal as . Based on the augmented state constraint in [8, 9], (18) holds if there exist positive-definite matrices such thatThe estimation error can be represented by the polyhedral set; that is, , where , are the vertices of ; is the total number of vertices related to . Hence, similar to the polyhedral estimation error constraint in [8, 9], (30) is guaranteed byFurthermore, the estimation error set can also be confined by an ellipsoid; that is, , where is known at time . In this case, based on the ellipsoidal estimation error set in [8], (30) is guaranteed by

Lemma 7. Suppose at time that (22), (23), (30), and (31) are satisfied. If (34) is satisfied,then (19) is guaranteed.

Proof. Since (22), (23), (30), and (31) are satisfied, it can be inferred that . Define as the th row of the -order identity matrix. ThenIf (36) is satisfied,then . By applying the Schur complement, (10), and the convexity of the polytopic description of system (11), it is shown that (36) is equivalent toBy pre- and postmultiplying the left-hand side of (37) with and , respectively, and then applying (25) and ,  , one can obtain (34).

Lemma 8. Suppose at time that (22), (23), (30), and (31) are satisfied. If there exists a positive-definite matrix such thatwhere is the th diagonal element of , then (20) is satisfied.

Proof. Define as the th row of the -order identity matrix. Then,where . If there exists a symmetric matrix satisfying (39) such thatthen . By applying the Schur complement, it is shown that (41) is equivalent toBy pre- and postmultiplying the left-hand side of (42) with and , respectively, considering the convexity of the polytopic description of system (11), and then applying (25), (38) is obtained.

By applying “Proposition  2” in [24], when the complexity parameter in [24] is , the sufficient conditions for (22) and (38) to hold areWith larger , less conservative conditions can be obtained, but the computational burden is heavier. By applying “Proposition  2” in [24] and considering the convex representation of the input saturation, the number of LMIs incurred by (43) is and that by (44) is .

By summarizing the above derivations, it is shown that when the estimation error is represented by the polyhedral set , problem (11)–(15) can be approximated byWhen the estimation error is confined by the ellipsoidal set , problem (11)–(15) can be approximated byProblems (45) and (46) can be solved by an LMI tool for the fixed scalar over . When problem (45) or (46) is solved at time , the controller parameters are calculated by (10) and (24). The representation of saturated output feedback controller by the parameter-dependent form of both the actual output feedback controller and the auxiliary feedback controller increases the decision variables and also the LMI conditions to describe constraints. Then the computational burden is increased (see Section 6).

3.3. The Recursive Feasibility

At each time , based on Proposition  3 in [8], when problem (45) or (46) is solved, the optimal solution for problem (45) or (46) can be obtained. At time , by choosingthe feasible solution for problem (46) can be constructed.

Furthermore, in problem (45) or (46), (43) guarantees that augmented state , will converge within the ellipsoidal invariant set . Since (43) considers all the possible realizations of polytopic model parametric uncertainty and bounded disturbance, leads to ; that is,Based on the selected parameters in (48), it can be obtained that (50) is equivalent tothat is, based on the invariance of the augmented closed-loop system, the estimation error set at time is confined by an ellipsoidal set . Note that the true estimation error must satisfy (51), which has considered all the possible realizations of polytopic model parametric uncertainty and bounded disturbance. Hence, it may be conservative for the estimation error constraint at time . If we can obtain the bounds of the polyhedral set and satisfy , then (32) is less conservative than (33) for the estimation error constraint at time . Solving problem (45) at time can obtain better control performance than solving problem (46). If , then exceeds the bounds of . In this case, solve problem (46) at time , and apply a central symmetric polyhedral set for outer approximating . By this way, the estimation error set at time is also a polyhedral set. In the following section, we apply the zonotopic set-membership estimation approach for refreshing the polyhedral estimation error set.

4. The Refreshment of Estimation Error Set

At each time , for quasi-LPV system (2), the model parameters are known at the current time . The estimated state is calculated by (5). When the control input is calculated by (6) and implemented to system (2), then, according to (2) and the definition of estimation error, it is shown thatBy applying (3), (13), and zonotopic computation properties  1 and  2 in [9], it can be obtained thatFrom (53), the estimation error set at time can be represented by -order zonotope , whereIt can be obtained from (53) that, with the evolution of time, the order of zonotope will increase. Accordingly, the generators and the vertices of zonotope will increase, which may lead to more LMIs in the main optimization problem to describe the estimation error constraints and then increase the on-line computational burden of the main optimization problem. Hence, we impose the restriction on the order of by . When , Property  3 in [9] is applied to construct a low-order -zonotope for outer approximating the high-order zonotope ; that is, , where, and are calculated according to (2) and (3) in [9], respectively. From the above derivations, at time , the estimation error set is of one of the following forms:When the estimation error set is obtained, it can be described by vertices; that is, .

At each time , when is obtained by (56) or (57) and is obtained by (48), it can be obtained that if the following problem is feasible,then it will satisfy the fact that . If problem (58) is infeasible, then . For this case, we apply the following method in [8] for outer approximation of by a polyhedral set.

Define the polyhedral set , where is a transformation matrix and ; , . Both and are off-line prespecified. Let . Then the set can perform the outer approximation of the set . In [8], is a central symmetric polyhedral set with the center of origin. Hence, it can be represented by a zonotope. By applying the zonotopic computation method in PnPMPC Toolbox [25], we denote it as .

Remark 9. In [79], the main optimization problem considers an ellipsoidal bounded disturbance, while in the refreshment of estimation error set, the outer approximation of the ellipsoidal bounded disturbance by the polyhedral set is utilized, which may introduce conservatism. However, the main optimization problem and the refreshment of estimation error set in the present paper consider the componentwise bounded disturbance, which is represented by the zonotope (see (3)). Hence, the outer approximation of bounded disturbance is not involved in the refreshment of estimation error set.

5. The Overall Solution to the Saturated Dynamic OFRMPC

The overall solution includes solving problem (45) or (46) and calculating the controller parameters by (10) and (24) then refreshing the estimation error set. We summarize the overall algorithm as follows.

Algorithm 10 (the overall solution to the saturated dynamic OFRMPC). Prespecify a scalar , , and . Choose a feasible , and , where . Select the maximal order of zonotope related to the estimation error set as . Let . At each time , perform the following steps: (1)If , then solve problem (45); else solve problem (46).(2)Calculate the controller parameters by (10) and (24). Implement the control input by (6). Obtain and set by (5) and (48), respectively.(3)Calculate by (53). If , then ; else .(4)If problem (58) is feasible, then let ; else let , perform outer approximation of by , represent it by zonotope , and let .

Theorem 11. For uncertain system (2), suppose that for all and . By applying Algorithm 10, if problem (45) is feasible at time , then the recursive feasibility is guaranteed. The augmented state will converge to a neighborhood of ; the input and the output constraints are satisfied for all .

Proof. As mentioned in Section 3.3, based on Proposition  3 in [8], when problem (45) or (46) is solved at time , the optimal solution is or . Based on (47), (48), and (49), the feasible solution for problem (46) at time can be obtained. Hence, if problem (45) is feasible at time , then, at time , consider the following two cases: (1)If problem (58) is feasible, then let , and solve problem (45).(2)If problem (58) is infeasible, then let , solve problem (46), and let .For the above two cases, when problem (45) or (46) is solved at time , it can guarantee that, in the main optimization problem, is always satisfied. Hence, based on the invariance of the augmented closed-loop system, the augmented state constraint at time is satisfied. By further choosing the feasible solution as (47) and (49), the recursive feasibility is guaranteed. In problem (45) or (46), (43) guarantees thatBy applying the definition of bounded disturbance and (28), the following condition is satisfied: By summing (60) from to and applying , it achievesHence, at each time , is a feasible bound of performance cost at time . By reoptimization of problem (45) or (46) at time , it must result in . By induction, will not increase with time .
Furthermore, consider the following uncorrupt system (i.e., the system is not corrupted by disturbance):Equations (23) and (43) also guaranteeBy summing (63) from to and then applying (62), it achieves Hence, is also an upper bound of performance cost for the uncorrupt system (62). means that ; that is, both and are bounded as a function of . The fact that tends to a constant value means that, with the evolution of time, will converge to , and will converge to a neighborhood of . The satisfaction of the input and output constraints is due to the input saturation and (39) and (44), respectively.

6. The Computational Burden for the Compared Algorithms

We compare Algorithm 10 with Algorithm  8 in [9] and Algorithm  1 in [8]. The complexity analysis for the compared algorithms is listed in Table 1. The complexity analysis for solving the main optimization problem by an LMI tool is polynomial time, which (regarding the fastest interior-point algorithms) is proportional to , where is the number of scalar LMI variables and is the number of scalar LMI rows (see [26]). From Table 1, it can be seen that the computational burden of the main optimization problem for Algorithm  8 in [9] and Algorithm  1 in [8] is different in the numbers of LMIs to describe the augmented state constraint. However, of the main optimization problem in Algorithm 10 is increased due to the introduced auxiliary controller. The saturated input is considered in the main optimization problem, which is related to the convex hull of the parameter-dependent form of the actual dynamic output feedback controller and the introduced auxiliary controller. Hence, the increase of LMIs to describe the constraints in the main optimization problem leads to the increased (see the terms with underline in Table 1).

7. Numerical Example

Consider the nonlinear continuous stirred tank reactor (CSTR) model, which has been widely studied (see, e.g., [79]). The nonlinear CSTR model is represented as the quasi-LPV system around the operation point, which considers the physical constraints and bounded disturbance. The details corresponding to system (2) are as follows: ; Furthermore, based on Lemma  2 in [20], we choose if ; else ; .

For comparing Algorithm 10 with Algorithm  8 in [9] and Algorithm  1 in [8], we choose , , , and . The initial state , where ; . The disturbance sequence satisfies (), where , are randomly generated from the interval , which can be represented by a zonotope . For Algorithm  8 in [9] and Algorithm  1 in [8], we choose and ; that is, (set in the present paper is the outer approximation of bounded disturbance in [8, 9]). For Algorithm  1 in [8] and Algorithm 10, we select , . Furthermore, for Algorithm  1 in [8], we let ; . For Algorithm 10 and Algorithm  8 in [9], we restrict ().

The simulation results are shown in Figures 1 and 2. Figure 1 compares the trajectories of the augmented closed-loop system. The control input signals shown in Figure 2 illustrate that the input constraints are satisfied. The time spent on the simulations is  s (41.13 s and 42.58 s) for Algorithm 10 (Algorithm  8 in [9] and Algorithm  1 in [8]). From the compared algorithms, Algorithm 10 in the present paper has the best control performance, but the computational burden is increased. The LMI Toolbox of Matlab 7.14 and zonotopic computation method in PnPMPC Toolbox [25] (AMD Phenom II Processor 2.70 GHz, 8 G Memory) are utilized for the simulations.

8. Conclusions

This paper considers a synthesis approach of dynamical OFRMPC for the quasi-LPV system with input saturation. The saturated control input is represented by a convex hull involving the parameter-dependent form of the actual output feedback controller and the introduced auxiliary controller. The main optimization problem is solved as convex optimization and reduces the conservatism of the dynamic output feedback controller design. The refreshment of estimation error set not only obtains the bounds of the estimation error set at the next sampling time, but also guarantees the recursive feasibility of the optimization problem. Compared with the previous results, the control performance is improved at the cost of higher computational burden.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (NSFC, 61403297) and the Fundamental Research Funds for the Central Universities (JB140417).