W,+0)fXq:H@8Vol.7, No.7, 48/52J2008K"3~VlL==,)fKhk~AoNmP9H9Y)f†P 6 > ' &Y@ o &f e∗<∗∗∗∗Robust Temperature Control of a Reformer by Using Stable Continuous-Time Generalized PredictiveControl
Naohiko Ishibashi∗ , Mingcong Deng∗∗ and Akira Inoue∗∗

A temperature controller design using a stable continuous-time anti-windup generalized predictive control(CAGPC) system is applied to a reformer of a fuel cell system in this paper. The controlled process includes
uncertain time delay and input constraints. The controller is designed by using coprime factorization representation and Youla-Kucera parametrization to be robust to the uncertain time delay and the input constraints. The
effectiveness of the proposed design scheme is confirmed by a real application to the fuel cell system temperaturecontrol.
Key Words: stable continuous-time predictive control, anti-windup, robustness, reformer, fuel cell

1.

O8aKa/, Jq;&CO2 SPLo:N,W-+i3<8'M
l<7gs79F`NaZ,|T5lF*j, HmQ*hS.,OH3QN.,O,nK*$FO, GNb,RA3AE
SJPolymer Electrolyte Fuel Cell: PEFCK/E79F`
Na|BQ=,>^lF$k.
3AESOeGH@Gr=X?~5;er8.9kaxG/E9k?aK$/j<sJ(Mk.<r@k3H,G-k.
7+7$9~VN/Er@k?aKO"3*JeGN!k_w,,WHJj$=uGOB4-*hS$sUiNdj+i=,J,ON!k_wN_VO$qG"k. =N?a$lL
K.L7F$k:=eGO3AJTT,9$LPG$t}JIK
reGjCAJ~A,9K~A9k=X?~oJ~AoKr
V/3HG$9~VKo?k/Er@F$k. 33G$~A,9fNH.rBjK]D3H,Bj7?/E*hS PEFC
NQW-N]KR,k,$=N?aKO~AobK<67?
(^9Yr\89YKBj*K)f7$=X?~rBj5;
k,W,"k. ~AoN)fKO (1) /E+O0N)f$(2)/E+OeN)fN 2 o`"k,$:lb$~AobGN~AWm;9OTN+J`@~Vr}A$"/Ae(<?N0nOOKO)B,"k. D^j~O)srbC?`@~VW†IEEE International Conference on Industrial Technology(200612 )∗( )∗∗∗∗∗G/=!
/ n0)E)H t,3gX

Mitsubishi Heavy Industries, LTD.
Okayama University
Received January 4, 2008

