Modeling Cellular Networks in Fading Environments with Dominant Specular Components

Modeling Cellular Networks in Fading Environments with Dominant Specular Components


Stochastic geometry (SG) has been widely accepted as a fundamental tool for modeling and analyzing cellular networks. However, the fading models used with SG analysis are mainly confined to the simplistic Rayleigh fading, which is extended to the Nakagami-m fading in some special cases. However, neither the Rayleigh nor the Nakagami-m accounts for dominant specular components (DSCs) which may appear in realistic fading channels. In this paper, we present a tractable model for cellular networks with generalized two-ray (GTR) fading channel. The GTR fading explicitly accounts for two DSCs in addition to the diffuse components and offers high flexibility to capture diverse fading channels that appear in realistic outdoor/indoor wireless communication scenarios. It also encompasses the famous Rayleigh and Rician fading as special cases. To this end, the prominent effect of DSCs is highlighted in terms of average spectral efficiency.

1 Introduction

Modern cellular networks have evolved from the ubiquitous hexagonal grid to irregular multi-tier structure that randomly changes from one geographical location to another [1, 2, 3]. To cope with such evolution, efforts are spent to develop tractable stochastic geometry (SG) models that account for large-scale spatial randomness as well as different sources of uncertainties that emerge in modern cellular networks, such as multipath fading, shadowing, and power control [4]. The last decade witnessed significant progress developing all aspects of the SG models, except for the fading environments. In the context of multi-path fading, the Rayleigh fading assumption is almost a common factor in the literature. In addition to the Rayleigh fading case, there are some proposals that incorporate Nakagami-m fading into tractable SG analysis [5, 6].

Both the Rayleigh and Nakagami-m assumptions are favorable because they lead to a desirable exponential expressions for the conditional (i.e., conditioned on network geometry) signal-to-interference-plus-noise-ratio (SINR) performance metrics, which enables averaging via the moment generating function (MGF) of the interference. Although [4, Sec.III], [7, 8] show alternative techniques to circumvent the necessity of such representation, there are yet no SG models including state-of-the-art fading models that may include individually resolvable multipath components, which are denoted as multiple dominant specular components (DSCs) [9]. As argued in [9, 10, 11, 12, 13], both the Rayleigh and Nakagami-m models may fail to capture realistic fading environments. The Rayleigh fading ignores the line-of-sight (LOS) component in the received signal, which is prominent in outdoor cellular communication. It is also a single parameter fading model that is not flexible enough to model complex indoor fading environment. The Nakagami-m fading can neither aptly capture fading with LOS component nor accurately model indoor fading environment as argued in [9, 10, 12, 11]. Furthermore, the tractability of SG analysis with Nakagami-m fading necessities an integer value for , which decreases the model’s flexibility to capture complex fading environments.

In order to develop SG models in terms of multipath fading, it is required to incorporate a flexible fading model that can span different fading channels including those with multiple DSCs. The state-of-the-art studies in multi-path fading (see [9, 10, 11, 12, 13] and the references therein) show that fading channels with a finite number of DSCs with dominant powers in addition to a diffuse component appear in several practical indoor/outdoor scenarios. A simple example with single DSC appears cellular networks where one LOS path exists between the BS and a mobile user. In this case, the resulting fading environment follows the well-known Rician fading model. Furthermore, multiple DSCs may be created by reflections from metal objects (e.g., light-posts, cars) in close proximity to the BSs and/or user, in which the reflections would have comparable amplitude to the dominant LOS path due to the negligible power absorption factor of metals. A fading channel with a small number of specular components without a diffuse component may appear in millimeter wave (mmW) communication where high directional antennas are used for short-range communication [14]. Also, severe NLOS fading environments can be modeled via a small number of specular components without diffuse components [13, 11, 12]. Last but not least, field measurements for indoor multipath fading channels confirm the existence of multiple DSCs in addition to the diffuse component [15]. Therefore, there is a need to incorporate new flexible fading models that capture realistic multipath fading environments with multiple DSCs.

