# Output Feedback Stabilization of Unstable Heat Equations with Time Delay in Boundary Observation

Authors: Kun-Yi Yang; Lan Chen
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

In this paper, an unstable heat equation where the observation signal suffers from a given time delay has been considered. In order to remove the unstable boundary condition, the original system has been converted into the stabilizing system by backstepping transformation under which the form of the stabilizing controller can be obtained. For the sake of overcoming instability produced by the time delay, firstly, we construct the Luenberger observer and design the predictor. Then, the estimated output feedback controllers have been obtained based on the observer and predictor systems. It is shown that the closed-loop system is stable exponentially. At last, numerical simulations have been given to illustrate the effectiveness of the stabilizing output feedback controllers.


---

## Body

## 1. Introduction

In the past decades, much effort has been made for the problems of unstable or antistable infinite-dimensional systems which can be applied to practical engineering (see [1–13]). For example, reference [14] relates the loss of the heat to the surrounding medium and the destabilizing heat generation inside the rod. The backstepping method introduced into the PDE systems in a few years has systematically been applied to stabilize some unstable or antistable hyperbolic and parabolic equations (see [6–8, 11, 13, 15, 16]). For example, Kang and Guo [6] consider boundary stabilization for a cascade of unstable heat PDE systems where an output feedback is designed by the backstepping transformation. Moreover, reference [13] has applied the backstepping transformation to overcome the destabilizing boundary condition for a wave equation. The predictor-based feedback control may be an effective method to make the unstable reaction-diffusion equation compensate the delay (see [7, 11]). There are some of other references related to stability of parabolic equations, for example [17]. In this paper, we consider a heat equation which is controlled from one end and contains instability at the other end. The output feedback controllers designed by the backstepping method use the Volterra transformation to map an unstable PDE into a stable target PDE. It is shown that the controller based on the observer and predictor is also available for the unstable parabolic equation. In this sense, it is meaningful for this paper to consider the unstable heat equation with time delay in boundary observation.It is a common phenomenon of the practical systems with time delay which means that the observation signals need time to achieve the controller ([18, 19]). Unfortunately, even a small amount of time delay may make the originally stable system unstable [20]. Actually, for distributed parameter control systems, stabilization of the system with observation or control suffering from time delay represents difficult mathematical challenges, as was mentioned in [21]. In recent years, references [22–24] have introduced the separation principle to make the wave, beam, and Schrödinger equations where observation signals suffer from the given time delay stabilized. This paper has been devoted to the parabolic heat equation with the delayed observation signal.In this paper, we consider stabilization of the heat equation with the observation subject to a given time delay where an unstable boundary condition is coupled in the boundaryx=0, under the Neumann boundary control:(1)wtx,t=wxxx,t,0<x<1,t>0,wx0,t=−qw0,t,t≥0,wx1,t=ut,t≥0,yt=w0,t−τ,t>τ,wx,0=w0x,0≤x≤1where wx,t is the state, u is the input or control, y is the observation or output, τ is the given time delay, and q>0 is in the boundary condition which may be large enough to emphasize the instability effect. The objective of this paper is to design output feedback controls based on the estimated state in order to stabilize systems (1) for any given time delay.This paper is organized as follows: In Section2, we present the controller design, observer and predictor construction, and stability analysis for unstable Neumann actuation problems. Then, in the next Appendix section, by the main techniques involved in the previous section, we briefly analyze the stability of the unstable Dirichlet actuation problem. Section 4 gives simulation results.


## 2. Stabilization with Neumman Boundary Control

In this section, we apply the backstepping transformation to obtain the stabilizing feedback for system (1) under the Neumman boundary control. And then, the observer system has been constructed while the predictor system has been designed. At last, the output feedback control based on the estimated state has been shown to stabilize system (1) for any given time delay.


### 2.1. Feedback Control via Backstepping Transform

System (1) is considered in the energy space ℋ=L20,1 and the input(or output) space U=Y=ℂ. The energy of system (1) is(2)Et=12∫01w2x,tdx.For the original system (1), make the invertible change of variable(3)vx,t=wx,t+c0+q∫0xeqx−ηwη,tdη,0≤x≤1,t≥0,c0>0,q>0,whose equivalent transformation is(4)wx,t=vx,t−c0+q∫0xe−c0x−ηvη,tdη,0≤x≤1,t≥0,c0>0,q>0.In order to convert (1) into the exponentially stable system [25](5)vtx,t=vxxx,t,0<x<1,t>0,vx0,t=c0v0,t,t≥0,c0>0,vx1,t=0,t≥0,the feedback control in (1) may be proposed as that(6)ut=−c0+qw1,t−qc0+q∫01eq1−ηwη,tdη,t≥0,c0>0,q>0,under which system (1) becomes that [8](7)wtx,t=wxxx,t,0<x<1,t>0,wx0,t=−qw0,t,t≥0,q>0,wx1,t=−c0+qw1,t−qc0+q∫01eq1−ηwη,tdη,t≥0,c0>0,q>0.Based on the exponential stability of system (5) and the equivalence between systems (5) and (7), we have the following theorem [8].Theorem 1.
For any initial valuew0∈ℋ, under feedback (6), the closed-loop system (7) is exponentially stable.


