# Energy-efficient target tracking in wireless sensor networks: a quantized measurement fusion framework.

1. IntroductionTarget tracking through wireless sensor networks (WSNs) is a problem with a large spectrum of applications [1, 2], such as surveillance [3], natural disaster relief [4], traffic monitoring [5], and pursuit evasion games. Although it can have several advantages, target tracking in a WSN is a challenging task since each sensor node typically has very limited power supply and communication bandwidth [6, 7]. Besides, target tracking in this circumstance needs collaborative communication and computation among multiple sensors since information generated by a single sensor node is usually incomplete or inaccurate (e.g., [8]). Collaborative data processing and in-network processing have been extensively studied to reduce redundant communication in target tracking applications. Clustering or group-based techniques such as information-driven sensor query (IDSQ) [9] and dynamic convoy tree-based collaboration (DCTC) [10] are well-known solutions for collaborative target tracking. These schemes explicitly select and manage sensors that participate in target tracking.

Moreover, as mentioned above, each sensor node in the network has very limited energy/power source and communication bandwidth. Thus, by using quantized message for storage or transmission in the place of the original measurement, considerable savings in storage or transmission bandwidth can be realized. The fusion center (FC) will combine the quantized message from local sensors to produce a final estimate of the parameter. The problem of decentralized estimation and tracking based on quantized measurements has been studied in earlier works such as [11,12]. Recently, a distributed estimation approach based on the sign of innovations (SOI) is developed in [13], where only the transmission of a single bit per observation is required. In [14], particle filtering algorithms for tracking a single target using data from binary sensors are presented with the posterior Cramer-Rao lower bound (CRLB) derived. The posterior CRLBs on the mean squared error (MSE) of target tracking in WSNs with quantized range-only measurements are discussed in [15]. Both constant velocity (CV) and constant acceleration (CA) models for target dynamics and general range-only measuring model for local sensors are included, while the sensor nodes are assumed to be randomly distributed. Very recently, the problem of target tracking in a WSN with quantized measurements using adaptive thresholds is considered in [16]. Attention is focused on the design of measurement quantizer with adaptive thresholds. Based on the probability density function (PDF) of the signal amplitude measured at a random location, an adaptive design method for quantization thresholds is proposed by maximizing the entropy. However, there is little research on collaborative target tracking with nonlinear dynamics and/or observing equation from the quantized measurement fusion perspective.

In this paper, we consider the quantized measurement fusion problem for collaborative target tracking with focus on the energy and bandwidth scheduling strategy. Consequently, a complete framework of quantized measurement fusion for target tracking in a WSN is proposed in this paper. To be more specific, the following problems are considered.

(1) After receiving the quantized measurements from local node, the FC fuses the quantized measurements, for which we consider two fusion approaches. One is the augmented approach, which merges all the noises corrupted quantized messages into a vector. The other is the weighted approach under quasi-best linear unbiased estimation (quasi-BLUE) fusion rule. Note that the latter approach has a lower computational load since the observation vector dimension kept unchanged regardless of the number of the sensors deployed [17]. This is highly desirable in densely deployed WSN since the computational burden will keep low regardless of the number of the sensors activated.

(2) With the focuses on tradeoff between energy/power consumption and fusion accuracy for quantized measurements, we give closed-form solution to the optimization problem for bandwidth scheduling for both approaches. In the optimal bandwidth scheduling problem, the total energy consumption is minimized subject to a constraint on quantization covariance or the overall MSE incurred by weighted fusion.

(3) Due to the fact that generally passive sensors are employed in the sensor networks, we consider the nonlinear Gaussian discrete-time system model. The sequential importance resampling (SIR) particle filtering strategy is used for estimating the target state. The posterior CRLBs (see, e.g., [18, 19] and the references therein) for target tracking using quantized measurements are also derived.

The rest of the paper is organized as follows. In Section 2, the measurement fusion problem using quantized measurements will be formulated, and the quantization strategy is introduced. In Section 3, two quantized measurement fusion approaches are proposed, namely, the augmented approach and the weighted approach. The bandwidth scheduling problems for both approaches are explored, respectively, in Section 4 and Section 5. The tracking algorithm and the performance analysis are provided in Section 6. A numerical example is given in Section 7, which is followed by Section 8 that involves the conclusion and future work.

2. Problem Formulation

Considering the state estimation problem of the following nonlinear discrete-time system in a network with N sensors deployed

x(k+1) = f(x(k),k) + a(k) (1)

y, (k) = [h.sub.i] (x(k),k) + [v.sub.i] (k), i = 1,2,3,... ,N, (2)

