On the Performance of Decode-and-Forward Equal-Gain-Combining Relay Systems over Weibull Fading Channels

Relay-assisted wireless communications, where both the relay and the final destiny employ diversity-combining techniques, represent a compelling strategy for improving the signal-to-noise ratio (SNR) for mobile terminals, mainly at millimeter-wave (mmWave) frequency bands. In this sense, this work considers a wireless network that employs a dual-hop decode-and-forward (DF) relaying protocol, in which the receivers at the relay and at the base station (BS) use an antenna array. Moreover, it is considered that the received signals are combined at reception using equal-gain-combining (EGC). Recent works have enthusiastically employed the Weibull distribution so as to emulate the small-scale fading behavior in mmWave frequencies, which also motivates its use in the present work. For this scenario, exact and asymptotic expressions for the system’s outage probability (OP) and average bit error probability (ABEP) are derived in closed form. Useful insights are gained from these expressions. More precisely, they illustrate how the system and fading parameters affect the performance of the DF-EGC system. Monte Carlo simulations corroborate the accuracy and validity of the derived expressions. Furthermore, the mean achievable rate of the considered system is also evaluated via simulations. Useful insights regarding the system performance are obtained from these numerical results.


Introduction
Wireless communication systems are affected by diverse phenomena, among them, multipath fading produces a random attenuation over the received signals. This phenomenon can be emulated using different statistical distributions depending on the operating scenario. For instance, the Rayleigh and Rice distributions are typically considered for non-line-of-sight (NLOS) and line-of-sight (LOS) propagation scenarios, respectively [1]. However, these distributions do not adequately emulate the fading behavior at millimeterwave (mmWave) frequencies. As a consequence, other more general distributions, such as the Weibull distribution, have been considered for these scenarios [2][3][4]. In addition, it has been shown that the Weibull distribution is a versatile fading model due to its flexibility and excellent fit to experimental data [5,6]. Moreover, the Weibull fading model encompasses, as special cases, the Rayleigh, exponential and one-sided Gaussian distributions. Interestingly, the authors in [7] have shown that, unlike other distributions, the Weibull fading model allows emulating propagation scenarios worse (more severe) than those generated by Rayleigh fading, mainly for vehicle-to-vehicle (V2V) communications and wireless sensor network (WSN) scenarios. These applications are encompassed within the ultra-reliable low-latency communications (URLLC) and the massive machine-type communications (mMTC) use cases, respectively, which are established for fifth-generation (5G) and beyond 5G (B5G) mobile networks [8,9].
Several works have focused on enhancing the coverage and performance of wireless systems by counteracting the detrimental effects of fading. Diversity and combining techniques have been extensively studied in the literature [10][11][12][13][14]. Such techniques focus on processing various replicas of the transmitted signal, which are affected by independent fading [15,16]. The analysis of diversity-combining techniques involves the study of the statistics of the sum of random variables [17]. Thus, the probability density function (PDF) and the cumulative distribution function (CDF) of the sum of random variables are fundamental for evaluating the performance of wireless systems employing diversity techniques. In [18], the authors found an expression to approximate the PDF of the sum of independent and identically distributed (i.i.d.) Weibull random variables. More recently, exact expressions for the PDF and the CDF of the sum of i.i.d. Weibull random variables were derived in [19].
Among the different combination techniques proposed, maximal-ratio-combining (MRC) and equal-gain-combining (EGC) stand out [20][21][22]. In particular, MRC is the best combination technique, but it requires the estimation of the channel gains, that is, the receiver requires the fading amplitude for each diversity branch, and then phase compensation is performed on the signals to add them coherently. On the other hand, EGC has quite acceptable performance and lower implementation complexity than MRC since it only involves a coherent summation of the received signals [23]. To do so, the EGC receiver performs a phase detection and then a compensation stage. Conversely, the MRC receiver requires the same structures as EGC and, additionally, it employs pilot symbols to estimate the fading amplitude. By the above, EGC constitutes an interesting alternative considering the trade-off between performance and implementation complexity, especially for future mobile systems in which tens or hundreds of receiving antennas will be used.
In [24,25], the performance of MRC and EGC systems in generalized fading channels and non-identical Weibull fading channels was evaluated in terms of the outage probability (OP) [26], respectively. More specifically, in [24], the authors derived expressions to evaluate the OP in EGC systems considering Nakagami-m fading. Moreover, in [25], the authors indicated that there was no available analytical expression for analyzing the momentgenerating-function of the output signal-to-noise ratio (SNR) for EGC receivers operating over Weibull fading channels. As a consequence, they derived approximate expressions for the OP in this scenario. In [27], the authors proposed a multi-tag EGC scheme over Nakagami-m fading channels. In particular, this work analyzed a radio frequency identification (RFID) communication scenario. OP expressions were derived using the characteristic function method. Moreover, in [28], the OP for L-branch EGC diversity receivers was evaluated via simulations using the concept of importance sampling (IS) for Rayleigh and Rician fading channels. In [29], the OP and error rate performance of EGC over correlated Beaulieu-Xie fading channels were evaluated, where closed-form asymptotic expressions (lower and upper bounds) were derived for these performance indicators.
Several relaying techniques have been proposed in the literature to improve the performance of wireless communication systems. In [30], the authors developed lowcomplexity cooperative communication protocols that combat multipath-induced fading. Among them, we highlight the decode-and-forward (DF) protocol, which considers a relay that demodulates the signals and then modulates and re-transmits them to the final receiver (destination) [31]. This process increases the capacity and reliability of user terminals over bad or severe propagation conditions.
The DF protocol has been studied in different scenarios. In [32,33], the performance of wireless systems considering DF relays and Rayleigh fading channels was analyzed in terms of the OP and the average bit error probability (ABEP), respectively. In particular, in [32], an exact closed-form expression for the OP was derived assuming dissimilar fading parameters. Moreover, in [33], the authors derived ABEP expressions for binary differential phase-shift keying, as well as OP expressions for noise-limited systems considering Rayleigh fading channels. In [34], the performance of a DF relay system over Rayleigh fading channels, in which the receivers employ the selection combining technique, was analyzed. For this scenario, an exact closed-form expression to evaluate the OP was derived. In [31], a dualhop DF relaying system was considered in which the relay and destination nodes were subject to co-channel interference in the presence of Weibull fading. For this scenario, OP expressions based on moment estimators were derived, which were obtained as an infinite summation of Meijer's G-functions [35] ([Equation (5.3-1)]). In [36], the performance of dual-hop DF cognitive relay networks over independent non-identically distributed (i.n.i.d.) Weibull fading channels was evaluated in terms of the OP, the symbol error probability (SEP), and the ergodic capacity. Integral form expressions were derived for evaluating these key performance indicators (KPIs). In [37], a dual-hop DF system, in which the receivers employ a single antenna, was considered, where one hop was subject to η-µ fading and the other hop was subject to κ-µ fading. Here, SEP, OP, and outage capacity expressions were obtained in terms of infinite series. In [38], the OP of antenna selection schemes in dual-hop DF relay networks over Nakagami-m fading channels was investigated, considering that the source, relay, and destination were equipped with multiple antennas. Exact expressions in terms of nested summations were found to evaluate the OP. In [39], a power beacon-assisted cooperative decode-amplify-forward relaying network was analyzed, where the relaying transmission could be selected to remain silent or transmit information employing a DF or amplify-and-forward (AF) protocol based on the channel quality among the source, relay, and destination. For this scenario, closed-form expressions of the OP and the throughput over Rayleigh fading channels were derived. Recently, in [40], the authors considered non-orthogonal multiple access (NOMA) systems that employed dual-hop hybrid optical wireless and radio frequency communications with the aid of a DF relay. In particular, the optical wireless channel underwent exponential-generalized Gamma (EGG) fading, whereas the radio frequency channel underwent Rayleigh fading. Exact closed-form expressions of the OP were obtained in terms of the Mejier-G function, and approximated closed-form expressions of the ergodic capacity were derived in terms of the incomplete Gamma function [41] ([Equation (6.5.3)]). Finally, in [42], an exact closed-form expression for the SNR at the output of NOMA-EGC receivers was derived using a Laplace transform approach for a generalized κ − µ fading channel. Unfortunately, the derived expression was only valid for κ → 0, which is equivalent to Nakagami-m fading. Table 1 summarizes the key contributions of the aforementioned works. Based on the literature review, and to the best of the authors' knowledge, a wireless system that employs dual-hop DF relaying and EGC receivers has not previously been analyzed over Weibull fading channels. This is the primary motivation for this work, where a dual-hop cooperative network operates in which the relay node (R) assists the transmission from a source (S) to the destination (D) using the DF protocol. The relay operates in half-duplex mode, and S does not have a direct link to D due to the presence of deep shadowing. In addition, S and D employ EGC with a different number of receiving antennas. By considering this scenario, the main contributions of this paper are summarized as follows: • The PDF and the CDF of the SNR at the output of EGC receivers are calculated considering Weibull fading channels. • Novel exact and asymptotic closed-form expressions to evaluate the OP and the ABEP of dual-hop DF-EGC relaying receivers are also derived, which are easy to manipulate and do not rely on functions that are not built-in in common mathematical software. • Monte Carlo simulations in different scenarios show the accuracy of the derived expressions. Furthermore, the achievable rate of the considered system is also evaluated via simulations. Useful insights regarding the system performance are obtained from these numerical results. Cooperative communications X X X X X [P] = Proposal presented in this work.
The remaining sections are organized as follows: Section 2 presents some preliminary results related to the sum of Weibull random variables. Section 3 presents the system and channel models. Exact expressions for calculating the OP and the ABEP of DF-EGC systems are derived in Section 4, where an asymptotic analysis is also carried out. The numerical results are presented in Section 5. Finally, the main conclusions of this work are summarized in Section 6.
In the following, f X (X) and F X (x) denote the PDF and the CDF of the random variable X, respectively. In addition, E[·] denotes the expectation, P(·) represents the probability operator, R + denotes the set of positive real numbers, and the min(a, b) function returns the argument with the lowest value.

