Next Article in Journal
Novel Authentication Protocols Based on Quadratic Diophantine Equations
Previous Article in Journal
Distributed Optimization Algorithm for Composite Optimization Problems with Non-Smooth Function
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mobile Sensor Networks for Finite-Time Distributed H Consensus Filtering of 3D Nonlinear Distributed Parameter Systems with Randomly Occurring Sensor Saturation

1
School of Internet of Things, Wuxi Vocational College of Science and Technology, Wuxi 214028, China
2
School of Internet of Things Engineering, Jiangnan University, Wuxi 214122, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(17), 3134; https://doi.org/10.3390/math10173134
Submission received: 2 August 2022 / Revised: 22 August 2022 / Accepted: 29 August 2022 / Published: 1 September 2022
(This article belongs to the Topic Distributed Optimization for Control)

Abstract

:
This paper is concerned with designing a distributed bounded H consensus filter to estimate an array of three-dimensional (3D) nonlinear distributed parameter systems subject to bounded perturbation. An optimization framework based on mobile sensing is proposed to improve the performance of distributed filters. The measurement output is obtained from a mobile sensor network, where a phenomenon of randomly occurring sensor saturation is taken into account to reflect the reality in a mobile networked environment. A sufficient condition is established by utilizing operator-dependent Lyapunov functional for the filtering error system to be finite-time bounded. Note that the velocity law of each mobile sensor is included in this condition. The effect from the exogenous perturbation to the estimation accuracy is guaranteed at a given level by means of H consensus performance constraint. Finally, simulation examples are presented to demonstrate the applicability of the theoretical results.

1. Introduction

Many physical phenomena in the real world, such as chemical reactions [1], heat conduction [2], fluid flow [3], and other spatially distributed processes, have states that depend not only on time, but also on space. In general, such phenomena can be described by partial differential equations, called distribution parameter systems (DPS). Robotic arms and aerial refueling tubes with flexible structures, temperature distribution of catalytic reaction rods in chemical processes, large-scale population migration, information dissemination, and population consumption in social economy can all be modeled as distributed parameter systems [4]. Since the 1960s, the development of modern partial differential equations and functional analysis has established a strict theoretical basis for distributed parameter systems. The well-posed issue as well as stability, controllability [5], observability, and optimal control of distributed parameter systems have been intensively studied. On the other hand, the control issues of distributed parameter systems have been an active area of research [6,7].
The filtering or state estimation issues play a critical role to signal processing and control engineering areas [8]. For distributed parameter systems, there are essentially two technical routes to address the challenge of filtering. One is to discretize the distributed parameter system into a finite dimensional system using a differential approach or finite element method. Then, finite dimensional filters are designed to achieve the state estimation of the system. A robust Kalman filter was given after discretizing the gas pipeline transient flow equation in [9]. Moreover, a discrete-time Kalman filter [10] was designed for linear infinite dimensional systems that are discretized in time rather than spatially approximated under the structure and energy preserving Crank-Nicolson framework. In [11], robust filtering of discrete processes was studied by designing recursive filters in a finite horizon with the use of two-dimensional Riccati-like difference equations. The advantage of this technical route is that various filter design schemes in finite dimensional systems can be introduced to address the filtering issue of distributed parameter systems. However, the approximation from discretization may lead to bias in the filtering results. On the flip side, an advanced research topic is the design of filters directly for distributed parameter systems. The functional analysis LMI(LOI) approach [12,13] becomes an effective tool to directly solve the filtering issue of distributed parameter systems. An extended Kalman filter for semilinear infinite dimensional systems was proposed in [14] by utilizing functional analysis. Robust H filtering is required when the system has uncertain and unknown perturbations [15]. A network-based H filter [16] was designed for a semilinear N-D diffusion equation under distributed in space measurements. A reliable H filtering issue [17] was developed for switched parabolic systems. These issues regarding H filtering have been resolved simply and effectively by adopting the LMI approach. Most of the existing filtering strategies are based on a centralized way, where information can be collected from all sensor nodes. A natural choice to enhance effectiveness and save energy is to introduce a distributed way. Each sensing node communicates only with its neighboring nodes in sensor networks. Many existing results deal with the distributed filtering issue for networked control systems such as nonlinear systems [18] and stochastic systems [19]. However, less research has been conducted on distributed filtering of distributed parameter systems. Some studies have shown that distributed estimation is quite effective in reaching consensus. In recent years, special attention has been paid to the distributed H consensus filtering of various network systems, such as linear systems [20], cyber-physical systems [21], and delay systems [22].
In the practical controller design, the system control should be finished in a given amount of time. Peter Dorato [23] introduced the short-time stability concept for this reason in 1961. Due to its importance in the study of transient performance, the subject of finite time stability has received a lot of attention recently. Finite-time stability differs from the usual Lyapunov stability in that it is concerned with the evolution of the system state over a given finite time [0, T]. In [24], the distributed H filtering in a finite-time horizon was studied for a class of Takagi–Sugeno fuzzy systems. For singular systems with Markovian jumping, finite-time H filtering was also investigated [25]. Noticeably, distributed consensus H filtering of distributed parameter systems is still lacking today. Moreover, as mentioned before, most studies on filtering of distributed parameter systems have used fixed static sensor networks. Therefore, this paper aims to develop a mobile sensing approach to finite-time convergence analysis of distributed H consensus filters for a class of 3D parabolic systems.
It is well known that the nonlinearity caused by environmental circumstances is a widely occurring factor in engineering. As of now, the processing of nonlinear terms in the system has been concerned with allowing nonlinear functions to satisfy specific conditions, including the Lipschitz condition [26], the Lipschitz-like condition based on one-norm [27], and the sector bounded constraints [28]. Moreover, it is also enabled to use the decomposition technique of nonlinearity, which divides the nonlinear function into a linear part and a high-order nonlinear part, as proposed in [29]. In complex spatio-temporal distribution processes, the two parts of this nonlinearity may appear in a probabilistic way, namely, their type or intensity may have changed randomly. This model description has been little studied, and this paper will introduce this new view to discuss the nonlinearity.
For accurate estimation of the state of the actual network, measurement outputs need to be collected and processed to minimize the effects of possible noise and incomplete information. However, technical limitations restrict the sensor from outputting a signal with infinite amplitude, which makes the output subject to “sensor saturation”. Saturation occurs not only to degrade the estimated performance, but may also lead to undesired oscillations or even unstable behavior. In actual network environments, probabilistic sensor saturation is possible due to intermittent saturation brought on by intermittent sensor failures, changeable saturation level brought on by aging sensors, and abrupt changes in the network environment. This phenomenon is referred to as randomly occurring sensor saturation. Due to the nonlinear characteristics of saturation, the randomly occurring sensor saturation is actually a scenario where the linear and higher order nonlinear parts of the nonlinearity occur in a probabilistic way.
Guiding by the arguments above, we focus on building a novel distributed H consensus filter for a class of 3D nonlinear distributed parameter systems such that the related filtering error system is finite-time stable while taking into account the noise attenuation level γ . The main contributions of this paper can be summarized as follows:
(1) A distinguishing feature of the research issue addressed is that the distributed filtering issue of 3D distributed parameter systems is studied by using an effective mobile sensing approach.
(2) The physical plant is under considered, nonlinear, bounded perturbation, randomly occurring sensor saturation, which makes more practical significance of our current study.
(3) The nonlinear function used to describe the uncertainty of a complex environment is decomposed into a part that satisfies the sector condition and a part that has the saturation property, and these two parts appear in a random way. To the authors’ knowledge, this nonlinear decomposition model is seldom seen in the processing of nonlinear terms in system filtering and control problems.
(4) The proposed finite-time distributed H consensus filtering technique not only achieves state estimation in a short time, but also the estimation error satisfies the bounded H consensus performance, which serves as a fast converging robust filtering technique suitable for complex spatially distributed processes.
(5) The system and filtering techniques are discussed in some special scenarios so that the associated filtering issues can be solved more easily.
The rest of the manuscript is structured as follows. In Section 2, the 3D nonlinear distributed parameter systems with a network of m mobile sensors is presented; some preliminaries are given. In Section 3, a mobile sensing approach is used to solve the distributed bounded H consensus filter design issue. The discussion of some scenarios related to the main results is available in Section 4. A numerical example is provided in Section 5 to show the effectiveness of the proposed approach. In Section 6, conclusions are drawn.
Notation 1.
The notation is standard, except as otherwise specified. ζ ˙ ( t ) means the derivative of the function ζ ( t ) with respect to time t. Prob { x } denotes the probability that event x will occur. E { x } stands for the expectation of stochastic variable x. I denotes the identity matrix with appropriate dimension; diag { x 1 , x 2 , , x m } stands for an m-dimensional diagonal matrix. If M is a matrix, M T denotes its transpose and λ max ( M ) means the largest eigenvalue of M. Moreover, matrices without explicit dimensions are assumed to be the appropriate dimension ones.
Let H be a Hilbert space with inner product · , · and corresponding induced norm | · | . Let B be a reflexive Banach space with norm denoted by · . It is assumed that B is embedded densely and continuously in H . Let B * denotes the conjugate dual of B with induced norm · * . It follows B H B * with both embedding dense and continuously, the result of which is that we have | g | ϱ g , g B , for any constant ϱ > 0 [30].

