An Energy-Efficient Collaborative Target Tracking Framework in Distributed Wireless Sensor Networks

Energy consumption and tracking accuracy are two significant issues for collaborative tracking in distributed wireless sensor networks (DWSNs). To obtain a benefit from those issues, most of the recent work tends to reduce the spatial redundancy, while ignoring utilizing the attribute of time redundancy. In this paper, a novel energy-efficient framework of collaborative signal and information fusion is proposed for acoustic target tracking. The proposed fusion algorithm is based on neural network aggregation model and Gaussian particle filtering (GPF) estimation. And the neural network based aggregation (NNBA) can reduce spatial and time redundancy. Furthermore, a fresh cluster head (CH) selection method demanding less task handover is also presented to decrease energy consumption. The analyzed framework coupled with simulations demonstrates its excellent performance in tracking accuracy and energy consumption.


Introduction
With the fast development of microelectromechanical system (MEMS), digital signal processing, and embedded computing, sensor nodes are becoming smaller and cheaper.Therefore, it is possible to deploy large-scale distributed wireless sensor networks (DWSNs) to "achieve quality through quantity" in many complex applications [1,2].
However, challenges and difficulties also exist in target tracking DWSNs [3][4][5]: (1) sensor resources limitation [6] (especially communication bandwidth, energy, sensing range, and processing capabilities); (2) deployment and coordination of a mountain of sensors; (3) redundant data caused by similar measurements of adjacent sensors.Accomplishing the tracking task requires not only the methods for data fusion and collaborative signal and information processing (CSIP), but also an efficient scheme of energy consumption.In the literature, various researches have been proposed for target tracking in DWSNs [7].
In [8,9], the nodes detect the mobile target and send measurements directly to the sink node so that it can estimate the target location using all data.But the centralized schemes, regarded as the first CSIP (CSIP1), depend on the sink node to collect and process measurements, as the sink node will consume a lot of energy and be apt to die quickly with constraint battery; the network lifetime is consequently shortened.
Energy aware probabilistic target localization algorithm for cluster-based WSN is proposed in [10].Sensors detecting a target report a short message to the cluster head (CH).Then the CH localizes and determines the subset of sensors in the vicinity of target.This scheme can reduce communication bandwidth and save lots of energy.Nevertheless, the CH undertakes a heavy tracking task, which shortens its lifetime.The work in [11] proposed an adaptive dynamic clusterbased tracking (ADCT) protocol in which CHs are selected dynamically; meanwhile sensor nodes surrounding the target are waken up to construct a cluster with the target's moving in the network.Such scheme can be called the second CSIP (CSIP2).However, the election of the first CH is achieved by competition, which causes an additional burden in communication bandwidth.
A distributed predictive tracking (DPT) protocol for tracking moving objects with separate algorithms for nodes and CHs is presented in [12].The CH adopts the target descriptor to identify target and predicts its next location.Then it informs predicted area CH who activates  sensors 2 International Journal of Distributed Sensor Networks to sense the target.But most of the efforts have been spent on selecting the node nearest to the predicted target position for the CH.This scheme can be considered as the third CSIP (CSIP3).
Moreover, the work in [13] showed three algorithms for tracking a mobile target in WSNs utilizing cluster-based architecture, namely, adaptive head, static head, and selective static head.In [14], an intensive review was conducted to improve energy efficiency in collaborative target tracking in DWSNs.
Previous work was more focused on the spatial redundancy of the distributed sensor network for target tracking.However, the study on the time redundancy is still important.The purpose of this paper is to design a novel energy-efficient collaborative signal and information fusion processing (CSIFP) framework for acoustic target tracking in DWSNs, to collaboratively estimate the position and velocity of a moving target.In our framework, sensor nodes detecting a target are activated and dynamically form a cluster around the target.Each member node transfers the deredundant data to the CH.CH aggregates and processes these messages to estimate the current state of the target and reports the result to the sink.The CH is reselected relying on certain rule in Sections 4.1 and 4.2.To estimate the target position, the Gaussian particle filter (GPF) algorithm [15] is adopted.Compared with the particle filter (PF) [16][17][18], it omits resampling to reduce the computational complexity and avoid transmitting all particles and weights when task handover happens.Furthermore, the neural network based aggregation (NNBA) model relies on the cooperation of the nodes, which can reduce the spatial and time redundancy.
The paper is organized as follows.Section 2 presents the system assumptions and the system models.Section 3 proposes the framework of the distributed collaborative target tracking framework is illustrated.Section 4 discusses a neuron network based aggregation model for target tracking.Section 5 describes the simulation environment and Section 6 shows and analyzes the results.Finally, Section 7 summarizes the work and gives a conclusion.(v) As the paper focuses on the energy consumption of target tracking in DWSNs, the energy consumed by sensor nodes in the sleeping state is not considered.However, it does not mean that the energy consumption of idle sensor nodes can always be ignored.