### 2.2. Design of Observer and Predictor

This section has been devoted to constructing the observer system and designing the predictor system. In the first step, the observer has been constructed to estimate the statewx,s,0≤x≤1,0≤s≤t−τ,t>τ from the known observation ys+τ,0≤s≤t−τ,t>τ stated as the following system:(8)w^sx,s=w^xxx,s,0<x<1,0<s<t−τ,t>τ,w^x0,s=−qys+τ+k1w^0,s−ys+τ,0≤s≤t−τ,t>τ,k1>0,w^x1,s=us,0≤s≤t−τ,t>τ,w^x,0=w^0x,0≤x≤1,where w^0 is the arbitrarily assigned initial value of the observer.Set the error for the observer system(9)εx,s=w^x,s−wx,s,0≤x≤1,0≤s≤t−τ,t>τ.Then, by (1) and (8), εx,s satisfies the following:(10)εsx,s=εxxx,s,0<x<1,0<s<t−τ,t>τ,εx0,s=k1ε0,s,0≤s≤t−τ,t>τ,k1>0,εx1,s=0,0≤s≤t−τ,t>τ,εx,0=w^0x−w0x,0≤x≤1.System (10) can be written as follows:(11)ddsε·,s=Aε·,s,where the operator A is defined as follows:(12)Aφ=φ″,∀φ∈DA,DA=φ∈H20,1|φ′0=k1φ0k1>0,φ′1=0.It is known thatA generates an exponentially stable C0-semigroup on ℋ such that(13)eAsℋ≤Me−ωs,∀s>0,for some positive constants M,ω ([25]). Thus, for any w0∈ℋ and w^0∈ℋ, there exists a unique solution of (12) such that(14)ε·,sℋ≤Me−ωsε·,0ℋ,∀0≤s≤t−τ,t>τ.In the second step, the predictor system has been designed for the statewx,s,0≤x≤1,t−τ≤s≤t,t>τ as follows:(15)w^sx,s,t=w^xxx,s,t,0<x<1,t−τ<s<t,t>τ,w^x0,s,t=−qw^0,s,t,t−τ≤s≤t,t>τ,q>0,w^x1,s,t=us,t−τ≤s≤t,t>τ,w^x,t−τ,t=w^x,t−τ,0≤x≤1,t>τ,Let the error for the predictor system(16)εx,s,t=w^x,s,t−wx,s,0≤x≤1,t−τ≤s≤t,t>τ.Then, by (1), (15), and (16), εx,s,t satisfies the following system:(17)εsx,s,t=εxxx,s,t,0<x<1,t−τ<s<t,t>τ,εx0,s,t=−qε0,s,t,t−τ≤s≤t,t>τ,q>0,εx1,s,t=0,t−τ≤s≤t,t>τ,εx,t−τ,t=εx,t−τ,0≤x≤1,t>τ.It can be expressed as that(18)ddsε·,s,t=ℬε·,s,t,where(19)ℬφ=φ″,∀φ∈Dℬ,Dℬ=φ∈H20,1||φ′0=−qφ0q>0,φ′1=0.Next, we define a new inner product equivalent as the energy in (2)(20)φ1,φ2ℋ1=∫01coshα1−xφ1xφ2xdx,for some given α>0 satisfying that 2qcoshα−αsinhα<0.Thus, for anyφ∈Dℬ, we have that(21)φ,φℋ1=∫01coshα1−xφx2dx.By simple computations of integration by parts, it results that(22)ℬφ,φℋ1=φ″,φℋ1=∫01coshα1−xφ″xφxdx,=qcoshα−αsinhα2φ02+∫01coshα1−xα2φx2/2−φ′x2dx,≤α2/2·φ,φℋ1.Thus, system (17) generates a strongly continuous semigroup eℬs which indicates that(23)eℬsℋ≤M1eω1s,∀s≥0,for some positive constants M1 and ω1 [26]. Therefore, the solution of system (17) can be represented as follows:(24)ε·,s,t=eℬs−t+τε·,t−τ,t−τ≤s≤t,t>τ,which together with (14) and (23) result that(25)ε·,t,tℋ≤MM1eω+ω1τε·,0ℋe−ωt,∀t>τ.The estimated state has naturally been chosen as follows:(26)w˜x,t=w^x,t,t,∀t>τ.assured by the theorem below.Theorem 2.
For∀t>τ, we have that(27)w·,t−w˜·,tℋ≤MM1eω+ω1τw^0−w0ℋe−ωt,where w^0 is the initial value of the observer (8), w0 is the initial value of the original system (1), and M,ω and M1,ω1 are constants in (13) and (23), respectively.


