International Journal of Distributed Sensor NetworksVolume 2012 (2012), Article ID 460924, 14 pagesdoi:10.1155/2012/460924

Research Article

The Effects of Motion on Distributed Detection in Mobile Ad Hoc Sensor Networks

Xusheng Sun and Edward J. Coyle

School of Electrical and Computer Engineering, Georgia Institute of Technology, Atlanta, GA 30332-0250, USA

Received 22 July 2011; Revised 1 November 2011; Accepted 4 November 2011

Academic Editor: Frank Ehlers

Copyright © 2012 Xusheng Sun and Edward J. Coyle. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We consider a set of mobile wireless sensors that collect observations about a brief, localized event. As they continue to move about, one of them processes its observations, decides that an event of interest occurred, and wants to determine if other sensors confirm its results. This sensor assumes the role of a Cluster Head (CH) and requests that all other sensors that collected observations at that time/location reply to it with their decisions. The motion of the sensors since the observation time determines how many wireless hops their decision must cross to reach the CH. We analyze the effect of this motion in the 1D case by modeling each sensor's motion as a Correlated Random Walk (CRW). We also account for measurement errors and communication or processing errors in each wireless hop. Quantities, such as the error probability of the final decision at the CH and the minimum energy required to collect the local decisions from all relevant sensors, can then be directly calculated as functions of time and the parameters of the CRW, the measurement noise and the channel noise. These results allow rapid characterization of the time-dependence of distributed detection algorithms that are executed in mobile sensor networks.

1. Introduction

The systems-level work reported in this paper is motivated by a security scenario that has arisen in the context of the eStadium project [1–3], which develops smartphone applications/games, wireless communication networks, and wireless sensor networks for football games and other large-scale events. One security scenario of great interest occurs in the hour immediately before or immediately after a large event. At these times, large numbers of fans, stadium personnel, and concession staff are typically walking toward or away from the stadium. This is a time when they are very vulnerable to deliberate or accidental exposure to hazardous chemical or biological agents.

Suppose that the smartphones that the fans are carrying, or possibly the tickets that they have bought, have sensors embedded in them that can detect these agents. Capabilities like this are the goal of the US Department of Homeland Security's “Cell-All” program [4]. These sensors would necessarily be small and inexpensive and might thus have high false-alarm rates or low rates of correct detection. It is thus important that as many of their individual decisions as possible be fused with each other to ensure a correct final decision. They must therefore be able to communicate with each other or with a designated Cluster Head (CH) to accomplish this goal. These mobile, battery-powered, wireless sensors can thus be modeled as a mobile ad hoc network that is supporting an application that performs distributed detection.

Our first goal is to develop algorithms to efficiently calculate the performance of distributed detection algorithms in this dynamic and complex scenario. These numerical algorithms should account for as many factors as possible, including the motions that are typical of people in crowds, measurement errors in the sensors, errors made during wireless communication, and fusion algorithms with low enough complexity that they can be executed on very low-power processors or on processors that are shared by many different applications.

This paper is organized as follows. Section 2 describes the scenario that has motivated the work in this paper and that is analyzed as an example in Section 8. Prior work in this area is summarized in Section 3. In Section 4, the transient behavior of semi-infinite correlated random walks is reviewed. The solution to the finite case, which is a good model for the motion of people in the scenario of interest to us, is provided in Section 5. Its statistical properties that have the greatest impact on detection applications are analyzed in Section 6. In Section 7, the performance of detection tasks in mobile wireless sensor networks is determined. Numerical and simulation results are provided in Section 8 for the scenario described in Section 2. They demonstrate (i) the use of the mobility model in calculating the minimum energy required to collect data from each sensor participating in the data fusion process, and (ii) the decision error probability at the cluster head once the decisions from all participating sensors have been gathered. These two quantities are functions of both time and the parameters of the mobility model and we demonstrate how the results in this paper provide the ability to rapidly and accurately calculate or simulate them. Some conclusions and a discussion of future work are provided in Section 9.

2. Scenario: Distributed Detection in a Mobile Ad Hoc Network

To evaluate the tools developed in this paper, we assume a situation like that shown in Figure 1. The wireless nodes shown in the figure are carried by people or vehicles that are moving. They have been moving for some time, generally toward the stadium entrance, which is at position at the right end of each subfigure. Their motion to the right is thus strongly dominant over motion to the left, which is shown by arrows to the right that are larger than arrows pointing to the left. Each node may occasionally switch its direction of motion and briefly move against the general flow, just as people do in a real crowd.

(a) The nodes at the time at which an event occurs. The nodes that observed the event are blue. All nodes are moving, generally toward the stadium entrance at the right,

(b) The nodes at time , the time when some node is the first to process its observations, decide an event occurred and request decisions from other nodes. It floods the network with a message declaring itself to be the CH, organizes the nodes into a multihop cluster, and requests decisions from all nodes with relevant data

(c) The nodes at time . All nodes affected by the event have, on request, completed the processing of their observations and are reporting their local decision back through the cluster to the CH

Figure 1: A security scenario that arises when large crowds are moving toward a stadium or other venue to participate in an event.

At time we assume without loss of generality that the sensors are in a configuration like that shown in Figure 1(a). We assume that the sensors that are colored blue have each collected observations/samples that have been affected by a chemical or biological agent. It takes time for the sensors to process these observations/samples and the processors managing the sensors may have to finish other tasks before processing them. During this processing delay, the people carrying the sensors continue to move toward the stadium.

At time one of the sensors shown in Figure 1(b) has finished processing its observations. The results satisfy some criterion—a preliminary positive detection—that requires that sensor to seek other sensor's results so they can be fused with its own to make a highly reliable decision. This sensor declares itself to be a cluster head (CH) and floods the network of sensors—shown in Figure 1(b) as one sensor transmitting at time —with a request for other nodes to send it their results if they were at approximately the same location and time that the CH acquired its samples. This first request for data will suppress or supersede any other requests by any other sensors for data generated for this location and time. Thus, all other sensors with data about this event will forward it to the CH. For our analysis, we will consider two situations: (i) the CH remains stationary, and (ii) the CH continues to move with the crowd.

At time , the CH's request for reports has reached all other sensors, including those with nothing to report, as shown in Figure 1(c). Sensors with data to report, shown in blue, begin transmitting their results to the CH. In some cases, the CH may be multiple wireless hops away, so intermediate nodes will relay these transmissions. We assume that no local fusion takes place as the nodes' data is relayed to the CH.

It is during this data collection phase that the behavior becomes quite interesting. From the moment it starts, we want to determine the best decision that the CH can make given the multi-hop network that exists at time . If is close to , the sensors would not have spread out very much, so each one may be within one wireless hop of the CH. If is significantly greater than , then the nodes may be either very spread out or, if they have all gotten close to the stadium entrance, have begun to bunch up again. In any of these cases, there will be communication errors, but more errors will occur when the nodes are spread out because additional transmissions are needed to relay the data. The algorithms we use for distributed detection account for these situations and can even determine how much energy each node should dedicate to sending its data—where the energy is measured in terms of the number of bits the node sends.

We develop analytical models of this scenario and numerical methods to calculate the best possible performance that can be expected at any time . We are thus interested in the behavior of the network as it changes over time, where these changes are due to the motion of the crowd. The transient behavior of stochastic models of this motion is critical to understanding this scenario and eventually determining the optimal time at which to make a decision.

We must first model how sensors organize into a cluster, which must be done very rapidly in mobile, event-driven scenarios. We thus assume the existence of a very fast, low-complexity algorithm for clustering that is triggered by an event—in this case the request by one node for other nodes' data. The algorithm we choose is the one defined and analyzed in [5, 6]. We will only consider the single-cluster case that is created in the scenario in this paper; the more general case of hierarchical clustering is relevant if we need to collect results from all sensors participating in an event. In this more general case, many thousands of sensors may be involved, one for each person attending the event.

In the single-cluster scenario, the communication architecture is some variation of the 2-hop cluster shown in Figure 2. The sensor in the center of the cluster, called the CH, is the one that requested the decisions from other sensors. The sensors shown in the first ring are one wireless hop from the CH; the ones in the second ring are two hops away, and so forth. Three sensors in a sector of ring 2 are shown forwarding their decisions to one sensor in ring 1. Each sensor's observation is affected by measurement noise and all communications throughout the network will be affected by channel noise, fading, and transceiver errors. Because the cluster is multi-hop, decisions forwarded from the outer rings will suffer repeated exposure to the sources of communication error and more energy must be expended to get them to the CH.