JKm;9H@$9(k3H,G-k.
\@8GO$/E+O0N~O)sr<$TN+JWm;
9KP7$Bj"3~VlL==,)fJCAGPCK79F`r_W9k. CK$~AoK,Q9k?aK$>h$L9KM!5lF-?$1) `@~VNTN+5KmP9HJ)f
O$2) ~O)sKmP9HJ)fON2DNj!rH_go;k. sF9k_Wj!O{s,r==H Youla QiaHj<<7gs ’ rQ$$D.N2J,_Wj!rsF9k.
(1)`@~VNTN+5KmP9HBjJ)fOraak?
aK$~O)sr5k7?~A)fOKP7, 3sHm<
iOBjG"j , +DDk<WObBj$9JoA$/
BjJ=,)for_W9k.
(2)Dk<WC-K*$F, ~O)sKmP9HKBjJd~o ("sAo$I"CW3sHm<i ) HJk?aNr
oraa$(1) Gaa?)foKrorUC9k.
\j!GO$e- 1) NTN+J`@~VKP9kmP9H-KP7FO$`@~V`Naw rQ$$`@~VNTN+5Hawm9r)fONTN+5H7$TN+JWm;9KP9k/Bj"3~VlL==,)fN_Wj!rQ$k.
e- 2) N~O)sXNmP9H-O$~O)srTN+J
s~AWGH_J7$3sHm<i,/BjNH-TN+5
rbC?Dk<WOOmP9HBjG"k3H Khj$(1)G@?roKUC5lkroH7Faak.
3) 8)

5)

6)

4)

6)

2.

bGk=Hdj_j

3AESQ~Ao79F`NJ1J=$r Fig.1 K(9. ~Ao (Reformer) O(^r<67?(^X (Reforming CatTR 0007/08/707–0048 c 2008 SICE

W,+0)fXq:H@8 h 7 , h 7 f 2008 /49u(t) = σ(u1 (t))J2K)f~O u(t) OD.N)sru1kbNH9k.
σ(v) = umaxv

uminif v > umax ≥ 0if umin ≤ v ≤ umax

J3K

if v < umin ≤ 033G, u (t) O~O)sru1k0N.fG"j3sHm<iNPOr(9. Fig. 2 K*$F$w O2H~OG"
j$P O)fP]r=9. \*O Fig.1 G T r\89YK)f9k3sHm<i C r_W9k3HG"k.
111Fig. 1 The reformer sub-systemalyst) K$6AH7F:=eGO3AHer!k7$eGjC
AJ~A,9 (H2 Rich Reformed Gas) K~A9k. ~A,
9O3AES (Fuel Cell) NeGK&Kwil$@GKKw
il?u$H&K$/ELK~8F~A,9fNeG,Cq
5l$DC?eG,9O3Fo (Combustor) G3Fu$H
HbK3F5l$SP5lk. ~A,9H.fNl@=:G(CO) r 10[ppm] J<K)f7J1lP3AESN/EPO
Oc<9k?a$~A,9fNH.OBj7F)f5lF$
J1lPJiJ$.
~A,9NH.)fGO (1) /E+O0N)f+i$(2) /E+OeN)fKZj9(FTJ&. CK$(1) N/E+O0N)fGO$3AESK~A,9r!k9k0KO~Aob
N~A?~9Y,\89YKBj7F)f5lF$k,W,"j$/E+O0K*1k~Aob9YN\89YXN0j~V;LO$3AESN/E+O~VN;LKDJ,j$/0~K*1k?>z(N~erO+k3H,G-$EWG"k.
hCF$\@8GO/E+O0N~AoN9Y)fdjr
M!9k. 9JoA$~A,9O4L,3FoK!k5lF$
kuVrM(k. Fig.1 K(9H*j~AobN9Y T r$)fo C rp7F3A!kLKhCF)f7$e.LOl
jH9k. ~AobN~A?~O[.?~G"k?a$?~K,WJ.LO3FoGN3F.KhjEolk. 7?,C
F$C KhCF3A!kLr)f9k3HKhCF3F.,)f5l$~Ao9Y T r)f9k3H,G-k.
!3A!kLH~Ao9YVKO!0K(9H*j$TN+J`@~Vr<C? n !YlWm;9H>j9k.
C[s]B[s]e}U (s) +V (s) J1KY (s) = {A[s]A[s]A [s]33G, A[s], B[s] OBjJ?`0H9k. Y (s), U (s) *hS V (s) O=l>lPO (Fig.1 fN9Y T ), ~O (P<J
3A!kL), 0pr=9. Y (s)$U (s) O=l>l y(t), u(t)NiWi9XtG"k. A[s], B[s] O=l>l, n !, m !N{s?`0G"k. T OTN+J$NN`@~VG"k,,0 < T ≤ T ≤ T Ne<BM T , T O{NH9k. ^?, deg(C[s]) = n HJkh&K C[s] rjak . A [s] O`@~Vaw+iI87?`@~V`G"j, !aGaak.
33G$-fH7F$A[s] O s N?`0r(7 A(s) O s N-}Xtr(9.

BjGmP9HJ3sHm<iN_Wj!

3.

BjGmP9HJ3sHm<iN"k4j:`

3. 1

Fig. 2 The stable control systemr{NN`@~VNN_JkMH
7F!0K(9 ℓ ! Pade awrQ$k .
B [s]
J4Ke
≈
A [s]
TD (0 < T − ≤ T ≤ T + )

4)

−TD s

1

A1 [s] =(−1)l sl + (−1)l−1 q1 sl−1 + · · · + ql
B1 [s] =sl + q1 sl−1 + · · · + ql
(l + i)!
qi =(T )−i , l = 1, 2, · · ·(l − i)!i!

11P =10001000−01+0−+12)

1J5K

33G,1−T s1=

B[s] −T se
A[s]B[s]B1 [s] + (A1 [s]B[s]e−Ts − B[s]B1 [s])A[s]A1 [s]J6K

0+iN_JkbGkO!0Nh&K?(ilk.
B[s]B [s]P̂ (s)=, ∆N (s) = P (s) − P̂ (s)J7K
A[s]A [s]sF9k3sHm<iO Fig.2 K(9 Youla-Kucera Qia
Hj<<7gsKhCF?(ilk.
R(s)= (Y (s) − Q(s)N (s)) K(s)J8K
U (s)= X(s) + Q(s)D(s)
J9K
V (s)= Y (s) − Q(s)N (s)J10K
33G, Q(s) ∈ RH O/BjJU#<IPC/3sHm<iN_WQia<? G"j, I>-= K(s) O"3~V GPC! KhCF_W5lk. X $Y $N $D$K, Q N_WO!
(4)

1

1−1∞3)

12)

W,+0)fXq:H@8 h 7 , h 7 f 2008 /50

NH*j?(ilk. GiK, X(s) ∈ RH H Y (s) ∈ RH
r!0N Bezout y0r~-9kh&Kaak.
X(s)N (s) + Y (s)D(s) = 1
J11K
D.K, ~O)sr}?J$Dk<WONBjJC-?`0T [s] rjA7, eRN{s,rrJ<Nh&Kjak ’ .
∞∞3) 10)0gC[s] + F0 [s],C[s]K(s) = g

Y (s) =X(s) =C[s] + G0 (s)
C[s]
12J K

iNPOO/69klgb"j , ~O)sKP~7?"sAo$sI"CW3sHm<i,,WG"k. \@8GO$~O)srTN+Js~AWGKV-9(F$3Ns~AW
GKmP9HJ3sHm<iraak. 9JoA$D.Nk
L ’ rQ$k3HKhCFmP9Hroraak3H,G-k. Fig.2 G(7?s~AOO Fig.3 G(7?OH_J9
3H,G-k. 33G, φ OD.NX8r~-9k ’ ’ .
7)

3) 8)

