Next Article in Journal
Fuzzy Portfolio Selection in COVID-19 Spreading Period Using Fuzzy Goal Programming Model
Previous Article in Journal
Multi-Attribute Online Decision-Making Driven by Opinion Mining
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Boundary Stabilization of Heat Equation with Multi-Point Heat Source

School of Mathematics and Statistics, Shandong Normal University, Jinan 250014, China
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(8), 834; https://doi.org/10.3390/math9080834
Submission received: 25 March 2021 / Revised: 8 April 2021 / Accepted: 9 April 2021 / Published: 12 April 2021

Abstract

:
In this paper, we consider boundary stabilization problem of heat equation with multi-point heat source. Firstly, a state feedback controller is designed mainly by backstepping approach. Under the designed state controller, the exponential stability of closed-loop system is guaranteed. Then, an observer-based output feedback controller is proposed. We prove the exponential stability of resulting closed-loop system using operator semigroup theory. Finally, the designed state and output feedback controllers are effective via some numerical simulations.

1. Introduction

The parabolic partial differential equations (PDEs) have played a crucial role in the real-world practice over the past several decades. At early stages, Masi et al. [1] investigated the microscopic fluctuations of interacting particle systems on a lattice by the nonlinear parabolic PDE. In [2], Price et al. showed the analysis and enhancement on special types of images based on parabolic PDEs. Afterwards, under measuring the temperature changes of ovens’ processing area on different points, the heating capability of forced convection reflow ovens and the effect of the ovens’ construction on it, were considered in [3]. By the experimental data, the parabolic PDE was used to describe the specific characteristics. Besides, Pardo et al. [4] applied the three-dimensional parabolic PDE to the non-destructive evaluation of minefields. Recently, Gao et al. [5] proposed a method of parabolic equation modeling to address the problem of predicting electromagnetic wave propagation caused by atmospheric dust and rough sea surfaces in the maritime environment. Regarding the study of the parabolic PDEs, there are some other results and references [6,7,8,9,10] therein.
In the control engineering filed, many scholars have been keen on boundary stabilization problems of parabolic PDEs during the past decades. For boundary control problem of a class of unstable scalar parabolic PDEs, ref. [11] employed a gradient-based optimization method to parameterize the feedback kernel as a second-order polynomial, and made the optimized kernel generate closed-loop stability with restricting on the kernel coefficients. Liu [12] proposed the successive approximation to address the boundary stabilization problem of unstable parabolic PDE. Wu et al. [13] made use of fuzzy control approach to solve the stabilization problem of nonlinear parabolic PDEs. Vazquez and Krstic [14] adopted finite-dimension feedback linearization method to transform parabolic PDEs with Volterra nonlinearities into a stable system by the Volterra series nonlinear operators. Actually, the proposed finite-dimension feedback linearization method was infinite dimensional extension of the backstepping approach. It is well known that backstepping transformation is mainstream method to deal with boundary stabilization problems of parabolic PDEs, such as, linear parabolic PDEs [15,16,17,18,19], nonlinear parabolic PDEs [20,21,22], quasi-linear parabolic PDEs [23], coupled parabolic PDEs and ODE [24,25,26,27].
To the best of our knowledge, heat equation is an extremely classic class of parabolic PDEs. The boundary control results of heat equation with unstable term or source term have been proposed and references [19,28,29,30,31,32,33,34] therein. Krstic first proposed a separation principle result based on the passive and swapping identifiers and combined backstepping method to solve adaptive boundary control problem for heat equation with term λ u ( x , t ) in [33], where λ is unknown parameter. Afterward, based on expressly parametrized control formula, adaptive controllers for parabolic PDEs with the term λ u ( x , t ) were designed in [31], where λ is an unknown constant parameter. Additionally, Baccoli et al. [28] employed backstepping method to represent boundary stabilization result of the coupled reaction-diffusion processes with term Λ Q ( x , t ) , where Λ R n × n is a real-valued square matrix. In [19], heat equation system is considered by
Γ t ( ξ , t ) = Γ ξ ξ ( ξ , t ) + τ 0 Γ ( ξ 0 , t ) , 0 < ξ < 1 , t ( 0 , + ) , Γ ξ ( 0 , t ) = 0 , t 0 , Γ ( 1 , t ) = U ^ ( t ) , t 0 , Γ ( ξ , 0 ) = Γ 0 ( ξ ) , 0 ξ 1 ,
where U ^ ( t ) is boundary control, Γ 0 is initial value, τ 0 is a constant, ξ 0 ( 0 , 1 ) . Zhou and Guo designed a state feedback controller based on backstepping approach for problem (1). Motivated by [19], we consider extension of heat system (1)
Γ t ( ξ , t ) = Γ ξ ξ ( ξ , t ) + i = 1 Z τ i Γ ( ξ i , t ) , 0 < ξ < 1 , t ( 0 , + ) , Γ ξ ( 0 , t ) = 0 , t 0 , Γ ( 1 , t ) = U ^ ( t ) , t 0 , Y = { Γ ( ξ 1 , t ) , Γ ( ξ 2 , t ) , . . . , Γ ( ξ Z , t ) } , Γ ( ξ , 0 ) = Γ 0 ( ξ ) , 0 ξ 1 ,
where 0 < ξ i < 1 , i = 1 , 2 , , Z , Z [ 2 , + ) is a positive integer, τ i ( i = 1 , 2 , , Z ) is a constant parameter, U ^ ( t ) is control, Γ 0 is initial value, i = 1 Z τ i Γ ( ξ i , t ) is non-local term. As a matter of convenience, let τ = i = 1 Z τ i . In our paper, owing to the existence of our multi-term τ , normal backstepping transformation is no longer applicable. A new backstepping transformation is constructed in our paper. Unlike the unstable terms mentioned above, term i = 1 Z τ i Γ ( ξ i , t ) in system (1) is intermediate points related. In addition, i = 1 Z τ i Γ ( ξ i , t ) can be seen as 0 1 τ i δ ( s ξ i ) Γ ( s , t ) d s (a nice explanation in [19]). When i = 1 , system (2) becomes system (1). Actually, [35] adopted backstepping method to solve the output feedback problem of transport equation with non-local term μ u ( x 0 ) . Besides that there are other results about hyperbolic PDEs with integral term by backstepping approach in [36,37].
The rest of paper is organized as follows. In Section 2, we design a state feedback controller by backstepping approach to stabilize system (2). In Section 3, an output feedback controller is constructed. The resulting closed-loop system is proved to be exponentially stable. Section 4 illustrates the effectiveness of the proposed controller on the basis of some simulations. Finally, the concluding remarks are introduced in Section 5.

2. The State Feedback Controller Design by Backstepping