Figure 2: A 2-hop, 2D cluster. The mobility of the sensors produces this temporary, multi-hop cluster. Each sensor's decision may be incorrect because of measurement noise; transmitted packets may suffer bit errors because of noisy communication channels. If the number of hops in the cluster increases with time, the energy required for communication will increase and detection performance will decrease.

The mobility model governing each sensor and the time since the event of interest will clearly affect the number of hops, and thus both the probability of detection at the CH and the energy consumed to gather the local decisions. It will also affect the network's performance in terms of coverage, maximum throughput, and throughput-delay trade-offs.

3. Prior Work on Mobility Models and Distributed Detection

Many mobility models have been proposed for analyzing the behavior of mobile ad hoc networks like the one described previously. In [7, 8], those models are categorized into four classes: random models, models with temporal dependency, models with spatial dependency, and models with geographic restrictions. For instance, random walks and random waypoint models are random models; Gauss-Markov mobility models and Smooth Random models are models with temporal dependency; group mobility models [7] are models with spatial dependency; Pathway mobility models and Obstacle mobility models are models with geographic restrictions. These models are useful as analytical or simulation tools but none of them can account for all of the constraints or features of real systems.

We thus use Correlated Random Walks (CRWs) as the model for the motion of the sensors in our security scenario. They can account for time dependency, geographical restrictions, and nonzero drift. They are more general than random walks but are still amenable to analysis, sometimes in closed-form. The limiting distribution for a discrete-time 1D doubly-infinite CRW on the integers was derived in [9]. The probabilities of being at any lattice point at the th step for a discrete-time 1D CRW in different cases were found in [10–13]. The absorbing probability and expected duration of a discrete-time CRW was found in [14]. The transient behavior of the continuous-time CRW on a semi-infinite, 1D state space, is solved, in closed-form in some cases, in [15]. The behavior of multiple sensors whose motion is modeled by discrete-time processes related to CRWs has been studied in [16]. In this paper, we use continuous-time, 1D CRWs on finite state spaces to calculate the behavior of mobile sensors in the security scenario described previously.

Ultimately, our goal is to determine—by numerical analysis instead of simulation—the effect of the mobility of the nodes on the detection problem that exists at the application layer of the mobile, ad hoc network. The performance of these applications will vary with time because of the nodes' motions. Understanding the significance of these transient effects, and being able to precisely determine such quantities as a lower bound on the probability of detection, is critical to developing a strategy for reaching a decision meeting certain criteria in the minimum amount of time.

In [17], an optimal distributed detection strategy in a single-hop network was studied. In [18], a decentralized detection problem under bandwidth-constrained communication was investigated. The noise at different sensors was assumed to be independent, but the statistics of the noise were assumed to be unknown to the CH, so it treats all received detection results equally. It is shown in [19] how the performance of detection algorithms is improved by knowledge of the channel. The limits of detection performance in a one-hop sensor cluster with nonideal channels were determined in [20].

In [21], the performance of distributed detection in a random sensor field is analyzed. The works of [22, 23] discuss the distributed detection problem for Gaussian signals under communication constraints. In [24], the authors consider the detection and localization problem of material releases with sparse sensor configurations. In [25], the sensors adopt robust binary quantizers for distributed detection. Censoring and sequential tests have also been adopted in distributed detection to achieve the same detection accuracy with less energy [26, 27], and they have been combined to further save energy [28, 29].

4. Basic Results on Transient Analysis of a Correlated Random Walk on

In this section, we briefly review results on the transient probability distributions of continuous-time, 1D CRWs on [15]. In this case, the sensor moves according to the following rules.

It takes a step in the same direction as its previous step with probability or depending on whether its previous step was in the positive or negative direction, respectively. It takes a step in the opposite direction of its previous step with probability or depending on whether its previous step was in the positive or negative direction, respectively. On reaching a (reflecting) boundary, it takes a step in the opposite direction with probability one.

The time at which the sensor takes its next step is governed by a Poisson process of intensity .

A CRW on with a reflecting boundary at 0 can be modeled as a quasi-birth-death (QBD) process [30] with the state transition diagram shown in Figure 3. The state at any level is entered when the sensor moves to location from location . The state at any level is entered when the sensor moves to location from location .

