#### Abstract

Stochastic orders on point processes are partial orders which capture notions like being larger or more variable. Laplace functional ordering of point processes is a useful stochastic order for comparing spatial deployments of wireless networks. It is shown that the ordering of point processes is preserved under independent operations such as marking, thinning, clustering, superposition, and random translation. Laplace functional ordering can be used to establish comparisons of several performance metrics such as coverage probability, achievable rate, and resource allocation even when closed form expressions of such metrics are unavailable. Applications in several network scenarios are also provided where tradeoffs between coverage and interference as well as fairness and peakyness are studied. Monte-Carlo simulations are used to supplement our analytical results.

#### 1. Introduction

Point processes have been used to describe spatial distribution of nodes in wireless networks. Examples include randomly distributed nodes in wireless sensor networks or ad hoc networks [1–4] and the spatial distributions for base stations and mobile users in cellular networks [5–10]. In the case of cognitive radio networks, locations of primary and secondary users have been modeled as point processes [11–14]. Random translations of point processes have been used for modeling of mobility of networks in [15]. Stationary Poisson processes provide a tractable framework, but suffer from notorious modeling issues in matching real network distributions. Stochastic ordering of point processes provides an ideal framework for comparing two deployment/usage scenarios even in cases where the performance metrics cannot be computed in closed form. These partial orders capture intuitive notions like one point process being more dense or more variable. Existing works on point process modeling for wireless networks have paid little attention to how two intractable scenarios can be nevertheless compared to aid in system optimization.

Recently stochastic ordering theory has been used for performance comparison in wireless networks which are modeled as point processes [16–20]. Directionally convex (DCX) ordering of point processes and its integral shot noise fields have been studied in [16]. The work has been extended to the clustering comparison of point processes with various weaker tools, including void probabilities and moment measures, than DCX ordering in [17]. In [18], usual stochastic ordering of random variables capturing carrier-to-interference ratio has been established in cellular systems. Ordering results for coverage probability and per user rate have been shown in multiantenna heterogeneous cellular networks [19]. In [20], Laplace functional (LF) ordering of point processes has been introduced and used to study interference distributions in wireless networks. Several examples of the LF ordering of specific point processes have been also introduced in [20], including stationary Poisson, mixed Poisson, Poisson cluster, and binomial point processes.

While limited examples of the LF ordering have been provided in [20], in this paper we apply the LF ordering concept to several general classes of point processes such as Cox, homogeneous independent cluster, perturbed lattice, and mixed binomial point processes which have been used to describe distributed nodes of wireless systems in the literature. We also investigate the preservation properties of the LF ordering of point processes with respect to independent operations such as marking, thinning, random translation, and superposition. We prove that the LF ordering of original point processes still holds after applying these operations on the point processes. To the best of our knowledge, there is no study on LF ordering of general classes of point processes and their preservation properties in the literature. Through the preservation properties of LF ordering with respect to independent operations on a point process such as marking, thinning, random translation, and superposition which is done in our manuscript, we can consider several effects of real systems such as propagation effects over wireless channels, multiple access schemes, heterogeneous network scenarios, and mobile networks and compare performances without having to obtain closed form results for a wide range of performance metrics such as coverage probability, achievable rate, and resource allocation under these real system effects. In addition to the performance comparison, the stochastic ordering of point processes provides insights into the system design and evaluation such as network deployment and user selection schemes as we summarize in this paper.

The paper is organized as follows: In Section 2, we introduce mathematical preliminaries. Section 3 introduces ordering of point processes. In Section 4, we show the preservation properties of LF ordering. Sections 5.1 and 5.2 introduce applications of stochastic ordering of point processes in wireless networks. We also elaborate the qualitative insights of the stochastic ordering of point processes on the system design and evaluation in Section 5.3. Section 6 presents simulations to corroborate our claims. Finally, the paper is summarized in Section 7.

#### 2. Mathematical Preliminaries

##### 2.1. Stochastic Ordering of Random Variables

Before introducing ordering of point processes, we briefly review some common stochastic orders between random variables, which can be found in [21, 22].

###### 2.1.1. Usual Stochastic Ordering

Let and be two random variables (RVs) such thatThen is said to be smaller than in the* usual stochastic order* (denoted by ). Roughly speaking, (1) says that is less likely than to take on large values. To see the interpretation of this in the context of wireless communications, when and are distributions of instantaneous SNRs due to fading, (1) is a comparison of outage probabilities. Since are positive in this case, is sufficient in (1).

###### 2.1.2. Laplace Transform Ordering

Let and be two nonnegative RVs such thatThen is said to be smaller than in the* Laplace transform* (LT)* order* (denoted by ). For example, when and are the instantaneous SNR distributions of a fading channel, (2) can be interpreted as a comparison of average bit error rates for exponentially decaying instantaneous error rates (as in the case for differential-PSK (DPSK) modulation and Chernoff bounds for other modulations) [23]. The LT order is equivalent tofor all* completely monotonic* (c.m.) functions [22, pp. 96]. By definition, the derivatives of a c.m. function alternate in sign: , for , and . An equivalent definition is that c.m. functions are positive mixtures of decaying exponentials [22]. A similar result to (3) with a reversal in the inequality states thatfor all that have a completely monotonic derivative (c.m.d.). Finally, note that . This can be shown by invoking the fact that is equivalent to whenever is an increasing function [22] and that c.m.d. functions in (4) are increasing.

##### 2.2. Point Processes and Random Measures

Point processes have been used to model large-scale networks [1, 2, 13, 24–29]. Since wireless nodes are usually not colocated, our focus is on* simple* point processes, where only one point can exist at a given location. In addition, we assume the point processes are locally finite; i.e., there are finitely many points in any bounded set. Unlike [20], stationary and isotropic properties are not necessary in this paper. In what follows, we introduce some fundamental notions that will be useful.

###### 2.2.1. Campbell’s Theorem

It is often necessary to evaluate the expected sum of a function evaluated at the point process . Campbell’s theorem helps in evaluating such expectations. For any nonnegative measurable function which runs over the set of all nonnegative functions on ,The intensity measure of in (5) is a characteristic analogous to the mean of a real-valued random variable and defined as for bounded subsets . Thus, is the mean number of points in . If is stationary, then the intensity measure simplifies as for some nonnegative constant , which is called the intensity of , where denotes the dimensional volume of . For stationary point processes, the right side in (5) is equal to . Therefore, any two stationary point processes with the same intensity lead to equal average sum of a function (when the mean value exists).

A random measure is a function from Borel sets in to random variables in . The* Laplace functional * of random measure is defined by the following formulaThe Laplace functional completely characterizes the distribution of the random measure [25]. A point process is a special case of a random measure where the measure takes on values in the nonnegative integer random variables. In the case of the Laplace functional of a point process, can be written as in (6). As an important example, the Laplace functional of Poisson point process of intensity measure isIf the Poisson point process is stationary, the Laplace functional simplifies with .

###### 2.2.2. Laplace Functional Ordering

In this section, we introduce the Laplace functional stochastic order between random measures which can also be used to order point processes.

*Definition 1. *Let and be two random measures such thatwhere runs over the set of all nonnegative functions on . Then is said to be smaller than in the Laplace functional (LF) order (denoted by ).

In this paper, we focus on the LF order of point processes unless otherwise specified. Note that the LT ordering in (2) is for RVs, whereas the LF ordering in (8) is for point processes or random measures. They can be connected in the following way:Hence, it is possible to think of LF ordering of point processes as the LT ordering of their aggregate processes. Intuitively, the LF ordering of point processes can be interpreted as the LT ordering of their aggregate interferences. The Laplace transform of interference is particularly convenient for the determination of outage probabilities under certain fading distributions (e.g., Rayleigh) [30] and is also useful in performance comparison. The LF ordering of point processes also can be translated into the ordering of coverage probabilities and spatial coverage which will be discussed in detail in Section 5. The probability generating functional is another mathematical tool for studying point processes [1]. The probability generating functional will provide the same results in our scenarios. However, we will focus on the Laplace functional of point process in our study since it is easy to relate the Laplace functional to interference metrics as mentioned above.

###### 2.2.3. Voronoi Cell and Tessellation

The Voronoi cell of a point of a general point process consists of those locations of whose distance to is not greater than their distance to any other point in ; i.e.,The Voronoi tessellation (or Voronoi diagram) is decomposition of the space into the Voronoi cells of a general point process.

#### 3. Ordering of General Classes of Point Processes

The examples for LF orderings of some specific point processes have been provided in [20]. In this section, we introduce the LF ordering of general classes of point processes.

##### 3.1. Cox Processes

A generalization of the Poisson process is to allow for the intensity measure itself being random. The resulting process is then Poisson conditional on the intensity measure. Such processes are called* doubly stochastic Poisson processes* or* Cox processes*. Consider a random measure on . Assume that, for each realization , an independent Poisson point process of intensity measure is given. The random measure is called the driving measure for a Cox process. The LF ordering of Cox processes depends on their driving random measures.

Theorem 2. *Let and be two Cox processes with driving random measures and , respectively. If , then .*

*Proof. *The proof is given in Appendix A.

The mixed Poisson process is a simple instance of a Cox process, where the random measure is described by a positive random constant so that . Since the Laplace functional of the mixed Poisson process can be expressed as , using (7), and because and the c.m. property of , the LF ordering of mixed Poisson processes has the following relationship: if , then .

##### 3.2. Homogeneous Independent Cluster Processes

A general cluster process is generated by taking a parent point process and daughter point processes, one per parent, and translating the daughter processes to the position of their parent. The cluster process is then the union of all the daughter points. Denote the parent point process by , and let be the number of parent points. Further let , be a family of finite points sets, the untranslated clusters or daughter processes. The cluster process is then the union of the translated clusters:If the parent process is a lattice, the process is called a* lattice cluster process*. Analogously, if the parent process is a Poisson point process, the resulting process is a* Poisson cluster process*.

If the parent process is stationary and the daughter processes are finite point sets which are independent of each other, are independent of , and have the same distribution, the procedure is called* homogeneous independent clustering*. In this case, only the statistics of one cluster need to be specified, which is usually done by referring to the* representative cluster*, denoted by which is distributed the same as any . In this class of point processes, the LF ordering depends on the parent process and the representative process as follows.

Theorem 3. *Let and be two homogeneous independent cluster processes having representative clusters and , respectively. Also, let and be the parent point processes of two homogeneous independent cluster processes and , respectively. If and , then .*

*Proof. *The proof is given in Appendix B.

##### 3.3. Perturbed Lattice Processes with Replicating Points

Lattices are deterministic point processes defined aswhere is a matrix with , the so-called generator matrix. The volume of each Voronoi cell is and the intensity of the lattice is [31]. The perturbed lattice process is a lattice cluster process. Denote the lattice point process by , and let be the number of lattice points. Further let , be untranslated clusters. In each cluster, the number of daughter points are a random variable , independent of each other and identically distributed. Moreover, these points are distributed by some given spatial distribution. The entire process is then the union of the translated clusters as in (11). If the replicating points are uniformly distributed in the Voronoi cell of the original lattice, the resulting point process is a stationary point process and called a* uniformly perturbed lattice process*. If, moreover, the numbers of replicas are Poisson random variables, the resulting process is a stationary Poisson point process [32]. Now, we can define the following LF ordering of such point processes.

Theorem 4. *Let and be two uniformly perturbed lattice processes with numbers of replicas being nonnegative integer valued random variables and , respectively, and with the same mean . If , then .*

*Proof. *The proof is given in Appendix C.

Based on Theorem 5.A.21 in [21], the* smallest* and* biggest* LT ordered random variables can be defined as follows: Let be a random variable such that and let be a random variable degenerate at . Let be a nonnegative random variable with mean . ThenFrom Theorem 4 and (13), the uniformly perturbed lattice processes with replicating points with nonnegative integer valued distribution and will be the* smallest* and* biggest* LF ordered point processes, respectively, among uniformly perturbed lattice processes with the same average number of points. The smallest LF ordered uniformly perturbed lattice process exhibits clustering since some Voronoi cells have 2 points but other cells do not have any point. This observation is in line with the intuition that clustering diminishes point processes in the LF order.

##### 3.4. Mixed Binomial Point Processes

In binomial point processes, there are a total of fixed points uniformly distributed in a bounded set . The density of the process is given by where is the volume of . If the number of points is random, the point process is called a mixed binomial point process. As an example, with Poisson distributed , the point process is called a finite Poisson point process. The intensity measure of mixed binomial point processes is . In these point processes, one can show the following.

Theorem 5. *Let and be two mixed binomial point process with nonnegative integer valued random distribution and , respectively. If , then .*

*Proof. *The proof is given in Appendix D.

Similar to Theorem 4, Theorem 5 enables LF ordering of two point processes whenever an associated discrete random variable is LT ordered.

#### 4. Preservation of Stochastic Ordering of Point Processes

In what follows, we will show that the LF ordering between two point processes is preserved after applying independent operations on point processes such as marking, thinning, random translation, and superposition of point processes.

##### 4.1. Marking

Consider the dimensional Euclidean space , , as the state space of the point process. Consider a second space , called the space of marks. A marked point process on (with points in and marks in ) is a locally finite, random set of points on , with some random vector in attached to each point. A marked point process is said to be independently marked if, given the locations of the points in , the marks are mutually independent random vectors in , and if the conditional distribution of the mark of a point depends only on the location of the point it is attached to.

Lemma 6. *Let and be two point processes in . Also let and be independently marked point processes with marks with identical distribution in . If , then .*

*Proof. *The proof is given in Appendix E.

##### 4.2. Thinning

A thinning operation uses a rule to delete points of a basic process , thus yielding the* thinned point process *, which can be considered as a subset of . The simplest thinning is *-thinning*: each point of has probability of suffering deletion, and its deletion is independent of locations and possible deletions of any other points of . A natural generalization allows the retention probability to depend on the location of the point. A deterministic function is given on , with . A point in is deleted with probability and again its deletion is independent of locations and possible deletions of any other points. The generalized operation is called *-thinning*. In a further generalization the function is itself random. Formally, a random field is given which is independent of . A realization of the thinned process is constructed by taking a realization of and applying -thinning to , using for a sample of the random field . Given and given , the probability of in also belonging to is . As long as an independent thinning operation regardless of , and is applied on point processes, the LF ordering of the original pair of point processes is retained.

Lemma 7. *Let and be two point processes in , and let and be independently thinned point processes both with any identical independent thinning operation which could be either , or -thinning on both and . If , then .*

*Proof. *The proof is given in Appendix F.

Since the thinned point process is a locally finite random set of points on , with a binary random variable in attached to each point, independent thinning can be considered as the independent marking operation on a point process as discussed in the previous section.

##### 4.3. Random Translation

In this section, the stochastic operation that we consider is random translation. Each point in the realization of some initial point process is shifted independently of its neighbors through a random vector in where are independent of each other and the conditional distribution of a random vector of a point depends only on the location of the point . The resulting process is . The random translation preserves the LF ordering of point process as follows.

Lemma 8. *Let and be two point processes in , and let and be the translated point processes with common distribution for the translation . If , then .*

*Proof. *The proof is given in Appendix G.

Similar to the independent thinning operation, since the random translated point process is a locally finite random set of points on , with some random vector in attached to each point, the random translation can be considered as the independent marking operation on a point process.

##### 4.4. Superposition

Let and be two point processes. Consider the unionSuppose that with probability one the point sets and do not overlap. The set-theoretic union then coincides with the superposition operation of general point process theory. The superposition preserves the LF ordering of point processes as follows.

Lemma 9. *Let and be mutually independent point processes and and be the superposition of point processes. If for , then .*

*Proof. *The proof is given in Appendix H.

#### 5. Applications to Wireless Networks

In the following discussion, we will consider the applications of stochastic orders to wireless network systems.

##### 5.1. Cellular Networks

In this section, the comparisons of performance metrics will be derived based on the LF ordering of point processes for spatial deployments of base stations (BSs) and mobile stations (MSs).

###### 5.1.1. System Model

We consider the downlink cellular network model consisting of BSs arranged according to some point process in the Euclidean plane. For the deployment of BSs, a deterministic network such as lattice points or stochastic network such as a Poisson point process may be considered. Consider an independent collection of MSs, located according to some point process which is independent of . Figure 1 shows an example of cellular network consisting of stationary Poisson point processes with different intensities for BSs and MSs, respectively. For a traditional cellular network, assume that each user associates with the closest BS, which would suffer the least path-loss during wireless transmission. It is also assumed that the association between a BS and a MS is carried out in a large time scale compared to the coherence time of the channel. The cell boundaries are defined through the Voronoi tessellation of the BS process. Our goal is to compare performance metrics such as total cell coverage probability through stochastic ordering tools. The spatial coverage of cellular networks is also compared based on the LF order of the BS point processes.

In order for the total cell coverage probability to be compared, the signal to interference plus noise ratio (SINR) of a user at should be quantified. The effective channel power between a user and its associated typical BS is , which is a nonnegative RV. The SINR with additive noise power is given bywhere is the path-loss function which is a continuous, positive, nonincreasing function of the Euclidean distance from the user located at to the typical BS . The following is an example of a path-loss model [1, 26, 33, 34]:for some , and , where is called the path-loss exponent, determines whether the path-loss model belongs to a singular path-loss model () or a nonsingular path-loss model (), and is a compensation parameter to keep the total received power normalized regardless of the values of path-loss exponent. In (15), is the accumulated interference power at a user located at given bywhere denotes the set of all BSs which is modeled as a point process and is a positive random variable capturing the (power) fading coefficient between a user and the interfering BS. Moreover, are i.i.d. random variables and independent of and .

###### 5.1.2. Ordering of Performance Metrics in a Cellular Network

In the following discussion, we will introduce performance metrics involving the stochastic ordering of aggregate process in the cell which is associated with the BS . By studying spatial character of networks and investigating the spatial distributions of mobile users, we can compare system performances using stochastic ordering approach without actual system performance evaluation. In addition, the preservation properties of the LF order in Lemmas 6–9 guarantee that the performance comparison results based on the LF ordering of point processes are not changed with respect to any identical independent random operation on or such as marking, thinning, translation, and superposition.

*(a) Total Cell Coverage Probability*. In multicast/broadcast scenarios, multiple users receive a common signal from their associated BS. Therefore, the probability that SINRs of all served users are greater than a minimum threshold is an important measure to ensure the signal reception quality of every user and it is called* total cell coverage probability*. This metric can be ordered, if the underlying point processes are LF ordered.

Theorem 10. *Let be an arbitrary fixed BS deployment. Also, let and be two point processes for MS deployments, and let and be two user location sets belonging to the typical Voronoi cell . and are random variables that are functions of and , respectively. Let in (15) be exponentially distributed independent RVs, and let in (17) be independent RVs that are also independent of , , and . If , then for where is over all RVs, , , , and .*

*Proof. *The proof is given in Appendix I.

Note that Theorem 10 holds for any arbitrary as long as it is identical under both scenarios. Equation (18) shows that LF ordering of user point processes implies ordering of the multivariate complementary cumulative distribution functions (CCDFs) of the SINRs.

*(b) Network Spatial Coverage*. The network spatial coverage is an important performance metric to design BS deployment in cellular networks. We assume that BSs are distributed by a point process and each BS has a fixed radius of coverage . Denote the random number of BSs covering a fixed location bywhere is a -dimensional ball of radius centered at the point . Denote the probability generating function of the random number of BSs covering location by . Since , is a c.m. function with . Note that represents the probability whether the location is covered by at least one BS from the definition of the probability generation function. Thus, if , then from the property of LT ordering in (3) and consequently . This means the probability that any arbitrary point in is covered by at least one BS with the cell deployment by is always less than the probability with . Due to random effects of real systems such as shadowing, different transmission power per each BS, and obstacles, the range of coverage can be a nonnegative random variable. Since the random range of coverage can be considered as independent marking in Section 4.1, the ordering of spatial coverage probabilities still holds from Lemma 6 under the assumption that the random range of coverage is independent of the BS deployment . The ordering of spatial coverage probabilities also holds with random ellipse or square instead of a ball from Lemma 6.

With performance metrics such as network spatial coverage and network interference, our study can provide design guidelines for network deployment to increase spatial coverage of networks or provide less interference from networks. As an example, consider the effects of clustering as in Section 3.2. In the case that the daughter clusters are assumed to be Poisson, the clustering of nodes diminishes a point process in the LF order. From an interference point of view, the clustering of interfering nodes causes less interference in the LT order between interference distributions. This translates into an increased coverage probability and improved capacity for the system. However, the clustering of nodes also causes less spatial coverage. Therefore, the proper point process for network deployment should be studied for balancing between interference and spatial coverage.

##### 5.2. Cognitive Networks

In the following discussion, we will consider the applications of stochastic orders to cognitive network systems where there is an increasing interest in developing efficient methods for spectrum management and sharing.

###### 5.2.1. System Model

Let us consider an underlay cognitive radio network which contains a primary user (PU) and many secondary users (SUs) with an average interference power constraint . The PU is located at the origin. The SUs are uniformly randomly located in a certain area [35–37]. It is assumed that there is a BS to coordinate the SUs’ transmission. In order to satisfy the interference constraint, the BS selects only active number of users out of total users through user selection schemes and the selected active SUs are allowed to transmit their signals to the BS. User selection creates a thinned point process for the active SUs which has uniformly distributed points over the area . The resulting thinned point process can be considered as a mixed binomial point process as in Section 3.4. The system model of the cognitive radio network is illustrated in Figure 2. Under the system model, based on Theorem 5, one can design user selection schemes which guarantee the same average interference power from the active SUs to the PU and the same average sum rate for the active SUs. On the other hand, the* distribution* of the instantaneous interferences from the active SUs to the PU is such that they are LT ordered and cause ordered performance metrics such as a coverage probability and an achievable rate of the PU according to the user selection schemes. In other words, the smaller LT ordered interference distribution guarantees the better coverage probability and achievable rate of the PU. We now discuss these in detail.

*(a) Average Interference Power Constraint*. The instantaneous aggregate interference from the active SUs to the PU can be expressed aswhere is a point process for the active SUs and is a positive random variable capturing the (power) fading coefficient between an active SU located at and the PU. From Campbell’s theorem, the average of aggregate interference from the active SUs, , is the same as long as the average number of the active SUs is fixed to regardless of the distribution for (for a proof, please see Appendix J). Therefore, we need to select in order to satisfy the average interference power constraint .

*(b) Average Sum Rate*. On the other hand, if there is no interference between the active SUs by adopting code-division multiple access (CDMA), the instantaneous sum of achievable rates of the active SUs, , is also random and can be expressed aswhere is an effective fading channel between an active SU located at and the BS for SUs located at , is the interference power from the PU to the SU located at , and is an additive noise power. From Campbell’s theorem, the average of sum of achievable rates of active SUs, , is the same regardless of the point processes for the active SUs as long as the average number of the active SUs is equal to (proof is the same as the average interference case in Appendix J with a different function).

###### 5.2.2. Randomized User Selection Scheme for Secondary Users

We now discuss different ways to perform randomized user selection which corresponds to different ways of thinning the SU point process. One can apply the stochastic ordering tool to design user selection schemes based on Theorem 5 which produce mixed binomial point processes for active SUs. Given which satisfies the interference constraint , the smaller LT ordered random number of active SUs provides the smaller LF ordered point process ; that is, from Theorem 5. Consequently, the resulting mixed binomial point processes cause LT ordered aggregate interferences to the PU, . The LT ordered aggregate interferences from the active SUs to the PU yield ordered performance metrics for the PU such as a coverage probability and an achievable rate due to the LT ordered interferences and the coverage probability and the achievable rate having c.m. property with respect to the interferences [20]. On the other hand, the average sum of achievable rates of the active SUs and the average of interference power remain the same, and as discussed in Sections 5.2.1(a) and 5.2.1(b).

We now give examples of a user selection scheme using the stochastic ordering approach. If the active SUs are chosen among total SUs with a probability independently, the number of active SUs is a binomial random variable, . The number of active SUs can follow discrete distributions other than binomial if different modes of operation are adopted. For another example, if the SUs are selected before a predetermined failure numbers with a probability occuring, then the number of selected active SUs follows a negative binomial distribution with parameter and denoted as . Since [38], from Theorem 5. Therefore, the aggregate interferences from the active SUs are LT ordered , while and . For these operations, the BS for SUs only needs to know the number of total SUs and the average number of active SUs .

From the discussion in Sections 5.2.1(a) and 5.2.1(b), it is noted that the smaller LT ordered random number of active SUs causes less interference to the PU, while the average interference and average sum rate are the same. From Theorem 5.A.21 in [21] and the discussion in Section 3.3, the smallest LT ordered random number is and the biggest is the fixed . Even though the causes the smallest LT ordered interference to the PU, when , it causes large instantaneous interference to the PU at this moment. Otherwise, the provides balanced traffic from the active SUs since the fixed number of random set of users is always selected. However, it causes the bigger LT ordered interference to the PU than any other distributions for active random SUs . Therefore, the proper distribution for random number of active SUs should be studied for balancing peakyness of instantaneous interference power and fairness of active SUs’ traffic.

##### 5.3. Discussion and Tradeoffs

###### 5.3.1. Coverage versus Interference

As shown in the previous sections, the stochastic ordering of point processes can be applied to performance comparison and system design. Since the performance metrics depend on the stochastic ordering of point processes for spatial distribution of nodes, we can study spatial character of networks and investigate the spatial distributions of transmitting nodes in order to evaluate system performances using stochastic ordering approach. Then our investigations can be applied to the comparison of performance metrics such as number of served users and total cell coverage probability without having closed form expressions for the metrics. With performance metrics such as network spatial coverage and network interference [20], our study can provide design guidelines for network deployment to increase spatial coverage of networks or provide less interference from networks. As an example, consider the effects of clustering. The clustering of nodes reduces a point process in the LF order. From an interference point of view, the clustering of interfering nodes causes less interference in the LT order between interference distributions. This translates into an increased coverage probability and improved capacity for the system [20]. However, the clustering of nodes also causes less spatial coverage. Therefore, the proper point process for network deployment should be studied for balancing between interference and spatial coverage.