where x(k) [member of] [R.sup.n] and [y.sub.i](k) [member of] [R.sup.m] are the state of the system and the ith measurement at the time step k, respectively; f and h; are known, possibly nonlinear functions of the state x(k) and time step k. [omega](k) [member of] [R.sup.p] is the noise process caused by disturbances and modeling errors; [v.sub.i](k) [member of] [R.sup.m] are additive measurement noise vectors of the ith (i = 1, 2, 3,... ,N) sensor. It is assumed that the noise vectors w(k) and [v.sub.i](k) are zero mean and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)

The initial state x(0) with mean [x.sub.0] and variance [P.sub.0] is independent of [omega](k) and [v.sub.i](k), i = 1,2,3,... ,N. To ease the analysis, we also assume that all sensors are synchronised and have the same measurement rate.

Due to bandwidth and energy limitations, each sensor quantizes its measurement into a [b.sub.i](k)--bit message (where [b.sub.i](k) is the bandwidth to be determined later) and transmits this locally processed data to a FC. Then the FC estimates the state vector of the target according to (1) and the quantized messages. The purpose is to collaboratively estimate the target state using quantized measurements fusion strategy while meeting constraints on total energy consumption. More specifically, we consider the following aspects: (i) how to fuse the quantized messages; (ii) how to optimize bandwidth scheduling among sensors; and (iii) how to estimate the target state using these quantized messages. Before going further, we first describe the quantization scheme used in this paper.

2.1. Probabilistic Quantization Strategy. In the following, we present a brief review of the quantization scheme in this paper. Suppose that the measurement [y.sub.i](k) is bounded to a finite interval [-W,W], where W is a known parameter decided by the sensors' dynamic range. Due to bandwidth and power limitations, each sensor quantizes its observation into a [b.sub.i]--bit message, where bt is the bandwidth to be determined later. We, therefore, have [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] quantization points given by the set [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The points are uniformly spaced such that [DELTA] = [[tau].sub.j+1] - [[tau].sub.j] for j [member of] [1,2,...,[L.sub.i] - 1}. It follows that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Now suppose [y.sub.i](k) [member of] [[[tau].sub.j], [[tau].sub.j+1]) and let [z.sub.i](k) = p([y.sub.i](k)), where p(x) denotes the probabilistic operation. Then [y.sub.i](k) is quantized in a probabilistic manner:

Pr {[y.sub.i] (k) = [[tau].sub.j+1]} = r, Pr {[y.sub.i](k) = [[tau].sub.j]} = 1 - r, (4)

where r = ([y.sub.i](k) - [[tau].sub.j])/[DELTA]. The following lemma, adopted from [16], discusses two important properties of probabilistic quantization.

Lemma 1. Suppose interval [-W, W] bounded measurement [Y.sub.i](k) [member of] [[[tau].sub.j], [[tau].sub.j+1]), and let [z.sub.i](k) be a [b.sub.i]--bit quantization of [y.sub.i](k). The message [z.sub.i](k) is an unbiased representation of [y.sub.i](k); that is,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (5)

Assuming that the sensor messages are perfectly received by the FC, then the quantized message at the ith sensor can be modelled as

[z.sub.i] (k) = [y.sub.i] (k) + [q.sub.i], i = 1,2,3,... ,N, (6)

where qt is the zero mean quantization error with a variance less than [DELTA]2/4. Moreover, as performed locally at each sensor without coordination, the quantized noises atkth time iteration are mutually independent.

3. Two Quantized Measurement Fusion Approaches

By Lemma 1, the relationship between the original signal of ith sensor and the data received by the FC can be reformulated by

[z.sub.i](k) = [h.sub.i] (x(k),k) + [n.sub.i] (k), i = 1,2,3,... ,N, (7)

where [n.sub.i](k) = [v.sub.i] (k) + [q.sub.i](k) are mutually independent with zero mean and variance:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (8)

By quantized measurement fusion, the N local sensor observation models in (7) can be integrated into the following vector model:

z(k) = h(x(k) ,k) + n(k), (9)

where n(k) are zero mean fused noise with variance R(k) that will be determined as follows.

Commonly, there exist two approaches to obtain the quantized measurement fusion (see, e.g., [17, 20] and the references therein). One is the augmented approach that merges all the received data into a vector form as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

Then by approximately considering the noise [n.sub.i](k) to be white noises, we can use nonlinear filtering techniques to get the state vector estimate (see, e.g., [21]).

The other approach first combines the quantized measurements based on the BLUE rule (we will refer to weighted approach thereafter), and then traditional filtering schemes are employed to estimate the state of the target. Specifically, the fusion measurement can be obtained by weighted observation using the quasi-BLUE scheme:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (11)

The incurred MSE is thus

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (12)