In this section, we introduce the backstepping transformation as below
P ( ξ , t ) = Γ ( ξ , t ) i = 1 Z a i ( ξ ) 0 ξ i b i ( s ) Γ ( s , t ) d s 0 ξ c ( ξ , s ) Γ ( s , t ) d s ,
where a i , b i ( i = 1 , 2 , , Z ) and c are kernels and to be determined later.
Under transformation (3), we convert system (2) into the following target system
P t ( ξ , t ) = P ξ ξ ( ξ , t ) + σ P ( ξ , t ) , 0 < ξ < 1 , t ( 0 , + ) , P ξ ( 0 , t ) = 0 , t 0 , P ( 1 , t ) = 0 , t 0 ,
where σ ( , 0 ] is any given number.
Now, we prove the kernels are unique. Firstly, define the triangular domain Λ
Λ : = ( ξ , s ) R 2 | ( ξ , s ) [ 0 , 1 ] × [ 0 , ξ ] .
For transformation (3), we take derivative of it with regard to ξ , there is
P ξ ( ξ , t ) = Γ ξ ( ξ , t ) i = 1 Z a i ( ξ ) 0 ξ i b i ( s ) Γ ( s , t ) d s c ( ξ , ξ ) Γ ( s , t ) 0 ξ c ξ ( ξ , t ) Γ ( s , t ) d s .
Take derivative of Equation (5) with regard to ξ
P ξ ξ ( ξ , t ) = Γ ξ ξ ( ξ , t ) i = 1 Z a i ( ξ ) 0 ξ i b i ( s ) Γ ( s , t ) d s d c ( ξ , ξ ) d ξ Γ ( ξ , t ) c ( ξ , s ) Γ ξ ξ ( ξ , t ) c ξ ( ξ , ξ ) Γ ( ξ , t ) 0 ξ c ( ξ , s ) Γ ( s , t ) d s .
Next, we calculate the derivation of (3) with regard to t
P t ( ξ , t ) = Γ t ( ξ , t ) i = 1 Z a i ( ξ ) 0 ξ i b i ( s ) Γ t ( s , t ) d s 0 ξ c ( ξ , s ) Γ t ( s , t ) d s = Γ ξ ξ ( ξ , t ) + i = 1 Z τ i Γ ( ξ i , t ) i = 1 Z a i ( ξ ) [ b i ( ξ i ) Γ ξ ( ξ i , t ) b i ( 0 ) Γ ξ ( 0 , t ) b i ( ξ i ) Γ ( ξ i , t ) + b i ( 0 ) Γ ( 0 , t ) + 0 ξ i b ( s ) Γ ( s , t ) d s + 0 ξ i b i ( s ) d s j = 1 Z τ j Γ ( ξ j , t ) ] c ( ξ , ξ ) Γ ξ ( ξ , t ) + c s ( ξ , ξ ) Γ ( ξ , t ) c s ( ξ , 0 ) Γ ( 0 , t ) 0 ξ c s s ( ξ , s ) Γ ( s , t ) d s 0 ξ c ( ξ , s ) d s i = 1 Z τ i Γ ( ξ i , t ) .
According to P t ( ξ , t ) = P ξ ξ ( ξ , t ) + σ P ( ξ , t ) in (4), one has
i = 1 Z 0 ξ i ( a i ( ξ ) b i ( s ) a i ( ξ ) b i ( s ) + σ a i ( ξ ) b i ( s ) ) Γ ( s , t ) d s + 0 ξ ( c ξ ξ ( ξ , s ) c s s ( ξ , s ) + σ c ( ξ , s ) ) Γ ( s , t ) d s + ( 2 d d ξ c ( ξ , ξ ) σ ) Γ ( ξ , t ) + i = 1 Z Γ ( ξ i ) [ τ i + a i ( ξ ) b i ( ξ i ) τ i j = 1 Z 0 ξ j a j ( ξ ) b j ( s ) d s τ i 0 ξ c ( ξ , s ) d s ] i = 1 Z a i ( ξ ) b i ( ξ i ) Γ ξ ( ξ i ) i = 1 Z a i ( ξ ) b i ( 0 ) + c s s ( ξ , 0 ) Γ ( 0 , t ) = 0 .
In addition, by P ξ ( 0 , t ) = 0 and Γ ξ ( 0 , t ) = 0 , we get
P ξ ( 0 , t ) = Γ ξ ( 0 , t ) i = 1 Z a i ( 0 ) 0 ξ i b i ( s ) Γ ( s , t ) d s c ( 0 , 0 ) Γ ( 0 , t ) .
By the (8) and (9), there is
a i ( ξ ) b i ( s ) a i ( ξ ) b i ( s ) + σ a i ( ξ ) b i ( s ) = 0 , 2 d d ξ c ( ξ , ξ ) σ = 0 , c ξ ξ ( ξ , s ) c s s ( ξ , s ) + σ c ( ξ , s ) = 0 , i = 1 Z a i ( ξ ) b i ( 0 ) + c s ( ξ , 0 ) = 0 , τ i + a i ( ξ ) b i ( ξ i ) τ i j = 1 Z 0 ξ j a j ( ξ ) b j ( s ) d s τ i 0 ξ c ( ξ , s ) d s = 0 , b i ( ξ i ) = a i ( 0 ) = 0 .
For our main result, the following assumption is given by
S 1 = 2 + i = 1 Z τ i ξ i 2 0 , if τ = 0 ; S 2 = i = 1 Z τ i cosh ( τ ξ i ) 0 , if τ > 0 ; S 3 = i = 1 Z τ i cos ( τ ξ i ) 0 , if τ < 0 .
Lemma 1.
Under assumption (11), for (10), there exist classical solutions a i ( · ) C 2 ( [ 0 , 1 ] ) , b i ( · ) C 2 ( [ 0 , ξ i ] ) , and c ( · , · ) C 2 ( Λ ) .
Proof of Lemma 1.
For convenience, we split (10) into
a i ( ξ ) b i ( s ) a i ( ξ ) b i ( s ) + σ a i ( ξ ) b i ( s ) = 0 , τ i + a i ( ξ ) b i ( ξ i ) τ i j = 1 Z 0 ξ i a j ( ξ ) b j ( s ) d s τ i 0 ξ c ( ξ , s ) d s = 0 , b i ( x i ) = a i ( 0 ) = 0 ,
and
2 d d ξ c ( ξ , ξ ) σ = 0 , c ξ ξ ( ξ , s ) c s s ( ξ , s ) + σ c ( ξ , s ) = 0 , i = 1 Z a i ( ξ ) b i ( 0 ) + c s ( ξ , 0 ) = 0 .
Firstly, we prove that for (12), there exist classical solutions a i ( · ) C 2 ( [ 0 , 1 ] ) , b i ( · ) C 2 ( [ 0 , ξ i ] ) . According to the first equation of (12), one has
a i ( ξ ) + σ a i ( ξ ) a i ( ξ ) = a i ( ξ ) β 2 a i ( ξ ) a i ( ξ ) = b i ( ξ ) b i ( ξ ) = τ ,
where σ = β 2 . For computing (12) conveniently, we make | τ | = λ 2 . The proof is split the following three steps.
Step 1: τ = 0 .
There are the following solutions for Equation (14)
a i ( ξ ) = η i c h β ξ , b i ( s ) = ζ i ( s ξ i ) .
From the second of (12), one gives
a i ( ξ ) b i ( ξ i ) τ i = a 1 ( ξ ) b 1 ( ξ 1 ) τ 1 ,
where i = 2 , 3 , , Z . Combining (15) and (16), we have
η i ζ i = η 1 ζ 1 τ i τ 1 .
From (13), one has
d 2 d ξ 2 0 ξ c ( ξ , s ) d s = d d ξ [ 0 ξ c ξ ( ξ , s ) d s + c ( ξ , ξ ) ] = 0 ξ c ξ ( ξ , s ) d s + c ξ ( ξ , ξ ) + d d ξ c ( ξ , ξ ) = i = 1 Z a i ( ξ ) b i ( 0 ) σ 0 ξ c ( ξ , s ) d s + σ .
By (15) and (17), there is
i = 1 Z a i ( ξ ) b i ( 0 ) = η 1 ζ 1 τ 1 i = 1 Z τ i c h β ξ = 0 .
Then,
d 2 d ξ 2 0 ξ c ( ξ , s ) d s = α 0 ξ c ( ξ , s ) d s + σ .
For Equation (20), we obtain
0 ξ c ( ξ , s ) d s 1 = cosh β ξ .
Now, let
F i ( ξ ) = τ i + a i ( ξ ) b i ( ξ i ) τ i j = 1 Z 0 ξ j a j ( ξ ) b j ( s ) d s τ i 0 ξ c ( ξ , s ) d s .
By simple calculation, we can get
F i ( ξ ) = a i ( ξ ) b i ( ξ i ) τ i j = 1 Z 0 ξ j a j ( ξ ) b j ( s ) d s τ i 0 ξ c ξ ( ξ , s ) d s σ 2 τ i ξ ,
and
F i ( ξ ) = a i ( ξ ) b i ( ξ i ) τ i j = 1 Z 0 ξ j a j ( ξ ) b j ( s ) d s τ i 0 ξ c ξ ξ ( ξ , s ) d s τ i c ξ ( ξ , s ) α 2 τ i .
According to (14), (15), (17) and (24), one has F i ( ξ ) = 0 . That means F i ( ξ ) is constant. Since a i ( 0 ) = 0 , there is F i ( 0 ) = 0 , which means F i ( ξ ) is constant. Suppose that F i ( ξ ) = 0 , when ξ = 0 , we get
F i ( 0 ) = τ i + a i ( 0 ) b i ( ξ i ) τ i j = 1 Z 0 ξ i a j ( 0 ) b j ( s ) d s = 0 .
From (15), (17), (25), and the first equation of assumption (11), one has
η i ζ i = 2 τ 1 2 + j = 1 Z τ j ξ j 2 = 2 τ 1 S 1 .
Hence,
a i ( ξ ) b i ( s ) = 2 τ 1 ( s ξ i ) cosh β ξ 2 + j = 1 Z τ j ξ j 2 = 2 τ 1 ( ξ i s ) cosh β ξ S 1 .
Step 2: τ > 0 .
By (14), one has
a i ( ξ ) = η i cosh ( λ 2 + β 2 ξ ) , b i ( s ) = ζ i ( sinh ( λ s ) tanh ( λ ξ i ) cosh ( λ s ) ) .
Similar calculation to step 1, we get
η i ζ i = η 1 ζ 1 τ i cosh ( λ ξ i ) τ 1 cosh ( λ ξ 1 ) ,
and
0 ξ c ( ξ , s ) d s 1 = i = 1 Z a i ( ξ ) b i ( 0 ) λ 2 .
Substituting (14), (28), (29), (30) into (24), similarly, there is F i ( ξ ) = 0 . By the same some steps as step 1, we can also deduce F i ( ξ ) is constant. Supposed that F i ( ξ ) = 0 , when ξ = 0 , it gives
F i ( 0 ) = τ i + a i ( 0 ) b i ( ξ i ) τ i j = 1 Z 0 ξ j a j ( 0 ) b j ( y ) d y = τ i + 4 η 1 ζ 1 τ i λ τ 1 cosh ( λ ξ 1 ) τ i j = 1 Z η 1 ζ 1 τ j cosh ( λ ξ j ) τ 1 cosh ( λ ξ 1 ) ( 4 cosh ( λ ξ j ) 1 ) = τ i + τ i j = 1 Z η 1 ζ 1 τ j cosh ( λ ξ j ) τ 1 cosh ( λ ξ 1 ) = 0 .
From (31) and the second equation of assumption (11), one has
η 1 ζ 1 = λ τ 1 cosh ( λ ξ 1 ) j = 1 N τ j cosh ( λ ξ j ) = λ τ 1 cosh ( λ ξ 1 ) S 2 .
By (29) and (32), it can be obtained that
η i ζ i = λ τ i cosh ( λ ξ i ) S 2 .
Hence,
a i ( ξ ) b i ( s ) = 1 S 2 λ τ i cosh ( λ ξ i ) cosh ( λ 2 + β 2 ξ ) ( sinh ( λ s ) tanh ( λ ξ i ) cosh ( λ s ) ) .
Step 3: τ < 0 .
There is
a i ( ξ ) + σ a i ( ξ ) a i ( ξ ) = a i ( ξ ) β 2 a i ( ξ ) a i ( ξ ) = b i ( ξ ) b i ( ξ ) = λ 2 .
By (35), one has
a i ( ξ ) = η i cos ( λ 2 β 2 ξ ) , if λ 2 β 2 > 0 , a i ( ξ ) = η i , if λ 2 β 2 = 0 , a i ( ξ ) = η i cosh β ( β 2 λ 2 ξ ) , if λ 2 β 2 < 0 , b i ( s ) = ζ i sin ( λ ( ξ i s ) ) .
Similar to step 1, we can get
η i ζ i = η 1 ζ 1 τ i τ 1 ,
and
0 ξ c ( ξ , s ) d s 1 = i = 1 Z a i ( ξ ) b i ( 0 ) λ 2 .
Substituting (14), (36), (37), (38) into (24), we can still ensure F i ( ξ ) = 0 and F i ( ξ ) is constant. Supposed that F i ( ξ ) = 0 , when ξ = 0 , it gives
F i ( 0 ) = τ i + a i ( 0 ) b i ( ξ i ) τ i j = 1 Z 0 ξ i a j ( 0 ) b j ( s ) d s = τ i + η 1 ζ 1 τ i λ τ 1 + τ i j = 1 Z η 1 ζ 1 τ j λ τ 1 ( 1 cos ( λ ξ j ) ) = 0 .
From (39) and the last equation of assumption (11), there is
η 1 ζ 1 = λ τ 1 j = 1 Z τ j cos ( λ ξ j ) = λ τ 1 S 3 ,
which means
η i ζ i = λ τ i S 3 .
Hence, we have
a i ( ξ ) b i ( s ) = λ τ i cos ( λ 2 β 2 ξ ) sin ( λ ( ξ i s ) ) S 3 , if λ 2 β 2 > 0 , a i ( ξ ) b i ( s ) = λ τ i sin ( λ ( ξ i s ) ) S 3 , if λ 2 β 2 = 0 , a i ( ξ ) b i ( s ) = λ τ i cosh ( β 2 λ 2 ξ ) sin ( λ ( ξ i s ) ) S 3 , if λ 2 β 2 < 0 .
Now, we show that for (13), there exists classical solution c ( · , · ) C 2 ( Λ ) . Let H ( ξ ) = i = 1 Z a i ( ξ ) b i ( 0 ) , and
G ( x ^ , y ^ ) = c ( ξ , s ) , govern by x ^ = ξ + s , y ^ = ξ s .
After a simple calculation, G ( x ^ , y ^ ) can be written as
G ( x ^ , y ^ ) = σ ( x ^ + y ^ ) 4 0 y ^ H ( s ) d s σ 4 0 y ^ 0 α G ( α , s ) d s d α σ 4 0 y ^ s x ^ G ( α , s ) d α d s .
Next, we only need to show that there exists a unique solution G ( · , · ) C ( [ 0 , 2 ] × [ 0 , 1 ] ) for (44) under assumption (11). Define norm
G = max ( x ^ , y ^ ) [ 0 , 2 ] × [ 0 , 1 ] e a ^ x ^ e b ^ y ^ | G ( x ^ , y ^ ) | ,
where σ 2 a ^ b ^ 1 . It is easy to see that norm G is completed normed space. Let
T G = σ ( x ^ + y ^ ) 4 0 y ^ H ( s ) d s σ 4 0 y ^ 0 α G ( τ , s ) d s d α σ 4 0 y ^ s x ^ G ( τ , s ) d α d s .
Thus,
| T G 1 ( x ^ , y ^ ) T G 2 ( x ^ , y ^ ) | = | σ 4 0 y ^ 0 α | G 1 ( α , s ) G 2 ( α , s ) | d s d α σ 4 0 y ^ s x ^ | G 1 ( α , s ) G 2 ( α , s ) | d α d s | σ 2 0 y ^ 0 x ^ | G 1 ( τ , s ) G 2 ( α , s ) | d s d α = σ 2 0 y ^ 0 x ^ | G 1 ( τ , s ) G 2 ( α , s ) | d s d α = σ 2 0 y ^ 0 x ^ e a ^ s e b ^ τ e a ^ s e b ^ τ | G 1 ( τ , s ) G 2 ( α , s ) | d s d α σ 2 0 y ^ 0 x ^ e a ^ s e b ^ α e a ^ s e b ^ α d s d α G 1 ( x ^ , y ^ ) G 1 ( x ^ , y ^ ) σ 2 a ^ b ^ e a ^ s e b ^ α G 1 ( x ^ , y ^ ) G 1 ( x ^ , y ^ ) .
Hence, we have e a ^ x ^ e b ^ y ^ | T G 1 ( x ^ , y ^ ) T G 2 ( x ^ , y ^ ) | σ 2 a ^ b ^ G 1 ( x ^ , y ^ ) G 1 ( x ^ , y ^ ) , then
T G 1 ( x ^ , y ^ ) T G 2 ( x ^ , y ^ ) σ 2 a ^ b ^ G 1 ( x ^ , y ^ ) G 1 ( x ^ , y ^ ) .
By contraction mapping principle, unique solution G ( · , · ) C ( [ 0 , 2 ] × [ 0 , 1 ] ) for (44) holds. □
Remark 1.
From the above proof, we cannot explain that there are unique solutions a i ( · ) C 2 ( [ 0 , 1 ] ) , b i ( · ) C 2 ( [ 0 , ξ i ] ) for Equation (12). However, a i ( ξ ) b i ( s ) is unique, which indicates transformation (3) is unique.
Next, we show that transformation (3) is reversible. Assume that inverse backstepping transformation is governed by
Γ ( ξ , t ) = P ( ξ , t ) + i = 1 Z d i ( ξ ) 0 ξ i h i ( s ) P ( s , t ) d s + 0 ξ k ( ξ , s ) P ( s , t ) d s .
Analogously, from (6) to (7), we get
Γ ξ ξ ( ξ , t ) = P ξ ξ ( ξ , t ) + i = 1 Z d i ( ξ ) 0 ξ i h i ( s ) Γ ( s , t ) d s d d ξ ( k ( ξ , ξ ) ) Γ ( ξ , t ) + k ( ξ , ξ ) Γ ξ ( ξ , t ) k ξ ( ξ , ξ ) Γ ξ ( ξ , t ) 0 ξ k ξ ξ ( ξ , s ) Γ ( s , t ) d s ,
and
Γ t ( ξ , t ) = P t ( ξ , t ) + σ P ( ξ , t ) + i = 1 Z d i ( ξ ) 0 ξ i h i ( s ) ( P s s ( s , t ) + σ P ( s , t ) ) d s 0 ξ k ( ξ , s ) ( P s s ( s , t ) + σ P ( s , t ) ) d s = P t ( ξ , t ) + σ P ( ξ , t ) + i = 1 Z d i ( ξ ) [ h i ( ξ i ) P ξ ( ξ i , t ) h i ( ξ i ) P ξ ( ξ i , t ) h i ( 0 ) P ξ ( 0 , t ) + 0 ξ i h i ( s ) P ( s , t ) d s + σ 0 ξ i h i ( s ) P ( s , t ) d s ] + k ( ξ , ξ ) P ξ ( ξ , t ) k s ( ξ , ξ ) P ( ξ , t ) + k s ( ξ , 0 ) P ( 0 , t ) + 0 ξ k s s ( ξ , s ) P ( s , t ) d s + σ 0 ξ k ( ξ , s ) P ( s , t ) d s .
Owing to Γ t ( ξ , t ) = Γ ξ ξ ( ξ , t ) + i = 1 Z τ i Γ ( ξ i , t ) , one has
i = 1 Z [ 0 ξ i ( d i ( ξ ) h i ( s ) + d i ( ξ ) h i ( s ) + σ d i ( ξ ) h i ( s ) j = 1 Z τ i d j ( ξ i ) h j ( s ) τ i k ( ξ i , s ) ) P ( s , t ) d s ] + 0 ξ ( k s s ( ξ , s ) k ξ ξ ( ξ , s ) + σ k ( ξ , s ) ) P ( s , t ) d s ( 2 d d ξ k ( ξ , ξ ) σ ) P ( ξ , t ) + i = 1 Z P ( ξ i ) ( τ i + d i ( ξ ) h i ( ξ i ) ) + ( i = 1 Z d i ( ξ ) h i ( 0 ) + k s ( ξ , 0 ) ) P ( 0 , t ) + i = 1 Z d i ( ξ ) h i ( ξ i ) P ξ ( ξ i , t ) = 0 .
Additionally,
Γ ξ ( 0 , t ) = P ξ ( 0 , t ) i = 1 Z d ( 0 ) 0 ξ i h i ( s ) P ( s , t ) d s k ( 0 , 0 ) P ( 0 , t ) .
By (51) and (52), we get
k s s ( ξ , s ) k ξ ξ ( ξ , s ) + σ k ( ξ , s ) = 0 , 2 d d ξ k ( ξ , ξ ) σ = 0 , i = 1 Z d i ( ξ ) h i ( 0 ) + k s ( ξ , 0 ) = 0 , τ i + d i ( ξ ) h i ( ξ i ) = 0 , d i ( ξ ) h i ( s ) + d i ( ξ ) h i ( s ) + σ d i ( ξ ) h i ( s ) τ i k ( ξ i , s ) j = 1 Z τ i d j ( ξ i ) h j ( s ) = 0 , h i ( ξ i ) = d ( 0 ) = 0 .
Lemma 2.
Under assumption (11), for (53), there exist classical solutions d i ( · ) C 2 ( [ 0 , 1 ] ) , h i ( · ) C 2 ( [ 0 , ξ i ] ) , and k ( · , · ) C 2 ( Λ ) .
Proof of Lemma 2.
Similarly, let (53) split into the following
τ i + d i ( ξ ) h i ( ξ i ) = 0 , d i ( ξ ) h i ( s ) + d i ( ξ ) h i ( s ) + σ d i ( ξ ) h i ( s ) τ i k ( ξ i , s ) j = 1 Z τ i d j ( ξ i ) h j ( s ) = 0 , h i ( ξ i ) = d ( 0 ) = 0 ,
and
k s s ( ξ , s ) k ξ ξ ( ξ , s ) + σ k ( ξ , s ) = 0 , 2 d d ξ k ( ξ , ξ ) σ = 0 , i = 1 Z d i ( ξ ) h i ( 0 ) + k s ( ξ , 0 ) = 0 .
Firstly, we prove that for (54), there exist classical solutions d i ( · ) C 2 ( [ 0 , 1 ] ) , h i ( · ) C 2 ( [ 0 , ξ i ] ) . By τ i + d i ( ξ ) h i ( ξ i ) = 0 in (54), one has
d i ( ξ ) = τ i h i ( ξ i ) .
Combining the first equation of (54) and (56), we have
h i ( s ) + ( σ τ ) h i ( s ) + h i ( ξ i ) k ( ξ i , s ) = 0 , h i ( ξ i ) = 0 .
We discuss (57) by the following three cases.
Case 1: τ > σ
By (57), we have
h i ( s ) = h i ( ξ i ) [ 1 2 τ σ ( e τ σ s τ σ ξ i e τ σ s + τ σ ξ i ) e τ σ s 2 τ σ ξ i s e τ σ s k ( ξ i , s ) d s + e τ σ s 2 τ σ ξ i s e τ σ s k ( ξ i , s ) d s ] .
Then, it gives
h i ( s ) d i ( ξ ) = τ i [ 1 2 τ σ ( e τ σ s τ σ ξ i e τ σ s + τ σ ξ i ) e τ σ s 2 τ σ ξ i s e τ σ s k ( ξ i , s ) d s + e τ σ s 2 τ σ ξ i s e τ σ s k ( ξ i , s ) d s ] .
Case 2: τ = σ
According to (57), we obtain
h i ( s ) = h i ( ξ i ) ( s ξ i ) h i ( ξ i ) ξ i s ξ i y k ( ξ i , x ) d x d y .
So we have
h i ( s ) d i ( ξ ) = τ i ( s ξ i ) + τ i ξ i s ξ i y k ( ξ i , x ) d x d y .
Case 3: τ < σ
Similarly, by (57), one has
h i ( s ) = h i ( ξ i ) 1 σ τ [ ( sin ( σ τ ξ i ) cos ( σ τ s ) + sin ( σ τ s ) cos ( σ τ ξ i ) ) + cos ( σ τ s ) x i s sin ( σ τ x ) k ( ξ i , x ) d x sin ( σ τ s ) x i s cos ( σ τ x ) k ( ξ i , x ) d x ] .
Then,
h i ( s ) d i ( ξ ) = τ i 1 σ τ [ ( sin ( σ τ ξ i ) cos ( σ τ s ) + sin ( σ τ s ) cos ( σ τ ξ i ) ) + cos ( σ τ s ) x i s sin ( σ τ x ) k ( ξ i , x ) d x sin ( σ τ s ) x i s cos ( σ τ x ) k ( ξ i , x ) d x ] .
Now, we prove that there exists classical solution k ( · , · ) C 2 ( Λ ) for (55). We firstly show that kernel k ( ξ , s ) is independent of h i ( s ) . In fact, we obtain h i ( x i ) = h i ( 0 ) from Case 1–Case 3. For convenience, let A = h i ( 0 ) h i ( ξ i ) . By (58), we get
h i ( 0 ) = h i ( ξ i ) [ 1 2 ( e τ σ ξ i + e τ σ ξ i ) + 1 2 0 ξ i e τ σ x k ( ξ i , x ) d x + 1 2 0 x i e τ σ x k ( ξ i , x ) d x ] .
Define
f ( ξ ) = 0 ξ i e τ σ x k ( ξ i , x ) d x , g ( ξ ) = 0 ξ i e τ σ x k ( ξ i , x ) d x .
Then, combining (55), we obtain
f ( x ) τ f ( ξ ) = α e τ σ x k s ( ξ , 0 ) τ σ k ( ξ , 0 ) , f ( 0 ) = 0 , f ( 0 ) = 0 ,
and
g ( x ) τ g ( ξ ) = α e τ σ ξ k s ( ξ , 0 ) + τ σ k ( ξ , 0 ) , g ( 0 ) = 0 , g ( 0 ) = 0 .
Let F ( ξ ) = f ( ξ ) + g ( ξ ) , one has
F ( ξ ) = 2 k s ( ξ , 0 ) τ k ( ξ , s ) e τ σ ξ i e τ σ ξ i .
Then,
A = k s ( ξ , 0 ) τ .
Hence, A is any constant in this case, moreover, k s ( ξ , 0 ) = A τ . By (60), we have
h i ( 0 ) = h i ( ξ i ) + h i ( ξ i ) 0 ξ i k ( ξ i , x ) d x .
By (66) with τ = σ , we can get 0 ξ i k ( ξ i , x ) d x = 0 . So A = 1 holds by (66). By Equation (63), there is
h i ( 0 ) = h i ( ξ i ) cos ( σ τ ξ i ) + 1 σ τ 0 ξ i cos ( σ τ x ) k ( ξ i , x ) d x .
Let
γ ( ξ ) = 0 ξ cos ( σ τ x ) k ( ξ , x ) d x .
From (72), one has
γ ( ξ ) τ γ ( ξ ) = σ cos ( σ τ ξ ) k s ( ξ , 0 ) , γ ( 0 ) = 0 , γ ( 0 ) = 0 .
There is a unique solution in (73), which is written by
γ ( ξ ) = cos ( σ τ ξ ) + k s ( ξ , 0 ) τ .
Combining (72) and (74) at ξ = ξ i , i = 1 , 2 , , Z , we obtain
γ ( ξ i ) = 0 ξ i cos ( σ τ x ) k ( ξ i , x ) d x = cos ( σ τ ξ i ) + k s ( ξ i , 0 ) τ .
Then, substituting (75) into (71), we have
A = k s ( ξ , 0 ) τ .
We can see A is also any constant in this case.
Next, we prove that A = 1 . Substituting (48) into (3), we make the kernels satisfying
a i ( ξ ) [ b i ( s ) + s ξ i b i ( x ) k ( x , s ) d x ] = h i ( s ) [ d i ( ξ ) j = 1 Z a j ( ξ ) 0 ξ i b j ( s ) d j ( s ) d s 0 ξ c ( ξ , x ) d j ( s ) d x ] ,
and
k ( ξ , s ) = c ( ξ , s ) + s ξ c ( ξ , x ) k ( x , s ) d x .
From (77), we have
h i ( 0 ) h i ( ξ i ) = 1 b i ( ξ i ) ( b i ( 0 ) + 0 ξ i b i ( x ) k s ( x , 0 ) d x ) .
For simplifying calculation, we only give the result under case τ > 0 . In fact, case τ > 0 is the same as τ < 0 . Substituting b i ( s ) = ζ i ( sinh ( λ s ) tanh ( λ ξ i ) cosh ( λ s ) ) of (28) into (79), one has A = 1 .
Finally, we can see that (55) is the same as (13) except for the last equations. However, it can be transformed into form “ H ( ξ ) ” for the last of Equation (55). So the remaining steps are similar from (43) to (47). □
Remark 2.
Similar to Lemma 1, from the above proof, we cannot explain that there are unique solutions d i ( · ) C 2 ( [ 0 , 1 ] ) , h i ( · ) C 2 ( [ 0 , ξ i ] ) for Equation (54). However, d i ( ξ ) h i ( s ) is unique, which indicates transformation (48) is unique.
We design the state feedback controller
U ^ ( t ) = i = 1 Z a i ( 1 ) 0 ξ i b i ( s ) Γ ( s , t ) d s + 0 1 c ( 1 , s ) Γ ( s , t ) d s .
The closed-loop system (2) corresponding with controller (80) is described by
Γ t ( ξ , t ) = Γ ξ ξ ( ξ , t ) + i = 1 Z τ i Γ ( ξ i , t ) , 0 < ξ < 1 , t ( 0 , + ) , Γ ξ ( 0 , t ) = 0 , t 0 , Γ ( 1 , t ) = i = 1 Z a i ( 1 ) 0 ξ i b i ( s ) Γ ( s , t ) d s + 0 1 c ( 1 , s ) Γ ( s , t ) d s , t 0 , Γ ( ξ , 0 ) = Γ 0 ( ξ ) , 0 ξ 1 .
The system (81) is considered in the state space H = L 2 ( 0 , 1 ) . We define system operator A : D ( A ) H for closed-loop system (81)
A ϕ = ϕ + i = 1 Z τ i ϕ ( ξ i , t ) , ϕ A , A = { ϕ H 2 ( 0 , 1 ) A ϕ H , ϕ ( 1 ) = i = 1 Z a i ( 1 ) 0 ξ i b i ( s ) Γ ( s , t ) d s + 0 1 c ( 1 , s ) Γ ( s , t ) d s } .
Theorem 1.
Under assumption (11), for each initial value Γ 0 ( · ) H , the closed-loop system (81) admits a unique solution Γ ( · , t ) C ( [ 0 , + ) ; H ) . Moreover, the closed-loop system generates an exponentially stable C 0 -semigroup such that
e A t L A e ρ t , t 0 ,
where L A and ρ two positive constants.
Proof. 
Define linear operator A 0 : D ( A 0 ) H for the target system (4)
A 0 φ = φ + σ φ ; φ D ( A 0 ) , D ( A 0 ) = { φ H 2 ( 0 , 1 ) φ ( 0 ) = 0 , φ ( 1 ) = 0 } .
For any ω D ( A 0 ) ,
Re A 0 φ , φ = Re φ + σ φ , φ φ 2 + σ φ 2 0 ,
which means A 0 is dissipative on H . On the other hand, for any φ ^ H , A 0 ( φ ) = φ ^ , by φ ^ = φ + σ φ , we have
φ ( ξ ) = 1 σ [ coth ( σ ) 0 1 cosh ( σ s ) φ ^ ( s ) d s 0 1 sinh ( σ s ) φ ^ ( s ) d s ] cosh ( σ ξ ) + 1 σ sinh ( σ ξ ) 0 ξ cosh ( σ s ) φ ^ ( s ) d s 1 σ cosh ( σ ξ ) 0 ξ sinh ( σ s ) φ ^ ( s ) d s ,
which shows that A 0 1 L ( H ) is compact on H . By the Lumer Phillips theorem [38], A 0 generates a compressed C 0 -semigroup e A 0 on H . That is to say,
| | e A 0 t | | L A 0 e δ t ,
where L A 0 are δ two positive constants. Define Lyapunov function
E ( t ) = 1 2 0 1 P 2 ( ξ , t ) d ξ .
Find the derivative of (85), and the calculus of the derivatives gives as the following quantity
E ˙ ( t ) = 0 1 P ξ 2 d t + 2 σ E ( t ) 2 σ E ( t ) ,
which means
E ( t ) e 2 σ t E ( 0 ) .
Thus, the target system (4) is exponentially stable. Based on the transformation (3) and reversible transformation (48), we define the following bounded invertible operator K : H H
K ϕ = φ , φ = ϕ i = 1 Z a i ( 1 ) 0 ξ i b i ( s ) Γ ( s , t ) d s 0 1 c ( 1 , s ) Γ ( s , t ) d s , K 1 φ = ϕ , ϕ = φ + i = 1 Z d i ( 1 ) 0 ξ i h i ( s ) P ( s , t ) d s + 0 1 k ( 1 , s ) P ( s , t ) d s .
So there exists a bounded reversible operator K satisfying A = K 1 A 0 K . Hence, the operator A yields an exponentially stable C 0 -semigroup on H , that is,
e A t = K 1 e A 0 t K .
The Theorem 1 holds by the exponential stability of e A 0 t . □