### 2.3. Stabilization of the Closed-Loop System

For the estimated state (24) and (7), we naturally obtain the estimated output feedback control law as follows:(28)u∗t=0,0≤t≤τ,−c0+qw˜1,t−qc0+q∫01eq1−ηw˜η,tdη,t>τ,c0>0,q>0.Then, the closed-loop system becomes a system of partial differential equations as given below:(29)wtx,t=wxxx,t,0<x<1,t>0,wx0,t=−qw0,t,t≥0,q>0,wx1,t=u∗t,t≥0,wx,0=w0x,0≤x≤1,(30)w^sx,s=w^xxx,s,0<x<1,0<s<t−τ,t>τ,w^x0,s=−qys+τ+k1w^0,s−ys+τ,0≤s≤t−τ,t>τ,q>0,k1>0,w^x1,s=u∗s,0≤s≤t−τ,t>τ,w^x,0=w^0x,0≤x≤1,(31)w^sx,s,t=w^xxx,s,t,0<x<1,t−τ<s<t,t>τ,w^x0,s,t=−qw^0,s,t,t−τ≤s≤t,t>τ,q>0,w^x1,s,t=u∗s,t−τ≤s≤t,t>τ,w^x,t−τ,t=w^x,t−τ,0≤x≤1,t>τ.which is equivalent to the system below:(32)wtx,t=wxxx,t,0<x<1,t>τ,wx0,t=−qw0,t,t>τ,q>0,wx1,t=−c0+qw1,t+ε1,t,t−qc0+q∫01eq1−ηwη,t+εη,t,tdη,t>τ,q>0,wx,0=w0x,0≤x≤1,(33)εsx,s=εxxx,s,0<x<1,0<s<t−τ,t>τ,εx0,s=k1ε0,s,0≤s≤t−τ,t>τ,k1>0,εx1,s=0,0≤s≤t−τ,t>τ,εx,0=ε0x,0≤x≤1,(34)εsx,s,t=εxxx,s,t,0<x<1,t−τ<s<t,t>τεx0,s,t=−qε0,s,t,t−τ≤s≤t,t>τ,q>0,εx1,s,t=0,t−τ≤s≤t,t>τ,εx,t−τ,t=εx,t−τ,0≤x≤1,t>τ.Theorem 3.
For∀t>τ, anyw0∈ℋandw^0∈ℋ, there exists a unique solution of systems (32)–(34) such thatw·,t∈Cτ,∞;ℋ,ε·,s∈C0,t−τ;ℋ,ε·,s,t∈Ct−τ,t×[τ,∞;ℋ), and for anyε·,0∈DAwhereAis defined in (11), system (32) decays exponentially in the meaning that(35)w·,tℋ≤M0e−ω0t−τw0ℋ+c0+qL0M0Meω0τC1′2+q2e2qC2′2ωAε·,0ℋ·e−ω0t/2+e−ωt/2·eωτ,where L0 is defined in (54) and (55) and M0 and M are defined in (13) and (57), respectively.Proof.
Due to (21) and (23), for any ε·,t−τ∈ℋ determined by (33), there exists a unique solution of (34) such that(36)ε·,s,tℋ≤M1eω1s−t+τε·,t−τℋ,∀t−τ≤s≤t,t>τ.which together with (14) tell us that(37)ε·,s,tℋ≤MM1e−ω+ω1t−τeω1sε·,0ℋ,∀t−τ≤s≤t,t>τ.Thus, ε·,s,t∈Ct−τ,t×[τ,∞);ℋ.
Now for system (34), because of (18) and (19), the solution of (34) can be represented as [23](38)ε·,s,t=∑n=0∞anteλnsφn·,where(39)λn=−nπ2+2q+On−1,φnx=cosnπx+On−1.
As a result, we have that(40)ε1,t,t=∑n=0∞anteλntφn1,∀t>τ,εx,t,t=∑n=0∞anteλntφnx,∀0≤x≤1,t>τ.
Thus, by Cauchy–Schwartz inequality,(41)|ε1,t,t|≤∑n=0∞λnanteλnt−τ21/2·∑n=0∞λn−1eλnτ21/2,∀t>τ.
Similarly, we obtain that(42)|εx,t,t≤∑n=0∞λnanteλnt−τ21/2·∑n=0∞λn−1eλnτ21/2,∀0≤x≤1,t>τ.
Moreover, from the fact that(43)Aε·,t−τℋ=∑n=0∞λnanteλnt−τ21/2,∀t>τ,where A is defined in (11), there exists some positive constant C independent of t such that(44)ε1,t,t≤CAε·,t−τℋ,∀t>τ,εx,t,t≤Aε·,t−τℋ,∀0≤x≤1,t>τ,where(45)C=∑n=0∞λn−1eλnτ21/2.
Moreover, by (14) and C0-semigroup theory, the following inequality holds:(46)Aε·,t−τℋ≤Me−ωt−τAε·,0ℋ,for any t>τ, where Mandω are given in (13).
In conclusion, based on the equations from (37) to (43), we have that(47)ε1,t,t≤CMe−ωt−τAε·,0ℋ,∀t>τ,εx,t,t≤CMe−ωt−τAε·,0ℋ,∀0≤x≤1,t>τ.
From the second inequality of (44), it is shown that(48)∫01eq1−ηεη,t,tdη≤eqCMe−ωt−τAε·,0ℋ,∀t>τ,q>0.
On the other hand, the equivalent system of (32) can be illustrated easily by the transformation (2) as follows:(49)vtx,t=vxxx,t,0<x<1,t>τ,vx0,t=c0v0,t,t>τ,c0>0,vx1,t=−c0+qε1,t,t−qc0+q∫01eq1−ηεη,t,tdη,t>τ,q>0.
It can be written as follows:(50)ddtv·,t=A0v·,t+ℬ0c0+qε1,t,t+q∫01eq1−ηεη,t,tdη,where(51)A0φ=φ″,∀φ∈DA0,DA0=φ∈H10,1|φ′0=φ′1=0,(52)ℬ0=−δx−1.
A direct computation shows(53)ℬ0A0−1f=−∫011c0+ηfηdη,∀f∈ℋ.which means that ℬ0A0−1 is bounded.
For the Lyapunov function of system (46),(54)E0t=12∫01v2x,tdx.
The direct computation shows that(55)E˙0t=−c0v20,t−∫01vx2x,tdx,which means that(56)c0∫0Tv20,tdt≤E00.for any T>0. This inequality together with (50) illustrates that ℬ0 is admissible for eA0t. Therefore, there exists a unique solution of (46) such that v·,t∈Cτ,∞;ℋ. Since the operator A0 generates an exponentially stable C0-semigroup [25], it follows from (44) and (45) that [27, 28](57)∫τt/2eA0t/2−sℬ0c0+qε1,s,s+q∫01eq1−ηεη,s,sdηdsℋ≤c0+qL0ε1,·,·+q∫01eq1−ηεη,·,·dηL2τ,t/2≤c0+qL02∫τt/2ε21,s,s+q2∫01e2q1−ηε2η,s,sdηds1/2≤c0+qL0MC1′2+q2e2qC2′2/ωAε·,0ℋ,(58)∫t/2teA0t−sℬ0c0+qε1,s,s+q∫01eq1−ηεη,s,sdηdsℋ≤c0+q∫0teA0t−sℬ00⋄t/2ε1,s,s+q∫01eq1−ηεη,s,sdηℋ≤c0+qL0ε1,·,·+q∫01eq1−ηεη,·,·dηL2t/2,t≤c0+qL02∫t/2tε21,s,s+q2∫01e2q1−ηε2η,s,sdηds1/2.≤c0+qL0MeωτC1′2+q2e2qC2′2/ωAε·,0ℋ·e−ωt/2,for some constants L0>0 independent of(59)ε1,t,t+q∫01eq1−ηεη,t,tdη,where(60)u⋄τ0vt=ut0≤t≤τ0,vtt>τ0.
Furthermore, the solution of system (46) can be represented as(61)v·,t=eA0t−τv·,τ+c0+q∫τteA0t−s−τℬ0ε1,s,s+q∫01eq1−ηεη,s,sdηds=eA0t−τv·,τ+c0+qeA0t/2−τ∫τt/2eA0t/2−sℬ0ε1,s,s+q∫01eq1−ηεη,s,sdηds+c0+qe−A0τ∫t/2teA0t−sℬ0ε1,s,s+q∫01eq1−ηεη,s,sdηds.
SinceA0 generates an exponentially stable C0-semigroup, there exist two positive constants M0,ω0 such that(62)eA0tℋ≤M0e−ω0t,which together with (54), (55), and (58) show that(63)v·,tℋ≤M0e−ω0t−τv·,τℋ+c0+qL0M0Meω0τC1′2+q2e2qC2′2ωAε·,0ℋ·e−ω0t/2+e−ωt/2·eωτ,∀t>τ.
At last, since system (32) and (46) are equivalent to each other based on the transformation (2), naturally we have that(64)w·,tℋ≤M0e−ω0t−τw0ℋ+c0+qL0M0Meω0τC1′2+q2e2qC2′2ωAε·,0ℋ·e−ω0t/2+e−ωt/2·eωτ,∀t>τ,where L0 is defined in (54) and (55), and M0 and M are defined in (55) and (13), respectively.


