Next Article in Journal
Forecasting the Volatility of the Cryptocurrency Market by GARCH and Stochastic Volatility
Next Article in Special Issue
Optimal Location and Sizing of DGs in DC Networks Using a Hybrid Methodology Based on the PPBIL Algorithm and the VSA
Previous Article in Journal
Nanofluid Flow on a Shrinking Cylinder with Al2O3 Nanoparticles
Previous Article in Special Issue
A Control Method for IPMSM Based on Active Disturbance Rejection Control and Model Predictive Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Fault Estimation for Hyperbolic PDEs

1
School of Automatin, Central South University, Changsha 480013, China
2
Department of Chemical & Materials Engineering, University of Alberta, Edmonton, AB T6G 2V4, Canada
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(14), 1613; https://doi.org/10.3390/math9141613
Submission received: 28 June 2021 / Revised: 5 July 2021 / Accepted: 6 July 2021 / Published: 8 July 2021
(This article belongs to the Special Issue Model Predictive Control and Optimization for Cyber-Physical Systems)

Abstract

:
The new adaptive fault estimation scheme is proposed for a class of hyperbolic partial differential equations in this paper. The multiplicative actuator and sensor faults are considered. There are two cases that require special consideration: (1). only one type of fault (actuator or sensor) occurs; (2). two types of faults occurred simultaneously. To solve the problem of fault estimation, three challenges need to be solved: (1). No prior information of fault type is known; (2). Unknown faults are always coupled with state and input; (3). Only one boundary measurement is available. The original plant is converted to Observer canonical form. Two filters are proposed and novel adaptive laws are developed to estimate unknown fault parameters. With the help of the proposed update laws, the true state of the faulty plant can be estimated by the proposed observers composed of two filters. By selecting a suitable Lyapunov function, it is proved that under unknown external disturbance, the estimation errors of state parameters and fault parameters decay to arbitrarily small value. Finally, the validity of the proposed observer and adaptive laws is verified by numerical simulation.

1. Introduction

Process monitoring involving the fault detection and estimation is crucial to process safety. The existing technology can be divided into three categories: based on mathematical model, based on data-driven empirical model and fault detection technology combining empirical model and plant prior information. Data driven method can refer to system identification areas [1,2,3,4], etc. While the mathematical model method has better extrapolation, and empirical model is more convenient to design. Many works have made important contributions in the aspect of fault detection: see [5,6,7,8,9,10,11,12]. Those contributions contain adaptive observer [13,14], sliding mode observer [15]. The recent works [16,17] present full-order and reduced-order sliding mode observer (SMO) for Markov jump systems. In [18], the linear matrix inequality (LMI) theory is used to realize the fault estimation of synchronous actuator and sensor for Markov jump system. On this basis, the stochastic process model is established, the approximate fault estimation is made, and the time to failure (TTF) estimation is achieved in [19,20].
In practical application, robust fault estimation is the key to realize real-time monitoring, diagnosis and fault-tolerant control. Many contributions have been made to solve the related problems. In [21], the augmented system is first configured to contain the plant state of interest and the augmented state of the fault. An optimization method of unknown input observer (UIO) based on LMI is proposed to minimize the stability of the estimation error dynamic and the influence of disturbance. Aiming at the problem of state estimation and fault estimation for discrete time systems, in [22], a robust estimation method based on linear matrix inequality (LMI) technique and generalized systems theory is proposed. An observer design method is proposed in [23] for Takagi-Sugeno fuzzy systems in the presence of process uncertainties and unexpected failures. Concretely, a robust state estimation and fault estimation for all controlled objects are realized by using UIO technology, augmented system theory and sliding mode control method.
A lot of industrial systems, such as fluid flow and chemical reaction processes, all exhibit spatio-temporal dynamics, that is, the state of the system is spatial-varying and time-varying. Hence, the temporal dynamic representation is not suitable for displaying their behavior, and the observer design is a challenging problem. Moreover, the fault detection and estimation of space-time systems are more challenging. Most of the existing fault detection methods are based on modal analysis. For example, [24] uses the modal analysis technique based method to develop a fault detection observer. Similarly, the work [25] designs a robust observer to solve the problem of fault detection and state estimation.
The above techniques make use of the model approximation, which results in false alarm due to model reduction. In [26], an adaptive PDE observer is developed for positive system to realize fault accommodation. In our recent work [27], two filters are proposed based on the Luenberger observer, and then novel parameter adaptation laws were proposed to address boundary fault estimation issues.
In this paper, we propose a scheme for state and fault estimation of a linear hyperbolic PDE system with unknown disturbances based on the plant observer canonical form. Unlike [28], where additive faults are considered and some real positive conditions need to be satisfied, this article considers multiplicative faults. Especially, when the actuator and sensor fail at the same time, the estimation problem becomes more challenging.
The contributions are as follows:
  • The problem of state and fault estimation in two cases is solved: (1). only one type of fault (sensor or actuator) occurs; (2). the sensor and actuator failed simultaneously.
  • By proposing novel transformation, the observer canonical form of the system of interest is obtained. Consequently, two auxiliary filters are designed, and based on which the novel adaptation laws are proposed.
  • We show that in sense of Lyapunov stability fault parameter estimation converges exponentially fast.
This paper is organized as follows: Section 2 introduces the system considered and multiplicative boundary faults. Section 3 describes the design of fault parameter adaptation law based on the observer canonical form. In Section 4, numerical simulations are performed to illustrate the validity of the proposed results.

Notation

Given the vector variables v ( ξ ) , ξ [ 0 , 1 ] which is spatially varying and continuous, the following operator is defined:
H β v ( ξ ) = 0 1 e β ξ v ( ξ ) d ξ ,
where β is a real number, with the derived norm
v β 2 = H β v T v = 0 1 e β ξ v T ( ξ ) v ( ξ ) d ξ .
The integral operator (1) has the following property
2 H β v v ξ = e β v 2 ( 1 ) v 2 ( 0 ) β v β 2 ,
where v ξ represents the derivative of v ( ξ ) w.r.t the spatial variable ξ . In subsequent sections, we will write t v ( ξ , t ) (or v t ( ξ , t ) ) and ξ v ( ξ , t ) (or v ξ ( ξ , t ) ) to stand for v ( ξ , t ) t and v ( ξ , t ) ξ , respectively.
The norm v β is equivalent to the standard L 2 -norm, namely, there exist positive constants κ 1 , κ 2 such that
κ 1 v β v 0 κ 2 v β .
For brevity, the argument in time and space are often omitted, i.e., u = u ( t ) , v = v ( ζ , t ) .

2. Problem Statement

We consider the following PDEs:
v t ( ξ , t ) = v ξ ( ξ , t ) + g ( ξ ) v ( 0 , t ) + ϕ ( t ) + 0 ξ f ( ξ , ζ ) v ( ζ , t ) d ζ + ξ 1 h ( ξ , ζ ) v ( ζ , t ) d ζ ,
v ( 1 , t ) = u ( t ) ,
y ( t ) = v ( 0 , t )
defined on the domain t R + , ξ 0 , 1 , where g , f and h denote system coefficient functions, v ( ξ , t ) L 2 ( 0 , 1 ) , ( ξ , t ) [ 0 , 1 ] × R + stands for the system state, L 2 ( 0 , 1 ) denotes a real Hilbert space with the inner product v 1 , v 2 = 0 1 v 1 ( z ) v 2 ( z ) d z , v 1 , v 2 L 2 ( 0 , 1 ) . As shown in Figure 1, and the sensor output is denoted by y ( t ) , and u ( t ) stands for actuator input. ϕ ( t ) is an external unknown disturbance satisfying:
Assumption 1.
The disturbance ϕ ( t ) is bounded:
| ϕ ( t ) | ϕ ¯ , t R +
where ϕ ¯ is a known constant.
From Figure 1, the boundary condition and the output can be defined according to fault occurrence (at t = T 0 ):
  • Actuator Fault:
    v ( 1 , t ) = u ( t ) , t < T o , θ a u ( t ) , t T o ,
  • Sensor Fault:
    y ( t ) = v ( 0 , t ) , t < T o , θ s v ( 0 , t ) , t T o ,
    where θ a and θ s denote actuator and sensor faults, respectively, θ s R + and θ a R .
In practice, due to the obstacle of unknown faults, the information of boundary conditions v ( 1 , t ) and v ( 0 , t ) is not available.
In this work, the objective is to estimate the unknown fault during the following fault occurrences:
  • [p1]. Actuator fault occurrence (9) after time t = T o ;
  • [p2]. Sensor fault occurrence (10) after time t = T o ;
  • [p3]. Simultaneous sensor and actuator fault occurrence at time t = T o .
Assumption 2.
Fault parameters are assumed to be time-invariant and bounded:
θ a θ ¯ a m
θ ̲ s l θ s θ ¯ s m ,
with θ ¯ a m > 0 , θ ̲ s l > 0 , θ ¯ s m > 0 as known parameters.

3. Adaptive Estimation Laws

Since the available information is the sensor output y ( t ) and the input u ( t ) , it is challenging to address the fault estimation problems. To this end, different strategies will be proposed based on the cases [p1], [p2] and [p3].
We define the triangles
T L = z , ζ 0 , 1 × 0 , 1 , z ζ
T U = z , ζ 0 , 1 × 0 , 1 , z ζ
and the spaces V L = C T L ; R and V U = C T U ; R .
The considered real process system (called as plant) (5)–(7) is first converted into the following form:
t v o ( ξ , t ) = ξ v o ( ξ , t ) + Φ ( ξ ) v o ( 0 , t ) + Ψ ( ξ ) v o ( 1 , t ) + T ϕ t , v o ( 1 , t ) = θ a u ( t ) , y ( t ) = θ s v o ( 0 , t )
by performing the following transformation and inverse transformation:
v ( ξ , t ) = Q p o , q o + I L 2 v o ( · , t ) ξ
and
v o ( z , t ) = Q k o , l o + I L 2 v ( · , t ) ( ξ )
with k o , p o V L and l o , q o V U , I L 2 denoting the identity operator on L 2 ( 0 , 1 ) and the operator Q ψ 1 , ψ 2 defined by:
Q ψ 1 , ψ 2 m ( ξ ) = 0 ξ ψ 1 ( ξ , ζ ) m ( ζ ) d ζ + ξ 1 ψ 2 ( ξ , ζ ) m ( ζ ) d ζ
for all m L 2 ( 0 , 1 ) , and ξ [ 0 , 1 ] given ψ 1 V L and ψ 2 V U .
In (13), T ϕ t = I L 2 + Q k o , l o ϕ t ( 1 ) , and Φ ( ξ ) and Ψ ( ξ ) give the flexibility to make transformation successful. Here, variables k o , l o satisfy the following set of kernel equations:
ξ k o ( ξ , ζ ) + ζ k o ( ξ , ζ ) = f ( ξ , ζ ) Φ ( ξ ) l ( 0 , ζ ) + ξ 1 l o ( ξ , s ) f ( s , ζ ) d s + ζ ξ k o ( ξ , s ) f ( s , ζ ) d s + 0 ζ k o ( ξ , s ) h ( s , ζ ) d s ξ , ζ [ 0 , 1 ]   s . t . ζ ξ , ξ 1 ,
ξ l o ( ξ , ζ ) + ζ l o ( ξ , ζ ) = h ( ξ , ζ ) Φ ( ξ ) l o ( 0 , ζ ) + ζ 1 l o ( ξ , s ) f ( s , ζ ) d s + ξ ζ l o ( ξ , s ) h ( s , ζ ) d s + 0 ξ k o ( ξ , s ) h ( s , ζ ) d s ξ , ζ [ 0 , 1 ]   s . t . ξ ζ
with boundary conditions:
l o ( 0 , ζ ) = 0 , k o ( 1 , ζ ) = 0 .
From (15) and (19), it follows that
v ( 0 , t ) = v o ( 0 , t ) , v ( 1 , t ) = v o ( 1 , t ) .
During the transformation, the terms Φ ( ζ ) and Ψ ( ζ ) can be computed:
Ψ ( ζ ) = l o ( ξ , 1 ) ,
Φ ( ξ ) = I L 2 + Q k o , l o [ g ( · ) ] ( ξ ) k o ( ξ , 0 ) .
Remark 1.
Transformations in (14) and (15) is bounded, and let us define k ¯ o = max ξ , ζ [ 0 , 1 ] k o ( ξ , ζ ) and l ¯ o = max ξ , ζ [ 0 , 1 ] l o ( ξ , ζ ) , then we can obtain T ϕ ( t ) T ¯ ϕ , where T ¯ ϕ = 1 + k ¯ o + l ¯ o ϕ ¯ .

