Journal of Analytical Methods in Chemistry

Journal of Analytical Methods in Chemistry / 2009 / Article

Research Article | Open Access

Volume 2009 |Article ID 164568 | 11 pages |

A Two-Time Scale Decentralized Model Predictive Controller Based on Input and Output Model

Academic Editor: Peter Stockwell
Received11 May 2009
Accepted02 Jul 2009
Published30 Sep 2009


A decentralized model predictive controller applicable for some systems which exhibit different dynamic characteristics in different channels was presented in this paper. These systems can be regarded as combinations of a fast model and a slow model, the response speeds of which are in two-time scale. Because most practical models used for control are obtained in the form of transfer function matrix by plant tests, a singular perturbation method was firstly used to separate the original transfer function matrix into two models in two-time scale. Then a decentralized model predictive controller was designed based on the two models derived from the original system. And the stability of the control method was proved. Simulations showed that the method was effective.

1. Introduction

Applications of singular perturbation in control theory can be traced back to 1970s [1]. When there are both fast and slow dynamics in a system, the perturbation conception can be adopted to simplify the system. Phillips first combined the optimal control theory with a singular perturbation system and presented a two-stage design of linear feedback control [2]. Later, nonlinear singular perturbation systems [3], large scale systems [4], and high-gain feedback systems [5] were investigated with optimal control theory, too. Recently, process control researchers took notice of the mechanical reasons for two-time scale characteristics of systems. And the existence of fast and slow dynamics was studied in several previous papers. Yi and Luyben analyzed the dynamic characteristic of coupled reactor/column systems in a series of papers and illuminated the mechanisms of fast and slow dynamics in some special systems [6โ€“8]. Contou-Carrere and Daoutidis focused on the coexistence of fast and slow dynamics in integrated process networks [9]. They explained that large flow rates brought in fast dynamic and led to a time scales separation of dynamics. They also designed a precompensator for a distillation model. Sequentially, Kumar and Daoutidis further analyzed the dynamics of process with material and energy recycle, and introduced a controller design framework consisting of properly coordinated controllers in fast and slow time scales [10]. Vora and Daoutidis introduced a nonstandard form singular perturbation method and analyzed the existence of fast and slow dynamics in a nonlinear system [11]. Kumar and Christofides dealt with two-time scale chemical processes and modeled them by using nonlinear ordinary differential equations with large parameters of the form 1/๐œ€ [12]. They obtained a standard singularly perturbed representation. All those methods focused on process thermodynamic models for which took much work to get precise data.

In actual processes, model predictive control (MPC) was regarded as โ€œthe only advanced control methodology which has made a significant impact on industrial control engineeringโ€ [13]. And MPC has no limit of the model form. The input-output model which can be obtained easily by identification is usually used. However, all the methods mentioned above is from the point of optimal control which is an offline method. In addition, all the above methods regard the steady state output of fast dynamics as the input or known quantity of slow dynamics but pay little attention on the influence of slow dynamics on fast dynamics. MPC adopts online rolling optimization and uses feedback to correct prediction. The challenge of applying MPC in singular perturbation system becomes the compromise between the control interval and the predictive horizon. Fast dynamic needs small control interval, while slow dynamic needs large predictive horizon. There are few papers about two-time scale MPC. Buescher and Baum introduced a two-time scale approach to nonlinear model predictive control, and used a โ€œgappingโ€ method to smooth the control quality [14]. But they did not consider the dynamic characteristics of the model in their control algorithm which may be not suitable for special systems.

Therefore, we focused on designing an MPC for systems with two-time scale characteristic in this paper. First, we introduced the background of this field. Then we described the two-time scale decomposition of a transfer function matrix with different dynamics in different channels. At part 3, we presented a kind of two-time scale decentralized MPC algorithm step by step and proved its stability. At last, we gave several simulations to test the validity of the two-time scale decentralized MPC algorithm.

2. Two-Time Scale Decomposition of a Transfer Function Matrix

In some systems, the dynamics varies with different channels. And the response speeds of those different channels vary so much even in different time scales. The characteristics can be got from the transfer function matrix intuitively. We simply took a two-in-two-out first-order transfer function matrix, for example and gave the following definition.

Definition 1. A Two-in-two-out first-order transfer function matrix ๐บ(๐‘ ) with two-time scale characteristic is presented below, where ๐บ(๐‘ )=๎‚ธ๐บ11(๐‘ )๐บ12(๐‘ )๐บ21(๐‘ )๐บ22(๐‘ )๎‚น๐บ11(๐‘ )=๐‘Ž11๐‘11๐‘ +1๐‘’๐œ11๐‘ ,๐บ12(๐‘ )=๐‘Ž12๐‘12๐‘ +1๐‘’๐œ12๐‘ ,๐บ21(๐‘ )=๐‘Ž21๐‘21๐‘ +1๐‘’๐œ21๐‘ ,๐บ22(๐‘ )=๐‘Ž22๐‘22๐‘ +1๐‘’๐œ22๐‘ ,(1) when ๐‘Ž11, ๐‘Ž12, ๐‘Ž21, ๐‘Ž22 are in the same order, ๐‘11,๐‘12โ‰ช๐‘21,๐‘22, ๐œ11,๐œ12,๐œ21,๐œ22 are in the same order.
Considering the system
๐‘Œ=๐บ(๐‘ )๐‘ˆ,(2) the response speed of output ๐‘ฆ1 is much faster than that of output ๐‘ฆ2. Then a traditional central controller cannot satisfy the demands of ๐‘ฆ1 and ๐‘ฆ2 simultaneously. For example, for the fast channel a very short control interval is needed to provide enough dynamic characteristics, while for the slow channel a very large predictive horizon is needed to ensure the stability of the controller. Thus a normal model predictive controller cannot satisfy both demands. Therefore we designed a decentralized controller based on a two-time scale method.
Some papers [9, 11, 15] introduced a singular perturbation method to obtain two-time scale models which are based on state space model. In order to transform the transfer function matrix into a two-time scale form, the transfer function matrix form should be transferred into a state space form. But the delay terms cannot be expressed in a state space form. Thus we firstly took a transfer function matrix without delay to illustrate the two-time scale decomposition, and then discussed the situation with delay terms.

2.1. Without Delay Terms

Based on the method mentioned in literature[16], system ๐บ๎…ž(๐‘ )=๎‚ธ๐บ๎…ž11(s)๐บ๎…ž12(s)๐บ๎…ž21(s)๐บ๎…ž22(s)๎‚น,(3) where ๐บ๎…ž11(s)=๐‘Ž11/๐‘11๐‘ +1, ๐บ๎…ž12(s)=๐‘Ž12/๐‘12๐‘ +1, ๐บ๎…ž21(s)=๐‘Ž21/๐‘21๐‘ +1, ๐บ๎…ž22(s)=๐‘Ž22/๐‘22๐‘ +1, ๐‘Ž11, ๐‘Ž12, ๐‘Ž21, ๐‘Ž22 are in the same order, ๐‘11,๐‘12โ‰ช๐‘21,๐‘22, can have this form