Target Dynamic Model and Measurement
Model.Target tracking aims to estimate the real-time state of a moving object, which is mainly based on the target dynamic models.The appropriate models can not only describe the evolution of the target state with respect to time, but also represent a large number of observations, reducing the data transmission in the DWSNs [19].
The target state is denoted by a four-dimensional vector in 2D Cartesian coordinate plane: where   and   denote the target position in -axis and -axis at time   and ξ and η denote the corresponding velocities.
We describe the target dynamic model with addictive noise as the state-space model in the following form: where x  , z  , and u  are the target state, observation, and control input vectors, respectively, at the discrete time   , w  is the system noise at   , and k  is the observation noise, assumed to be statistically independent from w  .  and ℎ  are state transition and observation functions, which are vector valued (possibly time-varying).
For nearly constant velocity (CV) or white noise acceleration model [19], u  is approximated to zero, and (2) can be rewritten by where Note that  is the sampling period.[20][21][22].Energy based acoustic sensors are widely used in DWSNs due to low cost and universality.The intensity of the source attenuates at a rate that is inversely proportional to the distance from the source.Given   sensors, the signal intensity measured on the th sensor is given by

Measurement Model
where  is the strength of the acoustic emitted by the target during .The vectors r  and r  denote the coordinates of the target and the th sensor at time step   , respectively.‖r  − r  ‖ is the Euclidean distance between the th sensor node and the target.  is sensor gain factor of the th sensor, which can be measured during sensor calibration. is the intensity decay factor.In [21], it is shown that the effective decay, , is approximately 2 and this value is assumed for this paper.
V , is white Gaussian measurement noise, V , ∼ (0,  2 V ), independent of   , and independent of other sensors' noise.In practice,  2  V can be estimated based on the measurements when the acoustic source is absent.For concise notation,  2 V is assumed to be identical for all sensor nodes,  2 V =  2 V for 1 ≤  ≤   , and it is assumed to be known.

Distributed Collaborative Target Tracking Framework
In this paper, we develop an energy-efficient distributed collaborative target tracking framework.This effective framework involves four modules: the tracking task initialization, CH election, node selection, and tracking task handover.Before specifying the aforementioned contents, several messages are listed in Table 1.

Tracking Task Initialization.
Assume that a network with   sensors is deployed on a 2D area.For a given acoustic energy sensor node, it can detect a target once the acoustic intensity measurement exceeds a threshold, th, corresponding to the calculated strength at the maximum sensing distance   .Sensor nodes detecting the target are activated, whereas the other nodes remain in the sleeping state for saving energy.Considering the influence of the measurement noise, we define a set of candidate nodes whose measurements of acoustic intensity are not lower than  th / and remaining energy is greater than zero.Only the candidate nodes could participate in the target tracking process and transmit the measurements to the CH.The candidate nodes set  is denoted as follows: where  , is the measurement of sensor node  at time   and the energy threshold  th is  ⋅ (/ 2  ).Parameter  ∈ (0, 1] is customized, and a bigger  means more candidate nodes will be chosen.  is the remaining energy of the node , and   is the communication radius.
At initial time step as well as the target detection restart phase, each candidate sensor node starts a timer   and will broadcast CHInfoMesg after the timer expires.The length of   is According to (6), the sensor with the largest measurement is nearest to the target; thus,   is the smallest.So the sensor node nearest to the target will first send its CHInfoMesg.The other nodes which have received this message register the broadcaster as their CH, terminate the timer, and give up TargetDegradeMesg CH transmits the result to the sink broadcasting message.Through the above process, the candidate closest to the target is elected to be the CH.The cluster member nodes transmit their observation to CH, and the CH aggregates the collected information and performs the tracking task using GPF.