## 2.1. Feedback Control via Backstepping Transform

System (1) is considered in the energy space ℋ=L20,1 and the input(or output) space U=Y=ℂ. The energy of system (1) is(2)Et=12∫01w2x,tdx.For the original system (1), make the invertible change of variable(3)vx,t=wx,t+c0+q∫0xeqx−ηwη,tdη,0≤x≤1,t≥0,c0>0,q>0,whose equivalent transformation is(4)wx,t=vx,t−c0+q∫0xe−c0x−ηvη,tdη,0≤x≤1,t≥0,c0>0,q>0.In order to convert (1) into the exponentially stable system [25](5)vtx,t=vxxx,t,0<x<1,t>0,vx0,t=c0v0,t,t≥0,c0>0,vx1,t=0,t≥0,the feedback control in (1) may be proposed as that(6)ut=−c0+qw1,t−qc0+q∫01eq1−ηwη,tdη,t≥0,c0>0,q>0,under which system (1) becomes that [8](7)wtx,t=wxxx,t,0<x<1,t>0,wx0,t=−qw0,t,t≥0,q>0,wx1,t=−c0+qw1,t−qc0+q∫01eq1−ηwη,tdη,t≥0,c0>0,q>0.Based on the exponential stability of system (5) and the equivalence between systems (5) and (7), we have the following theorem [8].Theorem 1.
For any initial valuew0∈ℋ, under feedback (6), the closed-loop system (7) is exponentially stable.