3) 7) 8)33G, g, F [s], G [s], L [s], !0N E [s], F [s], G [s],H [s], L [s] O8% ’ KhCFjak. ^?,s C[s]F [s]=+ E [s]J13K
A[s]A [s]A[s]A [s]G [s]B[s]B [s]E [s]=+ H [s]
J14K
C[s]
C[s]
0k

00kkk3) 10)

kkkk

1

11kkk

J15K
GeK, Q(s) O Q(s) r^`$(9), (10) 0N U (s), V (s), *hS$S(s) = V (s) U (s) rQ$F$!0NTy0r~-9kh&K*V.
sk B[s]B1 [s]
Lk [s]= Hk [s] +A[s]A1 [s]A[s]A1 [s]Fig. 3 The equivalent diagram of the proposed control system−1k"−S(s)1#(D + N S)

−1k k"∆N0#k<1J16K
3NroOTBjJDk<WO,/BjHJk.2$sj}
H7FNilF$k . 9JoA$3Nror~?9H-B
jJ3sHm<i S(s) KhCF~O)sr}?J$Wm;
9OBjHJk. 33G, ∆N O (7) 0GjA5l$D.N
0r~7 $e0N=jKxQ9k3H,G-k.
∞∞11)

u2 = φ(u1 ) = u1 − 2u
V + U P = Y − QN + (X + QD)(

A1 [jω]B[jω]|{|P̃o (jω)|T2 (ω)
A[jω]+ |P̃e (jω)|T1 (ω)}, ω ∈ [0, ∞)T1 (ω)=

T2 (ω)=(|e−jωT+J21K
Fig.3 G (w, u ) +i (u , y) XN~POVX8 G O,
2"+G12 (s)G21 (s)G22 (s)

#"u2w

#J22K

G22 (s) = (I + V + U P )−1 P V R

NDk<WOO (19) 0Nror~99YFN φ KP
7, D.Nror~?93sHm<i S(s) KhjmP9HB
jHJk.
1. S(s) = V (s) U (s) ,Bj.
Fig.2< ω ≤ T2π++ 1|, if T +2π+T −2, if ω > T2π+−1

2. G11 (s) ∈ RHg,/5B

HJk. q , · · ·, q O (5) 0G?(ilF$k. (16) 0r~9h&K Q r*VH, GeK, u = u H7, (9) 0, (10) 0+i$N`@~VKmP9H/BjJ3sHm<i S(s) = V U
r_W9k3H,G-k.
l