Preliminaries
The performance analysis of EGC receivers involves the sum of N fading envelopes, where N is the diversity order, or equivalently, the number of antennas in the array, i.e., where X n ∈ R + is a Weibull random variable with PDF given by [43] ([Equation (2)]) where k > 0 and λ > 0 represent the shape and scale parameters of the Weibull distribution, respectively. Assuming that {X n } N n=1 is a set of i.i.d. Weibull random variables, the PDF of Z can be obtained from [19] ([Equation (4)]) as follows: where Γ(·) is the gamma function [41] ([Equation (6.1.1)]) given by and the coefficients δ i are obtained recursively employing the following expressions [19] ([Equation (5)]) Although (3) has an infinite summation, according to [19], it is enough to use a suitable finite number of terms to guarantee quite accurate results.

System and Channel Models
Consider a dual-hop cooperative network where R assists the transmission from S to D using the DF protocol, as shown in Figure 1. The relay operates in half-duplex mode, and S does not have a direct link to D due to the presence of deep shadowing. In addition, the nodes R and D employ N r and N d antenna elements, respectively. The impinging received signals at nodes R and D undergo i.i.d. Weibull fading and they are combined using the EGC technique. Moreover, φ r and φ d denote the instantaneous SNR at the outputs of the EGC receivers at R and D, respectively. The diversity branches are affected by additive white Gaussian noise (AWGN). Thus, each noise sample can be modeled by a zero-mean complex Gaussian random variable with variance where N 0 is the unilateral noise power spectral density and T s is the symbol duration.