As a step forward to account for fading channels with multiple DSCs, this paper incorporates the generalized two-ray (GTR) fading model into tractable stochastic geometry analysis. The GTR model is adopted due to its mathematical elegance and practical significance1. The GTR fading channel is characterized by four tunable parameters that explicitly account for two DSCs plus the diffuse component, which gives high flexibility to capture a diverse spectrum of fading channels. With the proper adjustment of the GTR parameters, it captures several fading channels as special cases such as deterministic, Rician, Rayleigh, hyper-two ray, and hyper-Rayleigh [13, 12]. Hence, accounting for GTR fading in a SG model leads to unified analysis for a diverse spectrum of fading channels. It is worth noting that the two DSCs, which are explicitly accounted for by the GTR model, are not necessarily LOS components. However, they may also be considered as two NLOS paths that experienced favorable propagation when compared to the others in the diffuse part. Using this interpretation, the severe multi-path fading channels captured by the GTR model can be justified [13, 12].

To the best of the authors’ knowledge, this work is the first to incorporate such a flexible fading model into a tractable SG analysis. It is worth mentioning that we incorporate the GTR fading model into a simplistic network setup to observe its explicit effects. Analyzing more advanced network setups with GTR fading is postponed to future work.

2 System Model

We consider a downlink single-tier cellular network with single-antenna BSs that are deployed according to the PPP with intensity . Each element belongs to and denotes the location of the i BS. Single-antenna users’ equipment (UEs) are spatially distributed according to a stationary point process with intensity such that . All BSs transmit with a constant transmit power of P. Users associate to the BSs according to the average radio signal strength (RSS) rule. Universal frequency reuse is adopted with no intra-cell interference.

It is assumed that the signal power decays according to the power-law with the distance , where is the path-loss exponent. We focus on the performance of a test UE located at the origin. According to Slivnyak’s theorem, there is no loss of generality with this assumption. For the sake of simple exposition, we define the set that contains the ordered distances from the test user to the BSs in . Following the RSS rule, and are the distance between the test UE and, respectively, his serving BS and the i nearest interfering BS.

2.1 The GTR Fading Model

We consider a GTR model with two specular components plus a diffuse component. The received baseband signal can be represented as:


where and are respectively the intended and i BS interfering symbols which are assumed to be drawn from Gaussian codebooks, and are the (amplitude, phase) of the two DSCs of the intended signal, and are the zero mean Gaussian distributed intended signal diffuse components with total power of , and are the (amplitude, phase) of the two DSCs of the interfering signal from the i BS, and are the zero mean Gaussian distributed i BS interfering signal diffuse components with total power of , and is the complex Gaussian noise. Let be the phase difference between the two DSCs and denotes the probability density function (PDF) of , then according to [13], the GTR model is defined by four parameters, namely the diffuse power (DP) , the specular-component-to-diffuse-power-ratio (SDPR) , peak-to-average-specular-component-power-ratio (PASPR) , and .

Observing and , we notice that the GTR model reduces to Rayleigh fading for and reduces to Rician fading for . It also captures the two-ray, hyper-two ray, and hyper-Rayleigh by letting , , and properly choosing . The distribution function discriminates between different types of GTR fading models. Particularly, is uniformly distributed in the range of in the GTR-U model, whereas is uniformly distributed in the truncated range of in the GTR-T model. In the GTR-V model, follows the Von Mises distribution. Note that the GTR-T and GTR-V capture the cases when the angles of the two specular components are correlated, which appears in some practical scenarios [13].

2.2 Methodology of Analysis

For a given realization of network geometry and channels’ gains, the interference term in (2.1) has a Gaussian distribution. Hence, treating interference as noise, the instantaneous SINR for the test user can be expressed as:


where is excluded from in (2) because the serving BS does not contribute to the interference. In SG analysis, we are interested in spatially averaged performance metrics, which requires the conditional (i.e., conditioning on ) PDFs of and , in addition to the PDF of . In a PPP network with RSS association, the PDF of is known to be , [1]. Also, the conditional PDF of is straightforward to obtain from the PDF of . However, the aggregate interference is usually characterized via its Laplace transform (LT)2. Hence, only performance metrics that are expressed in terms of the LT of can be evaluated.


For the sake of organized presentation, we devote Section 3 to derive the conditional LT of the aggregate interference with different GTR fading environments. In Section 4, we first express the averaged spectral efficiency (i.e., ) in terms of the conditional LTs derived in Section 3. Then, the averaging step over is done to obtain the spatially average spectral efficiency.

3 Interference Characterization in GTR Fading Environment

We assume that all fading channels are independent from the BSs’ locations, independent from each other, and are identically distributed (i.i.d.) according to the GTR model with parameters: , , , and . Let represents the expected power gain of the fading channel, then from the definition of , we have . For general , the LT of the aggregate interference in GTR fading environment is given by the following lemma.

Lemma 1.

In a GTR fading environment with parameters , , , and , the LT of the conditional aggregate interference (conditional on ) in a PPP cellular network with intensity is given by equation (2.2), where , and


Refer to Appendix .1. ∎

The expression in (2.2) contains two integrals, which may increase the computational complexity. Therefore, we obtain a simpler lower-bound on (2.2) which is given by the following lemma.

Lemma 2.

The LT given in equation (2.2) is lower bounded by (6), where is the confluent hypergeometric function, is given by equation (4), and


Refer to Appendix .2. ∎

It is worth emphasizing that is independent of the LT variable and is a function of the fading parameters only. Hence, for a known fading parameters, is a constant w.r.t. the interference LT. The accuracy of the lower-bound in (6) is validated in Section 5. Based on Lemma 1 and Lemma 2, the interference for different GTR fading models can be characterized. For instance, the interference with Rician faded channels is characterized by the following corollary,

Corollary 1.

In a Rician fading environment with parameters and , the LT of the conditional aggregate interference in a PPP cellular network with intensity is given by equation (2.2), and lower-bounded by (6), where and,


The corollary is obtained from Lemma 1 and Lemma 2 by setting . ∎

Interestingly, and for Rician fading is obtained in closed forms. The interference in GTR-U fading case is characterized via the following corollary,

Corollary 2.

In a GTR-U fading environment with parameters , , and , the LT of the conditional aggregate interference in a PPP cellular network with intensity is given by equation (2.2), and lower-bounded by (6), where and,


The corollary is obtained from Lemma 1 and Lemma 2 by setting . ∎

The GTR-U assumes that the phases of the two DSCs (i.e., and ) are independent, and hence, the value of the phase difference is uniformly distributed from . In some cases, correlation between and may exist, which limits the range that spans. The scenarios where and are correlated are captured via the GTR-T and GTR-V cases. In the GTR-T, the distribution of is assumed to be uniformly distributed in the range of , where is the truncation parameter that can be manipulated to capture different correlation scenarios. In the GTR-V case, is assumed to follow the Von Mises, given by


where is the phase difference distribution parameter that can be manipulated to concentrate the PDF of around a certain value in the range. It is worth mentioning that both the GTR-T and GTR-V reduce to the GTR-U for and , respectively. The LT of the conditional aggregate interference in the cases of GTR-T and GTR-V is given by the following corollaries

Corollary 3.

In a GTR-T fading environment with parameters , , , and , the LT of the conditional aggregate interference in a PPP cellular network with intensity is given by equation (2.2), and lower-bounded by (6), where and,


The corollary is obtained from Lemma 1 and Lemma 2 by setting for , and otherwise. ∎

Corollary 4.

In a GTR-V fading environment with parameters , , , and , the LT of the conditional aggregate interference in a PPP cellular network with intensity is given by equation (2.2), and lower-bounded by (6), where and,


The corollary is obtained from Lemma 1 and Lemma 2 by using (12) for the distribution of . ∎

The cases provided in Lemma 1, Lemma 2, and Corollaries 1 to 4 represent the case in which the DSCs arrive in addition to a diffuse component. This implicitly implies that interference signals have experienced favorable fading channels. Another case of special interest is when interference links experience severe fading, which can be captured via a two ray model without a diffuse component [13]. In this case, we set and such that is kept constant. For a general the LT of the aggregated interference in severe fading environment is given in the following lemma,

Lemma 3.

The LT of the aggregate interference in a GTR fading environment with and such that is constant, is given by


such that,


where denotes the lower incomplete gamma function and is given by equation (4).


Refer to Appendix C. ∎

For each of the special cases of GTR fading, in (18) can be evaluated via the appropriate . For GTR-U, GTR-T, and GTR-V, expressions for are given, respectively, by


It is worth highlighting that the severe fading channels are captured by the constructive and destructive interference of the two rays. Consequently, the two DSCs should be explicitly consider to account for such severe fading cases. In contrast, the Rician fading accounts only for one specular component plus the diffuse component, and hence, setting and leads to a deterministic channel gain.

4 Performance Analysis

In this section, the conditional LTs of the aggregate interference obtained in Section 3 are utilized to evaluate the spatially averaged spectral efficiency.

4.1 Average Spectral Efficiency

Starting form the Lemma introduced in [16], which states that


where .

Exploiting (21) the spectral efficiency can be directly expressed in terms of the LT of interfering and useful links. Following [13], the LT of the desired signal power while conditioning on is given by equations (22), (23), and (24) for GTR-U, GTR-T, and GTR-V fading channel, respectively, where .

Theorem 1.

The average spectral efficiency for a PPP cellular network with intensity in a GTR fading environment with parameters and for the desired and interfering fading channels is given by the following equation


where is given by equation (2.2) for a general , lower bounded by equation (6) and when it reduces to equation (17).


Follows from substituting (2) in (21) and averaging over . ∎

5 Results

Figure 1: Average spectral efficiency vs. truncation parameter (p) assuming GTR-T, Rician and Rayleigh for the desired channel and GTR-U for interfering channels. Solid lines, the diamonds and the X’s represent the results obtained analytically, by simulations and by simulations assuming Rayleigh fading for all the fading channels except the nearest 5 BSs. The chosen values are the same for the two ray model where , , and such that the average received power is kept constant for all cases.

This section presents numerical results for the spectral efficiency for GTR fading environments. All results are validated by independent system level simulations. The match between the analysis and simulation of all cases confirms the validity of our expressions. Unless otherwise stated, we used the following parameters’ values: , , , , , and .

Fig. 1 shows the spectral efficiency obtained via Theorem 2 for GTR-U fading on the interfering links and different fading models on the useful link. For the sake of fair comparison, we keep a constant for both links. For GTR-T on the useful link, the performance is highly affected by the value of . Lower implies higher correlations between the DSCs phases which lead to a high probability of destructive interference, and vice versa. The figure also shows that Rician fading on the useful link has higher performance than the Rayleigh fading due to the LOS path, and higher performance than the GTR-T fading due to the absence of destructive DSCs interference.

Figure 2: Ergodic rate assuming Rayleigh fading for the desired channel and GTR-T for the interferes channel. Solid lines, dashed lines, and the diamonds represent the rate obtained by the exact expression, the lower bound, and simulations.

