Interference in Poisson Networks with Isotropically Distributed Nodes
Abstract
Practical wireless networks are finite, and hence nonstationary with nodes typically nonhomogeneously deployed over the area. This leads to a locationdependent performance and to boundary effects which are both often neglected in network modeling. In this work, interference in networks with nodes distributed according to an isotropic but not necessarily stationary Poisson point process (PPP) are studied. The resulting link performance is precisely characterized as a function of (i) an arbitrary receiver location and of (ii) an arbitrary isotropic shape of the spatial distribution. Closedform expressions for the first moment and the Laplace transform of the interference are derived for the path loss exponents and , and simple bounds are derived for other cases. The developed model is applied to practical problems in network analysis: for instance, the accuracy loss due to neglecting border effects is shown to be undesirably high within transition regions of certain deployment scenarios. Using a throughput metric not relying on the stationarity of the spatial node distribution, the spatial throughput locally around a given node is characterized.
I Introduction
Stochastic geometry, in particular the theory of point processes, has recently attracted much attention in the field of interference modeling and performance analysis for wireless networks with many uncertainties such as mobile/dynamic user locations and channel fading. In a nutshell, the locations of the nodes are modeled as a realization of a stochastic point process rather than assuming a fixed spatial configuration. Since the emitted signals undergo a distantdependent path loss, the interference experienced by a given node becomes random. Its statistical properties, moreover, depend on several factors including the law of the spatial distribution of nodes.
Ia Related Work and Motivation
Interference modeling and network analysis using tools from stochastic geometry have become a multifaceted research field [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. Early works on interference modeling assumed a stationary PPP for the interferer locations, cf. [1, 2]. Using a similar model, the spatial throughput of decentralized networks with Aloha medium access control (MAC) was analyzed in [3, 4]. Following these works, the node locations have mostly been modeled by stationary point processes which typically leads to locationindependent statistical properties of the considered performance quantities, e.g., interference, outage probability or throughput. Among these advances, the transmission capacity framework [5, 6, 7, 9, 8] substantially contributed to a better understanding of the interactions between the basic system parameters of a wireless network. Besides, stationary models with nonhomogeneous node deployments, e.g., PoissonCluster [10] and Matérn hardcore models [11, 12], were also investigated as such models are wellsuited for studying more sophisticated medium access control (MAC) schemes. Treated as general motioninvariant, these models were further analyzed in [15, 13, 14] in a unifying way.
Stationarity of the spatial distribution of nodes is a desirable property since it allows for analytic tractability and, even more important, it represents a key requirement for applying certain performance metrics such as the transmission capacity metric [5]. In practice though, wireless networks exhibit a nonstationary spatial node distribution; for instance, because of a finite network area with boundary regions. Consequently, performancerelevant quantities such as the experienced interference typically vary across the network area, thereby complicating modeling and system design. Besides this simple example, more complex deployments are often found in practice, e.g., wireless sensor networks created by airdrop [21] or spontaneous formation of hot spots [22]. The spatial configuration of such hotspots is typically dictated by user motion and by geometric constraints as illustrated in [23] for the example of a campuswide WiFi deployment. Furthermore, there is a growing need for cellular operators to better understand not only the temporal variations in user traffic demands but also its spatial dependence; for instance, the optimal interplay between smallcell deployments and WiFi offloading–a promising approach to boost capacity in dense areas–requires carefully pinpointing areas of peaktraffic demands [24] and identifying locations for deployment so to reach the mobile users [25]. Hence, analytic tools for quantifying the network performance while taking into account user mobility as well as hardtopredict spatial configurations are of eminent importance.
The need for nonstationary models for characterizing more complex node deployments was reported for instance in [26]. The authors discussed techniques that generate nonuniform node distributions for the purpose of efficient network simulations. In [16], a nonstationary and isotropic node distribution was assumed for analyzing multiantenna receivers in the presence of interference. The analysis showed that the shape of the spatial distribution has a considerable impact on link performance. In [27], a first attempt was made towards analyzing the link performance at an arbitrary receiver location and for an arbitrary isotropic node distribution.
The shortcomings associated with the stationarity assumption are summarized below:
Infiniteness and boundary effects: Stationarity implies that the network consists of infinitely many nodes spread over an infinitely large region. However, the number of nodes as well as the network area is finite in practice. Boundary effects are ignored although they play a critical role in realworld networks because of unequal performance among the nodes in terms of, e.g., local topology, interference/noiselimited performance, etc. Ensuring the qualityofservice (QoS) level targeted before deployment hence becomes difficult.
Model artifacts for path loss exponent : In the stationary case and with a path loss exponent of , the sum interference is almost surely (a.s.) infinite [5]. More specifically, stationary models lose their accuracy as the path loss exponent tends to since interference then becomes dominated by the infinite number of far nodes.
Lack of local throughput metric: Nonstationary models prevent the use of certain throughput metrics such as the transmission capacity. This is because it is no longer possible to infer the global performance from the local analysis as the local performance is locationdependent in nonstationary deployments.
IB Contributions and Outcomes
We extend prior work by modeling the node locations by an isotropic PPP, with stationarity being a special case. The network model is explained in detail in Section II. The contributions are summarized below.
Interference characterization: The first moment and distribution function are studied for both arbitrary and Rayleigh fading channels. In the case of arbitrary fading, we derive in closedform the exact first moment as a function of the spatial shape of the node distribution and of the arbitrary receiver location for the path loss exponents and in Section IIIA. Using these results, an upper bound on the tail probability of the interference is derived. A corresponding lower bound that is not limited to the above values of is also presented for suitable spatial shapes. In Section IIIB, we derive the Laplace transform of the interference for the Rayleigh fading model. We also address the case , which eluded a meaningful analysis in the stationary model. One important insight is that for , one can find situations in which a.s. infinitely many nodes contribute to the interference while the interference remains finite a.s. This result sharpens prior statements about the nature of the interference for and suggests that there exists a transition between sparse and dense networks.
Outage probability and model accuracy: The locationdependent outage probability is characterized in Section IVA as a function of the spatial shape. We demonstrate the use of the developed model in Section IVB by comparing it to a stationary model that uses a local approximation capturing the nonhomogeneity in the spatial nodes distribution. We show that, depending on the spatial shape in question, large approximation errors are avoided by the developed model, particularly in transition areas where boundary effects come into play.
Applications: We propose a metric that is capable of quantifying the local throughput in nonstationary networks in Section V and demonstrate through an example how this metric can be used. In this example, we refine an existing result on code division multiple access (CDMA) systems in decentralized networks for the case . We also show how the model accuracy for carriersense based MACs in decentralized networks can be significantly increased by means of the developed model.
Notation: Sansserifstyle letters () represent random variables while serifstyle letters () represent deterministic variables or constants. The imaginary unit is given by and denotes the real part of a complexvalued number. We define by and a disc and a circle, respectively, centered at with radius . The origin is denoted by .
Ii Mathematical Model
We consider a packetbased wireless network with identicallyequipped nodes isotropically distributed in . The nodes are assumed to be slotsynchronized. In a randomly chosen slot, some nodes wish to transmit a packet. We assume that the locations of these transmitters follow an isotropic PPP on with intensity (equivalently, density) being defined on . Throughout this work we will denote by to the random location of the th node as well as to the th node itself. Due to isotropy of , is rotationinvariant [28] and depends solely on the distance to the origin, i.e., for all .^{1}^{1}1In a very few situations, we will switch between Cartesian coordinates () and the corresponding polar coordinates () when appropriate. We do not expect any confusion thereof. For notational convenience, we define . The next definition is a consequence of the fact that can be described as the resulting intensity after locationdependent thinning of a stationary PPP of some constant intensity [11].
Definition 1.
The shape function , reflecting the spatial shape of , is defined on by the relation , where is some intensity scaling constant.
We assume that each transmitter has an intended receiver located at fixed distance . The fixed distance assumption, which can be seen as a target transmission distance dictated by the network protocol, is commonly accepted, see for example [5]. In order to measure the (spatiallyaveraged) link performance in the network we define a reference link, cf. [7]; the reference link consists of a reference receiver placed at an arbitrary location and of an associated reference transmitter placed at , where lies somewhere on the circle . Note that neither the receiver nor the transmitter are part of the point process . By the SlivnyakMecke Theorem [28], the statistics of are not affected by the addition of the reference link.
We consider a path loss plus blockfading channel with independent and identically distributed (i.i.d.) fading coefficients. The power path loss between two locations is given by the path loss function with path loss exponent , where ensures boundedness of . The power fading coefficient between a transmitter at and the reference receiver is given by , where for all . When appropriate we will drop the index in .
We assume that all nodes transmit with the same fixed transmit power and at a common information rate. The sum interference power at the reference receiver is then given by
(1) 
Treating interference as white noise, the instantaneous signaltointerferenceplusnoise ratio at the reference receiver is given by
(2) 
where is the average signaltonoise ratio. We assume that the nodes employ strong channel coding such that the outage event is a steep function of the . Focusing on the case where explicit transmitter coordination as well as CSI feedback is not possible, the outage probability is a useful metric for characterizing the link performance.
Definition 2.
The outage probability of the reference link is given by
(3) 
where is a modulation and coding specific threshold.
Iia Spatial shapes chosen for illustrations
We next introduce four exemplary spatial shapes used for illustrations and numerical evaluations. They are chosen such to roughly characterize typical scenarios in wireless networks and to help increasing the reader’s intuition about the results. An exact validation of the chosen spatial models through comparison with realworld deployments is outside the scope of this paper. The spatial shapes are depicted in Fig. 2.

Scenario a), finite network: This scenario reflects the basic property all practical networks share: the network area is finite, or equivalently, the node density tends to zero for sufficiently large distances to the network center. It is assumed that the density first remains constant over a wide range as in the stationary case. At the network boundary the density then starts to decay rapidly until it becomes zero.