3. Output Feedback Controller Design

In this section, an observer-based output feedback control for system (2) is designed as follows
Γ ^ t ( ξ , t ) = Γ ^ ξ ξ ( ξ , t ) + i = 1 Z τ i Γ ( ξ i , t ) , 0 < ξ < 1 , t ( 0 , + ) , Γ ^ ( 1 , t ) = U ^ ( t ) , t 0 , Γ ^ ξ ( 0 , t ) = 0 , t 0 , Γ ^ ( ξ , 0 ) = Γ ^ 0 ( ξ ) , 0 ξ 1 .
Let Γ ¯ = Γ ^ Γ . Error system Γ ¯ is given by
Γ ¯ t ( ξ , t ) = Γ ¯ ξ ξ ( ξ , t ) , 0 < ξ < 1 , t ( 0 , + ) , Γ ¯ ( 1 , t ) = 0 , t 0 , Γ ¯ ξ ( 0 , t ) = 0 , t 0 , Γ ¯ ( ξ , 0 ) = Γ ¯ 0 .
It is well known that system (88) yields a unique exponentially stable solution in state space H . We propose the following output feedback controller
U ^ ( t ) = i = 1 Z a i ( 1 ) 0 x i b i ( s ) Γ ^ ( s , t ) d s + 0 1 c ( 1 , s ) Γ ^ ( s , t ) d s ,
where a i , b i and c are determined by (10).
The closed-loop system composed of system (2) and system (87) under observer (89), is governed by
Γ t ( ξ , t ) = Γ ξ ξ ( ξ , t ) + i = 1 Z τ i Γ ( ξ i , t ) , 0 < ξ < 1 , t ( 0 , + ) , Γ ( 1 , t ) = i = 1 Z a i ( 1 ) 0 x i b i ( s ) Γ ^ ( s , t ) d s + 0 1 c ( 1 , s ) Γ ^ ( s , t ) d s , t 0 , Γ ξ ( 0 , t ) = 0 , t 0 , Γ ^ t ( ξ , t ) = Γ ^ ξ ξ ( ξ , t ) + i = 1 Z τ i Γ ( ξ i , t ) , 0 < ξ < 1 , t ( 0 , + ) , Γ ^ ( 1 , t ) = i = 1 Z a i ( 1 ) 0 ξ i b i ( s ) Γ ^ ( s , t ) d s + 0 1 c ( 1 , s ) Γ ^ ( s , t ) d s , t 0 , Γ ^ ξ ( 0 , t ) = 0 , t 0 , Γ ( ξ , 0 ) = Γ 0 ( ξ ) , Γ ^ ( ξ , 0 ) = Γ ^ 0 ( ξ ) , 0 ξ 1 .
Theorem 2.
Under assumption (11), for any initial datum ( Γ ( · , 0 ) , Γ ^ ( · , 0 ) ) H × H , the closed-loop system (90) admits a unique solution ( Γ ( · , t ) , Γ ^ ( · , t ) ) C ( [ 0 , + ) ; H × H ) . Moreover, the closed-loop system (90) is exponentially stable, that is, for any t > 0 , there exist L A ^ > 0 , δ ^ > 0 such that
( Γ ( · , t ) , Γ ^ ( · , t ) ) H L A ^ e δ ^ t ( Γ ( · , 0 ) , Γ ^ ( · , 0 ) ) H .
Proof. 
To prove our result, we aim to transform system (90) into an equivalent one. Then, we divide the equivalent system into two sub-systems which are proved to be well-posed and exponentially stable.
Owing to
Γ ^ Γ ¯ = 0 1 1 1 Γ Γ ^ ,
the closed-loop system (90) is equivalent to PDEs as follows
Γ ^ t ( ξ , t ) = Γ ^ ξ ξ ( ξ , t ) + i = 1 Z τ i Γ ( ξ i , t ) , x ( 0 , 1 ) , t ( 0 , + ) , Γ ^ ( 1 , t ) = i = 1 Z a i ( 1 ) 0 ξ i b i ( s ) Γ ^ ( s , t ) d s + 0 1 c ( 1 , s ) Γ ^ ( s , t ) d s , t 0 , Γ ^ ξ ( 0 , t ) = 0 , t 0 , Γ ¯ t ( ξ , t ) = Γ ¯ ξ ξ ( ξ , t ) , ξ ( 0 , 1 ) , t ( 0 , + ) , Γ ¯ ( 1 , t ) = 0 , t 0 , Γ ¯ ξ ( 0 , t ) = 0 , t 0 , Γ ^ ( ξ , 0 ) = Γ ^ 0 ( ξ ) , Γ ¯ ( ξ , 0 ) = Γ ¯ 0 , 0 ξ 1 .
Next, we prove that system(93) is exponentially stable. Notice that “ Γ ¯ -part” of system (93) is special case at σ = 0 for system (4), so “ Γ ¯ -part” has a unique solution which is exponentially stable. Now, we only need to prove that “ Γ ^ -part” of system (93) has a unique solution and is exponentially stable. Similar to transformation (3), we give the following transformation
p ( ξ , t ) = Γ ^ ( ξ , t ) i = 1 Z a i ( ξ ) 0 ξ i b i ( s ) Γ ^ ( s , t ) d s 0 ξ c ( ξ , s ) Γ ^ ( s , t ) d s ,
under which “ Γ ^ -part” of system (93) is mapped to
p t ( ξ , t ) = p ξ ξ ( ξ , t ) + σ p ( ξ , t ) + i = 1 Z Γ ¯ ( ξ i , t ) [ 1 + j = 1 Z a j ( ξ ) 0 ξ j b j ( s ) d s + 0 1 c ( ξ , s ) d s ] , ξ ( 0 , 1 ) , t ( 0 , + ) , p ( 1 , t ) = 0 , t 0 , p ξ ( 0 , t ) = 0 , t 0 .
Now, we are in a position to prove that system (95) admits a unique solution and is exponentially stable. Based on (84), the solution of system (95) can be obtained as
p ( · , t ) = e A 0 t p ( · , 0 ) + 0 t e A 0 ( t s ) F p ( ξ ) i = 1 Z Γ ¯ ( ξ i , s ) d s ,
where F p ( ξ ) = 1 + j = 1 Z a j ( ξ ) 0 ξ j b j ( s ) d s + 0 1 c ( ξ , s ) d s . Define the following energy function E 1 ( t ) of system (88)
E 1 ( t ) = 1 2 0 1 ( Γ ¯ ( ξ , t ) ) 2 d ξ .
A simple computation of derivative of (97) with respect to t shows that
E ˙ 1 ( t ) = 0 1 ( Γ ¯ ξ ( ξ , t ) ) 2 d ξ 0 1 ( Γ ¯ ( ξ , t ) ) 2 d ξ = 2 E 1 ( t ) ,
which means
E 1 ( t ) e 2 t E 1 ( 0 ) .
For any t 2 > t 1 > 0 , one has
t 1 t 2 | Γ ¯ ξ ( ξ , s ) | 2 d s = t 1 t 2 | ξ 1 Γ ¯ ξ ( ξ , s ) d ξ | 2 d s t 1 t 2 | 0 1 Γ ¯ ξ ( ξ , s ) d ξ | 2 d s t 1 t 2 0 1 | Γ ¯ ξ ( ξ , s ) | 2 d ξ d s = t 1 t 2 E ˙ 1 ( s ) d s = E 1 ( t 1 ) E 1 ( t 2 ) E 1 ( 0 ) ( e 2 δ 1 t 1 + e 2 δ 1 t 2 ) .
Owing to the proof process of Theorem 1, there is
e A 0 t p ( · , 0 ) e A 0 t p ( · , 0 ) L A 0 e δ t p ( · , 0 ) .
Simultaneously,
0 t e A 0 ( t s ) F p ( ξ ) i = 1 Z Γ ¯ ( ξ i , s ) d s 0 t e A 0 ( t s ) F p ( ξ ) i = 1 Z Γ ¯ ( ξ i , s ) d s .
Additionally, from (99), for any i = 1 , 2 , , Z , we obtain(see [39])
0 t e δ 1 ( t s ) | Γ ¯ ( ξ i , s ) | d s = 0 t 2 e δ 1 ( t s ) | Γ ¯ ( ξ i , s ) | d s + t 2 t e δ 1 ( t s ) | Γ ¯ ( ξ i , s ) | d s = 0 t 2 e δ 1 ε | Γ ¯ ( ξ i , t ε ) | d ε + t 2 t e δ 1 ( t s ) | Γ ¯ ( ξ i , s ) | d s [ 0 t 2 e 2 δ 1 ε d ε ] 1 2 [ 0 t 2 | Γ ¯ ( ξ i , t ε ) | 2 d ε ] 1 2 + [ t 2 t e 2 δ 1 ( t s ) d s ] 1 2 [ t 2 t | Γ ¯ ( ξ i , s ) | 2 d s ] 1 2 1 2 δ 1 e 2 δ 1 t [ 0 | Γ ¯ ( ξ i , ε ) | 2 d ε ] + 1 2 δ 1 E 1 ( 0 ) e δ 1 t 4 .
Then, it can be obtained that
0 t e δ 1 ( t s ) | i = 1 Z Γ ¯ ( ξ i , s ) | d s Z e 2 δ 1 t 2 δ 1 [ 0 | Γ ¯ ( ξ i , ε ) | 2 d ε ] + Z 2 δ 1 e δ 1 t 4 E 1 ( 0 ) .
Therefore, system (95) is exponentially stable, which means Theorem 2 holds. □

