Triple-Mode Model Predictive Control Using Future Target Information

: In this paper, we propose a triple-mode model predictive control (MPC) algorithm that uses future target information to improve tracking performance. To explicitly take into account the future target information in the MPC optimization, the proposed triple-mode control law encompasses three parts: (i) the future target information feedforward, (ii) the output feedback, and (iii) the extra degrees of freedom for constraint satisfaction. The ﬁrst two parts of the control law are off-line designed through unconstrained MPC, and the optimal future trajectory horizon is obtained by golden section search based on the integral of squared error (ISE) criterion. The ﬁnal part is calculated by the on-line MPC algorithm aiming to satisfy constraints. Furthermore, we analyze the feasibility and convergence properties of the proposed algorithm. The method is demonstrated by the simulation of the shell fundamental control problem and also tested on the coordinated control problem in the power plant. The test results show that the proposed algorithm can increase tracking performance dramatically due to the proper selection of future trajectory horizon.


Introduction
Concern about multivariable constrained system control has become a central issue in the control domain. As a result, much research in recent years has focused on the development of advanced control algorithms that are process friendly. Model predictive control (MPC) has been widely and successfully applied in the process industry, primarily due to its superior performance and the capability of dealing with a multivariate constrained problem [1][2][3]. Since the step-response model is more convenient to obtain and more suitable for the process industry, this model has been successfully applied to numerous industrial processes by Aspen Technology, Honeywell Hi-Spec, and others. There has been an increase in the adoption of dynamic matrix control in MPC commercial software packages such as DMC+, SMC, RMPCT, HIECON [1].
The key idea of model predictive control is to use the historical information and model prediction to calculate the optimal input trajectory by minimizing the deviation between the process output and the desired target in a receding horizon way. There exists a rich supporting theory [4][5][6] for its success in industry. Moreover, the applications of model predictive control have already gone beyond petrochemical or chemical industries and extended to power plants [7], traffic, robots [8] and any other fields.
The research to date has tended to focus on the regulation problem (steer the system to the origin) [4] for facilitating the analysis and design. The regulation problem assumes the hypothesis that the future setpoint is fixed or changes instantaneously, thus future target information is rarely discussed in the current literature. Considering the setpoint changes, Limón et al. [9] adopted a modified cost function and a stabilizing extended terminal constraint in MPC optimization to obtain the feasible setpoints in case that the setpoints are not feasible. Their work is still preserving no consideration of using the future target information. The MPC optimization have the capability of considering the future target information [10], which may be helpful in control since it allows the controller to react to the setpoint changes in advance [11]. Carrasco and Goodwin [12] proposed a feedforward model predictive control, in which the feedforward was designed individually before the feedback design. However, how much future target information (so-called "preview" in their work) should be used in feedforward design is not discussed in their work. Rossiter and Grinnell [13] developed an extended input horizon generalized predictive controller (EIHGPC) algorithm, in which they parameterized the degrees of freedom to spread the control increments over the prediction horizon and adopted a heuristic algorithm to select a proper future target horizon, hence improves the tracking performance. Their method only pertains to the unconstrained case. In [14], the authors proposed a dual-mode constrained MPC algorithm, which was composed of a linear state feedback and extra degrees of freedom. Under such a control law, the feedforward could not be designed individually and the future target horizon was determined through a trial and error process.
However, the inappropriate use of future target information may cause the deterioration of control performance [13,15]. There remains a need for an effective method that can determine a proper future target horizon. Moreover, the current implementations of MPC fail to separately design the feedforward, feedback and constraint satisfaction, which could not yield the desirable performance from each part. In this study, a triple-mode constrained model predictive control algorithm is developed, which has the merit of independent design for each mode of the control law. The main ingredients of the triple-mode control law are: (i) the feedforward design using future target information, (ii) the output feedback, and (iii) the extra degrees of freedom for constraint satisfaction. In addition, a golden section search method is proposed to optimize the future trajectory horizon, which can be used to design the feedforward to improve the tracking performance.
The outline of this article is as follows. First, we describe the typical dynamic matrix control. Next, we give the future target information feedforward structure derived from the unconstrained control law and propose an optimization method for the selection of future trajectory horizon. We then introduce the triple-mode model predictive control using future target information. Finally, we provide two illustrative simulation examples.

