Journal of Control Science and Engineering

Journal of Control Science and Engineering / 2012 / Article
Special Issue

Model Predictive Control

View this Special Issue

Research Article | Open Access

Volume 2012 |Article ID 630645 |

He De-Feng, Yu Li, "Pole Placement-Based NMPC of Hammerstein Systems and Its Application to Grade Transition Control of Polypropylene", Journal of Control Science and Engineering, vol. 2012, Article ID 630645, 6 pages, 2012.

Pole Placement-Based NMPC of Hammerstein Systems and Its Application to Grade Transition Control of Polypropylene

Academic Editor: Baocang Ding
Received01 Jul 2011
Accepted11 Sep 2011
Published24 Nov 2011


This paper presents a new nonlinear model predictive control (MPC) algorithm for Hammerstein systems subject to constraints on the state, input, and intermediate variable. Taking no account of constraints, a desired linear controller of the intermediate variable is obtained by applying pole placement to the linear subsystem. Then, actual control actions are determined in consideration of constraints by online solving a finite horizon optimal control problem, where only the first control is calculated and others are approximated to reduce the computational demand. Moreover, the asymptotic stability can be guaranteed in certain condition. Finally, the simulation example of the grade transition control of industrial polypropylene plants is used to demonstrate the effectiveness of the results proposed here.

1. Introduction

Hammerstein systems consist of the cascade connection of a static (memoryless) nonlinear function followed from a linear dynamic system. Under certain assumptions (such as fading memory), the Hammerstein model could approximately represent nonlinear dynamics of real systems and has been successfully applied to many kinds of industrial processes such as pH neutralization [1], distillation [2], and polymerization transition [3]. In recent years, constrained control of Hammerstein systems has become one of the most needed and yet very difficult tasks for the process industry.

Model predictive control (MPC) is an effective control algorithm for handling constrained control problems, and various MPC algorithms have been proposed for control of Hammerstein systems with constraints [4ā€“11]. Making use of the geometry structure of Hammerstein model, [4ā€“9] developed two-step MPC schemes for Hammerstein systems with input constraints, where the intermediate variable was firstly obtained by linear MPC and then the actual control was calculated by solving nonlinear algebraic equations, desaturation, and so forth. Although two-step MPC algorithms possess more light computational burden than that where the nonlinearities were incorporated into optimal control problems directly [10, 11], solving nonlinear algebraic equations will inevitably have error and the restricted control is generally different from the desired one for constrained Hammerstein systems [6]. Therefore, the performance and stability properties may be deteriorated in the presence of constraints. Moreover, to the best of authorsā€™ knowledge, the conversional two-step MPC schemes usually have limited ability to dealing with input constraints. So, [12] presented a new two-step MPC scheme for Hammerstein systems with overall constraints, where the actual controller was determined by solving a finite horizon optimization problem that was described as tracking quadratic, optimal trajectories of the intermediate variable. However, the online computational demand of optimization problem of the scheme is still an open question.

In this paper, we consider the problem of the grade transition control of industrial polypropylene plants and propose an efficient NMPC scheme for Hammerstein systems with overall constraints based on the pole placement method [13]. We firstly obtained desired trajectories of Hammerstein systems by using pole placement to the linear subsystem and then calculated actual control actions by solution of finite horizon optimal control problems. In order to reduce the online computational burden of the optimization problem, only the first control action is calculated and the remaining is approximated by solving nonlinear algebraic equations over predictive horizon. Thus, the performance and stability properties can be guaranteed in the face of overall constraints, with moderate increment of the computational demand. Finally, an example of the grade transition control of propylene plants is used to demonstrate the effectiveness of the proposed algorithm.

2. Pole Placement-Based NMPC Algorithm