CH Selection.
It is expectable that CHs will consume much more energy than others because they undertake more burdensome tasks.For balancing the energy consumption, each node needs to take turn as CH by certain scheduling algorithm.However, CH rotation process consumes a lot of energy for transmitting message between CHs.Therefore, we propose a novel CH selection scheme considering these two factors.
At time  −1 , it is assumed that the CH is CH −1 .According to the predicted target position at time   and cluster members information, CH −1 builds a candidate CH set   : where  ∈ [0, 1] is customized and a smaller  means CH changes more frequently.  is the remaining energy of node , and  th is the threshold energy.If   is lower than  th , a node is not suitable for being a cluster head.The current CH checks itself whether it belongs to the candidate set   or not.If so, it will continue taking charge of the tracking task; otherwise, it will hand over the CHHan-doverMesg including tracking information and necessary parameters to CH  .In order to select an appropriate new CH, we define the attraction force   of the node , exerting on the current CH.The larger value means the node  has greater probability to be chosen as next CH: Note that  dir, ,  dis, , and  power, are the direction attraction component force, distance attraction component force, and power attraction component force exerting on CH −1 by node , respectively., , and  are three nonnegative constants which adjust the ratios of the above three component forces, and  +  +  = 1.Consider the following.
(i) Direction attraction component force  dirt, : the direction from the current cluster head CH −1 to node International Journal of Distributed Sensor Networks  and the direction from CH −1 to the predicted target position x are more closer, which means that the node  is nearly on the direction of the target motion and has a high probability to be the cluster head at time step   .We define the cosine value of this angle to indicate the force  dirt, : where r CH −1 is a vector from CH −1 to the predicted target position x and r  is vector from the current cluster head CH −1 to node .⟨r CH −1 , r  ⟩ is the inner product.‖r CH −1 ‖ and ‖r  ‖ represent the length of r CH −1 and r  , respectively.
(ii) Distance attraction component force  dis, : for the sake of reducing the frequency of CH reselection, we choose the candidate node as far as the current CH on the premise of an acceptable tracking accuracy. dis, is defined as follows: where dis (, CH −1 ) is the Euclidian distance between the node  and the current cluster head CH −1 and   is the communication radius.
(iii) Power attraction component force  power, :  power, is defined to indicate the node battery level, which is given by where  0 is the initial energy of the node .
Finally, CH  is selected according to the following rule: Noting that it is possible for multiple candidate nodes to have the same maximum attraction force, in this case, we will select one node randomly among these nodes as the new CH.

Member Node Selection.
The node selection determines the total energy consumption and data fusion accuracy, which has a significant impact on the performance of the whole sensor network.
To implement an optimal cluster with member nodes, more nodes should be selected to improve the tracking accuracy.On the other hand, fewer nodes should be selected to reduce the energy consumption for transmission.So we need to make a compromise.Considering the influence of measurement noise, we choose the nodes located in the communication range of the current CH from the set of candidate nodes   .The set of member nodes are denoted by where dis (, CH) is the distance between the node  and the position of the current CH at   .The constraint dis (, CH) ≤   is resistant to the measurement noise.The current CH selects a member node from   as the new CH, CH  , which is along the direction of estimated target motion as well as far away from the current CH.Then CH −1 sends a CHHandoverMesg message to CH  .This message contains the current target state estimation as well as some necessary parameters.After receiving the CHHandoverMesg from CH −1 , CH  will continue to take charge of the target tracking task and broadcast CHInfoMesg to inform its cluster members.