ฬ‡โ€Œ๐‘‹=๐ด๐‘‹+๐ต๐‘ˆ,๐‘Œ=๐ถ๐‘‹,(4) where ๐ด=diag๎€บ๐ด11,๐ด12,๐ด21,๐ด22๎€ป,๐ต=โŽกโŽขโŽขโŽขโŽฃ๐ต11๐ต12๐ต21๐ต22โŽคโŽฅโŽฅโŽฅโŽฆ,๐ถ=๎‚ธ๐ถ11๐ถ12๐ถ21๐ถ22๎‚น;(5)(๐ด๐‘–๐‘—,๐ต๐‘–๐‘—,๐ถ๐‘–๐‘—) is a one-dimension state space form of ๐บ๎…ž๐‘–๐‘—(s):


So we can rewrite (4) in this form:

โŽกโŽขโŽขโŽขโŽฃ๐œ€ฬ‡๐‘ฅ1๐œ€ฬ‡๐‘ฅ2ฬ‡๐‘ฅ3ฬ‡๐‘ฅ4โŽคโŽฅโŽฅโŽฅโŽฆ=โŽกโŽขโŽขโŽขโŽขโŽขโŽขโŽฃ๐œ€๐ด110๐œ€๐ด12๐ด210๐ด22โŽคโŽฅโŽฅโŽฅโŽฅโŽฅโŽฅโŽฆโŽกโŽขโŽขโŽขโŽฃ๐‘ฅ1๐‘ฅ2๐‘ฅ3๐‘ฅ4โŽคโŽฅโŽฅโŽฅโŽฆ+โŽกโŽขโŽขโŽขโŽฃ๐œ€๐ต11๐œ€๐ต12๐ต21๐ต22โŽคโŽฅโŽฅโŽฅโŽฆ๎‚ธ๐‘ข1๐‘ข2๎‚น,๎‚ธ๐‘ฆ1๐‘ฆ2๎‚น=๎‚ธ๐ถ11๐ถ12๐ถ21๐ถ22๎‚นโŽกโŽขโŽขโŽขโŽฃ๐‘ฅ1๐‘ฅ2๐‘ฅ3๐‘ฅ4โŽคโŽฅโŽฅโŽฅโŽฆ,(7) where ๐œ€ is a very small positive constant.

This form can be regarded as

ฬ‡โ€Œ๐‘‹=๐ด11๐‘‹+๐ด12๐‘+๐ต1๐‘ˆ,๐œ€ฬ‡โ€Œ๐‘=๐ด21๐‘‹+๐ด22๐‘+๐ต2๐‘ˆ,๐‘Œ=๐ถ1๐‘‹+๐ถ2๐‘+๐ท๐‘ˆ,(8) where ๐ด11=๎‚ธ๐ด2100๐ด22๎‚น,๐ด12=๎‚ธ0000๎‚น,๐ด21=๎‚ธ0000๎‚น,๐ด22=๎‚ธ๐œ€๐ด1100๐œ€๐ด12๎‚น,๐ต1=๎‚ธ๐ต2100๐ต22๎‚น,๐ต2=๎‚ธ๐œ€๐ต1100๐œ€๐ต12๎‚น,๐ถ1=๎‚ธ00๐ถ21๐ถ22๎‚น,๐ถ2=๎‚ธ๐ถ11๐ถ1200๎‚น,๐ท=๎‚ธ0000๎‚น.(9) And the system in this form was discussed in several literatures [9, 15, 17], too. We got the slow model

๐‘‘๐‘‹๐‘ ๐‘‘๐‘ก=๐ด๐‘ ๐‘‹๐‘ +๐ต๐‘ ๐‘ˆ,๐‘Œ=๐ถ๐‘ ๐‘‹๐‘ +๐ท๐‘ ๐‘ˆ,(10) and the fast model

๐‘‘๐‘๐‘“๐‘‘๐œ=๐ด๐‘“๐‘๐‘“+๐ต๐‘“๐‘ˆ,๐‘Œ=๐ถ๐‘“๐‘๐‘“+๐ท๐‘“๐‘ˆ,(11) where ๐ด๐‘ =๐ด11โˆ’๐ด12๐ดโˆ’122๐ด21,๐ต๐‘ =๐ต1โˆ’๐ด12๐ดโˆ’122๐ต2,๐ถ๐‘ =๐ถ1โˆ’๐ถ2๐ดโˆ’122๐ด21,๐ท๐‘ =๐ทโˆ’๐ถ2๐ดโˆ’122๐ต2,๐ด๐‘“=๐ด22,๐ต๐‘“=๐ต2,๐ถ๐‘“=๐ถ2,๐ท๐‘“=๐ท.(12)

We denoted the transfer function of the slow model ๐บ๐‘ (๐‘ ) by

๐บ๐‘ (๐‘ )=๐ถ๐‘ ๎€ท๐‘ ๐ผโˆ’๐ด๐‘ ๎€ธโˆ’1๐ต๐‘ +๐ท๐‘ =๎‚ธ๐‘Ž11๐‘Ž12๐บ๎…ž21(๐‘ )๐บ๎…ž22(๐‘ )๎‚น(13) and the transfer function of the fast model ๐บ๐‘“(๐‘ ) in ๐œ time scale by

๐บ๐‘“(๐‘ )=๐ถ๐‘“๎€ท๐‘ ๐ผโˆ’๐ด๐‘“๎€ธโˆ’1๐ต๐‘“+๐ท๐‘“=๐ถ2๎‚€๐‘ ๐ผโˆ’๐ด22๎‚โˆ’1๐ต2+๐ท=๐ถ2๎‚€๐‘ ๐ผโˆ’๐ด22๎‚โˆ’1๐ต2=โŽกโŽขโŽฃ๐œ€๐‘Ž11๐‘11๐‘ +๐œ€๐œ€๐‘Ž12๐‘12๐‘ +๐œ€00โŽคโŽฅโŽฆ,(14) and the fast model in t-time scale is

๐บ๐‘Ÿ(๐‘ )=๐บ๐‘“(๐œ€๐‘ )=โŽกโŽขโŽฃ๐‘Ž11๐‘11๐‘ +1๐‘Ž12๐‘12๐‘ +100โŽคโŽฅโŽฆ.(15)