Remark 2. Because the first approach makes use of all the raw quantized messages, it most likely provides a better estimate. However, the latter approach has a lower computational load since the observation vector dimension is kept unchanged regardless of the number of the sensors deployed [17]. This is especially preferable in densely deployed wireless sensor network since the number of sensors deployed in a snapshot will be large, which makes the dimension of observation space high. The cost of calculating inverse matrix, which is involved in computing the gain of the filter, is in proportion to cubic dimension of system observation space; therefore, the computational burden will increase when using the augmented approach as (10).

Remark 3. It is easy to see that the accuracy of fused quantized measurement is better if the variance of the quantized noise is smaller. We can make its upper bound small, which means that more bandwidths need to be supplied. However, the sensor power and transmission bandwidth are limited in a WSN. In the following, we will consider the bandwidth scheduling problem that solves the tradeoff between energy/power constraints and fusion performance for both quantized measurement fusion approaches.

4. Bandwidth Scheduling for Augmented Approach

In this section, we consider the bandwidth scheduling problem for the quantized measurement fusion via augmented approach. First, a practical energy model is introduced. Then, based on the energy model and the quantization scheme, we set up the optimization problem for bandwidth scheduling. Finally the explicit optimal solution is obtained through the Lagrange technique.

4.1. The Energy Model. We assume that the channel between the ith sensor and the FC experiences a path loss proportional to [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].where [d.sub.i] is the transmission distance between the ith sensor and the FC. Then the consumed energy of the ith sensor at time step k is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is the energy density [22], in which [rho] is a constant depending on the noise profile, and [P.sub.b] is the target bit error rate assumed common to all sensor-to-FC links.

4.2. Optimization Problem Setup. Our primary goal is to minimize the total energy consumption while meeting the constraints on fusion performance. Thus, we consider the following optimization:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (14)

where D(k) > 0 is a given performance constraint. Without illegibility, we will drop out the time step index k thereafter.

