Outage probability of fixed-gain dual-hop relay selection channels with heterogeneous fading

In this paper, we develop a framework that can be used to analyze the outage probability of a dual-hop fixed-gain amplify-and-forward relay system performing single relay selection. We consider three selection strategies: (1) first-hop selection, whereby the relay is chosen to maximize the signal-to-noise ratio (SNR) at the relay irrespective of the second-hop channels; (2) second-hop selection, in which the relay is chosen such that the second-hop SNR is maximized; and (3) dual-hop selection, where the end-to-end SNR is maximized. The proposed analytical framework is capable of treating systems operating in heterogeneous fading conditions, where all or a subset of the source-relay and relay-destination channels has non-identical distributions or even experience completely different fading processes. We apply the framework to calculate new exact series representations of the outage probability for the three aforementioned cases of relay selection when all channels adhere to the Nakagami-m model. Our analysis is corroborated by simulations. Finally, we provide a discussion of how the proposed framework can be applied to analyze other fading configurations using the techniques detailed herein.


Introduction
The paradigm of cooperative communication has been shown to provide significant advantages in cellular networks (and many other applications) for more than a decade [1]. With the deployment of LTE networks across the world in recent years [2] and the shifting of focus in the research community to 5G cellular technology [3], the ideas surrounding cooperative communication are evolving. Of particular interest in the discussion of 5G is relayaided communication. The heterogeneous structure of 5G networks will consist of layers of macrocells, small cells, relays, and device-to-device (D2D) networks [4]. In areas lacking a wired backhaul infrastructure, relays are, and will continue to be, deployed to act as picocell base stations to local user equipments (UEs) and to mimic (for the most part) a UE in the view of the macrocell base station. Within the D2D networking paradigm, relays will form a key component in the control infrastructure of a 5G network, linking D2D networks to the macrocell base station [4]. This idea can be further extended to connect D2D *Correspondence: justin.coon@eng.ox.ac.uk 2 Department of Engineering Science, University of Oxford, Parks Road, OX1 3PJ Oxford, UK Full list of author information is available at the end of the article clusters to the macrocell by using one of the constituent UEs in the cluster as a relay [5,6].
Fixed-gain amplify-and-forward (FGAF) relay systems have attracted a lot of attention recently due to their low complexity in practical implementation. To date, numerous efforts have been devoted to the performance analysis of dual-hop FGAF relay systems (see [7] and the references therein), including error rate, capacity, and outage calculations. With regard to the latter, closed-form outage probability expressions for dual-hop FGAF were derived in [8] for Rayleigh fading channels, in [7,9,10] for Nakagamim channels, and in [11] for Rician channels. In [12], a moment generation function (MGF) approach was used to analyze the performance of dual-hop FGAF systems operating in arbitrary, or generalized, fading conditions. Heterogeneous fading conditions were also considered in [13] in the context of systems employing maximum ratio combining and transmit antenna selection at the relay. Many of these reported results were recently generalized and incorporated into a multi-hop framework in [14].
Relay selection, the process by which one relay is chosen among all the available relays to forward the data from the source to the destination [15], has been shown to be a simple and effective scheme that offers the same diversity order as cooperative schemes that use all available relays [16]. Recently, this concept has been considered within the 5G framework [17]. Three main selection schemes have been considered in the literature: (1) first-hop selection (FHS), in which the selected relay maximizes the signal-to-noise ratio (SNR) at the first hop; 2) second-hop selection (SHS), in which the second-hop SNR is maximized; and (3) dual-hop selection (DHS), where the the end-to-end SNR is the objective to be optimized 1 . Several published works have analyzed the outage probability for dual-hop AF relay selection systems [18][19][20][21]. However, the outage probability analysis carried out to date for such systems presents three main limitations: 1. Much of the published literature considers FHS or DHS, neglecting the case of SHS despite the fact that the latter is an important practical model. 2. Source-relay (S-R) and relay-destination (R-D) channels have largely been assumed to follow Rayleigh, Nakagami-m, or Rician distributions in relay selection systems, with some important fading distributions having been largely ignored (e.g., Weibull and Hoyt). 3. Heterogeneous fading conditions have generally not been studied for relay selection channels, although such scenarios may arise in a plethora of applications, such as a cellular base station transmitting to a relay located in an enterprise zone, which then conveys the signal to users inside an office block.
In this paper, we aim to overcome the aforementioned limitations by presenting a framework for analyzing the outage probability of FGAF systems employing relay selection in heterogeneous channel conditions. The proposed framework applies generally for arbitrary configurations of S-R/R-D channels, e.g., Nakagami-m/Rician, Hoyt/Weibull, etc. In order to facilitate exposition, we focus on Nakagami-m fading in this work and derive the outage probability for systems operating in heterogeneous channels for FHS, SHS, and DHS selection strategies. We also provide complete diversity results for each selection mechanism and qualitatively discuss the application of the framework to other fading configurations.
The rest of the paper is organized as follows. The dualhop FGAF system model is presented in the "System model" section. The analytical framework is detailed in the "Analytical outage probability framework" section, which is then applied in the "Analysis for Nakagami-m" section to study heterogeneous Nakagami-m scenarios. The analysis is corroborated with numerical results in the "Simulation results and discussion" section. Extending the framework to other more general fading configurations is discussed and concluding remarks are given in the "Conclusions" section.

System model
Consider a dual-hop transmission system where N ≥ 1 relays are available to forward information from a source node to a destination node using an FGAF protocol (see Fig. 1). There is no direct S-D link. Suppose the nth relay is selected among the N available relays. A complex symbol with unit average power is transmitted from the source to the selected relay over a flat fading channel h 1,n . The received signal at this relay is then amplified with a fixed gain G, then transmitted from the relay to the destination over a flat fading channel h 2,n . We assume that h 1,n and h 2,n are independent but not necessarily identically distributed random variables.
Defineγ as a reference average received SNR parameter. Clearly, the average received SNR values corresponding to different hops will be different in practice. This is accounted for in the model by denoting the variances of the zero mean additive white Gaussian noise at the nth relay and the destination as σ 2 1,n and σ 2 2 , respectively, and defining parameters ρ 1,n and ρ 2 such thatγ = ρ 1,n /σ 2 1,n for n = 1, . . . , N andγ = ρ 2 /σ 2 2 . Thus, the SNR at the nth relay of the first hop is given by X n = |h 1,n | 2 ρ 1,nγ , and the SNR corresponding to the local channel at the destination (i.e., not the end-to-end SNR) is Y n = |h 2,n | 2 ρ 2γ . For the relay selection process, a single relay is chosen among the N available relays according to the FHS, SHS, or DHS selection criteria. For FHS, the selection is made at the source to maximize the received SNR at the relay. The chosen relay can be notified of its selection through, for example, a dedicated control channel. Exact details of how control signaling might be performed are beyond the scope of this paper. Consequently, under the FHS approach, the kth relay is selected such that k = arg max{X n }. Analogously, for SHS, we have k =

Analytical outage probability framework
Suppose the kth relay is selected to convey the source message to the destination. Irrespective of the selection criteria, we can write the end-to-end SNR as where C k =γ /(G 2 ρ 1k ). For DHS, the SNR outage probability 2 can be written as where z is the outage threshold. For both FHS and SHS, we let ℘ h,k denote the probability that the kth relay is selected according to criteria h ∈ {FHS, SHS}. Applying the law of total probability, we can write Note that we have slightly abused notation here since P(γ k < z) is conditioned on the selection mechanism (FHS or SHS) and the index of the selected relay k . Nevertheless, the important point is that γ k is a function of two independent variates, X k and Y k . This enables us to write a general formula for the cumulative distribution function (c.d.f.) of γ k , which we encompass in the following lemma (the proof is given in the Appendix).

Lemma 1.
Let γ k be defined as in (1). Then, the c.d.f. of γ k can be written as where c is a real constant 3 , and Although the lemma seems somewhat inaccessible due to the integrals, we note that it provides an exact expression as a functional of the first-and second-hop fading density functions, regardless of the underlying distributions. The only stipulation is that the channels corresponding to the two hops are statistically independent, a condition that is often, if not always, satisfied in practice. We will revisit that and apply this lemma to a great effect in the next section.
Turning our attention to FHS and SHS, it remains to calculate ℘ h,k . For FHS, this quantity is given by where it is assumed that {X n } are statistically independent. The expression for ℘ SHS,k is the same, but with Y replacing X.
Applying Lemma 1, we can now write the outage probability for DHS as and the outage probability for FHS and SHS as Consequently, to apply this framework, we must be able to calculate the functional λ k (s, z) as well as the contour integral Since λ k (s, z) is determined by the channel distributions, it can often be calculated easily. For the contour integral, we will exploit the pole structure of λ k (s, z) and apply the residue theorem from complex analysis to solve the integral, thus yielding series expressions of the outage probabilities in general. These technical details will be the focus of the next section.

Analysis for Nakagami-m
For the purpose of presenting how the analytical framework can be applied, we assume the channels for the two hops adhere to the same class of distribution, but may have different parameters. The same approach can be used to study arbitrary combinations of fading distributions for the different channels. Suppose all channels adhere to a Nakagami-m smallscale fading model. In this case, the probability density function (p.d.f.) of X n is given by and the p.d.f. of Y n is defined similarly. The corresponding c.d.f. of X n is given by where γ (a, x) = x 0 t a−1 e −t dt is the lower incomplete gamma function [22]. The parameter m 1,n > 1/2 is the shape parameter of the distribution (denoted by m 2,n for the p.d.f. of Y n ) while θ 1,n = E[ |h 1,n | 2 ] is effectively the scale parameter (correspondingly, θ 2,n for Y n ).

FHS
To calculate P o (z) for the FHS scheme, let us first consider the probability that the kth relay is selected. Using the absolutely convergent series expansion along with (6) enables us to write where the summation is (N − 1)-fold, with the indices being the set L k = {l 1 , . . . , l k−1 , l k+1 , . . . , l N }, with l n = 0, 1, 2, . . . for all n. In addition, β 1 = N n=1 ρ 1,n /θ 1,n and m 1,k = N n=1 m 1,n + N n=1, n =k l n . Now consider the calculation of λ k (s, z) in Lemma 1. The second integral in (5) is a gamma function. The calculation of the first integral follows from an application of the theory of order statistics (cf. [12, eq. (17)]) to obtain an expression for f X k . Using the series expansion for the lower incomplete gamma function mentioned above, the result can be written as where is the confluent hypergeometric function (Tricomi's solution) [23].
We are now in a position to calculate the integral ι k (z) and thus the outage probability, which is made much easier than it may first appear through the use of the residue theorem [24]. More specifically, we note that for a fixed third argument, U(·, ·, ·) is entire in its first two arguments [23]. Furthermore, for some complex number a, (x + a) has simple poles at x = −a − j for j = 0, 1, . . .. It follows from the Mellin inversion theorem [25] that the path of integration that defines ι k (z) is a line extending from −i∞ to i∞ just to the right of the imaginary axis, i.e., c > 0 in (9). By noting that we can construct a path in the left half splane such that as |s| → ∞ then | λ k (s, z) ds| → 0 along this path 4 , we can close the path of integration over which ι k (z) is defined with a semi-circular arc traversing this half plane counter-clockwise from i∞ to −i∞, thereby encircling the poles of λ k (s, z). Finally, we can employ the residue theorem to obtain a series expression for ι k (z). For the case where m 2,k is not an integer 5 , we end up with the expression Applying (16) to (8) gives the outage probability for FHS in Nakagami-m channels.
Although the outage probability is expressed as an infinite series, we will verify later that truncating the series to only a few terms provides an excellent approximation. Moreover, the hypergeometric function U(·, ·, x) is well behaved for small arguments 6 , and thus, it is possible to construct an asymptotic series of Poincaré type forγ → ∞. Discussion of such an activity is rather involved and would detract from the focus of this paper; hence, we defer it to a future contribution.
We also note that the outage expression given above applies to the general case where each of the N available channels at a given hop is subject to non-identical Nakagami-m fading with different scale and shape parameters, i.e., m 1,n = m 1,q for n = q and so on. The result is therefore more general than that given in [19], where the outage probabilities are obtained for the case where m 1,n = m 1 and m 2,n = m 2 for all n.
To conclude this subsection, we provide a brief note on the case where m 2,k is an integer. Under this condition, λ k (s, z) has second-order poles at s = −m 2,k , −m 2,k − 1, . . .. The residue theorem can still be applied, but the calculations become cumbersome, and careful bookkeeping must be enforced in order to ensure the resulting series expansion is accurate. This is particularly true wheñ m 1,n has certain properties (such as being an integer), as this will cause the structure of the expansion to vary. For our practical purposes, however, it is much more intuitive and efficient to evaluate the outage probability using noninteger values of m 2,k to approximate integers. We will show through simulations that this simple approximation works very well in practice.

SHS
The outage probability for SHS can be obtained by following a similar approach as was outlined above for FHS. Without presenting the details of the calculations, we directly present the outage probability expression to be where β 2 andm 2,k are defined in a similar manner to β 1 andm 1,k .

Single relay
Letting N = 1 corresponds to the case where there is a single relay, and thus, selection does not take place, i.e., the system model is the canonical dual-hop FGAF topology with no direct link. The outage probability can be calculated from Lemma 1 where, clearly, k = 1 and f X k and f Y k take the form given in (10). The resulting expression for the outage probability is where θ q and m q are the scale and shape parameters for the Nakagami-m channel at the qth hop, and ρ q = σ 2 qγ where σ 2 1 and σ 2 2 are the noise variances at the relay and the destination, respectively.
Although an outage expression for dual-hop FGAF systems with a single relay operating in Nakagami-m fading channels was provided in [7, eq. (20)], the result presented here is simpler and easier to evaluate numerically. To the best of our knowledge, (18) has not been reported in the literature.

DHS
The outage probability for DHS is obtained by substituting the expression for P o (z) given by (18) for the quantity in the brackets of (7), i.e., 1− 1 2πi c+i∞ c−i∞ λ k (s, z) = P o (z). In doing so, all parameters should be replaced by those corresponding to the kth path for k = 1, . . . , N. We omit the explicit expression for brevity.

Diversity analysis
As mentioned above, one may obtain asymptotic expansions for the outage probabilities corresponding to the abovementioned scenarios. Although a full asymptotic theory is rather involved and beyond the scope of this paper, it is relatively straightforward to obtain the diversity order for FHS, SHS, and DHS. This can be done in a rather cumbersome manner by invoking the series expansion representation of U(·, ·, x) (cf. [23]) and retaining the leading order inγ . It is, however, much more straightforward to use the theorem recently reported in [14]. For the two-hop relay selection schemes treated here, this theorem states that the outage probability of the system decays asymptotically like O (lnγ ) κ−1γ (s 0 ) where s 0 is the pole (of order κ) of the function that lies furthest to the right in the s-plane, with Z 1 and Z 2 respectively denoting the local first-and second-hop SNRs after selection. A consequence of this theorem is that the diversity order, defined in the usual way as is simply given by (see [14] for details) First, consider the FHS scheme, and suppose all shape parameters in the first hop are distinct 7 . For this case, let such that Z 1 = X n and Z 2 = Y n . The p.d.f. of Z 1 is where f X n (z) and P(X k ≤ z) are given in (10) and (11), respectively. Using (12), it is straightforward to show that where ζ = N n=1 m 1,n . Consequently, the pole that lies furthest to the right in the s-plane is s = −ζ . Moreover, it is easy to see that which has a dominant pole at s = −m 2,n . In systems employing selection diversity, it is well known that outage is dominated by the worst case channel. Since all secondhop channels are independent, this pertains to the case where the selection process results in the second channel having the minimum shape parameter. It follows from this fact and the calculations above that the diversity order for FHS is d = min To calculate the diversity order for DHS, it is easiest to substitute (2) into (20), which leads to These results are summarized in Table 1. One can confirm that these diversity orders specialize to the previous published results when the scale parameters and shape

Simulation results and discussion
In this section, we present numerical results in an effort to validate the analysis presented above. For all of the results presented below, we consider the case where there are two relays available, i.e., N = 2, and the ρ and θ parameters are all set to 1 if not otherwise stated. Various channel configurations were considered; Table 2 provides an overview, and the corresponding diversity orders are given in Table 3. For configurations Nak-1 and Nak-2, we considered the case where there is only a single relay, and the channels in both hops are non-identically or identically distributed, respectively. For configurations Nak-3, Nak-4, and Nak-5, we considered systems with relay selection, where Nak-3 and Nak-4 represent the cases where the N channels are identically distributed at the second hop and the first hop, respectively, and Nak-5 describes the case where all channels at all relays in both hops are identically distributed. Nak-6 considered the generalized case where the fading parameters are arbitrarily chosen such that channels corresponding to each relay for each hop are non-identically distributed, and the θ and ρ parameters for this case are also arbitrarily chosen to be θ 1,1 = 0.6, θ 1,2 = 0.8, θ 2,1 = 0.9, θ 2,2 = 1.1, ρ 1,1 = 1, ρ 1,2 = 0.9, and ρ 2 = 0.8. Note that for the Nak-2 configuration, we have deliberately chosen the parameters such that m 2 is an integer. The effect of second-order poles in this case was discussed in the "Analysis for Nakagami-m" section. Here, instead of considering the residues at the second-order poles, we simply use a non-integer value of m 2 to obtain approximate numerical results. It will be shown that such an approximation yields reasonably accurate results compared with numerical simulations. Figures 2, 3, and 4 depict the results for channels Nak-1 to Nak-5 with FHS, SHS, and DHS, respectively. It is shown that the analytical results agree well with the simulations for all SNR values, even though the series expressions were truncated to ten terms. In addition, the diversity orders given in Table 3 can be verified from the  Table 3 Diversity orders corresponding to configurations outlined in Table 2 Configuration FHS SHS DHS figures. In particular, it is observed in Fig. 4 that the outage probability corresponding to Nak-5 outperforms that for Nak-2 at SNR values greater than 10 dB due to the higher diversity gain associated with Nak-5. An additional interesting feature that is brought out in these results is the importance of properly choosing the selection scheme according to the prevailing statistical channel conditions. Observing Figs. 2 and 3, we see that the Nak-4 relay selection system performs worse than the single-relay Nak-1 system for SNR values greater than 10 dB when FHS is employed (see Fig. 2). This is due to the comparatively low diversity order resulting from the selection process (see Table 3). However, if SHS is employed in the Nak-4 system, the diversity advantage offered by the selection scheme is realized, as shown in Fig. 3. Figure 5 illustrates the outage probability for the Nak-6 configuration. Again, the numerical results are shown to agree with the simulation results for all selection schemes. In addition, it is observed that for this particular configuration, FHS significantly outperforms SHS owing to the fact that the diversity order offered by the former is 2.4, whereas it is 0.9 in the latter case. The disparity in the diversity orders due to the use of different selection strategies provides significant insights into practical design for relay selection systems.
In Fig. 6, we present the numerical results obtained by setting the shape parameters to 1 (i.e., Rayleigh fading) at the second hop or at both hops, thus depicting Nakagami-Rayleigh or Rayleigh-Rayleigh fading scenarios. The results in Fig. 6 were obtained by using m 1,1 = 1.3 and m 1,2 = 2.5 at the first hop for the Nakagami-Rayleigh configuration. It is observed that for all three selection methods, the Nakagami-Rayleigh configuration provides better performance than the Rayleigh-Rayleigh mode, which satisfies intuition since Rayleigh channels exhibit more severe fading than Nakagami-m channels. It is also interesting to observe that SHS yields an outage gain compared to FHS in the Nakagami-Rayleigh channel, indicating the advantage of performing selection at the hop that experiences more severe fading.

Conclusions
The framework presented herein can be applied to analyze many practical system configurations. The key components that are required are derived from Lemma 1, whence we see that density functions of the effective channel gains X k and Y k are required along with integrals that are closely linked to their Mellin transforms. Certainly, having this information available makes the calculation of ℘ h,k possible in re (6) as a series expansion if not in a simple closed form.
Taking a step back, it is clear that the example discussed in the previous section was made possible by the pole structure of λ k (s, z), which in turn arose from the Mellin for some shape and scale parameters m and θ, the poles of which lie at s = −m(1 + j) for j = 0, 1, . . .. Thus, we can see that problems involving Weibull distributions will also experience a rich pole structure, and the full power of the residue theorem will be applicable.
Other fading distributions are also worth considering within this framework. Rician and Hoyt densities involve Bessel functions, which lead to hypergeometric functions ( 1 F 1 in the case of Rician and 2 F 1 in the case of Hoyt) through the evaluation of the integral x s f X (x) dx in (5). But these functions are either entire or meromorphic in the arguments and thus can be dealt with by calculating the residues just as was done for Nakagami-m fading. Although the method is straightforward, the exact results are often challenging to interpret in written form given only a few pages in which to do so, and thus, we have omitted them from this discussion.
It should be clear by now that the fading distributions corresponding to each hop in the network need not belong to the same class, e.g., Rician, Weibull, etc. Indeed, the power of the framework is that it is relatively agnostic to the functions f X and f Y , as long as the related integrals given in (5) exist and the result facilitates contour integration vis-à-vis (4).
As a final note, we have liberally made use of infinite series representations of special functions through this analysis. For the examples discussed in the previous section, these series were convergent and thus well behaved under the operation of integration. However, if one is concerned mostly with asymptotics (e.g., high SNR), this condition need not be satisfied. In fact, very  accurate approximations to the outage probability of some relatively complex systems can be attained from truncated asymptotic series. Again, we leave the details to another forum due to space restrictions. 1 The term partial selection is also used in the literature to describe relay selection based on channel knowledge pertaining to a single hop. 2 As is customary, we consider SNR outage since there is a one-to-one relationship between this performance metric, mutual information outage, and many symbol-error-rate outage expressions encountered in practice when flat fading is assumed. 3 This constant arises from the use of the Mellin inversion theorem [25] and basically separates the poles of the integrand in a particular and well-defined way for a given function λ k . The interval in which c lies is explored for specific examples in the next section. 4 This follows from Stirling's formula [22] and an application of Jordan's lemma [26]. 5 We will return to the case where m 2,k is an integer later. 6 Technically, there is a singularity at x = 0, but for the purposes of our system analysis, this condition will never hold. 7 When the shape parameters are not distinct, the same approach that is outlined in this section can be taken, but the algebra and differentiation become cumbersome.