Performance Analysis
In this section, the performance of a DF-EGC relay system is evaluated in terms of the OP and the ABEP for binary modulations.
The methodology used for the analysis is described as follows: The PDF and the CDF of the SNR at the output of EGC systems are first calculated based on the preliminary results given in Section 2. Then, using the CDF expression and considering the operation scheme of DF relay systems, an exact expression to evaluate the OP of DF-EGC systems over Weibull fading is derived. Considering that the region of interest for the operation of wireless systems occurs in high SNR, an asymptotic analysis of the OP expression is performed. In the following, using the PDF-and CDF-derived expressions and using some preliminary results of the BEP for binary modulations over AWGN channels, an exact expression to calculate the ABEP of DF-EGC relay systems is derived. Then, an asymptotic closed-form analysis of this expression is also performed to provide more details about the system behavior. Some insights are obtained from all the derived expressions, which are also validated via Monte Carlo simulations in Section 5.

PDF and CDF of the SNR in EGC Systems
The instantaneous SNR at the output of an EGC receiver can be written as [26] ([Equation (9.51)]) where the factor 1/N normalizes the received power per antenna, and P is the received power per symbol.
From (1) and (8), the instantaneous SNR for an EGC system can be rewritten as it was employed (7), where the received energy per symbol is E s = PT s , and E s /N 0 is the normalized SNR per symbol. Since Z ≥ 0, (9) is a monotonic function, with the aid of [44] ([Equation (5.16)]) given by where W is the number of roots of the function y = g(x), X is a random variable, g(x) is a function of X and g (x) is the derivative of g(x), and, using (3), it can be shown that the PDF of Φ is given by (11) Figure 2 shows the PDF of the SNR for an EGC system, parameterized by different values of N and considering k = 1.5 in the Weibull fading. Figure 2a assumes λ = 0.5 and Figure 2b assumes λ = 1. This change in the λ parameter causes the PDF to have different behaviors. Thus, in the case of λ = 0.5, it can be observed that there is a greater probability that the random variable φ assumes small values. This is evidenced in the horizontal axis of the figures, where the horizontal axis of Figure 2a has values up to 10 and the horizontal axis of Figure 2b has values up to 50. In the figures, it should be noted also that the simulated PDF results agree with the theoretical PDF results in all the scenarios. In particular, the analytical PDFs are generated using 80 terms in the summation of (11).   From (11), the CDF of Φ can be obtained as