Let [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Then the optimization problem in (14) can be equivalently rewritten as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

The lower limit on Bt is given by the fact that [b.sub.i] [member of] [Z.sup.+.sub.0] and the definition of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Furthermore, we relax the integer Bt to be a real positive number so as to render the problem tractable and drop the last constraint on Bt as they become inactive at the optimum point. The major advantage of the alternative problem formulation (15) is that it admits the form of convex optimization and leads to a closed-form solution as shown next.

4.3. Optimal Solution. By leveraging the standard Lagrange technique, the optimal solution to (15) can be obtained as follows [23]. Without loss of generality, we first assume that [[beta].sub.1] [greater than or equal to] [[beta].sub.2] [greater than or equal to] ... [greater than or equal to] [[beta].sub.N] and define

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

Let 1 [less than or equal to] [K.sub.1] [less than or equal to] N be the unique integer such that f([K.sub.1]) < 1 and f([K.sub.1] - 1) [greater than or equal to] 1. Following a similar way as in [16], it is ready to show that such [K.sub.1] is unique unless f(K) < 1 for all 1 [less than or equal to] K [less than or equal to] N, in which case we simply set [K.sub.1] = 1. Then simple manipulations of (16) lead to the optimal solution given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (17)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (18)

Once the optimal real-value [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is computed, the associated bit loads can be readily obtained through upper integer rounding, as in [16, 22].

Remark 4. Recall from (12) that the energy consumption of each sensor is proportional to the path loss gain [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Hence, large value of the energy consumption corresponds to sensor deployed far away from the FC, usually with poor background channel gain. The most important observation to make from the optimal solution (18) is that the assigned bandwidth is inversely proportional to the logarithm of the energy consumption of this node. This is intuitively reasonable since sensors with better link conditions should be allocated with more bits to improve the estimation fusion accuracy.

5. Bandwidth Scheduling for Weighted Approach

By the same energy model adopted in Section 4, we consider the following optimization problem according to the MSE recurred by BLUE fusion in (11)-(12):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

where D(k) > 0 is a given performance constraint.

5.1. Alternative Formulation. To facilitate the analysis, we also relax the integer [b.sub.i](k) to be a real positive number and drop the last constraint on [b.sub.i](k) as they become inactive at the optimum point. However, even with this relaxation, the problem (19) remains nonconvex in [b.sub.i](k). Fortunately, we show next that the problem (19) can be reformulated as a convex one by performing a change of variables. Define = [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Then, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)

Therefore, the problem (19) can be transformed into a problem with variable r = ([r.sub.1], [r.sub.2],..., [r.sub.N]) shown as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (21)

The upper limit on r{ is given by the fact that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. It is easy to check that the optimization problem (21) is convex with respect to the new variable r = ([r.sub.1], [r.sub.2],...,[r.sub.N]).

5.2. Optimal Solution. By leveraging the standard Lagrange technique, we can derive the optimal solution as follows. Without loss of generality, we assume that [[beta].sub.1] [greater than or equal to] [[beta].sub.2] [greater than or equal to] ... [greater than or equal to] [[beta].sub.N] and define

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)

Let 1 [less than or equal to] [K.sub.1] [less than or equal to] N be the unique integer such that f([K.sub.1]) < 1 and f([K.sub.1] - 1) [greater than or equal to] 1; if f(K) < 1, for all 1 [less than or equal to] K [less than or equal to] N, then simply set [K.sub.1] = 1. Then simple manipulations lead to the optimal solution given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (23)

where [(x).sup.+] equals 0 when x < 0 and otherwise is equal to x. Recall the definition of [r.sub.i]; we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (24)

Thus, (23) implies that the optimal solution of [b.sub.t] is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)

Once the optimal real-value [b.sup.*.sub.i] is computed, the associated bit loads can be obtained through upper integer rounding. As mentioned above, large value of the energy consumption corresponds to sensor deployed far away from the FC, usually with poor background channel gain. In this point, the optimal solution (25) is similar to (18): for those active nodes, the assigned bandwidth is inversely proportional to the energy consumption of this node. In other words, to improve the fusion estimation accuracy, sensors with better link conditions should be allocated with more bits. Using the dynamic equation (1) and fused measurement (9), the state of the target can be estimated in the FC through nonlinear filtering technique, such as extended Kalman filter, unscented Kalman filter, and particle filter.

6. Particle Filtering and Posterior CRLBs for Tracking with Quantized Measurements

Bayesian sequential estimation, also known as Bayesian filtering, is the most commonly used framework for tracking applications. In Bayesian filtering, the tracking algorithm recursively estimates the target state [??](k | k) using the quantized measurement sequence z(1 : k), where z(1 : k) = [z(1),z(2),... ,z(k)}, namely, the posterior distribution p[x(k) | z(1 : k)}. It is well known that to recursively calculate the posterior distribution, one needs to have three distributions, namely, the initial state distribution p{x(0)}, the state transition model p{x(k) | x(k - 1)} that represents the state dynamics, and the likelihood function p{z(k) | x(k)} that depends on the measurement model [24].

Traditionally, Kalman filter provides an optimal solution to the Bayesian sequential problem for linear/Gaussian systems. In the case of nonlinear systems with Gaussian noises, extended Kalman filter can be used to provide a suboptimal solution by linearization of the nonlinear system equation and/or measurement model locally. Recently, however, it has been shown that even for linear/Gaussian systems, when the sensor measurements are quantized, extended Kalman filtering fails to provide an acceptable performance especially when the number of quantization levels is small [25]. For the problem of tracking with quantized measurements, in addition to the nonlinear mapping from the garget state to the sensor measurements (2), the final measurement fused model at the FC consists of quantized mapping resulting in a highly nonlinear and non-Gaussian system. Therefore, we adopt the sequential importance resampling (SIR) particle filter to solve the target state estimation problem (see, e.g., [24, 26,27] and the references therein). The advantage of the SIR particle filtering is that it is very easy to implement and computationally more efficient compared to other variants of particle filtering. As can be seen from the simulation in Section 7, the estimated performance by SIR particle filtering is quite satisfactory for tracking with quantized measurements, or in other words it performs very closely to its theoretical bound, namely, the posterior Cramer-Rao lower bound (CRLB) that will be derived in the second part of this section. Note that, for the three distributions mentioned above, the only remaining problem is the likelihood function p{z(k) | x(k)} which depends on both the sensor measurement model and the quantization strategy. With the help of the quantization model (4), the only thing we can infer from the quantized measurement [z.sub.i](k) = [[tau].sub.j] is that the original measurement [[tau].sub.j-1] < [y.sub.i](k) < [[tau].sub.j+1]. That is,

[[tau].sub.j-1] < (k) = [y.sub.i] (x (k) ,k) + [v.sub.i] (k) < [[tau].sub.j+1]. (27)

Therefore, for a given time step, we can use the sensor measurement model (9) and the Gaussian noise assumption with the inference in (27) to derive pizt(k) | x(k)}, the probability of a quantized sensor measurement taking a specific value [z.sub.i](k) = [[tau].sub.j] conditioned on the target state

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (28)

where [PHI] is a cumulative Gaussian distribution with mean zero and variance 1. Since sensor noises are assumed to be independent, the likelihood function at the FC can be formulated as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (29)

Using (29), the SIR particle filtering algorithm can now be applied to estimate the target state.

Now, we are in position to analyze the posterior CRLB for target tracking in WSNs with quantized measurements. Let [??](k I k) be an unbiased estimate of the state vector x(k) at time instant k, based on the quantized measurement sequence z(1 : k) and the prior knowledge of initial density p{x(0)}, where z(j) = [[[z.sub.1](j), [z.sub.2] (]), ... , [z.sub.N](j)].sup.T], j = 1,2, ...,k, are the augmented quantized measurement. The covariance of the estimation error is bounded by the posterior CRLB:

P(k | k) = E{[[??](k | k)-x (k)] [[[??](k | k)-x (k)].sup.T]} [greater than or equal to] [J.sup.-1.sub.k], (30)

where [J.sub.k] is the corresponding Fisher information matrix (FIM) and [J.sup.-1.sub.k] is the posterior CRLB. Tichavsky et al. provided an elegant method of computing the Fisher information matrix [J.sub.k] recursively [28]:

[J.sub.k+1] = [D.sup.22.sub.k] - [D.sup.21.sub.k] [([J.sub.k] + [D.sup.11.sub.k]).sup.-1][D.sup.12.sub.k], (31)

where

[D.sup.11.sub.k] = E {-[[DELTA].sup.x(k).sub.x(k)] log pix(k+1) | x (k)}} (32)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (33)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (34)

Here, [[DELTA].sup.x(k+1).sub.x(k)] = [[nabla].sub.x(k)][[nabla].sup.T.sub.t(k+1)] with V the first-order partial derivative operator. The recursion of (31) starts from an initial FIM [J.sub.0], which can be obtained from the a priori PDF p{x(0)}

[J.sub.0] = E{-[[DELTA].sup.x(0).sub.x(0)]) log p {x(0)}}. (35)

Note that the expectation E{x} in (32)-(33) is with respect to x(k) and x(k+ 1), whereas in (34) it is with respect to x(k), x(k +1), and z(k +1). From 31)-(34), it is not difficult to see that the quantized measurements only affect [D.sup.22,b.sub.k]. However, it is noticed that obtaining [D.sup.22,b.sub.k] analytically is not an easy task. A Monte Carlo simulation (see, e.g., [27]) is applied to approximate the theoretical posterior CRLB formulas, as can be found in Section 7.

7. A Simulation Example

Suppose S = 225 sensors are randomly deployed in the region of interest, which is 50 m x 50 m with the coordinate from (-25, -25) to (25, 25). Consider a target moving on noise circular trajectories as depicted in Figure 1. We use the CV to model the target dynamics [6, 7] with step-size T = 025. The measurement model of the ith sensor is [13]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (36)