Working Scheduling.
The distributed collaborative target tracking framework of the CH and member nodes is shown in Figure 1.Consider the following.
(i) The cluster head: CH is generated or reselected according to (14), and it broadcasts CHInfoMesg message to inform its cluster members.CH receives the DataMesg packet or DataDegradeMesg packet from its members, then chooses three messages with the biggest acoustic intensity measurement, and performs the GPF algorithm described in Section 4.  Note that, at any time, only several sensor nodes in the vicinity of the target need to collect the acoustic signal and the others keep sleeping.
The CH and member nodes run in the neural network based aggregation tracking model, which will be expatiated in Section 4.

Neural Network Based Aggregation Tracking Model
This section details the fusion processing algorithm based on neural network based aggregation tracking model.It combines the neural network of multilayer perceptron with clustering routing structure.This algorithm reduces time redundancy and spatial redundancy to reduce the traffic, save energy, and prolong the network lifetime.

4.1.
Collaborative Target Tracking Scheme.Figure 2 gives the diagram of NNBA model structure [23], which uses fourperceptron neural network.For DWSNs, each sensor node works as an underlying neuron; CH is like a central neuron, and the whole network is a complex neuron network system.Cluster member nodes are in the input layer and the first hidden layer.CH nodes are in the second and the third  hidden layer and output layer.For instance, each cluster with  members collects data of  species by different types.Consequently, this NNBA model has  ×  input layer nodes and the same number of first hidden layer neurons.Here, we only use the acoustic sensor, so  = 1.The quantity of neurons in the second hidden layer, the third hidden layer, and the output layer can be adjusted by the users according to the practical application.The hidden layers have no inevitable relationship with .It only deals with separate data of different types and need not be fully connected between the first and the second hidden layer.So do the second and the third hidden layer.But the treated data between the third hidden layer and output layer can be integrated, for they are full-connected.
According to this NNBA model, firstly sensor nodes collect the acoustic intensity measurements from the target and preprocess them according to the neuron function of the first hidden layer.Then the nodes send the results to the CH.The CH further processes to track the target using GPF algorithm detailed in Section 4.2.GPF is embedded in the second hidden layer, the third hidden layer, and output layer neurons.Finally, the CH transfers the estimated target state to the sink.Note that neural network model can do some modification according to the specific requirements in various applications.

Gaussian Particle
Filter.Recalling the theory of particle filtering [16][17][18], we estimate the target a posteriori distribution (x 0: | z 1: ) of x 0: , by approximating it with a random measure composed of particles x   and weights    , where  is a particle index.We denote it by Ξ = {x  0: ,    }  =1 with  being the number of particles.At each time step   , the particle filter carries out the following operations: (1) selection of most promising particle set, (2) particle propagation, (3) computation of particle weights, and (4) state estimation.
The Gaussian particle filter (GPF) [15] approximates the posterior (x  | z 0: ) by a Gaussian pdf, whose mean and covariance are calculated from the weighted particle set.In contrast to the standard PF [16][17][18], GPF algorithm does not require the process of resampling, which reduces the computational complexity.
Assume that the initial target state is known and its distribution accords with Gaussian distribution (x 0 ; x 0 ,  0 ).The steps of GPF algorithm are described as follows.

Process function
Output function (c) Normalize the weights as (d) Calculate the mean   and the covariance Σ  of approximate Gaussian distribution by  6), the sensor data is valuable when the distance ‖r  − r  ‖ between sensor node and the target is within certain limits, so the threshold  th / guarantees the accuracy of measurement.And the second function checks and decides whether to record this current input data.The output function is a specific threshold function, which is different from the standard threshold function with output 0 or 1.This function outputs 2() if 1() ≥ 1, else output none, to reduce the time redundancy.