In Fig. 2, we show the explicit effect of fading on the interfering links by fixing GTR-U fading on the useful link. Note that we fix the diffuse power for both the useful and interfering channels. In consistence with the previous results, Fig.2 manifests the effect of the interference fading model on the spectral efficiency. The Rayleigh fading contains the diffuse power only, and hence, has the highest spectral efficiency. In contrast, the Rician fading has the strongest channel gain due to the LOS connection, and hence, it represents the worst interference and lowest spectral efficiency. The GTR fading with multiple DSCs performance lies between the Rayleigh and Rician cases due to the probability of destructive DSCs interference. By decreasing (increasing) , the GTR-T approaches the Rayleigh (Rician) fading performance due to the higher (lower) probability of destructive interference. Last but not least, the figure confirms the tightness of the lower bound obtained by Lemma 2.

5.1 Discussion

These results emphasize the prominent effect of multiple DSCs on the network performance in terms of spectral efficiency. The results also show the flexibility of GTR-T model to capture different DSCs destructive interference scenarios, in which the correlation between the phase differences is captured by manipulating to control the range that spans.

While applying the GTR model to the useful link is quite intuitive, it might not be as much obvious why do we need to apply the GTR model to the interference links. To elaborate this point, we give the following justifications:

  • The interference is dominated by nearby interferers that may have LOS channels with the receivers especially in outdoor communication with macro BSs. As shown in Fig. 1, assuming GTR-U fading for the nearest 5 BSs only and Rayleigh fading for the other interferes does not have a noticeable effect on the averaged spectral efficiency.

  • In suburban and rural areas with parks and open spaces, LOS interference is most likely to happen.

  • The developed model is not restricted to LOS communication and may also model NLOS communication in which one or two of the NLOS paths have dominant powers.

  • For indoor environments, the specular component that go through doors and windows may have dominant powers w.r.t. the components that penetrate through walls with high attenuation coefficients. Hence, interfering signal coming from adjacent rooms/corridors may have multiple DSCs.

  • Beside the solid physical foundation of the GTR model, it offers a mathematically elegant and flexible model that spans different LOS and NLOS fading environments. For instance, it can capture extreme fading environments such as Rayleigh, the , hyper-Rayleigh and other channels [13, 12].

It is worth highlighting that the developed model parameters can be tuned to be a function of the interfering BS distance. More particularly, the fading severity can be adjusted to increase with the propagation distance, which complies with practice. Such modification is postponed to future work.

6 Conclusion

This paper presents a tractable stochastic geometry (SG) model with GTR fading channels that explicitly accounts for two dominant specular components (DSCs). Exact and lower bounds for the spectral efficiency are obtained. Depending on the phase difference between the two DSCs at the receiver, constructive/destructive interference may occur leading to high variability in the fading channel gains. To this end, the explicit effect of GTR models on the interference and useful links are investigated. The results confirm the prominent effect of DSCs fading models on the network performance in terms of spectral efficiency. As special case of the GTR, we consider the GTR-T in which the phase correlation between the two DSCs is captured by controlling the range that the phase difference spans. Finally, the model flexibility to capture several fading conditions, ranging from deterministic and favorable Rician to severe hyper-Rayleigh and hyper-two ray, in large-scale cellular networks is also highlighted.

Due to space constraints, only the outlines of the proofs are highlighted.

.1 Proof of Lemma 1

Following the same methodology in [1] and using the MGF of GTR fading in [13], the LT of is obtained as


where . By change of variables, , integration by parts, and some manipulations, (2.2) in Lemma 1 is obtained.

.2 Proof of Lemma 2

The integration in (2.2) can be only evaluated when the upper limit is . By definition, is always less than one. Since the integrand in (2.2) is positive, setting the integral upper limit to and using [17, Eq.(3.383)], the lower bound in (6) is obtained.

.3 Proof of Lemma 3

Setting and the average channel power is given by Starting form (26) and substituting by , by , setting . Then, by change of variables and integration by parts Lemma 3 is obtained.


  1. Models with more than two dominant specular components are not likely to appear in practice [13, 11].
  2. With slight abuse of notation, LT is used to denote the LT of a probability density function (PDF) of a random variable, which is equivalent to the moment generating function with negative argument.


  1. J. G. Andrews, F. Baccelli, and R. K. Ganti, “A tractable approach to coverage and rate in cellular networks,” IEEE Trans. Commun., vol. 59, no. 11, pp. 3122–3134, Nov. 2011.
  2. A. Guo and M. Haenggi, “Spatial stochastic models and metrics for the structure of base stations in cellular networks,” IEEE Trans. Wireless Commun., vol. 12, no. 11, pp. 5800–5812, Nov. 2013.
  3. W. Lu and M. D. Renzo, “Stochastic geometry modeling of cellular networks: Analysis, simulation and experimental validation,” CoRR, vol. abs/1506.03857, 2015. [Online]. Available:
  4. H. ElSawy, E. Hossain, and M. Haenggi, “Stochastic geometry for modeling, analysis, and design of multi-tier and cognitive cellular wireless networks: A survey,” IEEE Commun. Surveys Tuts., vol. 15, no. 3, pp. 996–1019, 2013.
  5. M. Di Renzo and W. Lu, “Stochastic geometry modeling and performance evaluation of MIMO cellular networks using the equivalent-in-distribution (EiD)-based approach,” IEEE Trans. Commun., vol. 63, no. 3, pp. 977–996, Mar. 2015.
  6. L. H. Afify, H. ElSawy, T. Y. Al-Naffouri, and M.-S. Alouini, “The influence of Gaussian signaling approximation on error performance in cellular networks,” IEEE Commun. Lett., Accepted 2015.
  7. M. D. Renzo and P. Guan, “Stochastic geometry modeling of coverage and rate of cellular networks using the Gil-Pelaez inversion theorem.” IEEE Commun. Lett., vol. 19, no. 9, pp. 1575–1578, Sep. 2014.
  8. M. D. Renzo, A. Guidotti, and G. E. Corazza, “Average rate of downlink heterogeneous cellular networks over generalized fading channels: A stochastic geometry approach,” IEEE Trans. Commun., vol. 61, no. 7, pp. 3050–3071, Jul. 2013.
  9. N. Beaulieu and X. Jiandong, “A novel fading model for channels with multiple dominant specular components,” IEEE Wireless Commun. Lett., vol. 4, no. 1, pp. 54–57, Feb 2015.
  10. M. Yacoub, “Nakagami-m phase-envelope joint distribution: A new model,” IEEE Trans. Veh. Technol., vol. 59, no. 3, pp. 1552–1557, March 2010.
  11. G. D. Durgin, T. S. Rappaport, and D. A. De Wolf, “New analytical models and probability density functions for fading in wireless communications,” IEEE Trans. Commun., vol. 50, no. 6, pp. 1005–1015, 2002.
  12. J. Frolik, “On appropriate models for characterizing hyper-rayleigh fading,” IEEE Trans. Wireless Commun., vol. 7, no. 12, pp. 5202–5207, Dec. 2008.
  13. M. Rao, F. Lopez-Martinez, M.-S. Alouini, and A. Goldsmith, “MGF approach to the analysis of generalized two-ray fading models,” IEEE Trans. Wireless Commun., vol. 14, no. 5, pp. 2548–2561, May 2015.
  14. S. Rangan, T. Rappaport, and E. Erkip, “Millimeter-wave cellular wireless networks: Potentials and challenges,” Proceedings of the IEEE, vol. 102, no. 3, pp. 366–385, March 2014.
  15. T. Rappaport, “Characterization of UHF multipath radio channels in factory buildings,” IEEE Trans. Antennas Propag., vol. 37, no. 8, pp. 1058–1069, Aug 1989.
  16. K. A. Hamdi, “A useful lemma for capacity analysis of fading interference channels,” IEEE Trans. Commun., vol. 58, no. 2, pp. 411–416, 2010.
  17. I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Seventh Edition.   Academic Press, 2007.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description