# Multichannel analysis of soft frequency reuse and user association in two-tier heterogeneous cellular networks

• 915 Accesses

• 3 Citations

## Abstract

The cell range expansion (CRE) is encouraged to be applied in the heterogeneous cellular networks (HCNs), to enhance the capacity by offloading macro users to small cells. However, the enhanced inter-cell interference coordination (eICIC) techniques are supposed to be used for mitigating the strong cross-tier interference suffered by the offloaded users and small cell edge users. To address this, a novel soft frequency reuse (SFR) scheme is adopted in this paper. We analyze multichannel downlink scenarios for the SFR scheme using the tools of stochastic geometry. In consideration of the random resource allocation and practical cell load model, the analytical results of coverage probability and average user rate are derived and validated through Monte Carlo methods. Furthermore, our results can reduce to simple closed-form under reasonable special case for modern urban cellular networks. The main evaluation of the performance in terms of average user rate is presented, and the optimal combination of association bias and parameters of the SFR scheme is also investigated. Numerical results show that the SFR scheme outperforms the frequency resource partitioning (FRP) scheme in any load condition. Moreover, the CRE with SFR scheme can improve the average user rate significantly.

## Introduction

Heterogeneous cellular networks (HCNs) are expected to be one of primary technologies for the emerging fifth generation (5G) mobile networks [1, 2], which can be considered as an efficient solution for dealing with the exploding data traffic demands of users. In a HCN, it is proved that the co-channel deployment of low-power base stations (BSs) (also known as small cells such as micro, pico, and femto BSs) with conventional macro BSs can yield the largest sum rate .

For a typical network, its user association rule that connects a user to a specific serving BS could substantially affect the network performance. In the existing conventional homogeneous networks, the maximum received signal strength (max-RSS) is widely adopted, which can achieve the anticipated performance. However, the max-RSS user association rule is not suitable for HCNs, since the large difference in transmit powers between small cells (e.g., pico BSs ≈30 dBm) and macro BSs(≈46 dBm) , and thereby most of users will be associated with the macro BSs. The unbalanced user load leads to overloaded macro BSs and inefficient resource utilization in small cells. To cope with this problem, the biased user association also known as cell range expansion (CRE) has been proposed , wherein the macro users are proactively offloaded to small cells. Nevertheless, the drawback of biased user association is that the offloaded macro users referred to as range-expanded small cell users are liable to experience severe interference from the nearby macro BSs. Naturally, the small cell edge users also are vulnerable users. In this context, the enhanced inter-cell interference coordination (eICIC) techniques  are expected to mitigate such strong interference.

### Motivation and related work

Previous works on the CRE in conjunction with suitable eICIC techniques can be divided into two general groups, namely, time-domain strategies  and frequency-domain strategies . One of time-domain strategies is the so-called almost blanking subframes (ABS) introduced in 3GPP Release 10, where macro BSs are periodically muted in order to mitigate interference to offloaded users . The ABS can be considered as a resource partitioning in the time domain. It can be seen in  that the user throughput can be enhanced through simulations. Using the tools of stochastic geometry, tractable expressions which can give clear insight into the performance gain are derived in [7, 9]. Furthermore, the work  analyzes the required number of ABS based on user throughput requirement, but the CRE is not captured. However, a larger association bias is required in a heavily loaded scenario, which causes macro BSs should be shut off about half the time. This is counter-intuitive and unrealistic. Hence, the reduced power subframes (RPS) are encouraged to be applied to address this issue . Instead of ABS, the RPS are allocated to macro interior users with lower transmit power. The capacity loss caused by ABS is thereby reduced because the sever interference to the small cell edge users is mitigated. The results in  show that the RPS provide better rate coverage than the ABS, with a given association bias.

We focus on the frequency-domain strategies in this paper. In [14, 15], the authors propose a frequency resource partitioning (FRP) scheme, wherein a certain fraction of frequency resources with full transmit power is preserved for offloaded users. Analytical expressions for coverage probability and rate distribution are presented in  using the tools of stochastic geometry. The optimal system performance can be achieved through jointly tune association bias and resource partitioning fraction. However, the full-load model that all the BSs are always active is not reasonable, especially for the network with CRE. This is because the assumption of full-load model cannot adequately reflect the enhancement in performance through load balancing. In detail, the interference from a BS is mainly dependent on its load. That is, the probability of a BS becomes the source of interference is directly proportional to its users. Taking such effects into consideration, the authors of [15, 19] propose a more practical cell load model that a BS is active on a given sub-channel only if the sub-channel is allocated for at least one user. In , both the coverage probability and average user rate of FRP scheme are dependent on user density and the resource partitioning fraction, apart from association bias. The results in  are more practical for design guidelines.

Nevertheless, the drawback of FRP scheme is the sacrifice of the spectral efficiency of macro users. Hence, the work  presents an evolved FRP scheme, wherein a certain fraction of frequency resources is allocated to not only offloaded macro users with full transmit power but also macro users with lower transmit power. By means of the transmit power reduction in macro tier, the spectral efficiency of macro users is guaranteed. However, the transmit power reduction is randomly applied to macro users, which probably leads to that the performance of macro users in poor situation gets worse. In addition, the small cell edge users also are vulnerable to interfered by the nearby macro BSs. Similar to the RPS, a novel soft frequency reuse (SFR) scheme proposed by [17, 18, 20] can mitigate the strong cross-tier interference by means of joint resource partitioning and transmit power reduction for macro interior users. The strong cross-tier interference suffered by offloaded macro users and pico edge users is commonly mitigated. More recently, the work  presents the spectral efficiency analysis for the proposed SFR scheme based on stochastic geometry. As a result, the system spectral efficiency can be significantly improved. Although the proportional fair resource allocation is adopted in , the analytical results cannot adequately capture the variation in performance because of the assumption of full-load model. To the best of our knowledge, none of the earlier works considered the impact of resource allocation and appropriate load model on the SFR based on stochastic geometry, which is fulfilled in this paper.

### Approach and contributions

Motivated by the works in [15, 19], we propose a general and tractable framework to analyze joint SFR, appropriate load model and CRE in a two-tier HCNs. Based on stochastic geometry, the locations of the BSs in each tier are modeled as a two-dimensional homogeneous Poisson Point Process (HPPP) , which have been proved to be as accurate as the hexagonal grid model. Furthermore, the Fractional Frequency Reuse (FFR) was analyzed using PPPs in order to tackle the strong interference from nearby macro BSs [24, 25]. In these works of , all analytical results have been verified via Monte Carlo methods. Note that all these works just adopt a simple resource allocation scheme, namely, the entire bandwidth of each BS is time-shared for its users. We consider a multichannel downlink based on orthogonal frequency division multiple access (OFDMA) technique, where each user is only served by one sub-channel. In this paper, we use the metric of average user rate, which can reflect the high spectral efficiency obtained from the proposed SFR scheme.

Based on the approach described above, the contributions of the paper are summarized as follows:

1. 1)

Based on the cell load model, we first derive the coverage probability of the network without SFR scheme, which can be used for the derivation of coverage probability and average rate of the network with SFR scheme.

2. 2)

Next, we derive the evaluation for network performance of the proposed SFR scheme is performed in terms of average user rate incorporated with the random resource allocation and cell load. For special case, the expressions of coverage probability are simple closed-form.

3. 3)

We compare the proposed SFR scheme with FRP scheme. Moreover, we comprehensively analyze the average user rate under different parameters by varying the association bias, resource partitioning factor, power control factor, and SINR thresholds. Then, we show the impact of the aforementioned parameters on the average user rate and investigate the optimal combination finally.

4. 4)

We show that the proposed SFR scheme is promising for improving the average user rate while considering the cell load.

The remainder of the paper is organized as follows: the system model, user association scheme, cell load analysis, and SFR scheme are presented in Section 2. Mathematical coverage probability and average user rate expressions are derived in Section 3. In Section 4, numerical results on the performance evaluation for the proposed SFR scheme are analyzed. Furthermore, the comparison with FRP scheme is presented. Finally, the paper is concluded in Section 5.

## System model

### Two-tier cellular network model

We consider a two-tier downlink cellular network based on OFDMA technique, i.e., the intra-cell interference is not considered. Without any loss of generality, let the first tier be macro (higher-power BS) tier, while let the second tier be pico (lower-power BS) tier. The locations of the BSs in jth tier are modeled as a two-dimensional HPPP Φ j with density λ j . Furthermore, the users are located according to another HPPP Φ u with density λ u , which is independent of {Φ j } j=1,2. For the co-channel deployment of the network, both the network tiers share the same set of available sub-channels (denoted by C). Then, let |C| denote the total number of available sub-channels. Moreover, every BS in jth tier transmits with same power $$P_{j}^{T}$$ and thus, the power per sub-channel of a BS in jth tier is kept constant at $${P_{j}} = {{P_{j}^{T}} \left / {\left | \mathbf{C} \right |}\right.}$$. For tractability, the standard power loss propagation model is applied in both tiers with the same path loss exponent α>2. As far as random channel fluctuations, Rayleigh fading with mean 1 (denoted as H x exp(1)) is applied at each channel. The noise is assumed to be additive with power σ 2. All macro and pico BSs are assumed to be open access in this paper. That is, the number of users served by each BS is unlimited.

### User association

We consider a user association that each user chooses its serving BS based on maximum biased-received-power (BRP) (termed biased user association) [79, 1517]. Denote R j as the distance of the typical user from its nearest BS of jth tier. A typical user at the origin is associated with the nearest BS in the kth tier if

$$k = \arg\max \left\{ {{P_{1}}R_{1}^{- \alpha },{P_{2}}BR_{2}^{- \alpha }} \right\},$$
(1)

where B≥0 dB is the association bias for tier 2 (pico tier). For simplicity, it is assumed that the association bias for tier 1 (macro tier) is unity in this paper. As the association bias B increases, more macro users will be offloaded to the corresponding pico BSs. Here, we let $${\mathcal {U}_{k}}$$ denote the set of users in the kth tier, which satisfies $${\mathcal {U}_{1}} \cup {\mathcal {U}_{2}} = \mathcal {U}$$.

Based on (1), the probabilities for the typical user to associate with macro and pico tiers, denoted by $${\mathcal {A}_{1}}$$ and $$\mathcal {A}_{2}$$, can be derived as in , i.e.,

$$\mathcal{A}_{1} = \frac{\lambda_{1}}{\lambda_{1} + \lambda_{2} \left(\frac{P_{2}B}{P_{1}}\right)^{2/\alpha}}\ \text{and}\ \mathcal{A}_{2} = \frac{\lambda_{2}}{\lambda_{1} \left(\frac{\mathrm{P}_{1}}{\mathrm{P}_{2}\mathrm{B}}\right)^{2/\alpha}+ \lambda_{2}}.$$
(2)

Using (2), the probability density function (PDF) of the distance between the typical user and its serving BS, denoted by $${f_{{X_{k}}}}\left (x \right)$$, is also given in  as

$$f_{X_{k}} (x) = \frac{2\pi\lambda_{k}x}{\mathcal{A}_{k}} \text{exp} \left(\frac{-\pi\lambda_{k}x^{2}}{\mathcal{A}_{k}}\right).$$
(3)

Since all BSs in every tier transmit with the same transmit power, each user belongs to macro tier or pico tier is always associated with the nearest macro BS or pico BS based on the maximum BRP. The two-tier cellular network coverage region will be constituted by two independent Voronoi tessellations. Furthermore, the size of each Voronoi cell is an i.i.d. random variable . On the basis of statistical property of the PPP, the probability mass function (PMF) of the number of users associated with a randomly chosen kth tier BS can be expressed as

{} \begin{aligned} \mathbb{P} \left(N_{k}= n \right) = \frac{3.5^{3.5}}{n!} \frac{\Gamma(n + 3.5)}{\Gamma(3.5)} \left(\frac{\lambda_{u}\mathcal{A}_{k}}{\lambda_{k}}\right)^{n} \left(3.5 + \frac{\lambda_{u}\mathcal{A}_{k}}{\lambda_{k}}\right)^{-(n + 3.5)}, \end{aligned}
(4)

where $$\Gamma \left (z \right) = \int _{0}^{\infty } {{t^{z - 1}}{e^{- t}}dt}$$ is the standard gamma function. The statistical property of the number of users is crucial for calculating a cell load.

### The main results of network without SFR

Before discussing the proposed SFR scheme, we firstly focus on the network’s SINR distribution without SFR, which is extremely essential for the classification of users. Following the analysis [22, 23], the main results for the biased user association are briefly presented here for the purpose of analytical evaluating the proposed SFR scheme.

#### Resource allocation and load statistics

For the sake of keeping simplicity and tractability of the PPP model, we adopt a simple random resource allocation scheme for both macro and pico users. That is, each BS, independent of the other BSs, randomly and uniformly selects one sub-channel for each of its users. If the number of users in a serving BS is greater than the number of available sub-channels, the resources can be equally allocated in time-division way. Hence, together with the use of OFDMA, the intra-cell interference can be ignored. With regard to other sophisticated resource allocation schemes like opportunistic and fair resource allocation, they are not analytically tractable and will be taken in consideration for our future work.

Different from the assumptions in  that all BSs of each tier are always active, we utilize a more practical load model that a BS is active when it has at least one user to serve. When a BS has no user to serve, its corresponding sub-channel set C will be left idle. Furthermore, it is also assumed that each BS has full buffer traffic downlink transmission for its each user. Hence, a BS is the source of interference when the BS must be active and simultaneously use the same sub-channel.

To obtain the SINR statistics on a given sub-channel, [15, 19] have derived the probability that a typical BS of each tier accesses a given sub-channel, which depends on the PMF of the number of users associated with that BS.

### Lemma 1

Let ρ k denote the probability that a typical BS in the kth tier accesses a given sub-channel from the set C. Then,

$${\rho_{k}} = 1 - \sum\limits_{n = 0}^{\left| \mathbf{C} \right| - 1} {\left({1 - \frac{n}{{\left| \mathbf{C} \right|}}} \right)} \mathbb{P} \left({{N_{k}} = n} \right).$$
(5)

Note that, we can term ρ k as the load of a typical BS in the kth tier. It also can be viewed as the probability that a typical BS becomes the source of interference because the BSs use the same sub-channel. Therefore, the macro BSs and pico BSs using the same sub-channel will form two independent homogeneous PPPs Ψ 1 of density λ 1 ρ 1 and Ψ 2 of density λ 2 ρ 2, respectively. In other words, the interfering sets Ψ 1 and Ψ 2 are independent thinning of the original PPPs Φ 1 and Φ 2, respectively, with retention probabilities ρ 1 and ρ 2, respectively [ 15 , 19 ].

#### SINR distribution

$${{\text{SINR}_{k}} = \frac{{{P_{k}}{H_{x}}{x^{- \alpha }}}}{{{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}} + {\sigma^{2}}} }}},$$
(6)

where $${I_{x,j}} = {P_{j}}{\sum \nolimits }_{y \in {\Psi _{j}}\backslash {b_{k}}} {{H_{y}}{{\left \| y \right \|}^{- \alpha }}}$$ is the cumulative interference from all the loaded BSs in the jth tier (except the user’s serving BS in the kth tier), and H x is the channel fading gain from the serving BS b k at a distance x.

Hence, the SINR distribution of a typical user can be thought of equivalently as the probability that the typical user can achieve a target SINR threshold T, i.e., the conditional coverage probability. The conditional coverage probability when a typical user $$u \in {\mathcal {U}_{k}}$$ is defined as

$${p_{c,k}}\left(T \right) \buildrel \Delta \over = \mathbb{P}\left({\text{SINR} > T\left| {u \in {\mathcal{U}_{k}}} \right.} \right).$$
(7)

### Theorem 1

Following the law of total probability, the coverage probability without SFR of a typical user $$u \in \mathcal {U}$$ is

$${p_{c}}\left(T \right) = {{\mathcal A}_{1}}{p_{c,1}}\left(T \right) + {{\mathcal A}_{2}}{p_{c,2}}\left(T \right),$$
(8)

where the conditional coverage probabilities without SFR are given by (9) and (10), $$Q\left ({a,b,c,d} \right) = {c^{{2 / b}}} + {a^{{2 / b}}} d\int _{{{\left ({\frac {c}{a}} \right)}^{{2 / b}}}}^{\infty } {\frac {{du}}{{1 + {u^{{b / 2}}}}}}$$, and $${\text {SNR}_{k}}\left (x \right) = \frac {{{P_{k}}{x^{- \alpha }}}}{{{\sigma ^{2}}}}$$.

