About the Authors:
Zhang Longge
* E-mail: [email protected]
Affiliation: Department of Mathematics and Physics, North China Electric Power University, Baoding, PR China
ORCID http://orcid.org/0000-0003-1575-1862
Yan Yan
Affiliation: Department of Mathematics and Physics, North China Electric Power University, Baoding, PR China
Abstract
In this paper, a new off-line model predictive control strategy is presented for a kind of linear parameter varying system with polytopic uncertainty. A nest of shrinking ellipsoids is constructed by solving linear matrix inequality. By splitting the objective function into two parts, the proposed strategy moves most computations off-line. The on-line computation is only calculating the current control to assure the system shrinking into the smaller ellipsoid. With the proposed formulation, the stability of the closed system is proved, followed with two numerical examples to demonstrate the proposed method’s effectiveness in the end.
Figures
Fig 4
Table 1
Fig 1
Fig 2
Table 2
Table 3
Fig 3
Fig 4
Table 1
Fig 1
Fig 2
Citation: Longge Z, Yan Y (2017) Robust shrinking ellipsoid model predictive control for linear parameter varying system. PLoS ONE 12(6): e0178625. https://doi.org/10.1371/journal.pone.0178625
Editor: Xiaosong Hu, Chongqing University, CHINA
Received: October 16, 2016; Accepted: May 16, 2017; Published: June 2, 2017
Copyright: © 2017 Longge, Yan. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: Support was provided by the National Nature Science Foundation of Hebei Provincial under Grants F2016502025 (ZL); Fundamental Research Funds for the Central Universities under Grant 9161015007 (ZL); National Natural Science Foundation of China under Grants 61403137 (ZL).
Competing interests: The authors have declared that no competing interests exist.
1. Introduction
Model predictive control (MPC), also known as receding or moving horizon control, is an effective control algorithm widely adopted in industry to deal with multivariable constrained control problem. MPC solves the constrained optimization problem at each sampling time and implements only the first element of the optimal control profile [1–6]. In recent ten years, some important branches is extended such as distributed MPC[7, 8], economic MPC[9, 10] and tube based MPC[11, 12] and so on.
Linear parameter varying (LPV) systems becomes a standard formalism in systems and control. It is introduced by Shamma [13] and it is an intermediate step between linear time-invariant (LTI) systems and non-linear plants. LPV systems can approximate many nonlinear systems and the gains can be automatically scheduled with respect to the parameters[14]. As the importance of LPV system, it has been widely investigated. Its common theme is to make the controller parameter dependent so that when the time-varying parameters are measured in real-time, the controller becomes self-scheduling and offers potential performance improvements over a fixed robust controller. Some researchers are concerned with the reducing the conservatism and improving the system’s performance. It is proposed that a state feedback MPC scheme based on a quasi-min-max algorithm [15], the first stage cost can be computed without any uncertainty, thus the first state cost can be determined separately from the rest of parameter changes. A robust dynamic output feedback MPC strategy is designed for a linear fractional representation represented systems[16]. Through the off-line designing robust state observer and on-line robust output feedback MPC controller designing, it is proposed a robust MPC scheme for LPV systems[17]. When the parameters have stochastic nature, via scenario optimization, a stochastic MPC for LPV systems is proposed[18]. In some practical applications, the varying parameters change rates are limited, and a novel algorithm is presented to estimate the future parameter variations which are predicted as a sequence of polytopic families [19, 20]. The designed feedback Robust MPC can improve the control performance. A robust MPC design method is proposed by Cao and Lin[21] to solve the influence of the actuator saturation, which degrades system performance and destroys the system’s stability. This method is improved by placing heavier weighting on the system corresponding to the actual linear feedback law[22].
On the other hand, the requirement of optimality leads to high on-line computation and this limits its application to relatively slow dynamics or small-scale process. To overcome this problem, some authors have proposed the off-line MPC. For example, a series of controllers corresponding to a sequence of nested asymptotically stable invariant ellipsoids is constructed off-line one with another[23]. This result is improved based on the nominal performance and followed with the improvement of the closed loop system’s feasibility and optimality[24].
In the LPV framework, it is assumed that the parameter is measureable or non-measurable. In the former case, the aforementioned researches have little use of the measured parameter vector. To solve this problem, this article aims to provide a formulation of decreasing the computation of the robust MPC for LPV system. The proposed method makes good use the information of the current states and parameter vector, and the future states can be described by the polytopic uncertain. The designed formulation employs the off-line designed shrinking ellipsoids which is first proposed in our previous work[25]. The shrinking ellipsoids have the followed character: at the current time the state , and at the next sampling time the state can be shrunk into the smaller ellipsoid with a controller, i.e. if , then . The on-line controller is calculated with low computation followed with the proof of the closed-loop system’s feasibility and stability.
The rest of this paper is organized as follows. Section 2 gives the problem description. In section 3, the shrinking ellipsoids are designed off-line, followed with the on-line control strategy. In section 4, two examples are presented to illustrate the effectiveness of the proposed strategy. Finally, we conclude the note in section 5.
The following notation will be used. Let Rn be the n-dimensional space of real valued vectors. For a matrix Q and a vector x ∈ Rn, xTQx will be denoted by . The matrix inequality A > B(A ≥ B) means that A and B are square symmetric matrices and A − B is (semi-) positive definite. The measured or actual value of variable x at real time k will be denoted by x(k) (or x(k | k)). x(k + i | k), i ∈ {0, 1, 2, …} is the predicted value of x at a future prediction time k + i predicted at real time k. The identity matrix with proper dimension is denoted by I. * denotes the corresponding transpose of the lower block part of symmetric matrices.
2. Problem description
Consider a linear discrete-time LPV system whose system matrices are affine functions of a parameter vector p(k):(1)where , and denote the state and input respectively. The time-varying parameter vector p(k) = [p1(k), p2(k), ⋯, pL(k) ∈ RL] belongs to a convex polytope , i.e., .
The aim of the research is to find a feedback control law(2)to achieve the followed performance cost(3)subject to state and output constraints(4)(5)where Q, R > 0 are weighting matrices. The MPC for the LPV system is transformed into a convex optimization problem using parameter dependent Lyapunov function[26]. It is less conservative as compared with the result of M.V. Kothare for existing the loosen variables G. However, the convex optimization involves L2 + 3L LMIs (L is the number of convexes polytope) and it requires prohibitive on-line computation.
Lemma 1[26]. Consider the System (1) at the sampling time k with unknown parameters. Given x(k | k), ymax, umax a state feedback control law is obtained by solving the following problem:(6)s.t.
(7)(8)(9)(10)
On the other hand, the speed of the closed-loop response can be influenced by specifying a minimum decay rate on the state x(||x(k)|| ≤ ρk ||x(0)||, 0 < ρ < 1) as follows:(11)for any [A(k + i), B(k + i)] ∈ Ω.
Lemma 2 [3]: For System (1), if it has a minimum decay rate of Eq (11), the following LMI must be satisfied:(12)
3. Shrinking ellipsoidal MPC for LPV system
In this section, it is assumed that both the parameter vector p(k) and the state x(k) are available in real-time. At sampling time k the system parameter vector p(k) is known exactly but unknown in the future. The designed strategy includes three stages: first, the method of seeking the minimum decay rate is presented for a fixed state; then a consequence of nested ellipsoids is constructed off-line based on the iterative method; finally the MPC algorithm is formulated. In the proposed method, most of the computation is moved off-line and only the low computation of calculating the input is left online.
Algorithm 3 (seeking the minimum decay rate). Consider the uncertain System (1) with Constraints (4) and (5). Let ρ = 1.
1. Investigate the feasibility of the following problem:
Problem 1: subject to Eqs (7) ~ (10) and (12).
2. If problem 1 is feasible, let ρ ≔ ρ − 0.01, go to step 1. Otherwise, return ρ.
Remark 4: The system is supposed to converge to the smaller invariant ellipsoid. The computed minimum decay rate is designed off-line and it is used to assure the stability of the closed-loop system.
Lemma 5 [25]. Suppose the matrices satisfy the followed condition 0 < ρ2Qi < Qi+1 < Qi(0 < ρ < 1), define and , then .
Algorithm 6 (off-line robust MPC). Consider the uncertain System (1) with Constraints (4) and (5). Given an initial feasible state x1, compute a sequence of minimizers as follows off-line. Let i ≔ 1.
1. Compute the minimum decay rate ρi at xi using algorithm 3.
2. Compute the optimizer γi,j, Qi,j, Xi,j, Yi,j, Zi,j at xi using lemma 1 with an additional constraint (ignored at time i = 1, j = 1, 2⋯, L), store γi,j, γi, Qi,j and , where
3. if i < N, choose a state xi+1 satisfying . Let i ≔ i + 1, go to step 1.
The state of the system is controlled into the smaller ellipsoid on-line. Suppose the state satisfying , we select the proper controller to drive the system into , that is . Consider the following objective function which is split into two parts:(13)
As the current state and parameter are unknown, so the first stage cost can be computed without uncertainty. This is the reason that the first control input can be separately from the rest of the future control law. Our on-line control strategy is based on the following two facts: the current state is in the ellipsoid , but not in , but in the following sampling time k + 1, the state is in , that is ; x(k + 1 | k) and p(k + 1 | k) cannot be measured on-line, so in the future sampling time, the system is described by polytope uncertainty, and the future control law can be designed off-line using the Algorithm 6. So in this case , then the objective function satisfied(14)
:(15)
The parameter vector p(k) and the current state x(k | k) are known, and only the input is free variable. The on-line work is only to compute the controller at the current sampling time(16)satisfied Eq (15) and(17)
The Conditions (15) and (17) are equivalent to the followed LMIs respectively(18)(19)
Theorem 7. The above optimization problem with the control law given by(20)can be solved by the following semi-definite programming(21)subject to Eqs (18) and (19).
Proof: Minimization of x(k | k)T Qx(k | k) + u(k | k)T Ru(k | k) + γi+1 is equivalent tosubject to x(k | k)T Qx(k | k) + u(k | k)T Ru(k | k) + γi+1 ≤ υ, using Schur complements, it is equivalent towhich proves Eq (19);
Substituting Eq (1) into Eq (15), one can get
Using Schur complements it can be expressed as
Remark 5: As the current parameter vector p(k) is known, there is no need to transform Eq (18) into the LMIs based on the convex of polytope.
The designed controller Eq (20) is composed of two parts: u(k | k) is the controller at the current sampling time and it can drive the state into the inner ellipsoid; is the off-line designed controller will be acted on the system in the following sampling time, and it can shrink the state into the smaller ellipsoid. The following theorem proves the closed-loop system’s stability.
Theorem 8. Given a dynamic System (1) and an initial state x(0) satisfying , the controller Eq (20) robustly asymptotically stabilizes the closed-loop system.
Proof: For the off-line minimization at xi, i = 2, …,N, the additional constraint is equivalent to . From lemma 5, it is known that the left part of the inequality implies that the constructed asymptotically stable invariant ellipsoid is inside , i.e. . So it is guaranteed that is monotonic decreasing with respect to the index i. The right part of the inequality can assure the algorithm’s feasibility, i.e., there exists in the feasible regions. From theorem 7, the control law Eq (20) is guaranteed to drive the state into the next nested ellipsoid. Lastly, the smallest is guaranteed to keep the state within and converge it to the origin.
Remark 10. As we all know, the fastest interior point algorithms show O(RS3) growth in computation where R is the total row size of the LMI system and S is the total number of scalar decision variables[23, 24, 27]. For LMI optimization Problem (6) the parameter S is given by and R is given by L(7nx + 2nu + 1) while for Problem (21) S is given by 2 and R is given by 1 + nx + nu. So the proposed strategy can reduce the on-line computation dramatically.
4. Numerical example
In this section, we present two examples to illustrate the effectiveness of the designed LPV system’s robust shrinking ellipsoid model predictive controller. The simulation is carried on Lenovo computer, and its processor is Intel® core™ [email protected].
Example 1. Consider the following LPV system given by[17](22)Where , , . The initial states of the System (22) is assumed as and the uncertain parameter α(k) belong to the following regions: α(k) ∈ [1, 5]. Then , where , , . We select some points on the x -axis, and get their minimum decay rate using the Algorithm 3 shown in Table 1.
[Figure omitted. See PDF.]
Table 1. Selected points and corresponding minimum decay of example 1.
https://doi.org/10.1371/journal.pone.0178625.t001
It is shown that the minimum decay rate is smaller than 0.8 from Table 1. In the simulation, the decay rate is set to be 0.8 and the nested ellipsoids is shown in Fig 1. The states of the controlled closed-loop system using the designed method (RSEMPC method) compared with the method in paper[26](PDLF method) are shown in Fig 2. The runtime of RSEMPC method is 0.001062 second while the PDLF method is 6.410566 second.
[Figure omitted. See PDF.]
Fig 1. Nested ellipsoids of example 1.
https://doi.org/10.1371/journal.pone.0178625.g001
[Figure omitted. See PDF.]
Fig 2. States of the closed-loop system of example 1.
https://doi.org/10.1371/journal.pone.0178625.g002
It is shown that the computation time of the RSEMPC method is much smaller than the PDLF method and most of the computation is left offline.
Example 2. Consider the following uncertain nonlinear model for non-isothermal control of a continuously stirred tank reactor (CSTR) where the exothermic reaction A → B takes place. A cooling coil is used to remove heat that is released in the exothermic reaction. The reaction rate constant k0 and the heat of reaction ΔHr×n are considered to be the uncertain parameters. The linearized model based on the component balance and the energy on the component balance is given as follows(23)where CA denotes the concentration of A in the reactor, T denotes the reactor temperature and TC denotes the temperature of coolant stream. The rate of reaction is first order with respect to component A. The relevant constants for the CSTR dynamic model are illustrated in Table 2.
[Figure omitted. See PDF.]
Table 2. Relevant constants for the CSTR dynamic model.
https://doi.org/10.1371/journal.pone.0178625.t002
Let , , and , where the subscript eq is used to denote the corresponding variable at equilibrium condition. The discrete-time model is obtained by discretizing using Euler first-order approximation with a sampling time of 0.15 min (F = 1 m3 / min, V = 1 m3, k0 = 109 ~ 5 × 109 min−1, E / R= 8330.1K, −ΔHr×n = 107 ~ 5 × 107 cal / kmol, ρ = 106 g / m3, UA = 5.34 × 106 cal / K, Cp = 1 cal / (gK), Ts = 394K and CAS = 0.265 kmol / m3)(24)where , , 1 ≤ α(k) = k0 / 109 ≤ 5 and 1 ≤ β(k) = −ΔHr×n / 107 ≤ 5. The polytopic Ω = Co{A1, A2, A3, A4}, , , , . It is the same as example 1, we can find the minimum decay rate is smaller than 0.8. In the simulation, the decay rate is set to be 0.8 and the nested ellipsoids constructed off-line by the prosed algorithm is shown in Fig 3. If the state of the system is in one ellipsoid at the present time k, then at the next sampling time k + 1, the state must be in the adjacent interior ellipsoid, so we call them as the shrinking ellipsoids. Table 3 shows the overall on-line numerical burdens compared with PDLF method. It is shown that the designed algorithm requires smaller on-line computation. Fig 4 shows the closed-loop responses of the system.
[Figure omitted. See PDF.]
Table 3. On-line numerical burdens in example 1.
https://doi.org/10.1371/journal.pone.0178625.t003
[Figure omitted. See PDF.]
Fig 3. Nested ellipsoids of example 2.
https://doi.org/10.1371/journal.pone.0178625.g003
[Figure omitted. See PDF.]
Fig 4. States of the closed-loop system of example 2.
https://doi.org/10.1371/journal.pone.0178625.g004
5. Conclusions
In this paper, a novel off-line MPC synthesis approach for a LPV system is presented. A sequence of state feedback gains corresponding to the sequences of nested ellipsoids is pre-computed. The on-line computation is only to calculate the input to control the states into the inner ellipsoids and it only left the current controller as a free variable, so most of the computation is moved off-line. The effectiveness of the propose method is illustrated by two simulation examples.
Supporting information
[Figure omitted. See PDF.]
S1 Fig. States and ellipsoids of example 1.
https://doi.org/10.1371/journal.pone.0178625.s001
(RAR)
S2 Fig. States and ellipsoids of example 2.
https://doi.org/10.1371/journal.pone.0178625.s002
(RAR)
Acknowledgments
This work was supported by the National Nature Science Foundation of Hebei Provincial under Grants F2016502025, Fundamental Research Funds for the Central Universities under Grant 9161015007, National Natural Science Foundation of China under Grants 61403137.
Author Contributions
1. Conceptualization: ZL.
2. Formal analysis: ZL.
3. Funding acquisition: ZL.
4. Methodology: ZL.
5. Software: YY.
6. Validation: YY.
7. Visualization: ZL.
8. Writing – original draft: ZL.
9. Writing – review & editing: YY.
Citation: Longge Z, Yan Y (2017) Robust shrinking ellipsoid model predictive control for linear parameter varying system. PLoS ONE 12(6): e0178625. https://doi.org/10.1371/journal.pone.0178625
1. Clarke DW, Mohtadi C, Tuffs P. Generalized predictive control—Part I. The basic algorithm. Automatica. 1987;23(2):137–48.
2. Garcia CE, Prett DM, Morari M. Model predictive control: theory and practice—a survey. Automatica. 1989;25(3):335–48.
3. Kothare MV, Balakrishnan V, Morari M. Robust constrained model predictive control using linear matrix inequalities. Automatica. 1996;32(10):1361–79.
4. Liu X, Feng S, Ma M. Robust MPC for the constrained system with polytopic uncertainty. International Journal of Systems Science. 2012;43(2):248–58.
5. Zou C, Manzie C, Nešić D. A Framework for Simplification of PDE-Based Lithium-Ion Battery Models. IEEE Transactions on Control Systems Technology. 2016;24(5):1594–609.
6. Zou C, Manzie C, Anwar S. Control-Oriented Modeling of A Lithium-Ion Battery for Fast Charging. IFAC Proceedings Volumes. 2014;47(3):3912–7. http://dx.doi.org/10.3182/20140824-6-ZA-1003.00829.
7. Camponogara E, Jia D, Krogh BH, Talukdar S. Distributed model predictive control. IEEE Control Systems. 2002;22(1):44–52.
8. Venkat AN, Hiskens IA, Rawlings JB, Wright SJ. Distributed MPC strategies with application to power system automatic generation control. IEEE transactions on control systems technology. 2008;16(6):1192–206.
9. Angeli D, Amrit R, Rawlings JB. On average performance and stability of economic model predictive control. IEEE transactions on automatic control. 2012;57(7):1615–26.
10. Diehl M, Amrit R, Rawlings JB. A Lyapunov function for economic optimizing model predictive control. IEEE Transactions on Automatic Control. 2011;56(3):703–7.
11. Limon D, Alvarado I, Alamo T, Camacho E. Robust tube-based MPC for tracking of constrained linear systems with additive disturbances. Journal of Process Control. 2010;20(3):248–60.
12. Mayne DQ, Kerrigan EC, Van Wyk E, Falugi P. Tube-based robust nonlinear model predictive control. International Journal of Robust and Nonlinear Control. 2011;21(11):1341–53.
13. Shamma JS. Analysis and design of gain scheduled control systems: Massachusetts Institute of Technology; 1988.
14. Chandra KPB, Alwi H, Edwards C. Fault Detection in Uncertain LPV Systems with Imperfect Scheduling Parameter Using Sliding Mode Observers. European Journal of Control. 2016.
15. Lu Y, Arkun Y. Quasi-min-max MPC algorithms for LPV systems. Automatica. 2000;36(4):527–40.
16. Lee S-M, Park JH. Output feedback model predictive control for LPV systems using parameter-dependent Lyapunov function. Applied Mathematics and Computation. 2007;190(1):671–6.
17. Park J-H, Kim T-H, Sugie T. Output feedback model predictive control for LPV systems based on quasi-min–max algorithm. Automatica. 2011;47(9):2052–8.
18. Calafiore GC, Fagiano L. Stochastic model predictive control of LPV systems via scenario optimization. Automatica. 2013;49(6):1861–6.
19. Casavola A, Famularo D, Franze G. A feedback min-max MPC algorithm for LPV systems subject to bounded rates of change of parameters. Automatic Control, IEEE Transactions on. 2002;47(7):1147–53.
20. Li D, Xi Y. The feedback robust MPC for LPV systems with bounded rates of parameter changes. Automatic Control, IEEE Transactions on. 2010;55(2):503–7.
21. Cao Y-Y, Lin Z. Min–max MPC algorithm for LPV systems subject to input saturation. IEE Proceedings-Control Theory and Applications. 2005;152(3):266–72.
22. Huang H, Li D, Lin Z, Xi Y. An improved robust model predictive control design in the presence of actuator saturation. Automatica. 2011;47(4):861–4.
23. Wan Z, Kothare MV. An efficient off-line formulation of robust model predictive control using linear matrix inequalities. Automatica. 2003;39(5):837–46.
24. Ding B, Xi Y, Cychowski MT, O’Mahony T. Improving off-line approach to robust MPC based-on nominal performance cost. Automatica. 2007;43(1):158–63.
25. Zhang L, editor Automatic Offline Formulation of Robust Model Predictive Control Based on Linear Matrix Inequalities Method. Abstract and Applied Analysis; 2013: Hindawi Publishing Corporation.
26. Wada N, Saito K, Saeki M. Model Predictive Control for Linear Parameter Varying Systems Using Parameter Dependent Lyapunov Function. Circuits and Systems II: Express Briefs, IEEE Transactions on. 2006;53(12):1446–50.
27. Gahinet P, Nemirovski A, Laub A, Chilali M. LMI Control Toolbox-for Use with Matlab, Natick, MA: The MATH Works. Inc; 1995.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2017 Longge, Yan. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
In this paper, a new off-line model predictive control strategy is presented for a kind of linear parameter varying system with polytopic uncertainty. A nest of shrinking ellipsoids is constructed by solving linear matrix inequality. By splitting the objective function into two parts, the proposed strategy moves most computations off-line. The on-line computation is only calculating the current control to assure the system shrinking into the smaller ellipsoid. With the proposed formulation, the stability of the closed system is proved, followed with two numerical examples to demonstrate the proposed method’s effectiveness in the end.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer