Meta Distribution of SIR in Ultra-Dense Networks with Bipartite Euclidean Matchings

10/29/2019 ∙ by Alexander P. Kartun-Giles, et al. ∙ 0

Ultra-dense networks maximise spatial spectral efficiency through spatial reuse. The way in which their dense limit is approached is difficult to understand mathematically, because the Euclidean combinatorial theory, despite being simple to state, is often difficult to work with, and often involves a stochastic model. Small changes in the way these models are defined can make significant differences to the qualitative properties of their predictions, such as concerning the meta distribution of the SIR, the data capacity, or their efficiency. In this paper we focus on triangular correlations and variable link distances in Euclidean space, addressing critical assumptions which are introduced in more elementary models. We study a bipartite Euclidean matching in order to study how incorporating geometric frustration leads to variation in the statistics of the meta distribution of the signal-to-interference ratio compared to these simplified models, comparing with the bipolar or downlink cellular networks studied by Haenggi and others.



There are no comments yet.


page 1

page 2

page 3

page 4

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

Geometry plays a fundamental role in the data capacity of dense wireless networks, since it limits their reuse of available spectrum [1]. The notion of spatial spectral efficiency (SSE) has therefore been introduced in order to measure the bits per second per Hz per unit volume the network can move from source to destination. On the whole, effective densification leads to greater SSE, up to some fundamental limit which is apparently independent of the communication strategy chosen [2]. See the highly cited works of Gupta and Kumar, and the following groundbreaking work of Francescetti, Di Renzo, Andrews, Haenggi, Baccelli, and others [3, 4, 5, 6, 7, 8, 9]. In an abstract setting, various ideas from “Euclidean combinatorial optimisation” [10], a field where combinatorial structures in -dimensional Euclidean space are optimised somehow, are therefore useful when theorising about this limit, and in other scenarios where spatially stochastic network models are used in wireless communications [11, 12, 13, 14, 15, 16, 17, 18].

In this paper, and following on from previous work [19, 20], we consider this problem in two-dimensional (2d) Euclidean space without boundary, from the perspective of optimal transport, which is an important topic surrounding transportation theory, measure theory, analysis and operations research, providing procedures for efficiently moving a set of ‘resources’ to a set of ‘destinations’, i.e. from mines to silos. With mines and silos, one pairs mines with silos so that each mine is assigned a unique silo, and each silo a unique mine [21]. Also, Mézard and Parisi, observed early links between Euclidean matching theory, and the physics of spin glass, a magnetic alloy important in condensed matter physics and complexity [22].

Consider user equipments (UEs) and base stations (BSs). Two binomial point processes in the Euclidean plane are paired-off in an edge-independent way known as a Euclidean matching, in particular a bipartite one, shown in Fig. 1 as studied by G. Sicuro [23], and we ask what optimal transport looks like in this setting. What is the maximum data rate in bits per second which the UE can pass to the BS [19]? This scenario is also relevant for machine-to-machine communication and sensor networks, where a massive number of devices form communication pairs, changing over time, aiming at distributing relevant data and measurement info over the full network. An upper bound on the maximum data rate concerns a functional defined on the space of bipartite matchings known as “transport plans”, which is minimised given the constraints each interferer places on its neighbouring transmitters. Certain transport plans will lead to some transmitter-receiver (Tx-Rx) pairs failing to meet an Signal-to-Interference Ratio (SIR) threshold with given reliability, so the introduction of the meta-distribution of the SIR is important in this setting. This then leads to the studies of spatial capacity [24, 25, 26]. We discuss the meta-distribution of the SIR in Section II.

Mathematically, the reliability

is defined as the conditional probability

, where is the operation threshold, and

, in our case, corresponds to the point process of transmitters and receivers in a bipartite matching. Assuming unit-mean exponential random variables for the power fading and Bernoulli variates with parameter

for the activity of each interferer, the reliability can be read as [5, Appendix]


where is the distance of the transmitter-receiver link whose reliability is computed, and the product is over the locations of interferers conditioned on the realization of .

Various Euclidean models have been studied from this perspective. For example, the bipolar model, where transmitter per unit volumes are distributed randomly in the infinite plane, and then each is assigned a receiver at a fixed distance in a uniformly random direction, has a certain fraction of Tx-Rx pairs who can achieve an SIR greater than with reliability greater than , while others fall below. The bipartite matching of Euclidean points is distinct since it incorporates Euclidean correlations, and variable link distances which shorten with . This has subtle but, as we show, very important effects on the statistics of the meta-distribution. The details of this follow in Sections II and III.

In Section II, we discuss the meta-distribution of the SIR. In Section III

we introduce the bipartite Euclidean matching, which is the unique setting of our results, and we derive the moments of the meta-distribution of the

SIR for it. While doing so, we compare the performance accuracy of various link distance models, including the popular bipolar model with nearest-neighbour communication. Finally, in Section IV, we discuss these results in light of spin glass and communication theory, and conclude.

Ii Meta Distribution

The classical analysis of wireless networks using stochastic geometry assumes suitable distributions for different parts of the system and averages the performance indicator (outage probability, data rate, etc.) over the ensemble of possible network states [27]. In equation (1

), we have already averaged over the distributions of the fading channel and the activity of interferers. If we average over the point process of interferers too, the result might not represent accurately the reliability of individual links in the network, especially when the standard deviation of the

Random Variable (RV) in (1) is large. Because of that, we may calculate the distribution (or at least the moments) of the RV in (1) instead [28].

The complementary Cumulative Distribution Function (CDF) of the RV , i.e., is usually referred to as the meta-distribution of the SIR. Given the operation threshold , this metric gives the percentage of random spatial realizations where the typical receiver meets the SIR target with reliability (probability calculated over fading and activity, see equation (1)) at least . For ergodic point processes, this percentage is equal to the fraction of links achieving a reliability higher than for each spatial realization (fixed but unknown locations).