{} \begin{aligned} {P}_{c,1} (T) &= \frac{2\pi\lambda_{1}}{\mathcal{A}_{1}} \int_{0}^{\infty}\ x\ \text{exp} \left(-\frac{T}{\text{SNR}_{1}(x)}\right)\\ &\quad\,\, \text{exp} \left\{-\pi\lambda_{1}Q(T,\alpha,1,\rho_{1})x^{2}\right.\\&\qquad\quad \,-\left. \pi\lambda_{2}(P_{2}/P_{1})^{2/\alpha}Q(T, \alpha, B, \rho_{2}) x^{2}\right\}dx \end{aligned}
(9)
{}{\begin{aligned} {P}_{c,2} (T) &= \frac{2\pi\lambda_{2}}{\mathcal{A}_{2}} \int_{0}^{\infty}\ x\ \text{exp} \left(-\frac{T}{\text{SNR}_{2}(x)}\right)\\&\quad\,\, \text{exp} \left\{- \pi\lambda_{1}(P_{1}/P_{2})^{2/\alpha}Q(T,\alpha,1/B,\rho_{1})x^{2}\right.\\&\qquad\quad\, -\left.\pi\lambda_{2}Q(T, \alpha, 1, \rho_{2}) x^{2}\right\}dx \end{aligned}}
(10)

### Proof

The proof is given in Appendix A. □

Unlike the results in , the coverage probability is dependent of ρ k , i.e., the load of each BS. Additionally, the conditional coverage probabilities without SFR are of special importance with respect to the proposed SFR scheme.

### SFR scheme

In this paper, the SFR scheme proposed by [17,18,20] is applied to our two-tier cellular network, as shown in Fig. 1.

#### Resource partitioning and power control

Similarly, both the network tiers are co-channel deployment. However for each cell, the entire available sub-channel set C is divided into two different subsets C 1 and C 2 with size |C 1| and |C 2|, respectively. The two different sunsets C 1 and C 2 have no intersection, i.e., C 1C 2= and C=C 1C 2. Therefore, a resource partitioning factor η is defined as η=|C 1|/|C|(0≤η≤1).

In the proposed SFR scheme, a typical BS classifies users with average SINR as two types of users: interior users and edge users. Instead of a geographic classification criterion, the SINR threshold can more adequately capture the randomness: the locations of the BSs and users [24,25]. Let T FR,k denote the SFR threshold of a typical BS in the kth tier. With SFR, each user calculates its SINR to the serving BS in the kth tier, and if it is less than the threshold T FR,k , then the user is an edge user, and otherwise the user is an interior user.

In this paper, we consider the following random resource allocation scheme. Let η be the fraction of resources (namely, C 1) allocated to the macro interior users and pico edge users. The remaining 1−η fraction of resources (namely, C 2) are allocated to the macro edge users and pico interior users. Nevertheless, the macro edge users and pico edge users, especially the offloaded macro users, usually suffer severe interference from the neighboring macro BSs. To accomplish this, a power control factor β(0<β<1) is introduced to the sub-channel set C 1 used by the macro BSs, i.e., $$P_{1}^{\mathrm {i}} = \beta {P_{1}}$$ and $$P_{1}^{\mathrm {e}} = {P_{1}}$$, where $$P_{1}^{\mathrm {i}}$$ is the transmit power of macro BSs for the interior users and $$P_{1}^{\mathrm {e}}$$ is the transmit power of macro BSs for the edge users. For all the pico users, the transmit power keeps full power transmission, the same as the macro edge users. Therefore, both high spectral efficiency and good user experience of the pico edge users can be achieved.

#### User association and load statistics

According to the proposed SFR scheme and user association scheme above, a typical user $$u \in \mathcal {U}$$ can lie in the following four disjoint sets:

\begin{aligned} u \in \left\{ \begin{array}{l} U_{1}^{\mathrm{i}}, \text{if}\ k = 1, P_{1} R_{1}^{-\alpha} \geq P_{2} B\, R^{-\alpha}_{2}\ \&\ \text{SINR}_{1} \geq T_{\text{FR}, 1}\\ U_{1}^{\mathrm{e}}, \text{if}\ k = 1, P_{1} R_{1}^{-\alpha} \geq P_{2} B\, R^{-\alpha}_{2}\ \&\ \text{SINR}_{1} < T_{\text{FR}, 1}\\ U_{2}^{\mathrm{i}}, \text{if}\ k = 2, P_{1} R_{1}^{-\alpha} < P_{2} B\, R^{-\alpha}_{2}\ \&\ \text{SINR}_{2} \geq T_{\text{FR}, 2}\\ U_{2}^{\mathrm{e}}, \text{if}\ k = 2, P_{1} R_{1}^{-\alpha} < P_{2} B\, R^{-\alpha}_{2}\ \&\ \text{SINR}_{2} < T_{\text{FR}, 2}\\ \end{array} \right., \end{aligned}
(11)

where $$\mathcal {U}_{1}^{\mathrm {i}}$$ is the set of interior users associated with the kth tier BSs and $$\mathcal {U}_{1}^{\mathrm {e}}$$ is the set of edge users associated with the kth tier BSs. Clearly, $${\mathcal {U}_{1}} \buildrel \Delta \over = \mathcal {U}_{1}^{\mathrm {i}} \cup \mathcal {U}_{1}^{\mathrm {e}}$$ is the set of macro users, $${\mathcal {U}_{2}} \buildrel \Delta \over = \mathcal {U}_{2}^{\mathrm {i}} \cup \mathcal {U}_{2}^{\mathrm {e}}$$ is the set of pico users, and $$\mathcal {U}_{1}^{\mathrm {i}} \cup \mathcal {U}_{1}^{\mathrm {e}} \cup \mathcal {U}_{2}^{\mathrm {i}} \cup \mathcal {U}_{2}^{\mathrm {e}} = \mathcal {U}$$.

For a randomly chosen user in Φ u , it will exactly belong to specific one of above four sets, according to the user association strategy in (11). Then, the probabilities that a randomly chosen user belongs to the sets $$\mathcal {U}_{k}^{\mathrm {i}}$$ and $$\mathcal {U}_{k}^{\mathrm {e}}$$ are $${{\mathcal A}_{k}}{p_{c,k}}\left ({{T_{\text {FR},k}}} \right)$$ and $${{\mathcal A}_{k}}\left ({1 - {p_{c,k}}\left ({{T_{\text {FR},k}}} \right)} \right)$$, which also can be interpreted as the average fraction of users belonging to the sets $$\mathcal {U}_{k}^{\mathrm {i}}$$ and $$\mathcal {U}_{k}^{\mathrm {e}}$$. Irrespective of the exact distribution of user locations, the numbers of interior and edge users in a typical BS are significant for characterizing the cell load. To make the proposed framework analytically tractable, the sets $$\mathcal {U}_{k}^{\mathrm {i}}$$ and $$\mathcal {U}_{k}^{\mathrm {e}}$$ can be equivalently modeled as independent homogeneous PPPs with densities $${{\mathcal A}_{k}}{p_{c,k}}\left ({{T_{\text {FR},k}}} \right)$$ and $${{\mathcal A}_{k}}\left ({1 - {p_{c,k}}\left ({{T_{\text {FR},k}}} \right)} \right)$$, respectively . Based on the same lines as the derivation in (4), their PMFs are given in the following lemma.

### Lemma 2

Let $$N_{k}^{\mathrm {i}}$$ and $$N_{k}^{\mathrm {e}}$$ denote the numbers of interior and edge users associated with a BS in the kth tier, respectively. Their PMFs are given by (12) and (13).

$$\begin{array}{@{}rcl@{}} {}\mathbb{P} \left({N_{k}^{\mathrm{i}} = n} \right) &=& \frac{{{{3.5}^{3.5}}}}{{n!}}\frac{{\Gamma \left({n + 3.5} \right)}}{{\Gamma \left({3.5} \right)}}{\left({\frac{{{\lambda_{u}}{{\mathcal A}_{k}}{p_{c,k}}\left({{T_{\text{FR},k}}} \right)}}{{{\lambda_{k}}}}} \right)^{n}}\\&& {\left({3.5 + \frac{{{\lambda_{u}}{{\mathcal A}_{k}}{p_{c,k}}\left({{T_{\text{FR},k}}} \right)}}{{{\lambda_{k}}}}} \right)^{- \left({n + 3.5} \right)}} \end{array}$$
(12)
$$\begin{array}{@{}rcl@{}} \mathbb{P} \left({N_{k}^{\mathrm{e}} = n} \right) &=& \frac{{{{3.5}^{3.5}}}}{{n!}}\frac{{\Gamma \left({n + 3.5} \right)}}{{\Gamma \left({3.5} \right)}}{\left({\frac{{{\lambda_{u}}{{\mathcal A}_{k}}\left({1 - {p_{c,k}}\left({{T_{\text{FR},k}}} \right)} \right)}}{{{\lambda_{k}}}}} \right)^{n}} \\&& {\left({3.5 + \frac{{{\lambda_{u}}{{\mathcal A}_{k}}\left({1 - {p_{c,k}}\left({{T_{\text{FR},k}}} \right)} \right)}}{{{\lambda_{k}}}}} \right)^{- \left({n + 3.5} \right)}} \end{array}$$
(13)

### Lemma 3

Let $$\rho _{k}^{\mathrm {i}}$$ and $$\rho _{k}^{\mathrm {e}}$$ denote the probabilities that a typical BS in the kth tier accesses a given sub-channel from the allocated subsets C 1 and C 2, respectively. Then,

\begin{aligned} \rho_{1}^{\mathrm{i}} = 1 - \sum\limits_{n = 0}^{\left| {{\mathbf{C}_{1}}} \right| - 1} {\left({1 - \frac{n}{{\left| {{\mathbf{C}_{1}}} \right|}}} \right)} \mathbb{P} \left({N_{1}^{\mathrm{i}} = n} \right), \end{aligned}
(14)
\begin{aligned} \rho_{1}^{\mathrm{e}} = 1 - \sum\limits_{n = 0}^{\left| {{\mathbf{C}_{2}}} \right| - 1} {\left({1 - \frac{n}{{\left| {{\mathbf{C}_{2}}} \right|}}} \right)} \mathbb{P} \left({N_{1}^{\mathrm{e}} = n} \right), \end{aligned}
(15)
\begin{aligned} \rho_{2}^{\mathrm{i}} = 1 - \sum\limits_{n = 0}^{\left| {{\mathbf{C}_{2}}} \right| - 1} {\left({1 - \frac{n}{{\left| {{\mathbf{C}_{2}}} \right|}}} \right)} \mathbb{P} \left({N_{2}^{\mathrm{i}} = n} \right), \end{aligned}
(16)
\begin{aligned} \rho_{2}^{\mathrm{e}} = 1 - \sum\limits_{n = 0}^{\left| {{\mathbf{C}_{1}}} \right| - 1} {\left({1 - \frac{n}{{\left| {{\mathbf{C}_{1}}} \right|}}} \right)} \mathbb{P} \left({N_{2}^{\mathrm{e}} = n} \right). \end{aligned}
(17)

As our explanation earlier, we can term $$\rho _{k}^{\mathrm {i}}$$ and $$\rho _{k}^{\mathrm {e}}$$ as the interior load and edge load of a typical BS in the kth tier, respectively. For a macro interior user served on a sub-channel from the subset C 1, it will be interfered by the interior loaded macro BSs and the edge loaded pico BSs. Therefore, its interfering sets $$\Xi _{1}^{\mathrm {i}}$$ and $$\Xi _{2}^{\mathrm {e}}$$ are independent thinning of the original PPPs Φ 1 and Φ 2, respectively, with retention probabilities $$\rho _{1}^{\mathrm {i}}$$ and $$\rho _{2}^{\mathrm {e}}$$, respectively. Similarly, the interfering sets $$\Xi _{k}^{\mathrm {i}}$$ and $$\Xi _{k}^{\mathrm {e}}$$ are two independent homogeneous PPPs with densities $${\lambda _{k}}\rho _{k}^{\mathrm {i}}$$ and $${\lambda _{k}}\rho _{k}^{\mathrm {e}}$$, respectively. The interfering sets of each user can be obtained as given in Table 1.

By employing the user association scheme and SFR scheme described above, the downlink received SINR of a typical user $$u \in \mathcal {U}$$ located at the origin is given by

$${} \begin{array}{l} \text{SINR} = \mathbf{1}\left({u \in \mathcal{U}_{1}^{\mathrm{i}}} \right)\frac{{\beta {P_{1}}{H_{x}}{x^{- \alpha }}}}{{\beta I_{x,1}^{\mathrm{i}} + I_{x,2}^{\mathrm{e}} + {\sigma^{2}}}} + \mathbf{1}\left({u \in \mathcal{U}_{1}^{\mathrm{e}}} \right)\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{I_{x,1}^{\mathrm{e}} + I_{x,2}^{\mathrm{i}} + {\sigma^{2}}}}\\ \qquad \quad + \mathbf{1}\left({u \in \mathcal{U}_{2}^{\mathrm{i}}} \right)\frac{{{P_{2}}{H_{x}}{x^{- \alpha }}}}{{I_{x,1}^{\mathrm{e}} + I_{x,2}^{\mathrm{i}} + {\sigma^{2}}}} + \mathbf{1}\left({u \in \mathcal{U}_{2}^{\mathrm{e}}} \right)\frac{{{P_{2}}{H_{x}}{x^{- \alpha }}}}{{\beta I_{x,1}^{\mathrm{i}} + I_{x,2}^{\mathrm{e}} + {\sigma^{2}}}} \end{array},$$
(18)

where 1(A) is the indicator function that takes the value 1 if the event A is true, H x is the channel fading gain from the serving BS b k at a distance x, and $$I_{x,j}^{\mathrm {i}} = {P_{j}}{\sum \nolimits }_{y \in \Xi _{j}^{\mathrm {i}}\backslash {b_{k}}} {{H_{y}}{{\left \| y \right \|}^{- \alpha }}}$$ and $$I_{x,j}^{\mathrm {e}} = {P_{j}}{\sum \nolimits }_{y \in \Xi _{j}^{\mathrm {e}}\backslash {b_{k}}} {{H_{y}}{{\left \| y \right \|}^{- \alpha }}}$$ are the cumulative interference from the interior and edge loaded BSs in the jth tier, respectively.

## Average user rate

This section is our main technical part. We first derive the general coverage probability for the proposed SFR scheme. Then, the methods of derivation are subsequently used for the average user rate. Moreover, we present a special case where α=4 and σ 2=0, representing for the modern cellular networks. For this case, the expressions of coverage probability reduce to simple closed-form, which can provide clear insight into the performance analysis of each user.

### Coverage probability

In the context of this paper, the conditional coverage probabilities $$p_{\text {FFR},k}^{\mathrm {i}}\left (T \right)$$ and $$p_{\text {FFR},k}^{\mathrm {e}}\left (T \right)$$ can be defined, respectively, as

$$p_{\text{FFR},k}^{\mathrm{i}}(T) \triangleq \mathbb{P} \left(\text{SINR} > T\left| u \in \mathcal{U}^{\mathrm{i}}_{k} \right.\right),$$
(19)
$$p_{\text{FFR},k}^{\mathrm{e}}(T) \triangleq \mathbb{P} \left(\text{SINR} > T\left| u \in \mathcal{U}^{\mathrm{e}}_{k} \right.\right).$$
(20)

Following from the law of total probability, the coverage probability of a typical user following the user association strategy in (11) is given by (21).

\begin{aligned} {p_{c,\text{FFR}}}\left(T \right) &= \sum\limits_{k = 1}^{2} {{\mathcal{A}}_{k}}\left({p_{c,k}}\left({{T_{\text{FR},k}}} \right)p_{\text{FFR},k}^{\mathrm{i}}\left(T \right)\right.\\ &\quad +\left. \left({1 - {p_{c,k}}\left({{T_{\text{FR},k}}} \right)} \right)p_{\text{FFR},k}^{\mathrm{e}}\left(T \right) \right) \end{aligned}
(21)
$$p_{\text{FFR}, 1}^{\mathrm{i}}(T) = \frac{\begin{array}{rl} \frac{2\pi\lambda1}{\mathcal{A}1} \int^{\infty}_{0} x\ \text{exp}\ \left(-\frac{T}{\text{SNR}_{1}(x,\beta)}-\frac{T_{\text{FR},1}}{\text{SNR}_{1}(x,1)}\right) \text{exp} \left\{-\pi\lambda_{1}x^{2}\left[1+2\xi\left(\rho^{\mathrm{i}}_{1},\rho_{1},T,T_{\text{FR},1},\alpha\right)\right]\right.&\\ \left.-\pi\lambda_{2}\left(\frac{P_{2}B}{P_{1}}\right)^{\frac{2}{\alpha}} x^{2} \left[1+2\xi\left(\rho^{\mathrm{e}}_{2}, \rho_{2}, (\beta B)^{-1} T, B^{-1} T_{\text{FR},1},\alpha\right)\right]\right\}&dx \end{array}}{p_{c,1}(T_{\text{FR},1})}$$
(22)
$$p^{\mathrm{i}}_{\text{FFR},1} (T)=\frac{\lambda_{1}\left(1 + \rho_{1}\sqrt{T_{\text{FR},1}} \ \text{arctan}(\sqrt{T_{\text{FR},1}})\right)+\lambda_{2}\sqrt{P_{2}/P_{1}}\left(\sqrt{B} + \rho_{2}\sqrt{T_{\text{FR},1}}\ \text{arctan}\ \left(\sqrt{T_{\text{FR},1}/B}\right)\right)}{\lambda_{1}\left(1+2\xi\left(\rho^{\mathrm{i}}_{1}, \rho_{1},T,T_{\text{FR},1},4\right)\right)+\lambda_{2}\sqrt{P_{2}B/P_{1}}\left(1+2\xi\left(\rho_{2}^{\mathrm{e}},\rho_{2},(\beta B)^{-1}T,B^{-1}T_{\text{FR},1},4\right)\right)}$$
(23)
$$\begin{array}{*{20}l} p_{\text{FFR}, 1}^{\mathrm{e}}(T) =& \frac{\frac{2\pi\lambda_{1}}{\mathcal{A}_{1}}\int^{\infty}_{0}\ x\ \text{exp}\left(-\frac{T}{\text{SNR}_{1}(x,1)}\right)\text{exp}\left\{-\pi\lambda_{1}Q\left(T,\alpha,1,\rho^{\mathrm{e}}_{1}\right)x^{2}-\pi\lambda_{2}\left(\frac{P_{2}}{P_{1}}\right)^{\frac{2}{\alpha}}Q\left(T,\alpha,B,\rho^{i}_{2}\right)x^{2}\right\}dx}{1{-p_{c,1}}(T_{\text{FR},1})}\\ &-\frac{\begin{array}{rl} \frac{2\pi\lambda1}{\mathcal{A}_1} \int^{\infty}_{0} x\ \text{exp}\ \left(-\frac{T+T_{\text{FR},1}}{\text{SNR}_{1}(x,1)}\right) \text{exp} \left\{-\pi\lambda_{1}x^{2}\left[1+2\xi\left(\rho^{\mathrm{e}}_{1},\rho_{1},T,T_{\text{FR},1},\alpha\right)\right]{\vphantom{\left(\frac{P_{2}B}{P_{1}}\right)^{\frac{2}{\alpha}}}}\right.&\\ \left.-\pi\lambda_{2}\left(\frac{P_{2}B}{P_{1}}\right)^{\frac{2}{\alpha}} x^{2} \left[1+2\xi\left(\rho^{\mathrm{i}}_{2}, \rho_{2}, B^{-1} T, B^{-1} T_{\text{FR},1},\alpha\right)\right]\right\}&dx \end{array}}{1{-p_{c,1}}(T_{\text{FR},1})} \end{array}$$
(24)
$$\begin{array}{ll} p_{\text{FFR},1}^{\mathrm{e}}\left(T \right) &= \frac{\left(\lambda_{1}+\lambda_{2}\sqrt{P_{2} B/P_{1}}\right)\left[\lambda_{1}\left(1+\rho_{1}\sqrt{T_{\text{FR},1}}\ \text{arctan}\left(\sqrt{T_{\text{FR},1}}\right) \right)+\lambda_{2}\sqrt{P_{2}/P_{1}}\left(\sqrt{B}+\rho_{2}\sqrt{T_{\text{FR},1}}\ \text{arctan}\left(\sqrt{T_{\text{FR},1}/B}\right)\right)\right]}{\left(1-\lambda_{1}-\lambda_{2}\sqrt{P_{2} B/P_{1}}\right)\left[\lambda_{1}\left(1+\rho_{1}^{\mathrm{e}}\sqrt{T} \ \text{arctan}\ (\sqrt{T})\right)+\lambda_{2}\sqrt{P_{2}/P_{1}}\left(\sqrt{B}+\rho^{\mathrm{i}}_{2}\sqrt{T}\ \text{arctan}\left(\sqrt{T/B}\right)\right)\right]}\\ &\quad-\frac{\left(\lambda_{1}+\lambda_{2}\sqrt{P_{2}B/P_{1}}\right)\left[\lambda_{1}\left(1+\rho_{1}\sqrt{T_{\text{FR},1}}\ \text{arctan}\left(\sqrt{T_{\text{FR},1}}\right) \right)+\lambda_{2}\sqrt{P_{2}/P_{1}}\left(\sqrt{B}+\rho_{2}\sqrt{T_{\text{FR},1}}\ \text{arctan}\left(\sqrt{T_{\text{FR},1}/B}\right)\right)\right]}{\left(1-\lambda_{1}-\lambda_{2}\sqrt{P_{2} B/P_{1}}\right)\left[\lambda_{1}\left(1+2\xi\left(\rho^{\mathrm{e}}_{1}, \rho_{1},T,T_{\text{FR},1},4\right)\right)+\lambda_{2}\sqrt{P_{2}B/P_{1}}\left(1+2\xi\left(\rho^{i}_{2},\rho_{2},B^{-1}T,B^{-1}T_{\text{FR},1},4\right)\right)\right]} \end{array}$$
(25)
$$p_{\text{FFR},{2}}^{\mathrm{i}}\left(T \right) = \frac{\begin{array}{l} \frac{2\pi {\lambda_{2}}}{{\mathcal A}_{2}}\int_{0}^{\infty} x\exp \left(-\frac{T + T_{\text{FR},2}}{\text{SNR}_{2}\left({x,1}\right)} \right)\exp \left\{ \vphantom{\left(\frac{P_{1}}{P_{2}B} \right)}\!\!- \pi {\lambda_{2}}{x^{2}}\left[1 + 2\xi \left({\rho_{2}^{\mathrm{i}},{\rho_{2}},T,{T_{\text{FR},2}},\alpha} \right)\right] \right. \\ \qquad \qquad \qquad \qquad\left. - \pi {\lambda_{1}}\left(\frac{P_{1}}{P_{2}B} \right)^{\frac{2}{\alpha}}{x^{2}}\left[1 + 2\xi \left(\rho_{1}^{\mathrm{e}},{\rho_{1}},BT,BT_{\text{FR},2},\alpha \right) \right] \right\}dx \end{array}}{p_{c,2}\left(T_{\text{FR},2}\right)}$$
(26)

As discussed in Section 2, we observe that the conditional coverage probability of a typical user depends on two SINR thresholds [24,25], first the SINR threshold (i.e., T FR,k ) on the allocated a sub-channel from the entire available sub-channel set C to determine its status (interior or not), and second the actual SINR threshold (i.e., T) on the newly allocated a sub-channel from the subset C 1 or C 2 to determine whether it is covered or not. Actually, these two SINR thresholds are correlated because the interference may be generated by the same set of BSs, which makes our analysis challenging. The following theorems give the conditional coverage probabilities for a typical user under different status.

### Theorem 2

(Macro tier, interior user): the coverage probability of the macro interior user is given by (22), where $$\text {SNR}_{k}\left ({x,a} \right) = \frac {{a{P_{k}}{x^{- \alpha }}}}{{{\sigma ^{2}}}}$$, p c,1(T FR,1) is given by (9), and $$\xi \left ({a,b,c,d,e} \right) = \int _{1}^{\infty } \left [ 1 - \left (1 - a\left ({1 - \frac {1}{{1 + c{v^{- e}}}}} \right) \right) \left (1 - b\left ({1 - \frac {1}{{1 + d{v^{- e}}}}} \right) \right) \right ]vdv$$.

### Proof

The proof is given in Appendix B. □

As we can see, ξ(a,b,c,d,e) is similar to ρ j Z(a,b,c) given by previous results in Theorem 1, they are different because of the dependence of two SINR thresholds. Moreover, ξ(a,b,c,d,e) can efficiently capture the disparity of the intra-tier and inter-tier interference before and after the proposed SFR scheme is applied.

Now, we turn our attention to the special case where α=4 and σ 2=0, which is significant in practice that is widely applied in lots of literatures . It is noted that the typical HCNs are interference-limited which the noise can be ignored compared to the interference. Furthermore, the special case can be considered as the scenario corresponding to an interference-limited urban cellular network [24,25,27], where FFR has been generally applied. In the special case, the expression (22) will be further simplified to a simple closed-form.

### Corollary 1

For α=4 and σ 2=0, the coverage probability of the macro interior user is given by (23), where $$\xi \left ({\rho _{1}^{\mathrm {i}},{\rho _{1}},T,{T_{\text {FR},1}},4} \right) =$$

$$\begin{array}{l} \frac{\rho_{1}^{\mathrm{i}}\sqrt{T} \arctan \left(\sqrt{T}\right)\left({T_{\text{FR},1}} - T - {\rho_{1}}{T_{\text{FR},1}}\right)}{2\left({{T_{\text{FR},1}} - T} \right)}\\ + \frac{{\rho_{1}}\sqrt {{T_{\text{FR},1}}} \arctan \left({\sqrt {{T_{\text{FR},1}}}} \right)\left({{T_{\text{FR},1}} - T + \rho_{1}^{\mathrm{i}}T} \right)}{2\left({{T_{\text{FR},1}} - T} \right)} \end{array},$$

and $$\xi \left ({\rho _{2}^{\mathrm {e}},{\rho _{2}},{{\left ({\beta B} \right)}^{- 1}}T,{B^{- 1}}{T_{\text {FR},1}},4} \right) =$$

$$\begin{array}{l} \frac{\rho_{2}^{\mathrm{e}}\sqrt {{{\left({\beta B} \right)}^{- 1}}T} \arctan \left({\sqrt {{{\left({\beta B} \right)}^{- 1}}T}} \right)\left({{T_{\text{FR},1}} - {\beta^{- 1}}T - {\rho_{2}}{T_{\text{FR},1}}} \right)}{2\left({{T_{\text{FR},1}} - {\beta^{- 1}}T} \right)}\\ + \frac{{\rho_{2}}\sqrt {{B^{- 1}}{T_{\text{FR},1}}} \arctan \left({\sqrt {{B^{- 1}}{T_{\text{FR},1}}}} \right)\left({{T_{\text{FR},1}} - {\beta^{- 1}}T + \rho_{2}^{\mathrm{e}}{\beta^{- 1}}T} \right)}{2\left({{T_{\text{FR},1}} - {\beta^{- 1}}T} \right)} \end{array}.$$

### Proof

When α=4, we have

$$Q\left({a,4,c,d} \right) = \sqrt{c} + \sqrt{ad} \arctan \left(\sqrt{a/c}\right).$$

Then, ξ(a,b,c,d,4)=

$${} \begin{array}{c} \int_{1}^{\infty} {\left[\! {1 - \left(\! {1 - a\left({1 - \frac{1}{{1 + c{v^{- 4}}}}} \!\right)} \right) \left({1 - b\left({1 - \frac{1}{{1 + d{v^{- 4}}}}} \right)} \right)} \!\right]vdv} \\ = \frac{{a\sqrt{c} \arctan \left({\sqrt{c}} \right)\left({d - c - bd} \right) + b\sqrt{d} \arctan \left({\sqrt{d}} \right)\left({d - c + ac} \right)}}{{2\left({d - c} \right)}} \end{array}.$$

Combining with (22) gives the desired result.

We see that the coverage probability of a macro interior user is a function of the SINR threshold T, the SFR threshold T FR,1, the power control factor β, the association bias B, and the load ρ 1, $$\rho _{1}^{\mathrm {i}}$$, ρ 2 and $$\rho _{2}^{\mathrm {e}}$$. Furthermore, the expression will have an indeterminate form when T=T FR,1. Let limit TT FR,1, the expressions will be simplified to $$\xi \left ({\rho _{1}^{\mathrm {i}},{\rho _{1}},T,{T_{\text {FR},1}},4} \right)=$$

$$\frac{{\left({2{\rho_{1}} + 2\rho_{1}^{\mathrm{i}} - {\rho_{1}}\rho_{1}^{\mathrm{i}}} \right)\sqrt {{T_{\text{FR},1}}} \arctan \left({\sqrt {{T_{\text{FR},1}}}} \right)}}{4} + \frac{{{\rho_{1}}\rho_{1}^{\mathrm{i}}{T_{\text{FR},1}}}}{{4\left({{T_{\text{FR},1}} + 1} \right)}},$$

and $$\xi \left ({\rho _{2}^{\mathrm {e}},{\rho _{2}},{{\left ({\beta B} \right)}^{- 1}}T,{B^{- 1}}{T_{\text {FR},1}},4} \right) =$$

$$\begin{array}{c} \frac{{\rho_{2}^{\mathrm{e}}\left({1 - {\beta^{- 1}} - {\rho_{2}}} \right)\sqrt {{{\left({\beta B} \right)}^{- 1}}{T_{\text{FR},1}}} \arctan \left({\sqrt {{{\left({\beta B} \right)}^{- 1}}{T_{\text{FR},1}}}} \right)}}{{2\left({1 - {\beta^{- 1}}} \right)}}\\ + \frac{{{\rho_{2}}\left({1 - {\beta^{- 1}} + \rho_{2}^{\mathrm{e}{\beta^{- 1}}}} \right)\sqrt {{B^{- 1}}{T_{\text{FR},1}}} \arctan \left({\sqrt {{B^{- 1}}{T_{\text{FR},1}}}} \right)}}{{2\left({1 - {\beta^{- 1}}} \right)}} \end{array}.$$

### Theorem 3

(Macro tier, edge user): the coverage probability of the macro edge user is given by (24).

### Proof

The proof is given in Appendix C. □

### Corollary 2

For α=4 and σ 2=0, the coverage probability of the macro edge user is given by (25), where $$\xi \left ({\rho _{1}^{\mathrm {e}},{\rho _{1}},T,{T_{\text {FR},1}},4} \right) =$$

\begin{aligned} &\frac{{\rho_{1}^{\mathrm{e}}\sqrt{T} \arctan \left({\sqrt{T}} \right)\left({{T_{\text{FR},1}} - T - {\rho_{1}}{T_{\text{FR},1}}} \right)}}{{2\left({{T_{\text{FR},1}} - T} \right)}}\\ & + \frac{{{\rho_{1}}\sqrt {{T_{\text{FR},1}}} \arctan \left({\sqrt {{T_{\text{FR},1}}}} \right)\left({{T_{\text{FR},1}} - T + \rho_{1}^{\mathrm{e}}T} \right)}}{{2\left({{T_{\text{FR},1}} - T} \right)}} \end{aligned},

and $$\xi \left ({\rho _{2}^{\mathrm {i}},{\rho _{2}},{{{B}}^{- 1}}T,{B^{- 1}}{T_{\text {FR},1}},4} \right) =$$

{} \begin{aligned} &\frac{{\rho_{2}^{\mathrm{i}}\sqrt {{B^{- 1}}T} \arctan \left({\sqrt {{B^{- 1}}T}} \right)\left({{T_{\text{FR},1}} - T - {\rho_{2}}{T_{\text{FR},1}}} \right)}}{{2\left({{T_{\text{FR},1}} - T} \right)}}\\ & + \frac{{{\rho_{2}}\sqrt {{B^{- 1}}{T_{\text{FR},1}}} \arctan \left({\sqrt {{B^{- 1}}{T_{\text{FR},1}}}} \right)\left({{T_{\text{FR},1}} - T + \rho_{2}^{\mathrm{i}}T} \right)}}{{2\left({{T_{\text{FR},1}} - T} \right)}} \end{aligned}.

