A mobile sensing approach to distributed consensus filtering of 2D stochastic nonlinear parabolic systems with disturbances

The focus of this paper is on distributed consensus filtering in two-dimensional (2D) space for stochastic nonlinear parabolic systems with disturbances. Distributed filters are used in mobile sensor networks to achieve consensus estimate. By utilizing a mobile sensing approach, the optimal framework for parabolic systems to improve filtering performance is provided. Sufficient conditions are created under the filtering error is boundedness by employing operator-dependent Lyapunov functional. The velocity law of each mobile sensor can also be used to guide the filtering error system toward rapid convergence. Finally, numerical examples are used to verify the effectiveness of the suggested approach.


Introduction
In the real world, most physical processes can be viewed as the evolution of parabolic systems.A number of applications employing this kind of system have started attracting growing interest, such as heat conduction, toxic gas diffusion, large-scale population movements and other processes in which evolution is a function of time and space.Observer-based control and signal processing are both fundamentally concerned with filtering parabolic systems.
As early as the early 1960s, the Kalman type filter, coupled with a weighted least squares approach (Meditch, 1971), variance method (Sakawa, 1972), and other techniques were used to analyse several filtering problems of linear parabolic systems to achieve optimal state estimation.After that, we extend the issue of filtering in parabolic systems to discrete cases, and we apply them to the study of air pollution (Omatu & Seinfeld, 1981).
Until date, the Kalman filter has remained one of the most effective methods for estimating measurement data.For networked control systems that gather information using a centralized structure, extended Kalman filtering (Hu et al., 2012), finite-horizon filtering (Wang et al., 2005) and H ∞ filtering (He et al., 2007) have all been extensively studied.Such a design of a centralized observer for one-dimensional heat conduction processes (Hidayat et al., 2011), even their fractional-order CONTACT Xueming Qian xmqian@protonmail.comprocesses (Chen et al., 2020), has also yielded meaningful results.
Unfortunately, the communication between nodes in sensor networks has limitations in practical applications, owing to the limitations of energy, bandwidth and range.Data fusion can be carried out in a distributed way, which means that each node can only communicate with its neighbours.For distributed strategies, switched system approach (Zhang et al., 2016), event-triggered mechanism (Li et al., 2020), fuzzy system method (Zhang et al., 2017) and Round-Robin type protocol (Ugrinovskii & Fridman, 2014) have been raised to study the distributed filtering issue for networked control systems.
Distributed estimation has now been the major tool for investigating filtering challenges involving stochastic parameters and nonlinearities in sensor networks with incomplete data.Sensor networks with incomplete data were used to investigate distributed filtering challenges for stochastic systems (Ding et al., 2012;Yu et al., 2013).In addition, Liu et al. (2015) addressed a recursive distributed filter using an event-triggered mechanism in which each sensor node communicates intermittently across sensor networks.
The consensus performance criterion was employed as a significant index to express the bounded consensus in terms of filtering error in Shen, Wang, and Hung (2010) and Yu et al. (2009).Distributed estimation is effective for consensus.Distributed consensus estimation combined with receding horizon techniques (Xu et al., 2022) has also been deeply investigated in distributed parameter systems.
Parabolic systems can be affected by noise for realistic physical processes.The disturbances might appear probabilistically.A randomly occurring nonlinearities (RONs) model is used to represent the nonlinear terms, which occur randomly, in the system.Numerous studies (Duan et al., 2014;Liu et al., 2011;Shen et al., 2011) on RONs have been suggested.
A networked environment is where data missing phenomena often arises due to sensor temporal failure or network transmission delays.In Nahi (1969), a method to deal with the data missing by using a Bernoulli distribution has been first presented and then investigated in Wei et al. (2009) and Shen, Wang, and Hung (2010) for the various filtering issues of networked control systems with probabilistic packet losses.Furthermore, signal quantization frequently occurs prior to transmission in a networked setting due to the constrained network bandwidth.It was noted that logarithmic quantization is suitable for communication with fewer bits.The sector bound technique has been extensively adopted for the quantization effects in Shen, Wang, Shu, et al. (2010).
As mentioned previously, most of the current filtering issues for systems adopt static sensor networks.Mobile sensor networks are receiving a lot of attention in the engineering industry since for their advantages of flexibility and energy reduction.For stochastic nonlinear parabolic systems using mobile sensor networks, the distributed consensus filtering challenge has not been well studied and is yet open.
The filtering issue for stochastic nonlinear parabolic systems with disturbances will be the subject of this study.Distributed consensus filters are constructed in two-dimensional (2D) space with the aid of mobile sensors.In this study, mobile sensing as a novel method will yield some innovative findings for the filtering issue of 2D parabolic systems.

Problem formulation and preliminaries
Consider a stochastic nonlinear parabolic system that describes a target process on a 2D plane.
where Q(t, θ) denotes the state of the 2D parabolic system at moment t ∈ R + , located at plane θ = (θ 1 , θ 2 ).The where γ ∈ [0, 1] are known positive constants.It is assumed that the stochastic variables γ (t) is independent of both σ (t) and the initial state of system.Hence, we have E{γ (t)} = γ .Outputs from N mobile sensors measured spatially.In this case, ith mobile sensor can provide output as follows: where d i (t) is external noise of ith mobile sensor, and ) refers to the spatial distribution of ith sensing device.It is usually assumed that distribution function ϕ i (θ ; θ i (t)) is non-negative and bounded and is expressed by where The distribution function shows that the distribution of each mobile sensor may also be different.Each moving sensing device is located at (θ i1 (t), θ i2 (t)).The timevarying location function changes with time, which characterizes the spatiotemporal trajectory of each mobile sensor.For every i, the stochastic variable where ᾱi ∈ [0, 1], i = 1, 2, . . ., N are known positive constants.It is assumed that the stochastic variables α i (t) is independent of both d i (t) and the initial state of system.Accordingly, we have In addition, the quantization effects are taken into account in this study.The quantizer q(•) is defined as would then be sent via the filter.The set of quantization levels for each q i (•) is provided by where ρ (i) represents the quantization density and is a constant.The logarithmic type quantizer q(•) is provided by After quantization, the measurements can be expressed as As a way of rewriting the 2D parabolic system (1) in the abstract form, a few symbols and their meanings are given here.
Let H be a Hilbert space.The abstract form of the 2D parabolic system is then as follows: where the state of system ) and its domain is given by D g, −Ag ≥ β 0 g 2 , for g ∈ B and constant β 0 > 0.
The measurement output can be expressed by The nonlinear function f satisfies the following Assumption 2.1.

Assumption 2.1:
There exists a constant l such that holds, for all s 1 , s 2 ∈ L 2 (S).
Assumption 2.3: It is to be noted that the observation operator is self-adjoint due to the symmetry of the spatial distribution ϕ i (θ ; θ i (t)), and its norm is determined by the embedding constant b and the measure of the spatial domain L(S).
Definition 2.1 is presented in the following to get the main results.

Definition 2.1 (Distributed consensus performance):
The filters are said to be distributed consensus filters if there exist η i > 0 and μ > 0 such that

Main results and proofs
The distributed consensus filter is designed for the ith mobile sensor in this paper: where Q i (t) is the ith mobile sensor's state estimation.ȳi (t), i = 1, 2, . . ., N is output measurement of the ith mobile sensor.As the adjoint of Φ * (θ (t)), Φ(θ(t)) performs the observation operator.k i > 0, i = 1, 2, . . ., N is the observer gain and G i , i = 1, 2, . . ., N is the consensus filter gain.In addition, It is obvious that B has an eigenvalue of zero, and all other eigenvalues are negative.
Letting e i (t) = Q i (t) − Q(t), the filtering error dynamical system can be expressed as follows: where Given that Φ(θ i (t)) is self-adjoint, it follows that the closed-loop operator A φ (θ i (t)) is also self-adjoint.Moreover, it can be known that the A φ (θ i (t)) is invertible when the boundedness and coercivity of A and Assumption 2.3 are considered combined.

Bounded consensus analysis
Theorem 3.1: Under Assumptions 2.1-2.3, the filters (14) are distributed consensus filters if and the velocity law of ith mobile sensor as follows: ) with ρ i1 > 0, ρ i2 > 0 are velocity gain of each mobile sensor.The estimated bound is presented as follows: Proof: Taking account into the properties of the linear operator A, the closed-loop operator A φ (θ i (t)) is easy to verify that it satisfies Consider the following Lyapunov functional, which is operator dependent: In the It ô formula, a stochastic differentiation of (20) results in A weak infinitesimal operator defined as L V(t) = lim →0 + E{V(t+ )|t}−V(t) along stochastic diffusion processes, according to filtering error system (15), yields In ( 22), a first term is deduced as follows: b ij e i (t), e j (t)) where C = G i = −A −1 φ (θ i (t)).By calculating the second term in ( 22), we obtain where ij are defined in Theorem 3.1.The choice deduces ( 24) negative definite and ρ ij > 0.
The last term of ( 22) gives the following: After combining ( 23)-( 26) the calculations above, we get It follows from the It ô formula that , then EV(t) − EV(0) < 0. This completes the proof.

Some special cases
Due to the generality of the conclusion given in Theorem 3.1, we consider below several special cases of its application.

Point measurement
From the distribution function (4), it follows that where When ε i1 and ε i2 small enough, we have Then, the spatial distribution of the mobile sensing device in (4) becomes a Dirac delta function.By choosing the measurement of the parabolic system changes from a distribution measurement to a point measurement.Point measurements are easier to implement in engineering than distributed measurements.
The following Theorem 3.2 is easily obtained from Theorem 3.1 and the spatial distribution of the mobile sensors (31), hence the proof is omitted.
Theorem 3.2: Under Assumptions 2.1-2.3, the filters ( 14) are distributed consensus filters if ( 16) holds and the velocity law of ith mobile sensor as follows: )) with ρi1 > 0, ρi2 > 0 are velocity gain of each mobile sensor.The estimated bound is presented as (19).

2D nonlinear parabolic systems with disturbances
If the nonlinear term and measurement output data missing are not random occur in the system (1) and output quantization phenomenon is not considered, the parabolic system can be simplified to the following form: The distributed consensus filter ( 14) can be reduced down to It is simple to derive the following corollary when Theorem 3.1 is taken into account.
Corollary 3.3: Under Assumptions 2. 1-2.3, the filters (36) are distributed consensus filters if l + κλ max (B) < 0, (37) and the velocity law of ith mobile sensor as follows: where i1 and i2 are defined in Theorem 3.1, ρi1 > 0, ρi2 > 0 are velocity gain of each mobile sensor.The estimated bound is presented as follows: Remark 3.1: For 2D nonlinear parabolic systems with disturbances, if a fixed sensor network is taken into consideration, the distributed consensus filters can be represented as Figure 1.The states of 2D parabolic system at four different moments.
while the system output can be expressed as yi (t) = ΦQ(t) + d i (t) dν i dt , i = 1, 2, . . ., N. The following conclusions are easy to draw: Under Assumptions 2.1-2.3, the filters (41) are distributed consensus filters if l + κλ max ( B) < 0, where B = B − KΦ.The estimated bound is presented as follows: .

2D stochastic nonlinear parabolic systems without disturbances
If noise or disturbances are not considered in the system (1) and its measured output, the parabolic system can be rewritten in the following form: Thus, the following corollary is easily introduced from Theorem 3.1.
Corollary 3.4: Under Assumptions 2.1-2.3, the filtering error system is mean square asymptotically stable, if ( 16) holds and the velocity law of ith mobile sensor as ( 17) and ( 18).The filter performance of mobile sensors is improved by the guidance scheme in the sense that the filtering error e i (t) converges to zero faster.

Numerical examples
To illustrate the benefit of the distributed consensus filter proposed in this paper, we provide a simulated scenario in this section.Take into consideration an array of 2D spatially distributed processes with Q(0, θ) = sin( π 2 θ) e −9θ 2 , θ = [0, 1] × [0, 1] as their initial condition.Dirichlet boundary conditions are followed by the system.Figure 1 shows the states of 2D spatial processes at four different times.
Figures 5-7 depict the evolution of the filtering error system for the three distributed consensus filters.
Figure 8 shows the evolution of L 2 norm for the three filtering error systems.

Conclusions
A distributed consensus filter is designed in this paper for the consensus estimation using mobile sensor networks.Mobile sensing has been proposed as a novel optimal framework for improving filter performance in 2D parabolic systems with stochastic nonlinear function and bounded disturbances.The measurement output of each mobile sensor includes external noise.These noises are Brownian motions with a stochastic feature that meet the bounded sector assumption.The convergence of the filtering error system is affected by stochastic nonlinear and disturbance terms.Filtering errors only converge within certain bounds.Sufficient criterion for the uniformly bounded of the filtering error system has been established which used the operator-dependent Lyapunov functional, and an estimated bound has been presented.In addition, the velocity law of each mobile sensor can be adjusted such that the filtering error converges more faster to its bounds.The optimization framework proposed in this study under the mobile sensing   approach allows for optimal performance of the filter design while explicitly incorporating the location of the mobile sensors into the process dynamics under this framework.

Figure 2 .
Figure 2. The measurement output of first mobile sensors with random data missing.

Figure 3 .
Figure 3.The measurement output of second mobile sensor with random data missing.

Figure 4 .
Figure 4.The measurement output of third mobile sensors with random data missing.

Figure 5 .
Figure 5.The evolution of filtering error systems for filter 1.

Figure 6 .
Figure 6.The evolution of filtering error systems for filter 2.

Figure 7 .
Figure 7.The evolution of filtering error systems for filter 3.

Figure 8 .
Figure 8.The evolution of L 2 norm for filtering error systems.

Figure 9 .
Figure 9.The trajectories of three mobile sensing devices.

Figure 10 .
Figure 10.The trajectory of first mobile sensor.

Figure 11 .
Figure 11.The trajectory of second mobile sensor.

Figure 12 .
Figure 12.The trajectory of third mobile sensor.
The inner product of H is •, • and its induced norm is | • |.Let B be a reflexive Banach space.• denotes the norm of B. It is assumed that B occupies embedded densely and continuously in H . Let B * denotes the conjugate dual of B. Notice that we have |ϕ| ≤ b ϕ , ϕ ∈ B, for any positive constant b.