Consider the discrete-time Hammerstein systems described byš‘„(š‘”+1)=š“š‘„(š‘”)+šµš‘£(š‘”),š‘£(š‘”)=š‘”(š‘¢(š‘”),š‘”),š‘”=0,1,ā€¦,(1) where š‘„āˆˆš‘…š‘›, š‘£āˆˆš‘…š‘š, and š‘¢āˆˆš‘…š‘Ÿ are the state, intermediate variable and control input, respectively. Vector field š‘” represents the nonlinear relationship between the input and intermediate variable, satisfying š‘”(0,š‘”)=0. It is assumed that (š“,šµ) is stabilizable and the state is available for feedback control. Without loss of generality, we also assume that the origin is the equilibrium of system (1).

Given constraints of system (1) as follows: š‘„(š‘”)āˆˆš‘€š‘„=ī€½š‘„āˆˆš‘…š‘›āˆ¶š‘„LBā‰¤š‘„ā‰¤š‘„UBī€¾,š‘”=0,1,ā€¦,š‘£(š‘”)āˆˆš‘€š‘£=ī€½š‘£āˆˆš‘…š‘šāˆ¶š‘£LBā‰¤š‘£ā‰¤š‘£UBī€¾š‘¢,š‘”=0,1,ā€¦,(š‘”)āˆˆš‘€š‘¢=ī€½š‘¢āˆˆš‘…š‘Ÿāˆ¶š‘¢LBā‰¤š‘¢ā‰¤š‘¢UBī€¾,š‘”=0,1,ā€¦,(2) where sets š‘€š‘„āŠ†š‘…š‘›, š‘€š‘£āŠ†š‘…š‘š, and š‘€š‘¢āŠ†š‘…š‘Ÿ are constraints on the state, intermediate variable, and input, respectively. The goal of this paper is to design an efficient NMPC algorithm for the Hammerstein system (1) with the constraint (2).

2.1. Design of NMPC Algorithm

Consider š‘› desired poles, denoted by {šœ†1,ā€¦,šœ†š‘›}, of the linear subsystem of (1), which represent some desired performance of system (1). Associated with this set of desired poles is a linear controller of the intermediate variableš‘£(š‘”)=āˆ’š¾š‘„(š‘”)(3) which is generated via Lyapunovā€™s method without any constraints [13]. Since the controller (3) is unable to be implemented by real control action, here it is labeled as š‘£(š‘”)š‘‘=āˆ’š¾š‘„(š‘”)š‘‘ with associated closed-loop states š‘„(š‘”)š‘‘.

Define a finite horizon objective function š½(š‘”) as follows:š½(š‘”)=š‘‡š‘āˆ’1ī“š‘–=0ī‚†ī€ŗš‘„(š‘”+š‘–+1āˆ£š‘”)āˆ’š‘„(š‘”+š‘–+1)š‘‘ī€»Tī€ŗƗš‘„š‘„(š‘”+š‘–+1āˆ£š‘”)āˆ’š‘„(š‘”+š‘–+1)š‘‘ī€»+š‘¢(š‘”+š‘–āˆ£š‘”)Tī‚‡,š‘…š‘¢(š‘”+š‘–āˆ£š‘”)(4) where the vector ā€¢(š‘”+š‘–āˆ£š‘”) is š‘–-step ahead prediction from time instant š‘”, integer 0<š‘‡š‘<āˆž is the predictive horizon, and matrices 0ā‰¤š‘„ and 0<š‘… are the weighting matrices of the state and input, respectively. Then a new two-step MPC algorithm of Hammerstein systems with overall constraints is presented below.

