Abstract

This paper is concerned with the problem of distributed estimation fusion over peer-to-peer asynchronous sensor networks with random packet dropouts. A distributed asynchronous fusion algorithm is proposed via the covariance intersection method. First, local estimator is developed in an optimal batch fashion by constructing augmented measurement equations. Then the fusion estimator is designed to fuse local estimates in the neighborhood. Both local estimator and fusion estimator are developed by taking into account the random packet losses. The presented estimation method improves local estimates and reduces the estimate disagreement. Simulation results validate the effectiveness of the proposed distributed asynchronous fusion algorithm.

1. Introduction

Distributed fusion and estimation problem over peer-to-peer sensor networks has attracted significant interest in the research community, because of its variety of applications, such as environmental monitoring, surveillance, and target tracking, [110]. In these applications, every sensor in the network does not only take measurements in a parallel manner but also acquires information from neighbors and processes it to get an estimate. Compared with the centralized fusion fashion, the main advantage of distributed fusion estimation is the computation burden alleviation and robustness enhancement.

In general, two main approaches are presented in the literature to solve the distributed fusion and estimation problem. The first approach is the consensus approach. The consensus approach was proposed in [11, 12], where local measurements are exchanged among neighbors to get local estimates, and then by using average consensus algorithms among neighbors every sensor in the network gets the same estimate in steady-state. In order to obtain the same estimate at every sensor in the network, the consensus approach may iterate several times for each new measurement. This is highly undesired when estimating the state of dynamic systems where new measurements need to be processed in a timely manner. The second approach is the diffusion approach. The diffusion approach was presented in [13], in which the estimates of local filters are calculated individually at each sensor by using the data from the neighborhood, and then the local estimates from the neighborhood are fused locally by a convex combination. Therefore, the diffusion approach is well suited for estimating the state of dynamic systems where new measurements are being taken in real time. Though the adaptive weights for the diffusion algorithm were presented in [14], the estimation error covariance information is not taken into account. In [15, 16], the covariance intersection (CI) algorithm was proposed, in which the coefficients of convex combination of the estimates are chosen by considering the error covariance information. Therefore, the CI algorithm was used to design the adaptive weights for the diffusion algorithm in [17]. Though the aforementioned works solved the distributed fusion and estimation problem in their specific perspectives, they did not take into consideration the packet dropouts which are usually unavoidable in networks.

Many results on control and state estimation in networked systems context are already presented for the problem of both unreliable communication links and measurement uncertainties. The interested readers may refer to [1823] and the references therein for further information. We will review only those works that are closely related to the state estimation in sensor networks with packet dropouts.

The problem of packet dropouts which are usually unavoidable in sensor networks has gained a lot of attention [2427]. Schenato addressed a minimum error covariance estimator for sensor networks with random packets losses in [28]. The critical data arrival rate which ensures the stability of Kalman filter is investigated in [29]. An optimal filter is proposed based on the probability of the packet dropouts in [30]. Sun proposed optimal linear estimators in the sense of minimum variance with respect to the packet dropout probability in [31]. Ma and Sun proposed a centralized linear optimal fusion estimator in the minimum variance sense in [32]. Though the aforementioned literatures deal with the random packet losses problem, the sensor asynchronous problem is not taken into account.

Most of the aforementioned works are developed under the synchronous assumption. And the real-world systems may usually be asynchronous. In regard to asynchronous sensor fusion, a series of linear weighted fusion (LWF) algorithms for two and more than two asynchronous sensors with and without feedback had been proposed separately in [3336]. By establishing state-space models at each sampling rate, a new fusion algorithm for asynchronous sensors had been presented in [37]. A minimum mean-squared error (MMSE) centralized asynchronous fusion algorithm had been presented in [38], which is suitable for arbitrary sensor number, sensor sampling rates, sensor initial sampling time instants, and fusion periods. By reconstructing the optimal centralized asynchronous fusion algorithm in [38], a so-called distributed asynchronous fusion algorithm is proposed in [39], which performs equivalently to the optimal centralized asynchronous fusion algorithm. Though the aforementioned works deal with the asynchronous sensor fusion problem, the random packet dropouts problem is not taken into account. Moreover, these works always need a fusion center to accomplish the estimation task.

