# Design ofH2/H∞ RMPC for Boiler Superheated Steam Temperature Based on Memoryless Feedback Multistep Strategy

Authors: Pu Han; Miao Liu; Dongfeng Wang; Hao Jia
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

The collection of superheated steam temperature models of a thermal power plant under different loads can be approximated to “multimodel” linear uncertain systems. After transformation, the tracking system was obtained from “multimodel” linear uncertain systems. For this tracking uncertain system, a mixedH2/H∞ robust model predictive control (HRMPC) based on a memoryless feedback multistep strategy is proposed. A multistep control strategy combines the advantages of predictive control rolling optimization with memoryless feedback control thoughts. It could effectively decrease the controller optimization parameter and ensure closed-loop system stability, and, at the same time, it also achieved acceptable control performance. Successful application to the superheated steam temperature system of a 300 MW thermal power plant verified the study of the HRMPC-P cascade controller design scheme in terms of feasibility and effectiveness.


---

## Body

## 1. Introduction

The superheated steam temperature of a thermal power plant directly affects the thermal efficiency, safe operation of the superheated pipe, and the steam turbine equipment. A superheated steam system is characterized by large inertia and a substantial amount of lag [1]. Most power plants continue to adopt a conventional steam temperature control system, such as a cascade PID scheme or double circuit control system with a lead steam temperature differential signal. In recent years, a number of advanced control algorithms were being researched for application in superheated steam temperature systems; these include adaptive fuzzy neural network control [2], adaptive predictive control [3], nonlinear generalized predictive control [4], network cascade control [5], and active disturbance rejection control [6].Model predictive control (MPC), also known as rolling time control, was first applied to the linear time-invariant (LTI) system [7, 8]. Despite the heavy computational burden incurred in its implementation, MPC is characterized by good robustness, which has been extended to some other important classes of systems, such as distributed systems [9] and nonlinear systems [10]. Its application in industrial contexts is more recent [11]. Uncertainty and disturbance are widespread for the actual system. Using the method of robust control treatment for an uncertain system as reference, robust model predictive control (RMPC) can effectively deal with model uncertainties and disturbance. Furthermore, the controlled system achieves asymptotic stability under the condition of meeting the feasibility. Therefore, as an important branch of model predictive control, RMPC attracts increasing attention [12–14].Linear matrix inequalities (LMI), just as the name implies, are linear in the matrix variables. Many control problems can be converted into a feasibility problem of an LMI system or a convex optimization problem with LMI constraints [15]. With the development of the LMI toolbox in MATLAB, the LMI approach has been widely used in the field of systems and robust control [16].Lately, many robust model predictive control algorithms have been proposed [13–17]. Although these algorithms enable the system to approach a stable state, they do not meet the required expectations. Robust tracking controls were investigated for uncertain discrete time systems and application in thermal power plant systems [18–23]. However, the designs were too complicated and some designs had to undergo online optimization control as independent optimization variables to guarantee the constraint conditions for the control law.In this paper, we designed a mixedH2/H∞ robust model predictive control (HRMPC) based on a memoryless feedback multistep strategy for superheated steam temperature “multimodel” or “multipacket” linear uncertain system of a thermal power plant. Section 2 states the problem to be solved and introduces some standard assumptions. Section 3 discusses the HRMPC method, which adopts a closed-loop multistep control strategy and the memoryless feedback control thought to ensure the stability of the closed loop. At the same time, it could also achieve good control performance. In addition, considering simplified design, robust stability, and practical application, we studied a HRMPC-P cascade controller design scheme with cascade control structure in Section 4. Simulation and experimental results are discussed in Section 5. Lastly, concluding remarks are presented in Section 6.Preliminaries are as follows:(1)
A∗SB=ASTSB.(2)
Schur Complement Lemma: consider the partitioned matrix(1)A=A11A12A21A22.
(a)
WhenA11 is nonsingular, A22-A21A11-1A12 is called the Schur complement of A11 in A.(b)
WhenA22 is nonsingular, A11-A12A22-1A21 is called the Schur complement of A22 in A.


## 2. Problem Statement