Algorithm 1 (modified two-step MPC, MTMPC). This algorithm comprises the following steps:Step 1. Given š‘› desired poles {šœ†1,ā€¦,šœ†š‘›} and off-line compute the gain matrix š¾ in (3) such that the eigenvalues of š“āˆ’šµš¾ are those specified in vector {šœ†1,ā€¦,šœ†š‘›}.Step 2. Set š‘‡š‘, š‘„, š‘… in the objective function in (4).Step 3. With š‘£(ā‹…)š‘‘=āˆ’š¾š‘„(ā‹…)š‘‘ and the current state š‘„(š‘”), solve on-line the nonlinear algebra equations (5) without constraints to obtain š‘‡š‘ desired control actions š‘¢(š‘”+š‘–)š‘‘š‘£(š‘”+š‘–)š‘‘ī€·āˆ’š‘”š‘¢(š‘”+š‘–)š‘‘ī€øš‘„,š‘”+š‘–=0,(š‘”+š‘–+1)š‘‘=(š“āˆ’šµš¾)š‘„(š‘”+š‘–)š‘‘,š‘„(š‘”)š‘‘=š‘„(š‘”),š‘–=0,1,ā€¦,š‘‡š‘āˆ’1.(5)Step 4. With the current state š‘„(š‘”), determine on-line the actual control action š‘¢(š‘”) as follows: š‘¢(š‘”)=argminš‘¢(š‘”āˆ£š‘”)š½(š‘”),s.t.š‘„(š‘”+š‘–+1āˆ£š‘”)=š“š‘„(š‘”+š‘–āˆ£š‘”)+šµš‘”(š‘¢(š‘”+š‘–āˆ£š‘”),š‘”+š‘–),š‘„(š‘”+š‘–+1)š‘‘=(š“āˆ’šµš¾)š‘„(š‘”+š‘–)š‘‘,š‘¢(š‘”+š‘–āˆ£š‘”)=š‘¢(š‘”+š‘–)š‘‘,š‘–=1,ā€¦,š‘‡š‘āˆ’1,š‘„(š‘”+1āˆ£š‘”)āˆˆš‘€š‘„,š‘”(š‘¢(š‘”āˆ£š‘”),š‘”)āˆˆš‘€š‘£,š‘¢(š‘”āˆ£š‘”)āˆˆš‘€š‘¢,š‘„(š‘”āˆ£š‘”)=š‘„(š‘”)š‘‘=š‘„(š‘”),(6) where š½(š‘”) is defined by (4).Step 5. Implement š‘¢(š‘”), set š‘”=š‘”+1, and go back Step 3.

Remark 2. For the algorithm here, only the first control action (i.e., š‘¢(š‘”āˆ£š‘”)) is optimized subject to constraints at each time. So, the online optimization problem has only š‘Ÿ decision variables versus š‘‡š‘š‘Ÿ decision variables for general nonlinear MPC, which results in a significant reduction of the online computational demand since the computational demand, in general, grows exponentially with the number of decision variables [14]. Thus, the predictive horizon š‘‡š‘ should be chosen reasonably large to ensure adequate performance. On the other hand, since only the first control is implemented in MPC scheme, it is expected that only enforcing constraints for the first prediction should give excellent results. In addition, the closed-loop performance of MPC systems clearly depends on how š‘› desired poles {šœ†1,ā€¦,šœ†š‘›} are arranged. The further work will focus on how to determine {šœ†1,ā€¦,šœ†š‘›} optimally.

It should be pointed out that the feasibility of optimization problem (6) may not be guaranteed theoretically with respect to arbitrary constraints in (2). However, it is usually ensured in practice via adjusting the intermediate variable constraint. Next, we show that the closed-loop MPC system (1) is asymptotically stable, provided that the feasibility of (6) is hold.

2.2. Analysis of Stability

In terms of Algorithm MTMPC, the MPC law of constrained Hammerstein system (1), (2) is defined as š‘¢(š‘”)MPC=š‘¢(š‘”āˆ£š‘”),š‘”=0,1,ā€¦(7) with the closed-loop systemī€·š‘„(š‘”+1)=š“š‘„(š‘”)+šµš‘”š‘¢(š‘”)MPCī€ø=ī€ŗī€·š‘¢,š‘”(š“āˆ’šµš¾)š‘„(š‘”)+šµš¾š‘„(š‘”)+š‘”(š‘”)MPC.,š‘”ī€øī€»(8)