When T=T FR,1, the limit TT FR,1 will simplify as $$\xi \left ({\rho _{1}^{\mathrm {e}},{\rho _{1}},T,{T_{\text {FR},1}},4} \right) =$$

$$\frac{{\left({2{\rho_{1}} + 2\rho_{1}^{\mathrm{e}} - {\rho_{1}}\rho_{1}^{\mathrm{e}}} \right)\sqrt {{T_{\text{FR},1}}} \arctan \left({\sqrt {{T_{\text{FR},1}}}} \right)}}{4} + \frac{{{\rho_{1}}\rho_{1}^{\mathrm{e}{T_{\text{FR},1}}}}}{{4\left({{T_{\text{FR},1}} + 1} \right)}},$$

and $$\xi \left ({\rho _{2}^{\mathrm {i}},{\rho _{2}},{{{B}}^{- 1}}T,{B^{- 1}}{T_{\text {FR},1}},4} \right) =$$

$$\frac{{\left({2{\rho_{2}} + 2\rho_{2}^{\mathrm{i}} - {\rho_{2}}\rho_{2}^{\mathrm{i}}} \right)\sqrt {{B^{- 1}}{T_{\text{FR},1}}} \arctan \left({\sqrt {{B^{- 1}}{T_{\text{FR},1}}}} \right)}}{4}+\frac{{{\rho_{2}}\rho_{2}^{\mathrm{i}}{{B^{- 1}}{T_{\text{FR},1}}}}}{{4\left({{B^{- 1}}{T_{\text{FR},1}} + 1} \right)}}.$$
$$p_{\text{FFR},{2}}^{\mathrm{i}}\left(T \right) = \frac{{\lambda_{2}}\left({1 + {\rho_{2}}\sqrt {{T_{\text{FR},2}}} \arctan \left({\sqrt {{T_{\text{FR},2}}}} \right)} \right) + {\lambda_{1}}\sqrt{{P_{1}} / {P_{2}}} \left(\sqrt{1 / B} + {\rho_{1}}\sqrt{T_{\text{FR},2}} \arctan \left({\sqrt {B{T_{\text{FR},2}}}} \right) \right)}{{\lambda_{2}}\left({1 + 2\xi \left({\rho_{2}^{\mathrm{i}},{\rho_{2}},T,{T_{\text{FR},2}},4} \right)} \right) + {\lambda_{1}}\sqrt {{{P_{1}}} / {{P_{2}}B}} \left({1 + 2\xi \left({\rho_{1}^{\mathrm{e}},{\rho_{1}},BT,B{T_{\text{FR},2}},4} \right)} \right)}$$
(27)
$$\begin{array}{l} p_{\text{FFR},2}^{\mathrm{e}}\left(T\right) = \frac{\frac{2\pi {\lambda_{2}}}{{\mathcal A}_{2}}\int_{0}^{\infty} x\exp \left(- \frac{T}{\text{SNR}_{2}\left({x,1}\right)} \right) \exp \left\{- \pi {\lambda_{2}}Q\left({T,\alpha,1,\rho_{2}^{\mathrm{e}}} \right){x^{2}} - \pi {\lambda_{1}}\left(\frac{P_{1}}{P_{2}} \right)^{\frac{2}{\alpha}}Q\left(\beta T,\alpha,{1 / B},\rho_{1}^{\mathrm{i}} \right){x^{2}} \right\}dx}{1 - {p_{c,2}}\left(T_{\text{FR},2}\right)}\\ \qquad \qquad - \frac{\begin{array}{l} \frac{2\pi{\lambda_{2}}}{{\mathcal{A}}_{2}}\int_{0}^{\infty} x\exp \left(- \frac{T + T_{\text{FR},2}}{\text{SNR}_{2}\left({x,1}\right)} \right)\exp \left\{\vphantom{\dot{\left[1 + 2\xi \left(\rho_{1}^{\mathrm{i}},{\rho_{1}},\beta BT,BT_{\text{FR},2},\alpha \right) \right]}\!}\!- \pi {\lambda_{2}}{x^{2}}\left[ 1 + 2\xi \left({\rho_{2}^{\mathrm{e}},{\rho_{2}},T,T_{\text{FR},2},\alpha} \right) \right] \right. \\ \left. \qquad \qquad \qquad \qquad - \pi {\lambda_{1}}\left(\frac{P_{1}}{P_{2}B}\right)^{\frac{2}{\alpha}}{x^{2}}\left[1 + 2\xi \left(\rho_{1}^{\mathrm{i}},{\rho_{1}},\beta BT,BT_{\text{FR},2},\alpha \right) \right] \right\}dx \end{array}}{1 - {p_{c,2}}\left({T_{\text{FR},2}}\right)} \end{array}$$
(28)
$$\begin{array}{l} p_{\text{FFR},2}^{\mathrm{e}}\left(T \right) = \frac{\left({\lambda_{1}}\sqrt{{P_{1}} / {P_{2}B}} + {\lambda_{2}} \right)\left[{\lambda_{2}}\left(1 + {\rho_{2}}\sqrt{T_{\text{FR},2}} \arctan \left(\sqrt{T_{\text{FR},2}} \right)\right) + {\lambda_{1}}\sqrt{{P_{1}} / {P_{2}}} \left(\sqrt{1 / B} + {\rho_{1}}\sqrt{T_{\text{FR},2}} \arctan \left(\sqrt{{BT}_{\text{FR},2}}\right) \right) \right]}{\left(1 - {\lambda_{1}}\sqrt{{P_{1}} / {P_{2}B}} - {\lambda_{2}} \right)\left[\lambda_{2}\left(1 + \rho_{2}^{\mathrm{e}}\sqrt{T} \arctan \left(\sqrt{T}\right)\right) + {\lambda_{1}}\sqrt{{P_{1}} / {P_{2}}} \left(\sqrt{1 / B} + \rho_{1}^{\mathrm{i}}\sqrt{{\beta}T} \arctan \left(\sqrt{{\beta} BT} \right)\right)\right]}\\ \qquad \qquad \qquad - \frac{\left(\lambda_{1}\sqrt{{P_{1}} / {P_{2}}B} + {\lambda_{2}} \right)\left[\lambda_{2}\left(1 + {\rho_{2}}\sqrt{T_{\text{FR},2}} \arctan \left(\sqrt{T_{\text{FR},2}}\right)\right) + {\lambda_{1}}\sqrt{{P_{1}} / {P_{2}}} \left(\sqrt{1 / B} + {\rho_{1}}\sqrt{T_{\text{FR},2}} \arctan \left(\sqrt{{BT}_{\text{FR},2}}\right)\right)\right]}{\left(1 - {\lambda_{1}}\sqrt{{P_{1}} / {P_{2}B}} - {\lambda_{2}} \right)\left[{\lambda_{2}}\left(1 + 2\xi \left({\rho_{2}^{\mathrm{e}},{\rho_{2}},T,{T_{\text{FR},2}},4} \right)\right) + {\lambda_{1}}\sqrt{{P_{1}} / {P_{2}B}} \left(1 + 2\xi\left(\rho_{1}^{\mathrm{i}},{\rho_{1}},\beta BT,{BT}_{\text{FR},2},4 \right)\right)\right]} \end{array}$$
(29)

### Theorem 4

(Pico tier, interior user): the coverage probability of the pico interior user is given by (26).

### Proof

A user $$u \in {\mathcal {U}_{2}}$$ with SINR>T FR,2 is allocated a sub-channel from available sub-channel subset C 2. Conditioning on its previous SINR and applying the Bayes’ rule, we have pFFR,2i(T)

\begin{aligned} &= \mathbb{P}\left({\frac{{{P_{2}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\hat I_{x,1}^{\mathrm{e}} + \hat I_{x,2}^{\mathrm{i}} + {\sigma^{2}}}} > T\left| {\frac{{{P_{2}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} > {T_{\text{FR},2}}} \right.} \right)\\ &= \frac{{\mathbb{P}\left({\frac{{{P_{2}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\hat I_{x,1}^{\mathrm{e}} + \hat I_{x,2}^{\mathrm{i}} + {\sigma^{2}}}} > T, \frac{{{P_{2}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} > {T_{\text{FR},2}}} \right)}}{{\mathbb{P}\left({\frac{{{P_{2}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} > {T_{\text{FR},2}}} \right)}} \end{aligned}

Following the method of Theorem 2 gives the desired result. □

### Corollary 3

For α=4 and σ 2=0, the coverage probability of the pico interior user is given by (27), where $$\xi \left ({\rho _{2}^{\mathrm {i}},{\rho _{2}},T,{T_{\text {FR},2}},4} \right) =$$

\begin{aligned} &\frac{{\rho_{2}^{\mathrm{i}}\sqrt T \arctan \left({\sqrt T} \right)\left({{T_{\text{FR},2}} - T - {\rho_{2}}{T_{\text{FR},2}}} \right)}}{{2\left({{T_{\text{FR},2}} - T} \right)}}\\ & + \frac{{{\rho_{2}}\sqrt {{T_{\text{FR},2}}} \arctan \left({\sqrt {{T_{\text{FR},2}}}} \right)\left({{T_{\text{FR},2}} - T + \rho_{2}^{\mathrm{i}}T} \right)}}{{2\left({{T_{\text{FR},2}} - T} \right)}} \end{aligned},

and $$\xi \left ({\rho _{1}^{\mathrm {e}},{\rho _{1}},BT,B{T_{\text {FR},2}},4} \right) =$$

\begin{aligned} &\frac{{\rho_{1}^{\mathrm{e}}\sqrt {BT} \arctan \left({\sqrt {BT}} \right)\left({{T_{\text{FR},2}} - T - {\rho_{1}}{T_{\text{FR},2}}} \right)}}{{2\left({{T_{\text{FR},2}} - T} \right)}}\\ &+ \frac{{{\rho_{1}}\sqrt {B{T_{\text{FR},2}}} \arctan \left({\sqrt {B{T_{\text{FR},2}}}} \right)\left({{T_{\text{FR},2}} - T + \rho_{1}^{\mathrm{e}}T} \right)}}{{2\left({{T_{\text{FR},2}} - T} \right)}} \end{aligned}.

When T=T FR,2, the limit TT FR,2 will simplify as $$\xi \left ({\rho _{2}^{\mathrm {i}},{\rho _{2}},T,{T_{\text {FR},2}},4} \right) =$$

\begin{aligned} \frac{{\left({2{\rho_{2}} + 2\rho_{2}^{\mathrm{i}} - {\rho_{2}}\rho_{2}^{\mathrm{i}}} \right)\sqrt {{T_{\text{FR},2}}} \arctan \left({\sqrt {{T_{\text{FR},2}}}} \right)}}{4} + \frac{{{\rho_{2}}\rho_{2}^{\mathrm{i}}{T_{\text{FR},2}}}}{{4\left({{T_{\text{FR},2}} + 1} \right)}}, \end{aligned}

and $$\xi \left ({\rho _{1}^{\mathrm {e}},{\rho _{1}},BT,B{T_{\text {FR},2}},4} \right) =$$

\begin{aligned} \frac{{\left({2{\rho_{1}} + 2\rho_{1}^{\mathrm{e}} - {\rho_{1}}\rho_{1}^{\mathrm{e}}} \right)\sqrt {B{T_{\text{FR},2}}} \arctan \left({\sqrt {B{T_{\text{FR},2}}}} \right)}}{4}+\frac{{{\rho_{1}}\rho_{1}^{\mathrm{e}}B{T_{\text{FR},2}}}}{{4\left({B{T_{\text{FR},2}} + 1} \right)}}. \end{aligned}

### Theorem 5

(Pico tier, edge user): the coverage probability of the pico edge user is given by (28).

### Proof

A user $$u \in {\mathcal {U}_{2}}$$ with SINR<T FR,2 is allocated a sub-channel from available sub-channel subset C 1. Conditioning on its previous SINR and applying the Bayes’ rule, we have pFFR,2e(T)

\begin{aligned} &= \mathbb{P}\left({\frac{{{P_{2}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\beta \hat I_{x,1}^{\mathrm{i}} + \hat I_{x,2}^{\mathrm{e}} + {\sigma^{2}}}} > T\left| {\frac{{{P_{2}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} < {T_{\text{FR},2}}} \right.} \right)\\ &= \frac{{\mathbb{P}\left({\frac{{{P_{2}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\beta \hat I_{x,1}^{\mathrm{i}} + \hat I_{x,2}^{\mathrm{e}} + {\sigma^{2}}}} > T,\frac{{{P_{2}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} < {T_{\text{FR},2}}} \right)}}{{\mathbb{P}\left({\frac{{{P_{2}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} < {T_{\text{FR},2}}} \right)}} \end{aligned}

Following the methods of Theorems 1 and 3 give the desired result. □

### Corollary 4

For α=4 and σ 2=0, the coverage probability of the pico edge user is given by (29), where $$\xi \left ({\rho _{2}^{\mathrm {e}},{\rho _{2}},T,{T_{\text {FR},2}},4} \right) =$$

$$\begin{array}{l} \frac{{\rho_{2}^{\mathrm{e}}\sqrt T \arctan \left({\sqrt T} \right)\left({{T_{{\text{FR}},2}} - T - {\rho_{2}}{T_{{\text{FR}},2}}} \right)}}{{2\left({{T_{{\text{FR}},2}} - T} \right)}}\\ + \frac{{{\rho_{2}}\sqrt {{T_{{\text{FR}},2}}} \arctan \left({\sqrt {{T_{{\text{FR}},2}}}} \right)\left({{T_{{\text{FR}},2}} - T + \rho_{2}^{\mathrm{e}}T} \right)}}{{2\left({{T_{{\text{FR}},2}} - T} \right)}} \end{array},$$

and $$\xi \left ({\rho _{1}^{\mathrm {i}},{\rho _{1}},\beta BT,B{T_{{\text {FR}},2}},4} \right) =$$

$$\begin{array}{l} \frac{{\rho_{1}^{\mathrm{i}}\sqrt {\beta BT} \arctan \left({\sqrt {\beta BT}} \right)\left({{T_{{\text{FR}},2}} - \beta T - {\rho_{1}}{T_{{\text{FR}},2}}} \right)}}{{2\left({{T_{{\text{FR}},2}} - \beta T} \right)}}\\ + \frac{{{\rho_{1}}\sqrt {B{T_{{\text{FR}},2}}} \arctan \left({\sqrt {B{T_{{\text{FR}},2}}}} \right)\left({{T_{{\text{FR}},2}} - \beta T + \rho_{1}^{\mathrm{i}}T} \right)}}{{2\left({{T_{{\text{FR}},2}} - \beta T} \right)}} \end{array}.$$

When T=T FR,2, the limit TT FR,2 will simplify as $$\xi \left ({\rho _{2}^{\mathrm {e}},{\rho _{2}},T,{T_{\text {FR},2}},4} \right) =$$

\begin{aligned} \frac{{\left({2{\rho_{2}} + 2\rho_{2}^{\mathrm{e}} - {\rho_{2}}\rho_{2}^{\mathrm{e}}} \right)\sqrt {{T_{\text{FR},2}}} \arctan \left({\sqrt {{T_{\text{FR},2}}}} \right)}}{4} + \frac{{{\rho_{2}}\rho_{2}^{\mathrm{e}}{{T_{\text{FR},2}}}}}{{4\left({{T_{\text{FR},2}} + 1} \right)}}, \end{aligned}

and $$\xi \left ({\rho _{1}^{\mathrm {i}},{\rho _{1}},\beta BT,B{T_{\text {FR},2}},4} \right) =$$

$$\begin{array}{c} \frac{{\rho_{1}^{\mathrm{i}}\left({1 - \beta - {\rho_{1}}} \right)\sqrt {\beta B{T_{\text{FR},2}}} \arctan \left({\sqrt {\beta B{T_{\text{FR},2}}}} \right)}}{{2\left({1 - \beta} \right)}}\\ + \frac{{{\rho_{1}}\left({1 - \beta + \rho_{1}^{\mathrm{i}}} \right)\sqrt {B{T_{\text{FR},2}}} \arctan \left({\sqrt {B{T_{\text{FR},2}}}} \right)}}{{2\left({1 - \beta} \right)}} \end{array}.$$

### Average user rate

Similar to definition of the conditional coverage probabilities, the average rate for a typical user is given as

\begin{aligned} {\mathcal R} = \sum\limits_{k = 1}^{2} {{{\mathcal A}_{k}}\left({{p_{c,k}}\left({{T_{\text{FR},k}}} \right){\mathcal R}_{\text{FFR},k}^{\mathrm{i}} + \left({1 - {p_{c,k}}\left({{T_{\text{FR},k}}} \right)} \right){\mathcal R}_{\text{FFR},k}^{\mathrm{e}}} \right)}. \end{aligned}
(30)

where $${\mathcal R}_{\text {FFR},k}^{\mathrm {i}}$$ and $${\mathcal R}_{\text {FFR},k}^{\mathrm {e}}$$ are the average rate of a typical user when it is an interior user and edge user associated with the kth tier BS, respectively.

The average user rate is computed in units of nats/sec/Hz (1 bit=ln(2)=0.693 nats), where it can represent the spectral efficiency of a user. Conditioning on a typical user at a distance x from its serving BS in the kth tier, $${\mathcal R}_{\text {FFR},k}^{\mathrm {i}}$$ and $${\mathcal R}_{\text {FFR},k}^{\mathrm {e}}$$ are thus defined as

$${\mathcal R}_{\text{FFR},k}^{\mathrm{i}} \buildrel \Delta \over ={\mathbb{E}_{x}}\left[ {{\mathbb{E}_{\text{SINR}}}\left[ {t_{k}^{\mathrm{i}}\ln \left({1 + \text{SINR}} \right)\left| {u \in \mathcal{U}_{k}^{\mathrm{i}}} \right.} \right]} \right],$$
(31)
$${\mathcal R}_{\text{FFR},k}^{\mathrm{e}} \buildrel \Delta \over ={\mathbb{E}_{x}}\left[ {{\mathbb{E}_{\text{SINR}}}\left[ {t_{k}^{\mathrm{e}}\ln \left({1 + \text{SINR}} \right)\left| {u \in \mathcal{U}_{k}^{\mathrm{e}}} \right.} \right]} \right].$$
(32)

where the expectation is taken with respect to the distance x, $$t_{k}^{\mathrm {i}}$$, and $$t_{k}^{\mathrm {e}}$$ are the faction of time the users $$u \in \mathcal {U}_{k}^{\mathrm {i}}$$ and $$u \in \mathcal {U}_{k}^{\mathrm {e}}$$ are served on a sub-channel, respectively. $$t_{k}^{\mathrm {i}}$$ and $$t_{k}^{\mathrm {e}}$$ can be considered as the reciprocal of average number of users sharing a given sub-channel, i.e., the load at a sub-channel.

### Lemma 4

The expectations of the faction of time on a given sub-channel $$\mathbb {E}\left [ {t_{k}^{\mathrm {i}}} \right ]$$ and $$\mathbb {E}\left [ {t_{k}^{\mathrm {e}}} \right ]$$ are given by (3336).

\begin{aligned} &\mathbb{E} \left[ {t_{1}^{\mathrm{i}}} \right] = \frac{{{\lambda_{1}}\left| {\mathbf{C}_{1}} \right|}}{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}\left[ {1 - {{\left({1 + \frac{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}{{3.5{\lambda_{1}}}}} \right)}^{- 3.5}}} \right]\\ &\qquad\quad + \sum\limits_{n = 0}^{\left| {\mathbf{C}_{1}} \right| - 1} {\left({1 - \frac{{\left| {\mathbf{C}_{1}} \right|}}{{n + 1}}} \right) \mathbb{P} \left({\tilde N_{1}^{\mathrm{i}} = n} \right)} \end{aligned}
(33)
{} \begin{aligned} &\mathbb{E} \left[ {t_{1}^{\mathrm{e}}} \right] = \frac{{{\lambda_{1}}\left| {\mathbf{C}}_{2} \right|}}{{{\lambda_{u}}{{\mathcal A}_{1}}\left({1 - {p_{c,1}}\left({{T_{\text{FR},1}}} \right)} \right)}}\left[ {1 - {{\left({1 + \frac{{{\lambda_{u}}{{\mathcal A}_{1}}\left({1 - {p_{c,1}}\left({{T_{\text{FR},1}}} \right)} \right)}}{{3.5{\lambda_{1}}}}} \right)}^{- 3.5}}} \right]\\ &\qquad\quad\!\! + \sum\limits_{n = 0}^{\left| {\mathbf{C}}_{2} \right| - 1} {\left({1 - \frac{{\left| {\mathbf{C}}_{2} \right|}}{{n + 1}}} \right)\mathbb{P}\left({\tilde N_{1}^{\mathrm{e}} = n} \right)} \end{aligned}
(34)
\begin{aligned} &\mathbb{E} \left[ {t_{2}^{\mathrm{i}}} \right] = \frac{{{\lambda_{2}}\left| {\mathbf{C}}_{2} \right|}}{{{\lambda_{u}}{{\mathcal A}_{2}}{p_{c,2}}\left({{T_{\text{FR},2}}} \right)}}\left[ {1 - {{\left({1 + \frac{{{\lambda_{u}}{{\mathcal A}_{2}}{p_{c,2}}\left({{T_{\text{FR},2}}} \right)}}{{3.5{\lambda_{2}}}}} \right)}^{- 3.5}}} \right]\\ &\qquad\quad+ \sum\limits_{n = 0}^{\left| {\mathbf{C}}_{2} \right| - 1} {\left({1 - \frac{{\left| {\mathbf{C}}_{2} \right|}}{{n + 1}}} \right) \mathbb{P} \left({\tilde N_{2}^{\mathrm{i}} = n} \right)} \end{aligned}
(35)
\begin{aligned} &\mathbb{E} \left[ {t_{2}^{\mathrm{e}}} \right] = \frac{{{\lambda_{2}}\left| {\mathbf{C}}_{1} \right|}}{{{\lambda_{u}}{{\mathcal A}_{2}}\left({1 - {p_{c,2}}\left({{T_{\text{FR},2}}} \right)} \right)}}\left[ {1 - {{\left({1 + \frac{{{\lambda_{u}}{{\mathcal A}_{2}}\left({1 - {p_{c,2}}\left({{T_{\text{FR},2}}} \right)} \right)}}{{3.5{\lambda_{2}}}}} \right)}^{- 3.5}}} \right]\\ &\qquad\quad+ \sum\limits_{n = 0}^{\left| {\mathbf{C}}_{1} \right| - 1} {\left({1 - \frac{{\left| {\mathbf{C}}_{1} \right|}}{{n + 1}}} \right) \mathbb{P} \left({\tilde N_{2}^{\mathrm{e}} = n} \right)} \end{aligned}
(36)

### Proof

The proof is given in Appendix D. □

Assuming the independence between the load and SINR [14,15,23], the following theorem gives the average rate for a typical macro interior user. Note that the average rate for other type of users can be obtained following the same procedure.

### Theorem 6

The average rate of a typical macro interior user is

$${\mathcal R}_{\text{FFR},{1}}^{\mathrm{i}} = \mathbb{E} \left[ {t_{1}^{\mathrm{i}}} \right] \times \int_{0}^{\infty} {p_{\text{FFR},1}^{\mathrm{i}}\left({{e^{t}} - 1} \right)} dt.$$
(37)

### Proof

Following the assumption of independence between the load and SINR, we have

$${} {\mathcal R}_{\text{FFR},{1}}^{\mathrm{i}} = \mathbb{E} \left[ {t_{1}^{\mathrm{i}}} \right] \times {\mathbb{E}_{x}}\left[ {{\mathbb{E}_{\text{SINR}}}\left[ {\ln \left({1 + \text{SINR}} \right)\left| {u \in \mathcal{U}_{1}^{\mathrm{i}}} \right.} \right]} \right].$$
(38)

Since $$\mathbb {E}\left [ \tau \right ] = \int _{0}^{\infty } {\mathbb {P}\left ({\tau > t} \right)dt}$$ for τ>0, we obtain □

\begin{aligned} &{\mathbb{E}_{x}}\left[ {{\mathbb{E}_{\text{SINR}}}\left[ {\ln \left({1 + \text{SINR}} \right)\left| {u \in \mathcal{U}_{1}^{\mathrm{i}}} \right.} \right]} \right]\\ &= \int_{0}^{\infty} {\mathbb{P}\left({\ln \left({1 + \text{SINR}} \right) > t\left| {u \in \mathcal{U}_{1}^{\mathrm{i}}} \right.} \right)} dt\\ &= \int_{0}^{\infty} {\mathbb{P}\left({\ln \left({1 + \frac{{\beta {P_{1}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\beta \hat I_{x,1}^{\mathrm{i}} + \hat I_{x,2}^{\mathrm{e}} + {\sigma^{2}}}}} \right) > t\left| {\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} > {T_{\text{FR},1}}} \right.} \right)} dt.\\ &= \int_{0}^{\infty} {\mathbb{P}\left({\frac{{\beta {P_{1}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\beta \hat I_{x,1}^{\mathrm{i}} + \hat I_{x,2}^{\mathrm{e}} + {\sigma^{2}}}} > {e^{t}} - 1\left| {\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} > {T_{\text{FR},1}}} \right.} \right)} dt \end{aligned}

Conditioning on x the distance to its serving macro BS and following the method of Theorem 2, we have

\begin{aligned} {\mathbb{E}_{x}}\left[ {{\mathbb{E}_{\text{SINR}}}\left[ {\ln \left({1 + \text{SINR}} \right)\left| {u \in \mathcal{U}_{1}^{\mathrm{i}}} \right.} \right]} \right] = \int_{0}^{\infty} {p_{\text{FFR},{1}}^{\mathrm{i}}\left({{e^{t}} - 1} \right)} dt. \end{aligned}

Plugging back into (38), the average rate (37) can be obtained.

## Simulation and numerical results

In this section, we present numerical results on the coverage probability and average user rate for the proposed SFR scheme. Furthermore, we compare the proposed SFR scheme with FRP scheme. The simulation parameters are in accordance with 3GPP technical reports . Unless otherwise stated, the transmit powers of a macro BS and a pico BS are $$P_{1}^{T} = 46 \ {\text {dBm}}$$ and $$P_{2}^{T} = 30 \ {\text {dBm}}$$. The densities of the two tiers are λ 1=1 BS/km2 and λ 2=5 BS/km2 with α = 4. The user density is λ u = 100 users/km2. For a typical LTE system with 10 MHz bandwidth, 50 sub-channels (i.e., |C| = 50) are available to each BS, each sub-channel bandwidth then is 200 kHz.

### Validation of analysis

In Fig. 2, the average user rate sweeping over a range of user density λ u is validated via Monte Carlo methods on a square window of 10 × 10 km2. There are three values for the combination of association bias, power control factor and resource partitioning factor (B,β,η) are shown. It can be obtained that the analytical results match quite well with the simulation results. Meanwhile, the simulation results give upper bound with less than a 0.26 dB gap from the analytical results. It also can be observed that the gap decreases with the increase in user density (namely, network load). In addition to the approximation for Voronoi cell areas distribution, the small gap is also due to the assumption of the independence between the load and SINR.

### Comparison with FRP scheme

Figure 3 compares the average user rate of the proposed SFR scheme and FRP scheme proposed by  in different load conditions. In this paper, the cell load model is dependent on user density and parameters of SFR scheme, apart from association bias. The average user rate of both schemes decreases with the increase in user density. This is because the number of users shared a common resource increases and so does the interference. From this figure, it can be seen that the proposed SFR scheme outperforms the FRP scheme in terms of average user rate in any load condition. The curves indicate that the proposed SFR scheme can enhance the average user rate by about 7.49% when λ u = 100 users/km2. Moreover, the average user rate can be further improved by properly setting up the parameters of SFR scheme.

### Average user rate: trends and discussion

The conditional coverage probabilities without SFR are significant for the classification of users in each tier. Figure 4 shows the variations in the previous cell load ρ k and p c,k (T FR,k ) with increasing association bias. It can be seen that the previous cell load of macro tier ρ 1 decreases with increasing association bias. So the coverage probability without SFR of macro users p c,1(T FR,1) increases because the interference from macro tier decreases. On the other hand, despite that the previous cell load of pico tier ρ 2 is directly proportional to the association bias, the coverage probability without SFR of pico users p c,2(T FR,2) firstly decreases and then increases with increasing association bias. This is because the offloaded users usually do not get the strongest received power and it results in degraded SINR, so the coverage performance of pico tier initially deteriorates. When beyond a certain association bias (10 dB in Fig. 4 b), the interference from macro tier for pico users becomes fairly low and also continues to decrease. Then, the aggregate interference decreases and it will results in improved coverage performance.

Figure 5 a shows that the conditional probabilities that a user is macro interior user and edge user (macro tier association probabilities) decrease with increasing association bias. Meanwhile, the conditional probability that a user is pico interior user increases with increasing association bias. However, the conditional probability that a user is pico edge user firstly increases and then decreases with increasing association bias. That is, the majority of offloaded users initially become pico edge users owing to their worse SINRs. When beyond a certain association bias, the interference from macro tier becomes fairly low, and thereby decreases its pico edge association probability. Figure 5 b shows the variations in the four types of cell loads with SFR with increasing association bias. It can be observed that the variations in the four types of cell loads with SFR are the same as that in Fig. 5 a. As η increases, more sub-channels become available for macro interior users and pico edge users, consequently their cell loads $$\rho _{1}^{\mathrm {i}}$$ and $$\rho _{2}^{\mathrm {e}}$$ decreases which are different from the cell loads $$\rho _{1}^{\mathrm {e}}$$ and $$\rho _{2}^{\mathrm {i}}$$.

For the sake that the proposed SFR scheme can be better understood, we analyze the conditional coverage probabilities and average rate of four types of users (i.e., the interior and edge users of macro and pico BSs). In Fig. 6, it can be seen that the coverage probability and average rate of macro users (both interior and edge users) increase with increasing association bias. This is because more macro users are offloaded to pico tier and more sub-channels become available for macro users, namely, the decrease in macro loads $$\rho _{1}^{\mathrm {i}}$$ and $$\rho _{1}^{\mathrm {e}}$$. Similar to , the variation in the performance with bias of users is mainly due to the change in the cell loads (i.e., $$\rho _{1}^{\mathrm {i}}$$, $$\rho _{1}^{\mathrm {e}}$$, $$\rho _{2}^{\mathrm {i}}$$, and $$\rho _{2}^{\mathrm {e}}$$), which are directly related to the interference from each tier. Moreover, there are obviously differences from previous works in terms of the coverage and rate trends of pico edge users, which their performance firstly decreases and then increases with increasing association bias. It is due to that the majority of offloaded users initially become pico edge users owing to their worse SINRs. Hence, the performance of pico edge users initially deteriorates, despite that the interference from macro tier decreases. When beyond a certain association bias, the interference from macro tier becomes fairly low, and thereby the performance of pico edge users increases. But, the performance of pico interior users has no significant change with respect to the increasing association bias. To mitigate the severe interference from macro tier suffered by pico edge users, especially the offloaded macro users, power control (transmit power reduction) is used for macro BSs. It can be obtained in Fig. 6 that the coverage probability and average rate of pico edge users are improved at the expense of the performance degradation of macro interior users. As η increases, more sub-channels become available for pico edge users, consequently their performance can be further improved. However, the performance of macro edge users and pico interior users remains relatively unchanged with regard to the power control.

#### Impact of SFR thresholds

For the proposed SFR scheme, SFR thresholds T FR,1 and T FR,2 are important design parameters and play a key role in controlling the gains of pico edge users. The effect of SFR threshold pair (T FR,1,T FR,2) on average user rate is shown in Fig. 7. It can be observed that the importance of properly dividing interior and edge regions on the average user rate. The optimal pair is (− 20 dB,4 dB) so that the average user rate attains a maximal value of 2.4304 nats/sec/Hz, which gives a significant increase (by around 47.7%) compared to the worst case that just gets 1.6454 nats/sec/Hz with (20 dB,− 20 dB). With the macro’s SFR threshold T FR,1 decreases gradually from 20 dB to − 20 dB and the matched pico’s SFR threshold T FR,2 increases from − 20 dB to 20 dB, the average user rate initially increases, but decreases beyond a certain SFR threshold pair and hence, the optimal SFR threshold pair exists. For the given value of association bias (B = 10 dB), the total number of pico users is much greater than the total number of macro users. Therefore, the variation in the overall average user rate is mainly contributed by pico users. As T FR,2 increases, the number of pico edge users increases while the number of pico interior users decreases. On the contrary, the number of macro interior users increase and the number of macro edge users decrease with the decrease in T FR,1. With the variation in SFR threshold pair, the average rate of a typical pico interior user decreases because the increase in shared resources and the decrease in interference from macro tier. Meanwhile, the average rate of a typical pico edge user varies inversely. Thus, the contribution of pico interior users towards the average rate initially dominates. But after a certain SFR threshold pair, the contribution of pico edge users towards the average rate eventually dominates. Furthermore, the main interference resulting in the degradation of network performance is the strong interference from macro tier. Hence, the macro interior regions tend to be larger because of the introduced power control can mitigate such strong interference.

#### Impact of power control

Figure 8 shows the effect of association bias and power control factor on average user rate. It can be seen that the average user rate can be significantly improved by introducing power control for macro interior users. With increasing power control factor, the average user rate initially increases, but decreases beyond a certain power control factor and hence, the optimal power control factor exists. Note that the average rate of macro edge users and pico interior users are invariant to power control factor. Meanwhile, with decreasing power control factor, the increase in average rate of pico edge users dominates in comparison with the decrease in average rate of macro interior users. Thus, the overall average rate is improved. But after a certain power control factor, the average rate of macro interior users is so small that the overall average rate decreases. It can be obviously observed that the optimal power control factor tends to be smaller because of the large transmit power of macro BS.

For a given power control factor, the average user rate initially increases as the association bias increases, but decreases beyond a certain association bias. Hence, the optimal association bias exists. This is because more macro users in poor situation is offloaded to pico tier with the increase in association bias and thus, the interference from macro tier decreases due to the macro cell load decreases. This leads to the average user rate initially increases. But after a certain association bias, the pico BS is overloaded so that the average user rate eventually decreases.

#### Impact of resource partitioning

Figure 9 shows the effect of association bias and resource partitioning factor on average user rate. For a given resource partitioning factor, the average user rate initially increases with increasing association bias, but decreases beyond a certain association bias. The reason of variation trend was already discussed above. For a given association bias, the average user rate initially increases with increasing resource partitioning factor, but decreases beyond a certain resource partitioning factor. This is because, more sub-channels become available for macro interior users and pico edge users as resource partitioning factor increases, then the corresponding cell load and interference decreases. Consequently, their average user rate increases. The average rate of macro edge users and pico interior users decreases because less sub-channels are allocated to them, resulting in their cell load and interference increases. Thus, the contribution of macro interior users and pico edge users towards the average rate initially dominates. But after a certain resource partitioning factor, the contribution of macro edge users and pico interior users towards the average rate eventually dominates. Furthermore, the resource partitioning factor is tightly related to association bias and SFR thresholds because of our practical cell load model.

#### Optimal average user rate

As discussed above, the association bias and parameters of the SFR scheme need to be carefully chosen for optimal average user rate. The average user rate is calculated from the combination of resource partitioning factor and power control factor. However, the results are obtained from varied cell load for each pair. That is, the association bias increases from 0 to 40 dB, and the threshold T FR,1 decreases gradually from 20 to − 20 dB and the matched threshold T FR,2 increases from − 20 to 20 dB. Hence, the optimal average user rate is 2.6217 nats/sec/Hz when B = 18 dB, η = 0.72, β = 0.03, T FR,1 = − 20 dB, and T FR,2 = 12 dB. Compared to the network without CRE that just gets 1.7704 nats/sec/Hz, the gain obtained from the CRE with SFR scheme can be as high as 48.1%.

## Conclusions

In this paper, we have presented an analytical framework to evaluate average user rate of HCNs with CRE and SFR scheme in a multichannel environment using the tools of stochastic geometry. Taking the random resource allocation and cell load into consideration, our network model can adequately capture the impact of CRE and SFR scheme on the performance. The numerical results have shown that the average user rate can be improved significantly, and further, the proposed SFR scheme outperforms the FRP scheme in any load condition. In addition, the optimal average user rate can be achieved by properly tuning the association bias and parameters of the SFR scheme. With the optimal combination, the gain can be as high as 48.1%.

## Appendix A

### Proof of Theorem 1

The conditional coverage probability without SFR of user $${u \in {\mathcal {U}_{k}}}$$ is

$${} {p_{c,k}}\left(T \right) = \int_{x = 0}^{\infty} {\mathbb{P}\left({\text{SINR} > T\left| {u \in {\mathcal{U}_{k}},{X_{k}} = x} \right.} \right)} {f_{{X_{k}}}}\left(x \right)dx.$$
(39)

Using the SINR expression in (6), we have

\begin{aligned} & \mathbb{P}\left({\text{SINR} > T\left| {u \in {\mathcal{U}_{1}},{X_{1}} = x} \right.} \right) = \mathbb{P}\left({\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}} + {\sigma^{2}}} }} > T} \right)\\ & = \mathbb{P}\left({{H_{x}} > {x^{\alpha} }P_{1}^{- 1}T\left\{ {{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}} + {\sigma^{2}}}} \right\}} \right)\\ & \overset{(a)}{=} {\mathbb{E}_{{I_{x,j}}}}\left[ {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}T\left\{ {{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}} + {\sigma^{2}}}} \right\}} \right)} \right]\\ & \overset{(b)}{=} \exp \left({ - \frac{T}{{\text{SNR}_{1}}(x)}} \right)\prod\limits_{j = 1}^{2} {{\mathbb{E}_{{I_{x,j}}}}\left[ {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}T{I_{x,j}}} \right)} \right]} \\ & = \exp \left({ - \frac{T}{{\text{SNR}_{1}}\left(x \right)}} \right)\prod\limits_{j = 1}^{2} {{\mathcal{L}_{{I_{x,j}}}}\left({{x^{\alpha} }P_{1}^{- 1}T} \right)} \end{aligned}
(40)

where $${\text {SNR}_{1}}\left (x \right) = \frac {{{P_{1}}{x^{- \alpha }}}}{{{\sigma ^{2}}}}$$, (a) follows that H x exp(1), (b) follows from the independence of I x,j and $${\mathcal {L}_{{I_{x,j}}}}\left (s \right)$$ is the Laplace transform of interference. Similarly

$$\begin{array}{l} \mathbb{P}\left({\text{SINR} > T\left| {u \in {\mathcal{U}_{2}},{X_{2}} = x} \right.} \right)\\ \qquad = \exp \left({ - \frac{T}{{{\text{SNR}}_{2}\left(x \right)}}} \right)\prod\limits_{j = 1}^{2} {{\mathcal{L}_{{I_{x,j}}}}\left({{x^{\alpha} }P_{2}^{- 1}T} \right)}. \end{array}$$
(41)

Considering the definition of the Laplace transform and the serving BS being b k , we have

{} \begin{aligned} {\mathcal{L}_{{I_{x,j}}}}\left(s \right) &= {\mathbb{E}_{{I_{x,j}}}}\left[ {\exp \left({ - s{I_{x,j}}} \right)} \right]\\ & = {\mathbb{E}_{{\Psi_{j}},\left\{ {{H_{y}}} \right\}}}\left[ {\exp \left({ - s{P_{j}}\sum\limits_{y \in {\Psi_{j}}\backslash {b_{k}}} {{H_{y}}{{\left\| y \right\|}^{- \alpha }}}} \right)} \right]\\ & \overset{(a)}{=} {\mathbb{E}_{{\Psi_{j}}}}\left[ {\prod\limits_{y \in {\Psi_{j}}\backslash {b_{k}}} {{\mathbb{E}_{{H_{y}}}}} \left[ {\exp \left({ - s{P_{j}}{H_{y}}{{\left\| y \right\|}^{- \alpha }}} \right)} \right]} \right]\\ & = {\mathbb{E}_{{\Psi_{j}}}}\left[ {\prod\limits_{y \in {\Psi_{j}}\backslash {b_{k}}} {{\mathcal{L}_{{H_{y}}}}\left({s{P_{j}}{{\left\| y \right\|}^{- \alpha }}} \right)}} \right]\\ & \overset{(b)}{=} \exp \left({ - 2\pi {\rho_{j}}{\lambda_{j}}\int_{{d_{jk}}\left(x \right)}^{\infty} {\left({1 - {\mathcal{L}_{{H_{y}}}}\left({s{P_{j}}{v^{- \alpha }}} \right)} \right)} vdv} \right)\\ & \overset{(c)}{=} \exp \left({ - 2\pi {\rho_{j}}{\lambda_{j}}\int_{{d_{jk}}\left(x \right)}^{\infty} {\frac{v}{{1 + {{\left({s{P_{j}}} \right)}^{- 1}}{v^{\alpha} }}}dv}} \right) \end{aligned}

where (a) is obtained from the independence of H y , (b) follows from the probability generating functional (PGFL)  of Ψ j and replacing v=y, and (c) follows that H y exp(1). Moreover, d jk (x) is the lower bound on distance of the jth tier when the typical user $$u \in {\mathcal {U}_{k}}$$, which can be obtained by using (1) as

\begin{aligned} \text{if} \ k = 1: {d_{11}} = x,{d_{21}} = \left({P_{2}B} / {P_{1}} \right)^{1 / \alpha }x\\ \text{if} \ k = 2: {d_{12}} = \left({P_{1}} / {P_{2}B} \right)^{1 / \alpha }x,{d_{22}} = x \end{aligned}.

Using the change of variables (sP j )−2/α v 2u, the integral can be simplified as

\begin{aligned} &\int_{{d_{jk}}\left(x \right)}^{\infty} {\frac{v}{{1 + {{\left({s{P_{j}}} \right)}^{- 1}}{v^{\alpha} }}}dv} \\ & \qquad \qquad = \frac{1}{2}\left({sP}_{j} \right)^{2 / \alpha }\int_{\left({sP}_{j}\right)^{- 2 / \alpha}{\left({{d_{jk}}\left(x \right)} \right)}^{2}}^{\infty} {\frac{du}{1 + {u^{\alpha / 2}}}} \\ & \qquad \qquad = \frac{1}{2}\left({sP}_{j}\right)^{2 / \alpha}Z\left(1,\alpha,\frac{\left(d_{jk}\left(x\right)\right)^{\alpha}}{s{P_{j}}} \right) \end{aligned}.

Hence, the Laplace transform of interference is

$$\begin{array}{@{}rcl@{}}{} {\mathcal{L}_{I_{x,j}}}(s) = \exp \left(\! - \pi {\rho_{j}}{\lambda_{j}}\left({sP}_{j}\right)^{2 / \alpha}{Z\left(\!1,\alpha,\frac{\left({d_{jk}}(x)\right)^{\alpha}}{{sP}_{j}} \!\right)} \!\right), \end{array}$$
(42)

where $$Z\left ({a,b,c} \right) = {a}^{2 / b}{\int \nolimits }_{\left ({\frac {c}{a}} \right)^{{2 / b}}}^{\infty }{\frac {du}{1 + {u^{b / 2}}}}$$.

Using (3) in (39) along with (4042), the conditional coverage probabilities given in Theorem 1 are obtained.

## Appendix B

### Proof of Theorem 2

A user $$u \in {\mathcal {U}_{1}}$$ with SINR>T FR,1 is allocated a sub-channel from available sub-channel subset C 1 and will experience new fading power $${\hat H_{x}}$$ and out-of-cell interference from both tiers, instead of H x and $${\sum \nolimits }_{j = 1}^{2} {{I_{x,j}}}$$. Conditioning on its previous SINR, the coverage probability of the macro interior user is

\begin{aligned} & p_{\text{FFR},1}^{\mathrm{i}}\left(T \right) = \\ & \qquad \mathbb{P}\left({\frac{{\beta {P_{1}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\beta \hat I_{x,1}^{\mathrm{i}} + \hat I_{x,2}^{\mathrm{e}} + {\sigma^{2}}}} > T\left| {\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} > {T_{\text{FR},1}}} \right.} \right) \end{aligned}
(43)

Following the Bayes’ rule, we have

\begin{aligned} & p_{\text{FFR},1}^{\mathrm{i}}\left(T \right) = \\ & \qquad \frac{{\mathbb{P}\left({\frac{{\beta {P_{1}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\beta \hat I_{x,1}^{\mathrm{i}} + \hat I_{x,2}^{\mathrm{e}} + {\sigma^{2}}}} > T,\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} > {T_{\text{FR},1}}} \right)}}{{\mathbb{P}\left({\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}}} > {T_{\text{FR},1}}} \right)}} \end{aligned}
(44)
$$\overset{(a)}{=} \frac{\begin{array}{l} \mathbb{E}\left[ {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}T\left({\hat I_{x,1}^{\mathrm{i}} + {\beta^{- 1}}\hat I_{x,2}^{\mathrm{e}} + {\beta^{- 1}}{\sigma^{2}}} \right)} \right) \times} \right.\\ \qquad \quad \left. {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}{T_{\text{FR},1}}\left({{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}} \right)} \right)} \right] \end{array}}{{\mathbb{E}\left[ {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}{T_{\text{FR},1}}\left({{I_{x,1}} + {I_{x,2}} + {\sigma^{2}}} \right)} \right)} \right]}}$$

