This article is part of the series MIMO Relays for Cooperative Wireless Networks.

Open Access Research

Impact of antenna correlation on the performance of partial relay selection

Nuwan S Ferdinand, Upul Jayasinghe, Nandana Rajatheva* and Matti Latva-aho

Author Affiliations

Department of Communications Engineering, , Oulu, Finland

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2012, 2012:261  doi:10.1186/1687-1499-2012-261

The electronic version of this article is the complete one and can be found online at: http://jwcn.eurasipjournals.com/content/2012/1/261


Received:13 April 2012
Accepted:12 July 2012
Published:16 August 2012

© 2012 Ferdinand et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 nsand ndantennas, 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 SRii 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,

| | h sm | | F = max 1 < i < L | | h si | | F (1)

where ||·||F denotes the Frobenius norm and hsi is the ns×1 channel vector between SR and the elements of hsi are modeled as mutually correlated Ralyeigh fading entries. Let ns×nscorrelation matrix at source be Φs, then Φ s = E [ h si h si ] , where 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 Rm and the received signal at Rmis given as,

y r = P s h sm w s x + v m (2)

where Ps is the transmitted power and hsm is given as in (1) and ws is the MRT weight vector which is defined as w s = h sm / | | h sm | | F . Additive white Gaussian noise (AWGN) component with Vmvariance at Rm is denoted as vm. Then Rm multiplies the received signal by gain G and transmits to the Dand the received signal at Dis given as,

y d = h md G P r ( P s h sm w s x + v m ) + v d (3)

where G is defined differently for the two relay schemes and is given in the next section. Pr is the transmitted power at Rm and 1×nd channel vector between RmDis hmdand its elements are mutually correlated such that the correlation matrix at D is Φ d = E [ h md h md ] . vd is noise vector at D and it elements are AWGN with Vdvariance. Now Dperforms MRC to obtain the signal as,

y d = w d h md G P r ( P s h sm w s x + v m ) + w d v d (4)

where w d = h md / | | h md | | F is the MRC weight vector. Now after some mathematical manipulations, we obtain the end-to-end SNR as,

γ e = P s V m | | h sm | | F 2 P r V d | | h md | | F 2 P r V d | | h md | | F + 1 G 2 V m (5)

Notation: Let ρ 1 = P s V m and ρ 2 = P r V d and define γ 1 = | | h sm | | F 2 ρ 1 and γ 2 = | | h md | | F 2 ρ 2 . Let the distinct eigenvalues of the correlation matrix at source Φsbe ϕ 1 , ϕ 2 , , ϕ n s and those of the correlation matrix at the destination Φdbe σ 1 , σ 2 , , σ n d .

Statistics of SNRs

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

p γ 2 ( z ) = u = 1 n d σ u n d 2 k = 1 , k u n d σ u σ k exp z ρ 2 σ u ρ 2 (6)

and the cumulative density function (cdf) of γ2can be derived using F γ 2 ( Θ ) = 0 Θ p γ 2 ( z ) dz with the help of ([27], Equation 2.321.2) as,

F γ 2 ( z ) = 1 u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) exp z ρ 2 σ u (7)

Now, we derive the cdf of γ1 as,

F γ 1 ( z ) = 1 w o = 1 L ( 1 ) w o + 1 L w o w 1 = 0 w o w 2 = 0 w 1 w n s 1 = 0 w n s 2 i = 1 , w n s = 0 n s × ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) w i 1 w i w i 1 w i × exp i = 1 , w n s = 0 n s ( w i 1 w i ) z ρ 1 ϕ i (8)

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,

G = 1 P s h sm F 2 + V m (9)

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

γ e = γ 1 γ 2 γ 1 + γ 2 + c γ 1 γ 2 γ 1 + γ 2 (10)

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 γedrops below a predefined threshold Θand it is mathematically represented as,

P out = Pr ( γ e < Θ ) = Pr γ 1 γ 2 γ 1 + γ 2 + c < Θ (11)

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

F γ e ( Θ ) = 1 2 w 0 = 1 L ( 1 ) w 0 + 1 L w 0 w 1 = 0 w 0 w 2 = 0 w 1 w n s 1 = 0 w n s 2 × i = 1 , w n s = 0 n s ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) w i 1 w i × w i 1 w i u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) × i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i ρ 2 σ u exp Θ ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) Θ ρ 1 ϕ i Θ ( Θ + c ) × K 1 2 i = 1 , w n s = 0 n s ( w i 1 w i ) Θ ( Θ + c ) ρ 1 ρ 2 ϕ i σ u (12)

where K1(·) 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,

F γ e ( Θ ) = 1 2 w 0 = 1 L ( 1 ) w 0 + 1 L w 0 w 1 = 0 w 0 w 2 = 0 w 1 w n s 1 = 0 w n s 2 × i = 1 , w n s = 0 n s ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) w i 1 w i × w i 1 w i u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) × i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i ρ 2 σ u Θ exp Θ ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) Θ ρ 1 ϕ i × K 1 2 Θ i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ρ 2 ϕ i σ u (13)

Average SER

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

P ser = E aQ ( 2 b γ e ) (14)

where E[·] is the expectation operator and Q(.) is the Gaussian Q function. The modulation schemes is defined by ab, mainly BPSK (a=1, b=1), M-ary PAM ( a = 2 ( M 1 ) M , b = 3 M 2 1 ) and M-PSK ( a = 2 , b = sin 2 Π M ) [28,29]. Carrying out integration by parts of (14), we obtain,

P ser = a b 2 Π 0 F γ e ( z ) z exp ( bz ) dz (15)

Now, we substitute F γ e ( Θ ) in (13) to (15) and perform the integration with the help of ([27], Equation 6.621.3), to obtain the closed form expression for the average SER of γeas,

P ser = a 2 4 a b w 0 = 1 L ( 1 ) w 0 + 1 L w 0 w 1 = 0 w 0 w 2 = 0 w 1 w n s 1 = 0 w n s 2 i = 1 , w n s = 0 n s ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) w i 1 w i w i 1 w i × u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i ρ 2 σ u Γ 5 2 Γ 1 2 1 ρ 2 σ u + i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i + b + β 1 5 2 × F 5 2 , 3 2 ; 2 ; 1 ρ 2 σ u + i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i + b β 1 1 ρ 2 σ u + i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ 1 + b + β 1 (16)

where

β 1 = 2 i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i ρ 2 σ u (17)

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 (Cerg) can be mathematically expressed as,

C erg = E [ R log 2 ( 1 + γ e ) ] (18)

Closed form expression for Cergcan be derived as,

C erg = R ln 2 χ 1 + χ 2 χ 3 (19)

where

χ 1 = w 0 = 1 L ( 1 ) w 0 + 1 L w 0 w 1 = 0 w 0 w 2 = 0 w 1 w n s 1 = 0 w n s 2 i = 1 , w n s = 0 n s × ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) w i 1 w i w i 1 w i e β 2 ρ 1 E 1 β 2 ρ 1 (20)

where E 1 ( z ) = z e x x dx is the exponential integral. χ2is given as,

χ 2 = u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) e 1 ρ 2 σ u E 1 1 ρ 2 σ u (21)

and

χ 3 = ξ 1 , 1 ρ 2 σ u β 2 ρ 1 ξ 2 , 1 ρ 2 σ u = β 2 ρ 1 (22)

ξ 1 = u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) w 0 = 1 L ( 1 ) w 0 + 1 L w 0 w 1 = 0 w 0 × w 2 = 0 w 1 w n s 1 = 0 w n s 2 1 ρ 2 σ u β 2 ρ 1 1 × i = 1 , w n s = 0 n s ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) w i 1 w i w i 1 w i × E 1 β 2 ρ 1 ρ 2 σ u e β 2 ρ 1 β 2 E 1 1 ρ 2 σ 2 ρ 1 e 1 ρ 2 σ u (23)

ξ 2 = u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) w 0 = 1 L ( 1 ) w 0 + 1 L w 0 w 1 = 0 w 0 × w 2 = 0 w 1 w n s 1 = 0 w n s 2 × i = 1 , w n s = 0 n s ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) w i 1 w i w i 1 w i × 1 + 1 1 ρ 2 σ u e 1 ρ 2 σ u E 1 1 ρ 2 σ u (24)

where

β 2 = i = 1 , w n s = 0 n s w i 1 w i ϕ i (25)

Proof

Appendix 2 □

High SNR analysis

Here we derive the high SNR expressions for the outage probability and the average SER. Let z = Θ ρ 1 and ρ2=μρ1.

High SNR outage probability

We rewrite the (7) by expanding the e z σ u using Maclaurin series as follows,

F γ 2 ( z ) = 1 u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) s = 0 z μ σ u s (26)

It is observed at high SNR the lower order terms ( z n ; n < n s ) sum to zero, hence, by collecting the higher order terms, we can express (28) in high SNR as,

F γ 2 ( z ) = u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) z μ σ u n d + o ( z n d + 1 ) (27)

We can further simplify (29) as,

F γ 2 ( z ) = z n d μ n d n d ! u = 1 n d σ u + o ( z n d + 1 ) (28)

Similarly, we can derive the F γ 1 ( z ) in high SNR as

F γ 1 ( z ) = z n s n s ! i = 1 n s ϕ i + o ( z n s + 1 ) L (29)

Simplification yields,

F γ 1 ( z ) = z n s L n s ! i = 1 n s ϕ i L + o ( z n s L + 1 ) (30)

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

F γ e ( z ) = ψ 1 z L n s + o ( z L n s + 1 ) , L n s < n d ψ 2 z n d + o ( z n d + 1 ) , L n s > n d ψ 3 z n + o ( z n + 1 ) , n s = L n d = n (31)

where,

ψ 1 = 1 n s ! i = 1 n s ϕ i L (32)

ψ 2 = 1 μ n d n d ! u = 1 n d σ u (33)

ψ 3 = ψ 1 + ψ 2 (34)

It is observed from (33) that the diversity gain Gd= minLnsnd. This shows that the performance of the system is dominated by one of the single links unless Lns=nd, hence, in order to fully utilize the resources, we like to keep Lns=nd. 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]

P ser = 2 t aψΓ t + 3 2 ( 2 b ρ 1 ) ( t + 1 ) Π + o ( ρ 1 ( t + 1 ) ) (35)

where a and b define the modulation scheme and ψis given as in (34)–(36). t=minLnsnd−1 and diversity gain Gd=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,

γ e = γ 1 γ 2 γ 2 + C (36)

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

G = E | | h sm | | F 2 1 P s | | h sm | | F 2 + V m (37)

then the constant Cgiven by

C = E γ 1 1 ( γ 1 + 1 ) 1 (38)

If the gain is selected as,

G = 1 P s E | | h sm | | F 2 | | h sm | | F 2 + V m (39)

the constant Cgiven by

C = E γ 1 ( γ 1 ) + 1 (40)

Outage probability

If outage probability is denoted by Pout, then

P out = P ( γ e < Θ ) = P γ 1 γ 2 γ 2 + C < Θ (41)

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

F γ e ( Θ ) = 1 2 w o = 1 L ( 1 ) w o + 1 L w o w 1 = 0 w o w 2 = 0 w 1 w n s 1 = 0 w n s 2 × v = 1 , w n s = 0 n s ϕ v n s 1 k = 1 , k v n s ( ϕ v ϕ k ) w v 1 w v × w v 1 w v u = 1 n d σ u n d 1 k = 1 , k i n d ( σ u σ k ) × i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i ρ 2 σ u ΘC × exp i = 1 , w n s = 0 n s ( w i 1 w i ) Θ ρ 1 ϕ i × K 1 2 i = 1 , w n s = 0 n S ( w i 1 w i ) ΘC ρ 1 ρ 2 ϕ i σ u (42)

where K1(.) is 1st order modified Bessel function of second type and C is the fixed gain.

Proof

Appendix 3 □

From Equation (8), cdf of γ1can be written as,

F γ 1 ( z ) = 1 α 4 exp z β 4 (43)

Where,

α 4 = w o = 1 L ( 1 ) w o + 1 L w o w 1 = 0 w o w 2 = 0 w 1 w n s 1 = 0 w n s 2 v = 1 , w n s = 0 n s × ϕ v n s 1 k = 1 , k v n s ( ϕ v ϕ k ) w v 1 w v w v 1 w v (44)

and

β 4 = exp i = 1 , w n s = 0 n s ( w i 1 w i ) z ρ 1 ϕ i (45)

Case 1

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

E [ γ 1 ] = α 4 β 4 (46)

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

C = α 4 β 4 + 1 (47)

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

Case 2

E 1 γ 1 + 1 = 0 1 γ 1 + 1 d F γ 1 dz d γ 1 = α 4 β 4 exp ( β 4 ) Ei ( β 4 ) (48)

Substituting to (40)

C = α 4 β 4 exp ( β 4 ) Ei ( β 4 ) 1 (49)

where Ei(·) is the exponential integral.

Average SER

Then substituting F γ e ( z ) from (44) to (15) and through the mathematical simplification with the help of ([27], Equation 6.614.5), the closed form expression for SER is obtained.

P ser = a 2 a b 4 w o = 1 L ( 1 ) w o + 1 L w o w 1 = 0 w o w 2 = 0 w 1 w n s 1 = 0 w n s 2 × v = 1 , w n s = 0 n s ϕ v n s 1 k = 1 , k v n s ( ϕ v ϕ k ) w v 1 w v w v 1 w v × u = 1 n d σ u n d 1 k = 1 , k i n d ( σ u σ k ) i = 1 , w n s = 0 n s ( w i 1 w i ) C ρ 1 ϕ i ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i + b 3 / 2 × exp 0 . 5 i = 1 , w n s = 0 n s ( w i 1 w i ) C ρ 1 ϕ i ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i + b × K 1 0 . 5 i = 1 , w n s = 0 n s ( w i 1 w i ) C ρ 1 ϕ i ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i + b K 0 0 . 5 i = 1 , w n s = 0 n s ( w i 1 w i ) C ρ 1 ϕ i ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i + b (50)

where K1(.) and K0(.) 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

E ( γ e h ) = h Γ ( h + 1 ) Γ ( h ) w o = 1 L ( 1 ) w o + 1 L w o w 1 = 0 w o w 2 = 0 w 1 × w n s 1 = 0 w n s 2 v = 1 , w n s = 0 n s ϕ v n s 1 k = 1 , k v n s ( ϕ v ϕ k ) w v 1 w v × w v 1 w v u = 1 n d σ u n d 1 k = 1 , k i n d ( σ u σ k ) × exp 0 . 5 i = 1 , w n s = 0 n s ( w i 1 w i ) C ρ 1 ϕ i ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i + b × i = 1 , w n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i h × W h , 0 . 5 0 . 5 i = 1 , w n s = 0 n s ( w i 1 w i ) C ρ 1 ϕ i ρ 2 σ u i = 1 , x n s = 0 n s ( w i 1 w i ) ρ 1 ϕ i (51)

where Wk,μ(.) 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,

C erg = 1 2 E γ e log 2 ( 1 + γ e ) (52)

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

C erg = 1 2 log 2 ( e ) ln ( 1 + E [ γ e ] ) E [ γ e 2 ] E [ γ e ] 2 2 ( 1 + E [ γ e ] ) 2 (53)

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,

F γ e ( z ) = 1 2 w o = 1 L ( 1 ) w o + 1 L w o w 1 = 0 w o w 2 = 0 w 1 w n s 1 = 0 w n s 2 × v = 1 , w n s = 0 n s ϕ v n s 1 k = 1 , k v n s ( ϕ v ϕ k ) w v 1 w v × w v 1 w v u = 1 n d σ u n d 1 k = 1 , k i n d ( σ u σ k ) × i = 1 , w n s = 0 n s ( w i 1 w i ) ϕ i σ u Dz μ × exp i = 1 , w n s = 0 n s ( w i 1 w i ) z ϕ i × K 1 2 i = 1 , w n s = 0 n s ( w i 1 w i ) Dz μ ϕ i σ u (54)

then (56) can be rewritten as,

F γ e ( z ) = 1 Ψ Dz μ exp β 5 z K 1 2 α 5 Dz μ (55)

where

Ψ = 2 w o = 1 L ( 1 ) w o + 1 L w o w 1 = 0 w o w 2 = 0 w 1 w n s 1 = 0 w n s 2 v = 1 , w n s = 0 n s × ϕ v n s 1 k = 1 , k v n s ( ϕ v ϕ k ) w v 1 w v w v 1 w v × u = 1 n d σ u n d 1 k = 1 , k i n d ( σ u σ k ) i = 1 , w n s = 0 n s ( w i 1 w i ) ϕ i σ u , (56)

β 5 = i = 1 , w n s = 0 n s ( w i 1 w i ) ϕ i (57)

and

α 5 = i = 1 , w n s = 0 n s ( w i 1 w i ) ϕ i σ u . (58)

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

F γ e ( z ) = 1 Ψ Dz μ S = 0 β 5 z S S ! 1 2 α 5 Dz μ 1 2 + k = 0 α 5 Dz μ 1 + 2 k 2 Υ k k ! ( k + 1 ) ! (59)

Now it is observed that z n , n < min ( L n s , n d ) sum to zero and after limiting to high order terms with some simplifications we have,

F γ e ( z ) = Ψ ( 1 ) N + 1 β 5 N 2 α 5 N ! k = 0 N 1 ( β 5 ) N k 1 α 5 1 + 2 k 2 ( N k 1 ) ! × D μ k + 1 Υ k k ! ( k + 1 ) ! z N + o z N + 1 (60)

where N=min(Lnsnd) and

Υ k = 1 2 ln α 5 Dz μ ψ ( k + 1 ) ψ ( k + 2 ) (61)

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],

P ser = 2 t aΛΓ t + 3 2 ( 2 b ρ 1 ) ( t + 1 ) Π + o ( ρ 1 ( t + 1 ) ) (62)

where

Λ = Ψ ( 1 ) N + 1 β 5 N 2 α 5 N ! k = 0 N 1 ( β 5 ) N k 1 α 5 1 + 2 k 2 ( N k 1 ) ! × D μ k + 1 Υ k k ! ( k + 1 ) ! (63)

whereΨ as in (58) and a and b define the modulation scheme. t=minLnsnd−1 and diversity gain Gd=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 Φsis ρ s | i j | and that of the correlation matrix Φdis ρ d | i j | . Without loss of generality we consider ρ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.

thumbnailFigure 1. Csi-assisted outage probability.

thumbnailFigure 2. Fixed Gain outage probability.

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.

thumbnailFigure 3. Symbol error rate, Csi-assisted systems.

thumbnailFigure 4. Symbol error rate, fixed gain systems.

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 ns=nd=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.

thumbnailFigure 5. Ergodic capacity, csi-assisted systems.

thumbnailFigure 6. Ergodic capacity, fixed gain systems.

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||hsi||Fand we rewrite (1) as,

γ 1 = ma x 1 < i < L γ i (64)

We find the pdf of γias [25,26],

p γ i ( z ) = i = 1 n s ϕ i n s 2 k = 1 , k i n s ϕ i ϕ k exp z ρ 1 ϕ i ρ 1 (65)

and cdf of γias,

F γ i ( z ) = 1 i = 1 n s ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) exp z ρ 1 ϕ i (66)

We assume that the relays are distributed homogeneously such that they have equal average SNR, further, we assume that the SRii channels are independent. Then we can derive the cdf of γ1as,

F γ 1 ( z ) = F γ i ( z ) L = 1 i = 1 n s ϕ i n s 1 k = 1 , k i n S ( ϕ i ϕ k ) × exp z ρ 1 ϕ i L (67)

Using multinominal theorem and after some simplifications, we derive F γ 1 ( z ) as in (8),

Appendix 2

Outage probability

Pdf of γ2 can be obtained from (6) and the cdf of γ1is derived in (8). Following the same procedure as mentioned in ([4], Appendix A), F γ e ( Θ ) can be expressed as,

F γ e ( Θ ) = 1 0 F ~ γ 1 Θ + Θ ( Θ + 1 ) z p γ 2 ( z + Θ ) dz (68)

where F ~ γ 1 ( z ) = 1 F γ 1 ( z ) . By substituting (8) and (6) to (70) and by mathematical simplifications we obtain the F γ e ( Θ ) as,

F γ e ( Θ ) = 1 w 0 = 1 L ( 1 ) w 0 + 1 L w 0 w 1 = 0 w 0 w 2 = 0 w 1 w n s 1 = 0 w n s 2 × i = 1 , w n s = 0 n s ϕ i n s 1 k = 1 , k i n s ( ϕ i ϕ k ) w i 1 w i w i 1 w i (69)

× u = 1 n d σ u n d 1 e Θ ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) Θ ρ 1 ϕ i k = 1 , k u n d ( σ u σ k ) ρ 2 σ u × 0 e z ρ 2 σ u i = 1 , w n s = 0 n s ( w i 1 w i ) Θ ( Θ + c ) z ρ 1 ϕ i dz (70)

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,

C erg = R ln 2 E ln 1 + γ 1 γ 2 γ 1 + γ 2 + 1 (71)

After some mathematical manipulations [33],

C erg = R ln 2 E ln ( 1 + γ 1 ) ( 1 + γ 2 ) 1 + γ 1 + γ 2 (72)

Now, we can rewrite (73)

C erg = R ln 2 E [ ln ( 1 + γ 1 ) ] + E [ ln ( 1 + γ 2 ) ] E [ ln ( 1 + γ 3 ) ] (73)

where γ3=γ1 + γ2. Now χ 2 = E [ ln ( 1 + γ 2 ) ] can be derived using pdf of γ2as,

χ 2 = 0 ( 1 + z ) u = 1 n d σ u n d 2 k = 1 , k u n d ( σ u σ k ) e z ρ 2 σ u ρ 2 dz (74)

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 χ 1 = E [ ln ( 1 + γ 1 ) ] as in (22). Performing convolution operation we obtain the pdf of γ3=γ1 + γ2as,

p γ 3 ( z ) = 0 z p γ 2 ( t ) p γ 1 ( z t ) dt (75)

We carry out the integration to obtain p γ 3 ( z ) as,

p γ 3 ( z ) = δ 1 ( z ) , 1 ρ 2 σ u β ρ 1 δ 2 ( z ) , 1 ρ 2 σ u = β ρ 1 (76)

where

δ 1 ( z ) = u = 1 n d σ u n d 1 k = 1 , k u n d ( σ u σ k ) w 0 = 1 L ( 1 ) w 0 + 1 L w 0