where a = 40 is the assumed known amplitude of the sound source (the target), [parallel](x(k), y(k)) - (xs(i),ys(i))[parallel] denotes the distance between the target and the ith sensor, and [v.sub.i](k) is the measurement noise with distribution p([v.sub.i](k)) = [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The sensor positions are assumed to be known in the FC.

The layout of the network is illustrated in Figure 1, where each little square stands for the location of a sensor. The nearest-neighbors multisensor scheduling scheme is used to perform the tracking task. In each time step, the sensors lying within the sensing radius are selected to form a temporary tasking group. All the member sensors in the group are assumed to perform the sensing task simultaneously. The leader is selected as the one with the minimum average communication distance from all the activated neighbors (meanDj) and the maximum energy residue ([E.sub.ri]). The selection can be based on arg min; [[delta] x mean [D.sub.i] + (1-d)[E.sup.-1.sub.ri]], where [delta] is a weight on the mean distance. In the simulation, we set [delta] = 05. Then, each local processor quantizes the mea surement to a finite-bit message according to the bandwidth allocation through (18) or 25) and transmits this message to the leader. Finally, the leader, which functions as a local FC in addition to its basic sensing function, estimates the target state according to the received quantized measurements. We assume that the sensor messages are perfectly received by the FC and that all of the sensors deployed into an area have the same characteristics, that is, the same sensing radius [[gamma].sub.s] = 8 m, the measurement noise covariance [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and the measuring amplitude range [0, W] with W = 20.

In Figure 1, the red solid curve is the track of the target traveling start from (15, -10), while the black and blue curves represent the estimate of the target state through the proposed augmented and weighted approaches, respectively. The initial distribution of the sound source is assumed to be Gaussian with mean [[13, -8].sup.T] with covariance matrix 5*[I.sub.2]. In the FC, the augmented and weighted approaches are implemented with optimization problems (14) and (19), respectively. Each activated sensor quantizes the local measurement according to the allocated bandwidth and transmits it to the FC. Then, the SIR particle filtering is used to estimate the state of the target according to (1)and(9). As can be seen from Figure 1,both fusion approaches are able to keep the trackvery close to the true target trajectory, while the weighted approach performs a little worse than the augmented approach. This is because the latter uses the quantized messages directly without any information while for the former the quantized messages are first weighted, which will bring some information loss.

We perform the simulation over 100 Monte Carlo runs each with 200 time steps. For performance analysis and comparison purpose, we first derive the posterior CRLB in this scenario.

7.1. The Posterior CRLBs in the Linear Dynamic Case. For the CV model and fused measurement model (9), (32)-(34) become

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (37)

The substitution of (37) into the recursion (31) and the application of the matrix inverse lemma yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (38)

In this case, the recursion starts from [J.sub.0] = P-1, where P0 is the initial state variance.

Remark 5. When unquantized measurements are used, it is not difficult to derive that [29]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (39)

where [nabla]h is the partial derivation of the nonlinear measurement model.

While in our quantized scenario, it is not easy to obtain [D.sup.22,b.sub.k] = E{-[[DELTA].sup.x(k+1).sub.x(k+1)] log p{z(k + 1) | x(k)}}. The following will show this in detail. Combining the measurement model (36) and the likelihood function at the FC in (29) with the properties of first-order Markov system [30], we are ready to formulate [D.sup.22,b.sub.k] as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (40)

where, for m, n = {1,2}, [x.sup.(m).sub.(k+1)] stands for the mth entry of the state vector x(k + 1) and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (41)

Hence, with the CV system model and measurement model (36), we have, for m = {1,2},

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (42)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (43)

Till now, we have shown how to calculate the FIM for the state estimation with quantized measurements. However, as mentioned above, obtaining this quantity analytically is not an easy task. Hence, a Monte Carlo simulation (see, e.g., [27]) is applied to approximate the theoretical posterior CRLB formulas. With the help of particle filtering (PF) techniques, each PDF in the expectation E{-[[DELTA].sup.x(k+1).sub.x(k+1)] log p{z(k +1) | x(k)}} can be represented by a set of samples (particles) with associated weights after an M-run Monte Carlo simulation; that is,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (44)

where s is the number of particles and [x.sup.(l).sub.(k+1)] and [z.sup.(l).sub.i(k+1)] are the Ith realization of x(k + 1) and zt(k + 1) during the jth Monte Carlo run, respectively. The posterior CRLB can be obtained by simply taking the inverse of the FIM.

7.2. Results and Discussion. The main performance criterion is the root mean square error(RMSE) of the target position estimates defined by [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. In Figure 2, the RMSEs of position for both augmented and weighted approaches are compared with the posterior CRLB using quantized measurements. As less tightly bound, the posterior CRLB for state estimation using unquantized measurements is also provided in Figure 2. The performances are based on 100 Monte Carlo runs. From Figure 2, we can see that the augmented approach performs very closely to the posterior CRLB using quantized measurements. It is clear that the augmented approach significantly outperforms the weighted approach since the raw quantized messages are used for the former. Besides, the difference between the augmented and weighted approaches lies in that the latter merges the quantized measurements into a large vector, which will bring large computational burden in the FC. This will be demonstrated immediately (cf. Table 1). The RMSE in y-direction performs similarly and, therefore, is omitted here.

Tracking accuracy does not denote everything. It is especially true for target tracking through WSNs because the energy and bandwidth are both strictly limited. Therefore, another purpose of our simulations is to observe how the percentage of activated sensors affects the tracking accuracy and the amount of energy saving. We varied the percentage of activated sensors lying in the sensing radius from 10% to 100%. For each scenario, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] sensors are activated, where P is the percentage, [N.sub.s] is the number of sensors lying within the sensing radius, and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is the operator of upper integer rounding. Furthermore, we define the mean of RMSE over sampling time as mean RMSE = (1 /K) [[summation].sup.K.sub.i=1] [RMSE.sub.i]. The mean RMSEs for both augmented and weighted approaches are shown in Figure 3. For comparison purpose, we also plotted the posterior CRLB with quantized measurements. Commonly, the mean RMSEs and CRLB decrease as the percentage of activated sensors increases. Furthermore, for all scenarios, the augmented approach also performs very closely to the CRLB and significantly outperforms the weighted approach.

Moreover, we compare the proposed power scheduling schemes (18) and 25) to the uniform quantization allocation scheme, in which each sensor quantizes the observation into the minimal same number of bits to achieve the target RMSE distortion [16]. The comparison of communication energy between the proposed scheduling strategies in (18) and (25) and the uniform quantization allocation scheme is shown in Figure 4, in which we also vary the percentage of activated sensors that lie in the sensing radius from 10% to 100%. The mean of communication energy reduction for (18) and (25) is 45.1% and 72.8% when the percentage of activated sensors is 100%, respectively. This will prolong the lifetime of the whole network abundantly for both fusion approaches. Besides, from Figure 4, it is obvious that the energy saving of (18) and (25) over uniform quantization allocation scheme increases with the percentage of sensors increase. In other words, the proposed bandwidth scheduling schemes save more energy in densely deployed sensor networks. Finally, computational complexity of the two approaches in the FC is compared in Table 1, which shows that the weighted approach has less computational burden than the augmented approach. Compared with the latter, the former saves more than 31% computational energy in the FC.