where (a) follows from the independence of $${\hat H_{x}}$$ and H x , and $${\hat H_{x}},{H_{x}} \sim \exp \left (1 \right)$$.

Conditioning on x, the distance to its serving macro BS, and focusing on the numerator of (44), we observe that the noise is independence of the interference and the expectation with respect to $$\hat I_{x,1}^{\mathrm {i}}$$, I x,1, $$\hat I_{x,2}^{\mathrm {e}}$$ and I x,2 is the joint Laplace transform $$\mathcal {L}\left ({{{\hat s}_{1}},{s_{1}},{{\hat s}_{2}},{s_{2}}} \right)$$ of $$\hat I_{x,1}^{\mathrm {i}}$$, I x,1, $$\hat I_{x,2}^{\mathrm {e}}$$, and I x,2 evaluated at $$\left ({{x^{\alpha } }P_{1}^{- 1}T,{x^{\alpha } }P_{1}^{- 1}{T_{\text {FR},1}}} \right.,{x^{\alpha } }{\left ({\beta {P_{1}}} \right)^{- 1}}T,\left. {{x^{\alpha } }P_{1}^{- 1}{T_{\text {FR},1}}} \right)$$.

Then, the joint Laplace transform is $$\mathcal {L}\left ({{{\hat s}_{1}},{s_{1}},{{\hat s}_{2}},{s_{2}}} \right)$$

\begin{aligned} & = \mathbb{E}\left[ {\exp \left({ - {{\hat s}_{1}}\hat I_{x,1}^{\mathrm{i}} - {s_{1}}{I_{x,1}} - {{\hat s}_{2}}\hat I_{x,2}^{\mathrm{e}} - {s_{2}}{I_{x,2}}} \right)} \right]\\ & \overset{(a)}{=} \mathbb{E}\left[ {\exp \left({ - {{\hat s}_{1}}\hat I_{x,1}^{\mathrm{i}} - {s_{1}}{I_{x,1}}} \right)} \right]\mathbb{E}\left[ {\exp \left({ - {{\hat s}_{2}}\hat I_{x,2}^{\mathrm{e}} - {s_{2}}{I_{x,2}}} \right)} \right], \end{aligned}
(45)

where (a) follows from the independence of $$\hat I_{x,1}^{\mathrm {i}} + {I_{x,1}}$$ and $$\hat I_{x,2}^{\mathrm {e}} + {I_{x,2}}$$.

Factoring out the first term of (45), we have the derivation of (46), where (a) follows that $${\hat H_{x}},{H_{x}} \sim \exp \left (1 \right)$$ and (b) is obtained from the probability generating functional (PGFL)  of Φ 1, replacing v=y and x is the lower bound on distance between the macro user and its own tier.