4. Simulation Results

In this section, some simulation results are presented to explain the effectiveness of proposed controller by the finite element method. For the open-loop system of (2) and the closed-loop system (90), we choose parameters
τ 1 = 1 , τ 2 = 1 , τ 3 = 3 ; ξ 1 = 0.3 , ξ 2 = 0.4 , ξ 3 = 0.5 ; σ = 0 ,
and initial values Γ 0 = 6 sin ( π ξ ) , Γ ^ 0 = 4 cos ( π ξ ) .
Figure 1 displays the solution of open-loop system (2). We can see that the solution of open-loop system (2) is growing fast and is unstable. Figure 2a,b display the solutions of “ Γ -part” and “ Γ ^ -part” of closed-loop system (90), respectively. Trajectory of the controller (89) is displayed in Figure 3. So it is clearly seen that the solution of closed-loop system (90) decays to zero and is stable under controller (89).

5. Concluding Remarks

The paper is presented the boundary output feedback stabilization of a heat equation with multi-point heat source mainly by backstepping approach. We first design the state feedback controller based on backstepping transformation. The exponential stability of closed-loop system is guaranteed. Secondly, the observer-based output feedback controller is constructed for infinite-dimensional systems. Furthermore, we prove that closed-loop system is exponentially stable. In the future, we will consider the result without assumption (11) and the case that τ i is a variable function or a more general continuous function hold.