Scenario b), urban with hotspot: In urban scenarios there may sometimes exist small areas with very high data traffic, i.e., communications hotspots. They are typically found in commercial areas or other public places, comprising many densely–and sometimes dynamically–deployed wireless architectures [29]. We model such a hotspot scenario by “adding a plateau of density” around the origin on top of an already existing level of density corresponding to the urban deployment. This level then decays to zero with increasing distance to the origin to reflect finiteness of the network.

Scenario c), scattered decentralized network: There are certain types of applications that preclude a detailed network layout for the reason of hostile environments or limited geographic access. For instance, large sensor networks are sometimes created by airdrop which results in a highly scattered spatial distribution of devices. We model such a behavior by an exponentially decreasing density around the origin.

Scenario d), carrier sensing in decentralized networks: This scenario is found in decentralized networks with transmitters employing carrier sensing to avoid excessive interference; suppose a transmitter (here located in the origin) is granted access to the medium. Consequently, other potential transmitters directly surrounding this transmitter defer their transmission as they sense the medium as busy, while other potential transmitters farther away sense the medium as free and therefore start to transmit, cf. Section VB for more details.
Iii Interference Analysis
We now study the interference statistics at the reference receiver . The analysis first focuses on the case of an arbitrary fading distribution. We derive the first moment of the interference and then use bounding techniques such as those used in [5] to characterize the interference distribution. For the Rayleigh model, we then derive the Laplace transform of the interference distribution.
Iiia Arbitrary Fading Model
IiiA1 First Moment of the Interference
The first moment of the interference measured at the reference receiver can in general be obtained by
(4) 
Obviously, one would have to know the distribution of , which unfortunately is known in closedform only for a few cases of stationary point processes [13]. A remedy to this problem is given by the Campbell Theorem [28], which allows us to derive the first result:
Theorem 1.
Let , , and . If for some , then
(5) 
where is given by
(6) 
A proof is given in Appendix B. The condition for some is necessary for to exist. The function in (6) has an interesting interpretation: it can be seen as the interferencedriving function as it determines the interference up to a scaling factor. Additionally, the first term in (6) can be interpreted as the interference field associated with the origin , while the second term effectively adds up the interference according to .
Remark 1.
If the reference receiver is located in the origin (), the asinhterm in has to be replaced by , cf. Identity 2 in Appendix A.
Glancing at the second term in (6), we note:
Corollary 1.
When and for all , can be interpreted as a complementary cumulative distribution function (CDF) with respect to a random distance to the origin, yielding
(7) 
Corollary 1 states that the integral in (6) can be seen as an averaging of the interference with respect to a random distance . Such a representation may be appropriate when analyzing networks with a priori unknown or fastvarying spatial configurations, for which a CDF is then used to model their spatial shape.
Corollary 2.
Let for some . Then, the expected number of interferers but the expected interference .
The intuition behind Corollary 2 is that, although the expected number of nodes contributing to the interference is unbounded, the network remains sufficiently sparse such that the first moment of the interference remains bounded. Note that for a PPP, if the expected number of interferers is unbounded, this implies that the number of interferers is a.s. infinite which can be verified by studying the Laplace transform of the PPP [28, 30]. Applying the Markov Inequality [31], this in turn means that the number of nodes is a.s. infinite while the interference remains a.s. finite. This particular finding is somewhat remarkable since it rearranges the commonlyaccepted perception, stating that whenever and the number of interferers is a.s. infinite, the interference is a.s. infinite as well [13, 7]. This perception indeed holds for stationary point processes but does not hold in general for nonstationary point processes, as demonstrated by Corollary 2.
Theorem 2.
Let for . Then, a.s.
A proof is given in Appendix B. Theorem 2 shows that whenever decays at most logarithmically, the interference is a.s. infinite. In particular, this includes the stationary case (since ) which is consistent with the literature [13]. Combining Corollary 2 and Theorem 2, we observe that for asymptotically decaying there exists a transition between sparse and dense networks. This transition determines whether or not the interference is a.s. finite in a nonstationary Poisson network with a.s. infinite number of interferers.
Remark 2.
By setting , the pure path loss case is also covered by the above results.
We now characterize the first moment of the interference for the case .
A proof is given in Appendix B. The term in (9) can be again interpreted as the interferencedriving function.
Corollary 3.
Let so that (stationary PPP). Then, by carefully taking the limits
the wellknown result for the stationary PPP is recovered [13].
IiiA2 Bounds on the Interference Distribution
We next treat the problem of bounding the tail of . A simple upper bound can be obtained using Markov’s inequality in combination with Theorem 1 and Theorem 3. For the construction of a lower bound, we first recall the definition of subharmonic functions.
Definition 3.
(Subharmonic functions [32, Ch. 2]): Let be an open set and let be a function twice continuously differentiable on . If on , then is called subharmonic on .
If is convex in a certain (onedimensional) region, then the intensity is subharmonic on the corresponding (twodimensional) region. Such a behavior may be often found at the network boundary, e.g., Scenario a) and b), or when the shape function exhibits a tail, e.g., Scenario c) and d). In this case a lower bound on the tail probability can be derived:
Theorem 4.
Let be subharmonic on and let . Denote by the maximum radius for which the closed ball is contained in , i.e., . Then,
(10) 
A proof is given in Appendix B. Note that subharmonicity includes the case of harmononicity. The construction of the lower bound in Theorem 4 basically builds on the socalled “dominant interferer” phenomenon introduced in [7], where it was also reported that the resulting bound is fairly tight. However, in our case the tightness of (10) strongly depends on the second derivative of and may not be guaranteed.
Corollary 4.
Let . For the pure path loss model (), (10) reduces to
(11) 
The restriction is necessary to allow the closest interferer to be dominant, otherwise this bounding technique would yield the trivial lower bound .
IiiB Rayleigh Fading Model
For the commonlyused Rayleigh fading model, it was demonstrated in [4] that the Laplace transform of the interference is useful for computing outage probabilities. We therefore derive the Laplace transform of , i.e., next for arbitrary and . Similar to Section IIIA, we focus again on and .
Theorem 5.
A proof is given in Appendix B. Note that for and for every , we have that for all . This in turn implies a.s. which is consistent with Theorem 2.
Remark 3.
Setting for all and , we recover the wellknown result for the homogeneous case with : .
Remark 4.
The case with Rayleigh fading may seem contradictory first: A path loss exponent equal to is typically observed in propagation environments without groundplane reflection [33]. In contrast, Rayleigh fading models the nonlineofsight (NLOS) case with many reflected paths impinging at the receiver. It turns out, however, that there may truly exist urban NLOS scenarios with considerably small path loss exponents () as demonstrated in [34]. Depending on the geometry of objects in the proximity of the receiver, the received signal may therefore still undergo severe smallscale fading. Hence, the case may serve as a theoretical limit of what can be expected roughly in Rayleigh fading environments with a small path loss exponent.
Iv Outage Probability and Model Accuracy
In this section, the outage probability is characterized for the underlying setting and the model accuracy is studied.
Iva Outage Probability
For arbitrary fading the bounds on the tail probability of can be used to bound the outage probability, see for example [7]. In the sequel, we will focus on the Rayleigh fading case and discuss the impact of the spatial shape on the resulting performance.
Corollary 5.
The outage probability at in the Rayleigh fading model is
(14) 
for the cases and , where is given by (13).
Corollary 5 now allows for measuring the exact outage probability at an arbitrary location and for an arbitrary spatial shape function . Fig. 3 shows vs. for and . The spatial shape function was chosen according to scenario a). In the noisefree case () the outage probability decreases monotonically with increasing distance to the origin. Furthermore, we observe that for the outage probability is higher and its slope is less steep than it is for . This is because for the individual interference contributions decay more slowly over distance than they do for . As a result, the interference is no longer dominated by only a few nearby interferers but it is determined by the large number of nodes nodes, including those relatively far away from . When receiver noise is considered ( dB) the behavior of changes considerably: while is on the same order as in the noisefree case around the center of the network, both curves converge to a constant outage probability level as the reference receiver eaves the center of the network. In fact, in this boundary region outage is primarily due to bad fades rather than to interference, thus rendering the performance noiselimited rather than interferencelimited. This transition–from the interferencelimited to the noiselimited regime–can be precisely tracked owing to the developed model; for instance, Fig. 2 suggests that the noiselimited regime commences somewhere around , while Fig. 3 reveals that this is not true at least for ( for the reason explained above.
IvB Exact vs. Approximate Model
Up to this point, it is not yet very clear how much can be gained by the interference model derived in this work. In order to quantify the gains, we compare our model to a simpler one that approximates the nonstationarityproperty of the interference locally around the reference receiver . In this simpler model, the interference field at is assumed to originate from a stationary PPP having constant intensity ; in other words, the networkwide spatial distribution of interferers is approximated locally by the density at location .
For the stationary case the outage probability for Rayleigh fading is wellknown [3, 4]. With the above approximation the approximate outage probability at is then given by
(15) 
The intuition behind (15) is that the exact outage probability is approximated by taking the outage probability expression corresponding to the stationary case and modulating the intensity by the spatial shape function at . To measure the difference between and we define the following metric.
Definition 4 (Logdivergence).
The logdivergence is defined as
(16) 
The logdivergence quantifies the ratio of the exact and approximate success probability () on the logarithmic scale for arbitrary . The normalization by is necessary to remove the dependency on so to measure the divergence resulting from the spatial shape only. For large positive the approximation overestimates the true outage probability, while for large negative outage probability is underestimated. The approximation is accurate whenever is small. The logdivergence can be computed using Corollary 5.
Corollary 6.
In the Rayleigh fading model and for , the logdivergence for the case is
(17) 
Fig. (a)a shows the logdivergence vs. for the shape functions introduced in Section IIA. It can be seen that the logdivergence exhibits an oscillatory behavior that depends on the degree of variability of . For example, in scenario a) and c) changes comparably slow, and thus the corresponding logdivergence shows only weak oscillations. The logdivergence for scenario a) remains low () over a wide range which suggests that the local approximation works well in this case. As for scenario c), however, the logdivergence is large around the origin which is due to the fact that the outage probability is highly overestimated as the (exponential) decay of to the righthand side is neglected by the approximation. A similar effect can be observed for scenario d): the outage probability is highly underestimated around the origin because the increasing density to the righthand side is neglected. The logdivergence for scenario b) exhibits rich oscillations due to a stronger varying shape function. As can be seen, these oscillations are high particularly in the transition region ().
To understand the impact of the logdivergence on outage probability, the relative error
(18) 
is shown in Fig. (b)b for the example . Fig. (b)b underlines the observations made in Fig. (a)a: the approximation works well for scenario a) while for the other three scenarios the approximation is relatively loose. Especially for the scenarios b) and d), the relative error is considerably high over a wide range (between 1% and 10% around the network center for scenario c)). For scenario b) the relative error is approximately 10% around the density midlevel. In case of scenario d) the approximation completely fails around the origin.
V Applications and Examples
In this section, the developed model is applied to problems in network modeling with nonstationary spatial node distributions.
Va Local Transmission Capacity
As argued in Section I, the transmission capacity metric cannot be applied to networks with nonstationary spatial node distributions. Based on the developed model, the definition of the transmission capacity can however be extended to account for nonstationarity and location dependency.
Definition 5 (Local transmission capacity).
The local transmission capacity is defined as
(19) 
and gives the maximal density of concurrent transmissions in an infinitesimal region around location subject to an outage probability constraint .
Since the local transmission capacity accounts for the spatial shape, this metric allows throughput bottlenecks to be spatially tracked and properly engineered, e.g., by balancing QoS among nodes irrespective of their location. For isotropic node distributions, in particular, the local transmission capacity depends on and can be computed/bounded by algebraic manipulations of the outage probability expressions derived in the previous sections.
Example: FHCDMA vs. DSCDMA in decentralized networks: In [7] it was shown that in the stationary PPP model, the transmission capacity gain of frequencyhopping (FH)CDMA compared to directsequence (DS)CDMA scales as , where is the processing gain. For very small path loss exponents () this result suggests that the gains of FHCDMA vanish irrespective of the processing gain . This observation, however, results from the stationarity assumption not being able to correctly capture the case . The next result rearranges this scaling result for the case for Rayleigh fading and the reference receiver located in the origin.
Corollary 7.
Let , and . In the Rayleigh fading model, the local transmission capacity gain of FHCDMA over DSCDMA at is
(20) 
VB Interference in Networks with TransmitterInhibition
In order to study decentralized networks with inhibition mechanisms such as carriersense medium access (CSMA) or local frequency division multiple access (FDMA) while ensuring analytic tractability, methods based on nonhomogeneous Poisson approximation have been proposed for stationary models [18, 12, 19, 20]. When such protocols are transmitterinitiated, e.g., transmitter sensing for CSMA, the resulting spatial distribution of interferers becomes nonhomogeneous but remains isotropic around the inhibiting transmitter since potential transmitters around the inhibiting transmitter are kept silent while others located farther away are likely to transmit. In contrast, the interference field at the associated receiver is not isotropic. To overcome this intractability, the receiver is assumed to be colocated with the inhibiting transmitter, thereby virtually rendering the interference field around the receiver isotropic at the cost of losing model accuracy. This loss depends on the distance between the inhibiting transmitter and the associated receiver. Using the developed model, we can now evaluate the accuracy loss resulting from assuming that transmitter and receiver are colocated for this CSMA modeling technique. We briefly summarize the basic ideas of this modeling technique and refer to [18, 19] for further details.
Nonhomogenous Poisson approximation for CSMA networks: Let the potential interferers be initially distributed according to a stationary PPP of density . Assume the that the reference transmitter and the reference receiver are located at and , respectively, and separated by . The inhibition mechanism is modeled in three steps: 1) The largescale density of active interferers is derived using a Matérntype 2 model [28], which captures the inhibition effect on a “macroscopic” level. The largescale density is then given by
(21) 
where is the sensing threshold. We then condition on the fact that is granted access to the channel and we are now seeking the statistical characterization of the interferers around after this conditioning. 2) This is where the nonhomogeneous Poisson approximation comes into play, yielding a “smallscale” density
(22) 
modeling the density of interferers around the reference transmitter . The term can be seen as the probability that an interferer at distance to does not sense the ongoing transmission of .The behavior of can be described by the spatial shape function of scenario d), cf. Fig. 2. 3) The density in (22) is then used to describe the interference around the reference receiver , although it reflects the interference experienced by the reference transmitter at . This simplification step increases analytic tractability at the cost of losing accuracy. The level of accuracy loss resulting from step 3) is next studied for the Rayleigh fading model with and . We use the same notion as in (18) to measure the relative accuracy loss as a function of , i.e.,
(23) 
where the second argument in now highlights the dependence on . Fig. 5 shows the accuracy loss vs. for different and dB. The density of potential transmitters is . It can be seen that depending on the value of , significant errors can occur. Interestingly, the relative accuracy loss decreases with which is due to an increased outage correlation at the locations and . The error remains small for small , i.e., shortrange communications, for the same reason as above. For very large , tends to zero which results from both and tending to zero.
Vi Conclusion
We extended prior work on interference modeling for wireless networks with isotropic but not necessarily stationary spatial distribution of nodes. The interference statistics and the outage probability were analyzed as a function of (i) an arbitrary receiver location inside the network and (ii) an arbitrary but isotropic node distribution. For the path loss exponents and closedform expressions were obtained while bounds were derived for other cases. The developed model led to some interesting insights that could not have been revealed previously due to limiting the analysis to stationary models only. The usefulness of the results was discussed and demonstrated through examples in network analysis related to outage probability, local throughput characterization and carriersensing mechanisms. It was found that the developed model increases model accuracy significantly and provides an adequate tool to describe locationdepended performance in networks with practical node distribution.
Appendix A Integral Identities
Identity 1.
If , , than
(24) 
where is the Legendre polynomial [35]. We will be using , leading to the Legendre polynomial given by .
Identity 2.
Let , where . Define , . By [35] and using the substitution , we have
(25) 
Identity 3.
Let , where . Then,
(26) 
where
(27) 
Proof:
The basic idea is to decompose the integrand into partial fractions and to apply Identity 1 and 2, yielding (26) after some algebraic manipulations. Note that according to [35], (24) and (25) hold only for realvalued parameters. However, we verified that they also hold for complexvalued parameters. ∎
Appendix B Proofs
Ba Proof of Theorem 1
We want to compute the expectation
(28) 
where the expectation is with respect to the interferer locations and the channel gains . Since the expectation operator linear, we can compute the expectation with respect to all first, i.e., . Applying the Campbell Theorem, yields
(29) 
Changing to polar coordinates and exploiting the isotropy property, we can rewrite (29) as
(30) 
We then apply Identity 1 to the inner integral of (30) to obtain
(31) 
Finally using product integration and applying Identity 2 to (31) and verifying the convergence of the upper limit using the constraint for some , yields the result.∎
BB Proof of Theorem 2
To prove that is a.s. infinite, we analyze its Laplace transform and check if it is zero for all . In the PPP case, the Laplace transform of the interference field is given by [11, 13]
(32) 
For (32) to become zero, the integral must not converge. We write
(33) 
where (a) follows from the inequality for [36]. Inserting the path loss function into the righthand side of (33) yields
(34) 
At the upper limit of the inner integral the integrand behaves as . So, the condition , where , is sufficient for the divergence of the integral. Because the Laplace transform of becomes zero in this case, this concludes the proof.∎
BC Proof of Theorem 3
BD Proof of Theorem 4
We follow the idea of dominant interferers which was introduced in [7] to bound the tail probability : let denote the set of all interferers, where each one taken by itself already results in the event . Clearly,
(35) 
where the set is denotes the elements outside on the same domain, i.e., . Note that both sums in (35) are nonnegative. Hence, we can now write
(36)  
where (a) follows from removing the “nondominant” part in (35), (b) follows from the definition of and (c) is a consequence of being Poisson distributed with mean . Using [11, Corollary 2.1.2], in (36) can be computed as
(37) 
We now translate by the vector to obtain a centric coordinate system, cf. Fig. 6. After switching to polar coordinates we obtain
(38) 
We now exploit the fact that is subharmonic around in a region . Let denote the maximal radius for which is contained in . Then, (38) can be bounded as
(39)  
where (a) follows from limiting the upper integration limit to . Inequality (b) is a consequence of subharmonicity [37, Ch.10], and (c) follows from Definition 1. Inserting (39) in (36) yields the result.∎
BE Proof of Theorem 5
We write
where (a) follows from algebraic manipulations and the i.i.d. property of the . (b) follows from the probability generating functional and the Laplace functional of a PPP [11]. After noting that for exponentially distributed , the integral is computed using Identity 1 and Identity 2 for and Identity 3 for .
Note that (a) in the proof holds for general point processes and some approximation techniques for computing the righthand side already exist [38]. The (b) part is for PPPs only.∎
BF Proof of Corollary 7
Solving (14) for and multiplying by yields the local transmission capacity
(40) 
after substituting . We are interested in the case . Using Identity 2 (case ) we find that
(41) 
Assuming that nodes employ pseudonoise sequences, FH effectively thins out the point process of interferers by (interference avoidance, ), while DS reduces the interference by a factor of (interference averaging, ), cf. [7]. Hence, using (40) and (41) the ratio can be written as
(42)  
where (a) follows from (14) and Remark 1 and (b) follows from algebraic manipulations. Now we show that the integral in (42) is finite. Assuming