Definition 3. Set š‘€āŠ†š‘…š‘› is called an attractive region of system (8) if for all š‘„(š‘”)āˆˆš‘€, the system trajectories š‘„(š‘ ;š‘„(š‘”))ā†’0 when š‘ ā†’+āˆž and satisfies š‘„(š‘”)āˆˆš‘€āŸ¹š‘„(š‘”+1)āˆˆš‘€,āˆ€š‘¢(š‘”)MPCāˆˆš‘€š‘¢š‘£,š‘”=0,1,ā€¦,(9) where subset š‘€š‘¢š‘£={š‘¢āˆˆš‘…š‘Ÿāˆ¶š‘”(š‘¢,ā‹…)āˆˆš‘€š‘£}āŠ†š‘€š‘¢.

Let a set of points {šœ†1,ā€¦,šœ†š‘›} be the desired poles of the linear subsystem š‘„(š‘”+1)=(š“āˆ’šµš¾)š‘„(š‘”). Then from inverse Lyapunovā€™s theorem [13], there is systematic matrices š‘ƒ>0 and š‘„>0 such that(š“āˆ’šµš¾)Tš‘ƒ(š“āˆ’šµš¾)āˆ’š‘ƒ<āˆ’š‘„.(10) In following, we present the result on our proposed algorithm.

Theorem 4. Consider the system (1), (2), there exists a nonempty region š‘€ such that the closed-loop system (8) is asymptotically stable if the inequality (11) is satisfied ī€·š‘‰(š‘„(š‘”+1))ā‰¤š‘‰š‘„(š‘”+1)š‘‘ī€ø,āˆ€š‘„(š‘”)āˆˆš‘€,š‘”=0,1,ā€¦,(11) where function š‘‰(š‘„)=š‘„Tš‘ƒš‘„ and š‘„(š‘”+1)š‘‘=(š“āˆ’šµš¾)š‘„(š‘”)š‘‘ with š‘„(š‘”)š‘‘=š‘„(š‘”). Moreover, the set š‘€ is an attractive region of system (8).

Proof. Consider the system (1), (2), define a level set š‘€ of function š‘‰(š‘„) as š‘€={š‘„āˆˆš‘…š‘›āˆ¶š‘‰(š‘„)ā‰¤š‘,š‘>0}āŠ†š‘€š‘„(12) such that š‘¢(š‘”)MPCāˆˆš‘€š‘¢š‘£, for all š‘„(š‘”)āˆˆš‘€, š‘”=0,1,ā€¦. Note that the set š‘€ is always nonempty since the constraint in (2) includes the origin as the internal point.
Let function š‘‰(š‘„) be a Lyapunov function candidate of system (8). Then, for any š‘„(š‘”)āˆˆš‘€, we have š‘‰(š‘„(š‘”+1))āˆ’š‘‰(š‘„(š‘”))=š‘„(š‘”+1)Tš‘ƒš‘„(š‘”+1)āˆ’š‘„(š‘”)Tš‘ƒš‘„(š‘”)=š‘„(š‘”)T(š“āˆ’šµš¾)Tš‘ƒ(š“āˆ’šµš¾)š‘„(š‘”)Tāˆ’š‘„(š‘”)Tš‘ƒš‘„(š‘”)+2š‘„(š‘”)T(š“āˆ’šµš¾)Tš‘ƒšµĪ”+Ī”TšµTš‘ƒšµĪ”=š‘„(š‘”)Tī€ŗ(š“āˆ’šµš¾)Tī€»š‘ƒ(š“āˆ’šµš¾)āˆ’š‘ƒš‘„(š‘”)+2š‘„(š‘”)T(š“āˆ’šµš¾)Tš‘ƒšµĪ”+Ī”TšµTš‘ƒšµĪ”,(13) where Ī”=š‘”(š‘¢(š‘”)MPC,š‘”)+š¾š‘„(š‘”)=š‘”(š‘¢(š‘”)MPC,š‘”)āˆ’š‘£(š‘”)š‘‘. Substituting (10) into (13) yields š‘‰(š‘„(š‘”+1))āˆ’š‘‰(š‘„(š‘”))=āˆ’š‘„(š‘”)T[]š‘„š‘„(š‘”)+2(š“āˆ’šµš¾)š‘„(š‘”)Tš‘ƒšµĪ”+Ī”TšµTš‘ƒšµĪ”=āˆ’š‘„(š‘”)Tš‘„š‘„(š‘”)+2š‘„(š‘”)Tš“Tš‘ƒšµĪ”+2(š‘”āˆ’Ī”)TšµTš‘ƒšµĪ”+Ī”TšµTš‘ƒšµĪ”=āˆ’š‘„(š‘”)Tš‘„š‘„(š‘”)+2š‘„(š‘”+1)Tš‘ƒšµĪ”āˆ’Ī”TšµTš‘ƒšµĪ”=āˆ’š‘„(š‘”)Tš‘„š‘„(š‘”)āˆ’š‘„(š‘”+1)š‘‘Tš‘ƒš‘„(š‘”+1)š‘‘+š‘„(š‘”+1)Tš‘ƒš‘„(š‘”+1)=āˆ’š‘„(š‘”)Tī€·š‘„š‘„(š‘”)āˆ’š‘‰š‘„(š‘”+1)š‘‘ī€ø+š‘‰(š‘„(š‘”+1)).(14) From the inequality (11), it is straight forward to obtain that š‘‰(š‘„(š‘”+1))āˆ’š‘‰(š‘„(š‘”))ā‰¤āˆ’š‘„(š‘”)Tš‘„š‘„(š‘”),āˆ€š‘„(š‘”)āˆˆš‘€.(15) Hence, function š‘‰(š‘„) is a Lyapunov function of system (8), and the system is asymptotically stable in š‘Ÿegion š‘€. Furthermore, from the definition of attractive regions and set š‘€ in (12), we derive that š‘€ is an attractive region of system (8).
This establishes the theorem.