Figure 3: A Correlated Random Walk (CRW) on .

Let be the row vector of probabilities that the chain is in any of the possible states at time , and define to be the Laplace transform of the vector . Letting denote the generator for this QBD process, then [31]

For this QBD process, For simplicity only, assume that the initial position of the sensor is 0; that is, . Define to be the transform of , the row-vector of transient probabilities for states on level of the process. Then, where Thus we get where

The boundary variables must also satisfy the following: where denotes the right eigenvector corresponding to the eigenvalue of whose magnitude is greater than or equal to 1 for all possible values of . For this CRW, is diagonalizable and its eigenvalues are given by where , for . The right eigenvectors , , corresponding to the eigenvalues and are

We find numerically that, for all possible values of , and , the magnitude of is greater than or equal to 1. Hence, where and , , .

Thus by recursion, where is after the mode whose eigenvalue is outside of the unit circle has been removed [31]. This ensures the stability of the recursion even when numerical errors occur in enforcing the orthogonality condition in (8).

5. New Results on Transient Analysis of Correlated Random Walks on

A CRW on with reflecting boundaries at 0 and can be modeled as a QBD process with the state transition diagram shown in Figure 4 and the following generator: We assume without loss of generality that the initial position of the sensor is 0; that is, . Then we get

Figure 4: A Correlated Random Walk on .

Thus, solving for in is reduced to solving where, if we define to be a identity matrix, and , and the matrix can be computed, via the Cayley-Hamilton theorem, as a linear combination of and . It is then straightforward to calculate the transient distributions of the process via stable algorithms for transform inversion [32–34].

As a check on the aforementioned results for the transient behavior of the CRW, we can use the previous results to find its limiting distribution. Since , it is easy to show that the limiting probability of the CWR on being at any location is [13] For the special case of , , and , for . Note, though, that the transient behavior is of interest to us in the rest of this paper.

Note that this approach to analyzing the transient behavior of a simple CRW can be extended to model more complex motion. The number of states in each level of the QBD—two states in the case of the CRW in Figure 4—can be increased. With four to five states per level, the walker can pause, can walk at different speeds at different times, and so forth. The only difference in the algorithms for calculating the transient distribution is that the eigenvalues and eigenvectors of may not be available in closed form. These techniques can also be extended to some classes of 2D CRWs [15].

6. Statistical Properties of Distances of Nodes from the CH

Suppose that the CH is motionless and there are a total of sensors in the cluster. As described in the introduction, they start at the the same location as the CH () and then move independently based on identical but independent CRW models. Define as the location of the th sensor at time . Further define and as the maximum distance and the minimum distance between any sensor in the cluster and the CH, respectively. Then,

Now assume that the CH is also mobile, that its position is governed by the same CRW as any sensor, and, as before, their motions are all independent. Then and become the maximum distance and the minimum distance between a specified sensor and any other sensor in the cluster. The probabilities remain the same whether the sensor is designated as the CH before or after motion of the sensors begins. Denote the position of the CH by . We have

The probability distribution of the distance between any two sensors is, for , , and, for , Define and as the maximum distance and the minimum distance between any two sensors in the cluster. Then, and for , where sensors are at the location and sensors are at the location , while all the others are in between. The calculation of the probability distribution of is much more complicated and, when , it is always zero.

Thus, for both the motionless- and mobile-CH cases, computing , , , and has been reduced to computing the transient distributions of CRWs. The computational techniques in Sections 2 and 3 thus enable fast and stable computation of statistics of interest in the structure of a cluster of mobile nodes as it evolves over time.

7. Detection in Wireless Sensor Networks

Here we assume a multi-hop sensor network in which each wireless hop is modeled as a Binary Symmetric Channel (BSC). The BSC cross-over probability captures the effect of channel errors on each individual decision bit and is easier to estimate than the full characteristics of the channel. The assumptions of symmetry and the same cross-over probability for each mote-to-mote channel are easily relaxed—they are used only to simplify the analysis so that important trends can be discerned. This paper also assumes that the error probabilities of the individual-received detection results are learned over time by the CH. Hence, the CH is assumed to know the optimal weights for the weighted median in the MAP detector [35] it uses to fuse the received detection results.