## 2.2. Design of Observer and Predictor

This section has been devoted to constructing the observer system and designing the predictor system. In the first step, the observer has been constructed to estimate the statewx,s,0≤x≤1,0≤s≤t−τ,t>τ from the known observation ys+τ,0≤s≤t−τ,t>τ stated as the following system:(8)w^sx,s=w^xxx,s,0<x<1,0<s<t−τ,t>τ,w^x0,s=−qys+τ+k1w^0,s−ys+τ,0≤s≤t−τ,t>τ,k1>0,w^x1,s=us,0≤s≤t−τ,t>τ,w^x,0=w^0x,0≤x≤1,where w^0 is the arbitrarily assigned initial value of the observer.Set the error for the observer system(9)εx,s=w^x,s−wx,s,0≤x≤1,0≤s≤t−τ,t>τ.Then, by (1) and (8), εx,s satisfies the following:(10)εsx,s=εxxx,s,0<x<1,0<s<t−τ,t>τ,εx0,s=k1ε0,s,0≤s≤t−τ,t>τ,k1>0,εx1,s=0,0≤s≤t−τ,t>τ,εx,0=w^0x−w0x,0≤x≤1.System (10) can be written as follows:(11)ddsε·,s=Aε·,s,where the operator A is defined as follows:(12)Aφ=φ″,∀φ∈DA,DA=φ∈H20,1|φ′0=k1φ0k1>0,φ′1=0.It is known thatA generates an exponentially stable C0-semigroup on ℋ such that(13)eAsℋ≤Me−ωs,∀s>0,for some positive constants M,ω ([25]). Thus, for any w0∈ℋ and w^0∈ℋ, there exists a unique solution of (12) such that(14)ε·,sℋ≤Me−ωsε·,0ℋ,∀0≤s≤t−τ,t>τ.In the second step, the predictor system has been designed for the statewx,s,0≤x≤1,t−τ≤s≤t,t>τ as follows:(15)w^sx,s,t=w^xxx,s,t,0<x<1,t−τ<s<t,t>τ,w^x0,s,t=−qw^0,s,t,t−τ≤s≤t,t>τ,q>0,w^x1,s,t=us,t−τ≤s≤t,t>τ,w^x,t−τ,t=w^x,t−τ,0≤x≤1,t>τ,Let the error for the predictor system(16)εx,s,t=w^x,s,t−wx,s,0≤x≤1,t−τ≤s≤t,t>τ.Then, by (1), (15), and (16), εx,s,t satisfies the following system:(17)εsx,s,t=εxxx,s,t,0<x<1,t−τ<s<t,t>τ,εx0,s,t=−qε0,s,t,t−τ≤s≤t,t>τ,q>0,εx1,s,t=0,t−τ≤s≤t,t>τ,εx,t−τ,t=εx,t−τ,0≤x≤1,t>τ.It can be expressed as that(18)ddsε·,s,t=ℬε·,s,t,where(19)ℬφ=φ″,∀φ∈Dℬ,Dℬ=φ∈H20,1||φ′0=−qφ0q>0,φ′1=0.Next, we define a new inner product equivalent as the energy in (2)(20)φ1,φ2ℋ1=∫01coshα1−xφ1xφ2xdx,for some given α>0 satisfying that 2qcoshα−αsinhα<0.Thus, for anyφ∈Dℬ, we have that(21)φ,φℋ1=∫01coshα1−xφx2dx.By simple computations of integration by parts, it results that(22)ℬφ,φℋ1=φ″,φℋ1=∫01coshα1−xφ″xφxdx,=qcoshα−αsinhα2φ02+∫01coshα1−xα2φx2/2−φ′x2dx,≤α2/2·φ,φℋ1.Thus, system (17) generates a strongly continuous semigroup eℬs which indicates that(23)eℬsℋ≤M1eω1s,∀s≥0,for some positive constants M1 and ω1 [26]. Therefore, the solution of system (17) can be represented as follows:(24)ε·,s,t=eℬs−t+τε·,t−τ,t−τ≤s≤t,t>τ,which together with (14) and (23) result that(25)ε·,t,tℋ≤MM1eω+ω1τε·,0ℋe−ωt,∀t>τ.The estimated state has naturally been chosen as follows:(26)w˜x,t=w^x,t,t,∀t>τ.assured by the theorem below.Theorem 2.
For∀t>τ, we have that(27)w·,t−w˜·,tℋ≤MM1eω+ω1τw^0−w0ℋe−ωt,where w^0 is the initial value of the observer (8), w0 is the initial value of the original system (1), and M,ω and M1,ω1 are constants in (13) and (23), respectively.


## 2.3. Stabilization of the Closed-Loop System

For the estimated state (24) and (7), we naturally obtain the estimated output feedback control law as follows:(28)u∗t=0,0≤t≤τ,−c0+qw˜1,t−qc0+q∫01eq1−ηw˜η,tdη,t>τ,c0>0,q>0.Then, the closed-loop system becomes a system of partial differential equations as given below:(29)wtx,t=wxxx,t,0<x<1,t>0,wx0,t=−qw0,t,t≥0,q>0,wx1,t=u∗t,t≥0,wx,0=w0x,0≤x≤1,(30)w^sx,s=w^xxx,s,0<x<1,0<s<t−τ,t>τ,w^x0,s=−qys+τ+k1w^0,s−ys+τ,0≤s≤t−τ,t>τ,q>0,k1>0,w^x1,s=u∗s,0≤s≤t−τ,t>τ,w^x,0=w^0x,0≤x≤1,(31)w^sx,s,t=w^xxx,s,t,0<x<1,t−τ<s<t,t>τ,w^x0,s,t=−qw^0,s,t,t−τ≤s≤t,t>τ,q>0,w^x1,s,t=u∗s,t−τ≤s≤t,t>τ,w^x,t−τ,t=w^x,t−τ,0≤x≤1,t>τ.which is equivalent to the system below:(32)wtx,t=wxxx,t,0<x<1,t>τ,wx0,t=−qw0,t,t>τ,q>0,wx1,t=−c0+qw1,t+ε1,t,t−qc0+q∫01eq1−ηwη,t+εη,t,tdη,t>τ,q>0,wx,0=w0x,0≤x≤1,(33)εsx,s=εxxx,s,0<x<1,0<s<t−τ,t>τ,εx0,s=k1ε0,s,0≤s≤t−τ,t>τ,k1>0,εx1,s=0,0≤s≤t−τ,t>τ,εx,0=ε0x,0≤x≤1,(34)εsx,s,t=εxxx,s,t,0<x<1,t−τ<s<t,t>τεx0,s,t=−qε0,s,t,t−τ≤s≤t,t>τ,q>0,εx1,s,t=0,t−τ≤s≤t,t>τ,εx,t−τ,t=εx,t−τ,0≤x≤1,t>τ.Theorem 3.
For∀t>τ, anyw0∈ℋandw^0∈ℋ, there exists a unique solution of systems (32)–(34) such thatw·,t∈Cτ,∞;ℋ,ε·,s∈C0,t−τ;ℋ,ε·,s,t∈Ct−τ,t×[τ,∞;ℋ), and for anyε·,0∈DAwhereAis defined in (11), system (32) decays exponentially in the meaning that(35)w·,tℋ≤M0e−ω0t−τw0ℋ+c0+qL0M0Meω0τC1′2+q2e2qC2′2ωAε·,0ℋ·e−ω0t/2+e−ωt/2·eωτ,where L0 is defined in (54) and (55) and M0 and M are defined in (13) and (57), respectively.Proof.
Due to (21) and (23), for any ε·,t−τ∈ℋ determined by (33), there exists a unique solution of (34) such that(36)ε·,s,tℋ≤M1eω1s−t+τε·,t−τℋ,∀t−τ≤s≤t,t>τ.which together with (14) tell us that(37)ε·,s,tℋ≤MM1e−ω+ω1t−τeω1sε·,0ℋ,∀t−τ≤s≤t,t>τ.Thus, ε·,s,t∈Ct−τ,t×[τ,∞);ℋ.
Now for system (34), because of (18) and (19), the solution of (34) can be represented as [23](38)ε·,s,t=∑n=0∞anteλnsφn·,where(39)λn=−nπ2+2q+On−1,φnx=cosnπx+On−1.
As a result, we have that(40)ε1,t,t=∑n=0∞anteλntφn1,∀t>τ,εx,t,t=∑n=0∞anteλntφnx,∀0≤x≤1,t>τ.
Thus, by Cauchy–Schwartz inequality,(41)|ε1,t,t|≤∑n=0∞λnanteλnt−τ21/2·∑n=0∞λn−1eλnτ21/2,∀t>τ.
Similarly, we obtain that(42)|εx,t,t≤∑n=0∞λnanteλnt−τ21/2·∑n=0∞λn−1eλnτ21/2,∀0≤x≤1,t>τ.
Moreover, from the fact that(43)Aε·,t−τℋ=∑n=0∞λnanteλnt−τ21/2,∀t>τ,where A is defined in (11), there exists some positive constant C independent of t such that(44)ε1,t,t≤CAε·,t−τℋ,∀t>τ,εx,t,t≤Aε·,t−τℋ,∀0≤x≤1,t>τ,where(45)C=∑n=0∞λn−1eλnτ21/2.
Moreover, by (14) and C0-semigroup theory, the following inequality holds:(46)Aε·,t−τℋ≤Me−ωt−τAε·,0ℋ,for any t>τ, where Mandω are given in (13).
In conclusion, based on the equations from (37) to (43), we have that(47)ε1,t,t≤CMe−ωt−τAε·,0ℋ,∀t>τ,εx,t,t≤CMe−ωt−τAε·,0ℋ,∀0≤x≤1,t>τ.
From the second inequality of (44), it is shown that(48)∫01eq1−ηεη,t,tdη≤eqCMe−ωt−τAε·,0ℋ,∀t>τ,q>0.
On the other hand, the equivalent system of (32) can be illustrated easily by the transformation (2) as follows:(49)vtx,t=vxxx,t,0<x<1,t>τ,vx0,t=c0v0,t,t>τ,c0>0,vx1,t=−c0+qε1,t,t−qc0+q∫01eq1−ηεη,t,tdη,t>τ,q>0.
It can be written as follows:(50)ddtv·,t=A0v·,t+ℬ0c0+qε1,t,t+q∫01eq1−ηεη,t,tdη,where(51)A0φ=φ″,∀φ∈DA0,DA0=φ∈H10,1|φ′0=φ′1=0,(52)ℬ0=−δx−1.
A direct computation shows(53)ℬ0A0−1f=−∫011c0+ηfηdη,∀f∈ℋ.which means that ℬ0A0−1 is bounded.
For the Lyapunov function of system (46),(54)E0t=12∫01v2x,tdx.
The direct computation shows that(55)E˙0t=−c0v20,t−∫01vx2x,tdx,which means that(56)c0∫0Tv20,tdt≤E00.for any T>0. This inequality together with (50) illustrates that ℬ0 is admissible for eA0t. Therefore, there exists a unique solution of (46) such that v·,t∈Cτ,∞;ℋ. Since the operator A0 generates an exponentially stable C0-semigroup [25], it follows from (44) and (45) that [27, 28](57)∫τt/2eA0t/2−sℬ0c0+qε1,s,s+q∫01eq1−ηεη,s,sdηdsℋ≤c0+qL0ε1,·,·+q∫01eq1−ηεη,·,·dηL2τ,t/2≤c0+qL02∫τt/2ε21,s,s+q2∫01e2q1−ηε2η,s,sdηds1/2≤c0+qL0MC1′2+q2e2qC2′2/ωAε·,0ℋ,(58)∫t/2teA0t−sℬ0c0+qε1,s,s+q∫01eq1−ηεη,s,sdηdsℋ≤c0+q∫0teA0t−sℬ00⋄t/2ε1,s,s+q∫01eq1−ηεη,s,sdηℋ≤c0+qL0ε1,·,·+q∫01eq1−ηεη,·,·dηL2t/2,t≤c0+qL02∫t/2tε21,s,s+q2∫01e2q1−ηε2η,s,sdηds1/2.≤c0+qL0MeωτC1′2+q2e2qC2′2/ωAε·,0ℋ·e−ωt/2,for some constants L0>0 independent of(59)ε1,t,t+q∫01eq1−ηεη,t,tdη,where(60)u⋄τ0vt=ut0≤t≤τ0,vtt>τ0.
Furthermore, the solution of system (46) can be represented as(61)v·,t=eA0t−τv·,τ+c0+q∫τteA0t−s−τℬ0ε1,s,s+q∫01eq1−ηεη,s,sdηds=eA0t−τv·,τ+c0+qeA0t/2−τ∫τt/2eA0t/2−sℬ0ε1,s,s+q∫01eq1−ηεη,s,sdηds+c0+qe−A0τ∫t/2teA0t−sℬ0ε1,s,s+q∫01eq1−ηεη,s,sdηds.
SinceA0 generates an exponentially stable C0-semigroup, there exist two positive constants M0,ω0 such that(62)eA0tℋ≤M0e−ω0t,which together with (54), (55), and (58) show that(63)v·,tℋ≤M0e−ω0t−τv·,τℋ+c0+qL0M0Meω0τC1′2+q2e2qC2′2ωAε·,0ℋ·e−ω0t/2+e−ωt/2·eωτ,∀t>τ.
At last, since system (32) and (46) are equivalent to each other based on the transformation (2), naturally we have that(64)w·,tℋ≤M0e−ω0t−τw0ℋ+c0+qL0M0Meω0τC1′2+q2e2qC2′2ωAε·,0ℋ·e−ω0t/2+e−ωt/2·eωτ,∀t>τ,where L0 is defined in (54) and (55), and M0 and M are defined in (55) and (13), respectively.