8. Conclusion and Future Work

The energy-constrained bandwidth scheduling problem for collaborative target tracking in WSNs has been investigated from the quantized measurement fusion perspective. Quantized measurement fusion via both augmented approach and weighted approach has been investigated with the focus on tradeoff between total energy consumption and fusion performance. The closed-form solution to the quantized-measurement-fusion-based bandwidth scheduling has been given for both approaches. Posterior CRLB on the tracking accuracy using quantized measurements is also derived. Future work will be concentrated on the tracking algorithm and the performance analysis under transmitting uncertainties.

http://dx.doi.org/10.1155/2014/682032

Conflict of Interests

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

Acknowledgments

The work was supported by the National Natural Science Foundation of China (61104210, 61100140, 61372049, and 61104072), China Scholarship Council (CSC), and the construct program of the key discipline in Hunan Province.

References

[1] I. F. Akyildiz, W. Su, Y. Sankarasubramaniam, and E. Cayirci, "Wireless sensor networks: a survey," Computer Networks, vol. 38, no. 4, pp. 393-422, 2002.

[2] J. Yick, B. Mukherjee, and D. Ghosal, "Wireless sensor network survey," Computer Networks, vol. 52, no. 12, pp. 2292-2330, 2008.

[3] M. Valera and S. A. Velastin, "Intelligent distributed surveillance systems: a review," IEE Proceedings--Vision, Image and Signal Processing, vol. 152, no. 2, pp. 192-204, 2005.

[4] Q. Wang, W.-P. Chen, R. Zheng, K. Lee, and L. Sha, "Acoustic target tracking using tiny wireless sensor devices," Proceedings of the International Workshop on Information Processing in Sensor Networks, vol. 2634, pp. 642-657, 2003.

[5] X. Li, W. Shu, M. Li, H.-Y. Huang, P.-E. Luo, and M.-Y. Wu, "Performance evaluation of vehicle-based mobile sensor networks for traffic monitoring," IEEE Transactions on Vehicular Technology, vol. 58, no. 4, pp. 1647-1653, 2009.

[6] J.-J. Xiao, S. Cui, Z.-Q. Luo, and A. J. Goldsmith, "Power scheduling of universal decentralized estimation in sensor networks," IEEE Transactions on Signal Processing, vol. 54, no. 2, pp. 413-422, 2006.

[7] Y. Zhou, J. Li, and D. Wang, "Collaborative target tracking in wireless sensor networks using quantized innovations and sigma-point Kalman filtering," in Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE '09), pp. 942-947, Seoul, Korea, July 2009.