Remark 5. From the proof of the above theorem, we know that the stability property of closed-loop system (8) is guaranteed by the linear controller of the intermediate variable, provided that the feasibility of optimization problem (6) and condition (11) holds. This suggests that the stability and optimality can be separated to some extent. Therefore, the design parameter of the objective function in (4) can be tuned freely to achieve more satisfactory performance, regardless the stability of the closed-loop system. In addition, it is pointed out that the attractive region š‘€ obtained here is merely a theoretical concept and the computation of š‘€ is not an easy task.

Remark 6. Instead of conventional two-step MPC algorithms, where the actual control was determined by solving nonlinear algebraic equations based on desired controllers of the intermediate variable resulted from linear MPC, actual control actions here are derived by solving a finite horizon optimal control problem (i.e., nonlinear MPC), which is defined as tracking the desired state trajectories driven by pole placement state feedback controller. Thus, overall constraints of Hammerstein systems are taken into account when designing the predictive controller, and then performance and stability properties are guaranteed. Meanwhile, the ability to handling constraints is enhanced in the algorithm proposed here. Note that the linear controller of intermediate variables in (3) can be any one that stabilizes the linear subsystems of Hammerstein models with certain desired performance.

3. Simulation Example

Taking an example of polypropylene (PP) grade transition control, we illustrate the effectiveness of the proposed algorithm.