Author Contributions

Methodology, F.-F.J.; investigation, B.-Q.Y.; formal analysis, Q.-Q.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the NSFC of China (62073203).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Masi, A.D.; Ferrari, P.A.; Lebowitz, J.L. Reaction-diffusion equations for interacting particle systems. J. Stat. Phys. 1986, 44, 589–644. [Google Scholar] [CrossRef]
  2. Price, C.B.; Wambacq, P.; Oosterlinck, A. Applications of reaction-diffusion equations to image processing. In Proceedings of the Third International Conference on Image Processing and its Applications, Warwick, UK, 18–20 July 1989. [Google Scholar]
  3. Illés, B.; Harsányi, G. Investigating direction characteristics of the heat transfer coefficient in forced convection reflow oven. Exp. Therm. Fluid Sci. 2009, 33, 642–650. [Google Scholar] [CrossRef]
  4. Pardo, F.; López, P.; Cabello, D.; Balsi, M. Efficient software-hardware 3D heat equation solver with applications on the non-destructive evaluation of minefields. Comput. Geosci. 2009, 35, 2239–2249. [Google Scholar] [CrossRef]
  5. Gao, Y.; Shao, Q.; Yan, B.Z.; Li, Q.F.; Guo, S.X. Parabolic equation modeling of electromagnetic wave propagation over rough sea surfaces. Sensors 2019, 19, 1252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Cheng, H.M.; Yuan, R. The spreading property for a prey-predator reaction-diffusion system with fractional diffusion. Frac. Calc. Appl. Anal. 2015, 3, 565–579. [Google Scholar] [CrossRef]
  7. Guo, W.; Jin, F.F. Adaptive error feedback regulator design for one-dimensional heat equation with unknown harmonic disturbance anticollocated with control. IEEE Trans. Autom. Control 2020, 2, 824–830. [Google Scholar] [CrossRef]
  8. Jin, F.F.; Guo, B.Z. Performance boundary output tracking for one-dimensional heat equation with boundary unmatched disturbance. Automatica 2018, 96, 1–10. [Google Scholar] [CrossRef]
  9. Kumar, S.; Cao, J.D.; Li, X.D. A numerical method for time-fractional reaction-diffusion and integro reaction-diffusion equation based on quasi-wavelet. Complexity 2020, 2020, 3291723. [Google Scholar] [CrossRef]
  10. Ma, T.F.; Yan, B.Q. Positive global solutions of nonlocal boundary value problems for the nonlinear convection reaction-diffusion equations. Bound. Value Probl. 2017, 2017, 9. [Google Scholar] [CrossRef] [Green Version]
  11. Ren, Z.G.; Xu, C.; Zhou, Z.C.; Wu, Z.Z.; Chen, T.H. Boundary stabilization of a class of reaction-advection-diffusion systems via a gradient-based optimization approach. J. Frankl. Inst. 2019, 356, 173–195. [Google Scholar] [CrossRef]
  12. Liu, W. Boundary feedback stabilization of an unstable heat equation. SIAM J. Control Optim. 2003, 42, 1033–1043. [Google Scholar] [CrossRef]
  13. Wu, H.N.; Wang, J.W.; Li, H.X. Exponential stabilization for a class of nonlinear parabolic PDE systems via fuzzy control approach. IEEE Trans. Fuzzy Syst. 2012, 20, 318–329. [Google Scholar] [CrossRef]
  14. Vazquez, R.; Krstic, M. Control of 1-D parabolic PDEs with Volterra nonlinearities, Part I: Design. Automatica 2008, 44, 2778–2790. [Google Scholar] [CrossRef]
  15. Boškoviā, D.; Balogh, A.; Krstic, M. Backstepping in infinite dimension for a class of parabolic distributed parameter systems. Math. Control Signal Syst. 2003, 16, 44–75. [Google Scholar] [CrossRef]
  16. Deutsche, J.; Gabrie, J. Fredholm backstepping control of coupled linear parabolic PDEs with input and output delays. IEEE Trans. Autom. Control 2014, 65, 3128–3135. [Google Scholar] [CrossRef]
  17. Kerschbaum, S.; Deutscher, J. Backstepping control of coupled linear parabolic PDEs with space and time dependent coefficients. IEEE Trans. Autom. Control 2020, 65, 3060–3067. [Google Scholar] [CrossRef]
  18. Sano, H. Neumann boundary stabilization of one-dimensional linear parabolic systems with input delay. IEEE Trans. Autom. Control 2018, 63, 3105–3111. [Google Scholar] [CrossRef]
  19. Zhou, Z.C.; Guo, C.L. Stabilization of linear heat equation with a heat source at intermediate point by boundary control. Automatica 2013, 49, 448–456. [Google Scholar] [CrossRef]
  20. Cerpa, E.; Coron, J.M. Rapid Stabilization for a Korteweg-de Vries equation from the left Dirichlet boundary condition. IEEE Trans. Autom. Control 2013, 58, 1688–1695. [Google Scholar] [CrossRef] [Green Version]
  21. Si, Y.C.; Xie, C.K.; Zhen, Z.Y.; Zhao, A.L. Local stabilization of coupled nonlinear parabolic equations by boundary control. J. Frankl. Inst. 2018, 355, 5592–5612. [Google Scholar] [CrossRef]
  22. Yu, X.; Xu, C.; Chu, J. Local exponential stabilization of Fisher’s equation using the backstepping technique. Syst. Control Lett. 2014, 74, 1–7. [Google Scholar] [CrossRef]
  23. Jadachowski, L.; Meurer, T.; Kugi, A. Backstepping observers for periodic quasi-linear parabolic PDEs. In Proceedings of the Proceedings of the 19th World Congress the International Federation of Automatic Control, Cape Town, South Africa, 24–29 August 2014. [Google Scholar]
  24. Benabdallah, A. Stabilization of a class of nonlinear uncertain ordinary differential equation by parabolic partial differential equation controller. Int. J. Robust Nonlinear Control 2020, 30, 3023–3038. [Google Scholar] [CrossRef]
  25. Tang, S.X.; Xie, C.K. Stabilization for a coupled PDE-ODE control system. J. Frankl. Inst. 2011, 348, 2142–2155. [Google Scholar] [CrossRef]
  26. Tang, S.X.; Xie, C.K.; Zhou, Z.C. Stabilization for a class of delayed coupled PDE-ODE systems with boundary control. In Proceedings of the 2011 Chinese Control and Decision Conference, Mianyang, China, 23–25 May 2011. [Google Scholar]
  27. Wang, J.; Krstic, M. Output feedback boundary control of a heat PDE sandwiched between two ODEs. IEEE Trans. Autom. Control 2009, 64, 4653–4660. [Google Scholar] [CrossRef]
  28. Baccoli, A.; Pisanoa, A.; Orlov, Y. Boundary control of coupled reaction-diffusion processes with constant parameters. Automatica 2015, 54, 80–90. [Google Scholar] [CrossRef]
  29. Dai, J.G.; Ren, B.B. UDE-based robust boundary control for an unstable parabolic PDE with unknown input disturbance. Automatica 2018, 93, 363–368. [Google Scholar] [CrossRef]
  30. Kang, W.; Guo, B.Z. Stabilization of unstable cascaded heat PDE system subject to boundary disturbance. IET Control Theo. Appl. 2016, 10, 1027–1039. [Google Scholar] [CrossRef] [Green Version]
  31. Krstic, M.; Smyshlyaev, A. Adaptive boundary control for unstable parabolic PDEs-Part I: Lyapunov design. IEEE Trans. Autom. Control. 2008, 53, 1575–1591. [Google Scholar] [CrossRef]
  32. Ragusa, M.A. Dirichlet problem associated to divergence form parabolic equations. Commun. Contemp. Math. 2005, 7, 551. [Google Scholar] [CrossRef] [Green Version]
  33. Smyshlyaev, A.; Krstic, M. Adaptive boundary control for unstable parabolic PDEs-Part II: Estimation-based designs. Automatica 2007, 43, 1543–1556. [Google Scholar] [CrossRef]
  34. Zhan, H.; Feng, Z. Stability of polytropic filtration equation with variable exponents. Adv. Differ. Equ. 2020, 25, 255–278. [Google Scholar]
  35. Wang, L.P.; Jin, F.F. Boundary output feedback stabilization of transport equation with non-local term. IMA J. Math. Control Inf. 2020, 37, 752–764. [Google Scholar] [CrossRef]
  36. Argomedo, F.B.; Krstic, M. Backstepping-forwarding control and observation for hyperbolic PDEs with fredholm integrals. IEEE Trans. Autom. Control 2015, 60, 2145–2160. [Google Scholar] [CrossRef] [Green Version]
  37. Bernard, P.; Krstic, M. Adaptive output-feedback stabilization of non-local hyperbolic PDEs. Automatica 2014, 50, 2692–2699. [Google Scholar] [CrossRef]
  38. Pazy, A. Semigroups of Linear Operators and Applications to Partial Differential Equations; Springer: New York, NY, USA, 1983. [Google Scholar]
  39. Su, L.L.; Guo, W.; Wang, J.M.; Krstic, M. Boundary stabilization of wave equation with velocity recirculation. IEEE Trans. Autom. Control. 2017, 62, 4760–4767. [Google Scholar] [CrossRef]