2. Problem Formulation and Preliminaries

In this paper, a class of 3D nonlinear distributed parameter systems is studied. Let η ( t , ζ ) denote the state of 3D distributed process at t [ 0 , + ) and at spatial location ζ = ( ζ 1 , ζ 2 , ζ 3 ) Ω , w h e r e   Ω = { ζ s | 0 ζ s l s , s = 1 , 2 , 3 } is a bounded region with smooth boundary Ω , and the measure of Ω is μ ( Ω ) > 0 . The adopted target object can be described by the following 3D nonlinear spatio-temporal distributed process.
η ( t , ζ ) t = s = 1 3 ζ s a s ( ζ s ) η ( t , ζ ) ζ s + H f ( t , η ( t , ζ ) ) + d ( ζ ) ν ( t ) z ( t ) = Ω b ( ζ ) η ( t , ζ ) d Ω
subject to the initial condition
η ( 0 , ζ 1 , ζ 2 , ζ 3 ) = η 0 ( ζ 1 , ζ 2 , ζ 3 ) ,
and having Dirichlet boundary conditions
η ( t , 0 ) = 0 , η ( t , · ) | Ω = 0 , t 0 .
The diffusion operators a s ( ζ s ) a ¯ s > 0 , s = 1 , 2 , 3 . d ( ζ ) denotes the spatial distribution of perturbation; ν ( t ) is external perturbation; b ( ζ ) is the spatial distribution of the output; H is a negative operator coefficient of nonlinear function f ( t , η ( t , ζ ) ) . It is assumed that f can be decomposed into a part f L ( t , η ( t , ζ ) ) that satisfies the sector condition and a part r ( t , η ( t , ζ ) ) that has the saturation property, and these two parts may occur in a probabilistic way.
f ( t , η ( t , ζ ) ) = ρ ( t ) f L ( t , η ( t , ζ ) ) + ( 1 ρ ( t ) ) r ( t , η ( t , ζ ) ) ,
where the stochastic variable ρ ( t ) R is a Bernoulli distributed white sequence taking values of 1 and 0 with
Prob { ρ ( t ) = 1 } = ρ ¯ Prob { ρ ( t ) = 0 } = 1 ρ ¯
where ρ ¯ [ 0 , 1 ] are known positive constant. It is assumed that the stochastic variable ρ ( t ) is independent mutually. We have
E { ρ ( t ) } = ρ ¯ and E { ( ρ ( t ) ρ ¯ ) 2 } = ρ ¯ ( 1 ρ ¯ ) .
Compared to the Lipschitz condition or sector condition, the decomposition method is another efficient method for dealing with nonlinearities. Also consider the occurrence of the decomposition term in a probabilistic way, the model is more suitable for the effect of environmental nonlinearity in practice.
The measurement output of m mobile sensors in 3D space are presented as follows.
y ( t ) = y 1 ( t ) y 2 ( t ) y m ( t ) = Ω ψ 1 ( ζ ; ζ 1 ( t ) ) η ( t , ζ ) d Ω Ω ψ 2 ( ζ ; ζ 2 ( t ) ) η ( t , ζ ) d Ω Ω ψ m ( ζ ; ζ m ( t ) ) η ( t , ζ ) d Ω ,
or ith measurement output can be stated as
y i ( t ) = Ω ψ i ( ζ ; ζ i ( t ) ) η ( t , ζ ) d Ω , i = 1 , 2 , , m ,
where ψ i ( ζ ; ζ i ( t ) ) denotes the spatial distribution of ith moving sensing device. The nonnegative function ψ i ( ζ ; ζ i ( t ) ) is bounded. The spatial distribution function of each mobile sensor, in the general sense, is given by
ψ i ( ζ s ; ζ i s ( t ) ) = φ i ( ζ ) if ζ Ω i 0 otherwise ,
where ζ 1 Ω i 1 = [ ζ i 1 δ i 1 , ζ i 1 + δ i 1 ] , ζ 2 Ω i 2 = [ ζ i 2 δ i 2 , ζ i 2 + δ i 2 ] , ζ 3 Ω i 3 = [ ζ i 3 δ i 3 , ζ i 3 + δ i 3 ] .
It is worth noting that the spatial distribution here implies that each mobile sensing device can have a different spatial distribution function. ( ζ i 1 ( t ) , ζ i 2 ( t ) , ζ i 3 ( t ) ) [ 0 , l 1 ] × [ 0 , l 2 ] × [ 0 , l 3 ] denotes the location of the time-varying center of mass of the ith sensing device. Thus, it draws the moving trajectory of the ith sensing device.
The saturation function σ ( · ) : R r R r is defined as
σ ( x ) = [ σ 1 T ( x 1 ) , σ 2 T ( x 2 ) , , σ r T ( x r ) ] T
with σ i ( x i ) = sign ( x i ) min { V i , max , V i } where V i , max is the ith element of the vector, and V max is the saturation level.
For every i, the stochastic variable ϵ i ( t ) R , i = 1 , 2 , , m is a Bernoulli distributed white sequence taking values of 1 and 0 with
Prob { ϵ i ( t ) = 1 } = ϵ ¯ i Prob { ϵ i ( t ) = 0 } = 1 ϵ ¯ i
where ϵ ¯ i [ 0 , 1 ] , i = 1 , 2 , , m are known positive constants. It is assumed that the stochastic variables ϵ i ( t ) are independent mutually in all i ( 1 i m ) . Accordingly, we have
E { ϵ i ( t ) ϵ ¯ i } = 0 and E { ( ϵ i ( t ) ϵ ¯ i ) 2 } = ϵ ¯ i ( 1 ϵ ¯ i ) .
The measurement output after considering the randomly occurring sensor saturation can be described as
y ˜ i ( t ) = ϵ i ( t ) y i ( t ) + ( 1 ϵ i ( t ) ) σ ( y i ( t ) ) , i = 1 , 2 , , m ,
In order to rewrite the 3D nonlinear distributed parameter system (1) as an evolution equation, the following notations are given, and their meanings are explained.
Consider a linear operator A : B B * satisfying the following assumptions:
Assumption 1.
A is bounded, that is
| g , A h | ε 0 g h ,
for g , h B and constant ε 0 > 0 .
Assumption 2.
A is coercive, that is
g , A g ϑ 0 g 2 ,
for g B and constant ϑ 0 > 0 .
The operator A is referred to as state operator. Moreover, the perturbation operator D : R B * is provided by D ν ( t ) = d ( ζ ) ν ( t ) , which satisfies g , D g μ d g , g . For a given time T f , 0 T f ν T ( t ) ν ( t ) d t d ¯ 2 , where μ d , d ¯ are positive constants. In addition, B : R B * is given by B g , h = Ω b ( ζ ) g ( ζ ) h ( ζ ) d Ω , which satisfies g , B g μ b g , g .
Then, it can be rewritten to obtain the corresponding compact form of the 3D nonlinear distributed parameter system (1).
η ˙ ( t ) = A η ( t ) + H F ( η ( t ) ) + D ν ( t ) z ( t ) = B η ( t )
where the state space is H = L 2 ( Ω ) , where η ( t , · ) = { η ( t , ζ ) : ζ Ω } is the state of the system. The space B is identified by the Sobolev space B = H 0 1 ( 0 , l ) = { h H 1 ( Ω ) | h | Ω = 0 } and its conjugate dual space B * is H 1 ( Ω ) . Let infinitesimal operator A = s = 1 3 ζ s a s ( ζ s ) η ( t , ζ ) ζ s and its domain be given by D ( A ) = { h L 2 ( Ω ) : h , h are absolutely continuous, h L 2 ( Ω ) and h | Ω = 0 } . The infinitesimal operator A generates a strongly continuous semigroup T ( t ) , t 0 , and the domain D ( A ) of the operator A is dense in H , [23].
The output can be expressed similarly as
y ˜ ( t ) = ϵ ( t ) y ( t ) + ( I ϵ ( t ) ) σ ( y ( t ) ) ,
where y ( t ) = Ψ ( ζ ( t ) ) η ( t ) , ϵ ( t ) = diag { ϵ 1 ( t ) , ϵ 2 ( t ) , , ϵ m ( t ) } and output operator Ψ ( ζ ( t ) ) is given by
Ψ ( ζ ( t ) ) g , h = Ω ψ 1 ( ζ ; ζ 1 ( t ) ) g ( ζ ) h ( ζ ) d Ω Ω ψ 2 ( ζ ; ζ 2 ( t ) ) g ( ζ ) h ( ζ ) d Ω Ω ψ m ( ζ ; ζ m ( t ) ) g ( ζ ) h ( ζ ) d Ω ,
where Ψ ( ζ ( t ) ) : B R 3 × R 3 × × R 3 m satisfying g , Ψ ( ζ ( t ) ) g μ c g , g .
To analyze finite-time H filtering, we consider the following assumptions and definitions.
Assumption 3.
There exists a constant l 0 such that
( s 1 s 2 ) T ( f L ( t , s 1 ) f L ( t , s 2 ) ) l 0 ( s 1 s 2 ) T ( s 1 s 2 )
Assumption 4.
r ( t , η ( t , ζ ) ) is bounded function satisfying
| r ( t , η ( t , ζ ) ) | r 0 ,
where r 0 is positive constant.
Assumption 5.
The saturation function σ ( x ) is sector bounded, i.e.,
σ T ( x ) σ ( x ) β
where β is a positive constant.
Definition 1
(Finite-time stability). For a given time constant T f , 3D filtering error systems with ν ( t ) = 0 are said to be finite-time stable with respect to ( α 1 , α 2 , T f , Q ) if
x ( t 0 ) , Q x ( t 0 ) α 1 2 x ( t ) , Q x ( t ) α 2 2
where α 2 > α 1 > 0 , Q is a positive definite operator.
Definition 2
(Finite-time boundedness). For a given time constant T f , 3D filtering error systems are said to be finite-time bounded with respect to ( α 1 2 , α 2 2 , T f , d ¯ , Q ) if condition (18) holds, where α 2 > α 1 > 0 , Q is a positive definite operator.
Definition 3
(Distributed bounded H consensus performance). The filters (20) are said to be distributed bounded H consensus filters if there exist w 0 > 0 such that the filtering error z ˜ i ( t ) satisfy the following inequalities
1 m i = 1 m | z ˜ i ( t ) | 2 γ 2 | ν ( t ) | 2 + w 0 α 2 2
where γ > 0 are some given disturbance attenuation level, for any i { 1 , 2 , , m } . If w 0 = 0 , then they are called distributed H consensus filters.
Lemma 1.
Let x , y be any m-dimensional real vectors, and let q be a positive scalar. Then the following inequality holds:
2 x , y q x , x + q 1 y , y .