Consider the following multipacket uncertain systems:(2)xk+1=Akxk+Bkuk+Bwωk,yk=Cxk,where x(k)∈Rn is the real-time state, u(k)∈Rm is the control input, y(k)∈Rl is the plant output, and ω(k)∈Rn is the unknown but bounded disturbance. A(k)∈Rn×n, B(k)∈Rn×m, and(3)Ak,Bk∈Ω=Ak,Bk∣Ak,Bk=∑i=1npλiAi,Bi,λi≥0,∑i=1npλi=1.A
i
,
B
i represent a state model of the system under a particular operating condition.In this paper, we consider Euclidean norm bounds and component-wise peak bounds on the inputuk+i∣k, given, respectively, as (4)uk+i∣k2≤umax,k,i≥0,(5)ujk+i∣k≤uj,max,k,i≥0,j=1,2,…,mwith a known bound umax>0.Assume thatA(k),B(k) are controllable, and the system external disturbance satisfying (6) is said to be an admissible disturbance(6)∑k=0∞wTkwk≤ϖwith a known upper bound ϖ>0.The RMPC system reference signalyr(k) is produced by the following system:(7)xrk+1=Arxrk+Brrk,yrk=Crxrk,where xr(k) are the reference states, r(k) are the expected input values, Ar is the Hurwitz matrix with an appropriate dimension, and Br, Cr are a constant matrices with an appropriate dimension.In order to ensure that the state variables of system (2) track the RMPC system reference signal yr(k), we designed a linear memoryless feedback controller as follows:(8)xck+1=xck+xrk-xk,uk=K1xk+K2xck,where xc(k)∈Rn is the controller state and K1, K2 are feedback controller gains that need to be solved.Substituting the memoryless feedback controller (8) into system (2), we have the close-loop system(9)xk+1=Ak+BkK1xk+BkK2xck+Bwωk,yk=Ckxk.By completing formulas (8) and (9), the tracking system can be obtained(10)x-k+1=A-k+B-kKx-k+B-wwk+Hxrk,yk=C-kx-k,where (11)x-k=xTkxcTkT,A-k=Ak0-II,B-wk=BwT0T,B-k=BTk0T,K=K1K2,H=0IT,C-=C0,A-k,B-k∈Ω-=A-k,B-k∣A-k,B-k=∑i=1npλiAik0-II,Bik0=∑i=1npλiA-ik,B-ik,λi≥0,∑i=1npλi=1.Therefore, the tracking system (10) is also a multipacket uncertain system. And the memoryless state feedback control law (8) can be rewritten as(12)uk=Kx-k.Control law (12) is required to make the tracking system (10) meet the following performance index.(a)
H
∞
Performance Index. Given the multiple packet uncertain time-varying discrete system (2) and a positive scalar γ>0, design a memoryless state feedback control law u(k) in the form of (12) such that the closed-loop tracking system (10) is asymptotically stable and satisfies (13)Twy∞≤γ,where Twy is the transfer function from ω(k) to y(k).(b)
H
2
Performance Index. Given the multiple packet uncertain time-varying discrete system (2) and a scalar α>0, design a memoryless state feedback control law u(k) in the form of (12) such that the output y(k) of the tracking system (10) is asymptotically stable and satisfies (14)y2≡∑k=0∞yTkyk≤α.


## 3. Design ofH2/H∞ Multistep Robust Predictive Controller

### 3.1.H2/H∞ Performance Index

