# In-Domain Control of a Heat Equation: An Approach Combining Zero-Dynamics Inverse and Differential Flatness

Authors: Jun Zheng; Guchuan Zhu
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

This paper addresses the set-point control problem of a one-dimensional heat equation with in-domain actuation. The proposed scheme is based on the framework of zero-dynamics inverse combined with flat system control. Moreover, the set-point control is cast into a motion planning problem of a multiple-input, multiple-output system, which is solved by a Green’s function-based reference trajectory decomposition. The validity of the proposed method is assessed through the analysis of the invertibility of the map generated by Green’s function and the convergence of the regulation error. The performance of the developed control scheme and the viability of the proposed approach are confirmed by numerical simulation of a representative system.


---

## Body

## 1. Introduction

Control of parabolic partial differential equations (PDEs) is a long-standing problem in PDE control theory and practice. There exists a very rich literature devoted to this topic, and it is continuing to draw a great attention for both theoretical studies and practical applications. In the existing literature, the majority of work is dedicated to boundary control, which may be represented as a standard Cauchy problem to which functional analytic setting based on semigroup and other related tools can be applied (see, e.g., [1–4]). It is interesting to note that, in recent years, some methods that were originally developed for the control of finite-dimensional nonlinear systems have been successfully extended to the control of parabolic PDEs, such as backstepping (see, e.g., [5–7]), flat systems (see, e.g., [8–14]), and their variations (see, e.g., [15, 16]).This paper deals with the output regulation problem for set-point control of a one-dimensional heat equation via pointwise in-domain (or interior) actuation. Notice that, due to the fact that the regularity of a pointwise controlled inhomogeneous heat equation is qualitatively different from that of boundary controlled heat equations, the techniques developed for boundary control may not be directly applied to the former case. This constitutes a motivation for the present work. The control scheme developed in this paper is based-on the framework of zero-dynamics inverse (ZDI), which was introduced by Byrnes and Gilliam in [17] and has been exploited and developed in a series of works (see, e.g., [18] and the references therein). It is pointed out in [19] that “for certain boundary control systems it is very easy to model the system’s zero dynamics, which, in turn, provides a simple systematic methodology for solving certain problems of output regulation.” Indeed, the construction of zero-dynamics for output regulation of certain interiorly controlled PDEs is also straightforward (see, e.g., [20]) and hence, the control design can be carried out in a systematic manner. Nevertheless, a main issue related to the application of ZDI method is that it leads to, in general, a dynamic control law. Thus, the implementation of such control schemes requires resolving the corresponding zero-dynamics, which may be very difficult for generic regulation problems, such as set-point control considered in the present work. To overcome this difficulty, we resort to the theory of flat systems [13, 21]. We show that, in the context of ZDI design, the control can be derived from the so-called flat output without explicitly solving the original dynamic equation. Moreover, in the framework of flat systems, set-point control can be cast into a problem of motion planning, which can also be carried out in a systematic manner. Note that it can be expected that the ZDI design is applicable to other systems, such as the interior control of beam and plate equations, as an alternative to the methods proposed in, for example, [11, 22, 23].The system model used in this work is taken from [20]. In order to perform control design based on the principle of superposition, we present the original system in a form ofparallel connection. As the control with multiple actuators located in the domain leads to a multiple-input, multiple-output (MIMO) problem, we introduce a Green’s function-based reference trajectory decomposition scheme that enables a simple and computational tractable implementation of the proposed control algorithm.The remainder of the paper is organized as follows. Section2 describes the model of the considered system and its equivalent settings. Section 3 presents the detailed control design. Section 4 deals with motion planning and addresses the convergence and the solvability of the proposed control scheme. A simulation study is carried out in Section 5, and, finally, some concluding remarks are presented in Section 6.


## 2. Problem Setting

