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

Authors: Jian Niu; Jun Zhao; Zuhua Xu; Jixin Qian
Journal: Journal of Automated Methods and Management in Chemistry
 (2009)

---

## Abstract

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.


---

## Body

## 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 matrixG(s) with two-time scale characteristic is presented below, where
(1)G(s)=[G11(s)G12(s)G21(s)G22(s)]G11(s)=a11b11s+1eτ11s,G12(s)=a12b12s+1eτ12s,G21(s)=a21b21s+1eτ21s,G22(s)=a22b22s+1eτ22s,
when a11,a12,a21,a22 are in the same order, b11,b12≪b21,b22,τ11,τ12,τ21,τ22 are in the same order.
Considering the system
(2)Y=G(s)U,
the response speed of output y1 is much faster than that of output y2. Then a traditional central controller cannot satisfy the demands of y1 and y2 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
(3)G′(s)=[G11′(s)G12′(s)G21′(s)G22′(s)],
where G11′(s)=a11/b11s+1, G12′(s)=a12/b12s+1, G21′(s)=a21/b21s+1, G22′(s)=a22/b22s+1, a11,a12,a21,a22 are in the same order, b11,b12≪b21,b22, can have this form(4)Ẋ=AX+BU,Y=CX,
where
(5)A=diag[A11,A12,A21,A22],B=[B11B12B21B22],C=[C11C12C21C22];(Aij,Bij,Cij) is a one-dimension state space form of Gij′(s):(6)∵b11,b12≪b21,b22∴|A11|,|A12|≫|A21|,|A22|.So we can rewrite (4) in this form:(7)[εẋ1εẋ2ẋ3ẋ4]=[εA110εA12A210A22][x1x2x3x4]+[εB11εB12B21B22][u1u2],[y1y2]=[C11C12C21C22][x1x2x3x4],
where  ε is a very small positive constant.This form can be regarded as(8)Ẋ=A̅11X+A̅12Z+B̅1U,εŻ=A̅21X+A̅22Z+B̅2U,Y=C̅1X+C̅2Z+D̅U,
where
(9)A̅11=[A2100A22],A̅12=[0000],A̅21=[0000],A̅22=[εA1100εA12],B̅1=[B2100B22],B̅2=[εB1100εB12],C̅1=[00C21C22],C̅2=[C11C1200],D̅=[0000].
And the system in this form was discussed in several literatures [9, 15, 17], too. We got the slow model(10)dXsdt=AsXs+BsU,Y=CsXs+DsU,
and the fast model(11)dZfdτ=AfZf+BfU,Y=CfZf+DfU,
where
(12)As=A̅11-A̅12A̅22-1A̅21,Bs=B̅1-A̅12A̅22-1B̅2,Cs=C1-C2A̅22-1A̅21,Ds=D̅-C2A̅22-1B̅2,Af=A̅22,Bf=B̅2,Cf=C̅2,Df=D̅.We denoted the transfer function of the slow modelGs(s) by(13)Gs(s)=Cs(sI-As)-1Bs+Ds=[a11a12G21′(s)G22′(s)]
and the transfer function of the fast model Gf(s) in τ time scale by(14)Gf(s)=Cf(sI-Af)-1Bf+Df=C̅2(sI-A̅22)-1B̅2+D̅=C̅2(sI-A̅22)-1B̅2=[εa11b11s+εεa12b12s+ε00],
and the fast model in t-time scale is(15)Gr(s)=Gf(εs)=[a11b11s+1a12b12s+100].Gf(s),Gr(s), and Gs(s) are all descriptions from different points of the real system G′(s). And they are all two-in-two-out systems but have less state variables than the real system G′(s). If we choose a long enough sample interval for G′(s) and  Gs(s), 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 G′(s) and Gr(s) 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 Gs(s) and  Gf(s), and also the transfer function of the original system G′(s):(16)lims→∞Gs(s)=Ds=lims→0Gf(s),(17)G′(s)=Gs(s)+Gf(εs)-Ds+O(ε).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 O(ε). Let(18)O(ε)=[00O21(ε)O22(ε)],
and let(19)Gt(s)=Gr(s)+O(ε)=[a11b11s+1a12b12s+1O21(ε)O22(ε)]=[G11′(s)G12′(s)O21(ε)O22(ε)].We had expressions of the fast modelGt(s) and the slow model Gs(s) about the model G′(s) 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 G(s) shows two-time scale characteristic, as an example. Let(20)X11(s)=G11(s)U1(s)=G11′(s)U1(s)eτ11s,X12(s)=G12(s)U2(s)=G12′(s)U2(s)eτ12s,X21(s)=G21(s)U1(s)=G21′(s)U1(s)eτ21s,X22(s)=G22(s)U2(s)=G22′(s)U2(s)eτ22s,HereGij′(s),i,j=1,2, denotes the transition process, and eτijs,i,j=1,2, denotes the delay time. When τ11,τ12,τ21,andτ22 are in the same order, Gij′(s), i,j=1,2, reflects the main dynamic characteristic. We can get the fast model in t-time scale(21)Gtd(s)=[G11′(s)eτ11sG12′(s)eτ12sO21(ε)O22(ε)].
and the slow model in t-time scale(22)Gsd(s)=[a11a12G21′(s)eτ21sG22′(s)eτ22s]Then we can design a decentralized controller based on characteristics of the fast modelGt(s)(Gtd(s)) and the slow model  Gs(s)(Gsd(s)).