Preliminary
Mathematical notations used throughout this paper are defined as follows. Given a vector x ∈ R n , the symbol ||x|| denotes the Euclidean 2-norm and ||x|| 2 P denotes the quadratic form of x. e.g., ||x|| 2 P = x T Px; A T denotes the transpose of A; Given two integers, l ≤ m, we define the set I l:m = {l, l + 1, ..., m − 1, m}. The superscript * denotes the optimal solution or cost according to the context.

Dynamic Matrix Control Algorithm
The design of the MPC algorithm in this paper follows the standard DMC approach proposed by Garcia and Morshedi [2]. Without loss of generality, we consider a MIMO linear dynamic system with n y ∈ R controlled variables and n u ∈ R manipulated variables. The DMC algorithm adopts a non-parametric step response model as the predictor, thus the ith output y i (t) ∈ R at time t can be formulated by the following step response model in which ∆u j (t) ∈ R denotes the jth input increment, i.e., ∆u j (t) := u j (t) − u j (t − 1). a i,j,k ∈ R denotes the kth element of the corresponding step response, and the step response sequence of the ith output to the jth input is {a i,j,1 , a i,j,2 , ..., a i,j,n , ...}. For a stable plant this sequence will asymptotically reach a constant value, i.e., a i,j,n+m ≈ a i,j,n , (m > 0). Based on the step response model, the output predictions over n future steps can be expressed as follows where a n y ,1,t a n y ,2,t · · · a n y ,n u ,t       , ∀t ∈ I 1:n ∆u(t) ∈ R n u is a vector of the future control moves at time t. Y(t) ∈ R n y ·n is the dynamic states of the system, and y(t + k|t − 1) ∈ R n y is the future output vector at time t + k assuming the input remain constant starting at time t − 1, which means ∆u(t + j) = 0 for j ≥ 0. The evolution of predictions can be obtained by translation of one-step prediction plus the contribution made by the latest control move ∆u(t).
To achieve the goal of offset-free tracking, a bias is introduced into the model to incorporate the feedback. It's assumed that the model prediction and the measured output at the current time is due to a constant step disturbance at the output. The constant output step disturbance assumption is standard in industrial applications [16]. The model bias e(t) ∈ R n y is based on the latest output measurement y(t) ∈ R n y and the predicted outputsŷ(t|t − 1) ∈ R n y e(t) =ỹ(t) −ŷ(t|t − 1) in whichŷ(t|t − 1) refers to the prediction of output at time t based on the measurements up to time t − 1. The bias is added to compensate the output prediction whereŶ(t|t − 1) denotes the estimate of Y(t) at time t before the latest measurements are obtained. Y(t|t) is the revised estimate of Y(t) on the basis of the feedback, and H := diag{I n y , ..., I n y } refers to the weighting matrix of the bias term. For open-loop stable systems, the future control moves ∆U M (t) ∈ R M·n u over the control horizon (M) is optimized to drive the model predicted outputsŶ P (t + 1|t) ∈ R P·n y as closely as possible to a desired future trajectory R(t + 1|t) ∈ R P·n y over prediction horizon (P). The DMC optimization problem is defined as where u max , u min , ∆u max , ∆u min are input constraints and increment constraints, respectively. The matrix A M P is the dynamic matrix with P × M block matrix and Λ P is the shift matrix with appropriate dimension. Q, W are the weighting matrices for outputs and the controller moves. The above notations readŶ At every sample time, controller evaluates its future input trajectories by solving the optimization problem (5) and sends the first step of its entire future input trajectory to its actuators. However, although the future target information R(t + 1|t) of the controlled variable were embedded in the optimization problem of DMC algorithm, the default choice in DMC algorithm was r(t + k|t) = r(t + 1|t), ∀k ∈ I 2:P . Little attention has been paid to the selection of an appropriate future trajectory horizon and the use of future target information.

