Modeling and Analysis of Cellular Networks using Stochastic Geometry: A Tutorial

• SG models spatial randomness in wireless networks and can incorporate fading, shadowing, and power control.Source: §I · p. 1
• Aggregate interference is generally non-Gaussian because distance-dependent path loss makes nearby interferers dominant, producing an impulsive, heavy-tailed α-stable distribution.Source: §I · p. 1
• Cellular-network SG modeling addresses spatial reuse and interference exclusion regions around receivers rather than idealized hexagonal deployments.Source: §I · p. 1, A. Using SG for Cellular Networks · p. 2
• The tutorial distinguishes itself through cellular-network focus, wireless-communication detail, and analytical treatment of Gaussian signaling approximation accuracy.Source: B. Motivation & Contribution · p. 2
• SG analysis relates cellular-network parameters and design variables to spatially averaged outage probability, error probability, and rate metrics.Source: §II · p. 3

Network model

• Base-station locations are modeled by a general two-dimensional point process, with repulsion capturing minimum-distance correlations from network planning.The process may be infinite for analytical simplicity, while finite-network modifications remain possible.Source: B. Network Abstraction · p. 4
• RSS-based association assigns the test user to the nearest base station x0, and the remaining ordered distances form the set ˜Ψ.The serving distance r0 is excluded from the interference set, and analysis initially conditions on a given r0.Source: B. Network Abstraction · p. 4

Signal model

• The received signal includes the intended transmission, aggregate interference from other base stations, and noise, with symbols and Rayleigh fading modeled as independent random variables.The test user is placed at an arbitrary origin for stationary point processes; otherwise, the analysis depends on location.Source: B. Network Abstraction · p. 4

Interference characterization

• Because interfering base-station locations and their number are uncertain, the analysis characterizes aggregate interference statistically through its PDF, characteristic function, or moments rather than instantaneously.Nearby interferers dominate because distance-dependent path loss prevents the central limit theorem from producing a Gaussian interference distribution.Source: B. Network Abstraction · pp. 4–5

SG tools

• Campbell’s theorem converts expectations of point-process sums into intensity integrals, while the PGFL converts random products into integrals useful for the interference characteristic function.Campbell’s theorem can obtain first moments and may extend to second moments, but higher-order moments and PGFL expressions are not generally tractable for arbitrary processes.Source: C. SG Analysis for iagg for General… · pp. 5–6

Tractability

• General point-process interference analysis is often analytically intractable, so the paper commonly approximates the interfering process by a PPP to retain tractability.The PPP approximation is motivated by the limited availability of closed-form higher-order moment and PGFL expressions for general point processes.Source: C. SG Analysis for iagg for General… · pp. 5–6

PPP approximation

• The PPP approximation is tractable because its PGFL, stationarity, and spatial ergodicity simplify interference analysis.Studies cited in the paper report close SINR agreement between PPP analysis and simulations using actual cellular topologies.Source: §IV · p. 6
• PPP models for repulsive cellular deployments should be parameterized by intensity λ and interference boundaries, especially the inner boundary.The outer boundary is usually taken as infinite because faraway base stations contribute negligibly.Source: §IV · p. 7

Interference characterization

• Using the PPP PGFL, the analysis derives the characteristic function of aggregate complex interference with an exclusion region around the receiver.The derivation uses the exclusion boundary implied by received-signal-strength association and is valid for η > 2.Source: A. Interference Characterization · p. 7
• The aggregate interference is circularly symmetric but neither Gaussian nor α-stable distributed.Without an exclusion region, it has an α-stable distribution with infinite moments; receiver protection produces finite moments under suitable conditions.Source: A. Interference Characterization · pp. 7–8
• Interference power is infinite at η = 2 or r0 = 0, whereas all cumulants and moments are finite for η > 2 and r0 > 0.For η > 2, interference power decays with exclusion radius as r0^(2−η) and increases linearly with λ and P.Source: A. Interference Characterization · p. 8
• For η > 2 and r0 > 0, positive finite kurtosis indicates an interference distribution with a heavier tail than the Gaussian distribution.The characteristic function, cumulants, expected power, and kurtosis expose these properties.Source: A. Interference Characterization · p. 8

PDF comparison

• The aggregate-interference PDF has a heavier tail than the Gaussian PDF, while smaller exclusion distance makes it approach the α-stable distribution.The comparison uses Gaussian and α-stable PDFs with matching parameters.Source: B. Numerical Results for iagg · p. 8

KS comparison

• The KS statistic compares the entire CDFs, so it does not capture deviations in tail probabilities.Fig. 6 uses relative KS distance to compare iagg with Gaussian and α-stable distributions.Source: B. Numerical Results for iagg · p. 8, C. Section Summary · p. 8
• The aggregate interference is neither Gaussian nor α-stable; as r0 increases, it deviates from α-stable behavior and approaches Gaussian behavior.The reported distributional relationship is summarized through the KS comparisons in Fig. 6.Source: B. Numerical Results for iagg · p. 8

Section summary

• The section concludes that PPP networks with receiver exclusion require interference models beyond Gaussian and α-stable assumptions.The analysis derives the aggregate-interference characteristic function and moments before turning to error probability.Source: C. Section Summary · p. 8

PDF comparison

• Fig. 5 plots the PDF of Re{iagg} for λ = 1 BS/km2, P = 10 W, η = 4, and r0 values of 250 m and 500 m.The PDFs are obtained by numerically inverting the characteristic function in (10).Source: B. Numerical Results for iagg · p. 8, Figure 5 · p. 9

Motivation

• Standard SEP expressions derived for Gaussian noise are not legitimate for cellular interference because aggregate interference is non-Gaussian.The issue arises when interference is treated under the same assumptions as AWGN.Source: §V · p. 9

Conditional Gaussian representation

• The EiD approach represents aggregate interference as a conditional Gaussian random variable and averages the conditional error probability to obtain the unconditional result.The representation is established by matching characteristic functions and therefore distributions.Source: §V · p. 9, A. Conditional Gaussian Representation… · pp. 9–10
• Conditioning on the random Gaussian-component variances makes interference plus noise circularly symmetric complex Gaussian with total variance equal to summed component variances plus N0.This conditional representation permits AWGN-based SNR formulas to be extended to cellular SINR analysis.Source: A. Conditional Gaussian Representation… · p. 10

ASEP analysis

• Exact ASEP analysis conditions on the Gaussian representation, applies the AWGN expression, and then deconditions over the non-Gaussian variables.The final expression involves an additional averaging step and remains conditioned on r0.Source: B. ASEP with Non-Gaussian Cellular… · p. 10
• Theorem 2 gives downlink ASEP for M-QAM in a PPP cellular network with Rayleigh fading, universal frequency reuse, and no intra-cell interference.The result is stated for a user at distance r0 from its serving base station.Source: B. ASEP with Non-Gaussian Cellular… · p. 11

Section summary

• The EiD pipeline characterizes interference with SG, constructs an equivalent Gaussian representation, performs conditional AWGN analysis, and deconditions to obtain ASEP.The exact expression is computationally intensive because it includes an integral involving a sum of hypergeometric functions.Source: C. Section Summary · p. 11

Gaussian signaling approximation

• Gaussian signaling abstracts interfering symbols by drawing them from a unit-variance complex Gaussian constellation, simplifying error-rate analysis.It is an approximation for interfering symbols from distinct constellations, not an approximation that makes aggregate interference Gaussian.Source: §VI · p. 11, A. Validation · p. 11
• Gaussian signaling directly yields the conditional Gaussian representation needed to use AWGN ASEP expressions.The paper states that this approximation avoids EiD complexity without compromising modeling accuracy.Source: §VI · p. 11

Validation

• The Gaussian signaling approximation preserves the aggregate-interference distribution and power of the original model.The characteristic functions and power expressions have equivalent forms, although higher-order even cumulants differ.Source: A. Validation · p. 12

Gaussian signaling approximation

• Conditioned on network geometry and channel gains, Gaussian interfering symbols make the aggregate interference conditionally complex Gaussian.This representation permits AWGN-based ASEP expressions for the received signal.Source: B. Approximate Error Probability… · p. 12

ASEP analysis

• The exact analysis, Gaussian approximation, and simulations are compared for ASEP across BS intensities, while Gaussian aggregate-interference modeling gives a loose estimate.The comparison confirms that the central limit theorem does not apply to aggregate interference in this setting.Source: B. Approximate Error Probability… · p. 13

Gaussian signaling approximation

• Gaussian signaling preserves the aggregate interference distribution, odd and second cumulants, and interference power for any constellation size.Differences from exact interference occur only in even cumulants above order two.Source: C. Section Summary · p. 13

ASEP analysis

• The approximation reduces the ASEP exponent from a sum of M hypergeometric functions to one hypergeometric function without sacrificing ASEP accuracy.For η = 4, the ASEP further reduces to a computationally simple inverse tangent function.Source: C. Section Summary · p. 13

Unified performance analysis

• The Gaussian approximation replaces the exact characteristic-function route with the Laplace transform of aggregate interference power, which is easier to derive and evaluate.The same transform supports ASEP, outage probability, and ergodic-rate calculations.Source: C. Section Summary · p. 13
• The SINR CDF is sufficient to characterize both outage probability and ergodic rate, providing a simpler alternative to the more involved ASEP expressions.Rate outage depends on interference and/or fading, and BER-based outage is also defined under the Gaussian signaling approximation.Source: A. Definition of Outage Probability… · p. 14

SINR CDF derivation

• The SINR CDF can be derived from the Laplace transform of aggregate interference because Rayleigh fading makes the useful-channel power exponentially distributed.This avoids directly computing the generally unavailable interference-power PDF.Source: B. SINR Distribution · p. 14

Scope and limitation

• The interference-power PDF is unavailable in closed form for practical cellular-network cases with an inner interference boundary at r0.Closed-form distributions occur only in special PPP cases whose interference boundaries extend from 0 to infinity.Source: B. SINR Distribution · p. 14

Performance metrics

• Theorem 4 gives downlink ergodic rate and outage probability for a PPP cellular network with Rayleigh fading, universal reuse, and no intra-cell interference.Both metrics are stated for a user at service distance r0.Source: B. SINR Distribution · p. 15
• The derivation obtains the SINR CDF from the interference Laplace transform and then uses it to compute ergodic rate and outage probability.This follows by substituting the transform expressions into the SINR result and applying the metric formulas.Source: B. SINR Distribution · p. 15
• Monte Carlo validation shows that the outage and ergodic-rate expressions capture effects of interferer intensity and interference boundaries.These metrics are presented as simpler alternatives for characterizing network behavior, with abstraction potentially hiding some underlying behavior.Source: B. SINR Distribution · p. 15

Section Summary

• ASEP, outage probability, and ergodic rate under Gaussian signaling all require the Laplace transform of aggregate interference power.For η = 4, the presentation focuses on simple transform evaluations, using a modulation-dependent argument for ASEP and threshold T for outage and rate.Source: C. Section Summary · p. 15, §VIII · p. 15

Baseline model

• The baseline model equates the interference exclusion distance with the service distance, while interference intensity and exclusion region distinguish network models.These parameters determine how the aggregate-interference transform changes across models.Source: §VIII · p. 15

Random Link Distance r0

• A random service distance requires an additional averaging step, producing a spatially averaged transform used for ASEP and the SINR CDF.The conditional transform is first obtained for r0, then averaged over its distribution.Source: A. Random Link Distance r0 · p. 16

Section Summary

• Outage probability and ergodic rate offer simpler network-behavior characterizations than ASEP, but their abstractions may hide detailed network behavior.The comparison is made between the simpler outage/rate expressions and the more involved ASEP formulas.Source: C. Section Summary · p. 15, A. Random Link Distance r0 · p. 16

Random Link Distance r0

• Noise changes the spatial averaging step because the noise term must be averaged together with the interference transform.The simplified averaging statement therefore does not apply directly in the prominent-noise case.Source: A. Random Link Distance r0 · p. 16

Load-aware networks

• Load-awareness is incorporated through an activity factor p that thins interfering BS intensity per channel while leaving the service-distance distribution based on λ.This reflects partially loaded networks in which some channels are idle.Source: B. Load-Aware Modeling · pp. 16–17

Multi-tier model

• Multi-tier cellular networks are modeled as independent PPP tiers, each with its own transmit power Pk, intensity λk, and path-loss exponent ηk.Users associate according to biased RSS controlled by tier-specific bias factors.Source: C. Multi-tier Cellular Networks · p. 17

Interference analysis

• Per-tier performance uses the aggregate interference from all tiers, with each tier contributing through its own intensity and interference boundary.Independence between tiers permits combining the per-tier interference transforms.Source: C. Multi-tier Cellular Networks · pp. 17–18

Association and biasing

• Bias factors control tier association and can increase small-cell coverage to offload users from macro BSs to small BSs.The figure illustrates this mechanism using different bias factors with common BS locations.Source: C. Multi-tier Cellular Networks · p. 17, Figure 11 · p. 18
• Tier-specific service-distance distributions depend on relative tier powers, bias factors, and path-loss exponents.For the common case ηk = 4, the distribution reduces to a simpler form.Source: C. Multi-tier Cellular Networks · p. 18
• Unbiased RSS association reduces SINR-dependent performance metrics to the single-tier case, independent of tier count, transmit powers, and BS intensities.This simplification follows from the expression for η = 4.Source: C. Multi-tier Cellular Networks · p. 19

Interference Coordination and Frequency Reuse

• The available spectrum is divided into Δ sub-bands, with reuse coordinated among base stations because PPP deployments do not support traditional hexagonal-grid reuse schemes.Each base station selects a sub-band through coordination with neighboring base stations.Source: D. Interference Coordination and… · p. 19
• User-centric coordination enlarges the interference-protection region to r_I = r_{Δ−1}, while association uses r_0 and their joint distribution.Averaging over the joint PDF of r_{Δ−1} and r_0 yields the spatially averaged interference Laplace transform.Source: D. Interference Coordination and… · p. 19
• The reuse analysis becomes a double-integral expression, but it remains efficiently evaluable in time and complexity relative to Monte Carlo simulations.The additional integration is the main analytical cost identified for coordination and frequency reuse.Source: D. Interference Coordination and… · p. 19
• Coordinated frequency reuse improves outage probability as the reuse factor Δ increases by changing both the interference boundary and interferer intensity.Figure 13 validates the corresponding interference-transform expression and shows this performance effect.Source: D. Interference Coordination and… · p. 19, Figure 13 · p. 19

Uplink Case Study

• In the uplink, per-UE power control maintains received power ρ at the serving base station, while each base station assigns a unique channel per user.The model assumes universal frequency reuse, no intra-cell interference, and a dense independent UE PPP.Source: E. Uplink Transmission · p. 20
• Uplink outage is higher than downlink outage because uplink transmission power is limited and association does not create geographical interference protection for uplink transmissions.Figure 15 verifies the uplink expression and the PPP approximation for interfering UEs.Source: E. Uplink Transmission · p. 20, Figure 15 · p. 21

General Fading

• Changing useful-link fading from exponential can prevent outage probability and ASEP from being expressed through the aggregate-interference Laplace transform, although ergodic rate remains expressible that way.General interfering-link fading preserves the framework but may make the interference transform more involved.Source: F. General Fading · pp. 20–21

Nakagami-m

• For integer Nakagami-m fading, the SINR CDF can be expressed using the Laplace transform of aggregate interference.The derivation uses the gamma CDF with integer shape parameter and the Laplace-transform identity.Source: 1) Nakagami-m · p. 21

Slow Fading

• With slow fading, tractability remains when association selects the base station offering the highest received signal strength.The displacement theorem captures shadowing by scaling the PPP intensity with a shadowing fractional moment.Source: 2) Additional Slow Fading · p. 21

MIMO Receive Diversity

• The tutorial presents receive diversity as a simple MIMO case study because a unified analytical framework for all MIMO configurations is difficult to provide.The treatment intentionally avoids the broader MIMO-system details available elsewhere.Source: G. Multiple Input Multiple Output… · p. 22
• With Nr receive antennas, the intended effective channel gain is gamma distributed with shape parameter Nr, while each interfering-link gain remains unit-mean exponential.This follows from summing Nr independent unit-mean exponential desired-link gains and independence between desired and interfering channels.Source: G. Multiple Input Multiple Output… · p. 22
• The receive-diversity SINR CDF uses the aggregate-interference Laplace transform with r_I = r_0, and Figure 17 validates the resulting expression.The example connects receive diversity to the earlier Nakagami-m analysis because the intended gain is gamma distributed.Source: G. Multiple Input Multiple Output… · p. 22, Figure 16 · p. 22
• MIMO analysis is more involved because precoding alters interfering signals and correlations may arise across antenna-branch interference.Even under Rayleigh fading, MIMO fading is no longer generally exponential.Source: G. Multiple Input Multiple Output… · p. 22

Network MIMO

• CoMP transmission serves each user from its nearest n base stations, whose signals therefore do not contribute to interference.The model considers a downlink single-tier cellular network with single-antenna base stations and user-centric CSI-agnostic cooperation.Source: H. Network MIMO · pp. 22–23
• Cooperation increases the geographical interference-protection region to rI = rₙ₋₁.The nearest n cooperating base stations are excluded from the interference field.Source: H. Network MIMO · p. 23

Discussion

• With two receive antennas and reuse factor 3, outage probability at T = 0 dB drops from almost 50% to below 5%.The combined use of receive diversity and frequency reuse produces the reported reduction.Source: §I · p. 23
• With an appropriate system model, PPP-based SG analysis can capture realistic network performance and provide acceptable performance characterization.The discussion attributes pessimistic results to assumptions such as universal frequency reuse, saturation, and peak base-station transmit power rather than to PPP itself.Source: §I · p. 23
• Simple system models can reveal parameter trends, but their results are illustrative and do not provide true numerical performance values.This scope boundary applies when trends are prioritized over absolute values.Source: §I · p. 23

Future directions

• PPP has limited ability to represent complex wireless topologies, motivating new tractable point processes for networks with D2D, V2V, and M2M communication.PPP is characterized by intensity and an interference boundary, which provide limited degrees of freedom for more complex correlations.Source: A. New Point Processes · p. 24

PPP distributions

• The appendix gives the distance distribution from a generic location in R² to the nearest point of a PPP with intensity λ.Source: Appendix I · p. 25
• The joint distance distribution to the nearest and nth points of a PPP is provided.Source: Appendix I · p. 25

PPP identities

• The probability generating functional is stated for a measurable function f and a PPP Φ.Source: Appendix I · p. 25

Voronoi geometry

• The appendix gives the area distribution of a generic PPP-Voronoi cell, using c = 3.575 for the R² Voronoi tessellation.Source: Appendix I · p. 25
• For independent PPPs of base stations and users, the appendix provides the probability mass function for the number of users in a generic base-station Voronoi cell.The point-process intensities are λb and λu.Source: Appendix I · p. 25

Auxiliary expressions

• The appendix includes the probability that a user associates with tier k and an averaging technique for a unit-mean Gamma variable, a random variable specified by its Laplace transform, and a constant.It also lists a PPP expression involving intensity, distance, path-loss exponent, and moments of the interfering symbol.Source: Appendix I · p. 25, APPENDIX II LEMMA 1 IN [165] · p. 26, APPENDIX III LT OF ζ · p. 26