๐บ๐‘“(๐‘ ), ๐บ๐‘Ÿ(๐‘ ), and ๐บ๐‘ (๐‘ ) are all descriptions from different points of the real system ๐บ๎…ž(๐‘ ). And they are all two-in-two-out systems but have less state variables than the real system ๐บ๎…ž(๐‘ ). If we choose a long enough sample interval for ๐บ๎…ž(๐‘ ) and ๐บ๐‘ (๐‘ ), the sample values can be similar, because the fast responses turn to steady states that can be regarded as constants in a very short time. If we choose a very short sample interval for ๐บ๎…ž(๐‘ ) and ๐บ๐‘Ÿ(๐‘ ) in a short enough period of time, the sample values can be similar, too. Because in such a short period of time the slow channel response is so slow that can be regarded as zero. But the output of the slow channel is not zero. Therefore, we modified the fast model in t-time scale when we designed a controller based on this fast model. Literature [18] proved the relationship between ๐บ๐‘ (๐‘ ) and ๐บ๐‘“(๐‘ ), and also the transfer function of the original system ๐บ๎…ž(๐‘ ):

lim๐‘ โ†’โˆž๐บ๐‘ (๐‘ )=๐ท๐‘ =lim๐‘ โ†’0๐บ๐‘“(๐‘ ),(16)๐บ๎…ž(๐‘ )=๐บ๐‘ (๐‘ )+๐บ๐‘“(๐œ€๐‘ )โˆ’๐ท๐‘ +๐‘‚(๐œ€).(17)

Equation (16) denotes that the initial value of the slow model equals to the final value of the fast model. And (17) denotes that the original system can be regarded as a sum of the slow model, the fast model, and a very little item ๐‘‚(๐œ€). Let

๐‘‚(๐œ€)=๎‚ธ00๐‘‚21(๐œ€)๐‘‚22(๐œ€)๎‚น,(18) and let

๐บ๐‘ก(๐‘ )=๐บ๐‘Ÿ(๐‘ )+๐‘‚(๐œ€)=โŽกโŽขโŽฃ๐‘Ž11๐‘11๐‘ +1๐‘Ž12๐‘12๐‘ +1๐‘‚21(๐œ€)๐‘‚22(๐œ€)โŽคโŽฅโŽฆ=๎‚ธ๐บ๎…ž11(๐‘ )๐บ๎…ž12(๐‘ )๐‘‚21(๐œ€)๐‘‚22(๐œ€)๎‚น.(19)

We had expressions of the fast model ๐บ๐‘ก(๐‘ ) and the slow model ๐บ๐‘ (๐‘ ) about the model ๐บ๎…ž(๐‘ ) without delay terms in t-time scale. Next we would consider the model with delay terms.

2.2. With Delay Terms

We took system (2), where ๐บ(๐‘ ) shows two-time scale characteristic, as an example. Let

๐‘‹11(๐‘ )=๐บ11(๐‘ )๐‘ˆ1(๐‘ )=๐บ๎…ž11(๐‘ )๐‘ˆ1(๐‘ )๐‘’๐œ11๐‘ ,๐‘‹12(๐‘ )=๐บ12(๐‘ )๐‘ˆ2(๐‘ )=๐บ๎…ž12(๐‘ )๐‘ˆ2(๐‘ )๐‘’๐œ12๐‘ ,๐‘‹21(๐‘ )=๐บ21(๐‘ )๐‘ˆ1(๐‘ )=๐บ๎…ž21(๐‘ )๐‘ˆ1(๐‘ )๐‘’๐œ21๐‘ ,๐‘‹22(๐‘ )=๐บ22(๐‘ )๐‘ˆ2(๐‘ )=๐บ๎…ž22(๐‘ )๐‘ˆ2(๐‘ )๐‘’๐œ22๐‘ ,(20)

Here ๐บ๎…ž๐‘–๐‘—(๐‘ ), ๐‘–,๐‘—=1,2, denotes the transition process, and ๐‘’๐œ๐‘–๐‘—๐‘ , ๐‘–,๐‘—=1,2, denotes the delay time. When ๐œ11,๐œ12,๐œ21,and๐œ22 are in the same order, ๐บ๎…ž๐‘–๐‘—(s), ๐‘–,๐‘—=1,2, reflects the main dynamic characteristic. We can get the fast model in t-time scale

๐บ๐‘ก๐‘‘(๐‘ )=๎‚ธ๐บ๎…ž11(๐‘ )๐‘’๐œ11๐‘ ๐บ๎…ž12(๐‘ )๐‘’๐œ12๐‘ ๐‘‚21(๐œ€)๐‘‚22(๐œ€)๎‚น.(21) and the slow model in t-time scale

๐บ๐‘ ๐‘‘(๐‘ )=๎‚ธ๐‘Ž11๐‘Ž12๐บ๎…ž21(๐‘ )๐‘’๐œ21๐‘ ๐บ๎…ž22(๐‘ )๐‘’๐œ22๐‘ ๎‚น(22)

Then we can design a decentralized controller based on characteristics of the fast model ๐บ๐‘ก(๐‘ )(๐บ๐‘ก๐‘‘(๐‘ )) and the slow model ๐บ๐‘ (๐‘ )(๐บ๐‘ ๐‘‘(๐‘ )).

3. Two-Time Scale Decentralized MPC

MPC is the only advanced control methodology which has made a significant impact on industrial control engineering [13]. And MPC is based on a predicted model. The MPC algorithm can be regarded as a combination of three parts: model prediction, roll optimization, and feedback rectification [19].

For the two-in-two-out system mentioned above with two-time scale characteristic, we designed a decentralized controller based on different time scales. We took the model without delay (๐บ๎…ž(๐‘ )) to illustrate the algorithm. The fast model provided abundant fast dynamic information to ensure the control quality. And the control interval was determined by the fast model. The slow model provided prediction horizon long enough to ensure the controllerโ€™s stability. In order to illustrate the algorithm, we defined that ๐‘ƒ๐‘  is the prediction horizon, ๐‘‡๐‘  is the sampling interval based on the slow model, ๐‘ƒ๐‘“ is the prediction horizon, ๐‘‡๐‘“ is the sampling interval based on fast model, and ๐‘€ is the manipulate horizon. ๐‘Ž๐‘–,๐‘—(๐‘ก) is the step response of ๐‘ฆ๐‘– from ๐‘ข๐‘— at ๐‘‡๐‘“ sample interval. We got the model vector ๐‘Ž๐‘–,๐‘—=[๐‘Ž๐‘–,๐‘—(1)โ‹ฏ๐‘Ž๐‘–,๐‘—(๐‘)]๐‘‡, ๐‘–=1,2, ๐‘—=1,2, and ๐‘ is a number large enough to fully reflect the fast and the slow part of the model. The two-time scale decentralized DMC algorithm was shown as follows.

Step 1. Model prediction based on slow model ๐บ๐‘ (๐‘ ).
On the t-time scale, the fast dynamic achieved a steady state. The model can be fully expressed by information not so necessary as ๐‘Ž๐‘–,๐‘—. Let ๐‘Ž๐‘ ,๐‘—=[๐‘Ž๐‘ ,๐‘—(1)โ‹ฏ๐‘Ž๐‘ ,๐‘—(๐‘ƒ๐‘ )]๐‘‡, ๐‘—=1,2, be the slow model vector, where ๐‘Ž๐‘ ,๐‘—(๐‘–)=[๐‘Ž2,๐‘—((๐‘–โˆ’1)โˆ—(๐‘‡๐‘ /๐‘‡๐‘“)+1)โ‹ฏ๐‘Ž2,๐‘—((๐‘–โˆ’1)โˆ—(๐‘‡๐‘ /๐‘‡๐‘“)+๐‘€)], ๐‘–=1,โ€ฆ,๐‘ƒ๐‘ . The slow predict model is ฬƒ๐‘ฆ๐‘ ,๐‘ƒ๐‘€(๐‘)=ฬƒ๐‘ฆ๐‘ ,๐‘ƒ0(๐‘)+๐ด๐‘ ฮ”๐‘ข,๐‘—=1,2,(23) where, ฬƒ๐‘ฆ๐‘ ,๐‘ƒ๐‘€(๐‘)โˆˆ๐‘…๐‘ƒ๐‘  is the predict value of future output and ฬƒ๐‘ฆ๐‘ ,๐‘ƒ0(๐‘)โˆˆ๐‘…๐‘ƒ๐‘  is the prime predict value of future output: ๐ด๐‘ =โŽกโŽขโŽขโŽขโŽขโŽขโŽฃ๐‘Ž๐‘ ,๐‘—(1)0โ‹ฎโ‹ฑ๐‘Ž๐‘ ,๐‘—(๐‘€)โ‹ฏ๐‘Ž๐‘ ,๐‘—(1)โ‹ฎโ‹ฎ๐‘Ž๐‘ ,๐‘—๎€ท๐‘ƒ๐‘ ๎€ธโ‹ฏ๐‘Ž๐‘ ,๐‘—๎€ท๐‘ƒ๐‘ โˆ’๐‘€+1๎€ธโŽคโŽฅโŽฅโŽฅโŽฅโŽฅโŽฆ.(24)

Step 2. Feedback Correction based on slow model ๐บ๐‘ (๐‘ ).
Let error vector be ๐‘’๐‘ (๐‘+1)=๐‘ฆ2(๐‘+1)โˆ’ฬƒ๐‘ฆ๐‘ ,๐‘ƒ๐‘€(๐‘+1๐‘),(25) and we can get ฬƒ๐‘ฆ๐‘ ,cor(๐‘+1)=ฬƒ๐‘ฆ๐‘ ,๐‘ƒ๐‘€(๐‘)+๐ป๐‘ ๐‘’๐‘ (๐‘+1),(26) where ๐ป๐‘ =โŽกโŽขโŽขโŽฃโ„Ž๐‘ (1)โ‹ฎโ„Ž๐‘ ๎€ท๐‘ƒ๐‘ ๎€ธโŽคโŽฅโŽฅโŽฆ,(27) Like in Step 4, at ๐‘+1 time point, the time origin changes from ๐‘ to ๐‘+1 time point, then the elements of vector ฬƒ๐‘ฆ๐‘Ÿ,cor(๐‘+1) should be moved, and the operation can be expressed by ฬƒ๐‘ฆ๐‘ ,๐‘0(๐‘+1)=๐‘†๐‘ ,0ฬƒ๐‘ฆ๐‘ ,cor(๐‘+1),(28) where ๐‘†๐‘ ,0=โŽกโŽขโŽขโŽขโŽฃ010โ‹ฎโ‹ฑโ‹ฑโ‹ฎ010โ‹ฏ01โŽคโŽฅโŽฅโŽฅโŽฆ.(29)

Step 3. Model prediction based on fast model ๐บ๐‘ก(๐‘ ).
On the ๐œ-time scale, the slow dynamic can be regarded as 0 which means that the slow dynamic changes very little. The model also can be fully expressed by information not so necessary as ๐‘Ž๐‘–,๐‘—. Let ๐‘Ž1,๐‘—=[๐‘Ž1,๐‘—(1)โ‹ฏ๐‘Ž1,๐‘—(๐‘๐‘Ÿ)]๐‘‡, ๐‘—=1,2. The fast predict model is ฬƒ๐‘ฆ๐‘Ÿ,๐‘ƒ๐‘€(๐‘˜)=ฬƒ๐‘ฆ๐‘Ÿ,๐‘ƒ0(๐‘˜)+๐ด๐‘Ÿฮ”๐‘ข,๐‘—=1,2,(30) where ๐ด๐‘Ÿ=โŽกโŽขโŽขโŽขโŽขโŽขโŽฃ๐‘Ž๐‘Ÿ,๐‘—(1)0โ‹ฎโ‹ฑ๐‘Ž๐‘Ÿ,๐‘—(๐‘€)โ‹ฏ๐‘Ž๐‘Ÿ,๐‘—(1)โ‹ฎโ‹ฎ๐‘Ž๐‘Ÿ,๐‘—๎€ท๐‘ƒ๐‘“๎€ธโ‹ฏ๐‘Ž๐‘Ÿ,๐‘—๎€ท๐‘ƒ๐‘“โˆ’๐‘€+1๎€ธโŽคโŽฅโŽฅโŽฅโŽฅโŽฅโŽฆ.(31)

Step 4. Feedback Correction based on fast model ๐บ๐‘ก(๐‘ ).
The input ๐‘ข(๐‘˜) was applied to the plant ๐บ๎…ž(s) at each ๐‘˜ time point, and (22) gave the predictive output ฬƒ๐‘ฆ๐‘Ÿ,๐‘ƒ๐‘€(๐‘˜+1โˆฃ๐‘˜). Let error vector be ๐‘’๐‘Ÿ(๐‘˜+1)=๐‘ฆ1(๐‘˜+1)โˆ’ฬƒ๐‘ฆ๐‘Ÿ,๐‘ƒ๐‘€(๐‘˜+1โˆฃ๐‘˜),(32) where ๐‘ฆ1(๐‘˜+1) is the sample value. We used this error vector to modify the infection by some unsure factors, and we gave the error a weight vector to modify the prediction of output: ฬƒ๐‘ฆ๐‘Ÿ,cor(๐‘˜+1)=ฬƒ๐‘ฆ๐‘Ÿ,๐‘ƒ๐‘€(๐‘˜)+๐ป๐‘Ÿ๐‘’๐‘Ÿ(๐‘˜+1),(33) where ๐ป๐‘Ÿ=[โ„Ž๐‘Ÿ(1)โ‹ฎโ„Ž๐‘Ÿ(๐‘ƒ๐‘“)]. At ๐‘˜+1 time point, the time origin changes from ๐‘˜ to ๐‘˜+1 time point, then the elements of vector ฬƒ๐‘ฆ๐‘Ÿ,cor(๐‘˜+1) should be moved, and the operation can be expressed by ฬƒ๐‘ฆ๐‘Ÿ,๐‘ƒ0(๐‘˜+1)=๐‘†๐‘Ÿ,0ฬƒ๐‘ฆ๐‘Ÿ,cor(๐‘˜+1),(34) where ๐‘†๐‘Ÿ,0=โŽกโŽขโŽขโŽขโŽฃ010โ‹ฎโ‹ฑโ‹ฑโ‹ฎ010โ‹ฏ01โŽคโŽฅโŽฅโŽฅโŽฆ.(35)

Step 5. Rolling horizon optimization based on fast model ๐บ๐‘ก(๐‘ ).
Let objective function be min๐ฝ(๐‘˜)=โ€–โ€–๐œ”1(๐‘˜)โˆ’ฬƒ๐‘ฆ๐‘Ÿ,๐‘ƒ๐‘€(๐‘˜)โ€–โ€–2๐‘„1+โ€–โ€–๐œ”2(๐‘)โˆ’ฬƒ๐‘ฆ๐‘ ,๐‘ƒ๐‘€(๐‘)โ€–โ€–2๐‘„2+โ€–ฮ”๐‘ข(๐‘˜)โ€–2๐‘…s.t.ฮ”๐‘ขminโ‰คฮ”๐‘ขโ‰คฮ”๐‘ขmax,๐‘ขminโ‰ค๐‘ขโ‰ค๐‘ขmax,(36) where ๐œ”๐‘–(๐‘˜), ๐‘–=1,2, are the reference value: ๐‘„1=diag๎€บ๐‘ž1(1)โ‹ฏ๐‘ž1๎€ท๐‘ƒ๐‘“๎€ธ๎€ป,๐‘„2=diag๎€บ๐‘ž2(1)โ‹ฏ๐‘ž2๎€ท๐‘ƒ๐‘ ๎€ธ๎€ป,๐‘…=blockโˆ’diag๎€ท๐‘…1,๐‘…2๎€ธ,๐‘…๐‘–=diag๎€บ๐‘Ÿ๐‘–(1)โ‹ฏ๐‘Ÿ๐‘–๎€ท๐‘€๐‘“๎€ธ๎€ป,๐‘–=1,2.(37) Without constraint, we can get the manipulated variables: ฮ”๐‘ข(๐‘˜)=๐ฟ๎‚†๎€ท๐ด๐‘‡๐‘Ÿ๐‘„1๐ด๐‘Ÿ+๐ด๐‘‡๐‘ ๐‘„2๐ด๐‘ +๐‘…๎€ธโˆ’1ร—๎€บ๐ด๐‘‡๐‘Ÿ๐‘„1๎€ท๐œ”1(๐‘˜)โˆ’ฬƒ๐‘ฆ๐‘Ÿ,๐‘ƒ0(๐‘˜)๎€ธ+๐ด๐‘‡๐‘ ๐‘„2๎€ท๐œ”2(๐‘)โˆ’ฬƒ๐‘ฆ๐‘ ,๐‘ƒ0(๐‘)๎€ธ๎€ป๎‚‡,(38) where ๐ฟ=โŽกโŽขโŽขโŽขโŽฃ10โ‹ฏ0๎„ฟ๎…€๎…€๎…€๎…€๎…€๎…€๎…€๎…ƒ๎…€๎…€๎…€๎…€๎…€๎…€๎…€๎…Œ๐‘€0010โ‹ฏ0๎„ฟ๎…€๎…€๎…€๎…€๎…€๎…€๎…€๎…ƒ๎…€๎…€๎…€๎…€๎…€๎…€๎…€๎…Œ๐‘€โŽคโŽฅโŽฅโŽฅโŽฆ.(39) Step 3 to step 5 form the inner circulate, and in every ๐‘‡๐‘“ time interval fast controller calculates a manipulate variable. Step 1 and Step 2 form the exterior circulate, and the slow controller provides the predictive value of the slow output to correct the fast controller every ๐‘›๐‘‡๐‘“ (n is a positive integer) time interval. So the framework of the control system was showed as shown in Figure 1.

4. Stability of the Two-Time Scale Decentralized MPC

We introduced the algorithm step by step in the above sections. And we would like to discuss the stability of the controller in this section. First we put forward a sufficient condition of the controller with one inner circulate.

Theorem 1. A two-time scale decentralized MPC with the control parameters ๐‘‡๐‘“, ๐‘‡๐‘ , and ๐‘ƒ๐‘  and one inner circulate is stable under the sufficient condition of a standard MPC, which has a control interval ๐‘‡๐‘“ and a predictive horizon ๐‘ƒ๐‘ (๐‘‡๐‘ /๐‘‡๐‘“).

Proof. Let the prediction horizon ๐‘ƒ๐‘ (๐‘‡๐‘ /๐‘‡๐‘“), be the manipulate horizon, let ๐‘€ be the objective function of the standard MPC be min๐ฝ๎…ž(๐‘˜)=โ€–โ€–๐œ”(๐‘˜)โˆ’ฬƒ๐‘ฆ๐‘ƒ๐‘€(๐‘˜)โ€–โ€–2๐‘„+โ€–ฮ”๐‘ข(๐‘˜)โ€–2๐‘…s.t.ฮ”๐‘ขminโ‰คฮ”๐‘ขโ‰คฮ”๐‘ขmax,๐‘ขminโ‰ค๐‘ขโ‰ค๐‘ขmax,(40)๐ฝ๎…ž(๐‘˜)=โ€–โ€–๐œ”(๐‘˜)โˆ’ฬƒ๐‘ฆ๐‘ƒ๐‘€(๐‘˜)โ€–โ€–2๐‘„+โ€–ฮ”๐‘ข(๐‘˜)โ€–2๐‘…=โ€–โ€–๐œ”1(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘˜)โ€–โ€–2๐‘„๎…ž1+โ€–โ€–๐œ”2(๐‘)โˆ’ฬƒ๐‘ฆ2,๐‘ƒ๐‘€(๐‘)โ€–โ€–2๐‘„๎…ž2+โ€–ฮ”๐‘ข(๐‘˜)โ€–2๐‘…=๐‘ƒ๐‘ (๐‘‡๐‘ /๐‘‡๐‘“)๎“๐‘–=1๐‘ž๎…ž1๎€ท๐œ”1(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘–+๐‘˜โˆฃ๐‘˜)๎€ธ2+๐‘ƒ๐‘ (๐‘‡๐‘ /๐‘‡๐‘“)๎“๐‘–=1๐‘ž๎…ž2๎€ท๐œ”2(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘–+๐‘˜โˆฃ๐‘˜)๎€ธ2+๐‘€๎“๐‘–=1๐‘Ÿฮ”๐‘ข(๐‘˜+๐‘–โˆฃ๐‘˜)2.(41) Because the response speed of ๐‘ฆ1 is very fast, the predictive value ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘–+๐‘˜โˆฃ๐‘˜) can be a fixed number when ๐‘– is larger than a certain number (๐‘ƒ๐‘“): ๐‘ƒ๐‘ (๐‘‡๐‘ /๐‘‡๐‘“)๎“๐‘–=1๐‘ž๎…ž1๎€ท๐œ”1(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘–+๐‘˜โˆฃ๐‘˜)๎€ธ2=๐‘ƒ๐‘“๎“๐‘–=1๐‘ž๎…ž1๎€ท๐œ”1(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘–+๐‘˜โˆฃ๐‘˜)๎€ธ2+Const.(42)
The response speed of ๐‘ฆ2 is very slow, so between a short time interval (๐‘‡๐‘ ) the predictive values ฬƒ๐‘ฆ2,๐‘ƒ๐‘€(๐‘–+๐‘˜โˆฃ๐‘˜), ๐ถโ‰ค๐‘–โ‰ค๐ถ+(๐‘‡๐‘ /๐‘‡๐‘“) (๐ถ is a positive number), can be linear correlation:
๐ถ+(๐‘‡๐‘ /๐‘‡๐‘“)๎“๐‘–=๐ถ๐‘ž๎…ž2๎€ท๐œ”2(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘–+๐‘˜โˆฃ๐‘˜)๎€ธ2=๐‘ž2๎€ท๐œ”2(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐ถ+๐‘˜โˆฃ๐‘˜)๎€ธ2,(43)๐‘ƒ๐‘ (๐‘‡๐‘ /๐‘‡๐‘“)๎“๐‘–=1๐‘ž๎…ž2๎€ท๐œ”2(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘–+๐‘˜โˆฃ๐‘˜)๎€ธ2=๐‘ƒ๐‘ ๎“๐‘—=1๐‘ž2๎€ท๐œ”2(๐‘˜)โˆ’ฬƒ๐‘ฆ1,๐‘ƒ๐‘€(๐‘—+๐‘˜โˆฃ๐‘˜)๎€ธ.(44) Combining (38), (40), (42), and (44), we got ๐ฝ๎…ž(๐‘˜)=๐ฝ(๐‘˜)+const.(45)
So min๐ฝ(๐‘˜) and min๐ฝ๎…ž(๐‘˜) had same answers, and the sufficient condition was proved.

When the inner circulate is large than one, that is, ๐‘›>1, the algorithm can maintain its stability if the slow dynamic changes little in the time interval of the inner circulate. But it is hard to find an upper limit for ๐‘›, because it is determined by the characteristic of the slow dynamic which can be quite different in different systems.

5. Case Study

A Model with Delay
We considered a two-in-two-out system. Two streams flow into a reactor, and ๐‘ข1 and ๐‘ข2 are the flow rates. The liquid level ๐‘ฆ1 and the temperature ๐‘ฆ2 are two controlled variables. ๐‘ฆ10=50cm and ๐‘ฆ20=295K are the initial stable states. The linear model of the system is ๎‚ธ๐‘ฆ1๐‘ฆ2๎‚น=๐บ(๐‘ )๎‚ธ๐‘ข1๐‘ข2๎‚น,(46) where ๐บ(๐‘ )=โŽกโŽขโŽขโŽขโŽฃ๐พ11๐‘‡11๐‘ +1๐‘’โˆ’๐œ11๐‘ ๐พ12๐‘‡12๐‘ +1๐‘’โˆ’๐œ12๐‘ ๐พ21๐‘‡21๐‘ +1๐‘’โˆ’๐œ21๐‘ ๐พ22๐‘‡22๐‘ +1๐‘’โˆ’๐œ22๐‘ โŽคโŽฅโŽฅโŽฅโŽฆ,(47)
The response of output ๐‘ฆ1 is much faster than that of ๐‘ฆ2. If a standard DMC controller is applied on this system, the sample interval is determined by the fast response and should be very small, and the predictive horizon is determined by the slow response and should be very large. In such a small interval, it is difficult to calculate the optimal manipulate variables, and the control quality may be bad. If we compromise the sample time interval of different channels, we can get the following control effect to track step signals. Liquid level ๐‘ฆ1 is set as 51โ€‰cm, and temperature is set as 296โ€‰K. Considering the uncertainty of the model, we chose the plant that each parameter above has 20% uncertainty to carry through the simulation. The plant is
๐‘ƒ(๐‘ )=โŽกโŽขโŽขโŽขโŽขโŽฃ๐พ11+0.2๐พ11๐œ€๐‘’โˆ’๎€ท๐œ11+0.2๐œ11๐œ€๎€ธ๐‘ ๎€ท๐‘‡11+0.2๐‘‡11๐œ€๎€ธ๐‘ +1๐พ12+0.2๐พ12๐œ€๐‘’โˆ’๎€ท๐œ12+0.2๐œ12๐œ€๎€ธ๐‘ ๎€ท๐‘‡12+0.2๐‘‡12๐œ€๎€ธ๐‘ +1๐พ21+0.2๐พ21๐œ€๐‘’โˆ’๎€ท๐œ21+0.2๐œ21๐œ€๎€ธ๐‘ ๎€ท๐‘‡21+0.2๐‘‡21๐œ€๎€ธ๐‘ +1๐พ22+0.2๐พ22๐œ€๐‘’โˆ’๎€ท๐œ22+0.2๐œ22๐œ€๎€ธ๐‘ ๎€ท๐‘‡22+0.2๐‘‡22๐œ€๎€ธ๐‘ +1โŽคโŽฅโŽฅโŽฅโŽฅโŽฆ,(48) where โˆ’1<๐œ€<1 is a random number. (1)๐‘‡=0.2, ๐‘ƒ=50.(2)๐‘‡=10, ๐‘ƒ=50.When the control interval is short, as seen from Figures 2 and 3, the maximum overshot of fast channel is too big and the response speed of slow channel is too slow. The reason is that the predictive horizon is not long enough, the prediction model is not fully used, and the feedback correction plays an important role in this set point tracking process. When the control interval is long, as seen from Figures 4 and 5, the slow channel shows good control quality, but the respond speed of fast channel is still a little slow. A long control interval means a low control frequency and the fast channel achieves its set point by several control steps, so the fast channel response is a little slow. And the fast channel can achieve a temporary stable state, so the response curve may have a stair shape. The compromise methods are not so perfect. Therefore, we designed a decentralized DMC controller, and let ๐‘‡๐‘“=0.2, ๐‘‡๐‘ =20, ๐‘›=5, ๐‘ƒ๐‘ =50, and ๐‘ƒ๐‘“=50. We used this decentralized controller to track step signals. As shown in Figures 6 and 7, both fast channel and slow channel showed very good control quality. This method combined fast dynamic information and stable state information. The model information was fully used. From the simulation, we found that the decentralized method showed better control quality than the two compromised method.