In propylene polymerization industry, melt index (MI, g/10ā€‰min) is usually used to identify grade indices of PP homopolymer and both of MI and concentration of ethylene (šøš‘”, %) in PP denote grade indices of total phase PP [12, 15ā€“18]. This implies that the nature of grade transition process is to operate the response of the variables MI and šøš‘” in polymer. Then according to the propylene polymerization mechanism, we have the following mathematical model which formulates the dynamics of the grade transition in polypropylene process [12]:āŽ”āŽ¢āŽ¢āŽ£Ģ‡š‘„1Ģ‡š‘„2āŽ¤āŽ„āŽ„āŽ¦=āŽ”āŽ¢āŽ¢āŽ£āˆ’1šœ010āˆ’šœāŽ¤āŽ„āŽ„āŽ¦āŽ”āŽ¢āŽ¢āŽ£š‘„1š‘„2āŽ¤āŽ„āŽ„āŽ¦+āŽ”āŽ¢āŽ¢āŽ£1šœ001šœāŽ¤āŽ„āŽ„āŽ¦āŽ”āŽ¢āŽ¢āŽ£š‘”1š‘”(š‘¢,š‘”)2āŽ¤āŽ„āŽ„āŽ¦,āŽ”āŽ¢āŽ¢āŽ£š‘”(š‘¢,š‘”)1š‘”(š‘¢,š‘”)2āŽ¤āŽ„āŽ„āŽ¦=āŽ”āŽ¢āŽ¢āŽ¢āŽ£š‘˜(š‘¢,š‘”)1+š‘˜2š‘¢1+š‘˜3ī€·š‘˜log4+š‘˜5š‘¢2+š‘˜6š‘¢3ī€øš‘˜7š‘¢2+1š‘˜8/š‘¢3+š‘˜9š‘¢3+š‘˜10āŽ¤āŽ„āŽ„āŽ„āŽ¦,āŽ”āŽ¢āŽ¢āŽ£āŽ¤āŽ„āŽ„āŽ¦š‘¦(š‘”)=1001š‘„(š‘”),š‘”ā‰„0,(16) where states š‘„1 and š‘„2 are the cumulative MI (MIš‘, g/10ā€‰min) and šøš‘” (šøš‘”š‘, %) of PP, respectively; intermediate variables š‘”1 and š‘”2 denote the instantaneous MI (MIš‘–, g/10ā€‰min) and šøš‘” (šøš‘”š‘–, %) of PP, respectively; input, š‘¢1, š‘¢2, š‘¢3 denote the reaction temperature š‘‡(K), concentration ratio of hydrogen to propylene š¶H2/š¶m (%), and concentration ratio of ethylene to propylene š¶m2/š¶m (%), respectively, and parameter šœ is the average polymer residence time (šœ = 2ā€‰h). The parameters of model š‘˜š‘– (š‘–=1,ā€¦,10) can be identified on-line by industrial data [15].

The grade transition sequence considered here is Aā†’Bā†’C, where grades A and B are PP homopolymer and C is total phase PP. It should be pointed out that the feed flow rate of ethylene is always zero (i.e., [š¶m2/š¶m]=0) in production of PP homopolymer grades A and B. The grade transition process is defined by the set-point of states and input/output constraints as shown in Table 1, where superscripts š‘, š‘–, and Ī” denote accumulation, instantaneous value, and increment constraints, respectively. Each grade transition is optimized over 10ā€‰h horizon with a sampling interval of 0.5ā€‰h.

GradeMI (g/10ā€‰min) šø š‘” (%) š‘‡ (K) š¶ H 2 (%) š¶ m 2 / š¶ m (%)

A ā†’ B [ 2 . 4 , 4 2 . 0 ] š‘ ā€” [ 3 4 1 . 1 5 , 3 4 5 . 1 5 ] [ 0 . 0 2 , 0 . 3 5 ] ā€”
[ 2 . 0 , 4 5 . 0 ] š‘– ā€”ā€” [ āˆ’ 0 . 1 , 0 . 1 0 ] Ī” ā€”
B ā†’ C [ 8 . 0 , 4 2 . 0 ] š‘ [ 0 . 0 , 2 . 5 ] š‘ [ 3 4 1 . 1 5 , 3 4 5 . 1 5 ] [ 0 . 2 0 , 0 . 3 5 ] [ 0 . 0 0 , 2 . 5 0 ]
[ 6 . 0 , 4 5 . 0 ] š‘– [ 0 . 0 , 2 . 6 ] š‘– ā€” [ āˆ’ 0 . 1 , 0 . 1 0 ] Ī” [ āˆ’ 0 . 1 , 0 . 1 0 ] Ī”