1−1~O)sr<&TN+JOKP9kmP9H_W~O)sN"kWm;9KP7FO$TBjJ3sHm<3. 2=G11 (s)G12 (s) = 2(I + V + U P )−1 V R− 1|, if 0 ≤ ω ≤ Tπ+ql + ql−2 s2 + · · ·Pee (s)= ls + q1 sl−1 + · · · + ql
ql−1 s + ql−3 s3 + · · ·Peo (s)= ls + q1 sl−1 + · · · + ql
1y1"

G21 (s) = −(I + V + U P )−1 P

2, if ω > Tπ+|e−jωTu1#

G11 (s) = −(I + V + U P )−1 (V + U P − I)

J17K

−jωT −
+ 1|, if 0 ≤ ω ≤ T +2π |e+T −

N (s) + ∆N (s))D(s)= D−1 (1 + U ∆N )4)|∆N (jω)| ≤ |

J18K
J19K
J20K1(z − φ(z))
2|φ(z)| ≤ |z|
σ(z) =

3. M D(s)M −1U (jω)∆N (jω)
| < 1, ω ∈ [0, ∞)4. |D(jω)3),8),M.

33G O5jMPQTsG"k hCF$0aNro (16) 0KC(F$eRN 4 DNror~93sHm<iraalPh$. 33G$ro 1 O (16) 0KhCF~?5lF$k.

W,+0)fXq:H@8 h 7 , h 7 f 2008 /4.

~AoN9Y)fn3kL\8M)f!t)fhV:N = 1 )fE_:λ = 0.1
G.=,hV:T = 0
Gg=,hV:T = 29_WQia<?:u = −0.001, u = 1,u12ndg = 0.2296, C[s] = s + 3,
F0 [s] = −17.5157s − 1.263, G0 [s] = −1.589,
L0 [s] = −27.6356s − 2.1979, Q(s) = −0.005
Table.2 Strongly stable controller design parameters:
T0 = 58.83s2 +7.0149s+0.2163

Dk<WONC-}x0/Bj Dk<WONK

3)−88.6885, −11.4746, −0.0358

D.N3AES/E79F`K*1k~AoNbGkrQ$
F, \sFj!NB!,QrTC?.
B[s] = 1

J23K
33GO`@~V,;$?a 1 !Naw (l = 1)$9JoA
T−= 2, T+= 3.7, TD = 3.5A1 [s] = 3.7s + 2,

B1 [s] = −3.7s + 2

rQ$?. n3GO n = 1 H7, \YJ_WQia<?O
Table.1, Table.2 K(9LjG"k.
)f'rsF!KpE-aa$
C + G0 = s + 1.4191,@?. D.K$(12) 0rQ$k.
−17.286s − 0.5742,
s+31113)I

