1. Introduction
In the last decade, sensor systems have been extensively studied in path planning [
1], environmental monitoring [
2], motor control [
3], and trajectory tracking [
4,
5], and so on. In multi-sensor systems, the accuracy and stability of the system are improved due to the joint data collection by multiple sensors. However, the impact of sensor failures or network attacks in the channel may lead to data transmission time-delay and random packet drops [
6,
7]. Therefore, the investigation of multi-sensor systems is of great importance.
Data processing in multi-sensor systems is performed in the form of fusion, and basic fusion methods include centralized [
8,
9] and distributed [
10,
11]. Centralized is ideally optimal, but when the number of sensors is large, fusion center data processing may be infeasible [
12,
13]. In contrast, the suboptimal distributed structure is more stable. As research goes further, adding an event-triggered transmission strategy to the system can reduce the energy consumption of sensors and decrease the communication burden. Ref. [
14] proposed a distributed event-triggered policy in which the subsystem only broadcasts state information to neighboring nodes when the local state error exceeds a specified threshold. Ref. [
15] proposed a data-driven transmission strategy based minimizing the volume of the non-transmission area. Ref. [
16] proposed a trigger decision based on the estimated variance, where a copy of the Kalman filter is run at the sensor node, and its measurement is transmitted only when the measurement prediction variance exceeds a certain threshold. The event-triggered transmission strategy in [
17] is based on a threshold-based strategy, where the event generator transmits a state measurement only when a signal exceeds a threshold value. A stochastic–deterministic dynamic event-triggered condition is proposed in [
18].
At the same time, the treatment of time-delay problems of systems has received much attention [
19]. The linear matrix inequality (LMI) [
20,
21] and partial differential equation (PDE) [
22,
23] methods are also commonly used in the time-delay treatment of systems. The state augmentation method in [
24] converts time-delay systems into non-time-delay systems with excellent results. The method in [
24] was used in [
25] for a multi-sensor system, but random packet drop was not considered.
State estimation is a pivotal research domain within industrial automation. Consequently, numerous estimation algorithms have been formulated, encompassing the likes of the Kalman filter, Wiener filter, and other notable methodologies. In the system modeling process, modelling errors are inevitable, so the estimator performance must have no sudden changes when the system parameters reasonably deviate from their nominal parameters [
26]. Those with this property are called robust state estimators, and many research methods are available [
27,
28,
29,
30]. A framework based on regularized least squares (RLS) is proposed in [
27], but the modelling errors are restricted to a specific form. A filter that compromises the nominal performance and uncertainty robustness is proposed in [
28]. A robust state estimator based on sensitivity penalty is proposed in [
29], which is not limited to structure-specific modelling errors. In addition, a robust state estimator based on the expectation minimization of estimation error is proposed in [
30]. The study [
31] presents an error estimator, which can be easily implemented in the code. Therefore, it is significant to employ robust state estimators in multi-sensor systems.
In this paper, we investigated the problem of robust fusion estimation for multi-sensor systems with uncertainty, restricted communication, random packet drops, state delay, and deterministic control inputs. A robust state estimator based on state augmentation and sensitivity penalty is used at the local scale. An analytic expression for the robust fusion estimator is derived based on event-triggered, and the pseudo-cross-covariance matrix of the fusion centers is updated. The consistent boundedness of the estimation error is proved. Several simulations verify the effectiveness of the fusion estimator.
The rest of this paper is briefly described below. The problem description and a brief description of the event-triggered transmission strategy are given in
Section 2. A robust fusion estimator for multi-sensor systems with state delays, deterministic control inputs, random packet drops, and communication constraints is derived in
Section 3. The boundedness of the fusion estimator is studied in
Section 4. Several sets of simulations are analyzed in
Section 5.
Section 6 concludes the paper.
2. Problem Formulation and Some Preliminaries
Consider the following discrete-time uncertain linear stochastic system with deterministic inputs and d-steps state delay
where
k represents the discrete-time and
i represents the sensor label. Furthermore,
is the state,
is the measurement,
represents the process noise,
is the deterministic control input, and
is the compound effect of measurement and communication errors. The following assumptions need to be made to guarantee the fitness of the state estimation problem.
(a)
and
are normally distributed with white noise,
,
, and
are mutually independent random variables.
where
,
, and
are known positive definite matrices and
denotes the Kronecker symbolic function.
(b) The elements in the matrices , , , and are known differentiable functions of the modelling errors, and the modelling errors consist of l mutually independent real-valued scalar bounded uncertainties , .
In the process of transmitting the measurement value Y from the sensor node to the fusion center, the channel may experience packet drops. A random variable r is defined to indicate the success or failure of the communication between the sensor node and the fusion center, taking the value of 1 for the successful transmission and 0 when the communication channel fails.
The aim of this paper is to develop a fusion algorithm based on local estimates from each sensor node for multi-sensor systems with parameter uncertainty, state delay, random packet drops, and communication rate limitations. To balance communication cost and estimation performance, an event-triggered transmission strategy like in [
15] is used in this paper.
Consider the following measurement channel
where
is the measurement output,
is the measurement matrix of the system,
represents the state, and
represents the measurement noise. A binary variable is denoted by
t, and when
indicates that the sensor node sends a measurement
Y and the other way around. The specific form of the event-triggered transmission strategy is as follows.
in which
and
are measurable sets. Generally, the center of mass of
is at the origin, that is,
. Note that the decision transmission in the event-triggered transmission strategy is actually when the difference between the measured value and the determined measured value is greater than a threshold value.
The transmission rate, for each sensor node is derived by . In addition, for any given desired transmission rate , a threshold can be easily determined.
Based on Lemma 1 in [
15], a virtual measure
is now defined, where it is uniformly distributed over, and is independent of,
X and
g. Suppose,
,
,
where
. Thus, the optimal transmission strategy is derived as
where
. The random variable
obeys the chi-square distribution with a degree of freedom m where
is the chi-square distribution function with a degree of freedom
m.
Remark 1. The considered multi-sensor system is shown in Figure 1. Each sensor node has state estimation performance with a state delay. Each sensor sends its local state estimate to the fusion center through an unreliable communication channel. All local estimates are fused with data at the fusion center through the best linear unbiased estimation criterion. 3. The Robust Fusion Estimation Procedure
Taking into account the impact of modelling errors on estimation performance, we adopt a robust state estimation algorithm based on sensitivity penalization [
29] to obtain local estimates for multi-sensor systems. A design parameter
,
, is defined to compromise between nominal estimation performance and performance deterioration due to modelling errors. Derived from the foundation of the Kalman filter, this robust state estimation algorithm utilizes sensitivity penalization of model uncertainty estimation errors. It shares a similar form and comparable computational complexity with the standard Kalman filter. When
, this estimator degenerates to the standard Kalman filter.
By introducing the augmentation matrix
X and augmenting the original system (
1) with states, the system becomes
in which,
As can be seen from the above transformation, the re-modeled system is a discrete linear uncertain system without state delay. Following the transformation of the system model from (
1) to (
2), it is evident that the system matrix dimension changes from
n to
.
Remark 2. In this paper, the system is considered only for constant state delays. Based on the state augmentation method, only the system matrix, input (control) matrix, and output matrix of the system need to be changed. The method transforms the original system into a non-time-delay system, but the system dimension will increase from the original n to n(n + d). The state augmentation method is simple and suitable when the delay step is low because the computational burden will increase when it is significant. However, the time delay step is generally manageable in practical production so the problem could be more influential.
To obtain the locally robust state estimate for the
th sensor node, we first define several important matrices
,
, and
, which play a key role in the parameter modification process, as follows:
Let . The detailed realization of the robust state estimation algorithm based on sensitivity penalty is given in Algorithm 1.
Here and are the pseudo-covariance matrices because and .
Based on the event-triggered transmission strategy in the second part, whether each sensor node sends a local state estimate to the fusion center is determined by
. The transmission strategy mentioned above can be expressed as
In order to guarantee the transmission rate , the vector , the positive definite weight coefficient matrix , and the positive real numbers must be chosen appropriately. is the local state estimate.
Notice that each local state estimate can be interpreted as a measurement
of the true state
collected through the virtual measurement channel defined as
where the estimation error
can be regarded as virtual measurement noise.
Algorithm 1: The local robust state estimation based on sensitivity penalty (Appendix A) |
|
Now, considering only the event-triggered transmission strategy, (
4) corresponds to the measurements received by the fusion center from sensor node i, that is,
. When sensor data are not transmitted, (
4) will be replaced by
Here,
is uniformly distributed within the ellipsoid mentioned in (
3) and is not correlated with the estimation error
.
According to the event-triggered transmission strategy, when there are packet drops in the communication channel from the estimator to the fusion center, the virtual measurement channel can be replaced with
where
is explicitly utilized in (
6) to indicate whether packet drop occurs in sensor transmission to the fusion center and
. The state of the multi-sensor system is shown in
Table 1. For simplicity, the event-triggered is abbreviated as ET, and the success of the transmission is simplified as PD. The
in (
6) is the predicted values of the
th sensor node for moment
k.
is the virtual measurement noise of the
th virtual channel for moment
k, which can be derived by
The fusion estimation with both random packet drops and event-triggered transmission strategies is investigated, and the following matrices are defined as
The information in the fusion center is obtained from the virtual measurement channel
In accordance with the best linear unbiased criterion (BLUE) in [
32], we can obtain the fusion estimate and its error covariance matrix.
In (
9),
is the covariance matrix of the virtual measurement noise, which is the global error covariance matrix of the estimation error. From
in (
8), the expression of
can be obtained as
in which
. The matrices
,
, and
in the formula are equal to
Then, we consider the state estimation errors of the following dynamic system.
The following relationships can be easily obtained
According to the above equation, the explicit expressions for the three pseudo mutual covariance matrices
,
, and
in (
11) can be derived as follows
in which
,
.
is a pseudo-covariance matrix in robust state estimation. Thus, there is
.
4. Some Properties of the Fusion Estimator
This section has the goal of investigating the steady-state properties of event-triggered robust fusion estimators for multi-sensor systems with deterministic inputs, random packet drops, and state delays. Assume that the modelling errors in this section are within the set , . The matrices , , and are denoted as , , and , respectively. In addition, the following assumptions need to be made.
(A) , , , , , , , , , , and are time-invariant.
(B) The uncertain linear system of (
1) is exponentially stable in the sense of Lyapunov and the matrices
,
,
,
,
,
,
,
are bounded whenever
and
.
(C) For every sensor node,
is detectable and the following matrix pair is detectable
where
.
Theorem 1 ([
15])
. Suppose that Assumptions (A), (B), and (C) hold and that each sensor transmits local estimate according to the event-triggered transmission strategy. If the weight matrix of the sensor node satisfies the conditionfor some positive real number , the estimation error is consistently bounded for any possible choice of , which means Proof of Theorem 1. Let
be the estimate obtained at time
k through
instead of
,
, which gives
so we have
Taking into account the first term on the right-hand side in (
15), since
is based on the vector
, the following inequality can be obtained
According to Assumptions (A), (B), and (C), then
is convergent, and
and
are also convergent [
33]. The estimation error has a bounded covariance matrix at each
k. This indicates that
is converged, and the estimation error covariance matrix is bounded.
From the inequality condition in Theorem 1 and the remainder of
, we can obtain
Hence, the uniform boundedness of
can be obtained by (
16). Now it is only necessary to prove that the second part of the right-hand side of inequality (
15) is uniform boundedness. Under the inequality condition in Theorem 1, it can be obtained as
When , it means that there is . Furthermore, it can be obtained that , then . The proof is done. □
To minimize the volume of the non-transported region,
and
can be appropriately denoted as
in which
. Two methods exist for determining the local prediction of
as per (
19). The first method utilized in this paper is a local estimation based on sensor nodes. This fusion estimation method does not necessitate broadcasting but requires each sensor node to retain past information. The second method is based on the
moment fusion estimation
.
5. Numerical Simulations
This section cites the tractor–car system detailed in [
34], shown in
Figure 2, and extends it to a multi-sensor system for sample simulations. The performance of the derived robust fusion estimator is demonstrated through comparison with the fusion estimator for the Kalman filter based on actual and nominal parameters using the same fusion method across two distinct sets of numerical simulations with modelling errors (fixed or not) and varying transmission rates and packet drop rates. This numerical simulation consists of two sensors. For each set, 500 time experiments were conducted, with 200 moments designated for each set, generating 200 input–output data pairs. In the simulations, the overall average estimated error variance
is computed for each moment, and the implementation of event-triggered and occurrence of packet drops are displayed.
Since the vehicle steering and directional angles in the tractor-car system are nonlinear, they can be linearized and expressed as
in which
,
,
,
,
, and
are the direction angle of the tractor, the direction angle of the car, the tractor steering angle, the process noise, d-step time-delay for state 1, and d-step time-delay for state 2, respectively.
is the state vector,
.
L,
k, and
v denote the length of the tractor, the sampling period, and the constant speed, respectively. Considering the system errors at linearization in the form of modelling errors
substituted into the system model, the matrix parameters are obtained as
In the numerical simulation, each parameter is taken as
,
s, and
, and a two-step state delay system was used. The matrix parameters are as follows
The packet drop process is assumed to be a stationary Bernoulli process. A constant value of is assigned to the filter design parameter .
In Case 1, the modeling errors
are assumed to be a fixed value of −0.8508. The transmission and packet drop rates for both sensors are set to 0.8 and 0.2, respectively.
Figure 3a illustrates the fusion estimation error over time, demonstrating that the robust fusion estimator proposed in this study outperforms the fusion estimator for the Kalman filter based on nominal parameters by approximately 7.800 dB.
Figure 3b,c depict the transmission of the two sensors and the packet drops of the communication channel, respectively. To clearly reflect the execution of the event-triggered,
is inverted, and
is treated similarly. Note that the plots of event-triggered realizations and packet drops here are from one of the 500 experiments used.
The modelling errors
are generated randomly and independently, conforming to a normal distribution with a truncation. The mean, standard variance, and truncation values of the normal distribution are set to
,
, and
, respectively.
Figure 4a illustrates that the derived estimator surpasses the performance of the Kalman filter based on nominal parameters, and it can be seen from the 200th moment that the estimator derived in this paper is 5.8600 (dB) lower than the nominal parameter-based Kalman filter.
Figure 4b,c show the realization of the sensor transmission and the channel packet drop over 200 moments, respectively.
In Case 2, the derived robust fusion estimator is tested using different transmission and packet drop rates. The modelling errors are the same as in Case 1 with a truncated normal distribution. Based on the analysis of
Figure 5, it is evident that the derived estimator exhibits effective and reliable operation even under diverse transmission rates generated by the employed event-triggered transmission strategy. However, variations in transmission rates give rise to disparities in estimation performance, a well-studied phenomenon. This can be attributed to the fact that higher transmission rates are associated with improved estimator performance. As the transmission rate increases, the fusion center receives a greater volume of estimation values, thereby leading to more accurate results. A reasonable analysis of
Figure 6 demonstrates that the derived estimator effectively maintains its reliability even under diverse packet drop rates. Nonetheless, differing packet drop rates introduce disparities in estimation performance, which is a valid observation. Higher packet drop rates correspond to inferior estimation performance. When compared to
Figure 5, it is apparent that the variation in estimation performance is greater for different packet drop rates than for different transmission rates.
As can be seen from the two sets of simulations, the proposed robust fusion estimator exhibits relatively better performance compared to the fusion estimator that ignores uncertainty. The derived robust fusion estimator is still applicable when the selection of modelling errors is not limited to the particular structure. The results show that the method is an effective multi-sensor fusion method in practical applications.