In the simulation running, the closed-loop desired poles of system (16) are placed at {šœ†1,šœ†2}={0.72,0.60}, and the corresponding state-feedback matrix š¾ in (3) is computed by š¾=[0.12,0;0,0.6]. Moreover, the parameters in the objective function (4) are chosen as penalty on the states, š‘„=š‘žš¼2, with š‘ž=1, penalty on the inputs, š‘…=š¼3 and a horizon length of š‘‡š‘=20. Assume that grade transition Aā†’B takes place at the tenth hour and Bā†’C at the thirtieth hour. Then some curves of the grade transition process being controlled by NMPC law are shown in Figures 1 and 2, where the solid line in Figure 1 is the set-point, the dotted line the instantaneous, and the dashed line the cumulative.

In order to calculate the quantity of off-specification polymer in grade transition process, we employ the criterion where on-specification polymer is defined as the polymer with MIš‘ and šøš‘”š‘ within a Ā±5% deviation of desired specification in Table 1. All other polymers beyond the criterion are categorized as off-specification production. Similarly, the transition time is defined as the time interval when the cumulative properties of polymer are out of the specification range Ā±5% of the previous grade and lasting until the cumulative properties enter and stay within the specification range Ā±5% of the new target grade. Therefore, the time of grade Aā†’B transition is 8.5 hours and the time of grade Bā†’C transition is 6.5 hours. Compare to the time of 10 hours in actual grade transition operation, this decreases the quantity of off-specification polymer in grade transition process and then increases the productive time of on-specification polymer. Also, Figure 2(a) suggests that there is no change in reaction temperature during the grade transition processes, which is desired in actual operation. Finally, it can be seen that the control laws, intermediate variables, and the states in the grade transition operation do not violate the constraints in Table 1.

4. Conclusion

Together with the pole placement method, the paper presented an efficient MPC algorithm for control of Hammerstein systems subject to constraints on states, inputs, and intermediate variables. Applying pole placement to the linear subsystem of Hammerstein models yielded linear controllers of the intermediate variable and then the actual control action was determined by online solving an optimization problem. The size of the online optimization problem depended only on the number of inputs, not on the predictive horizon, which greatly reduced the computational demand of the algorithm. With condition of feasibility of the optimization problem, the asymptotical stability of the closed-loop system with constraints was guaranteed and the results on simulation example of grade transition control of polypropylene plants showed the effectiveness of the results obtained here.