## 2.1. Without Delay Terms

Based on the method mentioned in literature[16], system
(3)G′(s)=[G11′(s)G12′(s)G21′(s)G22′(s)],
where G11′(s)=a11/b11s+1, G12′(s)=a12/b12s+1, G21′(s)=a21/b21s+1, G22′(s)=a22/b22s+1, a11,a12,a21,a22 are in the same order, b11,b12≪b21,b22, can have this form(4)Ẋ=AX+BU,Y=CX,
where
(5)A=diag[A11,A12,A21,A22],B=[B11B12B21B22],C=[C11C12C21C22];(Aij,Bij,Cij) is a one-dimension state space form of Gij′(s):(6)∵b11,b12≪b21,b22∴|A11|,|A12|≫|A21|,|A22|.So we can rewrite (4) in this form:(7)[εẋ1εẋ2ẋ3ẋ4]=[εA110εA12A210A22][x1x2x3x4]+[εB11εB12B21B22][u1u2],[y1y2]=[C11C12C21C22][x1x2x3x4],
where  ε is a very small positive constant.This form can be regarded as(8)Ẋ=A̅11X+A̅12Z+B̅1U,εŻ=A̅21X+A̅22Z+B̅2U,Y=C̅1X+C̅2Z+D̅U,
where
(9)A̅11=[A2100A22],A̅12=[0000],A̅21=[0000],A̅22=[εA1100εA12],B̅1=[B2100B22],B̅2=[εB1100εB12],C̅1=[00C21C22],C̅2=[C11C1200],D̅=[0000].
And the system in this form was discussed in several literatures [9, 15, 17], too. We got the slow model(10)dXsdt=AsXs+BsU,Y=CsXs+DsU,
and the fast model(11)dZfdτ=AfZf+BfU,Y=CfZf+DfU,
where
(12)As=A̅11-A̅12A̅22-1A̅21,Bs=B̅1-A̅12A̅22-1B̅2,Cs=C1-C2A̅22-1A̅21,Ds=D̅-C2A̅22-1B̅2,Af=A̅22,Bf=B̅2,Cf=C̅2,Df=D̅.We denoted the transfer function of the slow modelGs(s) by(13)Gs(s)=Cs(sI-As)-1Bs+Ds=[a11a12G21′(s)G22′(s)]
and the transfer function of the fast model Gf(s) in τ time scale by(14)Gf(s)=Cf(sI-Af)-1Bf+Df=C̅2(sI-A̅22)-1B̅2+D̅=C̅2(sI-A̅22)-1B̅2=[εa11b11s+εεa12b12s+ε00],
and the fast model in t-time scale is(15)Gr(s)=Gf(εs)=[a11b11s+1a12b12s+100].Gf(s),Gr(s), and Gs(s) are all descriptions from different points of the real system G′(s). And they are all two-in-two-out systems but have less state variables than the real system G′(s). If we choose a long enough sample interval for G′(s) and  Gs(s), 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 G′(s) and Gr(s) 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 Gs(s) and  Gf(s), and also the transfer function of the original system G′(s):(16)lims→∞Gs(s)=Ds=lims→0Gf(s),(17)G′(s)=Gs(s)+Gf(εs)-Ds+O(ε).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 O(ε). Let(18)O(ε)=[00O21(ε)O22(ε)],
and let(19)Gt(s)=Gr(s)+O(ε)=[a11b11s+1a12b12s+1O21(ε)O22(ε)]=[G11′(s)G12′(s)O21(ε)O22(ε)].We had expressions of the fast modelGt(s) and the slow model Gs(s) about the model G′(s) 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 G(s) shows two-time scale characteristic, as an example. Let(20)X11(s)=G11(s)U1(s)=G11′(s)U1(s)eτ11s,X12(s)=G12(s)U2(s)=G12′(s)U2(s)eτ12s,X21(s)=G21(s)U1(s)=G21′(s)U1(s)eτ21s,X22(s)=G22(s)U2(s)=G22′(s)U2(s)eτ22s,HereGij′(s),i,j=1,2, denotes the transition process, and eτijs,i,j=1,2, denotes the delay time. When τ11,τ12,τ21,andτ22 are in the same order, Gij′(s), i,j=1,2, reflects the main dynamic characteristic. We can get the fast model in t-time scale(21)Gtd(s)=[G11′(s)eτ11sG12′(s)eτ12sO21(ε)O22(ε)].
and the slow model in t-time scale(22)Gsd(s)=[a11a12G21′(s)eτ21sG22′(s)eτ22s]Then we can design a decentralized controller based on characteristics of the fast modelGt(s)(Gtd(s)) and the slow model  Gs(s)(Gsd(s)).