In the present work, we consider a scalar parabolic equation describing one-dimensional heat transfer with in-domain control, which is studied in [20]. Denote by u
(
x
,
t
) the temperature distribution over the one-dimensional space, x, and the time, t. The derivatives of u
(
x
,
t
) with respect to its variables are denoted by u
x and u
t, respectively. Consider m points x
j, j
=
1
,
…
,
m, in the interval (
0,1
) and assume, without loss of generality, that 0
=
x
0
<
x
1
<
x
2
<
⋯
<
x
m
<
x
m
+
1
=
1. Let Ω
≐
⋃
j
=
0
m
(
x
j
,
x
j
+
1
). The considered heat equation with in-domain control in a normalized coordinate is of the form (1a)
u
t
x
,
t
-
u
x
x
x
,
t
=
0
,
x
∈
Ω
,
t
>
0
,
(1b)
u
x
,
0
=
ϕ
x
,
(1c)
B
0
u
=
u
x
0
,
t
-
k
0
u
0
,
t
=
0
,
B
1
u
=
u
x
1
,
t
+
k
1
u
1
,
t
=
0
,
(1d)
u
x
j
+
=
u
x
j
-
,
j
=
1
,
…
,
m
,
(1e)
B
x
j
u
=
u
x
x
j
=
v
j
t
,
j
=
1,2
,
…
,
m
,where for a function f
(
·
) and a point x
∈
[
0,1
] we define(2)
f
x
=
f
x
+
-
f
x
-
,with x
- and x
+ denoting, respectively, the usual meaning of left and right hand limits to x. The initial condition is specified in (1b) with ϕ
(
x
)
∈
L
2
(
0,1
). It is assumed that, in system ((1a), (1b), (1c), (1d), and (1e)), we can control the heat flow at the points x
j for j
=
1
,
…
,
m; that is,(3)
v
j
t
=
u
x
x
j
=
u
x
x
j
+
,
t
-
u
x
x
j
-
,
t
.Note that, in ((1a), (1b), (1c), (1d), and (1e)), B
x
j
,
x
j
∈
(
0,1
), represents the pointwise control located in the domain.The space of weak solutions to system ((1a), (1b), (1c), (1d), and (1e)) is chosen to be H
1
(
0,1
). Note that system ((1a), (1b), (1c), (1d), and (1e)) is exponentially stable in H
1
(
0,1
) if B
0 and B
1 are chosen such that k
0
>
0 and k
1
>
0 [19].Denote a set of reference signals corresponding to the control support points of in-domain actuation by{
u
i
D
(
x
i
,
t
)
}
i
=
1
m, where u
i
D
(
x
i
,
t
)
∈
C
∞, i
=
1
,
…
,
m, for all t
∈
(
0
,
T
) and T
<
∞. Let e
i
(
t
)
=
u
(
x
i
,
t
)
-
u
i
D
(
x
i
,
t
) be the regulation errors. Let e
(
t
)
=
{
e
i
(
t
)
}
i
=
1
m and v
(
t
)
=
{
v
i
(
t
)
}
i
=
1
m.Problem 1.
The considered regulation problem for set-point control is to find a dynamic controlv
(
t
) such that the regulation error satisfies e
(
t
)
→
0 as t
→
∞.Note that although the model under the form ((1a), (1b), (1c), (1d), and (1e)) allows deducing easily the zero-dynamics, it is not convenient for motion planning and, in particular, for establishing the input-output map, which is essential for feedforward control design. For this reason, we introduce an equivalent formulation of the in-domain control problem described in ((1a), (1b), (1c), (1d), and (1e)) by replacing the jump conditions in (1e) by pointwise controls as source terms. The resulting system will be of the following form:(4a)
u
t
x
,
t
-
u
x
x
x
,
t
=
∑
j
=
1
m
δ
x
-
x
j
α
j
t
,
0
<
x
<
1
,
t
>
0
,
(4b)
u
x
,
0
=
ϕ
x
,
(4c)
B
0
u
=
u
x
0
,
t
-
k
0
u
0
,
t
=
0
,
B
1
u
=
u
x
1
,
t
+
k
1
u
1
,
t
=
0
,where δ
(
x
-
x
j
) is the Dirac delta function supported at the point x
j, denoting the position of control support, and α
j
:
t
↦
R, j
=
1
,
…
,
m, are the in-domain control signals.Lemma 2.
Considering weak solutions inH
1
(
H
1
0,1
,
[
0
,
T
]
), T
<
∞, system ((1a), (1b), (1c), (1d), and (1e)) and system ((4a), (4b), and (4c)) are equivalent if(5)
α
j
t
=
-
v
j
t
=
-
u
x
x
j
,
j
=
1
,
…
,
m
.Proof.
The proof follows the idea presented in [24]. Indeed, it suffices to prove “system ((1a), (1b), (1c), (1d), and (1e)) ⇒ system ((4a), (4b), and (4c)).” Let X
=
L
2
(
0,1
) be a Hilbert space equipped with the inner product 〈
v
,
w
〉
=
∫
0
1
v
(
x
)
w
(
x
)
d
x, for any v
,
w
∈
X. Let the operator A be defined by A
v
=
v
x
x, with domain D
(
A
)
=
{
v
∈
H
2
(
0,1
)
;
B
0
v
=
B
1
v
=
0
}. It is easy to see that A
∗, the adjoint of A, is equal to A. Let A
~ be an extension of A with the domain D
(
A
~
)
=
{
v
∈
X
;
v
∈
H
2
(
⋃
i
=
0
m
(
x
i
,
x
i
+
1
)
)
,
B
0
v
=
B
1
v
=
0
,
v
(
x
j
+
)
=
v
(
x
j
-
)
,
j
=
1
,
…
,
m
}. Let v
∈
D
(
A
~
), w
∈
D
(
A
∗
)
=
D
(
A
). Using integration by parts we obtain that(6)
A
~
v
,
w
=
v
,
A
v
+
∑
j
=
1
m
v
x
x
j
-
-
v
x
x
j
+
w
x
j
.Let X
-
1
=
D
A
∗
′, the dual space of D
(
A
). We need to define another extension for A. Let A
^
:
H
1
(
0,1
)
→
X
-
1 be defined by(7)
A
^
v
,
w
=
v
,
A
∗
w
∀
w
∈
D
A
∗
,with D
(
A
^
)
=
H
1
(
0,1
). Note that δ
(
·
-
x
j
) is not in X, but in the larger space X
-
1. It follows from (6), (7), and A
=
A
∗ that(8)
A
~
v
=
A
^
v
+
∑
j
=
1
m
v
x
x
j
-
-
v
x
x
j
+
δ
x
-
x
jin X
-
1. If v satisfies system ((1a), (1b), (1c), (1d), and (1e)), then v
˙
(
t
)
=
A
~
v
(
t
), which yields, considering (8), v
˙
(
t
)
=
A
^
v
+
∑
j
=
1
m
(
v
x
(
x
j
-
)
-
v
x
(
x
j
+
)
)
δ
(
x
-
x
j
). Finally, we can see that system ((1a), (1b), (1c), (1d), and (1e)) becomes system ((4a), (4b), and (4c)) with α
j
(
t
)
=
-
v
j
(
t
)
=
-
z
x
x
j, j
=
1
,
…
,
m, where we look for generalized solutions v
(
·
,
t
)
∈
D
(
A
^
)
=
H
1
(
0,1
) such that (8) is true in X
-
1.To establish in-domain control at every actuation point, we will proceed in the way ofparallel connection; that is, for every x
j
∈
(
0,1
), consider the following two systems:(9a)
u
t
x
,
t
-
u
x
x
x
,
t
=
0
,
x
∈
0
,
x
j
∪
x
j
,
1
,
t
>
0
,
(9b)
u
x
,
0
=
ϕ
j
x
,
(9c)
B
0
u
=
u
x
0
,
t
-
k
0
u
0
,
t
=
0
,
B
1
u
=
u
x
1
,
t
+
k
1
u
1
,
t
=
0
,
(9d)
u
x
j
+
=
u
x
j
-
,
(9e)
B
x
j
u
=
u
x
x
j
=
w
j
t
,
(10a)
u
t
x
,
t
-
u
x
x
x
,
t
=
δ
x
-
x
j
β
j
t
,
0
<
x
<
1
,
t
>
0
,
(10b)
u
x
,
0
=
ϕ
j
x
,
(10c)
B
0
u
=
u
x
0
,
t
-
k
0
u
0
,
t
=
0
,
B
1
u
=
u
x
1
,
t
+
k
1
u
1
,
t
=
0
,with ∑
j
=
1
m
ϕ
j
(
x
)
=
ϕ
(
x
). Similarly, systems ((9a), (9b), (9c), (9d), and (9e)) and ((10a), (10b), and (10c)) are equivalent provided u
∈
H
1
(
0,1
) and β
j
=
-
w
j
=
-
u
x
x
j. Let α
j
=
-
v
j
=
β
j
=
-
w
j
=
-
[
u
x
j
]
x
j for all j
=
1,2
,
…
,
m, where u
j denotes the solution to system ((10a), (10b), and (10c)). It follows that u
(
x
,
t
)
=
∑
j
=
1
m
u
j
(
x
,
t
) is a solution to system ((4a), (4b), and (4c)). Moreover,(11)
u
x
x
i
=
∑
j
=
1
m
u
x
j
x
i
=
u
x
i
x
i
=
v
i
,for all i
=
1,2
,
…
,
m. Hence u
(
x
,
t
)
=
∑
j
=
1
m
u
j
(
x
,
t
) is a solution to system ((1a), (1b), (1c), (1d), and (1e)). Therefore, throughout this paper, we assume α
j
=
-
v
j
=
β
j
=
-
w
j
=
-
[
u
x
j
]
x
j for all j
=
1,2
,
…
,
m. Due to the equivalences of systems ((1a), (1b), (1c), (1d), and (1e)) and ((4a), (4b), and (4c)), and systems ((9a), (9b), (9c), (9d), and (9e)) and ((10a), (10b), and (10c)), we may consider ((4a), (4b), and (4c)) and system ((9a), (9b), (9c), (9d), and (9e)) in the following parts.It is worth noting that, in general, the internal pointwise control of the heat equation cannot be allocated at arbitrary points. In particular, the approximate controllability may be lost if the support of the control is located on a nodal set of eigenfunctions (see, e.g., [25, 26]). Therefore, it is important to ensure that the controllability property of the system is insensitive to the location of control support, so that the expected performance can be achieved. Indeed, the loss of controllability will not happen to the considered system if the parameters k
0 and k
1 are chosen appropriately. Precisely, we call a point x
∈
(
0,1
) strategic if it is not located on a nodal set of eigenfunctions of the corresponding PDE [25, 26]. We have then the following.Proposition 3.
For system (4a) with the boundary conditions (4c), all the points x
∈
(
0,1
) are strategic for any k
0
>
0 and k
1
>
0. Furthermore, the approximate controllability of such a system is insensitive to the position of control support for all x
∈
(
k
1
/
(
k
0
+
k
1
)
,
1
).Proof.
The eigenfunctions of system (4a) with the boundary conditions (4c) are given by [27](12)
ψ
n
x
=
cos
⁡
μ
n
x
+
k
0
μ
n
sin
⁡
μ
n
x
,
n
=
1,2
,
…
,where μ
n are positive roots of the transcendental equation(13)
t
a
n
μ
=
μ
k
0
+
k
1
μ
2
-
k
0
k
1
,
k
0
>
0
,
k
1
>
0
.The solution of ψ
n
(
x
)
=
0 is given by(14)
tan
⁡
μ
n
=
-
μ
n
x
k
0
,
∀
x
∈
0,1
,
k
0
>
0
.Obviously, for any fixed x
∈
(
0,1
), the solution of (14) cannot be that of (13) for all n
=
1,2
,
…. Therefore, all the points x
∈
(
0,1
) are strategic for any k
0
>
0 and k
1
>
0. Furthermore, for μ to verify simultaneously (13) and (14), it must hold(15)
μ
2
=
k
0
k
1
-
x
k
0
k
0
+
k
1
.Therefore, μ is real-valued if and only if x
≤
k
1
/
(
k
0
+
k
1
)
<
1.Proposition3 implies that the approximate controllability of the considered system could hold in almost the whole domain by choosing k
0
≫
k
1.