3. Main Results

Consider the distributed consensus filter configuration with m mobile sensors as shown in Figure 1, where each sensor can receive information in a distributed way. Being aware of such a fact, we deployed the following filter structure on sensor node i.
η ^ ˙ i ( t ) = A η ^ i ( t ) + H F ( η ^ i ( t ) ) + Ψ * ( ζ i ( t ) ) k i [ y ˜ i ( t ) ϵ ¯ i Ψ ( ζ i ( t ) ) η ^ i ( t ) ] + G i κ j = 1 m π i j η ^ j ( t ) z ^ i ( t ) = B η ^ i ( t )
where η ^ i ( t ) is the state estimate of ith mobile sensor, and z ^ i ( t ) is the estimate for z ( t ) from the filter on ith mobile sensor; k i > 0 are the observer gains, and G i are the consensus filter gains. Moreover, κ denotes consensus strength coefficient; η ^ i ( 0 ) = η ^ i 0 η ( 0 ) for all i = 1 , 2 , , m . Π = ( π i j ) m × m is irreducible, π i j = π j i 0 , for i j and j = 1 m π i j = 0 , for all i = 1 , 2 , , m . It is obvious that Π has an eigenvalue of zero, and all other eigenvalues are negative.
Letting e i ( t ) = η ( t ) η ^ i ( t ) and z ˜ i ( t ) = z ( t ) z ^ i ( t ) , the filtering error system can be presented from (13) and (20) in the following.
e ˙ i ( t ) = A e i ( t ) + H ( F ( η ( t ) ) F ( η ^ i ( t ) ) ) Ψ * ( ζ i ( t ) ) k i [ ϵ i ( t ) y i ( t ) + ( 1 ϵ i ( t ) ) σ ( y i ( t ) ) ϵ ¯ i Ψ ( ζ i ( t ) ) η ^ i ( t ) ] G i κ j = 1 m π i j η ^ j ( t ) + D ν ( t ) z ˜ i ( t ) = B e i ( t )
Then, one can obtain
e ˙ i ( t ) = A ψ ( ζ i ( t ) ) e i ( t ) + H ( F ( η ( t ) ) F ( η ^ i ( t ) ) ) Ψ * ( ζ i ( t ) ) k i ( ϵ i ( t ) ϵ ¯ i ) Ψ ( ζ i ( t ) ) η ( t ) Ψ * ( ζ i ( t ) ) k i ( 1 ϵ i ( t ) ) σ ( y i ( t ) ) G i κ j = 1 m π i j η ^ j ( t ) + D ν ( t ) z ˜ i ( t ) = B e i ( t )
where A ψ ( ζ i ( t ) ) = A k i ϵ ¯ i Ψ * ( ζ i ( t ) ) Ψ ( ζ i ( t ) ) .

3.1. Finite-Time Bounded Analysis