Figure 1. Solution of open-loop system (2).
Figure 1. Solution of open-loop system (2).
Mathematics 09 00834 g001
Figure 2. Solution of closed-loop system (90): (a) “ Γ -part” of system (90); (b) “ Γ ^ -part” of system (90).
Figure 2. Solution of closed-loop system (90): (a) “ Γ -part” of system (90); (b) “ Γ ^ -part” of system (90).
Mathematics 09 00834 g002
Figure 3. Trajectory of the controller (89).
Figure 3. Trajectory of the controller (89).
Mathematics 09 00834 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hu, Q.-Q.; Jin, F.-F.; Yan, B.-Q. Boundary Stabilization of Heat Equation with Multi-Point Heat Source. Mathematics 2021, 9, 834. https://doi.org/10.3390/math9080834

AMA Style

Hu Q-Q, Jin F-F, Yan B-Q. Boundary Stabilization of Heat Equation with Multi-Point Heat Source. Mathematics. 2021; 9(8):834. https://doi.org/10.3390/math9080834

Chicago/Turabian Style

Hu, Qing-Qing, Feng-Fei Jin, and Bao-Qiang Yan. 2021. "Boundary Stabilization of Heat Equation with Multi-Point Heat Source" Mathematics 9, no. 8: 834. https://doi.org/10.3390/math9080834

APA Style

Hu, Q. -Q., Jin, F. -F., & Yan, B. -Q. (2021). Boundary Stabilization of Heat Equation with Multi-Point Heat Source. Mathematics, 9(8), 834. https://doi.org/10.3390/math9080834

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop