Performance analysis of joint relay selection and transmission schemes in multiple-antenna two-way relay systems

In this paper, we consider a two-way relay system with two multi-antenna sources and multiple single-antenna relays and study the performances for various transmission schemes with the assumption that perfect channel state information (CSI) is available at both the sources and the relays. Although a transmission scheme that combines source beamforming and relay selection (BF-RS) is known to improve a system performance, its performance has not been theoretically analyzed yet. In the paper, we focus on the performance analysis in terms of the symbol error probability (SEP) for the BF-RS scheme. The analytical upper and lower bounds of the end-to-end system SEP are derived in the closed form, and the asymptotic SEP expression is given in high SNR regime. By both analytical and simulation results, it is clearly shown that a full transmit and relay selection diversity gain can be achieved without code rate loss. The SEP curves from simulations show that our derived bounds can predict the performance accurately. Based on the derived analytical results, we propose an individual power allocation (IPA) scheme to save the total transmit power. We show that the IPA scheme can result in considerable energy saving with a comparable SEP performance.


Introduction
In order to achieve a higher throughput and spectral efficiency, protocols consisting of the multiple access (MAC) and broadcast (BC) phases in a two-way relay channel (TWRC) have attracted more research interests than conventional protocols using three or four phases. Since the physical-layer network coding (PNC) protocol can be easily applied to two-way relay systems, its performance has been widely studied [1].
For reliable transmissions, it is important to exploit diversity gain in multi-input multi-output (MIMO) wireless communications. Relay selection (RS) diversity can also be exploited in MIMO TWRC in conjunction with transmit diversity. For instance, To et al. [2] considered a two-way relay system using the Alamouti's code at the two sources while there is only one antenna at a relay. Gong et al. [3] studied a system consisting of two singleantenna sources and a single relay of two-antenna. Relay and antenna selection schemes are also addressed in many related works in MIMO TWRC. Chen et al. [4] proposed a RS and beamforming scheme for cooperative two-way relay systems with the amplify-and-forward (AF) protocol. Jia et al. [5] investigated the outage probability and ergodic capacity of the two-way network coding opportunistic relaying systems with the AF protocol. Jayasinghe et al. [6] proposed a robust joint precoder-decoder design scheme for a multiple-input multiple-output PNC-based two-way relay system.
Zhou et al. [7] proposed an opportunistic two-way relay scheme based on joint network coding and opportunistic relaying with the decode-and-forward (DF) protocol. The analysis and design of space-time trellis codes (STTCs) for TWRC with multi-antenna sources and single-antenna relays were studied with the AF protocol [8].
Although various approaches for efficient transmission schemes with RS are proposed, it is interesting to note that joint transmission and RS schemes for MIMO TWRC are not well studied. In other words, the individual analysis of either the transmission schemes or the RS schemes has been well investigated, but a best combination of the transmission and RS schemes under a certain performance criterion is not known, especially in asymmetrical multi-antenna two-way relay systems. In addition, it is proved that an opportunistic RS scheme with low complexity can achieve better performance than a fully distributed two-way relay scheme [7] and it is shown that the DF protocol outperforms the AF protocol in terms of symbol error probability (SEP) performance [9]. With the demand of reliable transmissions, combining the best known techniques in multiple-antenna two-way relay systems would lead us to the best joint scheme; however, its performance is not known yet.
In addition, it is well known that the power allocation (PA) plays an important role in wireless communications. With a total power constraint, the PA is widely used to improve the system performance. With individual power constraints, the PA can be used to save the system power consumption. Song et al. [10] proposed a PA scheme based on the asymptotic SEP in AF two-way relay systems. Zhang et al. [11] proposed a PA scheme to minimize the maximum individual outage probability. Talwar et al. [12] presented an optimal joint RS and PA scheme to maximize the worst received signal-to-noise ratio (SNR) in AF two-way relay systems. Do et al. [9] proposed a joint relay selection and PA scheme to improve the SEP performance in a two-way relay system with the DF protocol. To the best of our knowledge, all these works are based on the total transmission power constraint. However, the relays and the sources could be battery-powered, which means that individual transmission power constraints at nodes are more suitable for two-way relay systems with battery-powered nodes. Hence, under individual transmission power constraints, it would be needed to consider an individual power allocation (IPA) scheme to improve the system power efficiency in a multiple-antenna twoway relay system.
In this paper, as mentioned earlier, the combination of the best known techniques for a multiple-antenna twoway relay system consisting of two sources of multiple antennas and multiple relays of single antenna is considered. A transmission scheme that combines beamforming and RS (BF-RS) is mainly considered. We present an asymptotic SEP expression for the proposed BF-RS scheme and make comparison with other transmission schemes. Through analysis and simulation results, it can be shown that a full transmit and RS diversity gain can be achieved by the BF-RS scheme. Thus, we can confirm that the BF-RS scheme is optimal in terms of the diversity gain.
Furthermore, we consider a general case where two sources are equipped with M and N antennas and K relays participate in the transmission to derive closedform expressions for tight upper and lower bounds on SEP.
In order to improve the system power efficiency, an IPA scheme is also proposed to reduce the total transmission power cost. The performance of the BF-RS scheme and the IPA scheme is verified through simulations.
The main contributions of this paper are summarized as follows. (i) Despite that the BF-RS scheme is based on existing approaches, its performance has not been well investigated in a theoretical manner. We derive the upper and lower bounds of the end-to-end system SEP in an asymmetrical multi-antenna two-way relay system, where the source nodes are equipped with different numbers of antennas. The derived SEP expressions are tight and are verified with the simulated results. (ii) We conduct the asymptotic SEP performance analysis and evaluate the diversity order. We demonstrate that the BF-RS scheme can achieve a full transmit and RS diversity gain and outperforms other transmission and RS schemes. (iii) Under individual transmission power constraints, we propose an IPA scheme to improve the system power efficiency. The behaviors of the system power efficiency under different system configurations are investigated. Meanwhile, both simulated and analytical results show that the proposed IPA scheme can yield considerable energy saving with a comparable SEP performance to other PA schemes.
The rest of this paper is organized as follows. In Section 2, a system model of a multiple-antenna two-way relay system is introduced. In Section 3, we investigate the optimal beamforming and the RS criteria. In Section 4, an asymptotic SEP expression and the closed-form bounds are discussed. In Section 5, the IPA scheme for the BF-RS is given. Simulation results are shown to corroborate our theoretical results in Section 6. In Section 7, we draw the main conclusion.
The following notations are used in this paper. Lower and upper boldfaced letters are used for column vectors and matrices. (·) * and (·) T represent the conjugate and transpose, respectively. For a vector or matrix, · denotes the Frobenius norm. CN (μ, σ 2 ) denotes the circularly symmetric complex Gaussian distribution with mean μ and variance σ 2 . For a random variable (RV) X, E[ f (X)] represents the expectation of a function of X, f (X). f X (·), and F X (·) refer to the probability density function (PDF), and the cumulative density function (CDF) of X, respectively.

System model
In this paper, we consider a multiple-antenna two-way relay system where two multi-antenna source nodes, S 1 and S 2 , are equipped with M and N antennas, respectively, and K relay nodes, {R 1 , . . . , R K } are equipped with single antenna. S 1 and S 2 exchange their information with the help of one selected relay, say R k . This system model can be possessed in practice. For instance, two base stations equipped with multiple antennas might communicate via a selective single antenna relay, which is due to the size and complexity constraints. We assume that the channel state information (CSI) remains unchanged during a transmission time unit consisting of two phases (the first is the MAC phase and the second is the BC phase). The detailed system model is shown in Fig. 1.
In the MAC phase, the relays can receive signals from S 1 and S 2 . Based on the perfect CSI assumption, the best relay R k out of K relays is selected (the RS criteria will be discussed later). Then, the selected relay forwards the detected signals to both S 1 and S 2 in the BC phase. Using the PNC protocol, S 1 and S 2 can detect the desired information from each other. We assume that all the sources and relays operate in time division duplex mode, and perfect synchronization has been established between the sources and the relays before data transmission. In this section, we will consider a generic representation for various transmission schemes that can achieve a full transmit diversity such as beamforming, space time block coding (STBC), and so on.
We assume that the same signal constellation C is used at sources and relays throughout the paper. The size of C is denoted by C. In each transmission time unit, S 1 and S 2 will exchange their frames of L symbols, which are denoted by x = [x 1 , x 2 , . . . , x L ] and y = [y 1 , y 2 , . . . , y L ], respectively. Let the operator (·) represent the transmission scheme and the transmission matrices X = (x) where w i , i = 1, 2 are the beamforming weights. With X and Y, the received signal at R k (k = 1, 2, . . . , K) in the MAC phase can be expressed as ,k ] T are the channel vectors between relay R k and S 1 and relay R k and S 2 , respectively. k ∼ CN 0, σ 2 I L is the background noise vector at R k . P m is the transmit power of S m , m = 1, 2.
Since the CSI is assumed to be perfectly known at R k , under the maximum likelihood (ML) criterion, the estimate of x and y, denoted by x' and y', is given by where (X, For simplicity, we consider a modification of (2) in the following derivations, which does not change the final SEP results. The modified detection rule at R k is given by x', y' = arg min Fig. 1 The system model of a two-way relay system with multi-antenna sources and single-antenna relays At each relay, an XORed-like version of (x j , y j ), j = 1, 2, . . . , L, is to be detected, which is given by where l(x) is considered as the mapping function [1] and l −1 (x) represents the inverse operation of l(x). The XORed-like version of (x , y ) can also be written as z = In time slot j of the BC phase, S m receives where P r is the transmit power at R k , g 1k = T are the channel vectors between relay R k and S 1 and relay R k and S 3 , Finally, each source detects the frame sent by the other source in a similar way as (4):x = f (ẑ 2 , y) andŷ = f ẑ 1 , x wherex andŷ are the estimates of x and y, respectively.

Transmission scheme and RS criterion
In this paper, the main objective is to investigate the best performance among various transmission schemes and RS criteria in multiple-antenna two-way relay systems. By analyzing existing transmission schemes and RS criteria, we find 1 that the Alamouti code with RS schemes can only achieve a diversity gain of d = min(M, N) at a code rate of r = 1. The diagonal code with RS schemes can achieve a full diversity gain of d = K min(M, N) at a code rate of r = 1 min(M,N) . In other words, the Alamouti code with RS and the diagonal code with RS schemes suffer performance loss in either the diversity gain or the code rate. Among the transmission schemes, beamforming schemes can be easily combined with RS and achieve the best performance. The scheme can minimize the overall system SEP and achieve a full diversity gain at a code rate of r = 1. Thus, we focus on the performance analysis for the BF-RS scheme in this section. Prior to studying the performance of the BF-RS scheme, we first derive an expression of the system SEP.
Let P k MAC denote the SEP in the MAC phase at R k and P k BC,m denote the SEP in the BC phase for given h 1k and h 2k at S m . Then, we have It is noteworthy that a correct detection occurs if both the phases have erroneous decisions in a two-way transmission with the BPSK modulation. Therefore, averaging the SEPs over the two phases, the instantaneous SEP at source node S m , denoted by P k EtoE,m , is given by Thus, the instantaneous end-to-end SEP of both the source nodes S 1 and S 2 , defined by P k EtoE = 1 2 2 m=1 P k EtoE,m , is given by Note that in the rest of the paper, we confine ourselves to BPSK for tractable analysis as we will focus on the performance analysis of the BF-RS scheme.

The optimal beamforming at the two sources
We redefine the transmission matrices at the two sources as T are the beamforming vectors chosen by source S 1 and S 2 . For simplicity, the time index subscript j is omitted in the following derivation. We assume that the transmission power of the sources and the relays are identical, i.e., P 1 = P 2 = P r = P. We also denote the SNR byγ = P σ 2 , where σ 2 is the Gaussian noise variance. We focus on the case of g mk = h mk , m = 1, 2, by considering that the links between S m and R k are reciprocal. In addition, we assume that the h 1i,k 's, i = 1, 2, . . . , M and h 2i,k 's, i = 1, 2, . . . , N, k = 1, 2, . . . , K, are independent and h mi,k ∼ CN (0, 1), i.e., Rayleigh fading channels. At relay R k , we can rewrite (3) as Recall that the signals x and y are transmitted from S 1 and S 2 , respectively. The estimate of x and y is written as x and y , respectively. Based on pairwise error probability (PEP) analysis and assuming that all code u are equally likely, the average error rate is upper bounded by N are the equivalent system SNR of S 1 and S 2 , respectively, and Q(·) is the Gaussian Q-function ( [15], 2-1-97).
With power constraints at the sources, we minimize P k MAC subject to that the norms of w 1 and w 2 are bounded by constant, i.e., respectively, which leads to w 1 2 ≤ M, w 2 2 ≤ N. The optimal beamforming vectors can be obtained by maximizing the SNR of S1 and S 2 , respectively, as the objective is the sum of two monotonous Q-function expressions, and the constraints are independent of each other. Thus, the optimal beamforming vectors in the multiple-antenna two-way relay system can be easily obtained as Note that the results are similar to the maximum ratio transmission (MRT) [16]. It is proved that the MRT obtained from one-way relay systems is still valid in the two-way relay system [17,18].
Using the optimal beamforming vectors in (12), the lower bound on P k MAC can be obtained exactly in the same way as in [19]. Assuming a genie-aided source knows the message from the other before transmission, we have where β k min = min h 1k 2 , h 2k 2 and β k m = h mk 2 .

RS criterion
The opportunistic two-way relay scheme with a low complexity in [7] achieves the better performance compared to a fully-distributed two-way relay scheme. Thus, we focus on the schemes that only one best relay is selected out of K relays to decode-and-forward the signals in the second phase transmission. For tractable analysis, we use the equivalent SNR as an effective metric. When relay R k is selected, we define the equivalent SNR for the MAC phase at relay R k and the BC phase at S m , m = 1, 2 as SNR MAC k ≤γ β k min and SNR BC , m k =γ β k m , respectively. Using the optimal beamforming, we compare two RS methods: one is an optimal beamforming and RS (O-BF-RS) scheme and the other is a suboptimal beamforming and RS (S-BF-RS) scheme.
1) O-BF-RS: For O-BF-RS, the destination node is to select the best relay among all the relays, which is denoted by R R . The selected relay R R can achieve the minimum overall system SEP in (9). That is, 2) S-BF-RS: The overall system performance is limited by the worst link. Thus, as a suboptimal scheme, we select the relay which minimizes the maximum SEP of all links. This can be done by selecting the relay that maximizes the minimum equivalent SNR of all links. Then, the RS criterion is given as Actually, the performance difference between the optimal scheme and the suboptimal scheme is negligible as the overall performance is generally decided by the worst link. This result can be verified in Section 6. However, even the performance is almost the same, the computation complexity of the two schemes are quite different. The general SEP expression which is conditioned on the instantaneous received SNR can be expressed as SEP(γ ) = Q( √ (cγ )) [20], where γ is the received SNR, c is a constant which depends on the modulation mapping and Q(·) is the Gaussian Q-function ( [15], 2-1-97). As a result, computing the Q-function to obtain the SEP expression and then determining the best relay for the O-BF-RS is much more complicated than that of the S-BF-RS, in which only the received SNR is needed. Thus, the S-BF-RS is more practical. Note that, for convenience, we use the BF-RS to replace the S-BF-RS in the rest of this paper.
In the following section, we present the performance analysis for the BF-RS scheme in the multiple-antenna two-way relay system.

The bounds on the overall system SEP
Based on the optimal beamforming design in (12) and the RS criterion in (15), we begin to analyze the SEP and the diversity gain of the BF-RS scheme. The SEPs of the MAC and BC phases are denoted by P R MAC and P R BC , respectively. The SEPs are given by and Obviously, the bounds on the overall system SEP in (9) It is well-known that for two positive values a 1 and a 2 of similar order, we have a 1 Q( Based on the derivation in [9], the upper and lower bounds on P R EtoE can be obtained as

The distributions of the equivalent system SNR
Corresponding to (19), we introduce the equivalent system SNR as γ min R =γ β R min with the RS scheme. The distribution function of γ min R is required to analyze the error performance of the two-way relay system. Based on (15), we can obtain the following results.

Average end-to-end SEP expression
For BPSK modulation, the average SEP can be expressed as [20] where γ is the instantaneous SNR. (21), the upper and lower bounds of the average end-to-end SEP, denoted by P R EtoE , are given as

Asymptotic SEP of the BF-RS scheme
In this subsection, we derive an asymptotic SEP expression of the BF-RS scheme.

Lemma 2.
The PDF of the equivalent system SNR γ min R can be approximated by polynomial terms as From (24) and the expression in (19), an expression for the asymptotic average SEP of the BF-RS transmitting scheme can be found in the following proposition.

Proposition 2.
In high SNR regime, the asymptotic average SEP of the overall system is Proof. See Appendix 7.2.
According to Proposition 2, we can conclude that the diversity gain of the joint beamforming and RS scheme is K min (M, N). The BF-RS scheme achieves a full transmit and RS diversity gain.

Transmit power allocation
We have derived an asymptotic expression for the SEP of the BF-RS scheme in the last section, which was based on the assumption that the transmission power of the sources and relays are identical, i.e., P 1 = P 2 = P r = P. The scheme is called as equal PA (EPA) scheme. For the case that the relays and sources are wire-powered, the total power constraint is usually adopted. In [9], an optimal power allocation scheme (OPA) have been proposed to achieve the better SEP performance. However, for the case that the relays and sources are battery-powered, we need to consider individual transmission power constraints, i.e., 0 ≤ P 1 , P 2 , P r ≤ P. Thus, in this section, we propose an individual PA (IPA) scheme under the individual transmission power constraints.
We also denote the SNRs at the relay and the sources by γ m = P m σ 2 , m = 1, 2, and γ r = P r σ 2 , respectively. Hence, we can rewrite (18) as According to the derivation of (19), we can find that the SEP performance is mainly limited by the worst link in the two-way relay system. As a result, reducing the power consumption of the best link partly will hardly affect the system performance. Based on this observation, we propose an IPA scheme to save power consumption. The proposed IPA scheme is given by where η R m = β R m = h mR , m = 1, 2 or equivalently, 27) and P r = P.
As in the derivation of (19), it is obvious to draw the conclusion that the IPA scheme have the same bounds as the EPA scheme, which implies that their performances would be almost the same (this is later confirmed by simulations in Section 6). Meanwhile, it can be shown that considerable transmission power saving is achieved using the IPA compared to the EPA.
We use the proposed IPA scheme in the multipleantenna two-way relay system with M antennas sources and K relays. The average power consumption for the IPA, normalized by the EPA, is given by ≤ 1, the total transmission power consumption of the IPA is always less than that of the EPA. Since a closed-form solution of the total power consumption for the proposed IPA scheme is very complicated, we consider the dual case when M = N, since the system is more efficient on this antenna configuration. For notational convenience, we use the pair (M, K) to represent the system configuration, e.g., (2,4) represents that both the source nodes have M = 2 antennas and a single relay is selected out of K = 4 relays. We investigate the behaviors of E in Lemma 3 and Proposition 3 to gain insights into the performance analysis.

Lemma 3.
We consider the BF-RS scheme with the IPA in the two-way relay system. When (M, K) = (1, 1), for two RVs X = η R 1 , and Y = η R 2 , we have the CDF of Z = (30) Proof. See Appendix 7.3.
From Lemma 3, the expectation of Z is obtained by  As shown in Section 6, a better power saving can be achieved by decreasing the numbers of the sources antennas and the relays. It indicates the IPA is more suitable to be applied in the scenarios where the numbers of sources antennas and relays are limited.
At the end of this section, we consider the complexity order of the proposed IPA schemes in the multi-antenna two-way relay system. Note that the complexity order is dominated by complex multiplications in this paper. For Eq. (27), the calculations include the 2 norm operations, additions, comparisons, and divisions. As a result, the complexity order of the IPA with M-antennas and Nantennas at S 1 and S 2 , respectively, is

Simulation results
In this section, we evaluate the performances of the BF-RS and the IPA schemes. BPSK modulation is used and the i.i.d. Rayleigh flat fading channels are assumed in all the simulations. We assume that the transmission power of the sources and the relays are identical, i.e., P 1 = P 2 = P r = 1(0 dB), in the simulations of Figs. 2, 3, 4, 5 and 6. For convenience of the comparison, we assume that the total power constraint of each PA scheme is P T = 3 (= 4.77 dB). Furthermore, we assume that the sources and relays have the same noise variance.
Besides, we introduce the curves obtained by using (SNR − a) −d as the references to analyze the diversity gain, where a is a real number and d is a positive integer. If the SEP curves of the schemes tend to be parallel with (SNR − a) −d in high SNR region, it reveals that a diversity gain of d is achieved.

Simulated results
In Fig. 2, we compare the performance of the O-BF-RS and S-BF-RS schemes, where P 1 = P 2 = P r = 1. It is shown that the performance difference between the S-BF-RS and the O-BF-RS schemes is negligible. The intuitive explanation of the phenomenon is that the overall system performance of the two way relay system is limited by the worst link. In other words, once the relay which maximizes the worst SNR of the two links is selected, it is equivalent to minimize the overall system SEP. Thus, The qualitative explanation of the phenomenon is that once the SNR is low, more antennas is beneficial to achieving the better performance. Thus, the performance of (3,2) outperforms that of (2,4) in low SNR regions. However, since the diversity of (2,4) is larger than the diversity of (3,2), the performance of (2,4) outperforms that of (3,2) in high SNR regions, which further results in the intersection of the curves. Figure 3 shows the overall system SEP performance of different transmission schemes to validate the advantages of the BF-RS scheme. We compare the BF-RS scheme with  (2,2) and (2,4), we can find that the RS is only helpful for coding gain, but not diversity gain in the AC-RS scheme. Overall, the BF-RS scheme outperforms the AC-RS and DC-RS schemes in all the SNR regions.

Analytical results
In Fig. 4, we illustrate the diversity gain of the BF-RS scheme with different configurations of (M, K). The performance curves of the BF-RS scheme show the same decay trend with that of SNR −KM for all the configurations. The results show that the BF-RS scheme achieves a diversity gain of d = KM. In addition, we can find that the performance of configurations of (M, K) = (3, 2) and (4,2) are better than those of configurations of (2, 3) and (2,4), respectively. From this, we can draw the conclusion that increasing the number of the source antennas is more helpful than increasing the number of the relays on the overall system SEP performance. Figure 5 shows the SEP performance of a special case of the BF-RS scheme with (M, K, N) = (1, 2, 2), (2, 2, 2), and (3,2,2). The simulation result is within the curves of the theoretical upper and lower bounds. As shown in Fig. 5, the upper bound is very close to the simulated SEP curve. The accuracy of upper bound can indicate that the approximations in the theoretical derivation of (19) are accurate. Meanwhile, it can be noticed that the lower bound is parallel with the simulated SEP curve with a constant gap in all SNR region. In summary, our derived bounds can predict the performance accurately. Furthermore, the analytical asymptotic SEP and simulated SEP of the BF-RS scheme are presented as well. It is shown that, at the high SNR region, the asymptotic SEP given by (25) approaches the simulated SEP. The result has verified our analytical derivation. Figure 6 shows the SEP performance of the BF-RS scheme with the proposed IPA scheme. The results of the EPA and OPA [9] are also provided for comparison purposes. It can be observed that the SEP curves of the proposed IPA scheme are very close to the SEP curves of the EPA scheme. Meanwhile, it is shown that the gap between the performances of OPA and IPA is about 1-2 dB. This gain results from the fact that the more power consumption is used in OPA than that in IPA.