Theorem 1.
Under the distributed consensus filter (20) as given, the zero solution of filtering error system (21) is finite-time bounded with respect to ( α 1 2 , α 2 2 , T f , d ¯ , A ψ ( ζ i ( t ) ) , if there exist positive scalar λ 0 , nonnegative constant ϖ and two sets of positive constants p i , q i ( i = 1 , 2 , , m ) such that the following inequalities hold:
ρ ¯ l 0 + κ λ max ( Π ) < 0 ,
α 1 2 + d ^ λ 0 α 2 2 e ϖ T f
and the velocity law of mobile sensors as follows:
ζ ˙ i s ( t ) = c i s k i ϵ ¯ i Ξ i s , s = 1 , 2 , 3 .
where d ^ = μ d 2 d ¯ 2 2 i = 1 m q i ,
Ξ i s = ζ i s δ i s ζ i s + δ i s d φ i ( ζ s ) d ζ s φ i ( ζ s ) e i 2 ( t , ζ ) d ζ s + φ i 2 ( ζ i s δ i s + 0 ) e i 2 ( t , ζ i s δ i s ) φ i 2 ( ζ i s + δ i s + 0 ) e i 2 ( t , ζ i s + δ i s ) , s = 1 , 2 , 3 .
with c i s > 0 , s = 1 , 2 , 3 , i = 1 , 2 , , m is velocity gain of each mobile sensor. The guidance strategy for mobile sensors enhances the filter performance in the sense that the filtering error e i ( t ) converges to zero faster. The estimated bound is presented as follows:
lim t 1 m E i = 1 m | η ( t ) η ^ i ( t ) | 2 ( 1 ρ ¯ ) r 0 + ( 1 ρ ¯ ) 2 r 0 2 β ^ ( ρ ¯ l 0 + κ λ max ( Π ) ) ( ρ ¯ l 0 + κ λ max ( Π ) ) 2 ,
where β ^ = β μ c 2 2 m i = 1 m p i k i 2 ( 1 ϵ ¯ i ) 2 .
Proof. 
Closed-loop spatio-temporal operator A ψ ( ζ i ( t ) ) is easily verified that it meets
| i = 1 m g , A ψ ( ζ i ( t ) ) h | = | i = 1 m g , A k i ϵ ¯ i Ψ * ( ζ i ( t ) ) Ψ ( ζ i ( t ) ) h | ε 0 m g h + K ϵ ¯ Ψ ( ζ ( t ) ) g , Ψ ( ζ ( t ) ) h ε 0 m g h + λ max ( K ϵ ¯ ) ϱ 2 μ 2 ( Ω ) g h = ε 1 g h
where ε 1 = ε 0 m + λ max ( K ϵ ¯ ) ϱ 2 μ 2 ( Ω ) > 0 , ϵ ¯ = diag { ϵ ¯ 1 , ϵ ¯ 2 , , ϵ ¯ m } , K = diag { k 1 , k 2 , , k m } , Ψ ( ζ ( t ) ) = ϱ μ ( Ω ) .
i = 1 m g , A ψ ( ζ i ( t ) ) g = i = 1 m g , A k i ϵ ¯ i Ψ * ( ζ i ( t ) ) Ψ ( ζ i ( t ) ) g ϑ 0 m g 2 + K ϵ ¯ Ψ ( ζ ( t ) ) g , Ψ ( ζ ( t ) ) g ϑ 0 m g 2 + λ min ( K ϵ ¯ ) | Ψ ( ζ ( t ) ) g | 2 > ϑ 0 m g 2
where g , h B .
By (28) and (29), the linear operator A ψ ( ζ i ( t ) ) is negative definite and invertible. Thus, we consider the following Lyapunov functional that depends on the spatio-temporal operator A ψ ( ζ i ( t ) ) .
V ( t ) = 1 2 i = 1 m e i ( t ) , A ψ ( ζ i ( t ) ) e i ( t ) .
The infinitesimal operator L V is defined as L V ( t ) = lim Δ 0 + 1 Δ [ E V ( t + Δ ) | t V ( t ) ] , then along the solution of the filtering error system (22), we have
L V ( t ) = i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) e ˙ i ( t ) 1 2 i = 1 m E e i ( t ) , d A ψ ( ζ i ( t ) ) d t e i ( t ) .
The following result was simply inferred by taking into account (10)–(12), as well as by noting that the operator A ψ ( ζ i ( t ) ) is self-adjoint.
i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) e ˙ i ( t ) = i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) A ψ ( ζ i ( t ) ) e i ( t ) i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) H ( F ( η ( t ) ) F ( η ^ i ( t ) ) ) + i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) Ψ * ( ζ i ( t ) ) k i ( ϵ i ( t ) ϵ ¯ i ) Ψ ( ζ i ( t ) ) η ( t ) + i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) Ψ * ( ζ i ( t ) ) k i ( 1 ϵ i ( t ) ) σ ( y i ( t ) ) + i = 1 m e i ( t ) , A ψ ( ζ i ( t ) ) G i κ j = 1 m π i j e j ( t ) i = 1 m e i ( t ) , A ψ ( ζ i ( t ) ) D ν ( t )
By Assumptions 3 and 4 together with (4)–(6), the following holds:
i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) H ( F ( η ( t ) ) F ( η ^ i ( t ) ) ) = i = 1 m E e i ( t ) , F ( η ( t ) ) F ( η ^ i ( t ) ) ρ ¯ l 0 i = 1 m | e i ( t ) | 2 + 2 ( 1 ρ ¯ ) r 0 i = 1 m | e i ( t ) |
where H = A ψ 1 ( ζ i ( t ) ) .
Here, selecting the consensus filter gains G i = A ψ 1 ( ζ i ( t ) ) out of simplicity, we obtain
i = 1 m e i ( t ) , A ψ ( ζ i ( t ) ) G i κ j = 1 m π i j e j ( t ) = κ i = 1 m e i ( t ) , j = 1 m π i j e j ( t ) = κ i = 1 m j = 1 m π i j | e i ( t ) | 2
Moreover, it is simple to deduce that
i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) D ν ( t ) i = 1 m q i 2 D ν ( t ) , D ν ( t ) + i = 1 m 1 2 q i A ψ ( ζ i ( t ) ) e i ( t ) , A ψ ( ζ i ( t ) ) e i ( t ) i = 1 m q i 2 μ d 2 | ν ( t ) | 2 + i = 1 m 1 2 q i | A ψ ( ζ i ( t ) ) e i ( t ) | 2
By Lemma 1, the following holds:
i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) Ψ * ( ζ i ( t ) ) k i ( 1 ϵ i ( t ) ) σ ( y i ( t ) ) i = 1 m p i 2 E Ψ * ( ζ i ( t ) ) k i ( 1 ϵ i ( t ) ) σ ( y i ( t ) ) , Ψ * ( ζ i ( t ) ) k i ( 1 ϵ i ( t ) ) σ ( y i ( t ) ) ) + i = 1 m 1 2 p i E A ψ ( ζ i ( t ) ) e i ( t ) , A ψ ( ζ i ( t ) ) e i ( t ) μ c 2 i = 1 m p i k i 2 ( 1 ϵ ¯ i ) 2 2 σ ( y i ( t ) ) , σ ( y i ( t ) ) + i = 1 m 1 2 p i | A ψ ( ζ i ( t ) ) e i ( t ) | 2
Substituting (34)–(37) into (33) leads to
i = 1 m E e i ( t ) , A ψ ( ζ i ( t ) ) e ˙ i ( t ) i = 1 m | A ψ ( ζ i ( t ) ) e i ( t ) | 2 + ρ ¯ l 0 i = 1 m | e i ( t ) | 2 + 2 ( 1 ρ ¯ ) r 0 i = 1 m | e i ( t ) | + μ c 2 i = 1 m p i k i 2 ( 1 ϵ ¯ i ) 2 2 σ ( y i ( t ) ) , σ ( y i ( t ) ) + i = 1 m 1 2 p i | A ψ ( ζ i ( t ) ) e i ( t ) | 2 + κ i = 1 m j = 1 m π i j | e i ( t ) | 2 + i = 1 m q i 2 μ d 2 | ν ( t ) | 2 + i = 1 m 1 2 q i | A ψ ( ζ i ( t ) ) e i ( t ) | 2 i = 1 m 1 + 1 2 p i + 1 2 q i | A ψ ( ζ i ( t ) ) e i ( t ) | 2 + ( ρ ¯ l 0 + κ λ max ( Π ) ) i = 1 m | e i ( t ) | 2 + 2 ( 1 ρ ¯ ) r 0 i = 1 m | e i ( t ) | + μ c 2 i = 1 m p i k i 2 ( 1 ϵ ¯ i ) 2 2 β + i = 1 m q i 2 μ d 2 d ¯ 2
Furthermore, with respect to the second part of (32), we obtain
1 2 i = 1 m E e i ( t ) , d A ψ ( ζ i ( t ) ) d t e i ( t ) = 1 2 i = 1 m k i ϵ ¯ i e i ( t ) , d Ψ * ( ζ i ( t ) ) Ψ ( ζ i ( t ) ) d t e i ( t ) = s = 1 3 i = 1 m k i ϵ ¯ i ψ i ( ζ i s ( t ) ) e i ( t ) , ζ ˙ i s ( t ) ) d ψ i ( ζ i s ( t ) ) d ζ s e i ( t ) = s = 1 3 i = 1 m k i ϵ ¯ i ζ ˙ i s ( t ) ) 0 l i s φ i ( ζ s ) d φ i ( ζ s ) d ζ s e i 2 ( t ) d ζ s = i = 1 m k i ϵ ¯ i ζ ˙ i s ( t ) s = 1 3 Ξ i s
where Ξ i s , s = 1 , 2 , 3 have been defined in Theorem 1.
The choice
ζ ˙ i s ( t ) = c i s k i ϵ ¯ i Ξ i s , s = 1 , 2 , 3 .
Deduce (39) negative definite, c i s > 0 , s = 1 , 2 , 3 ; i = 1 , 2 , , m . Notice that ( ζ ˙ i 1 ( t ) , ζ ˙ i 2 ( t ) , ζ ˙ i 3 ( t ) ) is the velocity of ith moving sensing device.
Subsequently, from (38) and (39), it follows that:
L V ( t ) i = 1 m ϖ i | A ψ ( ζ i ( t ) ) e i ( t ) | 2 + d ^ + ρ ¯ l 0 + κ λ max ( Π ) i = 1 m | e i ( t ) | 2 + 2 ( 1 ρ ¯ ) r 0 i = 1 m | e i ( t ) | + β μ c 2 2 i = 1 m p i k i 2 ( 1 ϵ ¯ i ) 2
where ϖ i = 1 + 1 2 p i + 1 2 q i and d ^ is defined in Theorem 1.
Note that
i = 1 m | e i ( t ) | 2 = i = 1 m j = 1 m | e i ( t ) | | e j ( t ) | 1 2 i = 1 m j = 1 m | e i ( t ) | 2 + | e j ( t ) | 2 = m i = 1 m | e i ( t ) | 2
Then, it follows that
L V ( t ) ϖ 0 i = 1 m | A ψ ( ζ i ( t ) ) e i ( t ) | 2 + d ^ + ρ ¯ l 0 + κ λ max ( Π ) i = 1 m | e i ( t ) | 2 + 2 ( 1 ρ ¯ ) r 0 m i = 1 m | e i ( t ) | 2 + β μ c 2 2 i = 1 m p i k i 2 ( 1 ϵ ¯ i ) 2 < ϖ V ( t ) + d ^ + ρ ¯ l 0 + κ λ max ( Π ) m [ 1 m i = 1 m | e i ( t ) | 2 + 2 ( 1 ρ ¯ ) r 0 ρ ¯ l 0 + κ λ max ( Π ) 1 m i = 1 m | e i ( t ) | 2 β ^ ρ ¯ l 0 + κ λ max ( Π ) ] .
where ϖ 0 = max { ϖ i } , i = 1 , 2 , , m , ϖ = 2 ϖ 0 ϑ 0 m and β ^ is defined in Theorem 1. Let χ = 1 m i = 1 m | e i ( t ) | 2 and w ( χ ) = χ 2 + 2 ( 1 ρ ¯ ) r 0 ρ ¯ l 0 + κ λ max ( Π ) χ β ^ ρ ¯ l 0 + κ λ max ( Π ) . It is simple to observe that w ( x ) = 0 has two solutions
χ 1 , 2 = ( 1 ρ ¯ ) r 0 ± ( 1 ρ ¯ ) 2 r 0 2 β ^ ( ρ ¯ l 0 + κ λ max ( Π ) ) ρ ¯ l 0 + κ λ max ( Π )
where χ 1 < 0 and χ 2 > 0 . If χ χ 2 , then w ( χ ) 0 .
Thus, it is not difficult to draw the following conclusion
L V ( t ) < ϖ V ( t ) + d ^ .
From (45), we can conclude that
V ( t ) < e ϖ T f V ( 0 ) + m d ^ = e ϖ T f 1 2 i = 1 m e i ( 0 ) , A ψ ( ζ i ( t ) ) e i ( 0 ) + m d ^ < e ϖ T f m α 1 2 + d ^ .
From the definition of V ( t ) , there exists 0 < λ 0 1 2 such that
V ( t ) λ 0 i = 1 m e i ( t ) , A ψ ( ζ i ( t ) ) e i ( t )
Therefore, we obtain
i = 1 m e i ( t ) , A ψ ( ζ i ( t ) ) e i ( t ) < 1 λ 0 e ϖ T f m α 1 2 + d ^ m α 2 2 .
This completes the proof.  □