In order to reduce the conservativeness of the design, the multistep control setS0=x∣xTQ0-1x≤1 is utilized [11], and the control strategy (15) as the system control law in the future(15)ε=K01,…,KNu-11,K02,…,KNu-12,where Ki1, Ki2 are predicted feedback gains of the tracking system at sampling time k+i. Nu is control horizon.With respect to the input constraints, we can get(16)ul2=αlYiQi-1x-2=αlYiQi-1/2Qi-1/2x-2≤αlYiQi-1/2,l=1,…,m,where ul is the lth element of the control input corresponding to the feedback control gain of Ki1, Ki2 and αl is the lth row of m-dimensional identity matrix.For this control strategy, we select the Lyapunov function as follows:(17)Vx-k+i=x-Tk+iPix-k+i,i=0,…,N-1,where N is the prediction horizon and the length of N indicates N steps from time k, the predicted output value close to the expected value. Pi>0, when i>N, Pi=PN-1.Denote(18)βik=A-k+B-kK+L-C0.We can obtain that(19)Vx-k+i+1-Vx-k+i=x-k+i+1Pi+12-x-k+iPi2=x-Tk+iβiTkPi+1βikx-k+i+2x-Tk+iβiTkPi+1B-wwk+i+2x-Tk+iβiTkPi+1Hxrk+i+wTk+iB-wTPi+1B-wwk+i+2wTk+iB-wTPi+1Hxrk+i+xrTk+iHTPi+1Hxrk+i-x-Tk+iPix-k+i+yTk+iyk+i-yTk+iyk+i+γ2wTk+iwk+i-γ2wTk+iwk+i=x-Tk+ivTk+iWikx-Tk+ivTk+iT-yTk+iyk+i+γ2vTk+ivk+i.Denote(20)χ=βiTkPi+1βik-Pi+C-TC-;then(21)Wik=χβiTkPi+1B-wβiTkPi+1H∗B-wTPi+1B-wB-wTPi+1H∗∗HTPi+1H.Lemma 1.
For tracking system (10), denote Ki=Ki1H1+Ki2H2, where H1=I0, H2=0I, and then a memoryless state feedback control law (22)Ki=YiQi-1exists such that Twy∞≤γ if there exist γ>0, Qi∈R3n×3n and Yi∈R3m×3n satisfying(23)-Qi∗∗∗∗00∗∗∗00-αγ2I∗∗A-Qi+B-YiB-wH-Qi+1∗C-Qi000-αI≤0,(24)-1∗∗γ2ϖ-αγ2ϖ∗x-k0-Q0≤0,where i=0,…,N-1.Proof.
We sum both sides of equality (19) from i=0 to i=∞. Actually, when i>N-1 the system output is close to the expected value and limi→∞⁡V(x-(k+i))=0; then we obtain the following equation:(25)-x-TkPkx-k=∑i=0∞-yTk+iyk+i+γ2wTk+iwk+i+∑i=0∞x-Tk+iwTk+ixrTk+iWikx-Tk+iwTk+ixrTk+iT.Equation (25) can be rewritten as(26)∑i=0∞yTk+iyk+i=x-TkPkx-k+∑i=0∞γ2wTk+iwk+i+∑i=0∞x-Tk+iwTk+ixrTk+iWikx-Tk+iwTk+ixrTk+iT.
(1) Whilex(k)=0, the performance index (13) is equivalent to(27)∑i=0∞yTk+iyk+i-γ2wTk+iwk+i≤0.
Thus, from (26) it is known that if Wi(k)≤0, then (27) is satisfied. In other words, H∞ performance index is met.
(2) From (26), if Wi(k)≤0, satisfying the following conditions, then H2 performance index is met:(28)∑i=0∞yTk+iyk+i≤x-TkP0x-k+γ2ϖ≤α.
ForWi(k), (21) can be written as(29)Wik=χβiTkPi+1B-wβiTkPi+1H∗B-wTPi+1B-wB-wTPi+1H∗∗HTPi+1H=-Pi∗∗00∗00-γ2I+C-T00C-00+βiTkB-wTHTPi+1βikB-wH.
Further, using the Schur complement in (29), we have(30)-Pi∗∗∗∗00∗∗∗00-γ2I∗∗βikB-wH-Pi+1-1∗C-000-I≤0.
Multiply both sides of inequality (30) by using (31)diagα1/2Pi-1,α1/2I,α1/2I,α1/2I,α1/2I.Denote that Qi=αPi-1, Ki=YiQi-1; condition (23) can be obtained.
In the same way, using the Schur complement in (28) x-T(k)P0x-(k)+γ2ϖ≤α, we have(32)-α∗∗γ2ϖ-γ2ϖ∗x-k0-P0-1≤0.
Multiplying both sides of inequality (32) by using(33)diagα-1/2I,α1/2I,α1/2Ithen condition (24) can be obtained.


### 3.2. Constraint Condition