A Continuously Stirred Tank Reactor (CSTR) Model
Chen studied the nonlinearity of a CSTR and modeled the CSTR by the following nonlinear equations [20]: ๐‘‘๐ถ๐ด๐‘‘๐‘ก=ฬ‡โ€Œ๐‘‰๐‘‰๐‘…๎€ท๐ถ๐ด0โˆ’๐ถ๐ด๎€ธโˆ’๐‘˜1(๐‘‡)๐ถ๐ดโˆ’๐‘˜3(๐‘‡)๐ถ2๐ด,๐‘‘๐ถ๐ต๐‘‘๐‘ก=โˆ’ฬ‡โ€Œ๐‘‰๐‘‰๐‘…๐ถ๐ต+๐‘˜1(๐‘‡)๐ถ๐ดโˆ’๐‘˜2(๐‘‡)๐ถ๐ต,๐‘‘๐‘‡๐‘‘๐‘ก=ฬ‡โ€Œ๐‘‰๐‘‰๐‘…๎€ท๐‘‡0โˆ’๐‘‡๎€ธโˆ’1๐œŒ๐ถ๐‘๎€ท๐‘˜1(๐‘‡)๐ถ๐ดฮ”๐ป๐‘…๐ด๐ต+๐‘˜2(๐‘‡)๐ถ๐ตฮ”๐ป๐‘…๐ต๐ถ+๐‘˜3(๐‘‡)๐ถ2๐ดฮ”๐ป๐‘…๐ด๐ท๎€ธ+๐‘˜๐‘ค๐ด๐‘…๐œŒ๐ถ๐‘๐‘‰๐‘…๎€ท๐‘‡๐พโˆ’๐‘‡๎€ธ,๐‘‘๐‘‡๐พ๐‘‘๐‘ก=1๐‘š๐พ๐ถ๐‘๐พ๎€ท๐‘„๐พ+๐‘˜๐‘ค๐ด๐‘…๎€ท๐‘‡โˆ’๐‘‡๐พ๎€ธ๎€ธ.(49)
We chose ๐‘ข1=ฬ‡โ€Œ๐‘‰/๐‘‰๐‘… and ๐‘ข2=๐‘„๐พ as the manipulate variables, and ๐ถ๐ด and ๐‘‡๐พ as the controlled variables. Parameter values were given in Table 2.
We chose a steady state and identified the input-output model to design a model predictive controller.
The input-output model is
๎‚ธ๐ถ๐ด๐‘‡๐‘˜๎‚น=โŽกโŽขโŽขโŽฃ2.868๐‘ +42.850โˆ’3.633๐‘ +24.940.04559๐‘ +14.24โŽคโŽฅโŽฅโŽฆ๎‚ธ๐‘ข1๐‘ข2๎‚น.(50)
The respond speed of ๐‘‡๐‘˜ is a little faster than that of ๐ถ๐ด. And the sample frequency of the two kinds of sensors cannot be same due to the limits of the sensors. The temperature sensor can only be sampled in high frequency, while the concentration sensor can only be sampled in low frequency. In order to maintain high control frequency, soft-sensing methods are often used in standard MPC. Due to computational errors and other errors, the soft-sensing method is not a perfect way. The decentralized method presented in this paper can also deal with different sample frequencies in different channels. We compared 0.1 step tracing effect of the standard MPC and that of the decentralized method. Let the sample interval of the temperature sensor be 0.1โ€‰minute, and let the sample interval of the concentration sensor be 1โ€‰minute.
(1)Perfect soft-sensing, ๐‘‡=0.1โ€‰minute, ๐‘ƒ=50.(2)Without soft-sensing ๐‘‡=1 minute, ๐‘ƒ=50.(3)Decentralized MPC ๐‘‡๐‘“=0.1minute, ๐‘‡๐‘ =2minutes, ๐‘›=10.
For the systems with little dynamic differences in different channels, if the controlled variables can provide sufficient reliable information in high frequency, the standard MPC can be applied and good control quality can be achieved (Figures 8 and 9). If the controlled variables can only provide reliable information in low frequency, the tracking speed turns slow (Figures 10 and 11). Although the dynamic characteristics of the two channels have a few differences, the decentralized MPC can also achieve good control quality (Figures 12 and 13) because the decentralized method takes full advantage of the reliable information of different channels.

๐พ 1 1 ๐พ 1 2 ๐พ 2 1 ๐พ 2 2 ๐‘‡ 1 1 ๐‘‡ 1 2 ๐‘‡ 2 1 ๐‘‡ 2 2 ๐œ 1 1 ๐œ 1 2 ๐œ 2 1 ๐œ 2 2


๐ถ ๐ด 0 Feed concentration of species A5.1โ€‰mol/L
๐‘‡ 0 Feed temperature104.9โ€‰C
๐‘˜ 1 0 Collision factor for reaction 1: ๐‘˜ 1 ( ๐‘‡ ) = ๐‘˜ 1 0 ๐‘’ โˆ’ ๐ธ 1 / ๐‘‡ 1 . 2 8 7 ร— 1 0 1 2 h โˆ’ 1
๐‘˜ 2 0 Collision factor for reaction 2: ๐‘˜ 2 ( ๐‘‡ ) = ๐‘˜ 2 0 ๐‘’ โˆ’ ๐ธ 2 / ๐‘‡ 1 . 2 8 7 ร— 1 0 1 2 h โˆ’ 1
๐‘˜ 3 0 Collision factor for reaction 3: ๐‘˜ 3 ( ๐‘‡ ) = ๐‘˜ 3 0 ๐‘’ โˆ’ ๐ธ 3 / ๐‘‡ 9 . 0 4 3 ร— 1 0 9 ( m o l A ) โˆ’ 1 h โˆ’ 1
๐ธ 1 Normalized activation energy for reaction 1 โˆ’ 9758.3โ€‰K
๐ธ 2 Normalized activation energy for reaction 1 โˆ’ 9758.3โ€‰K
๐ธ 3 Normalized activation energy for reaction 1 โˆ’ 8560โ€‰K
ฮ” ๐ป ๐‘… ๐ด ๐ต Enthalpies of reaction 14.2โ€‰kj/mol A
ฮ” ๐ป ๐‘… ๐ด ๐ต Enthalpies of reaction 2 โˆ’ 11โ€‰kj/mol B
ฮ” ๐ป ๐‘… ๐ด ๐ต Enthalpies of reaction 3 โˆ’ 41.85โ€‰kj/mol A
๐‘˜ ๐‘ค Heat transfer coefficient for cooling jacket4.032โ€‰kj/(h m 2 K)
๐ด ๐‘… Surface of cooling jacket0.215โ€‰ m 2
๐‘‰ ๐‘… Reactor volume0.01โ€‰ m 3
๐‘š ๐พ Coolant mass5.0โ€‰kg
๐ถ ๐‘ ๐พ Heat capacity of coolant2.00โ€‰kj/(kg โ‹… K)
๐ถ ๐‘ Heat capacity3.01โ€‰kj/(kg โ‹… K)
๐œŒ Density0.9342โ€‰kg/L