## 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(G′(s)) 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 Ps is the prediction horizon, Ts is the sampling interval based on the slow model, Pf is the prediction horizon, Tf is the sampling interval based on fast model, and M is the manipulate horizon. ai,j(t) is the step response of yi from uj at Tf sample interval. We got the model vector  ai,j=[ai,j(1)⋯ai,j(N)]T,  i=1,2, j=1,2, and N 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 modelGs(s).
On the t-time scale, the fast dynamic achieved a steady state. The model can be fully expressed by information not so necessary asai,j. Let  as,j=[as,j(1)⋯as,j(Ps)]T,j=1,2, be the slow model vector, where  as,j(i)=[a2,j((i-1)*(Ts/Tf)+1)⋯a2,j((i-1)*(Ts/Tf)+M)], i=1,…,Ps. The slow predict model is
(23)ỹs,PM(p)=ỹs,P0(p)+AsΔu,j=1,2,
where, ỹs,PM(p)∈RPs is the predict value of future output and  ỹs,P0(p)∈RPs is the prime predict value of future output:
(24)As=[as,j(1)0⋮⋱as,j(M)⋯as,j(1)⋮⋮as,j(Ps)⋯as,j(Ps-M+1)].Step 2.
Feedback Correction based on slow modelGs(s).
Let error vector be(25)es(p+1)=y2(p+1)-ỹs,PM(p+1p),
and we can get
(26)ỹs,cor(p+1)=ỹs,PM(p)+Hses(p+1),
where
(27)Hs=[hs(1)⋮hs(Ps)],
Like in Step 4, at p+1 time point, the time origin changes from p to p+1 time point, then the elements of vector ỹr,cor(p+1) should be moved, and the operation can be expressed by
(28)ỹs,N0(p+1)=Ss,0ỹs,cor(p+1),
where
(29)Ss,0=[010⋮⋱⋱⋮010⋯01].Step 3.
Model prediction based on fast modelGt(s).
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 ai,j. Let a1,j=[a1,j(1)⋯a1,j(Nr)]T, j=1,2. The fast predict model is
(30)ỹr,PM(k)=ỹr,P0(k)+ArΔu,j=1,2,
where
(31)Ar=[ar,j(1)0⋮⋱ar,j(M)⋯ar,j(1)⋮⋮ar,j(Pf)⋯ar,j(Pf-M+1)].Step 4.
Feedback Correction based on fast modelGt(s).
The inputu(k) was applied to the plant G′(s) at each k time point, and (22) gave the predictive output ỹr,PM(k+1∣k). Let error vector be
(32)er(k+1)=y1(k+1)-ỹr,PM(k+1∣k),
where y1(k+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:
(33)ỹr,cor(k+1)=ỹr,PM(k)+Hrer(k+1),
where  Hr=[hr(1)⋮hr(Pf)]. At k+1 time point, the time origin changes from k to k+1 time point, then the elements of vector ỹr,cor(k+1) should be moved, and the operation can be expressed by
(34)ỹr,P0(k+1)=Sr,0ỹr,cor(k+1),
where
(35)Sr,0=[010⋮⋱⋱⋮010⋯01].Step 5.
Rolling horizon optimization based on fast modelGt(s).
Let objective function be(36)minJ(k)=∥ω1(k)-ỹr,PM(k)∥Q12+∥ω2(p)-ỹs,PM(p)∥Q22+∥Δu(k)∥R2s.t.Δumin≤Δu≤Δumax,umin≤u≤umax,
where ωi(k), i=1,2, are the reference value:
(37)Q1=diag[q1(1)⋯q1(Pf)],Q2=diag[q2(1)⋯q2(Ps)],R=block-diag(R1,R2),Ri=diag[ri(1)⋯ri(Mf)],i=1,2.
Without constraint, we can get the manipulated variables:
(38)Δu(k)=L{[ArTQ1(ω1(k)-ỹr,P0(k))+AsTQ2(ω2(p)-ỹs,P0(p))](ArTQ1Ar+AsTQ2As+R)-1×[ArTQ1(ω1(k)-ỹr,P0(k))+AsTQ2(ω2(p)-ỹs,P0(p))][ArTQ1(ω1(k)-ỹr,P0(k))+AsTQ2(ω2(p)-ỹs,P0(p))](ArTQ1Ar+AsTQ2As+R)1},
where
(39)L=[10⋯0︸M0010⋯0︸M].
Step 3 to step 5 form the inner circulate, and in every  Tf 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  nTf (n is a positive integer) time interval. So the framework of the control system was showed as shown in Figure 1.Figure 1
Framework of the control system.


## 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 parametersTf,  Ts,  and Ps and one inner circulate is stable under the sufficient condition of a standard MPC, which has a control interval Tf and a predictive horizon Ps(Ts/Tf).Proof.
Let the prediction horizonPs(Ts/Tf), be the manipulate horizon, let  M be the objective function of the standard MPC be
(40)minJ′(k)=∥ω(k)-ỹPM(k)∥Q2+∥Δu(k)∥R2s.t.Δumin≤Δu≤Δumax,umin≤u≤umax,(41)J′(k)=∥ω(k)-ỹPM(k)∥Q2+∥Δu(k)∥R2=∥ω1(k)-ỹ1,PM(k)∥Q1′2+∥ω2(p)-ỹ2,PM(p)∥Q2′2+∥Δu(k)∥R2=∑i=1Ps(Ts/Tf)q1′(ω1(k)-ỹ1,PM(i+k∣k))2+∑i=1Ps(Ts/Tf)q2′(ω2(k)-ỹ1,PM(i+k∣k))2+∑i=1MrΔu(k+i∣k)2.
Because the response speed of y1 is very fast, the predictive value ỹ1,PM(i+k∣k) can be a fixed number when i is larger than a certain number (Pf):
(42)∑i=1Ps(Ts/Tf)q1′(ω1(k)-ỹ1,PM(i+k∣k))2=∑i=1Pfq1′(ω1(k)-ỹ1,PM(i+k∣k))2+Const.
The response speed ofy2 is very slow, so between a short time interval (Ts) the predictive values ỹ2,PM(i+k∣k), C≤i≤C+(Ts/Tf) (C is a positive number), can be linear correlation:
(43)∑i=CC+(Ts/Tf)q2′(ω2(k)-ỹ1,PM(i+k∣k))2=q2(ω2(k)-ỹ1,PM(C+k∣k))2,(44)∑i=1Ps(Ts/Tf)q2′(ω2(k)-ỹ1,PM(i+k∣k))2=∑j=1Psq2(ω2(k)-ỹ1,PM(j+k∣k)).
Combining (38), (40), (42), and (44), we got
(45)J′(k)=J(k)+const.
SominJ(k) and minJ′(k) had same answers, and the sufficient condition was proved.When the inner circulate is large than one, that is,n>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  n, 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, andu1 and u2 are the flow rates. The liquid level y1 and the temperature y2 are two controlled variables. y10=50cm and y20=295K are the initial stable states. The linear model of the system is
(46)[y1y2]=G(s)[u1u2],
where
(47)G(s)=[K11T11s+1e-τ11sK12T12s+1e-τ12sK21T21s+1e-τ21sK22T22s+1e-τ22s],
The response of outputy1 is much faster than that of  y2. 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 y1 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
(48)P(s)=[K11+0.2K11εe-(τ11+0.2τ11ε)s(T11+0.2T11ε)s+1K12+0.2K12εe-(τ12+0.2τ12ε)s(T12+0.2T12ε)s+1K21+0.2K21εe-(τ21+0.2τ21ε)s(T21+0.2T21ε)s+1K22+0.2K22εe-(τ22+0.2τ22ε)s(T22+0.2T22ε)s+1],
where -1<ε<1 is a random number. (1)
T=0.2, P=50.(2)
T=10, P=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  Tf=0.2, Ts=20, n=5, Ps=50, and Pf=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.Figure 2
Response ofy1.Figure 3
Response ofy2.Figure 4
Response ofy1.Figure 5
Response ofy2.Figure 6
Response ofy1.Figure 7
Response ofy2.A Continuously Stirred Tank Reactor (CSTR) Model
Chen studied the nonlinearity of a CSTR and modeled the CSTR by the following nonlinear equations [20]:
(49)dCAdt=V̇VR(CA0-CA)-k1(T)CA-k3(T)CA2,dCBdt=-V̇VRCB+k1(T)CA-k2(T)CB,dTdt=V̇VR(T0-T)-1ρCp(k1(T)CAΔHRAB+k2(T)CBΔHRBC+k3(T)CA2ΔHRAD)+kwARρCpVR(TK-T),dTKdt=1mKCpK(QK+kwAR(T-TK)).
We choseu1=V̇/VR and u2=QK as the manipulate variables, and CA and TK 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
(50)[CATk]=[2.868s+42.850-3.633s+24.940.04559s+14.24][u1u2].
The respond speed ofTk is a little faster than that of CA. 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,T=0.1 minute, P=50.(2)
Without soft-sensingT=1 minute, P=50.(3)
Decentralized MPCTf=0.1minute, Ts=2minutes, n=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 (Figures8 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.Table 1
Values of parameters.
K11K12K21K22T11T12T21T22τ11τ12τ21τ2221.11.31.51210009000.50.7108Table 2
Parameter values.
VariableDefinitionValueCA0Feed concentration of species A5.1 mol/LT0Feed temperature104.9 Ck10Collision factor for reaction 1:k1(T)=k10e-E1/T1.287×1012h-1k20Collision factor for reaction 2:k2(T)=k20e-E2/T1.287×1012h-1k30Collision factor for reaction 3:k3(T)=k30e-E3/T9.043×109(molA)-1h-1E1Normalized activation energy for reaction 1-9758.3 KE2Normalized activation energy for reaction 1-9758.3 KE3Normalized activation energy for reaction 1-8560 KΔHRABEnthalpies of reaction 14.2 kj/mol AΔHRABEnthalpies of reaction 2-11 kj/mol BΔHRABEnthalpies of reaction 3-41.85 kj/mol AkwHeat transfer coefficient for cooling jacket4.032 kj/(hm2 K)ARSurface of cooling jacket0.215m2VRReactor volume0.01m3mKCoolant mass5.0 kgCpKHeat capacity of coolant2.00 kj/(kg·K)CpHeat capacity3.01 kj/(kg·K)ρDensity0.9342 kg/LTable 3
Stable states.
CACBTTku1u22.4308 mol/L1.0802 mol/L115.4559 C114.9944 C20min-1-400 kj/hFigure 8
Response ofCA.Figure 9
Response ofTk.Figure 10
Response ofCA.Figure 11
Response ofTk.Figure 12
Response ofCA.Figure 13
Response ofTk.


## 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 whenn=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.



---

*Source: 164568-2009-09-30.xml*