PA
In Fig. 7, we compare the power consumption by the IPA and EPA schemes. The power cost of the EPA is the same as that of the OPA, i.e., P T = 3. It can be observed that the proposed IPA scheme obtains almost 14 % power saving when (M, K) = (1, 1). The power cost increases with the number of source antennas and the number of relays. The results verify our derivation. In summary, the EPA is introduced as a benchmark in Figs. 6  show that the OPA outperforms the EPA in terms of the SEP once the total power constraint is adopted. Likewise, the IPA outperforms the EPA in terms of the power cost once the independent power constraint is adopted. In Fig. 8, we plot the curved surface of the expectation of Z in terms of different numbers of source antennas and relays. It is shown that the E[ Z] increases monotonically, which confirms that the effective power saving of the proposed IPA scheme is obtained when the numbers of antennas of the sources and the relays are limited.

Conclusions
In this paper, we have studied the performance of various transmission schemes for multiple-antenna two-way relay systems with RS in terms of SEP. Specifically, we have focused on the performance of the BF-RS scheme as it can outperform other schemes significantly. We derived the analytical upper and lower bounds of the end-to-end system SEP in the closed form as well as the asymptotic SEP expression in high SNR regime. The analytical performance predictions with the derived SEP expressions were verified by simulations. Based on individual power constraints, we also proposed an IPA scheme to save the total transmit power. Both analytical and simulated results indicated that the BF-RS scheme with IPA could save considerable amount of energy. Endnote 1, 2, . . . , M, k = 1, 2, . . . , K, is independent chi-square distributed with 2M degrees of freedom, the PDF of γ 1k Likewise, the PDF and CDF of γ 2k can be easily obtained as well. Define β k min = min(γ 1k , γ 2k ), the CDF of β k min is F β k Using order statistics, the CDF of the equivalent system SNR γ min R = max k=1,2,...,Kγ Consequently, using the binomial expansion, we can obtain the CDF of γ min R . The PDF of γ min R can be obtained by differentiating the CDF in (31). Applying Lemma 1 in (19) and (21), the analytical endto-end SEP in Proposition 2 can be obtained by the integral of [22, 2.
According to the Taylor series expansion formula and the limit theory, we can approximate the PDF by polynomial terms as  (19), the SEP of the overall system is upper bounded by Applying Lemma 2 in [23], the average SEP of the overall system is given by where (·) is the Gamma function.

Proof of Lemma 3
Since X = η R 1 and Y = η R 2 are RVs with generalized Rayleigh distribution [ can be expressed as For x, y ≥ 0, we have Due to the symmetry, it is straightforward to obtain I 2 from I 1 . Then, the CDF F Z (z) = I 1 + I 2 = 2I 1 .

Proof of Proposition 3
Based on the RS criterion as (15), when K → ∞, it is obvious to see lim 2F(y) . The expectation of Z = min(X,Y ) max(X,Y ) , when K → ∞, can be expressed as Thus, we have Let us now consider the expectation when M → ∞ and K ∈ N + . Since X = η R 1 and Y = η R 2 are generalized Rayleigh RVs, X 2 and Y 2 are chi-square distributed with 2M degrees of freedom. Using high-order moments of the generalized Rayleigh distribution [15, 2-1-138] and Stirling's formula [24], we find and E X 2 = 2 (M+1) (M) = 2M. Therefore, the coefficient of variation for X is obtained as lim It completes the proof.