## 3. Control Design and Implementation

In the framework of zero-dynamics inverse, the in-domain control is derived from the so-called forced zero-dynamics, or zero-dynamics for short. To work with theparallel connected system ((9a), (9b), (9c), (9d), and (9e)), we first split the reference signal as(16)
u
D
x
,
t
=
∑
j
=
1
m
γ
j
x
,
x
j
u
j
d
x
j
,
t
,
u
j
d
x
j
,
t
∈
H
1
0
,
T
,
for any
T
<
∞
,where γ
j
(
x
,
x
j
) will be determined in Theorem 7 (see Section 4). Denoting by ε
j
(
t
)
=
u
j
(
x
j
,
t
)
-
u
j
d
(
x
j
,
t
) the regulation error corresponding to system ((9a), (9b), (9c), (9d), and (9e)), the zero-dynamics can be obtained by replacing the input constraints in (9e) by the requirement that the regulation errors vanish identically; that is, ε
j
(
t
)
=
0. Thus, we obtain for a fixed index j
(17a)
ξ
t
x
,
t
=
ξ
x
x
x
,
t
,
x
∈
0
,
x
j
∪
x
j
,
1
,
t
>
0
,
(17b)
ξ
x
,
0
=
0
,
(17c)
ξ
x
0
,
t
-
k
0
ξ
0
,
t
=
0
,
ξ
x
1
,
t
+
k
1
ξ
1
,
t
=
0
,
(17d)
ξ
x
j
,
t
=
u
j
d
x
j
,
t
,where u
j
d
(
x
j
,
t
)
∈
H
1
(
0
,
T
) for any T
<
∞. It should be noticed that, by construction, we can always choose suitable initial data of u
d
(
x
,
0
) so that the zero-dynamics will have a homogenous initial condition. This will significantly facilitate control design. Then, we can get from (9e)(18)
w
j
=
u
x
j
x
j
=
ξ
x
j
x
j
,which shows that the in-domain control for system ((9a), (9b), (9c), (9d), and (9e)) can be derived from the solution of the zero-dynamics. Hence, ((17a), (17b), (17c), and (17d)) and (18) form a dynamic control scheme. This is indeed the basic idea of zero-dynamic inverse design. The convergence of regulation errors with ZDI-based control is given in following theorem.Theorem 4.
The regulation error corresponding to system ((9a), (9b), (9c), (9d), and (9e)), ε
j
(
t
)
=
u
j
(
x
j
,
t
)
-
u
j
d
(
x
j
,
t
), tends to 0 as t tends to ∞ for any u
j
d
(
x
j
,
t
)
∈
H
1
(
0
,
T
) and T
<
∞.The proof of Theorem4 can follow the development presented in Section III of [20] for the case of a heat equation with one in-domain actuator and hence, it is omitted. Note that a key fact used in the proof of this theorem is that the system given in ((9a), (9b), (9c), (9d), and (9e)) without interior control is exponentially stable for any initial data ϕ
j
(
x
)
∈
L
2
(
0,1
) if k
0
>
0 and k
1
>
0.To implement the dynamic control scheme composed of ((17a), (17b), (17c), and (17d)) and (18), we resort to the technique of flat systems [11, 13, 28, 29]. In particular, we apply a standard procedure of Laplace transform-based method to find the solution to ((17a), (17b), (17c), and (17d)). Henceforth, we denote by f
^
(
x
,
s
) the Laplace transform of a function f
(
x
,
t
) with respect to the time variable. Since ((17a), (17b), (17c), and (17d)) has a homogeneous initial condition, then for fixed x
j
∈
(
0,1
), the transformed equations of ((17a), (17b), (17c), and (17d)) in the Laplace domain read as(19a)
s
ξ
^
x
,
s
=
ξ
^
x
x
x
,
s
,
x
∈
0
,
x
j
∪
x
j
,
1
,
s
∈
C
,
(19b)
ξ
^
x
0
,
s
-
k
0
ξ
^
0
,
s
=
0
,
ξ
^
x
1
,
s
+
k
1
ξ
^
1
,
s
=
0
,
(19c)
ξ
^
x
j
,
s
=
u
^
j
d
x
j
,
s
.We divide ((19a), (19b), and (19c)) into two subsystems, that is, for fixed x
j
∈
(
0,1
), considering(20a)
s
ξ
^
x
,
s
=
ξ
^
x
x
x
,
s
,
0
<
x
<
x
j
,
s
∈
C
,
(20b)
ξ
^
x
0
,
s
-
k
0
ξ
^
0
,
s
=
0
,
(20c)
ξ
^
x
j
,
s
=
u
^
j
d
x
j
,
s
,
(21a)
s
ξ
^
x
,
s
=
ξ
^
x
x
x
,
s
,
x
j
<
x
<
1
,
s
∈
C
,
(21b)
ξ
^
x
1
,
s
+
k
1
ξ
^
x
,
s
=
0
,
(21c)
ξ
^
x
j
,
s
=
u
^
j
d
x
j
,
s
.Letξ
^
-
j
(
x
,
s
) and ξ
^
+
j
(
x
,
s
) be the general solutions to ((20a), (20b), and (20c)) and ((21a), (21b), and (21c)), respectively, and denote their inverse Laplace transforms by ξ
-
j
(
x
,
t
) and ξ
+
j
(
x
,
t
). The solution to ((17a), (17b), (17c), and (17d)) can be written as(22)
ξ
j
x
,
t
=
ξ
-
j
x
,
t
χ
0
,
x
j
+
ξ
+
j
x
,
t
χ
x
j
,
1
,where(23)
χ
x
Ω
j
=
1
,
x
∈
Ω
j
⊆
0,1
;
0
,
otherwise
.Then at each point x
i
∈
(
0,1
), by (18) and the argument of “parallel connection” (see Section 2), we have [
u
x
]
x
i
=
∑
j
=
1
m
[
u
x
j
]
x
i
=
[
u
x
i
]
x
i
=
[
ξ
x
i
]
x
i, i
=
1
,
…
,
m. Hence the in-domain control signals of system ((1a), (1b), (1c), (1d), and (1e)) can be computed by(24)
v
i
=
u
x
x
i
=
ξ
x
i
x
i
,
i
=
1
,
…
,
m
.In the following steps, we present the computation of the solution to system ((17a), (17b), (17c), and (17d)), ξ
j. Issues related to the generation reference trajectory u
D
(
x
,
t
) for system ((1a), (1b), (1c), (1d), and (1e)) will be addressed in Section 4.Note thatξ
^
-
j
(
x
,
s
) and ξ
^
+
j
(
x
,
s
), the general solutions to ((20a), (20b), and (20c)) and ((21a), (21b), and (21c)), are given by(25)
ξ
^
-
j
x
,
s
=
C
1
ϕ
1
x
,
s
+
C
2
ϕ
2
x
,
s
,
ξ
^
+
j
x
,
s
=
C
3
ϕ
1
x
,
s
+
C
4
ϕ
2
x
,
s
,with(26)
ϕ
1
x
,
s
=
sinh
⁡
s
x
s
,
ϕ
2
x
,
s
=
cosh
⁡
s
x
.We obtain by applying (20b) and (20c)(27)
C
1
ϕ
1
x
j
,
s
+
C
2
ϕ
2
x
j
,
s
=
u
^
j
d
x
j
,
s
,
C
1
-
k
0
C
2
=
0
,which can be written as(28)
ϕ
1
x
j
,
s
ϕ
2
x
j
,
s
1
-
k
0
C
1
C
2
=
u
^
j
d
x
j
,
s
0
.Let(29)
R
-
j
=
ϕ
1
x
j
,
s
ϕ
2
x
j
,
s
1
-
k
0
,
(30)
u
^
j
d
x
j
,
s
=
-
det
⁡
R
-
j
y
^
-
j
x
j
,
s
.We obtain(31)
C
1
C
2
=
adj
R
-
j
det
⁡
R
-
j
u
^
j
d
x
j
,
s
0
=
k
0
y
^
-
j
x
j
,
s
y
^
-
j
x
j
,
s
.Therefore, the solution to ((20a), (20b), and (20c)) can be expressed as(32)
ξ
^
-
j
x
,
s
=
k
0
ϕ
1
x
+
ϕ
2
x
y
^
-
j
x
j
,
s
.We may proceed in the same way to deal with ((21a), (21b), and (21c)). Indeed, letting(33)
R
+
j
=
ϕ
1
x
j
,
s
ϕ
2
x
j
,
s
ϕ
2
1
,
s
+
k
1
ϕ
1
1
,
s
s
ϕ
1
1
,
s
+
k
1
ϕ
2
1
,
s
,
(34)
u
^
j
d
x
j
,
s
=
det
⁡
R
+
j
y
^
+
j
x
j
,
s
,we get from ((21a), (21b), and (21c))(35)
C
3
C
4
=
s
ϕ
1
1
,
s
+
k
1
ϕ
2
1
,
s
y
^
+
j
x
j
,
s
-
ϕ
2
1
,
s
+
k
1
ϕ
1
1
,
s
y
^
+
j
x
j
,
s
,
(36)
ξ
^
+
j
x
,
s
=
s
ϕ
1
1
,
s
+
k
1
ϕ
2
1
,
s
ϕ
1
x
+
ϕ
2
1
,
s
+
k
1
ϕ
1
1
,
s
ϕ
2
x
y
^
+
j
x
j
,
s
.Applying the results from [30, 31], which are based on module theory, to (30) and (34), we may choose y
^
j
(
x
j
,
s
) as a basic output such that(37)
y
^
+
j
x
j
,
s
=
-
d
e
t
R
-
j
y
^
j
x
j
,
s
,
y
^
-
j
x
j
,
s
=
d
e
t
R
+
j
y
^
j
x
j
,
s
.It should be noticed that the concept of basic outputs (or flat outputs) plays a central role in flat system theory, because the system trajectory and the input can be directly computed from a basic output and its time derivatives [13, 21].Using the property of hyperbolic functions, we obtain from (32) and (36) that(38)
ξ
^
-
j
x
,
s
=
k
1
sinh
⁡
s
x
j
-
s
s
-
cosh
⁡
s
x
j
-
s
k
0
sinh
⁡
s
x
s
+
cosh
⁡
s
x
y
^
j
x
j
,
s
,
ξ
^
+
j
x
,
s
=
k
1
sinh
⁡
s
x
-
s
s
-
cosh
⁡
s
x
-
s
k
0
sinh
⁡
s
x
j
s
+
cosh
⁡
s
x
j
y
^
j
x
j
,
s
.Note that(39)
ξ
^
j
x
,
s
=
ξ
^
-
j
x
,
s
χ
0
,
x
j
+
ξ
^
+
j
x
,
s
χ
x
j
,
1is a solution to ((19a), (19b), and (19c)). Using the fact (40)
sinh
⁡
x
=
∑
n
=
0
∞
x
2
n
+
1
2
n
+
1
!
,
cosh
⁡
x
=
∑
n
=
0
∞
x
2
n
2
n
!
,we obtain the time-domain solution to ((17a), (17b), (17c), and (17d)), which is given by(41)
ξ
j
x
,
t
=
k
0
k
1
∑
n
=
0
∞
∑
k
=
0
n
x
2
k
+
1
x
j
-
1
2
n
-
k
+
1
2
k
+
1
!
2
n
-
k
+
1
!
y
j
n
-
k
0
∑
n
=
0
∞
∑
k
=
0
n
x
2
k
+
1
x
j
-
1
2
n
-
k
2
k
+
1
!
2
n
-
k
!
y
j
n
+
k
1
∑
n
=
0
∞
∑
k
=
0
n
x
2
k
x
j
-
1
2
n
-
k
+
1
2
k
!
2
n
-
k
+
1
!
y
j
n
-
∑
n
=
0
∞
∑
k
=
0
n
x
2
k
x
j
-
1
2
n
-
k
2
k
!
2
n
-
k
!
y
j
n
χ
0
,
x
j
+
k
0
k
1
∑
n
=
0
∞
∑
k
=
0
n
x
j
2
k
+
1
x
-
1
2
n
-
k
+
1
2
k
+
1
!
2
n
-
k
+
1
!
y
j
n
-
k
0
∑
n
=
0
∞
∑
k
=
0
n
x
j
2
k
+
1
x
-
1
2
n
-
k
2
k
+
1
!
2
n
-
k
!
y
j
n
+
k
1
∑
n
=
0
∞
∑
k
=
0
n
x
j
2
k
x
-
1
2
n
-
k
+
1
2
k
!
2
n
-
k
+
1
!
y
j
n
-
∑
n
=
0
∞
∑
k
=
0
n
x
j
2
k
x
-
1
2
n
-
k
2
k
!
2
n
-
k
!
y
j
n
χ
x
j
,
1
.Furthermore, by a direct computation we get(42)
ξ
^
x
j
x
j
=
k
0
k
1
s
+
s
sinh
⁡
s
+
k
0
+
k
1
cosh
⁡
s
y
^
j
x
j
,
s
=
k
0
k
1
∑
n
=
0
∞
s
n
2
n
+
1
!
+
k
0
+
k
1
∑
n
=
0
∞
s
n
2
n
!
+
∑
n
=
0
∞
s
n
+
1
2
n
+
1
!
y
^
j
x
j
,
s
.It follows from (24) that, in time domain, the control is given by(43)
v
j
t
=
ξ
x
j
x
j
=
k
0
k
1
∑
n
=
0
∞
y
j
n
x
j
,
t
2
n
+
1
!
+
k
0
+
k
1
∑
n
=
0
∞
y
j
n
x
j
,
t
2
n
!
+
∑
n
=
0
∞
y
j
n
+
1
x
j
,
t
2
n
+
1
!
.Finally, provided u
j
d
(
x
j
,
t
)
=
ξ
j
(
x
j
,
t
), for j
=
1
,
…
,
m, the reference trajectory u
D
(
x
,
t
) can be determined from (16) and (41).