The Second Hidden Layer.
The second hidden layer on the CH merges with the GPF algorithm (see Figure 4).CH chooses three maximum measurements and puts them in neurons, respectively, in the second hidden layer.According to GPF algorithm, CH draws samples {x  −1 }  =1 from the distribution (x −1 ) ∼ (x −1 ;  −1 , Σ −1 ), in which the mean  −1 and covariance Σ −1 can be computed by (16).In the process function, each neuron continues to complete the time update process, draws samples from the importance function (x  | z 0: ), and obtains the respective particle weights by (17) based on input parameter ,  −1 , Σ −1 and the same particle set {x  −1 }  =1 .At last, each neuron outputs its own normalization weights.

The Third Hidden Layer.
One neuron is used to complete the subsequent GPF algorithm in the third hidden layer.As shown in Figure 5, CH averages value of the input 2, the weight of the three particle sets.The process function is used ) ) Input function

Process function
Output function Input function

Process function
Output function Figure 5: NNBA model third hidden layer neuron function.
to estimate the target position x and covariance Σ  according to (19).The output of the function equals the estimated target position x represented by 3, which will be transferred to the output layer.
The GPF tracking algorithm is completely performed by the second and third hidden layer.Note that these two can be combined into a union, in which case the input data should be aggregated into a vector.

The Output Layer.
The output neuron model resembles the input neuron.Different from the input neuron, the output neuron model can be multi-input neuron, which is determined by the sensor type.In this study, only acoustic energy sensor is used, so there is one neuron in output layer.
As Figure 6 shows, CH compares the input   and (−1).If the difference exceeds the threshold 2, that is, Δ() ≥ 2, CH outputs   and sends the TargetToSinkMesg containing target state estimation to the sink; otherwise, the neuron outputs none, and CH sends TargetDegradeMesg message to the sink.The shortened messages reduce traffic and lower the energy consumption especially when the target stops or moves slowly.
Neurons in NNBA algorithm need some operation parameters (such as weight, thresholds, etc.) at runtime; as

Simulation Environment and Setup
In this section, we build a simulation platform to evaluate the performance of the proposed distributed collaborative target tracking framework.Firstly, we describe our simulation environment and parameters.Then we define the performance metrics to compare the proposed algorithm and others' work.We should also determine the thresholds  and , for they are playing an important role in selecting the candidate nodes and CH.The bigger the  is, the higher the tracking accuracy is, but the more energy the whole network consumes.The bigger the  is, the less frequently the CH hands over, and the more energy consumption the whole network has.On the other hand, the tracking accuracy is not closely tied with .According to large number of experiment results, a best tradeoff between good tracking accuracy and low energy consumption is achieved in the case of  = 0.64,  = 0.7.Simulation parameters are summarized in Table 2.

Tracking Performance.
Assume a target moves into the monitoring field at time  = 0 with initial state [50 6 0 6]  .The simulation tracking time is  sim = 100.The system process noise   is assumed as zero-mean Gaussian distribution with variance  2  = diag([0.1250.5 0.125 0.5]).The number of particles used in each algorithm is  = 600.
To evaluate the tracking performance for different algorithms, we calculate the root of mean square errors (RMSEs) International Journal of Distributed Sensor Networks of target position and velocity estimations at each time step under 50 Monte Carlo runs.The RMSEs metric [24] is defined as for position,

Energy Consumption.
The node energy depletion results from sensing module, communication module, and the processing module.The total energy consumption at time   can be summarized as follows: Among the three parts, wireless communication module plays the dominant role in the energy consumption.Many researches have been carried out on the low-energy radios.The first-order radio energy consumption model [25] is adopted for simulation.To transmit -bit message over a distance , the radio expends To receive -bit message, the radio expends where  -elec and  -elec , separately denote the energy to activate the baseband circuits when transmitting and receiving. -amp is the energy to run the RF module.The electronic energy,  elec , denotes the energy spent to activate the baseband circuits to transmit or receive 1 bit.And  fs and  mp denote the energy spent to run the RF module to transmit 1 bit in free space and multipath space separately. 0 is a threshold indicating that if the distance is less than  0 , the free-space model is used; otherwise, the multipath model is used.For the experiments, the communication energy parameters are set:  elec = 50 nJ/bit,  fs = 10 pJ/bit/m 2 , and  mp = 0.0013 pJ/bit/m 4 .From (23), we can derive  0 as The packet is defined in Table 3.In all algorithms, the HEADER of the listed messages in Table 3 is 8 bits.