For constraint (5), we have the following lemma.Lemma 2.
If there existQi∈R3n×3n and Yi∈R3m×3n, i=0,…,N-1 such that Lemma 1 is established, and there exist symmetric matrices X such that the following inequality is satisfied, then constraint (5) is met(34)-X∗YT-Q≤0,Xjj≤uj,max2,j=1,2,…,nm.Proof.
At sample timek, consider the Euclidean norm constraint (5).(35)uk+i∣k2≤umax,k,i≥0.
Following controller (12), we have(36)maxi≥0uk+i∣k22=maxi≥0YQ-1x-k+i∣k22≤maxz∈εYQ-1z22=λmaxQ-1/2YTYQ-1/2,where ε=z∣zTQ-1z≤1=z∣zTPz≤γ; using the Schur complement, we have(37)uk+i∣k22≤umax2,i≥0.If the inequality is established(38)-umax2I∗YT-Q≤0at sample time k(39)ujk+i∣k≤uj,max,i≥0,j=1,2,…,musing the Cauchy-Schwarz inequality [24], we have(40)maxi≥0ujk+i∣k2=maxi≥0YQ-1xk+i∣kj2≤maxz∈εYQ-1zj2≤YQ-1/2j22=YQ-1YTjj.
Thus, symmetric matrixX satisfies(41)-X∗YT-Q≤0,Xjj≤uj,max2,j=1,2,…,nm,where Xjj are diagonal elements of matrix X.
In conclusion, for tracking system (10), we have the memoryless feedback multistep H2/H∞ robust model predictive control algorithm from Lemmas 1 and 2. At sample time k, solve the following optimization problem:(42)minQi,Yi,Xi,ααSubjectto23,24,34.
The current system input of the controller is(43)uk=Kx-k.


### 3.3. Robust Stability

Theorem 3.
Consider the multipacket uncertain tracking system (14). If the optimization problem is feasible with the system state x-(k), then system (14) in close-loop is robustly stable.Proof.
At sample timek, optimization problem (42) is feasible with the optimal solution(44)Q0,…,QN-1,Y0,…,YN-1,X0,…,XN-1,γk.
While disturbance exists at sample timek+1, the optimization problem (42) has the optimal solution(45)Q1,…,QN-1,QN-1,Y1,…,YN-1,YN-1,X1,…,XN-1,XN-1,γ∗k,where γ(k+1)=γ∗(k), and by the principle of optimization, we have(46)γ∗k≥γk+1≥γ∗k+1.
When disturbance disappears, due to the fact thatWi(k)≤0, V(k+1)<V(k), it means that V(k) will decay to 0. Thus, the closed-loop system has robust stability.


## 3.1.H2/H∞ Performance Index