Motivated by the above discussion, we consider the distributed estimation fusion problem for an asynchronous sensor network with random packet dropouts. The sensor network is considered to be a peer-to-peer network with no fusion center, and the communication channels of the sensor network are unreliable. The main contribution of our paper is the development of the distributed asynchronous fusion algorithm in a peer-to-peer sensor network with random packet dropouts. Our algorithm is developed as follows. Firstly, a Bernoulli random variable is employed to describe the random packet dropouts. Then, by constructing a local augmented measurement with the random packet dropouts under consideration, an asynchronous local estimator is obtained at each sensor to compute local estimate. Further, the fusion estimators at different sensors are presented by utilizing a weighted fusion criterion based on covariance intersection method, and the error covariance information of local estimates is used to compute the fusion weights. The main advantage of the proposed distributed asynchronous fusion algorithm is the practicality when dealing with a practical asynchronous sensor network with unreliable communication links.

The rest of the paper is organized as follows. In Section 2, the asynchronous distributed estimation problem is formulated. Then, design procedures for the asynchronous distributed fusion algorithm are presented in Section 3. The effectiveness of the proposed algorithm is demonstrated by a simulation example in Section 4. Finally, the conclusion is provided in Section 5.

2. Problem Formulation

Consider a continuous-time dynamic system described by the following linear stochastic differential equation: where and is zero-mean white Gaussian noise with covariance matrix . Suppose that is the corresponding state transition matrix of , is the time when the fusion center executed the last fusion operation or the initial time if , is the oncoming fusion time instant, and is the fusion interval. Then, by discretizing the continuous-time linear system (1), we have where is zero-mean white Gaussian noise with covariance matrix .

Observations of system (1) are obtained by a sensor network which deployed spatially. The sensor network which consists of distributed sensors is considered to be a peer-to-peer network. The interconnection topology of the sensor network is denoted by . indicates that sensor and sensor are connected, otherwise implies that sensor and sensor are not connected. No fusion center exists in the sensor network, and every sensor acts also as an estimator which collects information only from its neighbors to generate local estimates. The set of sensors connected to the certain sensor is called the neighborhood of sensor and it is denoted by . Without loss of generality, assume that and the number of neighbors of sensor is given by the number of elements of , written as .

Suppose that the number of measurements provided by sensor during the th time interval is . So observations during the th time interval of system (1) can be modeled as where and are measurement matrix and zero-mean white Gaussian measurement noise with covariance matrix , respectively. corresponds to the th measurement of sensor during time interval . The set of measurements obtained by sensor during time interval is denoted by . Assume that , , and are mutually independent.

We also assume that the sensors in the network are time synchronized and time driven; therefore the sensors can easily obtain the time stamp through the digital communication network. The observations are exchanged between neighboring sensors once in every fusion interval. The packet-dropping process at every communication links is modeled as independent identically distributed (i.i.d.) Bernoulli random process. We use Bernoulli random variable (for all ) to indicate whether a packet is successfully transmitted from sensor to sensor at th time instant. If , the packet is successfully transmitted, otherwise, and the packet is dropped. Without loss of generality, it is assumed that .

The th measurement sensor got from sensor during time interval can be written as where is the packet arriving probability (PAP) and is the packet loss probability (PLP) [17].

3. Distributed Asynchronous Fusion Algorithm

3.1. Local Asynchronous Estimator

Without loss of generality, it is assumed that sensors in totally obtain a set of asynchronous measurements during the th time interval. Measurement is provided by sensor at time instant . After sensor collects all asynchronous measurements which are provided by sensors in , it reorders these measurements in time sequence. Suppose the th measurement at sensor during the th time interval which is denoted by corresponds to , then we have

As mentioned before, we can obtain the time stamp of current measurement easily. The packet is successfully transmitted, if the time stamp of the latest measurement is in the present fusion interval, that is, , otherwise . We can simplify the measurement as where , , and .