Simulation and Analysis
This section shows experimental studies for verifying the proposed framework.
6.1.Performance of Target Tracking.In order to demonstrate the GPF's satisfied performance of tracking accuracy, a comparison with PF [17] is carried out below.Here, we consider two cases: in PF-1, the mean and covariance for Gaussian noise are modeled precisely; the noise statistical parameters in PF-2 have large model error.Figure 7 displays a typical trajectory and the tracking obtained by three different filtering algorithms.Figure 8 demonstrates the comparison of RMSEs of the estimated position and velocity obtained by three different filtering algorithms at each time step under 50 Monte Carlo runs ( mc = 50).
As shown in Figures 7 and 8, we can see that the PF-1 algorithm has the best tracking accuracy, while PF-2 works badly.It is because that PF is greatly influenced by model parameters.
It is noted that the GPF algorithm outperforms the PF-2 with a close performance to PF-1 algorithm.Though GPF omits the resampling process, it preserves the desired accuracy of target estimation.Instead of handing over all particles to the new CH in a general PF algorithm, GPF only needs to transmit parameters (the mean and the covariance of the posterior distribution), which saves amount of energy.Further discussion will be detailed in Section 6.2.
To further validate the effectiveness of our proposed algorithm CSIFP, we compare it with those obtained from previous work [8,11,12], including CSIP1, CSIP2, and CSIP3 algorithms.All of them use GPF estimation method.
The tracking trajectories obtained by these algorithms are shown in Figure 9.It is clear that the tracking performance of CSIFP is comparable to the referenced methods (CSIP1, CSIP2, and CSIP3).CSIP1 performs best because of relying on the largest set of sensors' measurement.However, a rise in energy consumption appears.CSIP2 and CSIP3 approach a moderate precision due to their own strategies of CH selection and spatial fusion.It makes sense that the proposed CSIFP introduces not only spatial fusion but also time fusion, sacrificing the tracking accuracy in some degree.It is noteworthy that the slight performance degradation is worthy in contrast to the significant energy saving benefit in Section 6.2.

Performance of Network Energy Efficiency.
In this subsection, we will firstly explain why the GPF algorithm is adopted by comparing generic PF's and GPF's energy consumption.Then the performance of energy efficiency in our proposed CSIFP is validated as well as the CSIP1, CSIP2, and CSIP3.
From Figure 11, the cumulative energy consumptions in PF-1 algorithm and PF-2 algorithm are nearly the same.It is obvious that energy consumption is uncorrelated with prior model of system noise.The GPF, by contrast, consumes least energy of the three.The theoretical proof is derived.
In the PF algorithm,  particles and their weights are transferred when tracking task handover happens, and total quantity of data transmission is where  represents the dimension of the target state vector and 1 ≤  ≤ 6 for the 2-D target tracking. is the number of the particles, usually in the 10 2 order of magnitude or more.
In the GPF algorithm, the current CH transmits nothing but the Gaussian mean and covariance of the target estimation.
Considering the symmetry of the covariance, the quantity of the communication data is It can be seen from Figure 12, the traffic of general PF algorithm is much bigger than GPF, and the latter's traffic has nothing to do with the particle number.The accuracy of the algorithm can be improved without a heavier traffic load by increasing the number of particles, if the computation complexity is accepted.
Figure 13 presents the total communication energy consumption of four different collaborative processing schemes.We can see that the CSIP1 consumes the most energy because all the active nodes with the   range send their measurements to CH.The energy consumption of CSIP2 and CSIP3 is less than CSIP1 by choosing the candidate  To further verify the tracking performance of the proposed algorithm, the position RMSE and velocity RMSE are, respectively, illustrated in Figure 15   Moreover, we carry out statistics by average RMSE and total energy consumption, as illustrated in Table 4. CSIFP achieves nearly the same good tracking accuracy.It obtains a gain of 56.83% in the energy cost in contrast with CSIP1, with a sacrifice of 12.88% in precision.Compared to CSIP2 and CSIP3, CSIFP consumes the least energy.This gain benefits from not only the CH selection method but also the fusion algorithm based on NNBA.From the above results, a good tradeoff between tracking accuracy and energy cost is achieved.