In order to reduce the conservativeness of the design, the multistep control setS0=x∣xTQ0-1x≤1 is utilized [11], and the control strategy (15) as the system control law in the future(15)ε=K01,…,KNu-11,K02,…,KNu-12,where Ki1, Ki2 are predicted feedback gains of the tracking system at sampling time k+i. Nu is control horizon.With respect to the input constraints, we can get(16)ul2=αlYiQi-1x-2=αlYiQi-1/2Qi-1/2x-2≤αlYiQi-1/2,l=1,…,m,where ul is the lth element of the control input corresponding to the feedback control gain of Ki1, Ki2 and αl is the lth row of m-dimensional identity matrix.For this control strategy, we select the Lyapunov function as follows:(17)Vx-k+i=x-Tk+iPix-k+i,i=0,…,N-1,where N is the prediction horizon and the length of N indicates N steps from time k, the predicted output value close to the expected value. Pi>0, when i>N, Pi=PN-1.Denote(18)βik=A-k+B-kK+L-C0.We can obtain that(19)Vx-k+i+1-Vx-k+i=x-k+i+1Pi+12-x-k+iPi2=x-Tk+iβiTkPi+1βikx-k+i+2x-Tk+iβiTkPi+1B-wwk+i+2x-Tk+iβiTkPi+1Hxrk+i+wTk+iB-wTPi+1B-wwk+i+2wTk+iB-wTPi+1Hxrk+i+xrTk+iHTPi+1Hxrk+i-x-Tk+iPix-k+i+yTk+iyk+i-yTk+iyk+i+γ2wTk+iwk+i-γ2wTk+iwk+i=x-Tk+ivTk+iWikx-Tk+ivTk+iT-yTk+iyk+i+γ2vTk+ivk+i.Denote(20)χ=βiTkPi+1βik-Pi+C-TC-;then(21)Wik=χβiTkPi+1B-wβiTkPi+1H∗B-wTPi+1B-wB-wTPi+1H∗∗HTPi+1H.Lemma 1.
For tracking system (10), denote Ki=Ki1H1+Ki2H2, where H1=I0, H2=0I, and then a memoryless state feedback control law (22)Ki=YiQi-1exists such that Twy∞≤γ if there exist γ>0, Qi∈R3n×3n and Yi∈R3m×3n satisfying(23)-Qi∗∗∗∗00∗∗∗00-αγ2I∗∗A-Qi+B-YiB-wH-Qi+1∗C-Qi000-αI≤0,(24)-1∗∗γ2ϖ-αγ2ϖ∗x-k0-Q0≤0,where i=0,…,N-1.Proof.
We sum both sides of equality (19) from i=0 to i=∞. Actually, when i>N-1 the system output is close to the expected value and limi→∞⁡V(x-(k+i))=0; then we obtain the following equation:(25)-x-TkPkx-k=∑i=0∞-yTk+iyk+i+γ2wTk+iwk+i+∑i=0∞x-Tk+iwTk+ixrTk+iWikx-Tk+iwTk+ixrTk+iT.Equation (25) can be rewritten as(26)∑i=0∞yTk+iyk+i=x-TkPkx-k+∑i=0∞γ2wTk+iwk+i+∑i=0∞x-Tk+iwTk+ixrTk+iWikx-Tk+iwTk+ixrTk+iT.
(1) Whilex(k)=0, the performance index (13) is equivalent to(27)∑i=0∞yTk+iyk+i-γ2wTk+iwk+i≤0.
Thus, from (26) it is known that if Wi(k)≤0, then (27) is satisfied. In other words, H∞ performance index is met.
(2) From (26), if Wi(k)≤0, satisfying the following conditions, then H2 performance index is met:(28)∑i=0∞yTk+iyk+i≤x-TkP0x-k+γ2ϖ≤α.
ForWi(k), (21) can be written as(29)Wik=χβiTkPi+1B-wβiTkPi+1H∗B-wTPi+1B-wB-wTPi+1H∗∗HTPi+1H=-Pi∗∗00∗00-γ2I+C-T00C-00+βiTkB-wTHTPi+1βikB-wH.
Further, using the Schur complement in (29), we have(30)-Pi∗∗∗∗00∗∗∗00-γ2I∗∗βikB-wH-Pi+1-1∗C-000-I≤0.
Multiply both sides of inequality (30) by using (31)diagα1/2Pi-1,α1/2I,α1/2I,α1/2I,α1/2I.Denote that Qi=αPi-1, Ki=YiQi-1; condition (23) can be obtained.
In the same way, using the Schur complement in (28) x-T(k)P0x-(k)+γ2ϖ≤α, we have(32)-α∗∗γ2ϖ-γ2ϖ∗x-k0-P0-1≤0.
Multiplying both sides of inequality (32) by using(33)diagα-1/2I,α1/2I,α1/2Ithen condition (24) can be obtained.


## 3.2. Constraint Condition

For constraint (5), we have the following lemma.Lemma 2.
If there existQi∈R3n×3n and Yi∈R3m×3n, i=0,…,N-1 such that Lemma 1 is established, and there exist symmetric matrices X such that the following inequality is satisfied, then constraint (5) is met(34)-X∗YT-Q≤0,Xjj≤uj,max2,j=1,2,…,nm.Proof.
At sample timek, consider the Euclidean norm constraint (5).(35)uk+i∣k2≤umax,k,i≥0.
Following controller (12), we have(36)maxi≥0uk+i∣k22=maxi≥0YQ-1x-k+i∣k22≤maxz∈εYQ-1z22=λmaxQ-1/2YTYQ-1/2,where ε=z∣zTQ-1z≤1=z∣zTPz≤γ; using the Schur complement, we have(37)uk+i∣k22≤umax2,i≥0.If the inequality is established(38)-umax2I∗YT-Q≤0at sample time k(39)ujk+i∣k≤uj,max,i≥0,j=1,2,…,musing the Cauchy-Schwarz inequality [24], we have(40)maxi≥0ujk+i∣k2=maxi≥0YQ-1xk+i∣kj2≤maxz∈εYQ-1zj2≤YQ-1/2j22=YQ-1YTjj.
Thus, symmetric matrixX satisfies(41)-X∗YT-Q≤0,Xjj≤uj,max2,j=1,2,…,nm,where Xjj are diagonal elements of matrix X.
In conclusion, for tracking system (10), we have the memoryless feedback multistep H2/H∞ robust model predictive control algorithm from Lemmas 1 and 2. At sample time k, solve the following optimization problem:(42)minQi,Yi,Xi,ααSubjectto23,24,34.
The current system input of the controller is(43)uk=Kx-k.