3.2. Distributed H Consensus Performance Analysis

Next, we will concentrate on studying the H performance of the zero initial condition for the filtering error system (21).
Theorem 2.
Let the filter parameter k i and G i and the perturbation attenuation level γ > 0 be given. Then, the zero solution of the filtering error system (21) is finite-time bounded with respect to ( 0 , α 2 2 , T f , d ¯ , A ψ ( ζ i ( t ) ) ) , and z ˜ satisfies the bounded H consensus performance constraint (19) under the zero initial condition for all nonzero ν ( t ) , if under Assumptiosns 1–5 such that the following inequalities hold:
ρ ¯ l 0 + κ λ max ( Π ) + μ b 2 < 0 ,
γ + μ d 2 2 m i = 1 m q i < 0 ,
and the velocity law of mobile sensor is determined as (25) and (26). The estimated bound is given by
lim t 1 m E i = 1 m | η ( t ) η ^ i ( t ) | 2 ( 1 ρ ¯ ) r 0 + ( 1 ρ ¯ ) 2 r 0 2 β ^ l 1 l 1 2 .
Proof. 
It is simple to prove that (23)–(25) implies the filtering error system (21) is finite-time bounded. The H consensus performance of the closed-loop system will now be the subject of our discussion. The same functional Lyapunov candidate V ( t ) from Theorem 1 should be built. As in Theorem 1, similar line calculation yields
L V ( t ) ϖ m 2 α 2 2 + ρ ¯ l 0 + κ λ max ( Π ) m w ( χ ) + d ^ .
In order to address the system’s H consensus performance (19), we establish
J = E i = 1 m 0 T f | z ˜ i ( t ) | 2 γ 2 | ν ( t ) | 2 d t = E i = 1 m 0 T f z ˜ i ( t ) , z ˜ i ( t ) d t m γ 2 0 T f ν ( t ) , ν ( t ) d t .
From (45) to (46), with the zero initial condition, we have
J 0 T f i = 1 m B e i ( t ) , B e i ( t ) d t m γ 2 d ¯ 2 + L V ( t ) ϖ m 2 α 2 2 + l 1 m w ˜ ( χ ) + d ˇ ,
where l 1 = ρ ¯ l 0 + κ λ max ( Π ) + μ b 2 , w ˜ ( χ ) = χ 2 + 2 ( 1 ρ ¯ ) r 0 l 1 χ β ^ l 1 and d ˇ = ( γ 2 + μ d 2 2 m i = 1 m q i ) m d ¯ 2 .
It is easy to see that w ˜ ( χ ) = 0 has two solutions
χ 1 , 2 = ( 1 ρ ¯ ) r 0 ± ( 1 ρ ¯ ) 2 r 0 2 β ^ l 1 l 1 ,
where χ 1 < 0 and χ 2 > 0 . If χ χ 2 , then w ˜ ( χ ) > 0 . It is known with the condition (48) and (49) that l 1 m w ˜ ( χ ) + d ˇ < 0 .
Letting T f , we obtain that
i = 1 m 0 T f z ˜ i ( t ) , z ˜ i ( t ) d t < m γ 2 0 T f ν ( t ) , ν ( t ) d t + ϖ m 2 α 2 2 ,
namely,
1 m i = 1 m | z ˜ i ( t ) | 2 < γ 2 | ν ( t ) | 2 + ϖ 2 α 2 2 ,
which completes the proof of the theorem.  □
Remark 1.
A distributed consensus filter is built in (20) based on the mobile sensing approach. The filter can be adapted to the state estimation of complex environments. On the one hand, the nonlinear function in the system is assumed to have two parts, including the part that satisfies the sector condition and the part with saturation characteristics, and both parts occurred randomly. On the other hand, the randomly occurring saturation is considered at the output measurement. This leads to a better robustness of the filter.
Remark 2.
From (27), it follows that the minimal bound increases as r 0 or β increases. Thus, a very large filter gain should not be used if the bounds r 0 or β of the saturation is relatively high.

4. Discussion

Because of the generality of the conclusion given in Theorem 1, we discuss a few special scenarios of its application in the following.

4.1. Measurement Output with Saturation

If the measurement output is sensor saturation, the output can be expressed as
y ˇ i ( t ) = σ Ω ψ i ( ζ ; ζ i ( t ) ) η ( t , ζ ) d Ω , i = 1 , 2 , , m .
Its abstract form is
y ˇ ( t ) = σ ( Ψ ( ζ ( t ) ) η ( t ) ) .
Since the saturation function has obvious nonlinear characteristics, with the help of techniques for dealing with nonlinear terms in the system, we can decompose the saturation function into a sector nonlinear part and a high-order saturation part to facilitate our analysis.
σ ( y ( t ) ) = C 1 Ψ ( ζ ( t ) ) η ( t ) + R ( y ( t ) ) .
Considering Assumption 4, the high-order saturation part R ( y ( t ) ) satisfies
| R ( y ( t ) ) | r 1 .
For the nonlinear function f of the 3D system, take ρ ¯ = 1 ; hence f ( t , η ( t , ζ ) ) = f L ( t , η ( t , ζ ) ) satisfies Assumption 3.
To study the distributed H consensus filtering of the 3D nonlinear distributed parameter systems with sensor saturation, the distributed consensus filters (20) can be modified to the following form.
η ^ ˙ i ( t ) = A η ^ i ( t ) + H F ( η ^ i ( t ) ) + Ψ * ( ζ i ( t ) ) k i [ y ˇ i ( t ) C 1 Ψ ( ζ i ( t ) ) η ^ i ( t ) ] + G i κ j = 1 m π i j η ^ j ( t ) z ^ i ( t ) = B η ^ i ( t )
The filtering error system can be given from (13) and (59) as follows.
e ˙ i ( t ) = A ˜ ψ e i ( t ) + H ( F ( η ( t ) ) F ( η ^ i ( t ) ) ) Ψ * ( ζ i ( t ) ) k i R ( y i ( t ) ) G i κ j = 1 m π i j η ^ j ( t ) + D ν ( t ) z ˜ i ( t ) = B e i ( t )
where A ˜ ψ ( ζ i ( t ) ) = A k i C 1 Ψ * ( ζ i ( t ) ) Ψ ( ζ i ( t ) ) .
Similar to the proof of Theorem 1, let the Lyapunov function V ( t ) = 1 2 i = 1 m e i ( t ) , A ˜ ψ ( ζ i ( t ) ) e i ( t ) here. The following inequality
i = 1 m e i ( t ) , A ˜ ψ ( ζ i ( t ) ) Ψ * ( ζ i ( t ) ) k i R ( y i ( t ) ) r 1 i = 1 m ϑ ˜ i k i | e i ( t ) |
where ϑ ˜ i = ( ϑ 0 + k i l 1 μ 2 ( Ω ) ) μ ( Ω ) , are available when dealing with this Lyapunov functional.
Thus, it follows that
L V ( t ) ϖ m i = 1 m | A ˜ ψ ( ζ i ( t ) ) e i ( t ) | 2 + d ^ + l 0 + κ λ max ( Π ) i = 1 m | e i ( t ) | 2 r 1 λ max Θ ˜ K i = 1 m | e i ( t ) | + β ^ ,
where Θ ˜ = diag { ϑ ˜ 1 , ϑ ˜ 2 , , ϑ ˜ m } .
From the above, it is not difficult to draw the following corollary.
Corollary 1.
Let the filter parameter k i and G i and the perturbation attenuation level γ > 0 be given. Then, the zero solution of the filtering error system (60) is finite-time bounded with respect to ( 0 , α 2 2 , T f , d ¯ , A ψ ( ζ i ( t ) ) ) , and z ˜ satisfies the bounded H consensus performance constraint (19) under the zero initial condition for all nonzero ν ( t ) , if under Assumptions 1–5 such that the following inequalities hold:
l 0 + κ λ max ( Π ) + μ b 2 < 0 ,
γ + μ d 2 2 m i = 1 m q i < 0 ,
and the velocity law of the mobile sensor is determined as (25) and (26). The estimated bound is given by
lim t 1 m E i = 1 m | η ( t ) η ^ i ( t ) | 2 r 1 λ max Θ ˜ K + r 1 2 λ max 2 Θ ˜ K 4 β ^ l ^ 1 2 l ^ 1 2 .
where l ^ 1 = l 0 + κ λ max ( Π ) + μ b 2 .