In [35], a MAP approach to the distributed detection problem in a multi-hop cluster in a sensor network was developed. It considered a cluster with rings and motes in the th ring. Each mote makes a decision between two hypotheses, and , where “1” denotes that an event has occurred and “0” that it has not. The detection results by different motes are assumed to be i.i.d. Bernoulli random variables, each with a detection error probability . The noise processes in different wireless channels are assumed to be independent and white. Each hop in the network is modeled as a BSC with cross-over probability . The decisions made by motes in the outer rings are relayed by the motes in inner rings to the CH. Let the error probability of the detection results received by the CH from the th ring be . Then, for example, . Denote the detection result received by the CH from the th mote in the th ring by and arrange the detection results in the same ring in a vector . Let denote the event that a decision error happens at the CH.

In a one-hop cluster, suppose that the correct decision should be . Assume the prior probability and find a real number such that . Define . The MAP-based decision bit at the CH is 'th order statistic of . The decision error probability at the CH is It is shown in [35] that, as an estimate of the true decision bit, this weighted order statistic is biased but asymptotically unbiased. When the prior probability is , it simplifies into a majority logic operator.

Now we summarize key results in [35] for the multi-hop case. For a mote that is hops away from the CH, the detection result received by the CH after being relayed over these hops has error probability: We use the notation , which means that should be duplicated times. For simplicity, suppose that the prior probabilities are .

Theorem 1. Define , and assume that these 's can be scaled so that , where the 's are positive integers with . The MAP-based decision bit is then given by . The decision error probability at the CH is given by where is the number of occurrences of in the vector .

Because the weighted median is the MAP detector, we assume—even for sector-level fusion algorithms—that the CH calculates this weighted median after receiving inputs from all motes in the cluster. To understand what happens in these more complex cases, the asymptotic behavior of the weighted median must be determined for the single- and multi-hop cases; otherwise, it is too difficult to determine the probability of a clusterwide decision error at the CH.

The asymptotic behavior of the median filter in one-hop clusters has thus been studied using large deviation techniques in [36]:

We now find the error exponent for the decision error probability for the multi-hop case.

Theorem 2. In multi-hop sensor networks in which is the error probability of the individual detection results received by the CH from nodes in ring , the error probability of the MAP detector is upper bounded by [37]

Proof. In the context of the weighted median filter, the decision error probability at the CH is Since , let , and . We have For each term in the product, . Setting , we find

This error exponent is accurate and the bound is tight when the motes are densely deployed; that is, the number of sensors is large. Minimizing when the number of motes is finite is a very difficult combinatorial problem, so we can minimize this upper bound instead. Also note that the effect of each ring of motes on the decision error probability at the CH is apparent in this bound. It may thus be used to simplify many optimization problems in distributed detection in multi-hop scenarios.

Figure 5 compares the large deviation bound on the error probability with the error probability from simulations for a 3-ring cluster. It shows the desired result that the large deviation error exponent is accurate—the bound is parallel with the simulation result over the entire range of spatial densities. Of course, the bound is high by a multiplicative factor, which is often the case with large deviation techniques.

Figure 5: Comparison for a 3-ring cluster of the error probability obtained via simulation and its large deviations bound. The error exponent is clearly correct and the difference approaches zero asymptotically. Comparing the performance of complex strategies via their error exponents will thus correctly show which strategy is best when motes have large numbers of 1-hop neighbors. Direct calculations and simulations are used to confirm those results for small to moderate numbers of neighbors.

If the CH is static, then we have where

The expected energy consumption for collecting the data is then

If the CH is also mobile and its position is governed by a CRW, then we have where The expected energy consumption for collection of the data is

8. Numerical Results

Figure 6 shows the transient probability distribution of two CRWs on with different parameters. One CRW is asymmetric, which means that the walker tends, in the case shown, to prefer motion to the right. The other CRW is symmetric, which means that the walker has no preference for one direction or the other—but its motion is still correlated, unlike in a standard random walk.

(a) An asymmetric CRW with and . The walker starts at 0 and then drifts toward and tends to stay near 10

(b) A symmetric CRW with and . The walker starts at 0 and then it is eventually equally likely to be in each position except a boundary state

Figure 6: Transient probability distribution calculated for two CRWs on . The curve for a given shows the probability that the walker is at position at time .

Figure 7 shows the transient probability distribution of the distance from a static clusterhead to the furthest of five mobile sensors for each point in time. The static clusterhead is assumed to be at location 0. This distance is important for calculating the number of wireless hops between each sensor and the CH that is trying to gather data. If the transmission radius of each sensor is , then the number of hops to the CH is the ceiling of .

(a) All sensors move according to asymmetric CRWs with and

(b) All sensors move according to symmetric CRWs with and

Figure 7: Calculations of the transient probability distribution of the distance from a static CH to the furthest of five sensors that are moving according to independent CRWs on .

Figure 8 shows the transient probability distribution of the maximum distance between any of the five sensors and the CH when the mobile CH moves according to the same mobility model as the sensors. This is probably the most realistic case for mobile networks. It leads to an interesting phenomenon when the motion is on a finite grid and the CRWs are all asymmetric in the same direction. The sensors start at 0 when they gathered their observations. With time, they spread out but they all tend to move to the right. As they get close to the other barrier, they then tend to bunch up again. This effect is similar to a crowd moving in a given direction and then gathering around that destination. The distance here, when divided by the transmission radius of the sensors, tells the maximum size, in wireless hops, of the cluster.

(a) All sensors move according to asymmetric CRWs with and

(b) All sensors move according to symmetric CRWs with and

Figure 8: Calculations of the transient probability distribution of the distance from a mobile CH to the furthest of five sensors for from 0 to 10. The CH and the sensors are all moving independently and according to CRWs on .

Figure 9 shows the expected energy consumed—as a function of the time at which data collection starts—for a static or mobile CH to collect one packet of data from five mobile sensors. The six sensors all collected measurements at the same time/place but then continued to move. By the time a request to send in all data has been received, they are at locations in the state space with probabilities determined by their CRWs. The data collected and decisions made by the sensors that have wandered the furthest require the most energy to collect—that data must travel over the largest number of hops to reach the CH.

(a) All sensors move according to asymmetric CRWs with and

(b) All sensors move according to symmetric CRWs with and

Figure 9: Calculation of the expected energy consumed by the network when the collection of data from the five mobile sensors begins at time . The cases considered are when the transmission radius of each sensor is or and all sensors move according to symmetric or asymmetric CRWs.

Note that some curves in Figure 9 have a stepped appearance. This is due to the existence of the rings in the network that result from the fixed transmission radius for each node. When a node in ring moves further from the cluster head, it may cross from ring to ring . Its packets must then be relayed over one additional hop to reach the CH. This extra relay results in a significant jump in the energy that must be expended by the network.

Also note that the mobility of the CH has a significant impact on the energy consumed to collect data. If it is mobile and all sensors are moving according to a CRW that is asymmetric, then the energy required for collection first increases and then decreases. In this case, it is better to collect data either very quickly after the event, or, if that is not possible because of the time to process measurements, to wait until the CH and the other sensors bunch up again at their destination.

Figure 10 shows the large deviation bound of the error probability as a function of time for a static or mobile CH to collect one packet of data from five mobile sensors. It shows good agreement with the results in Figure 9, which shows the expected energy consumption as a function of time for a static or mobile CH to collect one packet of data from five mobile sensors.

(a) All sensors move according to asymmetric CRWs with and

(b) All sensors move according to symmetric CRWs with and

Figure 10: Calculation of the large deviation bound of the error probability when the collection of data from the five mobile sensors begins at time . The cases considered are when the transmission radius of each sensor is or with and and all sensors move according to symmetric or asymmetric CRWs.

Figure 11 shows the expected error probability as a function of time for a static or mobile CH to collect one packet of data from five mobile sensors. The same comments made previously about the effects of the mobility of the CH, and the effects of symmetric and asymmetric CRWs, apply here. The values of the probabilities of the nodes locations used to calculate these expected error probabilities were obtained in this case by simulations of the CRW.

(a) All sensors move according to asymmetric CRWs with and

(b) All sensors move according to symmetric CRWs with and

Figure 11: Expected error probability when the collection of data from the five mobile sensors begins at time . The cases considered are when the transmission radius of each sensor is or with and and all sensors move according to symmetric or asymmetric CRWs. The curves in this figure were generated from 5 million runs of a simulation of the detection process for each setting of the mobility parameters. The values of the probabilities of the nodes locations for the detection process were obtained by simulations of the CRW.