## 3.3. Robust Stability

Theorem 3.
Consider the multipacket uncertain tracking system (14). If the optimization problem is feasible with the system state x-(k), then system (14) in close-loop is robustly stable.Proof.
At sample timek, optimization problem (42) is feasible with the optimal solution(44)Q0,…,QN-1,Y0,…,YN-1,X0,…,XN-1,γk.
While disturbance exists at sample timek+1, the optimization problem (42) has the optimal solution(45)Q1,…,QN-1,QN-1,Y1,…,YN-1,YN-1,X1,…,XN-1,XN-1,γ∗k,where γ(k+1)=γ∗(k), and by the principle of optimization, we have(46)γ∗k≥γk+1≥γ∗k+1.
When disturbance disappears, due to the fact thatWi(k)≤0, V(k+1)<V(k), it means that V(k) will decay to 0. Thus, the closed-loop system has robust stability.


## 4. HRMPC-P Cascade Controller Design

The optimization goal of the inner loop is to quickly eliminate interference from the spraying system or burning system, and the optimization goal does not require difference. As a result, the proportional controller is often used in the inner loop, which is regarded as a quick follow-up system. In this paper, we study HRMPC-P controller design with a cascade control structure by considering a simplified design, robust stability, and practical application. A proportional controller is adopted for the inner loop. The HRMPC-P cascade control structure is as shown in Figure1.Figure 1
HRMPC-P control structure.Assume that the secondary superheated transfer functions are known, including both the leading and inertia segments. The HRMPC-P cascade control design steps are as follows:(1)
Select appropriate rolling optimization stepsN, control horizon Nu, and γ. The feasible regions of the controller broaden as Nu increases [11].(2)
Optimize the inner loop PID controller; in this work, the proportional controller simply needs to optimize the proportional gainδ.(3)
Set expectationsr and the reference system.(4)
At sample timek, when i=0, solve optimization problem (42), and obtain the optimal solution(47)Q0,…,QN-1,Y0,…,YN-1,X0,…,XN-1,γk.(5)
i=i+1,…, until i=N, and obtain the corresponding optimal solution(48)Q1,…,QN-1,QN-1,Y1,…,YN-1,YN-1,X1,…,XN-1,XN-1,γ∗k.(6)
Determine the current moment controller output(49)uk=Kkx-k.(7)
Takingu(k) into the inner loop discrete state space model (50), we have y2(k)(50)x2k+1=A2kx2k+B2kupk,y2k=Cx2k,wherex2(k) is the leading segment state, up(k) is the proportional control output, y2(k) is the leading segment output.(8)
Takingy2(k) into the tracking system (10), we have the state x-(k+1).(9)
Makek=k+1, jump to Step (4).


## 5. Simulation