###### 5.3.2. Burstiness versus Fairness

The stochastic ordering approach can also be applied to design user selection schemes for SUs in cognitive networks which satisfy the constraints of the system. Under the system setup described in Section 5.2, the smaller LT ordered active SUs causes less interference to the PU, while the average interference and average sum rate are the same. From Theorem 5.A.21 in [21] and the discussion in Section 3.3, the smallest and biggest LT ordered can be obtained as follows:where is any nonnegative integer random variable, is , and is the fixed . Even though the smallest LT ordered random number causes the smallest LT ordered interference to the PU, it causes more bursty traffic from the active SUs and more instantaneous interference when . Otherwise, the biggest LT ordered provides balanced traffic from the active SUs since the fixed number of random set of users is always selected. However, it causes the bigger LT ordered interference to the PU than any other distributions for active random SUs . Therefore, the proper distribution for active random SUs should be studied for balancing burstiness of interference and fairness of active SUs’ traffic.

###### 5.3.3. Random Effects

The system performances do depend on not only spatial deployment of a network, but also random effects of real systems such as random fading effects, random moving of mobile stations, coexisting of heterogeneous networks, and shadowing. From the preservation properties of LF ordering of point processes, it is noted that the random effects do not change the LF ordering of original point process when independent identical random effects are considered in two different spatial network deployment scenarios. As an example, consider two mobile networks where every node is moving. If the initial point processes of two mobile networks are LF ordered, then interferences raised from the mobile networks are LT ordered. Consequently the ordered interferences imply ordered performance metric results such as coverage probability and achievable rate. The ordered performance metric results are not changed regardless of the different mobility models such as constrained i.i.d. mobility, random walk, and Brownian motion. Since these random effects do not affect the performance comparison results and design guidelines using stochastic ordering approach, we need to focus the stochastic ordering of spatial deployments of nodes.

#### 6. Numerical Results

In this section, we now corroborate our theoretical results through numerical examples. Since the number of BSs is smaller than the number of MSs in general, we choose the smaller value of for BSs than for MSs.

##### 6.1. Cellular Networks

We show in Figure 3 the total cell coverage probabilities. It is assumed that the BS distribution is a stationary Poisson point process , and the compared user distributions follow also a Poisson point process and mixed Poisson process with the same intensity . Since the user distributions are LF ordered as , from Theorem 10, the total cell coverage probability of the users distributed by is greater than that of in the entire range of the SINR threshold as shown in Figure 3. Using the stochastic ordering approach, one can compare the coverage probabilities by investigating the spatial distributions of mobile users without actual system evaluation.

##### 6.2. Cognitive Networks

In Figure 4, the CDFs of the SINR of the primary user are shown where the spatial distributions of the active SUs are different mixed binomial processes. The LT ordered number of the active SUs of the negative binomial, Poisson, and binomial RVs ensures the ordering of the corresponding mixed binomial processes, from Theorem 5. Since the aggregate interferences from LF ordered point processes for spatial distributions of the active SUs are LT ordered, we observe which are the ratio between the PU’s effective fading channel and the interference from active SUs. This is because when the PU’s effective fading channel is exponentially distributed and interferences are LT ordered, the SIRs are reverse ordered in* usual stochastic order* [20]. The aforementioned relationship between SIR and interference distributions holds as long as the CCDF of the PU’s effective fading channel is a c.m. function by the property of LT ordering in (3). An example is the exponential distribution. Despite the ordering of the SIRs in Figure 4, the average of the sum of achievable rates of the active SUs in (21) is the same regardless of LF ordered point processes as shown in Figure 5.

In order to guarantee the best coverage probability for the PU, one can use the stochastic ordering approach to design the user selection schemes which cause less interference from active SUs to the PU. Since the negative binomial distribution for the number of active SUs causes the smallest LT ordered interference, one can design the user selection scheme with the minimum information such as the total number of SUs and the average number of active SUs which provides the negative binomial distribution as discussed in Section 5.2.

#### 7. Summary