4.2. Measurement Output with Packet Losses

If sensor saturation is not taken into account in the measurement output, then the stochastic variable ϵ i ( t ) in the output expression (65) implies the scenario of missing data in the measurement output.
y ¯ i ( t ) = ϵ i ( t ) Ω ψ i ( ζ ; ζ i ( t ) ) η ( t , ζ ) d Ω , i = 1 , 2 , , m .
Its abstract form can be expressed as
y ¯ ( t ) = ϵ ( t ) Ψ ( ζ ( t ) ) η ( t ) .
For the nonlinear function f of the 3D system, take ρ ¯ = 0 ; hence f ( t , η ( t , ζ ) ) = r ( t , η ( t , ζ ) ) satisfies Assumption 4.
Combined with the output (66), a modified distributed filter can be written as
η ^ ˙ i ( t ) = A η ^ i ( t ) + H F ( η ^ i ( t ) ) + Ψ * ( ζ i ( t ) ) k i [ y ¯ i ( t ) ϵ ¯ i Ψ ( ζ i ( t ) ) η ^ i ( t ) ] + G i κ j = 1 m π i j η ^ j ( t ) z ^ i ( t ) = B η ^ i ( t )
The filtering error system can be given from (13) and (67) in the following.
e ˙ i ( t ) = A ψ ( ζ i ( t ) ) e i ( t ) + H ( F ( η ( t ) ) F ( η ^ i ( t ) ) ) Ψ * ( ζ i ( t ) ) k i ( ϵ i ( t ) ϵ ¯ i ) Ψ ( ζ i ( t ) ) η ( t ) G i κ j = 1 m π i j η ^ j ( t ) + D ν ( t ) z ˜ i ( t ) = B e i ( t )
where A ψ ( ζ i ( t ) ) has been defined in (22).
Similar to the proof of Theorem 1, we can obtain the following corollary.
Corollary 2.
Let the filter parameter k i and G i and the perturbation attenuation level γ > 0 be given. Then, the zero solution of the filtering error system (21) is finite-time bounded with respect to ( 0 , α 2 2 , T f , d ¯ , A ψ ( ζ i ( t ) ) ) , and z ˜ satisfies the bounded H consensus performance constraint (19) under the zero initial condition for all nonzero ν ( t ) , if under Assumptions 1–5 such that the following inequalities hold:
κ λ max ( Π ) + μ b 2 < 0 ,
γ + μ d 2 2 m i = 1 m q i < 0 ,
and the velocity law of mobile sensor is determined as (25) and (26). The estimated bound is given by
lim t 1 m E i = 1 m | η ( t ) η ^ i ( t ) | 2 ( 1 ρ ¯ ) r 0 + ( 1 ρ ¯ ) 2 r 0 2 β ^ l 2 l 2 2 .
where l 2 = κ λ max ( Π ) + μ b 2 .

4.3. Measurement with Homogeneous Mobile Sensor Networks

In the following, the consensus filtering issue of system (13) is considered to be solved in a homogeneous mobile sensor network. The spatial distribution of the homogeneous mobile sensor network is given by
ψ i ( ζ s ; ζ i s ( t ) ) = ϕ if ζ Ω i 0 otherwise ,
where ζ 1 Ω i 1 = [ ζ i 1 δ i 1 , ζ i 1 + δ i 1 ] , ζ 2 Ω i 2 = [ ζ i 2 δ i 2 , ζ i 2 + δ i 2 ] , ζ 3 Ω i 3 = [ ζ i 3 δ i 3 , ζ i 3 + δ i 3 ] .
Furthermore, it is not difficult to establish that the system is mean square asymptotically stable if the condition 1 + 1 2 p i + 1 2 q i < ϖ i is attached in the (40) of proving the stability of the system.
The following corollary is easily obtained from Theorem 2; hence the proof is omitted.
Corollary 3.
Let the filter parameter k i and G i and the perturbation attenuation level γ > 0 be given. Then, the zero solution of the filtering error system (21) is mean square asymptotically stable, and z ˜ satisfies the H consensus performance constraint (19) under the zero initial condition for all nonzero ν ( t ) , if under Assumptions 1–5 such that the following inequalities hold:
1 + 1 2 p i + 1 2 q i < ϖ i , i = 1 , 2 , , m ,
ρ ¯ l 0 + κ λ max ( Π ) + μ b 2 < 0 ,
γ + μ d 2 2 m i = 1 m q i < 0 ,
and the velocity law of mobile sensor is given by
ζ ˙ i s ( t ) = c ^ i s k i ϵ ¯ i Ξ ^ i s , s = 1 , 2 , 3 .
where
Ξ ^ i s = ϕ 2 e i 2 ( t , ζ i s δ i s ) e i 2 ( t , ζ i s + δ i s ) , s = 1 , 2 , 3 .
with c ^ i s > 0 , s = 1 , 2 , 3 , i = 1 , 2 , , m is velocity gain of each mobile sensor. The estimated bound is determined as (50).

4.4. Point Measurement

In engineering, point measurements are easier to achieve than distributed measurements. Changing the spatial distribution of the mobile sensing device to Dirac delta function in (4), i.e., by selecting
ψ i ( ζ s ; ζ i s ( t ) ) = δ ( ζ ζ 1 ) δ ( ζ ζ 2 ) δ ( ζ ζ 3 ) ,
the measurement of the 3D distributed parameter system changes from distributed to point measurements.
The following Corollary 4 is easily obtained from Theorem 2 and the spatial distribution of the mobile sensors (78); hence the proof is omitted.
Corollary 4.
Let the filter parameter k i and G i and the perturbation attenuation level γ > 0 be given. Then, the zero solution of the filtering error system (21) is finite-time bounded with respect to ( 0 , α 2 2 , T f , d ¯ , A ψ ( ζ i ( t ) ) ) , and z ˜ satisfies the bounded H consensus performance constraint (19) under the zero initial condition for all nonzero ν ( t ) , if under Assumptions 1–5 such that the inequalities (48) and (49) hold, and the velocity law of mobile sensor is given by
ζ ˙ i s ( t ) = c ˇ i s k i ϵ ¯ i Φ i s , s = 1 , 2 , 3 .
where
Φ i s = δ i s e i 2 ( t , ζ i s δ i s ) e i 2 ( t , ζ i s + δ i s ) , s = 1 , 2 , 3 .
with c ˇ i s > 0 , s = 1 , 2 , 3 , i = 1 , 2 , , m is velocity gain of each mobile sensor. The estimated bound is determined as (50).
Remark 3.
If the stochastic variable in the measurement output is taken as ϵ ¯ i = 1 and the measurement output comes from a fixed sensor, then the output expression reduces to y ψ ( t ) = Ψ η ( t , ζ ) . The distributed consensus filters (20) can be simplified as
η ^ ˙ i ( t ) = A η ^ i ( t ) + H F ( η ^ i ( t ) ) + k i [ y ψ i ( t ) Ψ η ^ i ( t ) ] + G i κ j = 1 m π i j η ^ j ( t ) z ^ i ( t ) = B η ^ i ( t )
It is simple to deduce the following conclusions: Let the filter parameter k i and G i and the perturbation attenuation level γ > 0 be given. Then, the zero solution of the filtering error system (21) is finite-time bounded with respect to ( 0 , α 2 2 , T f , d ¯ , A ψ ) , and z ˜ satisfies the bounded H consensus performance constraint (19) under the zero initial condition for all nonzero ν ( t ) , if under Assumptions 1–5 such that the following inequalities (49) and ρ ¯ l 0 + κ λ max ( Π ˜ ) + μ b 2 < 0 where Π ˜ = Π K Ψ hold. The estimated bound is given by lim t 1 m E i = 1 m | η ( t ) η ^ i ( t ) | 2 ( 1 ρ ¯ ) r 0 + ( 1 ρ ¯ ) 2 r 0 2 β ^ l ˜ 1 l ˜ 1 2 where l ˜ 1 = ρ ¯ l 0 + κ λ max ( Π ˜ ) + μ b 2 .