$$\begin{array}{@{}rcl@{}} && \mathbb{E}\left[ {\exp \left({ - {{\hat s}_{1}}\hat I_{x,1}^{\mathrm{i}} - {s_{1}}{I_{x,1}}} \right)} \right]\\ && = \mathbb{E}\left[ {\exp \left({ - {{\hat s}_{1}}{P_{1}}\sum\limits_{y \in \Xi_{1}^{\mathrm{i}}} {{{\hat H}_{y}}{{\left\| y \right\|}^{- \alpha }} - {s_{1}}{P_{1}}\sum\limits_{y \in {\Psi_{1}}} {{H_{y}}{{\left\| y \right\|}^{- \alpha }}}} } \right)} \right] \\ && = \mathbb{E}\prod\limits_{y \in {\Phi_{1}}} {\exp \left({ - {{\hat s}_{1}}{P_{1}}{\mathbf{1}}\left({y \in \Xi_{1}^{\mathrm{i}}} \right){{\hat H}_{y}}{{\left\| y \right\|}^{- \alpha }} - {s_{1}}{P_{1}}{\mathbf{1}}\left({y \in {\Psi_{1}}} \right){H_{y}}{{\left\| y \right\|}^{- \alpha }}} \right)} \\ && = \mathbb{E}\prod\limits_{y \in {\Phi_{1}}} {\left({1 - \rho_{1}^{\mathrm{i}}\left({1 - \exp \left({ - {{\hat s}_{1}}{P_{1}}{{\hat H}_{y}}{{\left\| y \right\|}^{- \alpha }}} \right)} \right)} \right)} \left({1 - {\rho_{1}}\left({1 - \exp \left({ - {s_{1}}{P_{1}}{H_{y}}{{\left\| y \right\|}^{- \alpha }}} \right)} \right)} \right) \end{array}$$
(46)
$$\begin{array}{@{}rcl@{}} && \overset{(a)}{=} \mathbb{E}\prod\limits_{y \in {\Phi_{1}}} {\left({1 - \rho_{1}^{\mathrm{i}}\left({1 - \frac{1}{{1 + {{\hat s}_{1}}{P_{1}}{{\left\| y \right\|}^{- \alpha }}}}} \right)} \right)} \left({1 - {\rho_{1}}\left({1 - \frac{1}{{1 + {s_{1}}{P_{1}}{{\left\| y \right\|}^{- \alpha }}}}} \right)} \right) \\ && \overset{(b)}{=} \exp \left(- 2\pi {\lambda_{1}}\int_{x}^{\infty} \left[ 1 - \left({1 - \rho_{1}^{\mathrm{i}}\left({1 - \frac{1}{{1 + {{\hat s}_{1}}{P_{1}}{v^{- \alpha }}}}} \right)} \right) \left({1 - {\rho_{1}}\left({1 - \frac{1}{{1 + {s_{1}}{P_{1}}{v^{- \alpha }}}}} \right)} \right) \right] \text{vdv} \right) \\ && \qquad \quad \exp \left(- 2\pi {\lambda_{1}}{x^{2}}\int_{1}^{\infty} \left[ 1 - \left({1 - \rho_{1}^{\mathrm{i}}\left({1 - \frac{1}{{1 + T{v^{- \alpha }}}}} \right)} \right) \left({1 - {\rho_{1}}\left({1 - \frac{1}{{1 + {T_{\text{FR},1}}{v^{- \alpha }}}}} \right)} \right) \right] \text{vdv} \right) \times \\ && \exp \left(\! { - 2\pi {\lambda_{2}}{{\left(\! {\frac{{{P_{2}}B}}{{{P_{1}}}}} \!\right)}^{\frac{2}{\alpha }}}{x^{2}}\int_{1}^{\infty} {\left[\! {1 - \left(\! {1 - \rho_{2}^{\mathrm{e}}\left(\! {1 - \frac{1}{{1 + {{\left({\beta B} \right)}^{- 1}}T{v^{- \alpha }}}}} \!\right)}\! \right)\left(\! {1 - {\rho_{2}}\left(\! {1 - \frac{1}{{1 + {B^{- 1}}{T_{\text{FR},1}}{v^{- \alpha }}}}} \!\right)} \!\right)} \!\right] \text{vdv}}} \!\right) \end{array}$$
(47)
$$\begin{array}{@{}rcl@{}} && \frac{{2\pi {\lambda_{1}}}}{{{{\mathcal A}_{1}}}}\int_{0}^{\infty} x\exp \left(- \frac{T}{{\text{SNR}_{1}}\left({x,\beta} \right)} - \frac{T_{\text{FR},1}}{\text{SNR}_{1}\left({x,1} \right)} \right) \\ && \exp \left\{ - \pi {\lambda_{1}}{x^{2}}\left[ {1 + 2\xi \left({\rho_{1}^{\mathrm{i}},{\rho_{1}},T,{T_{\text{FR},1}},\alpha} \right)} \right] - \pi {\lambda_{2}}{{\left({\frac{{{P_{2}}B}}{{{P_{1}}}}} \right)}^{\frac{2}{\alpha }}}{x^{2}}\left[ {1 + 2\xi \left({\rho_{2}^{\mathrm{e}},{\rho_{2}},{{\left({\beta B} \right)}^{- 1}}T,{B^{- 1}}{T_{\text{FR},1}},\alpha} \right)} \right] \right\}dx \end{array}$$
(48)

Similarly, $$\mathbb {E}\left [ {\exp \left ({ - {{\hat s}_{2}}\hat I_{x,2}^{\mathrm {e}} - {s_{2}}{I_{x,2}}} \right)} \right ] =$$

\begin{aligned} \exp \left({ - 2\pi {\lambda_{2}}\int_{{{\left({\frac{{{P_{2}}B}}{{{P_{1}}}}} \right)}^{\frac{1}{\alpha }}}x}^{\infty} {\left[ {1 - \left({1 - \rho_{2}^{\mathrm{e}}\left({1 - \frac{1}{{1 + {{\hat s}_{2}}{P_{2}}{v^{- \alpha }}}}} \right)} \right)} \right.}} \right.\\ \qquad \qquad \qquad \qquad\left. {\left. { \times \left({1 - {\rho_{2}}\left({1 - \frac{1}{{1 + {s_{2}}{P_{2}}{v^{- \alpha }}}}} \right)} \right)} \right]{\text{vdv}}} \right). \end{aligned}

Hence, $$\mathcal {L}\left ({{x^{\alpha } }P_{1}^{- 1}T,{x^{\alpha } }P_{1}^{- 1}{T_{\text {FR},1}},{x^{\alpha } }{{\left ({\beta {P_{1}}} \right)}^{- 1}}T,{x^{\alpha } }P_{1}^{- 1}{T_{\text {FR},1}}} \right)$$ is given by (47).

Deconditioning on x, the numerator of (44) can is given by (48), where ξ(a,b,c,d,e)=

\begin{aligned} \int_{1}^{\infty} {\left[ {1 - \left({1 - a\left({1 - \frac{1}{{1 + {{\text{cv}}^{- e}}}}} \right)} \right)\left({1 - b\left({1 - \frac{1}{{1 + {{\text{dv}}^{- e}}}}} \right)} \right)} \right]{\text{vdv}}}. \end{aligned}

The denominator of (44) can be obtained from (9), we have the conditional coverage probability given in Theorem 2.

## Appendix C

### Proof of Theorem 3

A user $$u \in {\mathcal {U}_{1}}$$ with SINR < T FR,1 is allocated a sub-channel from available sub-channel subset C 2 and will experience new fading power $${\hat H_{x}}$$ and out-of-cell interference $$I_{x,1}^{\mathrm {e}}~+~I_{x,2}^{i}$$, instead of H x and $${\sum \nolimits }_{j = 1}^{2} {{I_{x,j}}}$$. The interference $$I_{x,1}^{\mathrm {e}}~+~I_{x,2}^{i}$$ is different from H x and $${\sum \nolimits }_{j = 1}^{2} {{I_{x,j}}}$$ because of the disparity of each BS’ load. Then, conditioning on its previous SINR, the coverage probability of the macro edge user is

$$\begin{array}{l} p_{\text{FFR},1}^{\mathrm{e}}\left(T \right) = \\ \qquad \mathbb{P}\left({\frac{{{P_{1}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\hat I_{x,1}^{\mathrm{e}} + \hat I_{x,2}^{\mathrm{i}} + {\sigma^{2}}}} > T\left| {\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}}} + {\sigma^{2}}}} < {T_{\text{FR},1}}} \right.} \right). \end{array}$$
(49)

Applying the Bayes’ rule, we have

\begin{aligned} &p_{\text{FFR},1}^{\mathrm{e}}\left(T \right)\\ &= \frac{{\mathbb{P}\left({\frac{{{P_{1}}{{\hat H}_{x}}{x^{- \alpha }}}}{{\hat I_{x,1}^{\mathrm{e}} + \hat I_{x,2}^{\mathrm{i}} + {\sigma^{2}}}} > T,\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}}} + {\sigma^{2}}}} < {T_{FR,1}}} \right)}}{{\mathbb{P}\left({\frac{{{P_{1}}{H_{x}}{x^{- \alpha }}}}{{{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}}} + {\sigma^{2}}}} < {T_{\text{FR},1}}} \right)}}\\ &\overset{(a)}{=} \frac{\begin{array}{l} \mathbb{E}\left[ {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}T\left({\hat I_{x,1}^{\mathrm{e}} + \hat I_{x,2}^{\mathrm{i}} + {\sigma^{2}}} \right)} \right) \times} \right.\\ \left. {\left({1 - \exp \left({ - {x^{\alpha} }P_{1}^{- 1}{T_{\text{FR},1}}\left({{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}}} + {\sigma^{2}}} \right)} \right)} \right)} \right] \end{array}}{{1 - \mathbb{E}\left[ {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}{T_{\text{FR},1}}\left({{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}}} + {\sigma^{2}}} \right)} \right)} \right]}} \end{aligned}
(50)

where (a) follows from the independence of $${\hat H_{x}}$$ and H x , and $${\hat H_{x}},{H_{x}} \sim \exp \left (1 \right)$$. The numerator of (50) can be decomposed as

\begin{aligned} &\mathbb{E}\left[ {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}T\left({\hat I_{x,1}^{\mathrm{e}} + \hat I_{x,2}^{\mathrm{i}} + {\sigma^{2}}} \right)} \right)} \right] - \\ & \qquad \mathbb{E}\left[ {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}T\left({\hat I_{x,1}^{\mathrm{e}} + \hat I_{x,2}^{\mathrm{i}} + {\sigma^{2}}} \right)} \right)} \right. \times \\ & \qquad \qquad \left. {\exp \left({ - {x^{\alpha} }P_{1}^{- 1}{T_{\text{FR},1}}\left({{\sum\nolimits}_{j = 1}^{2} {{I_{x,j}}} + {\sigma^{2}}} \right)} \right)} \right] \end{aligned}

The first term of the numerator represents the SINR on the newly allocated sub-channel and following the method of Theorem 1 gives (51).

$$\begin{array}{@{}rcl@{}} &&\frac{2\pi {\lambda_{1}}}{{\mathcal A}_{1}}\int_{0}^{\infty} {x\exp \left({ - \frac{T}{{\text{SNR}_{1}}\left(x \right)}} \right)} \exp \left\{ - \pi {\lambda_{1}}Q\left({T,\alpha,1,\rho_{1}^{\mathrm{e}}} \right){x^{2}} - \pi {\lambda_{2}}{\left({P_{2}} / {P_{1}} \right)}^{2 / \alpha}Q\left({T,\alpha,B,\rho_{2}^{\mathrm{i}}} \right){x^{2}} \right\}dx \end{array}$$
(51)
$$\begin{array}{@{}rcl@{}} &&\qquad \exp \left({ - 2\pi {\lambda_{1}}{x^{2}}\int_{1}^{\infty} {\left[ {1 - \left({1 - \rho_{1}^{\mathrm{e}}\left({1 - \frac{1}{{1 + T{v^{- \alpha }}}}} \right)} \right)\left({1 - {\rho_{1}}\left({1 - \frac{1}{{1 + {T_{\text{FR},1}}{v^{- \alpha }}}}} \right)} \right)} \right]{\text{vdv}}}} \right) \times \\ &&\exp \left({ - 2\pi {\lambda_{2}}{{\left({\frac{{{P_{2}}B}}{{{P_{1}}}}} \right)}^{\frac{2}{\alpha }}}{x^{2}}\int_{1}^{\infty} {\left[ {1 - \left({1 - \rho_{2}^{\mathrm{i}}\left({1 - \frac{1}{{1 + {B^{- 1}}T{v^{- \alpha }}}}} \right)} \right)\left({1 - {\rho_{2}}\left({1 - \frac{1}{{1 + {B^{- 1}}{T_{\text{FR},1}}{v^{- \alpha }}}}} \right)} \right)} \right]{\text{vdv}}}} \right) \end{array}$$
(52)
$$\begin{array}{@{}rcl@{}} &&\frac{2\pi {\lambda_{1}}}{{\mathcal A}_{1}}\int_{0}^{\infty} x\exp \left(- \frac{{T + {T_{\text{FR},1}}}}{\text{SNR}_{1}}\left({x,1} \right) \right) \\ &&\exp \left\{ - \pi {\lambda_{1}}{x^{2}}\left[ {1 + 2\xi \left({\rho_{1}^{\mathrm{e}},{\rho_{1}},T,{T_{\text{FR},1}},\alpha} \right)} \right] - \pi {\lambda_{2}}{{\left({\frac{{{P_{2}}B}}{{{P_{1}}}}} \right)}^{\frac{2}{\alpha }}}{x^{2}}\left[ {1 + 2\xi \left({\rho_{2}^{\mathrm{i}},{\rho_{2}},{B^{- 1}}T,{B^{- 1}}{T_{\text{FR},1}},\alpha} \right)} \right] \right\}dx \end{array}$$
(53)
$$\begin{array}{@{}rcl@{}} &&\mathbb{P}\left({\tilde N_{k}^{\mathrm{i}} = n} \right) = \frac{{{{3.5}^{3.5}}}}{{n!}}\frac{{\Gamma \left({n + 4.5} \right)}}{{\Gamma \left({3.5} \right)}}{\left({\frac{{{\lambda_{u}}{{\mathcal A}_{k}}{p_{c,k}}\left({{T_{\text{FR},k}}} \right)}}{{{\lambda_{k}}}}} \right)^{n}}{\left({3.5 + \frac{{{\lambda_{u}}{{\mathcal A}_{k}}{p_{c,k}}\left({{T_{\text{FR},k}}} \right)}}{{{\lambda_{k}}}}} \right)^{- \left({n + 4.5} \right)}} \end{array}$$
(54)
$$\begin{array}{@{}rcl@{}} &&\mathbb{P}\left({\tilde N_{k}^{\mathrm{e}} = n} \right) = \frac{{{{3.5}^{3.5}}}}{{n!}}\frac{{\Gamma \left({n + 4.5} \right)}}{{\Gamma \left({3.5} \right)}}{\left({\frac{{{\lambda_{u}}{{\mathcal A}_{k}}\left({1 - {p_{c,k}}\left({{T_{\text{FR},k}}} \right)} \right)}}{{{\lambda_{k}}}}} \right)^{n}}{\left({3.5 + \frac{{{\lambda_{u}}{{\mathcal A}_{k}}\left({1 - {p_{c,k}}\left({{T_{\text{FR},k}}} \right)} \right)}}{{{\lambda_{k}}}}} \right)^{- \left({n + 4.5} \right)}} \end{array}$$
(55)

Now conditioning on the second term of the numerator and following the method of Theorem 1, the joint Laplace transform $$\mathcal {L}\left ({{{\hat s}_{1}},{s_{1}},{{\hat s}_{2}},{s_{2}}} \right)$$ of $$\hat I_{x,1}^{\mathrm {e}}$$, I x,1, $$\hat I_{x,2}^{\mathrm {i}}$$, and I x,2 evaluated at $$\left ({{x^{\alpha } }P_{1}^{- 1}T,{x^{\alpha } }P_{1}^{- 1}{T_{\text {FR},1}},{x^{\alpha } }P_{1}^{- 1}T,} {{x^{\alpha } }P_{1}^{- 1}{T_{\text {FR},1}}} \right)$$ is