Choose two secondary superheated system transfer functions under different loads of 180 MW and 250 MW of a 300 MW power unit [21]. Design the controller to use the method discussed in Sections 4 and 5.(51)G1s=-120s+11.225s+14C°/%,G2s=-124s+1135s+14C°/%.The first part is a leading segment transfer function and the last part is the inertia segment transfer function for each transfer function.Using a sampling timeT=0.95 s, convert the transfer function mode into discrete state space mode. After discretization, the leading segment and the inertia segment are as in Table 1.Table 1
Discrete state space model.
Leading segment
Inertia segment
A
2
B
2
A
B
180 MW
0.9536
0
0
0
0.0464
0.9536
0
0
0
0.0464
0.9536
0
0
0
0.0464
0.95362
-
0.0464
0
0
0
0.9627
0
0
0
0.0373
0.927
0
0
0
0.0373
0.9627
0
0
0
0.0373
0.9627
0.0447
0
0
0
250 MW
0.9612
0
0
0
0.0388
0.9612
0
0
0
0.0388
0.9612
0
0
0
0.0388
0.9612
-
0.0388
0
0
0
0.9732
0
0
0
0.0268
0.9732
0
0
0
0.0268
0.9732
0
0
0
0.0268
0.9732
0.0268
0
0
0In each stateω(k)22≤0.01, ϖ=0.5,(52)C=0001,Bw=0000.1T,Ar=0.8270000.1730.8270000.1730.8270000.1730.827,Br=0.173000T,Cr=0001.The step responses of the system with different values ofγ are shown in Figure 2. Choose the transfer function G1(s) as an example, set the value of γ=0.1 and γ=2, and we have the step responses of the system. As is shown in the figure, the ability to restrain the interference of the system is increased with the decrease of γ.Figure 2
Step responses of the system for different values ofγ.The tracking performance of different controllers is shown in Figure5. Choose the transfer function G1(s) as an example. We compared the control quality of HRMPC-P cascade control with the conventional cascade PID control (PID-P) and dynamic matrix control [6] (DMC-P) through the system output tracking situation. The control laws and control structures are illustrated as follows.(1) PID Controller
(53)
u
P
I
D
t
=
K
P
e
t
+
K
i
∫
e
t
d
t
+
K
d
d
e
t
d
t
,
e
t
=
y
t
-
y
r
t
,where KP, Ki, and Kd all nonnegative, denoting the coefficients for the proportional, integral, and derivative terms; e(t) is the difference between a desired set point and a measured process variable. The PID-P control structure is shown in Figure 3.Figure 3
PID-P control structure.(2) DMC Controller
(54)
u
D
M
C
k
=
u
D
M
C
k
-
1
+
Δ
u
D
M
C
k
d
i
,
Δ
u
D
M
C
k
=
G
T
Q
G
+
λ
-
1
G
T
Q
Y
r
k
+
1
-
G
0
U
k
-
1
-
h
e
k
,where G is dynamic matrix; G0 is process matrix which is composed of dynamic coefficients; Q and h are the weighting matrices. di means the ith row of ΔuDMC(k), if the current moment i=1. The DMC-P control structure is shown in Figure 4.Figure 4
DMC-P control structure.Figure 5
Tracking performance of different controllers.(3) HRMPC-P Controller. The control law and control structure refer to (12) and Figure 1.As is shown in Figure5, the HRMPC-P cascade control enables the system to have good tracking performance such that the system output can meet the expectations faster. The inputs of the different controllers are shown in Figure 6(a). The inputs of HRMPC controller, DMC controller, and PID controller, respectively, are system states, control increment, and the difference between desired set points and measured process variables. The outputs of the different controllers are shown in Figure 6(b). The HRMPC-P controller has a smaller output range and faster adjusting speed and can complete the adjustment in about 200 s.Figure 6
Input and output of different controllers.
(a)
Input of different controllers
(b)
Output of different controllersAs is shown in Figure7, the load command changes from 180 MW to 250 MW when t=1000 s. Under this condition, when t=1000 s the system model translates from G1(s) to G2(s). The robustness of different controllers is shown in Figure 8. HRMPC-P cascade control still has an enhanced tracking effect, which can quickly adjust the system to the expected value while the system changes. Thus, HRMPC-P cascade control structure has stronger robustness.Figure 7
Load command.Figure 8
Robustness of different controllers.Figure9 shows the output of different controllers while the load command changes. The HRMPC-P controller also has a smaller output range and faster adjusting speed. The variation tendency of HRMPC controller gains K1 and K2 is shown in Figure 10. The variation range of these two gain groups is small, and they adjust quickly as well.Figure 9
Output of different controllers.Figure 10
The HRMPC controller gain.


## 6. Concluding Remarks

This paper presents a design strategy for a memoryless feedback multistepH2/H∞ RMPC controller for the superheated steam temperature of a thermal power unit. In this study, the collection of superheated steam temperature model under different load is approximated to “multimodel” linear uncertain systems. After transformation, the tracking system is obtained. Then through deduction, we have the HRMPC control flow. Moreover, successful application of the proposed HRMPC-P cascade control scheme to the superheated steam temperature system of a thermal power plant and improved control qualities are achieved. The HRMPC-P cascade control scheme is able to adapt to the whole operating range of the superheated steam temperature system of thermal power plants. Moreover, the decrease in the number of controller parameters may serve as a reference for actual engineering projects.



---

*Source: 2192393-2017-04-18.xml*