## 4. Motion Planning

For control purpose, we have to choose appropriate reference trajectories, or equivalently the basic outputs. Denote now byu
¯
D
(
x
) the desired steady-state profile. Without loss of generality, we consider a set of basic outputs of the form(44)
y
j
t
=
y
¯
x
j
φ
j
t
,
j
=
1
,
…
,
m
,where φ
j
(
t
) is a smooth function evolving from 0 to 1. Motion planning amounts then to deriving y
¯
(
x
j
) from u
¯
D
(
x
) and to determining appropriate functions φ
j
(
t
), for j
=
1
,
…
,
m.To this aim and due to the equivalence of systems ((1a), (1b), (1c), (1d), and (1e)) and ((4a), (4b), and (4c)), we consider the steady-state heat equation corresponding to system ((4a), (4b), and (4c)):(45a)
u
¯
x
x
x
=
∑
j
=
1
m
δ
x
-
x
j
α
¯
j
,
0
<
x
<
1
,
t
>
0
,
(45b)
u
¯
x
0
-
k
0
u
¯
0
=
0
,
u
¯
x
1
+
k
1
u
¯
1
=
0
.Based on the principle of superposition for linear systems, the solution to the steady-state heat equation ((45a) and (45b)) can be expressed as(46)
u
¯
x
=
∫
0
1
∑
j
=
1
m
G
x
,
ζ
δ
ζ
-
x
j
α
¯
j
d
ζ
=
∑
j
=
1
m
G
x
,
x
j
α
¯
j
,where G
(
x
,
ζ
) is the Green’s function corresponding to ((45a) and (45b)), which is of the form(47)
G
x
,
ζ
=
k
1
ζ
-
k
1
-
1
k
0
x
+
1
k
0
+
k
1
+
k
0
k
1
,
0
≤
x
<
ζ
;
k
1
x
-
k
1
-
1
k
0
ζ
+
1
k
0
+
k
1
+
k
0
k
1
,
ζ
≤
x
≤
1
.Indeed, it is easy to check that G
x
x
(
x
,
ζ
)
=
δ
(
x
-
ζ
) and G
(
x
,
ζ
) satisfies the boundary conditions, G
x
(
0
,
ζ
)
-
k
0
G
(
0
,
ζ
)
=
0 and G
x
(
1
,
ζ
)
+
k
1
G
(
1
,
ζ
)
=
0, the joint condition, G
(
ζ
+
,
ζ
)
=
G
(
ζ
-
,
ζ
), and the jump condition, [
G
x
(
x
,
ζ
)
]
ζ
=
1.Takingm distinguished points along the solution to ((45a) and (45b)), u
¯
(
x
1
)
,
…
,
u
¯
(
x
m
), we get(48)
u
¯
x
1
⋮
u
¯
x
m
=
G
x
1
,
x
1
⋯
G
x
1
,
x
m
⋮
⋱
⋮
G
x
m
,
x
1
⋯
G
x
m
,
x
m
α
¯
1
⋮
α
¯
m
.Note that, in (48), the matrix formed by the Green’s function defined an input-output map in steady-state, which is also called theinfluence matrix.Lemma 5.
The influence matrix chosen as in (48) is invertible. Thus,(49)
α
¯
1
⋮
α
¯
m
=
G
x
1
,
x
1
⋯
G
x
1
,
x
m
⋮
⋱
⋮
G
x
m
,
x
1
⋯
G
x
m
,
x
m
-
1
u
¯
x
1
⋮
u
¯
x
m
.Proof.
Form
=
1, since k
0
>
0,  k
1
>
0, and x
1
∈
(
0,1
), it follows that G
(
x
1
,
x
1
)
=
(
k
1
x
1
-
k
1
-
1
)
(
k
0
x
1
+
1
)
/
k
0
+
k
1
+
k
0
k
1
<
0. Hence it is invertible. We prove the claim for m
>
1 by contradiction. Suppose that the influence matrix is not invertible; then it is of rank m
-
1 or less. Without loss of generality, we may assume that, for some x
n
>
x
i with i
=
1
,
…
,
n
-
1, there exist n
-
1 constants l
1
,
l
2
,
…
,
l
n
-
1 such that(50)
G
x
j
,
x
n
=
∑
i
=
1
n
-
1
l
i
G
x
1
,
x
i
,
j
=
1
,
…
,
n
,where 1
<
n
≤
m and ∑
i
=
1
n
-
1
l
i
2
>
0. Let(51)
G
x
=
G
x
,
x
n
,
F
x
=
∑
i
=
1
n
-
1
l
i
G
x
,
x
i
.Equation (50) shows that F
(
x
)
=
G
(
x
) at every boundary point of x
1
,
x
2
,
x
2
,
x
3
,
…
,
[
x
n
-
1
,
x
n
]. Note that F
(
x
) is a linear function in x
1
,
x
2
,
x
2
,
x
3
,
…
,
[
x
n
-
1
,
x
n
] and that G
(
x
)
=
(
k
1
x
n
-
k
1
-
1
)
(
k
0
x
+
1
)
/
(
k
0
+
k
1
+
k
0
k
1
) in [
x
1
,
x
n
]; that is, G
(
x
) is a linear function in [
x
1
,
x
n
]. Hence F
(
x
)
≡
G
(
x
) in [
x
1
,
x
n
].
ByF
(
x
1
)
=
G
(
x
1
), we get(52)
k
1
x
n
-
k
1
-
1
=
∑
i
=
1
n
-
1
l
i
k
1
x
i
-
k
1
-
1
.By F
(
x
n
)
=
G
(
x
n
), we get(53)
k
0
x
n
+
1
=
∑
i
=
1
n
-
1
l
i
k
0
x
i
+
1
.Therefore(54)
∑
i
=
1
n
-
1
l
i
=
1
.By F
x
(
x
1
+
)
=
G
x
(
x
1
+
) and F
x
(
x
n
-
)
=
G
x
(
x
n
-
), we get(55)
k
0
k
1
x
n
-
k
1
-
1
=
k
0
k
1
∑
i
=
1
n
-
1
l
i
x
i
-
k
0
k
1
+
1
∑
i
=
2
n
-
1
l
i
+
l
1
k
1
=
k
0
k
1
∑
i
=
1
n
-
1
l
i
x
i
+
k
1
∑
i
=
1
n
-
1
l
i
.It follows that ∑
i
=
2
n
-
1
l
i
=
0, which yields, considering (54), l
1
=
1. By F
x
(
x
2
+
)
=
G
x
(
x
2
)
=
F
x
(
x
2
-
), we deduce(56)
l
1
k
1
k
0
x
1
+
1
+
l
2
k
1
k
0
x
2
+
1
+
k
0
∑
i
=
3
n
-
1
l
i
k
1
x
i
-
k
1
-
1
=
l
1
k
1
k
0
x
1
+
1
+
k
0
∑
i
=
2
n
-
1
l
i
k
1
x
i
-
k
1
-
1
,which gives l
2
=
0. Similarly, by F
x
(
x
j
+
)
=
G
x
(
x
j
)
=
F
x
(
x
j
-
)
(
j
=
3,4
,
…
,
n
-
2
), we obtain l
3
=
l
4
=
⋯
=
l
n
-
2
=
0
. Hence l
n
-
1
=
0. Then we deduce from (55) that(57)
k
0
k
1
x
n
-
k
1
-
1
=
k
1
k
0
x
1
+
1
.It follows from (53) that k
0
x
1
+
1
=
k
0
x
n
+
1. We conclude then by (57) that k
0
+
k
1
+
k
0
k
1
=
0, which is a contradiction to k
0
>
0 and k
1
>
0.In steady-state, we can obtain from (43) that(58)
v
¯
j
=
k
0
k
1
+
k
0
+
k
1
y
¯
x
j
=
-
α
¯
j
.Finally, y
¯
(
x
j
) can be computed by (49) and (58) for a given u
¯
D
(
x
).It is worth noting that (49) provides a simple and straightforward way to compute the static control from the prescribed steady-state profile. Indeed, a direct computation can show that applying (49) will result in the same static control obtained in [20] where aserially connected model is used.To ensure the convergence of (41) and (43), we choose the following smooth function as φ
j
(
t
):(59)
φ
j
t
=
0
,
if
t
≤
0
∫
0
t
exp
⁡
-
1
/
τ
1
-
τ
ε
d
τ
∫
0
T
exp
⁡
-
1
/
τ
1
-
τ
ε
d
τ
,
if
t
∈
0
,
T
1
,
if
t
≥
Twhich is known as Gevrey function of order σ
=
1
+
1
/
ε,  ε
>
0 (see, e.g., [13]).Lemma 6.
If the basic outputsφ
j
(
t
),  j
=
1
,
…
,
m, are chosen as Gevrey functions of order 1
<
σ
<
2, then the infinite series (41) and (43) are convergent.The claim of Lemma6 can be proved by following a standard procedure using the bounds of Gevrey functions (see, e.g., [13]):(60)
φ
k
+
1
t
≤
M
k
!
σ
K
k
,
∃
K
,
M
>
0
,
∀
k
∈
Z
≥
0
,
∀
t
∈
t
0
,
T
.Therefore, the details of proof are omitted.Theorem 7.
Assumek
0
>
0 and k
1
>
0. Let the basic outputs φ
j
(
t
),  j
=
1
,
…
,
m, be chosen as (59) with an order 1
<
σ
<
2. Let the reference trajectory of system ((1a), (1b), (1c), (1d), and (1e)) be given by (16) with(61)
γ
j
x
,
x
j
=
-
k
0
k
1
+
k
0
+
k
1
G
x
,
x
j
k
0
x
j
+
1
k
0
x
j
-
1
-
1
,
j
=
1
,
…
,
m
,where G
(
x
,
ζ
) is the Green’s function defined in (47). Then the regulation error of system ((1a), (1b), (1c), (1d), and (1e)) with the control given in (43) tends to zero; that is, e
i
(
t
)
=
u
(
x
i
,
t
)
-
u
i
D
(
x
i
,
t
)
→
0 as t
→
∞, for i
=
1,2
,
…
,
m.Proof.
By a direct computation we have(62)
e
i
t
=
x
i
,
t
-
u
i
D
x
i
,
t
=
u
x
i
,
t
+
∑
j
=
1
m
k
0
k
1
+
k
0
+
k
1
G
x
i
,
x
j
u
j
d
x
j
,
t
k
0
x
j
+
1
k
0
x
j
-
1
-
1
=
u
x
i
,
t
+
∑
j
=
1
m
k
0
k
1
+
k
0
+
k
1
G
x
i
,
x
j
ξ
j
x
j
,
t
k
0
x
j
+
1
k
0
x
j
-
1
-
1
≤
u
x
i
,
t
-
u
¯
x
i
+
u
¯
x
i
+
k
0
k
1
+
k
0
+
k
1
∑
j
=
1
m
G
x
i
,
x
j
y
¯
x
j
+
k
0
k
1
+
k
0
+
k
1
∑
j
=
1
m
G
x
i
,
x
j
ξ
j
x
j
,
t
k
0
x
j
+
1
k
0
x
j
-
1
-
1
-
k
0
k
1
+
k
0
+
k
1
∑
j
=
1
m
G
x
i
,
x
j
y
¯
x
j
.By (58) and (46), it follows u
¯
(
x
i
)
=
-
(
k
0
k
1
+
k
0
+
k
1
)
∑
j
=
1
m
G
(
x
i
,
x
j
)
y
¯
(
x
j
)
. Based on (41), (44), and the property of φ
j
(
t
) we have(63)
ξ
j
x
j
,
t
k
0
x
j
+
1
k
0
x
j
-
1
-
1
⟶
y
¯
x
j
as
t
→
∞
.Note that u
(
x
i
,
t
)
→
u
¯
(
x
i
)
a
s
t
→
∞
. Therefore e
i
t
→
0 as t
→
∞.Remark 8.
For anyx
∈
(
0,1
), replacing x
i by x in the proof of Theorem 7, we can get u
x
,
t
-
u
D
x
,
t
→
0 as t
→
∞, which shows that the solution u
(
x
,
t
) of system ((1a), (1b), (1c), (1d), and (1e)) converges to the reference trajectory u
D
(
x
,
t
) at every point x
∈
(
0,1
).