Outage Probability
The instantaneous SNR for dual-hop DF relay systems can be written as [31] where φ r and φ d were defined in Section 3.
Since φ r and φ d are i.i.d. random variables, the OP for dual-hop DF relay systems can be calculated as where we have used [44] ([Equation (6.81)]), given by where w, x, and y are random variables. In particular, x and y are independent random variables. Moreover, in (14), φ th is a threshold value for the SNR in order to guarantee an adequate system performance, and F Φ r (φ r ) and F Φ d (φ d ) are the CDF of the received SNR at R and D, respectively. By considering that R and D can employ a different number of receiving antennas, and as our system model considers EGC receivers at both R and D, then F Φ r (φ r ) and F Φ d (φ d ) are given by (12), and N is replaced by N r and N d , as appropriate. In addition, the subscripts r and d are used to differentiate the fading channel parameters related to the links S → R and R → D, respectively. Accordingly, the OP for the considered DF-EGC system can be written as where for x ∈ {r, d}. Now, an asymptote for the OP expression in the high SNR region is derived. It can be observed that, as E s /N 0 → ∞, the most dominant term in the infinite series of (16) appears for i = 0. Thus, by considering only this term, an asymptotic expression for the OP can be obtained as where the last term in the sum is negligible as E s /N 0 → ∞ because the exponent of the E s /N 0 is greater than that of the first two terms.
Moreover, the dominant term in (18) is the one in which the product k x N x is smallest, for x ∈ {r, d}. Therefore, the asymptote for the OP in the high E s /N 0 region can be written as OP where is named as the coding gain, and is the system diversity order.

Average Bit Error Probability
In (13), φ r and φ d are i.i.d. random variables; thus, from (15), it can be shown that the PDF of φ df , given by (13), is equal to Similarly to the OP analysis, F Φ r (φ r ) and F Φ d (φ d ) are given by (12), N is replaced by N r and N d , as appropriate, and the subscripts r and d are used to differentiate the fading channel parameters related to the links S → R and R → D, respectively. In addition, f Φ r (φ r ) and f Φ d (φ d ) are given by (11) and the same aspects as in the case of the CDF are considered. Therefore, (22) can be rewritten as where and for x ∈ {r, d}.
Then, the ABEP can be calculated as [1] where P b (φ df ) is the BEP in AWGN channels conditioned on the instantaneous value of φ df . For binary modulations, it is known that [26] ([Equation (9.2)]) where a depends on the modulation scheme: a = 1 is used for coherent binary-phase-shiftkeying (BPSK), a = 1/2 is employed for coherent orthogonal binary-frequency-shift-keying (BFSK), and a = 0.715 for coherent BFSK with minimum correlation. From (23)- (28), the ABEP is obtained as With the aid of [45] ([Equation (4.1.18)]), a change of variables, and after some algebraic manipulations, (29) can be rewritten as In (30), as E s /N 0 → ∞, the most dominant terms in the infinite series appear for i = 0 and j = 0. Thus, an asymptotic expression for the ABEP can be calculated considering only these terms; then, we obtain From (26), the terms in (31) with the product ζ r,0 ζ d,0 have the factor (E s /N 0 ) − 1 2 (k r N r +k d N d ) , which are negligible with respect to the terms with the factors (E s /N 0 ) − 1 2 k r N r or (E s /N 0 ) − 1 2 k d N d as E s /N 0 tends to infinity. Under this premise, we can take the most significative terms in (31) and obtaining Considering again (26), and as E s /N 0 → ∞, the most dominant term in (32) is the one in which the product k x N x is smallest for x ∈ {r, d}. Hence, after some algebraic manipulations, and using (24)- (26), the asymptote for the ABEP in the high E s /N 0 region can be written as where is the coding gain in the ABEP asymptote, where Υ x,i is given by (17), and ξ is the diversity gain, which is given by (21). Thus, the DF-EGC relay system diversity is not only determined by the receiver with the least number of antennas (R or D), but also by the Weibull distribution shape parameter k of each link, which plays an important role in the system performance.

Achievable Rate
The achievable rate in bps/Hz (bits per second/Hertz) of DF-EGC systems can be calculated considering Shannon's capacity expression for signal transmission over AWGN channels [46]; that is, where we have employed (13). The mean value of the achievable rate, i.e., E[C(φ df )], as well as the OP and the ABEP, are analyzed employing numerical results in the following section.