The state can be written by using (2) as where is the backward state transition matrix.

By using (6) and (7), can be rewritten as where , , and .

By denoting where , are identity matrices with proper dimensions.

We can get an augmented measurement equation at sensor as

From (2), (3), (6), and (8) we have

Suppose ; then Thus, the augmented measurement noise is zero-mean white Gaussian noise with covariance:where

From (2), (8), and (11), we have Further we have where .

It can be observed from (19) that the discrete-time process noise and the augmented measurement noise are cross correlated; this correlation comes from the common continuous-time process noise [40].

For sensor , the asynchronous local estimation problem in is converted to an estimation problem of the following system: where , , and .

Thus a local asynchronous estimation method is summarized in Theorem 1 in an optimal batch asynchronous data fusion fashion [38].

Theorem 1 (local asynchronous estimator). For sensor , after it gets asynchronous measurements from sensors in , the local asynchronous estimator at sensor is summarized as follow: given the local state estimate of system (1) and its error covariance matrix at time , the local estimate and its error covariance matrix at time are given by where and are the state estimate and one-step predicted estimate, respectively, and are the covariance matrices of filtering error and one-step prediction error, respectively, and the initial values of and are given, respectively, by and , where is a normally distributed random variable with mean equal to the initial truth ground of the state and covariance .

Proof. For neighborhood , the local asynchronous estimation problem at sensor is converted to the state estimation problem of system (20). This is a classical filtering problem with noise cross-covariance which can be solved by using the Kalman filter [40] for correlated process and measurement noises.

Remark 2. The random packets losses have been taken into account in Theorem 1. And in Theorem 1 each sensor in the sensor network generates a local estimate at the fusion time instant by using measurements only from its neighbors. These local estimates are suboptimal in the sense that not all the measurements in the sensor network are used. Therefore, a fusion criterion, which will improve each local estimation performance, weighted by matrices will be developed later in this paper. The main results will be presented in the following subsection.

3.2. CI-Based Fusion Estimator

In this subsection, a fusion criterion weighted by matrices is developed to generate fusion estimate at every sensor, and it takes into consideration the random packets dropouts. The criterion is first given in the following lemma in [16].

Lemma 3. Let denote unbiased estimates for the unknown state vector and let denote the corresponding estimated covariance matrices; the CI fusion estimator is given by the following convex combination: where and hold.

When local estimates calculated by the local estimators in Theorem 1 are available at the sensors in the sensor network, sensor collects them from its neighborhood to generate a fused estimate according to the fusion rule based on covariance intersection. For the packet-dropping process, it is also modeled as i.i.d. Bernoulli random process. Denote (for all ) as the packet transmit indicator. It is also assumed that .

Theorem 4 (CI-based fusion estimator). By gathering the local estimates provided by the sensors in neighborhood , the distributed weighted fusion filter at sensor is given by where are subject to , .

Proof. Theorem 4 follows directly from Lemma 3.

Remark 5. Different choices of can be used to optimize the update with respect to different performance criteria, such as minimizing the trace or determinant of . However such optimization is nonlinear and a high computation load is required for computing the optimal [41]. Thus, several fast CI algorithms that produce suboptimal solutions have been proposed in terms of trace or determinant minimization [41, 42]. In what follows, we choose to minimize the trace of . For sensor , a suboptimal noniterative weighting coefficients solution for fast CI fusion is given by

For clarity, the asynchronous distributed fusion algorithm is summarized in Algorithm 1.

(1) Initialization: , .
(2) Obtain local measurements .
(3) Broadcast message to neighbors.
(4) Receive messages , , from neighbors.
(5) Compute , , , and by using (9)–(19).
(6) Calculate the local estimate by using Theorem 1.
(7) Broadcast message to neighbors.
(8) Receive messages , , from neighbors.
(9) Calculate the fusion estimate by using Theorem 4 and (25).