Main Results
In this section, we first derive the explicit control law of the unconstrained MPC, which includes two parts: (i) the future target information feedforward, and (ii) the output feedback. Based on the future target information feedforward, we demonstrate that a poor future trajectory horizon has quite an impact on the system's performance. Then we propose an optimization method to determine the future trajectory horizon and the tripe-mode model predictive control using future target information. The triple-mode MPC is composed of the designed unconstrained control law and the extra degrees of freedom. The future target information is taken into account explicitly to improve tracking performance.

Illustrations of Poor Future Trajectory Horizon
Now we show that the unconstrained MPC control law will provide possibilities in the explicit expression of feedforward and feedback. Consider the performance objective J(t) in problem (5), one can derive the unconstrained law by minimizing it w.r.t the future control moves, according the extremum condition d( and making the following substitutions then we can get It can be seen that the future target information is involved in the first term of the control law: F r R(t + 1|t), which can be seen as the feedforward in the unconstrained control law. In addition, the second part of the control law can be seen as the output feedback based on the current output prediction.
For convenience, we define the future trajectory horizon in feedforward term for each output as f i ∈ R, ∀i ∈ I 1:n y , which indicates the length of future target information considered in the optimization. Any future targets beyond the future trajectory horizon will not be taken into account in the optimization problem, even if known. For a given f i , ∀i ∈ I 1:n y , the future target information adopted in MPC can be expressed as Furthermore, by combining the specific future trajectory horizon, the unconstrained law can be represented as follows: F r R(t + 1|t) = F 1 r 1 + ... + F n y r n y =F 1r1 + ... +F n yr n y (10) where F i r i , ∀i ∈ I 1:n y refers to the ith feedforward contribution, andF iri , ∀i ∈ I 1:n y is the reconstructed future target information feedforward component based on the designed future trajectory horizon f i .
, ∀i ∈ I 1:n y and F j , ∀j ∈ I 1:P·n y denotes the jth column vector of the feedforward matrix F r . It indicates that, for the ith output, the original P-dimension target vector r i has been replaced by a designed f i -dimension vectorr i . A major current focus in improving the control performance is how to select a proper f i , ∀i ∈ I 1:n y . It is known that the typical DMC algorithm considers no setpoint changes during the optimizaiton, which is equivalent to set the future trajectory horizon as f i = 1, ∀i ∈ I 1:n y .
Next, we examine the closed-loop performances of the MPC algorithm adopting different choices of future trajectory horizon. Consider a single-input, single-output system described by the following transfer function G(s) = 1 (10s + 1)(5s + 1) The sampling time is chosen as 5 s, a step change of setpoint happens at sampling instant 15 and the unconstrained MPC algorithm is applied to drive the system to its target. The arbitrary choices of future trajectory horizon are f = 1, 6, 15. The tuning parameters of DMC are: P = 15, M = 3 and Q = I, R = 5 * I in which I is the identity matrix.
As illustrated by Figure 1, with a default choice of future trajectory horizon in MPC ( f = 1), the system response is slugglish when the setpoint changes, in contrast, setting the future trajectory horizon as the prediction horizon could result in an aggressive response. A proper selection f = 6 deserves a desirable tracking performance. The integral of squared error (ISE) criterion is widely used to evaluate the system performance in control engineering, the discrete form of ISE reads in which error(t) refers to the difference between the expected value (target value) and the actual response value at time instant t.
The ISE values shown in Table 1 demonstrate that the standard MPC algorithm with no anticipation ( f = 1) owns large ISE value during the transition from one target to another. On the other hand, with large number of future trajectory horizon ( f = 15), the system would react far in advance and lead to bad performance too.

Future Trajectory Horizon Optimization Algorithm
Since the closed-loop response with smaller ISE value is better, we advocate use of the ISE value to select the proper future trajectory horizon. The ISE value can be regarded as a function of future trajectory horizon, i.e., ISE i = g( f i ), ∀i ∈ I 1:n y . A better f i deserves lower function value g( f i ), we take the idea of golden section search algorithm [17] to optimize the future trajectory horizon. The algorithm can be summarized as follows: According to Algorithm 1, we plot the ISE value subject to the iteration of a k , b k , λ k , µ k as Figure 2 shows. The optimal future trajectory horzion obtained is f = µ 6 = 6. Which testifies the result of last section.

Algorithm 1 The optimization of Future trajectory horizon based on ISE value
For any controlled variable y i , ∀i ∈ I 1:n y ; Step 1. Given the prediction horizon P, set a 1 = 1, b 1 = P, then calculate the initial point: and set k = 1; Step 2. Give a specified target, govern the process by using the control law Equation (8) and evaluate the ISE value g(·) over the entire runtime; Step and the optimal future trajectory horizon is obtained, break the optimization; Otherwise, if g(λ k ) > g(µ k ), go to Step 4, else if g(λ k ) ≤ g(µ k ), go to Step 5; Step 4. Set and go to Step 6; Step 5. Set and go to Step 6; Step 6. Set k = k + 1, go to Step 2.

Triple-Mode MPC Algorithm Using Future Target Information
For the conventional MPC algorithm, the optimization problem (5) is solved with n u · M degree of freedom ∆U M (t) to tackle constraints. However, the original optimization problem treats the input as decision variables, which breaks down the inherent future target information feedforward structure as described before. What's more, the default choice of future trajectory horizon f i = 1, ∀i ∈ I 1:n y in the MPC algorithm may lead to sluggish tracking performance. There exists no effective method to determine a proper future target horizon for MPC optimization.
With the proposed Algorithm 1 used, we next present a triple-mode MPC algorithm using future target information. Different from the typical DMC algorithm as described in preliminary section, the proposed triple-mode control law encompasses three parts: (i) the future target information feedforward, (ii) the output feedback, and (iii) the extra degrees of freedom. Each mode of the control law can be designed individually. The first two parts of the control law are off-line designed through unconstrained MPC, and the optimal future trajectory horizon is obtained by using Algorithm 1. The final part is calculated by the on-line MPC algorithm aiming to satisfy constraints.
Enhance the unconstrained MPC control law Equation (8) with an additional term C(t) ∈ R n u ·M in which C(t) = [c(t|t), c(t + 1|t), ..., c(t + M − 1|t)] T and c(t + i|t), i ∈ I 0:M−1 ∈ R N u representing degrees of freedom to get the constraints satisfied. It indicates that the original degree of freedom ∆U M (t) in constrained DMC optimization (5) is replaced by C(t), which not only provide the feedforward design from future target information but also guarantee the advantage of constraint satisfaction. At each sampling time, the controller solves the following optimization problem min C(t) where Equation (17f) refers to the terminal constraint used to ensure closed-loop stability when model predictive control is employed. After optimization, the optimal solution C * (t) = [c * (t|t), c * (t + 1|t), ..., c * (t + M − 1|t)] T can be obtained, and the optimal input increment is nothing but ∆u * (t) = d T (F r R(t + 1|t) − KŶ(t|t)) + c * (t|t) with d = [I n u , 0, ..., 0] T . Note that the proposed triple-mode MPC optimiazation problem (17) is a quadratic programming [18] with n u × M degrees of freedom. The degrees of freedom are the same as the original DMC optimization problem (5), thus the computation load of proposed triple-mode MPC is the same as which of DMC. A detailed implementation procedure of the proposed triple-mode MPC algorithm (Algorithm 2) is summarized as follows.

Algorithm 2 Triple-mode MPC algorithm using future target information
Offline computation: Carry out the Algorithm 1 to obtain the future trajectory horizon f i , ∀i ∈ R n y , the first two parts of the control law are designed through Equations (7) and (10); Online computation: Step 1. At sampling time t, controller receives measurement from the sensors and evaluate the modified estimate of predictions according to Equation (4); Step 2. Solve the optimization problem (17) to evaluate its future input trajectory; Step 3. Injects the optimal input into the plant; Step 4. Set t = t + 1, and go to Step 1.
In this section, the triple-mode MPC using future target information is presented. As has been previously stated, this triple-mode control law is based on the usage of the perturbations on the fixed dual-mode unconstrained control law. To this end, the following assumption is considered. Assumption 1. The step response sequence a i,j,k , ∀i ∈ I 1:n y , ∀j ∈ I 1:n u , ∀k ∈ N + is a Cauchy sequence: ∀ > 0, ∃N ∈ N, ∀n, ∀m, (m, n > N → |a i,j,n − a i,j,m | < ). Assumption 2. The weighting matrices Q ∈ R P·n y ×P·n y and R ∈ R M·n u ×M·n u are positive definite matrices. At least as many inputs as outputs is required: n u ≥ n y .
The future target is admissible if there exist a set of control sequence ∆U M (t) such that the system outputŶ P (t) can reach it and the constraints (e.g., ∆u min ≤ ∆U M (t) ≤ ∆u max , u min ≤ U M (t) ≤ u max ) can be satisfied during the stage. The following theorem follows immediately. Theorem 1. Consider that Assumption 1 and 2 hold, and the future target R(t + 1|t) is admissible. Then, for any feasible initial stateŶ(t|t), the proposed triple-mode controller ∆u * (t) = d T (F r R(t + 1|t) − KŶ(t|t)) + c * (t|t) asymptotically steers the system to R(t + 1|t).
The proof can be found in the Appendix A.

Applications
In this section, dynamic simulations are carried out to evaluate the performance of the proposed triple-mode controllers. This section presents two illustrative applications for the proposed methodology: (i) the shell fundamental control problem (case study I), and (ii) the coordinated control problem in the power plants (case study II). We aim at showing the better performance can be achieved by proper selection of future trajectory horizon. The simulation case study was performed on a PC with Intel Core i5 Quad 1.6GHz, 16GB memory and in the Matlab platform.

System I: A Heavy Oil Fractionator
In this example, consider a plant shown in Figure 3 presented in Prett and Morari [19]. in which a gaseous feed stream is separated by removing heat. Here we will only consider the two-by-two subsystem, which consists of two inputs: the top draw u 1 and the side draw u 2 ; two outputs: the top end point y 1 and the side end point y 2 . The model for this system is: The system will be represented using a 60-coefficient step response model with a 5-min sampling time. Set the initial operating point to be the origin and the upper bound of both output and input is 0.5, the lower bound of both output and input is −0.5. Three MPC algorithms are used to compare the system performance: the typical DMC, the proposed triple-mode MPC (Algorithm 2), and the DMC algorithm with an artificial selection of future target horizon (denoted as DMC-FF in the following). Tuning parameters of the MPC controllers are the same as follows: a prediction horizon of P = 60 and a control horizon of M = 8. The weighting matrices are set as identity matrix. In the DMC-FF algorithm, the future target horizon is selected as f 1 = 40, f 2 = 40. In the triple-mode MPC algorithm, the optimal future target horizon can be obtained as f 1 = 13, f 2 = 13 according to Algorithm 1. To compare the tracking performance, the evolution of setpoint is given as: the initial setpoint is coincident with the system's initial operating point; at sample time 20, change the setpoint sp 1 of output y 1 to its upper bound and the setpoint sp 2 of output y 2 is also changed to its upper bound when sample time 40; at sample time 60, set the setpoint of output y 1 to the origin and the setpoint of output y 2 is also set to the origin at sample time 80. Figures 4 and 5 present the output response and input response of the controlled system.  As it can be seen in Figure 4, both two outputs are sluggish under the typical DMC algorithm when the setpoint changed. On the other hand, using an artificial selection of future target horizon may cause the system to response to the future target information far in advance. The performance is degraded because the typical DMC algorithm overlooks the information of future target and the DMC-FF algorithm uses too much future target information. In contrast, the triple-mode MPC algorithm which taking into account the future target information properly can respond to subtlety, the performance of the system is improved dramatically.
Next, Table 2 shows the sum of the ISE value of the controlled system for each control strategy throughout the simulation horizon. The DMC-FF using the future target information can outperform the typical DMC, and the triple-mode MPC using an optimal future target horizon achieves the best tracking performance among the three algorithms.

System II: A Power Plant Coordinated Control System
Coordinated control system (CCS) is the system that coordinates the steam turbine and the boiler at the same time to respond to the various load demands quickly [20]. The boiler and the turbine possess different dynamic characteristics. The dynamic from the combustion rate to the main steam pressure is a slow process with large delay and great inertia. In contrast, the response from the turbine governor valve to the unit output is a relatively fast process.
The simplified coordinated control system can be shown in Figure 6. It can be seen that this system is a two-input, two-output control system, in which the manipulated variables are the fuel command u B (%) and the opening level of the turbine valve µ t (%), the controlled variables are the main steam pressure p t (MPa) and the output power N E (MW). Consider a 660MW unit in Datang Panshan Power Plant, China, in which the CCS is composed of the boiler of Type HG − 2023/17.6 − YM4 boiler and the turbine of Type N600 − 16.7/537/537 − I. The nonlinear model from the load 50% to load 100% was identified through the step test [21]. When the load N E = 550MW, the linearized model is: In this case study, it was established that the control system must always obey the speed restriction on the manipulated variables |∆µ t | ≤ 0.5%, |∆u B | ≤ 0.5%, and the system is initially at N E = 550 MW, p t = 16.9 Mpa. The sampling time is T = 5 s. Three MPC algorithms are used to compare the system performance: the typical DMC, the proposed triple-mode MPC (Algorithm 2), and the DMC algorithm with an artificial selection of future target horizon (denoted as DMC-FF in the following). Tuning parameters of the three MPC algorithms are the same as follows: the prediction horizon is P = 300 and the control horizon is M = 20; the weighting matrices are set as: the weight for N E is q 1 = 1 and for p t is q 2 = 100, the weight for inputs is W = 5 * I. In the DMC-FF algorithm, the future trajectory horizon for outputs N E and p t are set to f 1 = 80, f 2 = 120 respectively. In addition, the triple-mode MPC adopts f 1 = 22, f 2 = 26 to design the feedforward based on Algorithm 1. The main purpose of the CCS is to keep the output power responding to load demands quickly when the unit running in either constant or sliding pressure operation. The following two scenarios are considered:

Constant Pressure Operation
When the unit is running in constant pressure operation mode, the main steam pressure is maintained at the rated value regardless of any changes in load demand.
Assuming that the evolution of load demand is given as: the initial load demand is coincident with the system's initial operating point N E = 550 MW; at sample time 200, change the load demand to 600 MW then remain constant; at sample time 700, set the target load to 500 MW and then push it to the 550 MW at time instant 1500. In order to prevent the step load disturbance of the target load from impacting the entire control system of the unit, the setpoint of the output power is formed by limiting the load change rate after the target load is changed [22]. In general, the load change rate is determined by the operating personnel according to the operating conditions of the unit equipment.
Here we set the load change rate as 5 MW/min. Figure 7 presents the dynamics of the controlled outputs N E and p t . The figure shows the responses obtained by the typical DMC, the triple-mode MPC and the DMC-FF algorithms. The control actions implemented by the manipulated inputs µ t and u B that caused the responses in Figure 7, are shown in Figure 8. It can be observed the output power governed by triple-mode MPC can track the load demand more quickly while maintaining steam pressure at the constant value. In contrast, the DMC-FF algorithm uses too much future target information and responses to the load demand far in advance, which leads to large oscillations in the main steam pressure. The main steam pressure controlled by the typical DMC algorithm also oscillates a lot when the load demand changes. In this case, the typical DMC algorithm even outperforms the DMC-FF algorithm due to the improper selection of future target horizon in the DMC-FF algorithm. Table 3 reports the square errors between the output response and the desired setpoint value. In contrast, the proposed triple-mode MPC is expressed explicitly in terms of future information feedforward, and is hence expected to yield improved results for tracking performance.

Sliding Pressure Operation
We then proceed with the running of the unit concerning sliding pressure operation, and the optimal operating curve is plotted in Figure 9. The main steam pressure increases in line with the increase of the load when the load is less than 550 MW, and maintains constant when the load exceeds 550 MW. The load demands used for tracking are the same as before.  Figure 10 presents the dynamics of the controlled outputs N E and p t . The figure shows the responses obtained by the typical DMC, the proposed triple-mode MPC and the DMC-FF algorithms. The control actions implemented by the manipulated inputs µ t and u B that caused the responses in Figure 10, are shown in Figure 11. Analyzing the results in the first panel of Figure 10, one can notice that for the trajectory of the controlled output N E under the proposed triple-mode MPC algorithm is closer to the desired trajectory than which under the typical DMC algorithm. In addition, the trajectory governed by the DMC-FF algorithm can also respond to the load demand more quickly than the typical DMC algorithm. (The curves of triple-mode MPC and DMC-FF almost overlap in the left magnified figure) However, the DMC-FF algorithm responses to the load demand far too early, which leads to the loss of tracking performace as we can see in the right magnified figure.
On the other hand, as shown in the second panel of Figure 10, when the setpoint of the main steam pressure changes with the load demand, the accuracy of tracking deteriorates a lot under either the typical DMC or the DMC-FF algorithms. In contrast, using the triple-mode MPC can drive the main steam pressure p t to track its setpoint without too much performance loss.
A quantitative comparison of square errors for the whole simulation is also shown in Table 4. In the table, the power output and pressure errors of the triple-mode MPC algorithm are dramatically reduced compared to those of the typical DMC and the DMC-FF algorithms. Remark 1. In many practical processes such as power plants [23], wind tunnels [24], there exists future target information that can be used in feedforward design. When the MPC algorithm is applied to real systems, the closed-loop performance is related to the accuracy of the model and the future target information. If the model is consistent with the plant, and the future target information is correct, the control performance obtained is as good as which obtained in the case study.

Conclusions
In this work, we propose a triple-mode model predictive control algorithm using future target information. In order to explicitly take into account the future target information in the MPC optimization, the proposed triple-mode control law includes the future target information feedforward term, and the future trajectory horizon can be optimized off-line to design the feedforward. Furthermore, the extra degrees of freedom are added to append the control law such that the constraints can be met when the MPC is running. Experimental results on the shell fundamental control problem, and the coordinated control problem in the power plant demonstrate that, by incorporating future target information feedforward explicitly and using an optimal future target horizon, the proposed triple-mode MPC algorithm could better hedge against the future target information, mitigate the sluggish response, and improve the tracking performance.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Proof of Theorem 1
Proof of Theorem 1. A deviation model is obtained by computing an error model with respect to the future targets. To do so, consider a desired future trajectory described by R(t + 1|t), and define deviation variablesȲ that satisfy the dynamic modelȲ The zero regulation problem applied to the system in deviation variables finds ∆U M (t) that takes Y(t|t) to zero.