Numerical Results and Discussions
In this section, the performance of DF-EGC relay systems is analyzed in terms of the OP, the ABEP, and the mean achievable rate in different scenarios. For the OP and the ABEP, the exact analytical results are obtained employing 80 terms in the summations of (16) and (30), as appropriate. Monte Carlo simulations with 10 7 trials validate the analytical results. Figure 3 shows the OP as a function of the E s /N 0 ratio parameterized by the number of receiving antennas. In this first scenario, it is assumed that the number of antennas at R and at D are equal, i.e., N r = N d = N. Additionally, it is assumed that the channel conditions for the links S → R and R → D are also equal; therefore, k r = k d = 1.5 and, λ r = λ d = 1. In the figure, as E s /N 0 increases, the OP decreases, which is an expected result. In addition, when the number of antennas increases, the OP decreases since the additional number of antennas guarantees more diversity in the system. Moreover, it can be observed that the analytical results agree perfectly with the simulation results in all the scenarios. In particular, the exact analytical curve is obtained employing (16) and the asymptote is calculated via (19). Figure 4 shows the OP as a function of the E s /N 0 ratio parameterized by N r and N d . It is assumed that the link S → R is over different propagation conditions than the link R → D; thus, k r = 2.5 and k d = 2. In addition, it is considered that λ r = λ d = 1. In Figure 4a, where N r = 3, it can be observed that, when the number of antennas in D increases, the system performance improves. However, it can be observed that when D has more than 5 antennas, the OP does not decrease. From (19) and (21), the system performance is dictated by the link in which the product between the Weibull distribution parameter k and the number of antennas is smallest. Consequently, the link S → R dictates the performance in this scenario, and increasing the number of antennas in D does not reduce the OP. In counterpart, in Figure 3, it was noticed that increasing the number of antennas reduces the OP; however, in that scenario the number of antennas increased simultaneously in both R and D. Figure 4b is similar to Figure 4a, but N r = 5 is employed. When comparing both figures, it can be observed that lower OP values are obtained in Figure 4b. This occurs because the number of antennas in R increases from N r = 3 to N r = 5. Nevertheless, similarly to the previous scenario, the performance is dictated by the link S → R. Therefore, it can be observed that, from a number of antennas greater than N d = 7, the system performance does not improve.     Figure 5 shows the OP as a function of the E s /N 0 ratio parameterized by N d and N r . In this case, it is considered that k r = 2 and k d = 2.5. In Figure 5a, it can be observed that, when N d = 5 antennas, the performance improves slightly in the low SNR region when compared to the scenario N d = 3 antennas, but the OP does not improve in the high SNR region when compared to the same scenario. A similar behavior is observed in Figure 5b when N d = 7. This OP behavior was not observed in Figure 4 where k r > k d . Hence, these results evidence that the fading channel parameters (k and λ) have a significant impact on the DF-EGC system performance. Thus, increasing the number of receiving antennas at R or D does not necessarily guarantee that the OP is reduced based on the propagation conditions.    Figure 6 shows the OP as a function of N d , parameterized by N r , k r and k d , considering E s /N 0 = 10 dB and λ r = λ d = 1. It is interesting to observe that increasing the number of antennas in D only ensures a higher performance improvement when the number of antennas in R is also increased. Thus, it can be observed that, if the number of antennas in R is maintained fixed, and the number of antennas in D is increased, then the OP decreases up to a certain value, and then it no longer decreases even though N d increases. Therefore, the number of antennas in both EGC receivers must be properly selected in order to maintain an adequate trade-off between cost and performance. Finally, by comparing Figure 6a,b, it can be observed that, for the same number of antennas at R and D, the scenario of Figure 6a generates lower values of OP. For a better understanding of this result, consider the case where N r = 3 and N d = 4. In the scenario of Figure 6a, we have that N r k r = 7.5 and that N d k d = 8; as a consequence, from (21), the system diversity order is ξ = 3.75. On the other hand, in the scenario of Figure 6b, we have that N r k r = 6 and N d k d = 10, then ξ = 3. Hence, as the scenario of Figure 6a ensures a higher diversity order, it also ensures a lower OP.  (b) k r = 2, k d = 2.5. Figure 6. OP as a function of N d , parameterized by N r , k r and k d , considering E s /N 0 = 10 dB, and λ r = λ d = 1. Figure 7 shows the ABEP as a function of the E s /N 0 ratio, parameterized by N r considering N d = 5, k r = 1.5, k d = 0.9, λ r = 1.4, λ d = 1 and BPSK modulation, i.e, a = 1. Notice the accuracy of the exact and asymptotic theoretical results, obtained via (30) and (33), respectively, with the simulated results. In addition, it is observed that, as N r increases, the ABEP decreases. However, as the value of N r approaches the value of N d , the improvement in performance becomes smaller. In particular, when N r = N d = 5, the system performance is dictated by the link R → D, since the product k d N d is less than the product k r N r . Therefore, despite the value of N r being increased, the system diversity, or equivalently the system performance, is dictated by the product k d N d . Thus, increasing N r does not improve the ABEP. This can be verified analytically through (21) and (33), where it is determined that the system diversity order is equal to 2.25 when N d = 5. This produces a saturation in the ABEP despite the fact that N r increases.  Figure 8 shows the ABEP as a function of E s /N 0 , parameterized by N d considering N r = 5, k r = 1.5, k d = 0.9, λ r = 1.4, λ d = 1 and a = 1. In this case, the system performance for N r ≤ 5 is dictated by the link R → D since k d N d < k r N r . As a consequence, increasing N d up to 5 improves the system's performance; this system behavior will be maintained as long as k d N d < k r N r . Therefore, the results of Figures 7 and 8 indicate that increasing the number of antennas excessively in the relay or in the BS does not necessarily improve the system performance.  Figure 9 shows the ABEP as a function of E s /N 0 , parameterized by k r and a considering N r = 3, N d = 4, k d = 0.9, λ r = 1.4 and λ d = 1. Thus, this figure indicates the impact that the fading channel affecting the link S → R has on the system performance. In addition, a = 1 is used for BPSK and a = 1/2 is used for BFSK. It can be observed that, for both modulations, the theoretical results are in full agreement with the simulated results, and, as expected, BPSK has better performance than BFSK for the same value of k r . It can be observed that, as k r increases, the ABEP is reduced since the system diversity order increases. Interestingly, there is a greater improvement when k r increases from 0.5 to 1 than when it increases from 1 to 1.5. In the first case, the performance is dictated by the link S → R, given that k r N r < k d N d . On the other hand, when k r increases to 1.5, then the performance is now governed by the link R → D because k r N r > k d N d .