[8] J. Aslam, Z. Butler, F. Constantin, V. Crespi, G. Cybenko, and D. Rus, "Tracking a moving object with a binary sensor network," in Proceedings of the 1st International Conference on Embedded Networked Sensor Systems (SenSys '03), pp. 150-161, November 2003.

[9] F. Zhao, J. Shin, and J. Reich, "Information-driven dynamic sensor collaboration," IEEE Signal Processing Magazine, vol. 19, no. 2, pp. 61-72, 2002.

[10] W. Zhang and G. Cao, "DCTC: dynamic convoy tree-based collaboration for target tracking in sensor networks," IEEE Transactions on Wireless Communications, vol. 3, no. 5, pp. 1689-1701, 2004.

[11] J. A. Gubner, "Distributed estimation and quantization," IEEE Transactions on Information Theory, vol. 39, no. 4, pp. 1465-1467, 1993.

[12] W.-M. Lam and A. R. Reibman, "Design of quantizers for decentralized estimation systems," IEEE Transactions on Communications, vol. 41, no. 11, pp. 1602-1605, 1993.

[13] A. Ribeiro, G. B. Giannakis, and S. I. Roumeliotis, "SOI-KF: distributed Kalman filtering with low-cost communications using the sign of innovations," IEEE Transactions on Signal Processing, vol. 54, no. 12, pp. 4782-4795, 2006.

[14] P. M. Djuric, M. Vemula, and M. F. Bugallo, "Target tracking by particle filtering in binary sensor networks," IEEE Transactions on Signal Processing, vol. 56, no. 6, pp. 2229-2238, 2008.

[15] Y. Zhou, J. Li, and D. Wang, "Posterior Cramer-Rao lower bounds for target tracking in sensor networks with quantized range-only measurements," IEEE Signal Processing Letters, vol. 17, no. 2, pp. 157-160, 2010.

[16] Y. Zhou, J. Li, and D. Wang, "Target tracking in wireless sensor networks using adaptive measurement quantization," Science China Information Sciences, vol. 55, no. 4, pp. 827-838, 2012.

[17] Q. Gan and C. J. Harris, "Comparison of two measurement fusion methods for Kalman-filter-based multisensor data fusion," IEEE Transactions on Aerospace and Electronic Systems, vol. 37, no. 1, pp. 273-280, 2001.

[18] H. L. Van Trees and K. L. Bell, Bayesian Bounds for Parameter Estimation and Nonlinear Filtering & Tracking, Wiley-IEEE Press, 2007.

[19] H. Meng, M. L. Hernandez, Y. Liu, and X. Wang, "Computationally efficient PCRLB for tracking in cluttered environments: measurement existence conditioning approach," IET Signal Processing, vol. 3, no. 2, pp. 133-149, 2009.

[20] X. R. Li, Y. Zhu, and C. Han, "Unified optimal linear estimation fusion," in Proceedings of the 39th IEEE Conference on Decision and Control, pp. 10-25, Sydney, Australia, 2000.

[21] S. Sun, J. Lin, L. Xie, and W. Xiao, "Quantized Kalman filtering," in Proceedings of the 22nd IEEE International Symposium on Intelligent Control (ISIC '07), pp. 7-12, Singapore, October 2007

[22] S. Cui, A. J. Goldsmith, and A. Bahai, "Energy-constrained modulation optimization," IEEE Transactions on Wireless Communications, vol. 4, no. 5, pp. 2349-2360, 2005.

[23] S. Boyd and L. Vandenberghe, Convex Optimization, Cambridge University Press, Cambridge, UK, 2003.

[24] M. S. Arulampalam, S. Maskell, N. Gordon, and T. Clapp, "A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking," IEEE Transactions on Signal Processing, vol. 50, no. 2, pp. 174-188, 2002.

[25] Y. Ruan, P. Willett, A. Marrs, F. Palmier, and S. Marano, "Practical fusion of quantized measurements via particle filtering," IEEE Transactions on Aerospace and Electronic Systems, vol. 44, no. 1, pp. 15-29, 2008.

[26] N. J. Gordon, D. J. Salmond, and A. F. M. Smith, "Novel approach to nonlinear/non-Gaussian Bayesian state estimation," IEE Proceedings, Part F, vol. 140, no. 2, pp. 107-113, 1993.

[27] B. Ristic, S. Arulampalam, and N. Gordon, Beyond the Kalman Filtering: Particle Filtering for Tracking Application, Artech House, 2004.

[28] P. Tichavsky, C. H. Muravchik, and A. Nehorai, "Posterior cramer-rao bounds for discrete-time nonlinear filtering," IEEE Transactions on Signal Processing, vol. 46, no. 5, pp. 1386-1396, 1998.

[29] H. L. Van Trees, Detection, Estimation and Modulation Theory, John Wiley & Sons, New York, NY, USA, 1968.

[30] Y. Bar-Shalom, X. R. Li, and T. Kirubarajan, Estimation with Applications to Tracking and Navigation, John Wiley & Sons, New York, NY, USA, 2001.

Yan Zhou, (1,2) Dongli Wang, (1) Tingrui Pei, (1) and Yonghong Lan (1)

(1) College of Information Engineering, Xiangtan University, Xiangtan 411105, China

(2) Key Laboratory of Intelligent Computing & Information Processing of MOE, Xiangtan University, Xiangtan 411105, China

Correspondence should be addressed to Yan Zhou; yanzhou@xtu.edu.cn

Received 8 October 2013; Revised 27 December 2013; Accepted 6 January 2014; Published 27 February 2014

Academic Editor: Long Cheng

TABLE 1: The computational complexity. Algorithms Augmented Weighted Time (ms) 501.6 344.8 Ratio 1.0 0.6874

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Zhou, Yan; Wang, Dongli; Pei, Tingrui; Lan, Yonghong |

Publication: | International Journal of Distributed Sensor Networks |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 6775 |

Previous Article: | Histogram estimation for optimal filter skyline query processing in wireless sensor networks. |

Next Article: | Distributed relay-assisted retransmission scheme for wireless home networks. |

Topics: |