5. Simulation Example

In this section, simulations are given to illustrate how these results can be applied to achieve finite-time bounded in 3D filtering error system. With Dirichlet boundary conditions, we take into consideration a class of 3D nonlinear distributed parameter systems subject to the initial condition η ( 0 , ζ 1 , ζ 2 , ζ 3 ) = ( ζ 1 / l 1 ) 3 ( 1 ζ 1 / l 1 ) 3 ( ζ 2 / l 2 ) 3 ( 1 ζ 2 / l 2 ) 3 ( ζ 3 / l 3 ) 3 ( 1 ζ 3 / l 3 ) 3 , where l 1 = 1 , l 2 = 0.5 , l 3 = 0.8 . The system parameters are proposed as follows: a 0 = 10 5 and f ( t , η ( t , ζ ) ) = arctan ( 0.7 η ( t , ζ ) ) . The evolution of the 3D nonlinear distributed parameter system at four different times is shown in Figure 2.
The coordinate system in Figure 2 represents a three-dimensional space. Hence, the states of the system in three dimensions at time t = 0.25 , 1 , 2 , 5 s, respectively, need to be shown using four different spatial coordinate systems.
Two mobile sensing devices are taken into account, and their initial locations are selected as ( ζ 11 ( 0 ) , ζ 12 ( 0 ) , ζ 13 ( 0 ) ) = ( 0.18 , 0.10 , 0.25 ) and ( ζ 21 ( 0 ) , ζ 22 ( 0 ) , ζ 23 ( 0 ) ) = ( 0.56 , 0.14 , 0.32 ) . The sensor network coupling matrix Π = 1 1 1 1 . Their spatial distribution is defined by
ψ i ( ζ s ; ζ i s ( t ) ) = 1 if ζ Ω i 0 otherwise .
In this simulation, the probabilities are taken as ρ ¯ = 0.75 , ϵ ¯ 1 = 0.9 and ϵ ¯ 2 = 0.85 ; the perturbation attenuation level is γ = 0.95 ; σ ( y i ( t ) ) is a saturation function written in the following:
σ ( y i ( t ) ) = σ ( y i ( t ) ) = V i , max , if y i ( t ) > V i , max σ ( y i ( t ) ) = y i ( t ) , if V i , max y i ( t ) V i , max σ ( y i ( t ) ) = V i , max , if y i ( t ) < V i , max
The saturation values are taken as V 1 = V 2 = 0.07 . The measurement output of mobile sensors with randomly sensor saturation are depicted in Figure 3.
The initial condition for the distributed consensus filter is selected to be η ^ 1 ( 0 , ζ ) = η ^ 2 ( 0 , ζ ) = 0 . Gain for the distributed filters is determined by k 1 = 10 and k 2 = 30 . Figure 4 and Figure 5 illustrate how the filtering errors of the two distributed H consensus filters changed at four different times in the first 5 s. After 5 s, with increasing time, their filtering errors converge rapidly to the bounds. Figure 6 and Figure 7 show the estimation results of the two distributed filters at four different times, respectively.
Figure 8 displays the output estimation errors for the two filters. The estimation converges to the bounded range in a short time, as can be seen in Figure 8. The existence of the estimation error is caused by the bounded perturbation in the system. Once the perturbation decays or even disappears, then the estimation error converges rapidly to zero.
To reflect the enhanced filtering performance of the motion sensing approach, we introduce fixed sensors for comparison with mobile sensors. The two fixed sensors are located at ( ζ 11 , ζ 12 , ζ 13 ) = ( 0.18 , 0.10 , 0.25 ) and ( ζ 21 , ζ 22 , ζ 23 ) = ( 0.56 , 0.14 , 0.32 ) .These two locations are indicated by small circles in Figure 9. These two spatial locations are also the initial locations before the mobile sensor moves. As a comparison with the fixed sensor, the trajectories of the two mobile sensors under velocity constraints are also shown in Figure 9.
In Figure 10, the effect of the evolution of the L 2 norm of the filtering error for the two fixed sensors and the two mobile sensors is shown. The appreciable reduction of the state norm is observed for the mobile case.

6. Conclusions

This paper has addressed the issue of finite-time distributed H consensus filtering for a class of 3D nonlinear distributed parameter systems subject to bounded perturbation. The considered incomplete information in measured output is randomly occurring sensor saturation. Specialized conditions are used to treat the system and possible output saturation. A mobile sensing approach is proposed as a new framework for optimizing the filtering issue, thus enhancing the performance of the filter. Moreover, a novel nonlinear stochastic decomposition model is proposed. With a guaranteed H consensus disturbance rejection attenuation level and dynamics of the filtering error system that are finite-time bounded, a novel distributed H consensus filter has been developed. In light of the operator-dependent Lyapunov function, the finite-time bounded of the filtering error system can be achieved. For the simplified analysis, the spatial distribution of the selected mobile sensing devices is symmetric, and therefore, the closed-loop operator in the evolution equation is self-adjoint. In a more general case, the spatial distribution may not be symmetric, and even may not be derivative with respect to the location. Then the analysis of the closed-loop operator will be a challenge and is being considered by the authors. The filtering error system is able to converge faster, also dependent on conditions in the form of the velocity law for the mobile sensors. It has been demonstrated through numerical simulations that the proposed criteria are effective.

Author Contributions

Conceptualization, X.Q. and B.C.; methodology, X.Q.; software, X.Q.; formal analysis, X.Q.; resources, B.C.; writing—original draft preparation, X.Q.; writing—review and editing, X.Q.; supervision, B.C.; project administration, B.C.; funding acquisition, X.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Natural Science Foundation of the Jiangsu Higher Education Institutions of China (No. 17KJB510051), Qing Lan Project of the Jiangsu Higher Education Institutions (Young and Middle-aged Academic Leader (2022), Excellent Teaching Team (2020)), Soft Science Research Project of Wuxi (No. KX-22-B60), Advanced Research and Study Project for Academic Leaders of Jiangsu Higher Vocational Colleges (No. 2021GRFX068).

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.

Abbreviations

The following abbreviations are used in this manuscript:
3DThree-dimensional
DPSDistributed parameter systems
LMILinear matrix inequality
LOILinear operator inequality