Remark 6. In the distributed asynchronous fusion filter presented in Algorithm 1, each sensor in sensor network broadcasts its local measurement information and local estimate to its neighboring sensors; for every sensor, more information from different sensors, not only the neighboring ones, will be utilized to generate the fusion estimates, which helps improve the local estimate performance and reduce the disagreement of local estimates. Further, not all target information in the sensor network is used to generate the fusion estimate at one sensor; meanwhile the method we used to compute coefficient is suboptimal; therefore the fusion estimate is a suboptimal estimate, but it is an improved one compared with the local one.

4. Simulation

In this section, the effectiveness of the proposed estimator design method is demonstrated by presenting simulations of a target tracking system. Consider the following continuous-time dynamic system: where and is zero-mean white Gaussian noise with variance ; by discretizing the continuous-time linear system (26), we have where the process noise is a white Gaussian noise sequence with zero mean and covariance

The state of system (27) is observed by a sensor network with 12 sensor nodes, and the topology of the sensor network is shown in Figure 3. The observation equations of the sensors are given by (3) with measurement matrices , and ; measurement noise variances , and ; and sampling periods , , and , respectively. The fusion operation is performed with period (i.e., ) and initial time instant .

For performance comparison, root mean square errors (RMSE) of position and velocity given by [40] are used, where is the Monte-Carlo simulation times, is the ground truth of the th component of the state vector in the th run, and is the corresponding estimate; the disagreement potential of the estimates in the neighborhood is defined as where is the estimate of sensor at time instant , is the corresponding error covariance, and the simulation ends after 20 seconds.

We present the simulation results without packet losses firstly. As shown in Figure 1, sensors 2 and 12 are directly connected to sensor 4, thus neighborhood consists of sensors 2, 4, and 12, with sampling periods and , respectively. We take estimate performance of sensors in neighborhood to testify the effectiveness of the proposed algorithm. At each fusion instant , local estimate at sensor 4 is generated by collecting measurements from itself and sensors 2 and 12, and then sensor 4 collects local estimates from itself and sensors 2 and 12 to form fused estimates by using Theorem 4 and (25). Since trends of the simulation results are the same for both components, only the first component of the state vector is shown and analyzed below, and the simulation results of the proposed algorithm are compared with the optimal linear minimum mean-squared error centralized asynchronous fusion algorithm in [22].

Figure 2 compares estimate performance of the proposed local asynchronous estimator and fusion estimator with the centralized asynchronous fusion estimator proposed in [22]. Although the performance of local estimator and fusion estimator is close to the centralized one, the centralized asynchronous fusion algorithm outperforms them both. This verifies our conclusion that the proposed local asynchronous estimator and fusion estimator are suboptimal as we stated in Remarks 2 and 6.

Figure 3 compares the performance of the proposed local asynchronous estimator and fusion estimator. It can be observed that the RMSE of fusion estimates is smaller than that of the corresponding local estimates; this verifies that the fusion estimates improve the performance of local estimates.

Figure 4 compares the disagreement of estimates obtained by the proposed local asynchronous estimator and fusion estimator; the fusion estimator gains smaller disagreement potential than the corresponding local one. This verifies our conclusion that fusion estimator reduces the disagreement among different sensors.

Figures 5 and 6 shows the filtering performances of the sensors in with different PLPs. It can be observed that packet dropout degrades estimation performance. The fusion estimator performance still outperforms local estimator and both estimators RMSEs converge. This demonstrates the effectiveness of the proposed distributed asynchronous fusion algorithm.

5. Conclusion

A CI-based distributed asynchronous fusion algorithm has been proposed for decentralized asynchronous sensor networks with random packet dropouts in this paper. The proposed distributed asynchronous fusion algorithm is applicable to decentralized asynchronous sensor networks; it is capable of improving estimates and reducing disagreement in comparison with the local ones; meanwhile it is also able to deal with random packet losses. The simulation results demonstrate the effectiveness of the proposed distributed asynchronous fusion algorithm.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 61104186 and no. 61273076) and the Natural Science Foundation of Jiangsu Province (no. BK2012801).