## 5. Simulation Study

In the simulation, we implement system ((4a), (4b), and (4c)) with 12 actuators evenly distributed in the domain at the spot points {
1
/
13,2
/
13
,
…
,
12
/
13
}. The numerical implementation is based on a PDE solver, pdepe, in Matlab PDE Toolbox. In numerical simulation, 200 points in space and 100 points in time are used for the region [
0,1
]
×
[
0,0.5
]. The basic outputs φ
j
(
t
) used in the simulation are Gevrey functions of the same order. In order to meet the convergence condition given in Lemma 6, the parameter of the Gevrey function is set to ε
=
1.1. The feedback boundary control gains are chosen as k
0
=
10 and k
1
=
1. A perturbation of the form u
(
x
,
0
)
=
cos
⁡
π
x is applied at t
=
0 in the simulation.The desired steady-state temperature distribution is a piecewise linear curve, depicted in Figure1(a), which is a solution to ((45a) and (45b)). The corresponding static controls, α
¯
1
,
…
,
α
¯
12, are shown in Figure 1(b). Note that the dynamic control signals, α
i
(
t
), are smooth functions connecting 0 to α
¯
i for i
=
1
,
…
,
12. The evolution of temperature distribution with static and dynamic control, as well as the corresponding regulation errors with respect to the static profile defined as e
(
x
,
t
)
=
u
(
x
,
t
)
-
u
¯
D
(
x
), is depicted in Figure 2. The simulation results show that the system performs well with the developed control scheme. It can also be seen that the dynamic control provides a faster response time compared to the static one.Figure 1
Controller characteristics: (a) desired profile; (b) static control signals.
(a)
(b)Figure 2
Evolution of temperature distribution: (a) response with static control; (b) regulation error with static control; (c) response with dynamic control; (d) regulation error with dynamic control.
(a)
(b)
(c)
(d)


## 6. Conclusion

This paper presented a solution to the problem of set-point control of temperature distribution with in-domain actuation described by an inhomogeneous parabolic PDE. To apply the principle of superposition, the system is presented in aparallel connection form. The dynamic control problem introduced by the ZDI design is solved by using the technique of flat systems motion planning. As the control with multiple in-domain actuators results in a MIMO problem, a Green’s function-based reference trajectory decomposition is introduced, which considerably simplifies the control design and implementation. Convergence and solvability analysis confirms the validity of the control algorithm and the simulation results demonstrate the viability of the proposed approach. Finally, as both ZDI design and flatness-based control can be carried out in a systematic manner, we can expect that the approach developed in this work may be applicable to a broader class of distributed parameter systems.



---

*Source: 187284-2015-12-31.xml*