๐ถ ๐ด ๐ถ ๐ต ๐‘‡ ๐‘‡ ๐‘˜ ๐‘ข 1 ๐‘ข 2

2.4308โ€‰mol/L1.0802โ€‰mol/L115.4559โ€‰C114.9944โ€‰C20โ€‰ m i n โˆ’ 1 โˆ’ 400โ€‰kj/h

6. Conclusion

In this article, we focused on a kind of special system and designed a decentralized model predictive controller for it. This kind of system has different dynamics in different channels and exhibits two-time scale. A centralized MPC controller cannot satisfy the fast and the slow channels simultaneously. We used singular perturbation method to get the fast and the slow model from the original system. In actual processes, input-output models that can be obtained easily by identification were usually used to describe the real system. We demonstrated the singular perturbation method applying in transfer function matrix. Then we presented a decentralized model predictive controller based on the fast and the slow model and provided a sufficient condition for the algorithm stability when ๐‘›=1. Finally, the decentralized model predictive control algorithm was applied in two examples by simulation, and the validity of the control algorithm was tested. The simulation results proved that the two-time scale MPC is superior to the traditional MPC when the system had two-time scale characteristic.

The algorithm is based on the idea of fully using the information of the system. For the systems with two-time scale characteristics, the fast and slow channels are controlled, respectively, in the decentralized algorithm. This algorithm makes best use of the transition information of the fast channels and the slow channels and reduces the computation burden, which provides short control interval and increases the response speed. For those systems without two-time scale characteristics, this algorithm also works well. MPC has intensively been applied in the industrial process. The two-time scale MPC algorithm which is presented in this paper extends the applying scopes of MPC.


The authors gratefully acknowledge the financial support of 863 Program of China (no. 2007AA041402), National Key Scientific and Technical Project of China (no. 2007BAF22B05) and National Science Foundation of China (no. 60804023).


  1. P. V. Kokotovic, R. E. O'Malley Jr., and P. Sannuti, โ€œSingular perturbations and order reduction in control theoryโ€”an overview,โ€ Automatica, vol. 12, no. 2, pp. 123โ€“132, 1976. View at: Publisher Site | Google Scholar
  2. R. G. Phillips, โ€œA two-stage design of linear feedback controls,โ€ IEEE Transactions on Automatic Control, vol. 25, no. 6, pp. 1220โ€“1223, 1980. View at: Publisher Site | Google Scholar
  3. A. Saberi and H. K. Khalil, โ€œTwo time-scale feedback design of nonlinear singularly perturbed system,โ€ in Proceedings of the American Control Conference (ACC '83), pp. 441โ€“442, San Francisco, Calif, USA, June 1983. View at: Google Scholar
  4. P. V. Kokotovic, โ€œSubsystems, time scales and multimodeling,โ€ Automatica, vol. 17, no. 6, pp. 789โ€“795, 1981. View at: Publisher Site | Google Scholar
  5. H. K. Khalil and A. Saberi, โ€œDecentralized stabilization of nonlinear interconnected systems using high gain feedback,โ€ IEEE Transactions on Automatic Control, vol. 27, no. 1, pp. 265โ€“268, 1982. View at: Publisher Site | Google Scholar
  6. C. K. Yi and W. L. Luyben, โ€œDesign and control of coupled reactor/column systemsโ€”part 1: a binary coupled reactor/rectifier system,โ€ Computers and Chemical Engineering, vol. 21, no. 1, pp. 25โ€“46, 1996. View at: Publisher Site | Google Scholar
  7. C. K. Yi and W. L. Luyben, โ€œDesign and control of coupled reactor/column systemsโ€”part 2: more complex coupled reactor/column systems,โ€ Computers and Chemical Engineering, vol. 21, no. 1, pp. 47โ€“67, 1996. View at: Publisher Site | Google Scholar
  8. C. K. Yi and W. L. Luyben, โ€œDesign and control of coupled reactor/column systemsโ€”part 3: a reactor/stripper with two columns and recycle,โ€ Computers and Chemical Engineering, vol. 21, no. 1, pp. 69โ€“86, 1996. View at: Publisher Site | Google Scholar
  9. M.-N. Contou-Carrere and P. Daoutidis, โ€œDynamic precompensation and output feedback control of integrated process networks,โ€ in Proceedings of the American Control Conference (ACC '04), vol. 4, pp. 2909โ€“2914, Boston, Mass, USA, June-July 2004. View at: Google Scholar
  10. A. Kumar and P. Daoutidis, โ€œNonlinear dynamics and control of process systems with recycle,โ€ Journal of Process Control, vol. 12, no. 4, pp. 475โ€“484, 2002. View at: Publisher Site | Google Scholar
  11. N. Vora and P. Daoutidis, โ€œNonlinear model reduction of chemical reaction systems,โ€ in Proceedings of the American Control Conference (ACC '99), vol. 3, pp. 1583โ€“1587, San Diego, Calif, USA, June 1999. View at: Publisher Site | Google Scholar
  12. A. Kumar, P. D. Christofides, and P. Daoutidis, โ€œSingular perturbation modeling of nonlinear processes with nonexplicit time-scale multiplicity,โ€ Chemical Engineering Science, vol. 53, no. 8, pp. 1491โ€“1504, 1998. View at: Publisher Site | Google Scholar
  13. J. M. Maciejowski, Predictive Control with Constraints, Prentice Hall, London, UK, 2002.
  14. K. L. Buescher and C. C. Baum, โ€œA two-time-scale approach to nonlinear model predictive control,โ€ in Proceedings of the American Control Conference (ACC '95), vol. 3, pp. 2250โ€“2256, Seattle, Wash, USA, June 1995. View at: Google Scholar
  15. H. K. Khalil, โ€œOutput feedback control of linear two-time-scale system,โ€ IEEE Transactions on Automatic Control, vol. 32, no. 9, pp. 784โ€“792, 1987. View at: Publisher Site | Google Scholar
  16. S. F. Wang and J. H. Du, Illustrative Examples and Problems of Automatic Control Theory, Tsinghua University Press, Beijing, China, 2002.
  17. Y. G. Xi, Introduction on Dynamic Large-Scale Systems, National Defense Industry Press, Beijing, China, 1988.
  18. R. S. Madriz and S. S. Sastry, โ€œInput-output description of linear systems with multiple time-scales,โ€ International Journal of Control, vol. 40, no. 4, pp. 699โ€“721, 1984. View at: Publisher Site | Google Scholar
  19. Y. G. Xi, Predictive Control, National Defense Industry Press, Beijing, China, 1993.
  20. H. Chen, A. Kremling, and F. Allgower, โ€œNonlinear predictive control of a benchmark CSTR,โ€ in Proceedings of the 3rd European Control Conference (ECC '95), pp. 3247โ€“3252, Rome, Italy, 1995. View at: Google Scholar

Copyright © 2009 Jian Niu et al. 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.

More related articles

678ย Views | 485ย Downloads | 5ย Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.