Conclusions
In this paper, we have presented a high energy-effective strategy of collaborative signal and information fusion processing for acoustic target tracking applications in DWSNs.Within this framework, the whole network is assumed to be a neural network, in which the member nodes correspond to the neurons and the CH corresponds to the neural center.The nodes detect the target and transmit measurements to the CH.CH receives the information, responds to different neurons' stimulation, estimates the target state based on the GPF, and sends the estimation to the sink finally.

3. 4 .
Tracking Task Handover.Once the target is far away from the current CH and CH −1 ∉   , CH −1 cannot efficiently manage.CH −1 should hand over the tracking task to a new CH.This operation can not only shut down the old nonsignificant cluster nodes to save energy in time but also improve the tracking accuracy.

Figure 1 :
Figure 1: Flowchart of the distributed collaborative tracking scheme.

5. 1 .
Simulation Environment.We simulate a wireless sensor network with   = 200 randomly deployed acoustic sensors to monitor a 800 m × 800 m square field.All the sensor nodes within the monitoring region (excluding sink) have same structure and can detect the target within the sensing range.The measurements collected by sensor nodes are corrupted by zero-mean i.i.d.Gaussian noise with variance  2 V = 0.00005.The parameters of the measurement model in (6) are  = 1000,   = 1.Besides that, for the CH selection phase in the proposed framework, weight constants in(10) are  = 0.4,  = 0.3, and  = 0.3; energy threshold is set to 10% of the initial battery level; namely,  th = 0.1 ×  0 .
are the target position estimation at time step   in th Monte Carlo run, and  true  ,  true  are the true position at   .Similarly, ̂ξ   , ̂η   are the target velocity estimation at time step   in th Monte Carlo run, and ξ true  , η true  are the true velocity at   .
Figure 10 demonstrates the average result of 50 Monte Carlo runs from up to down: (a) the RMSE

Figure 11 :Figure 12 :
Figure 11: The cumulative communication energy consumption of PF and GPF after the entire tracking process.

Figure 13 :
Figure 13: The cumulative communication energy consumption of CSIFP, CSIP1, CSIP2, and CSIP3 after the entire tracking process.
and Figure 16.The position RMSE curve of CSIFP has obvious fluctuation.It is because that CSIFP reduces time redundancy by NNBA model.NNBA model will be inertial if a slight shift happens in contrast to the previous measurement.

Table 1 :
Messages used in the target tracking framework.
2 to estimate the current target state and predict the possible state at next time step.The CH also takes charge of sending the tracking results to the sink node.If the difference between estimated target position at   and the one at  −1 is greater than the threshold 2, CH .If so, the node transfers DataMesg including packet type, ID, and measurement to the CH.Otherwise it goes to sleep state.When the difference of the calculated distance between the node  and the target at   and  −1 is smaller than difference threshold 1, it sends a DataDegradeMesg, omitting measurement.The CH holds its previous measurement once receiving DataDegradeMesg.

Table 2 :
The set of the simulation parameters.

Table 4 :
Average result of performance.The former saves more energy due to selecting the farther member node as the new CH when CH changes, which reduces the handovers.The proposed CSIFP scheme consumes the least energy, because the CH selection method makes advantage of CSIP2, and the NNBA model is utilized to reduce the time redundancy as well.Consequently, the network payload approaches an expectable level.6.3.Performance of Robustness.Furthermore, for investigating the robustness of our proposed algorithm CSIFP, 6 independent Monte Carlo simulations with different trajectories are carried out.Figure14presents these random trajectories and tracking results obtained by these four algorithms.It is clear to us that CSIFP can keep good tracking accuracy when the trajectory changes.