3.1. Filters and Non-Adaptive Estimate

We design the following filters:
Λ t ( ξ , t ) = Λ ξ ( ξ , t ) , Λ ( ξ , 0 ) = Λ 0 ( ξ ) , ξ 0 , 1 ,
Λ ( 1 , t ) = y ( t )
and
Γ t ( ξ , t ) = Γ ξ ( ξ , t ) , Γ ( ξ , 0 ) = Γ 0 ( ξ ) , ξ 0 , 1 ,
Γ ( 1 , t ) = u ( t ) ,
where Λ 0 and Γ 0 are the initial conditions of the proposed filers. The explicit solutions to the filters (23)–(24) and (25)–(26) can be obtained:
Λ ( ξ , t ) = y ( t + ξ 1 ) ,
Γ ( ξ , t ) = u ( t + ξ 1 ) ,
Corresponding to different cases of fault occurrence, we consider the following non-adaptive state estimates:
c1. 
(Actuator fault estimation) For the problem [p1]:
v ¯ o ( ξ , t ) = θ a Γ ( ξ , t ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ , t d ζ + ξ 1 Φ ( ζ ) Λ 1 ζ ξ , t d ζ .
c2. 
(Sensor fault estimation) For the problem [p2]:
v ¯ o ( ξ , t ) = Γ ( ξ , t ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ , t d ζ + 1 θ s ξ 1 Φ ( ζ ) Λ 1 ζ ξ , t d ζ .
c3. 
(Simultaneous fault estimation) For the problem [p3]:
v ¯ o ( ξ , t ) = θ a Γ ( ξ , t ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ , t d ζ + 1 θ s ξ 1 Φ ( ζ ) Λ 1 ζ ξ , t d ζ .
Lemma 1.
We consider the plant (5)–(7) with different fault occurrence, the filters (23)–(26) and the state estimate v ¯ o in (29), (30) and (31). After a finite time T = 1 , with ϕ 0 , one will have
v ¯ o v o .
Proof. 
Let us define e μ = v o v ¯ o , then we have:
t e μ ( ξ ) = ξ e μ ( ξ ) , e μ ( 1 ) = 0 ,
the solution of which is
e μ ( ξ , t ) = e μ ( 1 ξ t , 0 ) , t 1 ξ 0 , t > 1 ξ ,
which directly shows that after T = 1 the e μ 0 , and therefore this completes the proof. □

3.2. Adaptive Laws

Motivated by the above results, the following adaptive observers are constructed:
c1. 
(Healthy State estimation) Healthy case:
v ^ o ( ξ , t ) = Γ ( ξ , t ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ , t d ζ + ξ 1 Φ ( ζ ) Λ 1 ζ ξ , t d ζ ,
y ^ ( t ) = v ^ o ( 0 , t ) ,
e y ( t ) = y ( t ) y ^ ( t )
where, e y ( t ) is the detection residual. v ^ o ( ξ , t ) and y ^ ( t ) denote estimates of v o ( ξ , t ) and y ( t ) , respectively.
c2. 
(Actuator fault estimation) For the problem [p1]:
v ^ o ( ξ , t ) = θ ^ a ( t ) Γ ( ξ , t ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ , t d ζ + ξ 1 Φ ( ζ ) Λ 1 ζ ξ , t d ζ ,
e y ( t ) = y ( t ) y ^ ( t ) = v o ( 0 , t ) v ^ o ( 0 , t )
where θ ^ a ( t ) is the estimate of θ a .
c3. 
(Sensor fault estimation) For the problem [p2]:
v ^ o ( ξ , t ) = Γ ( ξ , t ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ , t d ζ + 1 θ ^ s ( t ) ξ 1 Φ ( ζ ) Λ 1 ζ ξ , t , d ζ
e y ( t ) = y ( t ) y ^ ( t ) = θ s v o ( 0 , t ) θ ^ s ( t ) v ^ o ( 0 , t ) ,
where θ ^ s ( t ) denotes the estimate of θ s .
c4. 
(Simultaneous fault estimation) For the problem [p3]:
v ^ o ( ξ , t ) = θ ^ a ( t ) Γ ( ξ , t ) + z 1 Ψ ( ζ ) Γ 1 ζ ξ , t d ζ + 1 θ ^ s ( t ) ξ 1 Φ ( ζ ) Λ 1 ζ ξ , t d ζ ,
e y ( t ) = y ( t ) y ^ ( t ) = θ s v o ( 0 , t ) θ ^ s ( t ) v ^ o ( 0 , t ) .

3.3. Single Type of Fault Estimation

In this subsection, the adaptive law and the adaptive observer will be proposed to address the problem of the single type fault estimation, i.e, [p1] and [p2] respectively. In particular, the output residual e y ( t ) can be calculated: for actuator fault estimation [p1] e y ( t ) = v o ( 0 , t ) v ^ o ( 0 , t ) ; for sensor fault estimation [p2] e y ( t ) = θ s v o ( 0 , t ) θ ^ s ( t ) v ^ o ( 0 , t ) .
It is given that only one type of fault happens. In the problem [p1], the following adaptation law is employed:
θ ^ ˙ a ( t ) = proj [ θ ¯ a m , θ ¯ a m ] β G Γ ( t ) ( 0 ) e y ( t ) γ θ ^ a ( t ) 1 + d 2 , θ ^ a ( t ) , θ ^ a ( 0 ) [ θ ¯ a m , θ ¯ a m ]
with design parameters β > 0 and γ > 0 , where the operator G is defined by
G Γ ( t ) ( ξ ) = Γ ( ξ , t ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ , t d ζ
to estimate v o ( ξ , t ) and y ( t ) . In the problem [p2], the following adaptation law is applied:
θ ^ ˙ s ( t ) = proj [ θ ̲ s l , θ ¯ s m ] β G Γ ( t ) ( 0 ) e y ( t ) γ θ ^ s ( t ) 1 + d 2 , θ ^ s ( t ) , θ ^ s ( 0 ) [ θ ̲ s l , θ ¯ s m ]
to estimate the plant state, with
d 2 = μ 2 ( 0 ) + θ s m 2 μ 2 ( 0 ) + G Γ ( 0 ) 2
and the projection operator given as:
p r o j [ a , b ] η , ω = 0 if ω = a and η 0 , 0 if ω = b and η 0 , η otherwise .
Theorem 1.
The adaptive laws (44) and (46) corresponding to actuator and sensor faults, respectively, provide the following properties:
θ ^ a ( t ) θ ¯ a m , θ ̲ s l θ ^ s ( t ) θ ¯ s m , t 0 ,
θ ˜ a , θ ˜ s L ,
θ ^ ˙ a , θ ^ ˙ s L L 2 ,
ρ L L 2 ,
where θ ˜ a ( t ) = θ a θ ^ a ( t ) and θ ˜ s ( t ) = θ s θ ^ s ( t ) and
ρ = e y 1 + d 2
with d 2 defined in (47).
Proof. 
We define a new error signal for actuator and sensor faults, respectively:
μ ( ξ ) = v o ( ξ ) ξ 1 Φ ( ζ ) Λ 1 ζ ξ d ζ θ a Γ ( ξ ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ d ζ
and
μ ( ξ ) = v o ( ξ ) ξ 1 Φ ( ζ ) Λ 1 ζ ξ d ζ θ s Γ ( ξ ) ξ 1 Ψ ( ζ ) Γ 1 ζ ξ d ζ .
Therefore,
μ t ( ξ ) = μ ξ ( ξ ) + T ( ϕ ) , μ ( 1 ) = 0 .
(a) 
Proof part of Actuator fault:
Construct a Lyapunov function candidate:
V 1 = 1 ε μ b 2 + 1 2 η θ ˜ a 2 ,
where the design parameters b, ε and η are positive constants. Taking time derivative of (57) gives:
V ˙ 1 = 2 ε H b μ ξ μ + 2 ε H b μ T ( ϕ ) + θ ˜ a θ ˜ ˙ a η ,
where | T ( ϕ ) | T ¯ ϕ . Applying the property (3) and Young’s inequality, plugging the parameter adaptation law (44), and using θ ˜ proj [ θ ̲ , θ ¯ ] τ , ρ ^ θ ˜ τ yield
V ˙ 1 2 ε e b μ 2 ( 1 ) μ 2 ( 0 ) b μ b 2 θ ˜ a θ ^ ˙ a η + 2 ε H b μ T ¯ ϕ 2 ε μ 2 ( 0 ) 2 ε b μ b 2 + 1 ε μ b 2 + 1 ε T ¯ ϕ 2 β η θ ˜ a G [ Γ ( t ) ] ( 0 ) e y ( t ) 1 + d 2 + γ η θ ˜ a θ ^ a 1 + d 2
therein the identity θ ˜ ˙ a = θ ^ ˙ a is employed. Moreover, the inequality
2 ε H b μ T ¯ ϕ 1 ε μ b 2 + 1 ε T ¯ ϕ 2
was used. By combining (38)–(39) and (54), for the actuator fault, one gets:
e y = v o ( 0 ) v ^ o ( 0 ) = μ ( 0 ) + θ ˜ a G Γ ( 0 ) .
Based on the equation (60), one can replace θ ˜ a G Γ ( t ) ( 0 ) in (59), inserting θ ^ a = θ a θ ˜ a and using Cauchy-Schwartz inequality give:
V ˙ 1 2 ε μ 2 ( 0 ) 1 ε 2 b 1 μ b 2 β η 1 1 + d 2 e y 2 + β 2 η 1 1 + d 2 μ 2 ( 0 ) + β 2 η 1 1 + d 2 e y 2 γ 2 η 1 1 + d 2 θ ˜ a 2 + γ 2 η 1 1 + d 2 θ a 2 + 1 ε T ¯ ϕ 2 = 2 ε 1 1 + d 2 β 2 η μ 2 ( 0 ) 1 ε 2 b 1 μ b 2 1 1 + d 2 γ 2 η θ ˜ a 2 1 1 + d 2 β 2 η e y 2 ( t ) + 1 1 + d 2 γ 2 η θ a 2 + 1 ε T ¯ ϕ 2 2 c 1 V 1 + Δ 1
with the positive constant: c 1 = 1 2 min 2 b 1 , γ and
Δ 1 = 1 1 + d 2 γ 2 η θ ¯ a m 2 + 1 ε T ¯ ϕ 2 ,
where we choose the design parameters satisfying
2 ε > 1 1 + d 2 β 2 η , 2 b > 1 .
From (62), by choosing large values of η and ε , a small value of γ , then Δ 1 becomes an arbitrarily small positive constant. Based on the inequality (61), when V 1 ( t 0 ) > Δ 1 c 1 and V 1 ( t ) > Δ 1 c 1 for all t > t 0 , then V ˙ 1 ( t ) c 1 V 1 ( t ) , which means V 1 ( t ) V 1 ( t 0 ) e c 1 t . Hence, V 1 ( t ) is exponentially decreasing. Then, after a finite time t f = 1 c 1 ln c 1 V 1 ( t 0 ) Δ 1 , we get V 1 ( t ) Δ 1 c 1 , for all t > t f . Integrating (61) from zero to infinity gives ρ L 2 . From (60), we get for t > 1
e y 1 + d 2 = μ ( 0 ) 1 + d 2 + θ ˜ a G Γ ( 0 ) 1 + d 2 1 + θ ˜ a G Γ ( 0 ) 1 + d 2
which shows ρ L . From the adaptive law (44), we get
θ ^ ˙ a β e y 1 + d 2 G Γ ( 0 ) 1 + d 2 β ρ .
(b) 
Proof part of Sensor fault:
Construct the following Lyapunov candidate
V 2 = 1 ε μ b 2 + θ ˜ s 2 2 η θ ¯ s m 2
with parameters b > 0 , ε > 0 and η > 0 . Differentiating V 2 w.r.t time along the dynamics (56), one obtains:
V ˙ 2 = 2 ε H b μ ξ μ + 2 ε H b μ T ( ϕ ) + θ ˜ s θ ˜ ˙ s η θ ¯ s m 2 .
Applying the property (3) and Young’s inequality, substituting the parameter adaptation law (44), and applying θ ˜ proj [ θ ̲ , θ ¯ ] τ , ρ ^ θ ˜ τ , and inserting the parameter adaptation law (46), one gets:
V ˙ 2 2 ε e b μ 2 ( 1 ) 2 ε μ 2 ( 0 ) 2 ε b μ b 2 + 2 ε H b μ T ¯ ϕ θ ˜ s θ ^ ˙ s η θ ¯ s m 2 2 ε μ 2 ( 0 ) 1 ε 2 b 1 μ b 2 β θ ˜ s G Γ ( t ) ( 0 ) e y ( t ) η θ ¯ s m 2 1 1 + d 2 + γ θ ˜ s θ ^ s η θ ¯ s m 2 1 1 + d 2 + 1 ε T ¯ ϕ 2 ,
where θ ˜ ˙ s = θ ^ ˙ s is used. From (40)–(41) and (55), one obtains:
e y = θ s v o ( 0 ) θ ^ s v ^ o ( 0 ) = θ s μ ( 0 ) + θ ˜ s G Γ ( 0 ) .
Using (69) to replace the nonlinear term θ ˜ s ( t ) G Γ ( t ) ( 0 ) and employing Cauchy-Schwartz inequality give:
V ˙ 2 2 ε μ 2 ( 0 ) 1 ε 2 b 1 μ b 2 + 1 ε T ¯ ϕ 2 1 1 + d 2 β e y 2 ( t ) η θ ¯ s m 2 1 1 + d 2 β θ s μ ( 0 ) e y ( t ) η θ ¯ s m 2 + 1 1 + d 2 γ θ ˜ s θ s θ ˜ s η θ ¯ s m 2 2 c 2 V 2 + Δ 2
with the positive constant c 2 = 1 2 min γ , 2 b 1 and
Δ 2 = 1 1 + d 2 γ 2 η + 1 ε T ¯ ϕ 2 ,
where we choose the design parameters
2 ε > 1 1 + d 2 β 2 η , 2 b > 1 .
From the inequality (70) and for V 2 ( t 0 ) > K 2 Δ c 2 and V 2 ( t ) > K 2 Δ c 2 , t > t 0 , we have V ˙ 2 ( t ) V 2 ( t 0 ) e c 2 t , which implies that V 2 ( t ) is decreasing exponentially and there exists a finite time t f = 1 c 2 ln c 2 V 2 ( t 0 ) K 2 Δ such that V 2 ( t ) K 2 Δ c 2 , t > t f .
Integrating (70) from zero to infinity gives ρ L 2 . From (69), we get for t > 1
e y 1 + d 2 = θ ¯ s m μ ( 0 ) 1 + d 2 + θ ˜ s G Γ ( 0 ) 1 + d 2 1 + θ ˜ s G Γ ( 0 ) 1 + d 2
which shows ρ L . From the adaptive law (46), we get
θ ^ ˙ s β e y 1 + d 2 G Γ ( 0 ) 1 + d 2 β ρ .
This completes the proof. □

3.4. Simultaneous Faults Estimation

We turn to the most challenging estimation problem [p3]. The problem [p3] includes cases [p1] and [p2]: [p1] is equivalent to [p3] with θ s = 1 while [p2] can be regarded as [p3] with θ a = 1 . To solve the problem [p3], it is required to develop coupled parameter adaptation laws for estimating θ a and θ s simultaneously.
Let us consider that actuator and sensor faults simultaneously occur to the plant (5)–(7). We construct the adaptive observer (42) and (43) and the following parameter adaptation laws:
θ ^ ˙ a ( t ) = proj [ θ ¯ a m , θ ¯ a m ] β G Γ ( t ) ( 0 ) e y ( t ) γ a θ ^ a ( t ) 1 + d 2 , θ ^ a ( t ) , θ ^ a ( 0 ) [ θ ¯ a m , θ ¯ a m ] ,
θ ^ ˙ s ( t ) = proj [ θ ̲ s l , θ ¯ s m ] β θ ^ a ( t ) G Γ ( t ) ( 0 ) e y ( t ) γ s θ ^ s ( t ) 1 + d 2 , θ ^ s ( t ) , θ ^ s ( 0 ) [ θ ̲ s l , θ ¯ s m ]
with design parameters β > 0 , γ a > 0 and γ s > 0 .
Theorem 2.
(Simultaneous fault estimation) The adaptive laws (75) and (76) provide the following properties:
θ ^ a ( t ) θ ¯ a m , θ ̲ s l θ ^ s ( t ) θ ¯ s m , t 0 ,
θ ˜ a , θ ˜ s L ,
θ ^ ˙ a , θ ^ ˙ s L L 2 ,
ρ L L 2 .
Proof. 
By following the similar steps in Theorem 1, we introduce the following variable:
μ ( ξ ) = v o ( ξ ) ξ 1 Φ ( ζ ) Λ 1 ζ ξ d ζ θ s θ a Γ ( ξ ) + ξ 1 Ψ ( ζ ) Γ 1 ζ ξ d ζ .
Based the equations (42), (43) and (81), the error signal v ˜ ( ξ ) can be expressed:
v ˜ o ( ξ ) = v o ( ξ ) v ^ o ( ξ ) = μ ( ξ ) + θ ˜ s θ s θ ^ s ξ 1 Φ ( ζ ) Λ 1 ζ ξ d ζ + θ ˜ a G [ Γ ] ( ξ ) .
A direct calculation shows that μ also satisfies
μ t ( ξ ) = μ ξ ( ξ ) + T ( ϕ ) , μ ( 1 ) = 0 .
In addition, the detection residual e y ( t ) is expressed as:
e y = θ s v o ( 0 ) θ ^ s v ^ o ( 0 ) = θ s μ 0 + θ s θ ˜ a + θ ˜ s θ ^ a G [ Γ ] ( 0 ) .
We consider a Lyapunov candidate:
V 3 = 1 ε μ b 2 + θ s θ ˜ a 2 2 η Θ ¯ s m 2 + θ ˜ s 2 2 η Θ ¯ s m 2
with parameters b > 0 , η > 0 and ε > 0 , where Θ ¯ s m = max θ ¯ s m , θ ¯ a m . The sensor fault parameter θ s is positive and hence V 3 is positive definite. Taking the time derivative of V 3 gives:
V ˙ 3 = 2 ε H b μ ξ μ + 2 ε H b μ T ( ϕ ) + θ ˜ s θ ˜ ˙ s η Θ ¯ s m 2 + θ s θ ˜ a θ ˜ ˙ a η Θ ¯ s m 2 .
We use property (3) and Young’s inequality, the adaptation law (75) and (76), and θ ˜ proj [ θ ̲ , θ ¯ ] τ , ρ ^ θ ˜ τ to obtain:
V ˙ 3 2 ε μ 2 ( 0 ) 2 b 1 ε μ b 2 + 1 1 + d 2 θ s γ a θ ˜ a θ ^ a η Θ ¯ s m 2 + 1 1 + d 2 γ s θ ˜ s θ ^ s η Θ ¯ s m 2 1 1 + d 2 β θ s θ ˜ a + θ ˜ s θ ^ a G [ Γ ( t ) ] ( 0 ) e y ( t ) Θ ¯ s m 2 + 1 ε T ¯ ϕ 2 .
According to (84), one can replace θ s θ ˜ a ( t ) + θ ˜ s ( t ) θ ^ a ( t ) G [ Γ ( t ) ] ( 0 ) and apply Cauchy-Schwartz inequality to obtain:
V ˙ 3 2 ε μ 2 ( 0 ) 2 b 1 ε μ b 2 1 1 + d 2 θ s γ a θ ˜ a 2 η Θ ¯ s m 2 1 1 + d 2 γ s θ ˜ s 2 η Θ ¯ s m 2 + 1 1 + d 2 γ s θ ˜ s θ s η Θ ¯ s m 2 + 1 1 + d 2 θ s γ a θ ˜ a θ a η Θ ¯ s m 2 1 1 + d 2 β e y 2 ( t ) θ s μ ( 0 ) e y ( t ) Θ ¯ s m 2 + 1 ε T ¯ ϕ 2 2 c 3 V 3 ( t ) + Δ 3
with c 3 = 1 2 min 2 b 1 , γ a , γ s and Δ 3 = θ ¯ s m γ a 2 η 1 1 + d 2 + γ s 2 η 1 1 + d 2 + 1 ε T ¯ ϕ 2 , where θ ¯ s m / Θ ¯ s m 1 and θ ¯ a m / Θ ¯ s m 1 are used. The design parameters can be chosen:
2 ε > 1 1 + d 2 β 2 η , 2 b > 1 .
Repeating the similar steps as Theorem 1, we see that after a finite time t f , V 3 ( t ) Δ 3 c 3 holds. Hence, θ ˜ s ( t ) and θ ˜ a ( t ) become arbitrarily small. □

4. Simulation Results

In this section, we will examine the performance of the proposed method by simulating the occurrence of different faults [p1], [p2] and [p3].
We use the system in [29] as a representative system, which is a Korteweg-de Vries (KdV)-like equation typically describing the dynamics of shallow water wave. The functions in the system are given by
g ( ξ ) = q + p l o cosh ( o ) sinh ( o ( 1 ξ ) ) ,
f ( ξ , ζ ) = p l cosh ( o z ) cosh ( o ( 1 ζ ) ) cosh ( o ) + p l cosh ( o ( ξ ζ ) ) ,
h ( ξ , ζ ) = p l cosh ( o ξ ) cosh ( o ( 1 ζ ) ) cosh ( o )
with q = 1.25 , p = o = 0.1 and l = 10 . We formulate the boundary control:
u ( t ) = U s + 0 1 p s ( 1 , ζ ) v ( ζ ) d ζ
with U s = 50 , and the second term is a stabilization term. p s ( z , ζ ) satisfies the kernel Equation (14) in [30]. Moreover, the initial condition is v ( ξ , 0 ) = 10 sin ( π ξ ) for z [ 0 , 1 ] .
The actuator and sensor faults are as follows:
θ a = 3.5 , t > 8 , 1 + 2.5 ( 1 e 2 ( t 5 ) ) , t [ 5 , 8 ] , 1 , t 5 ,
θ s = 2 , t > 6 , 1 + ( 1 e 3 ( t 5 ) ) , t [ 5 , 6 ] , 1 , t 5 .
Different occurrences with different exogenous disturbance ϕ ( t ) are shown in Table 1, where ‘★’ indicates the fault occurrences and ‘-’ shows no faults.

4.1. Adaptive Fault Estimation

According to [p1], [p2] and [p3], we formulated the adaptive observers (38)–(43), respectively. Moreover, initial conditions are θ ^ a ( 0 ) = 1 and θ ^ s ( 0 ) = 1 . The design parameters are listed in the following table:

4.1.1. Actuator Fault Estimation

We simulate the system in the case of [p1]. The actuator fault occurs at T o = 5 s and becomes constant after T = 8 s . In Figure 2, the oscillation arises due to the perturbation ϕ ( t ) = ϕ 1 as shown in Table 1. Applying the adaptation law (44) ensures that in Figure 2 the estimate y ^ ( t ) ( 0 , t ) = v ^ (blue dotted line) converges to the measurement y ( t ) (solid line). As shown in Figure 3, the estimate θ ^ a ( t ) is exponentially convergent. We also use the common Lunberger observer (35)–(37), but the estimate y ^ x ( t ) (black dotted line) deviates from the measurement. Particularly, the mean square error (MSE) is calculated for the output estimates to show performance.

4.1.2. Sensor Fault Estimation

In this subsection, we investigate the performance of the proposed observer (40) and (41) and the adaptive law (46). We simulate the situation when the same system fails. Sensor fault occurs at T o = 5 s , and is kept unchanging after t = 6 s . A new external perturbation ϕ = ϕ 2 shown in Table 2 effects the system, where ‘★’ indicates the fault occurrences and ‘-’ shows no faults. A comparison between y ^ x ( t ) and y ^ ( t ) = θ ^ s ( t ) v ^ ( 0 , t ) , see Figure 4 and Figure 5, where y ^ x ( t ) is generated by the common Lunberger observer (35)–(37) and y ^ ( t ) is from the recommended observer (40) and (41). Apparently, the y ^ ( t ) (blue dashed line) exponentially converges to the measurement and the resulting residual e y ( t ) exponentially decays to a very small value in Figure 5. However, the common Lunberger observer (35)–(37) fails to give satisfactory results.

4.1.3. Simultaneous Fault Estimation

In this section, we carry out a simulation with the fault [P3]: both actuator and sensor faults occur at T o = 5 s and remain unchanged after T = 8 s . ϕ ( t ) = ϕ 3 is the external disturbance defined in Table 1. We apply the adaptive observer (42) and (43) and coupling parameter adaptive law (75) and (76). Expected state estimates v ^ ( ξ ) , θ ^ a and θ ^ s converge to v ( ξ ) , θ a and θ s . The results in Figure 6 and Figure 7 directly show the performance of the proposed method.The observation error generated by the proposed observer (42) and (43) v ˜ ( 0 , t ) = v ( z , t ) v ^ ( z , t ) and the residual e y ( t ) attenuates to a very small value around zero: see the dashed blue line in the Figure 6. However, the common observer (35)–(37) is unable to provide acceptable estimates as shown in the black dotted line in the Figure 6.

5. Conclusions

We propose a strategy to address the estimation problem of a class of linear hyperbolic PDE systems with single and simultaneous faults and unknown disturbance. The proposed observer implements state estimation and fault parameter estimation using the boundary measurement. The detailed observer design results are presented, and the unknown faults are estimated with parameter adaptation laws. Through numerical simulations, it is shown that the proposed parameter adaptation laws ensure arbitrarily small errors in the fault parameter estimation, despite the presence of unknown external disturbances. However, the common Lunberger observer fails to provide the exact state estimates when facing the faulty system. This method has the advantages of easy design and strong portability. Once the observer canonical form can be obtained, it is possible to extend the results of this work to more complex PDE systems. In this paper, we obtain that the error of fault parameters and state estimation is ultimately bounded.

Author Contributions

Software, Y.Y.; Writing—original draft, Y.Y.; Methodology, Y.Y.; Writing—review & editing, X.X. and S.D.; Supervision, S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Central South University grant number 202045010.

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. Yang, X.; Yin, S.; Kaynak, O. Robust identification of LPV time-delay system with randomly missing measurements. IEEE Trans. Syst. Man, Cybern. Syst. 2018, 48, 2198–2208. [Google Scholar] [CrossRef]
  2. Yang, X.; Yang, X. Local identification of LPV dual-rate system with random measurement delays. IEEE Trans. Ind. Electron. 2017, 65, 1499–1507. [Google Scholar] [CrossRef]
  3. Du, Y.; Budman, H.; Duever, T.A.; Du, D. Fault Detection and Classification for Nonlinear Chemical Processes using Lasso and Gaussian Process. Ind. Eng. Chem. Res. 2018, 57, 8962–8977. [Google Scholar] [CrossRef]
  4. Du, Y.; Du, D. Fault detection and diagnosis using empirical mode decomposition based principal component analysis. Comput. Chem. Eng. 2018, 115, 1–21. [Google Scholar] [CrossRef]
  5. Yin, S.; Gao, H.; Qiu, J.; Kaynak, O. Descriptor reduced-order sliding mode observers design for switched systems with sensor and actuator faults. Automatica 2017, 76, 282–292. [Google Scholar] [CrossRef]
  6. Chen, J.; Patton, R.J. Robust Model-Based Fault Diagnosis for Dynamic Systems. Springer Science & Business Media: New York, NY, USA, 2012; Volume 3. [Google Scholar]
  7. Ding, S.X. Model-Based Fault Diagnosis Techniques: Design Schemes, Algorithms, and Tools. Springer Science & Business Media: New York, NY, USA, 2008. [Google Scholar]
  8. Wang, Y.W.; Xiao, J.W.; Wang, H.O. Global synchronization of complex dynamical networks with network failures. Int. J. Robust Nonlinear Control. 2010, 20, 1667–1677. [Google Scholar] [CrossRef]
  9. Jiang, B.; Staroswiecki, M.; Cocquempot, V. Fault accommodation for nonlinear dynamic systems. IEEE Trans. Autom. Control. 2006, 51, 1578–1583. [Google Scholar] [CrossRef]
  10. Ragusa, M.A. Hölder Regularity Results for Solutions of Parabolic Equations, Variational Analysis and Applications. Book Ser. Nonconvex Optim. Appl. 2005, 79, 921–934. [Google Scholar]
  11. Wang, Y.; Wang, J.; Zhang, H. Reliable Control for Fuzzy Stochastic Hyperbolic Systems with Time Delay. In Proceedings of the 2nd IEEE-ASME International Conference on Mechatronics and Embedded Systems and Applications, Beijing, China, 13–16 August 2006; pp. 1–5. [Google Scholar]
  12. Xu, X.; Dubljevic, S. Actuator Fault Detection and Estimation for a Class of Hyperbolic PDEs Using Filter-Based Observer. In Proceedings of the 23rd International Conference on Methods & Models in Automation & Robotics (MMAR), Miedzyzdroje, Poland, 27–30 August 2018; pp. 148–153. [Google Scholar]
  13. Cristofaro, A.; Johansen, T.A. Fault tolerant control allocation using unknown input observers. Automatica 2014, 50, 1891–1897. [Google Scholar] [CrossRef]
  14. Yang, H.; Jiang, B.; Staroswiecki, M. Observer-based fault-tolerant control for a class of switched nonlinear systems. IET Control. Theory Appl. 2007, 1, 1523–1532. [Google Scholar] [CrossRef]
  15. Liu, M.; Shi, P. Sensor fault estimation and tolerant control for Itô stochastic systems with a descriptor sliding mode approach. Automatica 2013, 49, 1242–1250. [Google Scholar] [CrossRef]
  16. Yang, H.; Yin, S. Descriptor Observers Design for Markov Jump Systems with Simultaneous Sensor and Actuator Faults. IEEE Trans. Autom. Control. 2019, 64, 370–3377. [Google Scholar] [CrossRef]
  17. Yang, H.; Yin, S. Reduced order sliding mode observer-based fault estimation for Markov jump systems. IEEE Trans. Autom. Control. 2019, 64, 473–4740. [Google Scholar] [CrossRef]
  18. Li, X.; Ahn, C.K.; Lu, D.; Guo, S. Robust simultaneous fault estimation and nonfragile output feedback fault-tolerant control for Markovian jump systems. IEEE Trans. Syst. Man Cybern. Syst. 2019, 49, 1769–1776. [Google Scholar] [CrossRef]
  19. Wang, P.; Vachtsevanos, G. Fault prognostics using dynamic wavelet neural networks. AI EDAM 2001, 15, 349–365. [Google Scholar] [CrossRef]
  20. Chinnam, R.B.; Baruah, P. Autonomous diagnostics and prognostics through competitive learning driven HMM-based clustering. Proc. Int. Jt. Conf. Neural Netw. 2003, 4, 2466–2471. [Google Scholar]
  21. Gao, Z.; Liu, X.; Chen, M.Z. Unknown input observer-based robust fault estimation for systems corrupted by partially decoupled disturbances. IEEE Trans. Ind. Electron. 2015, 63, 2537–2547. [Google Scholar] [CrossRef] [Green Version]
  22. Gao, Z. Fault estimation and fault-tolerant control for discrete-time dynamic systems. IEEE Trans. Ind. Electron. 2015, 62, 3874–3884. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, X.; Gao, Z.; Zhang, A. Observer-based fault estimation and tolerant control for stochastic Takagi–Sugeno fuzzy systems with Brownian parameter perturbations. Automatica 2019, 102, 137–149. [Google Scholar] [CrossRef]
  24. Armaou, A.; Demetriou, M.A. Robust detection and accommodation of incipient component and actuator faults in nonlinear distributed processes. AIChE J. 2008, 54, 2651–2662. [Google Scholar] [CrossRef]
  25. Ghantasala, S.; El-Farra, N.H. Robust actuator fault isolation and management in constrained uncertain parabolic PDE systems. Automatica 2009, 45, 2368–2373. [Google Scholar] [CrossRef]
  26. Demetriou, M.A.; Ito, K.; Smith, R.C. Adaptive monitoring and accommodation of nonlinear actuator faults in positive real infinite dimensional systems. IEEE Trans. Autom. Control. 2007, 52, 2332–2338. [Google Scholar] [CrossRef]
  27. Xu, X.; Yuan, Y.; Wang, M.; Fu, M.; Dubljevic, S. Fault detection and estimation for a class of PIDE systems based on boundary observers. Int. J. Robust Nonlinear Control. 2019, 29, 5867–5885. [Google Scholar] [CrossRef]
  28. Zhu, J.W.; Yang, G.H.; Wang, H.; Wang, F. Fault Estimation for a Class of Nonlinear Systems Based on Intermediate Estimator. IEEE Trans. Autom. Control. 2016, 61, 2518–2524. [Google Scholar] [CrossRef]
  29. Bribiesca-Argomedo, F.; 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]
  30. Xu, X.; Dubljevic, S. Output regulation for a class of linear boundary controlled first-order hyperbolic PIDE systems. Automatica 2017, 85, 43–52. [Google Scholar] [CrossRef]
Figure 1. Actuator and sensor fault occurrences.
Figure 1. Actuator and sensor fault occurrences.
Mathematics 09 01613 g001
Figure 2. Comparison: the measurement y ( t ) , y ^ x ( t ) generated from (35)–(37)with MSE = 24 , and y ^ ( t ) from (38) and (39) with MSE = 0.14 .
Figure 2. Comparison: the measurement y ( t ) , y ^ x ( t ) generated from (35)–(37)with MSE = 24 , and y ^ ( t ) from (38) and (39) with MSE = 0.14 .
Mathematics 09 01613 g002
Figure 3. (Upper): The residual e x y ( t ) from (35)–(37), and residual e y ( t ) from (38) and (39); (Lower): θ a and the estimate θ ^ a ( t ) by (44).
Figure 3. (Upper): The residual e x y ( t ) from (35)–(37), and residual e y ( t ) from (38) and (39); (Lower): θ a and the estimate θ ^ a ( t ) by (44).
Mathematics 09 01613 g003
Figure 4. The y ( t ) (solid line), y ^ ( t ) (blue dash line) from (40) and (41) with MSE = 0.93 , and y ^ x ( t ) (black dash line) from (35)–(37) with MSE 58.3 .
Figure 4. The y ( t ) (solid line), y ^ ( t ) (blue dash line) from (40) and (41) with MSE = 0.93 , and y ^ x ( t ) (black dash line) from (35)–(37) with MSE 58.3 .
Mathematics 09 01613 g004
Figure 5. (Upper): The residual e x y ( t ) by (35)–(37), and residual e y ( t ) by (40) and (41); (Lower): the fault θ s and estimate θ ^ s ( t ) obtained by (46).
Figure 5. (Upper): The residual e x y ( t ) by (35)–(37), and residual e y ( t ) by (40) and (41); (Lower): the fault θ s and estimate θ ^ s ( t ) obtained by (46).
Mathematics 09 01613 g005
Figure 6. The residual e x y ( t ) using (35)–(37) with MSE = 143 , and e y ( t ) from (42) and (43) with MSE = 2.5 .
Figure 6. The residual e x y ( t ) using (35)–(37) with MSE = 143 , and e y ( t ) from (42) and (43) with MSE = 2.5 .
Mathematics 09 01613 g006
Figure 7. The faults θ a and θ s (solid line), and estimates θ ^ a ( t ) by (75) and θ ^ s ( t ) by (76).
Figure 7. The faults θ a and θ s (solid line), and estimates θ ^ a ( t ) by (75) and θ ^ s ( t ) by (76).
Mathematics 09 01613 g007
Table 1. Cases [p1], [p2], [p3] with disturbances ϕ ( t ) .
Table 1. Cases [p1], [p2], [p3] with disturbances ϕ ( t ) .
Fault θ a θ s ϕ ( t )
[p1]- ϕ 1 = 0.3 cos ( 2 π t ) + sin ( π t )
[p2]- ϕ 2 = cos ( π t ) + 0.2 sin ( 1.5 π t )
[p3] ϕ 3 = cos ( π t )
Table 2. Fault occurrences, proposed observers, adaptive laws (AL) and design parameters.
Table 2. Fault occurrences, proposed observers, adaptive laws (AL) and design parameters.
Fault Occurrence θ a θ s ObserverAL β γ
[p1]-(38)–(39)(44)910.03
[p2]-(40)–(41)(46)120.05
γ a γ s
[p3](42)–(43)(75)–(76)510.020.07
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yuan, Y.; Xu, X.; Dubljevic, S. Adaptive Fault Estimation for Hyperbolic PDEs. Mathematics 2021, 9, 1613. https://doi.org/10.3390/math9141613

AMA Style

Yuan Y, Xu X, Dubljevic S. Adaptive Fault Estimation for Hyperbolic PDEs. Mathematics. 2021; 9(14):1613. https://doi.org/10.3390/math9141613

Chicago/Turabian Style

Yuan, Yuan, Xiaodong Xu, and Stevan Dubljevic. 2021. "Adaptive Fault Estimation for Hyperbolic PDEs" Mathematics 9, no. 14: 1613. https://doi.org/10.3390/math9141613

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