## 3. Appendix: Stabilization with Dirichlet Boundary Control

For the system under the Dirichlet boundary control,(65)wtx,t=wxxx,t,0<x<1,t>0,wx0,t=−qw0,t,t≥0,w1,t=ut,t≥0,yt=w0,t−τ,t>τ,wx,0=w0x,0≤x≤1.By the similar design of estimated feedback control as that in the previous section, it can be stabilized for any given time delay. Firstly, we can apply the backstepping transformation to obtain the stabilizing feedback for system (62). And then, we design the observer and predictor systems. The output feedback control based on the estimated state can be shown to stabilize system (62) for any given time delay.


## 4. Numerical Simulation

In this section, numerical simulations have been given to show the effectiveness of the stabilizing output feedback controller. The space step is 0.025, and the time step is10−4. For the closed-loop system (29)–(31) with the Neumann boundary control, initial values have been chosen as that(66)w0x=x−1,ε0x=x2−2x−2,and the parameters have been chosen as that(67)k1=1,q=1,c0=1.while the time delay is τ=0.1. Then, it is found that the state of the closed-loop system is almost rest after a certain time shown in Figure 1.Figure 1
State of the closed-loop system (29)–(31) with Neumann boundary control.For the closed-loop system with the Dirichlet boundary control, we choose the initial values as follows:(68)w0x=x2−2x,ε0x=cosπx2,and the parameter has been chosen asq=1 while the time delay is τ=0.1. It is shown that the state of the closed-loop system is evidently stable shown in Figure 2.Figure 2
State of the closed-loop system with Dirichlet boundary control.


## 5. Conclusion

In this paper, we consider an unstable heat equation where there exists instability in the Neumann or Dirichlet boundary conditions. For the unstable Neumann boundary condition, we stabilize the heat equation where the observation signal suffers from a given time delay by designing the estimated state feedback controller based on the construction of the observer and predictor systems. Numerical simulations show the effectiveness of the stabilized controller. For the unstable Dirichlet boundary condition, a similar controller may be given to make the originally unstable system stable exponentially which can be shown in the simulation results of this paper. The future research direction may be the stabilization of the high dimensional or abstract parameter distributed systems with unstable effects [29, 30].



---

*Source: 6973646-2022-03-22.xml*