In this paper, Laplace functional ordering of broad classes of point processes was investigated. We showed the preservation of LF ordering of point process with respect to several operations, such as independent marking, thinning, random translation, and superposition. We introduced the applications of LF ordering of point processes to wireless networks such as cellular networks and cognitive networks, which provided guidelines for design of user selection schemes and transmission strategy for wireless networks. Tradeoffs between coverage and interference as well as fairness and peakyness were also discussed. The power of this approach is that network performance comparisons can be made even in cases where a closed form expression for the performances is not analytically tractable. We verified our results through Monte-Carlo simulations.

#### Appendix

#### A. Proof of Theorem 2

A Cox process is a Poisson point process conditional on the realization of the intensity measure. Therefore, the Laplace functional of the Cox process with driving random measure can be expressed as follows [25]:where is the Laplace functional of the Poisson process of intensity measure and is the set of intensity measures in (A.1). Equation (A.2) follows from the definition of Laplace functional of Poisson point process in (6), and (A.3) follows from the definition of Laplace functional of random measure in (8). Then, the proof follows from Definition 1.

#### B. Proof of Theorem 3

Let denote the Laplace functional of , i.e., the Laplace functional of the representative cluster translated by . It can be expressed as follows:From (B.1), the Laplace functional of the homogeneous independent cluster process with parent process can be expressed as follows [25]:From Definition 1 and (B.2), if and the same in both, then . Similarly, if and the same in both, then . Therefore, if and , then .

#### C. Proof of Theorem 4

The uniformly perturbed lattice process can be considered as the superpositions of independent finite point processes corresponding to each Voronoi cell. The finite point process in one of Voronoi cells consists of uniformly distributed random number of points . The Laplace functional of the finite point process can be expressed as follows:where is a region of the Voronoi cell. By denoting , is a c.m. function of since . Therefore, if , then we have the following relation:

From (C.2) the following LF ordering is obtained:Since the finite point processes are LF ordered, their superpositions are also LF ordered from Lemma 9. The proof of Theorem 4 is completed.

#### D. Proof of Theorem 5

Similar to Appendix C, the Laplace functional of a mixed binomial point process with random number of points can be expressed as follows:where is the bounded set of the point process. By denoting , is a c.m. function of since . Therefore, if , then we have the following relation:From (D.2) the following LF ordering is obtained:The proof of Theorem 5 is completed.

#### E. Proof of Lemma 6

The Laplace functional of an independently marked point process with a nonnegative function can be expressed as follows [39]:Here . Since , is a nonnegative function of . Then, the Laplace functional of marked point process follows fromTherefore, from (E.4) and the definition of LF ordering in (8), if , then .

#### F. Proof of Lemma 7

If is the Laplace functional of , then that of iswhere . From Definition 1, if , then . The -thinning is subset of -thinning. Analogous formula for -thinning follows by averaging with respect to the distribution of the random process . Since the inequality holds under every realization , their expectations also hold the inequality.

#### G. Proof of Lemma 8

Let denote the common distribution for the translations . For , the Laplace functional after random translation takes the form where . From Definition 1, if , then .

#### H. Proof of Lemma 9

Let and be mutually independent point processes and and be the superposition of point processes. The Laplace function of superposition of mutually independent point processes can be expressed as follows: converges if and only if the infinite sum of point processes is finite on bounded area . Therefore, from (H.1) and the definition of LF ordering in (8), if for , then .

#### I. Proof of Theorem 10

To prove Theorem 10, we first express for a generic Equation (I.3) follows from the assumption that and in (15) and (17) are independent under the given realizations of . From the assumption that is exponentially distributed (Rayleigh fading), (I.5) follows. In (I.5), is a nonnegative function of . Therefore, if , then (18) follows by (3) because is a c.m. function.

#### J. Proof of

From Campbell’s theorem in (5) and the intensity measure of the mixed binomial point processes , the average of interference can be expresses as follows:Conditioned on , is always equal. Therefore, the expectation of (J.2) with respect to is still equal. The proof is completed.

#### Data Availability

The simulation data used to support the findings of this study are available from the first author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work was supported in part by the National Science Foundation under Grant CCF 1117041. This work was also supported in part by the Institute for Information and Communications Technology Promotion (IITP) grant funded by the Korea government (MSIT) (No. 2017-0-01973 (Korea-Japan) International Collaboration of 5G mmWave Based Wireless Channel Characteristic and Performance Evaluation in High Mobility Environments).