### Abstract

Antenna correlation is generally viewed as an obstacle to realize the desired performance of a wireless system. In this article, we investigate the performance of partial relay selection in the presence of antenna correlation. We consider both channel state information (csi)-assisted and fixed gain amplify-and-forward (AF) relay schemes. The source and the destination are equipped with multiple antennas communicating via the best first hop signal-to-noise ratio (SNR) relay. We derived the closed form expression for outage probability, average symbol error rate (SER) for both schemes. Further, an exact expression is derived for the ergodic capacity in the csi-assisted relay case and an approximated expression is considered for the fixed gain case. Moreover, we provide simple asymptotic results and show that the diversity order of the system remains unchanged with the effect of antenna correlation for both types of relay schemes.

### Introduction

Two-hop amplify-and-forward (AF) relay networks have been investigated extensively in recent research [1-4]. The system with a source and a destination both equipped with multiple antennas communicating via a single antenna relay has received significant interest in most of the previous literature [5-11]. Different transmission and receive techniques were used and use of maximal ratio transmission (MRT) and maximal ratio combining (MRC) were among the most significant ones [5-9]. The analyzes in these cases were carried out with different fading channel environments for performance evaluation.

Antenna correlation is generally considered as a detrimental effect which degrades the system performance. To investigate this loss, several authors have studied the effect of antenna correlation in AF relay schemes. Authors in [7] have analyzed the channel state information (csi)-assisted AF relay network under antenna correlation with distinct eigenvalue distribution of correlation matrices and the fixed gain scheme has been considered in [12]. Then the general case of arbitrary distributed correlation matrix structures has been investigated by the authors in [9]. However, these evaluations are limited to the single source, relay and destination scenario.

It has been proven that the use of multiple relays with different selection methods can enhance the diversity and the performance [13-23]. There are several ways of selecting a relay for transmission. One method is referred to as the opportunistic relay selection [13,14] in which the relay with maximum instantaneous end-to-end signal-to-noise ratio (SNR) is considered. Synchronization is very important in this case. Another is the partial relay selection method, which can be carried out in two ways; by selecting either the first-hop relay [15,19,21] or the second hop relay [13,17,19] with the maximum instantaneous SNR. All these studies have been concentrated on the independent and identically distributed fading environments with some considering the effect of feedback delay.

#### Contribution

Although authors in the previous literature have studied the AF relay network under the effect of antenna correlation, all these works have been limited to single relay network. Hence, it motivated us to investigate the performance of partial relay selection with the effect of antenna correlation. We consider two types of AF relay schemes; csi-assisted and fixed gain relay. The exact closed form expressions for outage probability and average symbol error rate (SER) are derived for both schemes and an exact ergodic capacity expression is derived for the csi-assisted case and an approximation is found for the other case. Further, we study the system in high SNR and derive simple asymptotic expressions for outage probability and average SER for both cases. Our asymptotic analysis provide the depth in the system performance and it shows the variation of diversity gain. Finally, we give Monte Carlo simulations to verify our results.

### System model

Consider an AF relay network where a source (S) communicates with a destination (D) via the best relay (R). Both S and D are equipped with *n*_{s}and *n*_{d}antennas, respectively, and relays are equipped with a single antenna. Direct path
between source to destination is assumed to be unavailable due to heavy shadowing.
The csi is assumed to be available at S. When the csi is available at the transmitter, the optimal transmission scheme is
maximal ratio transmission (MRT) [24], hence, Suses MRT as the transmission scheme and destination uses MRC. We consider a system
where all the relays are homogeneously located having the same average SNR and we
further assume that S−R_{i}∀*i* channels are independent of each other. Source uses the csi to find the maximum SNR
relay from *L* number of relays in the first hop as,

where ||·||_{F} denotes the Frobenius norm and **h**_{si} is the *n*_{s}×1 channel vector between S−*R* and the elements of **h**_{si} are modeled as mutually correlated Ralyeigh fading entries. Let *n*_{s}×*n*_{s}correlation matrix at source be *Φ*_{s}, then
*E*[·] and (·)^{‡}denote the expectation operator and the Hermitian transpose, respectively. The communication
happens in two time slots as presented in numerous literature. During the first time
slot, S transmits the signal *x* to the selected relay R_{m} and the received signal at R_{m}is given as,

where *P*_{s} is the transmitted power and **h**_{sm} is given as in (1) and **w**_{s} is the MRT weight vector which is defined as
*V*_{m}variance at R_{m} is denoted as *v*_{m}. Then R_{m} multiplies the received signal by gain *G* and transmits to the Dand the received signal at Dis given as,

where *G* is defined differently for the two relay schemes and is given in the next section.
*P*_{r} is the transmitted power at R_{m} and 1×*n*_{d} channel vector between R_{m}−*D*is **h**_{md}and its elements are mutually correlated such that the correlation matrix at D is
**v**_{d} is noise vector at D and it elements are AWGN with *V*_{d}variance. Now Dperforms MRC to obtain the signal as,

where

*Notation:* Let
*Φ*_{s}be
*Φ*_{d}be

#### Statistics of SNRs

We can derive the probability density function (pdf) of *γ*_{2} as [25,26],

and the cumulative density function (cdf) of *γ*_{2}can be derived using

Now, we derive the cdf of *γ*_{1} as,

#### Proof

See Appendix 1. □

### Channel state information assisted relay

Here the relay uses the csi to amplify the received signal. It has been proven that
the csi-assisted relay outperforms the fixed gain relay in general. However, it has
a higher complexity in implementation when compared to the fixed gain one. In this
section, we derive exact expressions for outage probability, average SER and ergodic
capacity. First, we select the gain *G* for csi-assisted relay as,

Then the end-to-end SNR given in (5) can be rewritten as,

where *c*=1 for exact end-to-end SNR. The approximation holds for the medium to high SNR and
it provides a tight upper bound, we use the exact SNR to derive the outage probability
and ergodic capacity, and use the approximation for average SER.

#### Outage probability: csi-assisted relay

The outage probability is the probability that *γ*_{e}drops below a predefined threshold *Θ*and it is mathematically represented as,

We can derive the exact closed form expression for outage probability as,

where K_{1}(·) is the first order modified Bessel function of first kind.

#### Proof

Appendix 1. □

#### Outage probability: approximation

We simplify (12) by substituting *c*=0 to obtain the tight upper bound as given in (11) as,

#### Average SER

In this section, we derive the closed form expressions for average SER (*P*_{ser}) which is valid for several modulation schemes. According to technical literature
*P*_{ser} is defined as,

where *E*[·] is the expectation operator and *Q*(.) is the Gaussian Q function. The modulation schemes is defined by *a**b*, mainly BPSK (*a*=1, *b*=1), M-ary PAM (

Now, we substitute
*γ*_{e}as,

where

where *F*(*μ**ν*;*a*;*b*) is the Gauss hypergeometric function defined in ([27], Equation 9.10–9.13) and *Γ*(·) represents the Gamma function.

#### Ergodic capacity

The exact closed form expression for ergodic capacity has a significant importance
since it has not been derived even for a single user relay network with antenna correlation.
The ergodic capacity (*C*_{erg}) can be mathematically expressed as,

Closed form expression for *C*_{erg}can be derived as,

where

where
*χ*_{2}is given as,

and

where

#### Proof

Appendix 2 □

### High SNR analysis

Here we derive the high SNR expressions for the outage probability and the average
SER. Let
*ρ*_{2}=*μ**ρ*_{1}.

#### High SNR outage probability

We rewrite the (7) by expanding the

It is observed at high SNR the lower order terms (

We can further simplify (29) as,

Similarly, we can derive the

Simplification yields,

Following the same procedure as in ([30], Equations (A.09) and (A.10)), we can obtain the high SNR expression for

where,

It is observed from (33) that the diversity gain *G*_{d}= min*L**n*_{s}*n*_{d}. This shows that the performance of the system is dominated by one of the single
links unless *L**n*_{s}=*n*_{d}, hence, in order to fully utilize the resources, we like to keep *L**n*_{s}=*n*_{d}. It is observed here that the diversity gain of the system is not affected by the
antenna correlation. However, we have to note that this condition is only true for
the case where the correlation matrices have full rank.

#### High SNR average SER

Average SER in high SNR can be obtained using [31]

where *a* and *b* define the modulation scheme and *ψ*is given as in (34)–(36). *t*=min*L**n*_{s}*n*_{d}−1 and diversity gain *G*_{d}=*t* + 1.

### Fixed gain relay

In this section, we investigate the end to end performance of a dual-hop fixed gain network with multiple relays in the presence of antenna correlation at both ends. We derive closed-form expressions for outage probability, average SER, generalized moments of the end-to-end SNR and ergodic capacity. The asymptotic behavior of the system under high SNR is also considered. After performing some algebraic manipulations of (5), the end-to-end SNR for fixed gain can be expressed as,

According to literature, there are two common techniques for selecting a fixed gain
*G* in (38). If the gain is selected as,

then the constant *C*given by

If the gain is selected as,

the constant *C*given by

#### Outage probability

If outage probability is denoted by *P*_{out}, then

The closed form expression for outage probability can be derived as,

where *K*_{1}(.) is 1st order modified Bessel function of second type and *C* is the fixed gain.

#### Proof

Appendix 3 □

From Equation (8), cdf of *γ*_{1}can be written as,

Where,

and

#### Case 1

Using ([9], Equation 46) the expected value for random variable *γ*_{1} is calculated as,

Substituting to (42), the closed form expression for fixed gain can be derived as,

where *α*_{4}and *β*_{4} are as mentioned in (46) and (47), respectively.

#### Case 2

Substituting to (40)

where *Ei*(·) is the exponential integral.

#### Average SER

Then substituting

where *K*_{1}(.) and *K*_{0}(.) are the 1st and 2nd order modified Bessel functions of second type, respectively.

#### Generalized moments of SNR

In this section, we derive closed form expression for generalized moments of *γ*_{e} which is essential to the obtain ergodic capacity and the performance evaluation
of the system using the average output SNR and the degree of fading. Substituting
(44) into ([9], Equation 46) and after some mathematical manipulations with the help of ([27], Equation 6.643.3), we obtain

where *W*_{k,μ}(.) is the whittaker function defined in ([27], Equation 9.220.4).

#### Ergodic capacity

We derive a closed form expression for the ergodic capacity in multi relay network which is significant in determining the system performance especially in a correlated environment. Based on the literature ergodic capacity can be expressed as follows,

applying the expectation operator as in (54), an approximated result can be obtained as given below ([18], Equation 6)

By substituting (53) in to (55) for *h*=1, *h*=2 we can obtain approximated closed form expression for ergodic capacity.

#### High SNR analysis: outage probability

In this section, we analyze the fixed gain relay system in high SNR. Let *C*=*D**ρ*_{1}*ρ*_{2}=*μ**ρ*_{1}, and *z*=*Θ*/*ρ*_{1} then we can rewrite (44) as,

then (56) can be rewritten as,

where

and

Now, we expand the exponential function using Maclaurin series and Bessel function using ([27], Equation 8.446) in (57) to obtain,

Now it is observed that

where *N*=min(*L**n*_{s}*n*_{d}) and

where *ψ*(·) is Euler Psi function. It is observed from the fixed gain asymptotic outage expression
that the diversity gain of the system is similar to the csi-assisted relay scheme.

#### High SNR analysis: average SER

We can write the asymptotic average SER as [31],

where

where*Ψ* as in (58) and *a* and *b* define the modulation scheme. *t*=min*L**n*_{s}*n*_{d}−1 and diversity gain *G*_{d}=*t* + 1.

### Numerical analysis

Here we carry out the numerical analysis and verify our results using Monte Carlo
simulations. We use the exponential correlation matrix structure where (*i*,*j*)th element of the matrix *Φ*_{s}is
*Φ*_{d}is
*ρ*_{1}=*ρ*_{2}(*μ*=1) in all the cases shown in the figures. Fixed gain type 1 in (49) is used. Exponential
correlation matrices defined above have full rank. Hence, we obtain the desired diversity.

Figure 1 shows the outage probability variation with the average SNR of the first hop. Curves are plotted for different antenna configurations and correlation parameters. It is observed from the figure that the increase of the number of antennas, improves the outage probability. High SNR curves are also plotted where we can clearly see how diversity gain is varying. One can notice that the left three curves have a diversity gain of four and the two right most curves have a diversity gain of two. It is noticed that the increase of correlation decreases the performance. The outage probability variation for fixed gain relay is depicted in Figure 2. Here we see an improvement in the performance when the number of antennas and number of relays increase. However, as in the csi-assisted case we can observe that this improvement depends on the diversity gain. It is further observed that the increase of correlation decreases the performance. Moreover, one can notice that the csi-assisted relay outperforms the fixed gain one by approximately 3 dB. Moreover Monte Carlo simulation results exactly match with the analytical ones.

Average SER figures are depicted in Figures 3 and 4 for csi-assisted and fixed gain relay schemes, respectively. In both figures it is observed that the increase of number of antennas and the number of relays improve the average SER. Conversely, the increase of correlation parameters decreases it. Without loss of generality we have considered the BPSK, QPSK, and QAM schemes to demonstrate the average SER variation. High SNR curves are plotted and they are compatible with the exact ones in medium to high SNR and they show the diversity gain variation. Further, we can notice that the csi-assisted relay performs better than the fixed gain relay. Monte carlo simulations coincide with the analytical ones and it shows the accuracy of our results.

Figures 5 and 6 show the ergodic capacity variation for csi-assisted and fixed gain relay cases.
Without loss of generality, we fixed the number of antennas to be *n*_{s}=*n*_{d}=2 and *ρ*_{1}=*ρ*_{2}=2 dB. We have plotted the ergodic capacity variation against the number of relays
to demonstrate the fact that the ergodic capacity can be improved with the increase
of correlation for a higher number of relays. From the figures it is noticed that
the increase of correlation at the destination *ρ*_{d}, decreases the ergodic capacity , however, the increase of correlation parameter
at the source increases it. The reason for this behavior can be explained as follows;
the relay to destination is a point to point link, hence, the increase of correlation
decreases the ergodic capacity, however, the source to relay link is a point-to-multipoint
link, hence the increase of correlation parameter reduces the channel hardening effect
[26,32] which results in a higher capacity. Moreover, Monte Carlo simulations exactly coincide
with the analytical ones for csi-assisted one and are closely compatible with the
approximated fixed gain ergodic capacity.

### Conclusion

We have investigated the performance of a partial relay selection network with the effect of antenna correlation at the source and the destination. Two relay schemes; csi-assisted and fixed gain relay schemes have been considered and exact closed form expressions for outage probability, average SER and ergodic capacity have been derived. Our results can be used to quantify the effect of antenna correlation in partial relay selection. Further, we have provided an asymptotic analysis which can be used to obtain an insight of the system performance. In addition, we have showed that for a higher number of relays, the ergodic capacity can be improved with higher correlation at the source.

### Appendix 1

Let *γ*_{i}=*ρ*_{1}||**h**_{si}||_{F}and we rewrite (1) as,

We find the pdf of *γ*_{i}as [25,26],

and cdf of *γ*_{i}as,

We assume that the relays are distributed homogeneously such that they have equal
average SNR, further, we assume that the S−R_{i}∀*i* channels are independent. Then we can derive the cdf of *γ*_{1}as,

Using multinominal theorem and after some simplifications, we derive

### Appendix 2

#### Outage probability

Pdf of *γ*_{2} can be obtained from (6) and the cdf of *γ*_{1}is derived in (8). Following the same procedure as mentioned in ([4], Appendix A),

where

Performing the integration with the help of ([27], 3.471.9), we can obtain the closed form solution as in (12).

#### Ergodic capacity

We can rewrite (20) as,

After some mathematical manipulations [33],

Now, we can rewrite (73)

where *γ*_{3}=*γ*_{1} + *γ*_{2}. Now
*γ*_{2}as,

Performing the integration with the help of ([27], Equation 4.337.5) we obtain the closed form expression for *χ*_{2} as in (23). Similarly, we can derive
*γ*_{3}=*γ*_{1} + *γ*_{2}as,

We carry out the integration to obtain

where