9. Conclusions

In this paper, the solution to the finite state space CRW was provided and its statistical behavior was studied both analytically and numerically. As an illustration of the application-specific performance measures it can help address, we studied the impact of motion on the error probability of the final decision at the CH and the energy required to collect the decisions from all relevant sensors. Thus the temporary cluster head can decide the most appropriate time to call for data reports to help it make decisions within the allowed time frame and energy budget. In future research, we will study the impact of motion on the energy required to gather data and the error probability of the final decision in complex mobile wireless sensor networks.

References

X. Zhong, H. H. Chan, T. J. Rogers, C. P. Rosenberg, and E. J. Coyle, “The development and eStadium testbeds for research and development of wireless services for large-scale sports venues,” in Proceedings of the 2nd International Conference on Testbeds and Research Infrastructures for the Development of Networks and Communities (TRIDENTCOM '06), pp. 340–348, Barcelona, Spain, March 2006.

X. Zhong and E. J. Coyle, “eStadium: a wireless “living lab” for safety and infotainment applications,” in Proceedings of the 1st International Conference on Communications and Networking in China (ChinaCom '06), pp. 1–6, Beijing, China, October 2006.

A. Ault, J. V. Krogmeier, S. R. Dunlop, and E. J. Coyle, “Stadium: the mobile wireless football experience,” in Proceedings of the 3rd International Conference on Internet and Web Applications and Services (ICIW '08), Athens, Greece, June 2008.

“US DHS cell-all program,” http://www.dhs.gov/files/programs/gc_1268073038372.shtm

S. Bandyopadhyay, Q. Tian, and E. J. Coyle, “Spatio-temporal sampling rates and enery efficiency in wireless sensor networks,” IEEE/ACM Transactions on Networking, vol. 13, no. 6, pp. 1339–1352, 2005.

S. Bandyopadhyay and E. J. Coyle, “Minimizing communication costs in hierarchically-clustered networks of wireless sensors,” Computer Networks, vol. 44, no. 1, pp. 1–16, 2004.

F. Bai and A. Helmy, “A survey of mobility modeling and analysis in wireless ad hoc networks,” in Wireless Ad Hoc and Sensor Networks, Kluwer Academic Publishers, Boston, Mass, USA, 2004.

T. Camp, J. Boleng, and V. Davies, “A survey of mobility models for ad hoc network research,” Wireless Communications and Mobile Computing, vol. 2, no. 5, pp. 483–502, 2002.

S. Goldstein, “On diffusion by discontinuous movements, and on the telegraph equation,” Quarterly Journal of Mechanics and Applied Mathematics, vol. 4, no. 2, pp. 129–156, 1951.

E. Renshawn and R. Henderson, “The correlated random walk,” Journal of Applied Probability, vol. 18, pp. 403–414, 1981.

W. Böhm, “The correlated random walk with boundaries: a combinatorial solution,” Journal of Applied Probability, vol. 37, no. 2, pp. 470–479, 2000.

A. D. Proudfoot and D. G. Lampard, “A random walk problem with correlation,” Journal of Applied Probability, vol. 9, pp. 436–440, 1972.

R. Lal and U.N. Bhat, “Some explicit results for correlated random walks,” Journal of Applied Probability, vol. 27, pp. 757–766, 1989.

Y. L. Zhang, “Some problems on a one-dimensonal correlated random walk with various types of barriers,” Journal of Applied Probability, vol. 29, pp. 196–201, 1992.

S. Bandyopadhyay, E. J. Coyle, and T. Falck, “Stochastic properties of mobility models in mobile ad hoc networks,” IEEE Transactions on Mobile Computing, vol. 6, no. 11, pp. 1218–1229, 2007.

H. Cai and D. Y. Eun, “Toward stochastic anatomy of inter-meeting time distribution under general mobility models,” in Proceedings of the 9th ACM International Symposium on Mobile Ad Hoc Networking and Computing (MobiHoc '08), pp. 273–282, Hong Kong, China, May 2008.

K. Liu and A.M. Sayed, “Optimal distributed detection strategies for wireless sensor networks,” in Proceedings of the 42nd Annual Allerton Conference on Communication, Control and Computing, pp. 1651–1661, Monticello, Ill, USA, September 2004.

Z. Q. Luo, “Universal decentralized detection in a bandwidth constrained sensor network,” IEEE Transactions on Signal Processing, vol. 53, no. 8, pp. 2617–2624, 2005.

B. Chen and P. K. Willett, “On the optimality of the likelihood-ratio test for local sensor decision rules in the presence of nonideal channels,” IEEE Transactions on Information Theory, vol. 51, no. 2, pp. 693–699, 2005.

B. Chen, L. Tong, and P. K. Varshney, “Channel-aware distributed detection in wireless sensor networks,” IEEE Signal Processing Magazine, vol. 23, no. 4, pp. 16–26, 2006.

R. Niu and P. K. Varshney, “Performance analysis of distributed detection in a random sensor field,” IEEE Transactions on Signal Processing, vol. 56, no. 1, pp. 339–349, 2008.

W. Li and H. Dai, “Distributed detection in wireless sensor networks using a multiple access channel,” IEEE Transactions on Signal Processing, vol. 55, no. 3, pp. 822–833, 2007.

S. K. Jayaweera, “Bayesian fusion performance and system optimization for distributed stochastic Gaussian signal detection under communication constraints,” IEEE Transactions on Signal Processing, vol. 55, no. 4, pp. 1238–1250, 2007.

E. B. Fox, J. W. Fisher, and A. S. Willsky, “Detection and localization of material releases with sparse sensor configurations,” IEEE Transactions on Signal Processing, vol. 55, no. 5, pp. 2172–2181, 2007.

Y. Lin, B. Chen, and B. Suter, “Robust binary quantizers for distributed detection,” IEEE Transactions on Wireless Communications, vol. 6, no. 6, pp. 2172–2181, 2007.

S. Appadwedula, V. V. Veeravalli, and D. L. Jones, “Decentralized detection with censoring sensors,” IEEE Transactions on Signal Processing, vol. 56, no. 4, pp. 1362–1373, 2008.

Y. Mei, “Asymptotic optimality theory for decentralized sequential hypothesis testing in sensor networks,” IEEE Transactions on Information Theory, vol. 54, no. 5, pp. 2072–2089, 2008.

P. Addesso, S. Marano, and V. Matta, “Sequential sampling in sensor networks for detection with censoring nodes,” IEEE Transactions on Signal Processing, vol. 55, no. 11, pp. 5497–5505, 2007.

Y. Yao, “Group-ordered SPRT for distributed detection,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '08), pp. 2525–2528, Las Vegas, Nev, USA, April 2008.

M. F. Neuts, Matrix Geometric Solutions in Stochastic Models, Johns Hopkins University Press, London, UK, 1981.

J. Zhang and E. J. Coyle, “Transient analysis of quasi-birth-death processes,” Communications in Statistics, vol. 5, pp. 459–496, 1989.

J. Abate, G. L. Choudhury, and W. Whitt, “An introduction to numerical transform inversion and its application to probability models,” in Computational Probability, W. Grassman, Ed., pp. 257–323, Kluwer, Boston, Mass, USA, 1999.

J. Abate and W. Whitt, “The Fourier-series method for inverting transforms of probability distributions,” Queueing Systems, vol. 10, no. 1-2, pp. 5–87, 1992.

J. Abate and P. P. Valko, “Multi-precision Laplace transform inversion,” International Journal for Numerical Methods in Engineering, vol. 60, no. 5, pp. 979–993, 2004.

Q. Tian and E. J. Coyle, “Optimal distributed detection in clustered wireless sensor networks,” IEEE Transactions on Signal Processing, vol. 55, no. 7, pp. 3892–3904, 2007.

J. F. Chamberland and V. V. Veeravalli, “Asymptotic results for decentralized detection in power constrained wireless sensor networks,” IEEE Journal on Selected Areas in Communications, vol. 22, no. 6, pp. 1007–1015, 2004.

X. Sun and E. J. Coyle, “Low-complexity algorithms for event detection in wireless sensor networks,” IEEE Journal on Selected Areas in Communications, vol. 28, no. 7, Article ID 5555912, pp. 1138–1148, 2010.