References

  1. Cheng, Z.; Xu, F.; Luo, X. Analysis of spatial aggregation time lag in dynamic simulation of acetylene hydrogenation reactor distribution parameters and its influence on control performance. Chem. Autom. Instrum. 2022, 49, 274–279, 375. [Google Scholar]
  2. Makino, S.; Fukui, T.; Yoshida, T.; Hatsugai, Y. Edge states of a diffusion equation in one dimension: Rapid heat conduction to the heat bath. Phys. Rev. E 2022, 105, 024137. [Google Scholar] [CrossRef]
  3. Luo, K.; Wang, Q.; Li, Y.; Li, J.; Zhao, W. Research progress of magnetic fluid flow control based on high temperature gas effect. J. Mech. 2021, 53, 1515–1531. [Google Scholar]
  4. Liu, Z.; Ouyang, Y.; Song, Y.; He, W.; Wang, F. Parallel control of distributed parameter systems: From model-based control to data-driven intelligent control. J. Command. Control 2017, 3, 177–185. [Google Scholar]
  5. Liu, J.; Guo, B. Uniformly semidiscretized approximation for exact observability and controllability of one-dimensional Euler–Bernoulli beam. Syst. Control Lett. 2021, 156, 105013. [Google Scholar] [CrossRef]
  6. Li, Y.; Chen, C.; Ou, Y.; Bu, Z. Stochastic stability analysis of Markovian jump uncertain distributed parameter system. Control Eng. 2022, 29, 109–113. [Google Scholar]
  7. Bao, L.; Huang, P.; Fei, S. Exponential stability and L2 gain analysis of switched distributed parameter systems with time delay. J. Control Decis. 2020, 2, 385–393. [Google Scholar] [CrossRef]
  8. Guerrero-Sánchez, M.E.; Hernández-González, O.; Valencia-Palomo, G.; López-Estrada, F.R.; Rodríguez-Mata, A.E.; Garrido, J. Filtered observer-based ida-pbc control for trajectory tracking of a quadrotor. IEEE Access 2021, 9, 114821–114835. [Google Scholar] [CrossRef]
  9. Chen, L.; Jin, P.; Yang, J.; Li, Y.; Song, Y. Robust Kalman filter-based dynamic state estimation of natural gas pipeline networks. Math. Probl. Eng. 2021, 2021, 5590572. [Google Scholar] [CrossRef]
  10. Xie, J.; Dubljevic, S. Discrete-time Kalman filter design for linear infinite-dimensional systems. Processes 2019, 7, 451. [Google Scholar] [CrossRef]
  11. Wang, F.; Wang, Z.; Liang, J.; Liu, X. Robust finite-horizon filtering for 2-D systems with randomly varying sensor delays. IEEE Trans. Syst. Man Cybern. Syst. 2018, 50, 142–149. [Google Scholar] [CrossRef]
  12. Fridman, E.; Orlov, Y. An LMI approach to H boundary control of semilinear parabolic and hyperbolic systems. Automatica 2009, 45, 2060–2066. [Google Scholar] [CrossRef]
  13. Luo, Y.; Xia, W.; Liu, G.; Deng, F. W(1,2) (Ω)- and X(1,2) (Ω)-stability of delayed cellular neural network models with reaction diffusion terms. Chin. Sci. (Ser. F Inf. Sci.) 2009, 39, 608–616. [Google Scholar]
  14. Afshar, S.; Germ, F.; Morris, K.A. Extended Kalman filter based observer design for semilinear infinite-dimensional systems. arXiv 2022, arXiv:2202.07797. [Google Scholar]
  15. Gómez-Peñate, S.; López-Estrada, F.R.; Valencia-Palomo, G.; Rotondo, D.; Guerrero-Sánchez, M.E. Actuator and sensor fault estimation based on a proportional multiple-integral sliding mode observer for linear parameter varying systems with inexact scheduling parameters. Int. J. Robust Nonlinear Control 2021, 31, 8420–8441. [Google Scholar]
  16. Am, N.B.; Fridman, E. Network-based H filtering of parabolic systems. Automatica 2014, 50, 3139–3146. [Google Scholar]
  17. Yan, Z.; Guo, T.; Zhao, A.; Kong, Q.; Zhou, J. Reliable exponential H filtering for a class of switched reaction-diffusion neural networks. Appl. Math. Comput. 2022, 414, 126661. [Google Scholar] [CrossRef]
  18. Seo, J.; Yu, M.J.; Park, C.G.; Lee, J.G. An extended robust H filter for nonlinear constrained uncertain systems. IEEE Trans. Signal Process. 2006, 54, 4471–4475. [Google Scholar] [CrossRef]
  19. Gao, H.; Lam, J.; Wang, C. Robust energy-to-peak filter design for stochastic time-delay systems. Syst. Control. Lett. 2006, 55, 101–111. [Google Scholar] [CrossRef]
  20. Jiang, X.; Xia, G.; Feng, Z.; Jing, X. Distributed sampled-data nonfragile consensus filtering over sensor networks with topology switching and transmission delay. IEEE/ASME Trans. Mechatron. 2022, 27, 1379–1390. [Google Scholar] [CrossRef]
  21. Gao, R.; Yang, G.H. Distributed multi-rate sampled-data H consensus filtering for cyber-physical systems under denial-of-service attacks. Inf. Sci. 2022, 587, 607–625. [Google Scholar] [CrossRef]
  22. Zhu, F.; Liu, X.; Peng, L. Adaptive consensus-based distributed H filtering with switching topology subject to partial information on transition probabilities. Appl. Math. Comput. 2021, 411, 126534. [Google Scholar] [CrossRef]
  23. Dorato, P. Short-Time Stability in Linear Time-Varying Systems; Polytechnic Institute of Brooklyn: New York, NY, USA, 1961; pp. 33–39. [Google Scholar]
  24. Xiang, Z.; Sun, Y.N.; Mahmoud, M.S. Robust finite-time H control for a class of uncertain switched neutral systems. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1766–1778. [Google Scholar] [CrossRef]
  25. Guan, W.; Fu, L.; Ma, Y. Finite-time H filtering for discrete-time singular Markovian jump systems with time delay and input saturation. Complexity 2019, 2019, 4675397. [Google Scholar] [CrossRef]
  26. Cui, W.; Sun, S.; Fang, J.A.; Xu, Y.; Zhao, L. Finite-time synchronization of Markovian jump complex networks with partially unknown transition rates. J. Frankl. Inst. 2014, 351, 2543–2561. [Google Scholar] [CrossRef]
  27. Yang, X.; Lu, J. Finite-time synchronization of coupled networks with Markovian topology and impulsive effects. IEEE Trans. Autom. Control 2015, 61, 2256–2261. [Google Scholar] [CrossRef]
  28. Ma, N.; Liu, Z.; Chen, L. Robust and non-fragile finite time H synchronization control for complex networks with uncertain inner coupling. Comput. Appl. Math. 2018, 37, 5395–5409. [Google Scholar] [CrossRef]
  29. Wang, J.; Ding, C.; Pian, J.; Shi, Y. H model predictive control of discrete stochastic nonlinear systems with saturation. J. Harbin Univ. Technol. 2015, 20, 37–44. [Google Scholar]
  30. Curtain, R.F.; Zwart, H. An Introduction to Infinite-Dimensional Systems Theory: A State-Space Approach; Springer: New York, NY, USA, 2020; pp. 654–659. [Google Scholar]
Figure 1. The filtering issue in mobile sensor networks.
Figure 1. The filtering issue in mobile sensor networks.
Mathematics 10 03134 g001
Figure 2. The evolution of the 3D nonlinear distributed parameter systems.
Figure 2. The evolution of the 3D nonlinear distributed parameter systems.
Mathematics 10 03134 g002
Figure 3. Sensor measurement output.
Figure 3. Sensor measurement output.
Mathematics 10 03134 g003
Figure 4. Evolution of the filtering error system for the first filter.
Figure 4. Evolution of the filtering error system for the first filter.
Mathematics 10 03134 g004
Figure 5. Evolution of the filtering error system for the second filter.
Figure 5. Evolution of the filtering error system for the second filter.
Mathematics 10 03134 g005
Figure 6. Estimation of first filter.
Figure 6. Estimation of first filter.
Mathematics 10 03134 g006
Figure 7. Estimation of second filter.
Figure 7. Estimation of second filter.
Mathematics 10 03134 g007
Figure 8. Filter output estimation errors.
Figure 8. Filter output estimation errors.
Mathematics 10 03134 g008
Figure 9. Moving sensing device trajectory.
Figure 9. Moving sensing device trajectory.
Mathematics 10 03134 g009
Figure 10. Evolution of spatial L 2 norm.
Figure 10. Evolution of spatial L 2 norm.
Mathematics 10 03134 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qian, X.; Cui, B. Mobile Sensor Networks for Finite-Time Distributed H Consensus Filtering of 3D Nonlinear Distributed Parameter Systems with Randomly Occurring Sensor Saturation. Mathematics 2022, 10, 3134. https://doi.org/10.3390/math10173134

AMA Style

Qian X, Cui B. Mobile Sensor Networks for Finite-Time Distributed H Consensus Filtering of 3D Nonlinear Distributed Parameter Systems with Randomly Occurring Sensor Saturation. Mathematics. 2022; 10(17):3134. https://doi.org/10.3390/math10173134

Chicago/Turabian Style

Qian, Xueming, and Baotong Cui. 2022. "Mobile Sensor Networks for Finite-Time Distributed H Consensus Filtering of 3D Nonlinear Distributed Parameter Systems with Randomly Occurring Sensor Saturation" Mathematics 10, no. 17: 3134. https://doi.org/10.3390/math10173134

APA Style

Qian, X., & Cui, B. (2022). Mobile Sensor Networks for Finite-Time Distributed H Consensus Filtering of 3D Nonlinear Distributed Parameter Systems with Randomly Occurring Sensor Saturation. Mathematics, 10(17), 3134. https://doi.org/10.3390/math10173134

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