5.

gC + F0 = −17.2831s − 0.5742

X(s) =1

D−0.0988 + 0.1078j, −0.0988 − 0.1078j,

A[s] = 15.9s + 1,

K(7?)fo C K\sFj!r,Q7?H-Nn3kLr Fig.4 *hS Fig.5 K(9. Fig.4 O~Ao9Y
JFig.1 K*1k T KNW,9YQ=r(7F*j, B~OB
SMr$@~O\8Mr(7F$k. Fig.5 O Fig.1 K*1
k C NPOG"k3A!kLr(9. \n3K*$FOB
Q!K*$FHQ5lF$kj!$9JoA$a/\89Y
K~#5;k?a$^:$g-5 8[g/min] N9FCW~Or?($\89Yr6(kb$M 740 nK~#5;$=Ne$90sec eK, \89Yr 710 nH7F0j)fr+O7?.
l}$>h)f!H7F Fig.6 K(92+3Y)f! K
hk18i|roNbHGN)fkLr Fig.7 *hS Fig.8 K(9. )f+O0KO1MN`nrTC?. 3NH-N Fig.6fN)fQia<?O$\8MU#k? α = 0.4, β = 1.33,fc2$s K = 1.0, Q,~V T = 500[sec], y,~V
T = 0[sec], η = 0.8. Fig.7 O Fig.4 H1MK~Ao9Yr(7$Fig.8 O Fig.5 H1MK3A!kLr(9.
>B3K*$F, )f+O0O18roGB3rTC?,,Fig.4 GO 745.7 n+i)f,O^j, Fig.7 GO 743.0 n+iG"j, 2.7 nN9,"k,, 3N9O\n3,B!n3G"j, btN(^NuVNc$KhkbNG"k,, \89710 nKP7F 0.4%G"j, )fkLNfSKOFA7J$.
Fig.4$Fig.7 HbK*<P<7e<Hr7F$k,$=N}
3O, )f+O~@GO~AobN~AI9YO\89Yr
6(F$k,$~A(^O.FLr-7F=,K:95lF$J$?a, ex$~A?~J[.?~K,TolF*i:9
Ye:,3/?aG"k. 7+7, Fig.4 GO Fig.7 KfS7
F*<P<7e<H,.5/$s 700sec G0j7F$kNK
P7F$Fig.7 GOs 1300sec G0j7F$k3H+i$\sFj!KhjI>-,~e7F*j$\sFj!,-zG"k3H,(5lk.
Fig.1

Table.1 CGPC design parameters: w=1: Ny = 6CGPC51

Y (s) =s + 1.4191
s+358.83s2 + 35.5s + 2
58.83s2 + 7.0149s + 0.2613−3.7s + 2N (s) =58.83s2 + 7.0149s + 0.2613
D(s) =

hCF$)f'O (9) 0,(10) 0hj!0HJk.
gC[s] + F0 [s]+ Q(s)D(s)U (s) =X(s) + Q(s)D(s) =C[s]−17.2861s − 0.5742=
s+3
58.83s2 +35.5s+2− 0.005 ∗58.83s2 +7.0149s+0.2613
C[s] + G0 [s]− Q(s)N (s)V (s)=Y (s) − Q(s)N (s) =C[s]s + 1.4191−3.7s + 2=+ 0.005 ∗s+358.83s2 +7.0149s+0.2613

k@`@~VKTN+5*hS~O)BN2DNC-r-9k
3AES/E79F`K*1k~AoN9Y0j)fKBj
CAGPC 79F`N,QrTC?. ~O)s*hSTN+J`@~VNro<G, {s,r==H Youla-Kucera Qia
Hj<<7gsKhCF3sHm<ir_W7?. =NkL,
sFj!N-z-rB!n3rL8FN'G-?.
K

2 M 8 %

1 H. Demircioglu, ”Constrained continuous-time generalized
predictive control,” IEE Proc. Control Theory and Appl.,
146-5, 470/475, 1999.
2 M. Cannon and B. Kouvaritakis, ”Infinite horizon predictive control of constrained continuous-time linear systems,” Automatica, 36-7, 943/955, 2000.
3 M. Deng, A. Inoue, A. Yanou, and Y. Hirashima,”Continuous-time anti-windup generalized predictive control of non-minimum phase processes with input constraints,” Proc. of IEEE CDC, 2003, 4457/4462, 2003.
4 Deng,,, ””,, 35-7, 852/860,

KKK

df e\ TN+5rbDWm;9OKP9k1c,~)fO_W W,+0)fXq@88

52W,+0)fXq:H@8 h 7 , h 7 f 2008 /1999.
5 B. Kouvaritakis, J. A. Rossiter, and A. O. T. Chang, ”Stable generalized predictive control: an algorithm with guaranteed stability,” IEE Proc. Control Theory and Appl.,
139-4, 349/362, 1992.
6 M. V. Kothare, P. J. Campo, M. Morari, and C. N. Nett,”A unified framework for the study of anti-windup designs,” Automatica, 30-12, 1869/1883, 1994.
7 P. J. Campo, M. Morari, and C. N. Nett, ”Multivariableanti-windup and bumpless transfer: A general theory,”
Proc. of ACC, 1989, 1706/1711, 1989.
8,, ”Youla”,, 13-9, 395/402, 2000.
9,”,”,42-4, 330/333, 2003.
10, ,,,, ”YoulaParametrization
GPC”
8SICE, 214/215, 1999.
11 M. Deng, A. Inoue, K. Taketa, and Y. Hirashima,”Continuous-time anti-windup generalized predictive control of uncertain processes with input constraints and timedelays,” Proc. of the 43rd IEEE Conference on Decision
and Control (CDC), 5053/5058, 2004
12 H. Demircioglu and P. J. Gawthrop, ”Continuous-time
generalized predictive control(CGPC),” Automatica, 271, 55/74, 1991.
13,””,, 2000.

KKK

Fig. 4 Temperature change at T1 by CAGPCKYD k)
QiaHj<$7gsKpE/bGkI>?"sAo$sI"CW)fON>\_W 79F`)fpsXq@8o
K3\$gh Wm;9)fN=_$an$$h W,H)f
Kp< Y fe }D ?h {s,r==rQ$?
Khk"3~V NH% h s
fqYtXQViq@88
KK

Fig. 5 Control input by CAGPC

K-f G#8?kWu)f79F`NpCH~Q )H;QR

Nx T R pO

P 6 > 'J5qwK

1989 /e#gX)Xt4H. 1991 /e#gX)X&fJ=X!#)Xl6$N]x$;. 1/0)E)H ;Q\tK~R$=_Kjk. 2005/+i,3gXgX!+3JX&fJnN]x_
Xf. gH7F)f79F`N+/$Wm;9)fN&fK>v.

Fig. 6 Control System of Two-degree freedom Control

H@ oJ5qwK

f e<J5qwK

Fig. 7 Temperature change at T1 by Two degree freedomcontrol

Fig. 8 Control input by Two degree freedom control

1997 /'\gXgX!nN]x$;. 1/'\gXuj$2000 /$.j9 Exeter gX Research
Associate$2001 / 11 n NTT 3_eK1<7gsJXpC&fjj5<A"=7($H$2002 /11 n,3gX)XtK\k. =_1gX!+3J
X&fJZ5x. 79F`)f$NcGGJIN&fK>v. nNJXQK.

1968 /~TgXgX!$N]x$;. ~TgX
uj$+J@ Alberta gX Research Fellow$'\gXu5x$5x$1987 /5xH7F,3gX)XtK\k. =_1gX!+3JX&fJ5x.
=,)f${s,r$NcGGJIN&fK>v.
SICE Fellow. )XnN.