This work is supported by the National Natural Science Foundation of China under Grant no. 60904040, Specialized Research Fund for Doctoral Program of Higher Education of China under Grant no. 20093317120002 and Natural Science Foundation of Zhejiang Province under Grant no. Y1100911. Finally, the authors are grateful to the reviewsā€™ constructive comments. This paper includes revised and extended text of the first authorā€™s presentation at the 29th CCC, July 27, 2010, Beijiang, China.


  1. Z. Hasiewicz, ā€œHammerstein system identification by the Haar multiresolution approximation,ā€ International Journal of Adaptive Control and Signal Processing, vol. 49, pp. 691ā€“717, 1999. View at: Publisher Site | Google Scholar
  2. N. Bhandari and D. Rollins, ā€œContinuous-time Hammerstein nonlinear modeling applied to distillation,ā€ AIChE Journal, vol. 50, no. 2, pp. 530ā€“533, 2004. View at: Publisher Site | Google Scholar
  3. D. F. He, L. Yu, and M. S. Xue, ā€œOn operation guidance for grade transition in propylene polymerization plants,ā€ in Proceedings of the 19th Chinese Process Control Conference, pp. 339ā€“342, Beijiang, China, 2008, (Chinese). View at: Google Scholar
  4. Q. M. Zhu, K. Warwick, and J. L. Douce, ā€œAdaptive general predictive controller for nonlinear systems,ā€ IEE Proceedings D, vol. 138, no. 1, pp. 33ā€“40, 1991. View at: Google Scholar
  5. K. P. Fruzzetti, A. PalazoĒ§lu, and K. A. McDonald, ā€œNonlinear model predictive control using Hammerstein models,ā€ Journal of Process Control, vol. 7, no. 1, pp. 31ā€“41, 1997. View at: Google Scholar
  6. B. C. Ding, Y. G. Xi, and S. Y. Li, ā€œStability analysis on predictive control of discrete-time systems with input nonlinearity,ā€ Acta Automatica Sinica, vol. 29, no. 6, pp. 827ā€“834, 2003. View at: Google Scholar
  7. B. C. Ding, Y. G. Xi, and S. Y. Li, ā€œOn the stability of output feedback predictive control for systems with input nonlinearity,ā€ Asian Journal of Control, vol. 6, no. 3, pp. 388ā€“397, 2004. View at: Google Scholar
  8. B. Ding and Y. Xi, ā€œA two-step predictive control design for input saturated Hammerstein systems,ā€ International Journal of Robust and Nonlinear Control, vol. 16, no. 7, pp. 353ā€“367, 2006. View at: Publisher Site | Google Scholar
  9. H. T. Zhang, H. X. Li, and G. Chen, ā€œDual-mode predictive control algorithm for constrained Hammerstein systems,ā€ International Journal of Control, vol. 81, no. 10, pp. 1609ā€“1625, 2008. View at: Publisher Site | Google Scholar
  10. H. H. J. Bloemen and T. J. J. van den Boom, ā€œModel-based predictive control for Hammerstein systems,ā€ in Proceedings of the 39th IEEE Conference on Decision and Control, vol. 5, pp. 4963ā€“4968, Sydney, Australia, 2000. View at: Google Scholar
  11. H. H. J. Bloemen, T. J. J. van den Boom, and H. B. Verbruggen, ā€œModel-based predictive control for Hammerstein-Wiener systems,ā€ International Journal of Control, vol. 74, pp. 482ā€“495, 2005. View at: Publisher Site | Google Scholar
  12. D. F. He and L. Yu, ā€œNonlinear predictive control of constrained Hammerstein systems and its research on simulation of polypropylene grade transition,ā€ Acta Automatica Sinica, vol. 35, no. 12, pp. 1558ā€“1563, 2009 (Chinese). View at: Publisher Site | Google Scholar
  13. C. T. Chen, Linear system Theory and Design, Oxford University Press, New York, NY, USA, 3rd edition, 1999.
  14. A. Zheng, ā€œA computationally efficient nonlinear MPC algorithm,ā€ in Proceedings of the American Control Conference, pp. 1623ā€“1627, Albuquerque, NM, USA, June 1997. View at: Google Scholar
  15. J. R. Richards and J. P. Congalidis, ā€œMeasurement and control of polymerization reactors,ā€ Computers and Chemical Engineering, vol. 30, no. 10ā€“12, pp. 1447ā€“1463, 2006. View at: Publisher Site | Google Scholar
  16. H. S. Yi, J. H. Kim, C. Han, J. Lee, and S. S. Na, ā€œPlantwide optimal grade transition for an industrial high-density polyethylene plant,ā€ Industrial and Engineering Chemistry Research, vol. 42, no. 1, pp. 91ā€“98, 2003. View at: Google Scholar
  17. Y. Wang, H. Seki, S. Ohyama, K. Akamatsu, M. Ogawa, and M. Ohshima, ā€œOptimal grade transition control for polymerization reactors,ā€ Computers and Chemical Engineering, vol. 24, no. 2ā€“7, pp. 1555ā€“1561, 2000. View at: Publisher Site | Google Scholar
  18. T. Kreft and W. F. Reed, ā€œPredictive control of average composition and molecular weight distributions in semibatch free radical copolymerization reactions,ā€ Macromolecules, vol. 42, no. 15, pp. 5558ā€“5565, 2009. View at: Publisher Site | Google Scholar

Copyright Ā© 2012 He De-Feng and Yu Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.