Thus far, the moments of the meta-distribution of the SIR, , have been investigated for Poisson bipolar, cellular and heterogeneous wireless networks using the Probability Generating Functional (PGFL) of Poisson Point Process (PPP[5, 6]. The meta-distribution has also allowed to go beyond classical mean-field models [29], and maximize the actual density of active links in Poisson bipolar networks meeting a reliability constraint [24].

Assuming that the interferers follow a PPP of intensity , the integer moments of the meta-distribution can be calculated by raising (1) to the -th power and expressing the spatial average using the PGFL.


Next, we concentrate on the distribution of the link distance in Euclidean bipartite matchings.


Fig. 1: An example realisation of a bipartite Euclidean matching on the unit square with periodic boundary conditions, with pairs of points distributed uniformly at random [23]. One set of nodes transmit (red), while the other recieve (blue), though both sets of nodes may be transceivers. Isotropic radiation from the transmitter set may limit the capability of all links to achieve a signal-to-interference ratio threshold .

Iii System Model and Meta Distribtution

With and , consider the Poisson point process (PPP) of points, and also the PPP , also of points, drawn from the hypercube with periodic boundary conditions, which is equivalent to e.g. a torus when , with points distributed uniformly at random over its (flat) surface. Form a perfect bipartite matching of these points by assigning of the pairs with one end in each of and ‘active’ in such a way that every point is incident to exactly one active pair.

Call the Euclidean lengths of these edges . We reserve for the dimension of the hypercube. For , in [28], one assigns each edge its own data capacity based on the Shannon-Hartley theorem, since if the received signal power is well modelled by taking for the path loss exponent, then such a function of the edge length corresponds to the theoretical upper bound on the data capacity of a link over distance , given the bandwidth and noise are set to unity. For each matching, we therefore have a length and a capacity . The perfect matching which minimises is the shortest perfect matching, or just shortest matching. In this paper, we consider a perfect bipartite matching between the two sets minimizing the sum of Euclidean distances . We seek the proportion of links in the matching, which are able to achieve an SIR greater than with probability at least . Instead of analyzing the spatial performance of wireless ad hoc networks using a simple bipolar model as in [28, 5], we use the theory of Euclidean matchings for the geometry of the links in the network with the meta-distribution of SIR for its performance.

We briefly note that the bipartite matchings produced here are via the algorithm of Duff and Koster used in the MatLab (R2019b) ‘matchpairs’ toolbox [30]. They are often only close to optimal, due to the tractability of obtaining the unique shortest assignment.

Let us assume a realization of both and , of intensity , modeling the set of transmitters in a wireless ad hoc network deployed in the unit square . Given the number of transmitters, we deploy uniformly at random receivers, and in addition to that we place a receiver at the center of the square. Under periodic boundary conditions, all receivers in the network experience the same interference field. In addition, we assume that the point process of transmitters and receivers is ergodic. This is not as straightforward to prove as in Poisson bipolar and cellular networks in [28, 5], because nearby links experience correlated distances in a Euclidean matching. Under the ergodicity assumption (whose proof is deferred to a future study), the receiver at the origin is hereafter referred to as the typical receiver, and the meta-distribution of the SIR at the typical receiver will correspond to the fraction of links in the matching experiencing certain reliability.

In the performance evaluation of wireless ad hoc networks, the bipolar model assumes a fixed and known link distance  [28, 5]. Under the assumption that the nearest interferer can come arbitrarily close to the receiver, the moments of the meta-distribution of the SIR follow from equation (2).


where . Note that the above equation is equivalent to [5, Eq. (5)] using the diversity polynomials.

In the literature, the link distance is assigned a value either independently of the intensity of interferers or proportional to the nearest-neighbor distance  [31, Eq. (3.29)]. The latter agrees with the intuition, i.e., it takes into account the fact that for increasing density of users, a higher intensity of interferers should come along with shorter distances for the transmitter-receiver link. It is easy to see from (3) that the average success probability for nearest-neighbor communication becomes independent of the intensity . However, the scaling ignores the correlations of link distances inherent in a Euclidean matching.

We note that for , which describes the case of ultra-dense networks, the sum of link distances in bipartite Euclidean matchings in 2d space scales as  [32]. Recently, the pre-constant of the leading-order term has been also derived, showing that the scaling is proportional to  [33, 34, Eq. (6)]. Therefore we may assume that the mean link distance in a bipartite Euclidean matching scales as . One may now substitute this value instead of in (3) and proceed with the calculation of the moments. For instance, the first moment takes the following simple form:


Equation (4) can be seen as a complementary result to [31, Eq. (3.29)]. It indicates that the average success probability decreases for denser bipolar networks, while all other parameters remain fixed. Note that the new bipolar model assumes a fixed and known link distance, which is equal to the leading-order term of the mean link distance in a bipartite Euclidean matchings. As such it takes into account some of the correlation of link distances across the network through their mean. On the other hand, the classical bipolar model in [28, 5] specifies the link distance assuming that the link distances of Tx-Rx pairs are independent.

The new bipolar model in (4

) may offer an improved estimate for the link distance in a bipartite network, but it still neglects the fact that the link distances in the matching would naturally follow a distribution. Another popular distance model in wireless communications research is the Rayleigh distribution, inspired by the void probability of

PPP. In our system setup, we will consider the distribution , which has a mean equal to . The moments of the meta-distribution are calculated after substituting instead of in (3) and averaging over the above yielding


where .

Note that the Rayleigh distribution for the link distances is commonly associated with the performance evaluation of downlink cellular systems with nearest base station association. For a density of base stations and the typical receiver at the origin, the link distance follows the Rayleigh distribution  [27]. In addition, the nearest interfering base station must be located further than the serving base station. This is taken into account by setting the lower integration limit in the first line of (3) equal to instead of zero. The expression in (5) corresponds to a different system, i.e., the Rayleigh distribution has a different mean to reflect bipartite matchings in ultra-dense ad hoc networks instead of cellular systems, and the nearest interferer can also come closer to the receiver than its associated transmitter.

In bipartite Euclidean matchings the distance distribution for short links behaves linearly [22, Eq. (7)]

. As a result, a good candidate model might also be the Gamma distribution

with scale . The shape can be set to match the mean of the link distances, i.e., . In Fig. 2

, we see that the Gamma distribution provides better fit to the simulation than the Rayleigh. For comparison, we also include the exponential distribution, and the unit-step function pertinent to the bipolar models. We include the traditional bipolar model with the link distance equal to

, and the new bipolar model with link distance considering the Euclidean correlations in bipartite matchings. Clearly, both bipolar models are unable to reflect accurately the distribution of link gains in the network.


Fig. 2: Link distance distribution in bipartite Euclidean matchings with in unit square. We use the toolbox ’matchpairs’ in MatLab (R2019b) which implements an algorithm of Duff and Koster [30] to generate the transmitter-receiver pairs conditioned on the realization of the point processes . The link distances over spatial configurations are calculated and the CDF of their histogram is plotted. . The size of the bin is taken equal to . Hereafter, by ’Bipolar - literature’ we refer to the classical bipolar model with Poisson interferers of intensity , see equation (3), and transmitter-receiver link distance equal to .

For the exponential distance model, the moments of the meta-distribution are calculated by substituting instead of in (3) and averaging over an exponential distribution with rate .


where for brevity and is the complementary error function.

For the Gamma model, we average (3) over the distribution , where , yielding


In Fig. 3 we approximate the simulated meta-distribution using the various link distance models. In order to do that we calculate the first two moments for each model, see (3)(7

), and we map them to Beta distributions. We have checked that the Beta approximations provide a very good fit to the simulated meta-distributions for all the models. Note that due to the simplicity of (


) higher moments can be easily computed and used to numerically invert the characteristic function of the meta-distributions 

[5, Eq: (12)]. We see in Fig. 3 that the performance of the bipolar model is poor over the full range of the distribution and is clearly unsuitable to model the SIR in bipartite Euclidean matchings. The other models can capture the trend of the meta-distribution, with the Gamma model giving almost a perfect fit.

Let us assume that the network density increases while the product is kept fixed. For a bipolar model with link distance independent of the intensity as in [5], the average success probability should not change. By visual inspection of Fig. 3 and Fig. 4, we see that for bipartite Euclidean matchings the average success probability increases for denser networks with fixed . This is because the useful link distance (on average) decreases, while the mean interference level remains the same. Another remark is that denser networks augment the difference between the traditional and the new bipolar model. Finally, in Fig. 5 we plot the outage probability for various thresholds

. The Gamma model still outperforms all others. The bipolar model might be accurate at low thresholds, however, its prediction for the variance of the meta-distribution remains poor, see for instance Fig. 


Fig. 3: The meta-distribution of the SIR in bipartite Euclidean matchings with in a unit square. Activity , pathloss exponent , operation threshold . random spatial configurations are generated. For each configuration, we simulate the probability of successful reception (or reliability) at the typical receiver over random activity and fading realizations. The transmitter associated to the typical receiver is always active. The complementary CDF of the histogram of the reliability is finally plotted with bin size equal to .


Fig. 4: Intensity , activity . See the caption of Fig. 3 for the rest of the parameter settings.


Fig. 5: Outage probability with respect to the operation threshold . See the caption of Fig. 4 for the parameter settings.


Fig. 6: Meta-distribution at operation threshold . See the caption of Fig. 4 for the rest of the parameter settings.

Iv Discussion and Conclusions

The incorporation of a more sophisticated Euclidean model with variable link distance has been studied, and show to be significantly different from a simple model effectively trying to capture the same Tx-Rx pairing scenario. The careful introduction of Euclidean combinatorial theory into ultra-dense networks is crucial, as we see with such a large deviation in the qualitative properties of the meta-distribution statistics between the bipolar model and the bipartite Euclidean matchings.

A central development is to move beyond a random link assumption, incorporating the recent developments of Sicuro et al. in the theory of Euclidean matchings. In this paper, we have used just the leading-order term for the mean link distance in bipartite Euclidean matchings. Incorporating higher-order moments in our analysis would allow us to devise more accurate models for the performance of the network. Also, working on more general curved spaces may help understand the effect of geometry better. Nevertheless, the meeting of spin glass physics, and dense communication network theory, has significant potential for both fields.

V Acknowledgements

We thank G. Sicuro and Springer International Publishing AG for allowing reproduction of Fig. 1 from his PhD thesis [23]. We also thank Ginestra Bianconi, Justin Coon, Carl Dettmann, Harpreet Dhillon, Orestis Georgiou, Suhanya Jayaprakasam, Hyowon Kim, Eleonora Kreačić, Leo Laughlin and David Simmons for many helpful discussions.


  • [1] X. Ge, S. Tu, G. Mao, C. X. Wang, and T. Han, “5G Ultra-Dense Cellular Networks,” IEEE Wireless Communications, vol. 23, no. 1, pp. 72–79, March 2016.
  • [2] M. Franceschetti, M. D. Migliore, and P. Minero, “The capacity of wireless networks: Information-theoretic and physical limits,” IEEE Transaction on Information Theory, vol. 55, no. 8, August 2009.
  • [3] P. Gupta and P. R. Kumar, “The capacity of wireless networks,” IEEE Transactions on Infomation Theory, vol. 42, no. 6, pp. 388–404, March 2000.
  • [4] W. Lu and M. Di Renzo, “Stochastic geometry modeling of cellular networks: Analysis, simulation and experimental validation,” in Proceedings of the 18th ACM International Conference on Modeling, Analysis and Simulation of Wireless and Mobile Systems, ser. MSWiM ’15.   New York, NY, USA: ACM, 2015, pp. 179–188.
  • [5] M. Haenggi, “The meta distribution of the sir in poisson bipolar and cellular networks,” IEEE Transactions on Wireless Communications, vol. 15, no. 4, pp. 2577–2589, April 2016.
  • [6] ——, Stochastic Geometry for Wireless Networks, 1st ed.   New York, NY, USA: Cambridge University Press, 2012.
  • [7] A. K. Gupta, X. Zhang, and J. G. Andrews, “SINR and Throughput Scaling in Ultradense Urban Cellular Networks,” IEEE Wireless Communications Letters, vol. 4, no. 6, pp. 605–608, 2015.
  • [8] A. AlAmmouri, J. G. Andrews, and F. Baccelli, “A Unified Asymptotic Analysis of Area Spectral Efficiency in Ultradense Cellular Networks,” preprint available at arXiv:1712.07733, December 2017.
  • [9] Y. Wang, M. Haenggi, and Z. Tan, “The meta distribution of the sir for cellular networks with power control,” IEEE Transactions on Communications, vol. 66, no. 4, pp. 1745–1757, April 2018.
  • [10] J. E. Yukich, Probability Theory of Classical Euclidean Optimization Problems.   Springer, 1998.
  • [11] G. Knight, A. P. Kartun-Giles, O. Georgiou, and C. P. Dettmann, “Counting geodesic paths in 1-d vanets,” IEEE Wireless Communications Letters, vol. 6, no. 1, pp. 110–113, Feb 2017.
  • [12] A. P. Giles, O. Georgiou, and C. P. Dettmann, “Connectivity of soft random geometric graphs over annuli,” J. Stat. Phys., vol. 162, no. 4, p. 1068–1083, January 2016.
  • [13] A. P. Giles, O. Georgiou, and C. P. Dettmann, “Betweenness centrality in dense random geometric networks,” in 2015 IEEE International Conference on Communications (ICC), June 2015, pp. 6450–6455.
  • [14] A. P. Kartun-Giles, M. Barthelemy, and C. P. Dettmann, “Shape of shortest paths in random spatial networks,” Phys. Rev. E, vol. 100, p. 032315, Sep 2019.
  • [15] K. Koufos and C. P. Dettmann, “Temporal correlation of interference in bounded mobile ad hoc networks with blockage,” IEEE Communications Letters, vol. 20, no. 12, pp. 2494–2497, Dec 2016.
  • [16] ——, “Temporal correlation of interference and outage in mobile networks over one-dimensional finite regions,” IEEE Transactions on Mobile Computing, vol. 17, no. 2, pp. 475–487, Feb 2018.
  • [17] K. Koufos and C. P. Dettmann, “The meta distribution of the sir in linear motorway vanets,” IEEE Trans. Comms. (to appear), 2019.
  • [18] ——, “Distribution of cell area in bounded poisson voronoi tessellations with application to secure local connectivity,” Journal of Statistical Physics, vol. 176, no. 5, pp. 1296–1315, Sep 2019.
  • [19] A. Kartun-Giles, S. Jayaprakasam, and S. Kim, “Euclidean matchings in ultra-dense networks,” IEEE Communications Letters, vol. 22, no. 6, pp. 1216–1219, June 2018.
  • [20] E. Kreačić and G. Bianconi, “Statistical mechanics of bipartite z-matchings,” Europhysics Letters, vol. 126, no. 2, p. 28001, May 2019.
  • [21] C. Villani, Topics in Optimal Transport.   AMS, 2003.
  • [22] M. Mézard and G. Parisi, “The Euclidean Matching Problem,” J. Phys. France, vol. 49, no. 12, p. 2019–2025, December 1988.
  • [23] G. Sicuro, The Euclidean Matching Problem.   Springer Theses, Springer International Publishing AG, 2017.
  • [24] S. S. Kalamkar and M. Haenggi, “The spatial outage capacity of wireless networks,” IEEE Transactions on Wireless Communications, vol. 17, no. 6, pp. 3709–3722, June 2018.
  • [25] ——, “Per-link reliability and rate control: Two facets of the sir meta distribution,” IEEE Wireless Communications Letters, vol. 8, no. 4, pp. 1244–1247, Aug 2019.
  • [26] R. K. Ganti, J. G. Andrews, and M. Haenggi, “High-SIR Transmission Capacity of Wireless Networks with General Fading and Node Distribution,” IEEE Transactions on Information Theory, vol. 57, no. 5, pp. 3100–3116, 2011.
  • [27] J. G. Andrews, F. Baccelli, and R. K. Ganti, “A tractable approach to coverage and rate in cellular networks,” IEEE Transactions on Communications, vol. 59, no. 11, pp. 3122–3134, November 2011.
  • [28] R. K. Ganti and J. G. Andrews, “Correlation of link outages in low-mobility spatial wireless networks,” in 2010 Conference Record of the Forty Fourth Asilomar Conference on Signals, Systems and Computers, Nov 2010, pp. 312–316.
  • [29] S. P. Weber, X. Yang, J. G. Andrews, and G. de Veciana, “Transmission capacity of wireless ad hoc networks with outage constraints,” IEEE Transactions Information Theory, vol. 51, no. 12, pp. 4091–4102, Dec. 2005.
  • [30] I. S. Duff and J. Koster, “On algorithms for permuting large entries to the diagonal of a sparse matrix,” SIAM Journal on Matrix Analysis and Applications, vol. 22, no. 4, pp. 973–996, 2001.
  • [31] M. Haenggi and R. K. Ganti, “Interference in large wireless networks,” Foundations and Trends® in Networking, vol. 3, no. 2, pp. 127–248, 2009.
  • [32] M. Ajtai, J. Komlós, and G. Tusnády, “On optimal matchings,” Combinatorica, vol. 4, no. 4, pp. 259–264, Dec 1984.
  • [33] S. Caracciolo, C. Lucibello, G. Parisi, and G. Sicuro, “Scaling hypothesis for the euclidean bipartite matching problem,” Phys. Rev. E, vol. 90, p. 012118, Jul 2014.
  • [34] S. Caracciolo and G. Sicuro, “Scaling hypothesis for the euclidean bipartite matching problem. ii. correlation functions,” Phys. Rev. E, vol. 91, p. 062125, Jun 2015.