Mean Achievable Rate Results
Finally, in this section, the performance of DF-EGC systems is evaluated in terms of the mean achievable rate. Figure 10 shows the mean achievable rate as a function of the E s /N 0 ratio, parameterized by N r and N d considering k r = 1.5, k d = 0.9, λ r = 1.4, and λ d = 1. It can be observed that, as the E s /N 0 ratio increases, the achievable rate also increases, which is an expected result. Moreover, as N d increases, the mean achievable rate also increases. This occurs because the diversity at destiny increases. However, by comparing Figures 10a,b, it can be observed that the achievable rate increases more slowly when N r = 2. In fact, in this scenario, a saturation point of the achievable rate is observed; that is, from a certain value of N d , the achievable rate no longer increases (or increases very slowly) despite the increase in N d . Thus, similar to what was observed in the previous analysis, the performance of the DF-EGC system is given by the link in which the product k x N x is greater for x ∈ {r, d}. As a result, the saturation point of the rate when N d is increased is higher in the case of Figure 10b, when N r = 5. Figure 11 shows the mean achievable rate as a function of E s /N 0 , parameterized by λ d considering N r = N d = 4, k r = 1.5, k d = 0.9, and λ r = 1.4. Thus, this last figure shows the impact that the λ fading parameter has on the system performance. In particular, the parameter associated with the link R → D is considered. Although this parameter of the Weibull distribution is not directly related to the diversity of the system, it is observed that increasing λ d also increases the mean achievable rate. To explain this behavior, Equations (19) and (33) can be considered, where the factor Υ x,i appears, which depends on λ x , for x ∈ {r, d}. Therefore, an increase in λ d allows for improving the coding gain, which translates into better system performance. More specifically, this implies a lower OP, a lower ABEP, and, in the case of Figure 11, a higher achievable rate.  (b) N r = 5. Figure 10. Mean achievable rate as a function of E s /N 0 , parameterized by N r , and N d considering k r = 1.5, k d = 0.9, λ r = 1.4, and λ d = 1.  Figure 11. Mean achievable rate as a function of E s /N 0 , parameterized by λ d considering N r = N d = 4, k r = 1.5, k d = 0.9, and λ r = 1.4.

Conclusions
The performance of a DF-EGC dual-hop relay system was analyzed over the presence of Weibull fading channels, which is a scenario that has not been previously considered in the literature. In addition, unlike other previous works that considered Weibull fading in diverse scenarios and in which approximate or integral expressions were obtained to calculate some performance indicators, this work has presented exact and asymptotic closed-form expressions to evaluate the OP and the ABEP, which are a function of the fading parameters and the number of antennas at the relay and at the destiny. In addition, the derived expressions can be easily evaluated using widely available computing software, such as Matlab or Mathematica.
The analytical modeling showed that the system diversity is given by the link (S → R or R → D), where the product between the Weibull fading parameter k and the number of receiving antennas is smaller, which is a novel result. Thus, increasing the number of receiving antennas at R or D does not necessarily guarantee that the OP or the ABEP are reduced. For this reason, the expressions derived in this work are a helpful tool for designing DF-ECG systems and ensuring an adequate trade-off between cost and performance.
Finally, it is important to note that the ABEP for other modulation schemes, such as quadrature amplitude modulation (QAM), can be easily obtained based on the results presented in this work. In addition, an option for future research is the performance analysis of other techniques, such as by using the amplify-and-forward protocol over Weibull fading, or by considering other combination schemes, such as MRC.