$${} \begin{array}{l} \mathcal{L}\left({{{\hat s}_{1}},{s_{1}},{{\hat s}_{2}},{s_{2}}} \right)\\ = \mathbb{E}\left[ {\exp \left({ - {{\hat s}_{1}}\hat I_{x,1}^{\mathrm{e}} - {s_{1}}{I_{x,1}} - {{\hat s}_{2}}\hat I_{x,2}^{\mathrm{i}} - {s_{2}}{I_{x,2}}} \right)} \right]\\ = \mathbb{E}\left[ {\exp \left({ - {{\hat s}_{1}}\hat I_{x,1}^{\mathrm{e}} - {s_{1}}{I_{x,1}}} \right)} \right]\mathbb{E}\left[ {\exp \left({ - {{\hat s}_{2}}\hat I_{x,2}^{\mathrm{i}} - {s_{2}}{I_{x,2}}} \right)} \right] \end{array}$$

Hence, $$\mathcal {L}\left ({{x^{\alpha } }P_{1}^{- 1}T,{x^{\alpha } }P_{1}^{- 1}{T_{\text {FR},1}},{x^{\alpha } }P_{1}^{- 1}T,{x^{\alpha } }P_{1}^{- 1}{T_{\text {FR},1}}} \right)$$ is given by (52).

Deconditioning on x, the second term of the numerator can be expressed as (53), where ξ(a,b,c,d,e) is originally defined in Theorem 2.

The denominator of (50) also can be obtained from (9), we have 1−p c,1(T FR,1). Plugging back into (50), the conditional coverage probability given in Theorem 3 is obtained.

## Appendix D

### Proof of Lemma 4

Let $$\tilde N_{k}^{\mathrm {i}}$$ and $$\tilde N_{k}^{\mathrm {e}}$$ be the numbers of other interior and edge users, conditioned on the typical user being associated with that BS in the kth tier, respectively. Their PMFs are derived in a similar way in , given by (54) and (55). □

According to the proposed random resource allocation scheme, $$t_{k}^{\mathrm {i}}$$ and $$t_{k}^{\mathrm {e}}$$ satisfy

$$t_{1}^{\mathrm{i}} = \left\{ \begin{array}{ll} 1 & {\text{if}}\ {\tilde{N}}_{1}^{\mathrm{i}} + 1 \le \left| {\mathbf{C}_{1}} \right|\\ \left| {{\mathbf{C}_{1}}} \right| \left/ {\tilde{N}}_{1}^{\mathrm{i}} + 1\right. & {\text{otherwise}} \end{array} \right.,$$
$$t_{1}^{\mathrm{e}} = \left\{ \begin{array}{ll} 1 & {\text{if}}\ {\tilde{N}}_{1}^{\mathrm{e}} + 1 \le \left| {{\mathbf{C}_{2}}} \right|\\ \left| {{\mathbf{C}_{2}}} \right| \left/ {\tilde{N}}_{1}^{\mathrm{e}} + 1\right. & {\text{otherwise}} \end{array} \right.,$$
$$t_{2}^{\mathrm{i}} = \left\{ \begin{array}{ll} 1 & {\text{if}}\ {\tilde{N}}_{2}^{\mathrm{i}} + 1 \le \left| {{\mathbf{C}_{2}}} \right|\\ \left| {{\mathbf{C}_{2}}} \right| \!\left/ {\tilde{N}}_{2}^{\mathrm{i}} + 1\right. & {\text{otherwise}} \end{array} \right.,$$
$$t_{2}^{\mathrm{e}} = \left\{ \begin{array}{ll} 1 & {\text{if}}\ {\tilde{N}}_{2}^{\mathrm{e}} + 1 \le \left| {{\mathbf{C}_{1}}} \right|\\ \left| {{\mathbf{C}_{1}}} \right| \!\left/ {\tilde{N}}_{2}^{\mathrm{e}} + 1\right. & {\text{otherwise}} \end{array} \right..$$

Thus, $$\mathbb {E}\left [ {t_{1}^{\mathrm {i}}} \right ]$$ can be derived as

\begin{aligned} \mathbb{E}\left[ {t_{1}^{\mathrm{i}}} \right] &= \sum\limits_{n = 0}^{\left| {{\mathbf{C}_{1}}} \right| - 1} {\left({{\tilde{N}}_{1}^{\mathrm{i}} = n} \right)} + \sum\limits_{n = \left| {{\mathbf{C}_{1}}} \right|}^{\infty} {\frac{{\left| {{\mathbf{C}_{1}}} \right|}}{{n + 1}}\left({{\tilde{N}}_{1}^{\mathrm{i}} = n} \right)} \\ &= \sum\limits_{n = 0}^{\infty} {\frac{{\left| {{\mathbf{C}_{1}}} \right|}}{{n + 1}}\left({{\tilde{N}}_{1}^{\mathrm{i}} = n} \right)} + \sum\limits_{n = 0}^{\left| {{\mathbf{C}_{1}}} \right| - 1} {\left({1 - \frac{{\left| {{\mathbf{C}_{1}}} \right|}}{{n + 1}}} \right)\left({{\tilde{N}}_{1}^{\mathrm{i}} = n} \right)} \end{aligned}
(56)

Besides, the mean proportion of sub-channel allocated to the macro interior user is derived by (57), where (a) follows from the change of variables i = n + 1. Plugging back into (56), we have

$$\begin{array}{l} \mathbb{E}\left[ {t_{1}^{\mathrm{i}}} \right] = \frac{{{\lambda_{1}}\left| {{\mathbf{C}_{1}}} \right|}}{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}\left[ {1 - {{\left({1 + \frac{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}{{3.5{\lambda_{1}}}}} \right)}^{- 3.5}}} \right]\\ \qquad \qquad + \sum\limits_{n = 0}^{\left| {{\mathbf{C}_{1}}} \right| - 1} {\left({1 - \frac{{\left| {{\mathbf{C}_{1}}} \right|}}{{n + 1}}} \right)\mathbb{P}\left({{\tilde{N}}_{1}^{\mathrm{i}} = n} \right)} \end{array}$$

Following the same derivation procedure, the $$\mathbb {E}\left [ {t_{1}^{\mathrm {e}}} \right ]$$, $$\mathbb {E}\left [ {t_{2}^{\mathrm {i}}} \right ]$$, and $$\mathbb {E}\left [ {t_{2}^{\mathrm {e}}} \right ]$$ can be obtained.

\begin{aligned} & \mathbb{E} \left[ {\frac{1}{{{\tilde{N}}_{1}^{\mathrm{i}} + 1}}} \right] = \sum\limits_{n = 0}^{\infty} {\frac{1}{{n + 1}}\mathbb{P}\left({{\tilde{N}}_{1}^{\mathrm{i}} = n} \right)} \\ & \overset{(a)}{=} \sum\limits_{i = 1}^{\infty} \frac{{{{3.5}^{3.5}}}}{{i!}}\frac{{\Gamma \left({n + 3.5} \right)}}{{\Gamma \left({3.5} \right)}}{{\left({\frac{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}{{{\lambda_{1}}}}} \right)}^{i - 1}} \\ &\quad\times {{\left({3.5 + \frac{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}{{{\lambda_{1}}}}} \right)}^{- \left({i + 3.5} \right)}} \\ & = \frac{{{\lambda_{1}}}}{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}\sum\limits_{i = 1}^{\infty} \frac{{{{3.5}^{3.5}}}}{{i!}}\frac{{\Gamma \left({n + 3.5} \right)}}{{\Gamma \left({3.5} \right)}} \\&\quad\times{{\left({\frac{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}{{{\lambda_{1}}}}} \right)}^{i}}{{\left({3.5 + \frac{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}{{{\lambda_{1}}}}} \right)}^{- \left({i + 3.5} \right)}} \\ & = \frac{{{\lambda_{1}}}}{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}\left[ {\sum\limits_{i = 0}^{\infty} {\mathbb{P}\left({N_{k}^{\mathrm{i}} = i} \right) - \mathbb{P}\left({N_{k}^{\mathrm{i}} = 0} \right)}} \right]\\ & = \frac{{{\lambda_{1}}}}{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}\left[ {1 - {{\left({1 + \frac{{{\lambda_{u}}{{\mathcal A}_{1}}{p_{c,1}}\left({{T_{\text{FR},1}}} \right)}}{{3.5{\lambda_{1}}}}} \right)}^{- 3.5}}} \right] \end{aligned}
(57)

## References

1. 1

JG Andrews, S Buzzi, W Choi, SV Hanly, A Lozano, ACK Soong, JC Zhang, What will 5g be? IEEE J. Sel. Areas Commun.32(6), 1065–1082 (2014).

2. 2

E Hossain, M Hasan, 5g cellular: key enabling technologies and research challenges. IEEE Instrum. Meas. Mag.18(3), 11–21 (2015).

3. 3

X Lin, JG Andrews, A Ghosh, Modeling, analysis and design for carrier aggregation in heterogeneous cellular networks. IEEE Trans. Commun.61(9), 4002–4015 (2013).

4. 4

A Damnjanovic, J Montojo, Y Wei, T Ji, T Luo, M Vajapeyam, T Yoo, O Song, D Malladi, A survey on 3gpp heterogeneous networks. IEEE Wirel. Commun.18(3), 10–21 (2011).

5. 5

Kyocera, Potential performance of range expansion in macro-pico deployment (r1-104355). 3GPP TSG RAN WG1 Meeting-62, (Madrid, 2010).

6. 6

D Lopez-Perez, I Guvenc, Gdl Roche, M Kountouris, TQS Quek, J Zhang, Enhanced intercell interference coordination challenges in heterogeneous networks. IEEE Wirel. Commun. Mag.18(3), 22–30 (2011).

7. 7

S Mukherjee, I Guvenc, in Proceedings of 2011 Conference Record of the Forty Fifth Asilomar Conference on Signals, Systems and Computers (ASILOMAR): 6-9 Nov. 2011. Effects of range expansion and interference coordination on capacity and fairness in heterogeneous networks (IEEE, Pacific Grove, 2011), pp. 1855–1859.

8. 8

K Okino, T Nakayama, C Yamazaki, H Sato, Y Kusano, in Proceedings of 2011 IEEE International Conference on Communications Workshops (ICC): 5-9 June 2011. Pico cell range expansion with interference mitigation toward lte-advanced heterogeneous networks (IEEE, Kyoto, 2011), pp. 1–5.

9. 9

A Merwaday, S Mukherjee, I Guvenc, in Proceedings of 2013 IEEE Global Communications Conference (GLOBECOM): 9-13 Dec. 2013. On the capacity analysis of hetnets with range expansion and eicic (IEEE, Atlanta, 2013), pp. 4257–4262.

10. 10

M Cierny, H Wang, R Wichman, Z Ding, C Wijting, On number of almost blank subframes in heterogeneous cellular networks. IEEE Trans. Wirel. Commun.12(10), 5061–5073 (2013).

11. 11

Panasonic:Performance study on abs with reduced macro power (r1-113806). 3GPP TSG-RAN WG1 #67, (San Francisco, 2011).

12. 12

A Merwaday, S Mukherjee, I Guvenc, in Proceedings of 2014 IEEE Wireless Communications and Networking Conference (WCNC): 6-9 April 2014. Hetnet capacity with reduced power subframes (IEEE, Istanbul, 2014), pp. 1380–1385.

13. 13

H Hu, J Weng, J Zhang, Coverage performance analysis of feicic low-power subframes. IEEE Trans. Wireless Commun.15(8), 5603–5614 (2016).

14. 14

S Singh, JG Andrews, Joint resource partitioning and offloading in heterogeneous cellular networks. IEEE Trans. Wirel. Commun.13(2), 888–901 (2014).

15. 15

Y Dhungana, C Tellambura, Multichannel analysis of cell range expansion and resource partitioning in two-tier heterogeneous cellular networks. IEEE Trans. Wirel. Commun.15(3), 2394–2406 (2016).

16. 16

W Tang, S Feng, Y Liu, MC Reed, in Proceedings of 2015 IEEE Global Communications Conference (GLOBECOM): 6-10 Dec. 2015. Joint low-power transmit and cell association in heterogeneous networks (IEEE, San Diego, 2015), pp. 1–6.

17. 17

B Xie, BZ Zhang, RQ Hu, Y Qian, in Proceedings of 2016 IEEE International Conference on Communications (ICC): 22-27 May 2016. Spectral efficiency analysis in wireless heterogeneous networks (IEEE, Kuala Lumpur, 2016), pp. 1–6.

18. 18

Q Li, RQ Hu, Y Xu, Y Qian, Optimal fractional frequency reuse and power control in the heterogeneous wireless networks. IEEE Trans. Wirel. Commun.12(6), 2658–2668 (2013).

19. 19

H ElSawy, E Hossain, in Proceedings of 2013 IEEE Global Communications Conference (GLOBECOM): 9-13 Dec. 2013. Channel assignment and opportunistic spectrum access in two-tier cellular networks with cognitive small cells (IEEE, Atlanta, 2013), pp. 4477–4482.

20. 20

RQ Hu, Y Qian, An energy efficient and spectrum efficient wireless heterogeneous network framework for 5g systems. IEEE Commun. Mag.52(5), 94–101 (2014).

21. 21

HS Dhillon, RK Ganti, F Baccelli, JG Andrews, Modeling and analysis of k-tier downlink heterogeneous cellular networks. IEEE J. Sel. Areas Commun.30(3), 550–560 (2012).

22. 22

H-S Jo, YJ Sang, P Xia, JG Andrews, Heterogeneous cellular networks with flexible cell association: a comprehensive downlink sinr analysis. IEEE Trans. Wirel. Commun.11(10), 3484–3495 (2012).

23. 23

24. 24

T Novlan, RK Ganti, A Ghosh, JG Andrews, Analytical evaluation of fractional frequency reuse for ofdma cellular networks. IEEE Trans. Wirel. Commun.10(12), 4294–4305 (2011).

25. 25

T Novlan, RK Ganti, A Ghosh, JG Andrews, Analytical evaluation of fractional frequency reuse for heterogeneous cellular networks. IEEE Trans. Commun.60(7), 2029–2039 (2012).

26. 26

J-S Ferenc, Z Neda, On the size distribution of poisson-voronoi cells. Phys. A Stat. Mech. Appl.385(2), 518–526 (2007).

27. 27

TS Rappaport, Wireless communications principles and practice, 2nd Edition (Prentice Hall, New Jersey, 2002).

28. 28

3GPP TR36.814, Further advancements for e-utra physical layer aspects (v9.0.0) (2010).

29. 29

D Stoyan, W Kendall, J Mecke, Stochastic geometry and its applications (Wiley, West Sussex, 1996).

## Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grant No. 61271263).

## Author information

LG and SC conceived the proposed scheme. SC conducted the detailed derivation to evaluate the performance of the proposed scheme and wrote the manuscript. LG and ZS reviewed the manuscript. All authors have read and approved the final manuscript.

Correspondence to Shanya Cong.

## Ethics declarations

### Competing interests

The authors declare that they have no competing interests.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions
{"ft_lang_label":"__label__en","ft_lang_prob":0.86749357,"math_prob":0.999894,"size":59206,"snap":"2019-51-2020-05","text_gpt3_token_len":15599,"char_repetition_ratio":0.16671734,"word_repetition_ratio":0.15526684,"special_character_ratio":0.27748877,"punctuation_ratio":0.11785621,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.9999429,"pos_list":[0,1,2],"im_url_duplicate_count":[null,4,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2020-01-19T07:57:10Z\",\"WARC-Record-ID\":\"<urn:uuid:9ab796c3-dd9f-414d-b0fc-f6706a0cd581>\",\"Content-Length\":\"347469\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:c03a9c68-6c70-48fc-848d-be4a6b9e30f6>\",\"WARC-Concurrent-To\":\"<urn:uuid:f8c001ab-a6e5-49fa-af67-d97fb5ccf02f>\",\"WARC-IP-Address\":\"151.101.248.95\",\"WARC-Target-URI\":\"https://jwcn-eurasipjournals.springeropen.com/articles/10.1186/s13638-017-0954-8\",\"WARC-Payload-Digest\":\"sha1:AIQ3F7M4KMXCHOQOCIP7YBC3VFUUJY7B\",\"WARC-Block-Digest\":\"sha1:77EJFCEAGMQ3KUZW4KO4WPVKHNLPJ7AC\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2020/CC-MAIN-2020-05/CC-MAIN-2020-05_segments_1579250594333.5_warc_CC-MAIN-20200119064802-20200119092802-00409.warc.gz\"}"}