This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

A Hybrid Energy Harvesting Protocol for Cooperative NOMA: Error Performance Approach

Faical Khennoufa, Khelil Abdellatif, Ferdi Kara,  Hakan Kaya, Xingwang Li,   Khaled Rabie,   and Halim Yanikomeroglu F. Khennoufa and A. Khelil are with Laboratory of Exploitation and Valorization of the Saharan Energetics Resources (LEVRES), Department of Electrical Engineering, Echahid Hamma Lakhdar University, El-Oued, Algeria, email:{khennoufa-faical, abdellatif-khelil}@univ-eloued.dz.F. Kara is with the Department of Computer Engineering, Zonguldak Bulent Ecevit University, Zonguldak, Turkey, 67100, e-mail: [email protected]. Kaya is with the Electical and Electronics Engineering, Zonguldak Bulent Ecevit University, Zonguldak, Turkey, 67100, e-mail: [email protected]. Li is with School of Physics and Electronic Information Engineering, Henan Polytechnic University, Jiaozuo, China, e-mail: [email protected]. Rabie is with the Department of Engineering, Manchester Metropolitan University, Manchester M1 5GD, U.K. e-mail: [email protected]. Yanikomeroglu is with the Department of Systems and Computer Engineering, Carleton University, Ottawa, K1S 5B6, ON, Canada, e-mail: [email protected].
Abstract

Cooperative non-orthogonal multiple access (CNOMA) has recently been adapted with energy harvesting (EH) to increase energy efficiency and extend the lifetime of energy-constrained wireless networks. This paper proposes a hybrid EH protocol-assisted CNOMA, which is a combination of the two main existing EH protocols (power splitting (PS) and time switching (TS)). The end-to-end bit error rate (BER) expressions of users in the proposed scheme are obtained over Nakagami-mm fading channels. The proposed hybrid EH (HEH) protocol is compared with the benchmark schemes (i.e., existing EH protocols and no EH). Based on the extensive simulations, we reveal that the analytical results match perfectly with simulations which proves the correctness of the derivations. Numerical results also show that the HEH-CNOMA outperforms the benchmarks significantly. In addition, we discuss the optimum value of EH factors to minimize the error probability in HEH-CNOMA and show that an optimum value can be obtained according to channel parameters.

Index Terms:
BER, cooperative, EH, hybrid, NOMA.

I Introduction

Non-orthogonal multiple access (NOMA) is widely renowned for increasing the number of users and improving spectral efficiency in next-generation wireless networks [1]. Besides, NOMA can be easily integrated with other physical layer schemes. Hence, the cooperative-NOMA (CNOMA) has been one of the most investigated topics, recently due to its ability to increase spectral efficiency, user fairness and expand network coverage [2]. On the other hand, to increase the energy efficiency and recharge the battery-limited devices, the energy harvesting (EH) in CNOMA has also received significant attention [3]. In literature, there are two major protocols in EH: power splitting (PS) and time switching (TS). In the PS protocol, the EH and information detection will be performed simultaneously with a PS factor, whereas in TS protocol, the EH and information detection will be implemented at different times slots [4].

In order to show the performance gain of EH-NOMA over NOMA without EH, the ergodic capacity and outage probability of simultaneous wireless information and power transfer (SWIPT) aided CNOMA (SWIPT-CNOMA) have been evaluated [5, 6, 7, 8, 9]. Through examining the previous works, although most of the existing EH-CNOMA papers discuss investigating outage probability and ergodic capacity, the error performance (one of the most important key performance indicators (KPIs)) investigation is very limited. In [10, 11, 12], the error performance of CNOMA with EH (with PS protocol) amplify-and-forward (AF) relay has been performed. The authors of [13] evaluated the error performance of the mmth best relay selection in EH (with PS protocol) AF multiple relay-assisted CNOMA. Then in [14], the authors investigate the error performance of CNOMA with a EH relay (with TS protocol). As seen above discussion, the related EH-CNOMA papers deal with only one EH harvesting protocol. However, a comparison of the existing EH protocols is important to select the proper protocol to guarantee the best performance in CNOMA. Regardless superiority of EH to schemes without EH, it is possible to improve EH to reach the best performance. In order to improve the energy efficiency of EH, an optimized protocol called hybrid EH that combines PS and TS jointly has been proposed and evaluated in terms of throughput [15]. According to the findings in [15], the hybrid protocol outperforms both TS and PS protocols. The hybrid protocol can become one of the important manners to improve EH performance in CNOMA. The hybrid protocol has been evaluated in terms of ergodic capacity and outage probability for CNOMA in [16, 17, 18].

To the best of the authors’ knowledge, most of the existing papers investigate the EH-CNOMA schemes by using either TS or PS protocol and no work provides hybrid EH for CNOMA in terms of bit error rate (BER). Nevertheless, the hybrid protocol can improve the performance of EH [15], and it is considered as a practical scheme, which is a combination of the existing EH protocols. In addition, the BER performance of the EH-CNOMA schemes has not been studied well by considering different EH protocol. Motivated by aforementioned discussions, in this paper, we investigate a hybrid EH protocol in CNOMA (HEH-CNOMA) in terms of BER. The main contributions of the paper are given as follow.

  • We propose HEH-CNOMA scheme where the relay applies a hybrid EH protocol to improve the performance of the existing EH schemes in CNOMA networks.

  • We also present three benchmark schemes (i.e., PS, TS and no EH) to evaluate the performance of the proposed HEH-CNOMA. We derive the exact end-to-end BER expressions for HEH-CNOMA and three benchmarks over Nakagami-mm fading channels. The simulations validate the perfect match of the analysis. Based on the extensive simulations, we reveal that the proposed HEH-CNOMA outperforms all three benchmarks remarkably.

  • We also discuss the optimum values for EH parameters and power allocation coefficient in HEH-CNOMA and we present that optimum values can be obtained to improve performance of both users.

The rest of the paper is presented as follows. The considered HEH-CNOMA system model is introduced along with three benchmark schemes. In Section III we analyze the closed-forms BER of all schemes. The numerical results are presented in Section IV to validate analysis and compare with the benchmarks. Finally, Section V concludes the paper.

II System Model

We consider a CNOMA network consisting of a source (S), a decode-and-forward (DF) relay (R), and two users named user 1 (U1) and user 2 (U2). We assume that the R assists the users to receive the signal since the direct link between source and users does not exist due to large path-loss or obstacles. Each node is equipped with a single antenna. We suppose that the R node works in half-duplex (HD) mode. All communication links are exposed to a Nakagami-mm fading. The R node harvests energy from the radio frequency (RF) signal transmitted by the source (i.e., using three EH protocols defined below). In the second phase, the R node forwards the signal to two users using the harvested energy. The transmission block period is TT. The EH protocols, namely i) PS; ii) TS; and iii) hybrid (i.e., proposed) and their operations are described below.

  • PS: In the first T/2T/2 time slot, the received signal’s power is split into two parts at the R: The first part (ρ\rho) is used to harvest energy, while the second part (1-ρ\rho) is utilized to decode data in the decoder (where ρ\rho is the PS factor). The remaining time (i.e., T/2T/2) is devoted to delivering the data to the users.

  • TS: In the first βT\beta T time slot, the R harvests energy, where 0β10\leq\beta\leq 1 and β\beta is the TS factor. Then, in the first (1β)T/2(1-\beta)T\textfractionsolidus 2 of the remaining time, the S sends a superimposed coding (SC) information to the R. Then, in the last (1β)T/2(1-\beta)T\textfractionsolidus 2 time slot, the R uses the harvested energy to deliver this information signal to the users.

  • Hybrid (proposed): In the first βT\beta T time slot, the R harvests energy as being in TS mode, where 0β10\leq\beta\leq 1. Then, the S transmits an SC information to the R in the first (1β)T/2(1-\beta)T\textfractionsolidus 2 of the remaining time. Also, during this first (1β)T/2(1-\beta)T\textfractionsolidus 2 time period, the R harvest energy from the received information, where 0 ρ1\leq\rho\leq 1 (ρ\rho is the PS factor for energy harvesting during the first (1β)T/2(1-\beta)T\textfractionsolidus 2 at the R). Subsequently, the R uses the harvested energy from βT\beta T plus (1β)T/2(1-\beta)T\textfractionsolidus 2 time slots to forward the information signal to the users in the last (1β)T/2(1-\beta)T\textfractionsolidus 2.

  • No EH: The R is not capable of harvesting energy, and it uses its own energy source (e.g., battery). The S and the R cover half of the communication (T/2 time), and both consume equal power during these data transfers.

At the S, the messages of users are combined in an SC information signal with different power levels according to the users’ channel gain111In the NOMA system, the S or R combines the users’ messages in an SC signal with different power levels according to the channel gains for each user, where the user with higher channel gain will mostly have less power. Thus, we take into account that the U2 has a higher channel gain than U1, i.e., |h1|2<|h2|2|h_{1}|^{2}<|h_{2}|^{2} and the R receives this information signal in the first phase. Thus, the received information signal by the R node is given as

yr=Psϖ(α1s1+α2s2)hr+n,y_{r}=\sqrt{P_{s}\varpi}(\sqrt{\alpha_{1}}s_{1}+\sqrt{\alpha_{2}}s_{2})h_{r}+n, (1)

where ϖ\varpi is the coefficient that indicates how much of the source power is dedicated to information transfer, it will be (1ρ)(1-\rho) for PS and hybrid mode whereas equals to 1 for TS and no-EH modes. s1s_{1} and s2s_{2} are the U1 and U2 signals, respectively. α1\alpha_{1} and α2\alpha_{2} are the power allocation coefficients of U1 and U2, respectively in which α1>α2\alpha_{1}>\alpha_{2}, such that α1+α2=1\alpha_{1}+\alpha_{2}=1. PsP_{s} is the transmit power. hrh_{r} is the Nakagami-mm fading channel coefficient222The envelope of hrh_{r} and hkh_{k} follows Nakagami-mm fading with Ωr\Omega_{r} and Ωk\Omega_{k} spread and mrm_{r} and mkm_{k} shape parameters. between S-R. n𝒞𝒩(0,σ2)n\sim\mathcal{CN}(0,\sigma^{2}) is the additive white Gaussian noise (AWGN) at R. When the information signal is received, the RR node decodes the s1s_{1} firstly. Then, through the SIC, s2s_{2} is recovered.

For fairness, we consider the total consumed energy in all protocols is equal to ϵT\epsilon_{T}. Thus, ϵT\epsilon_{T} is consumed equally by the S and R in the no-EH mode, while it is consumed only by the S in EH protocols, the source power is expressed as

Ps=ϕPT,wherePT=ϵT/T,P_{s}=\phi P_{T},\ \text{where}\ P_{T}=\epsilon_{T}/T, (2)

where ϕ\phi is equal to 1 for no-EH, 2 for PS and 2β+1\frac{2}{\beta+1} for TS and hybrid. In EH protocols, the harvested energy during the power transfer in each protocol is given as

Eh={ηPsρT2|hr|2,PS,ηPsβT|hr|2,TS,ηPs(2β+ρ(1β)T)2)|hr|2,Hybrid,E_{h}=\left\{\begin{array}[]{lll}\eta P_{s}\rho\frac{T}{2}|h_{r}|^{2},&\mbox{PS},\\ \eta P_{s}\beta T|h_{r}|^{2},&\mbox{TS},\\ \eta P_{s}(\frac{2\beta+\rho(1-\beta)T)}{2})|h_{r}|^{2},&\mbox{Hybrid},\\ \end{array}\right. (3)

where η\eta (0<η<10<\eta<1) is the energy conversion efficiency factor. Thus, the transmit power of the relay through the harvested energy is given by

Pr={Ps|hr|2Ψ,PS, TS and Hybrid,Ps,no EH,P_{r}=\left\{\begin{array}[]{ll}P_{s}|h_{r}|^{2}\Psi,&\mbox{PS, TS and Hybrid},\\ P_{s},&\mbox{no EH},\end{array}\right. (4)

where

Ψ={ηρ,PS,2ηβ/(1β),TS,η(ρ+2β/(1β)),Hybrid.\Psi=\left\{\begin{array}[]{lll}\eta\rho,&\mbox{PS},\\ \nicefrac{{2\eta\beta}}{{(1-\beta)}},&\mbox{TS},\\ \eta(\rho+\nicefrac{{2\beta}}{{(1-\beta)}}),&\mbox{Hybrid}.\end{array}\right. (5)

Based on the received information signal, the R implements maximum likelihood detection (MLD)/successive interference cancellation (SIC), to detect s1s_{1} and s2s_{2} symbols respectively.Based on the signals decoded in the previous phase, the R implements a new SC signal and by using the harvested energy, it forwards this new SC signal to the users. The received information at the users is given by

yk=Prϖ(α1s1+α2s2)hk+n,k=1,2,y_{k}=\sqrt{P_{r}\varpi}(\sqrt{\alpha_{1}}s_{1}+\sqrt{\alpha_{2}}s_{2})h_{k}+n,\ k={1,2}, (6)

where hkh_{k}, k=1,2k={1,2} is the Nakagami-mm fading channel coefficient between RR-Uk. Through the information signal received by users, the U1 detects its symbols s1s_{1} directly, while the U2 use the SIC to detect its symbols s2s_{2}.

III Performance analysis

In this section, over Nakagami-mm fading channels, we derive the closed-form BER of the two phases. Then, we obtain the end-to-end (e2e) BER expressions for each user. At both S and R, a binary phase-shift keying (BPSK) modulation is utilized to obtain the s1s_{1} and s2s_{2} signals. Hence, two different energy levels with similar probabilities could be presented. By using a MLD, for all EH protocols and no EH, the R detects s1s_{1} firstly by treating s2s_{2} as noise. Thereafter, s2s_{2} is detected using the SIC. Then, the R implements an SC signal again for detected s1s_{1} and s2s_{2} symbols and forwards it to users by using the harvested energy (i.e., for the case of the EH protocols) or its own power (i.e., case of no EH). Through MLD, U1 detects its own symbols s1s_{1}, whereas U2 uses a SIC to identify its symbols s2s_{2}. By using the common probability of independent error event in two phases, the e2e BER at each user is given as in [19] by

Pe2e,k(e)=1(1Pr,k(e))(1Pk(e)),k=1,2,P_{e2e,k}(e)=1-\left.(1-P_{r,k}(e)\right.)\left.(1-P_{k}(e)\right.),\ k={1,2}, (7)

where Pr,k(e)P_{r,k}(e) is the BER of Uk at the relay. Pk(e)P_{k}(e) is the BER at Uk in the second phase.

III-A BER in the first phase

Based on the discussion above (i.e., for EH and no EH), R detects s1s_{1} in the first step. The BER of s1s_{1} at R, by following the steps in [20], can be given as

Pr,1(e)=12i=12Q(ϖγζi|hr|2),P_{r,1}(e)=\frac{1}{2}\sum_{{i}=1}^{2}\mathrm{Q}(\sqrt{\varpi\gamma\zeta_{i}|h_{r}|^{2}}), (8)

where γ=Psσ2\gamma=\frac{P_{s}}{\sigma^{2}}, ζi=[(α1+α2)2,(α1α2)2]\zeta_{i}=[(\sqrt{\alpha_{1}}+\sqrt{\alpha_{2}})^{2},(\sqrt{\alpha_{1}}-\sqrt{\alpha_{2}})^{2}]. The average BER (ABER) of (8) over a Nakagami-mm fading channel by using [21] is given as

Pr,1(e)=12i=1212πΓ(mr+0.5)Γ(mr+1)×(ϖγζiδ/2mr(1+ϖγζiδ/2mr)(mr+0.5))1(mr,ϖγζiδ),\begin{split}&P_{r,1}(e)=\frac{1}{2}\sum_{{i}=1}^{2}\frac{1}{\sqrt{2\pi}}\frac{\Gamma({m_{r}+0.5})}{\Gamma(m_{r}+1)}\\ &\times(\frac{\sqrt{\varpi\gamma\zeta_{i}\delta/2m_{r}}}{(1+\sqrt{\varpi\gamma\zeta_{i}\delta/2m_{r}})^{(m_{r}+0.5)}})\nabla_{1}(m_{r},\varpi\gamma\zeta_{i}\delta),\\ \end{split} (9)

where 1(mr,ϖγζiδ)=2F1(1,mr+0.5,mr+1,2mr2mr+ϖγζiδ)\nabla_{1}(m_{r},\varpi\gamma\zeta_{i}\delta)=_{2}F_{1}(1,m_{r}+0.5,m_{r}+1,\frac{2m_{r}}{2m_{r}+\varpi\gamma\zeta_{i}\delta}) and δ=Ωr\delta=\Omega_{r}. F12(.,.;.,.){}_{2}F_{1}(.,.;.,.) is the Gauss hypergeometric function.

We also need to analyze the BER performance of the correct and erroneous SIC as [20] to find the BER of s2s_{2} at R. Thus, the BER performance of s2s_{2} at R is given by

Pr,2(e)=12j=15νjQ(ϖγζj|hr|2),P_{r,2}(e)=\frac{1}{2}\sum_{{j}=1}^{5}\nu^{j}\mathrm{Q}(\sqrt{\varpi\gamma\zeta_{j}|h_{r}|^{2}}), (10)

where νj=[2,1,1,1,1]\nu^{j}=[2,-1,1,1,-1], ζj=[α2,(α1+α2)2,(2α1+α2)2,(α1α2)2,(2α1α2)2]\zeta_{j}=[\alpha_{2},(\sqrt{\alpha_{1}}+\sqrt{\alpha_{2}})^{2},(\sqrt{2\alpha_{1}}+\sqrt{\alpha_{2}})^{2},(\sqrt{\alpha_{1}}-\sqrt{\alpha_{2}})^{2},(2\sqrt{\alpha_{1}}-\sqrt{\alpha_{2}})^{2}]. By using [21], the ABER of (10) over the Nakagami-mm fading channel can be given as

Pr,2(e)=12j=15νj12πΓ(mr+0.5)Γ(mr+1)×(ϖγζjδ/2mr(1+ϖγζjδ/2mr)(mr+0.5))2(mr,ϖγζjδ),\begin{split}&P_{r,2}(e)=\frac{1}{2}\sum_{{j}=1}^{5}\nu^{j}\frac{1}{\sqrt{2\pi}}\frac{\Gamma({m_{r}+0.5})}{\Gamma(m_{r}+1)}\times\\ &(\frac{\sqrt{\varpi\gamma\zeta_{j}\delta/2m_{r}}}{(1+\sqrt{\varpi\gamma\zeta_{j}\delta/2m_{r}})^{(m_{r}+0.5)}})\nabla_{2}(m_{r},\varpi\gamma\zeta_{j}\delta),\\ \end{split} (11)

where 2(mr,ϖγζjδ)=2F1(1,mr+0.5,mr+1,2mr2mr+ϖγζjδ)\nabla_{2}(m_{r},\varpi\gamma\zeta_{j}\delta)=_{2}F_{1}(1,m_{r}+0.5,m_{r}+1,\frac{2m_{r}}{2m_{r}+\varpi\gamma\zeta_{j}\delta}).

III-B BER in the second phase

Using the harvested energy (i.e., using the power harvested of the three EH protocols as presented in section II), the R node forwards the SC signal to the users. The U1 detect its symbols directly using MLD, so the BER at U1 can be given by

P1(e)=12i=12Q(γζiψ|h1|2|hr|2).P_{1}(e)=\frac{1}{2}\sum_{{i}=1}^{2}\mathrm{Q}(\sqrt{\gamma\zeta_{i}\psi|h_{1}|^{2}|h_{r}|^{2}}). (12)

As seen in (12), the conditional BER includes multiplying of two random variables which belong to first and second phase link qualities. Although, the conditional BER in (12) defines the BER of second phase, it includes |hr|2|h_{r}|^{2} since the transmit power of the relay (i.e., PrP_{r}) is a function of it (4). Thus, the ABER of (12) can not be simply obtained by averaging Gamma distribution (i.e., the absolute square of Nakagami-mm). Firstly, the joint probability density function (PDF) of |h1|2|hr|2|h_{1}|^{2}|h_{r}|^{2} should be derived. By following steps in [22] through the PDF of the product of two Chi-Square distributions, the PDF can be obtained. After simplification and using [23, eq. (3.478.4)], we derive the ABER of (12) by using the Meijer-G function [23, eq.(9.311)] by following the steps of [22, 24], and the closed-form ABER at U1 is given as

P1(e)=12i=12Δ12πG3345(0,1u2,0.5u2,1u20.5u1,0.5u1,1u2,u2,0|Δ34),P_{1}(e)=\frac{1}{2}\sum_{{i}=1}^{2}\frac{\Delta_{1}}{2\sqrt{\pi}}G\begin{smallmatrix}3&3\\ 4&5\end{smallmatrix}\left(\begin{smallmatrix}0,1-u_{2},0.5-u_{2},1-u_{2}\\ 0.5u_{1},-0.5u_{1},1-u_{2},u_{2},0\end{smallmatrix}\middle|\frac{\Delta_{3}}{4}\right), (13)

where u1=mrm1,u2=0.5(mr+m1),Δ1=Δ2(m1ΩrmrΩ1)u12(1γζiψ)u2,Δ2=(mrΩr)mr(m1Ω1)m1Γ(mr)Γ(m1),Δ3=2m1mrΩ1Ωrγζiψu_{1}=m_{r}-m_{1},u_{2}=0.5(m_{r}+m_{1}),\Delta_{1}=\Delta_{2}(\frac{m_{1}\Omega_{r}}{m_{r}\Omega_{1}})^{\frac{u_{1}}{2}}(\frac{1}{\gamma\zeta_{i}\psi})^{u_{2}},\Delta_{2}=\frac{(\frac{m{r}}{\Omega_{r}})^{m{r}}(\frac{m{1}}{\Omega_{1}})^{m{1}}}{\Gamma(m_{r})\Gamma(m_{1})},\Delta_{3}=2\sqrt{\frac{m_{1}m_{r}}{\Omega_{1}\Omega_{r}\gamma\zeta_{i}\psi}}, mrm_{r} and mkm_{k} are the Nakagami fading parameters.

Also, we obtain the BER at U2 through the correct and erroneous SIC. Thus, the BER at U2 is given by

P2(e)=12j=15νjQ(γζjψ|h2|2|hr|2).P_{2}(e)=\frac{1}{2}\sum_{{j}=1}^{5}\nu^{j}\mathrm{Q}(\sqrt{\gamma\zeta_{j}\psi|h_{2}|^{2}|h_{r}|^{2}}). (14)

As discussed above, the ABER of (14) needs further modification. Thus, we find the ABER of (14) by using the same method of (13) and it is given as

P2(e)=12j=15νjΔ42πG3345(0,1u4,0.5u4,1u40.5u3,0.5u3,1u4,u4,0|Δ64),P_{2}(e)=\frac{1}{2}\sum_{{j}=1}^{5}\nu^{j}\frac{\Delta_{4}}{2\sqrt{\pi}}G\begin{smallmatrix}3&3\\ 4&5\end{smallmatrix}\left(\begin{smallmatrix}0,1-u_{4},0.5-u_{4},1-u_{4}\\ 0.5u_{3},-0.5u_{3},1-u_{4},u_{4},0\end{smallmatrix}\middle|\frac{\Delta_{6}}{4}\right), (15)

where u3=mrm2,u4=0.5(mr+m2),Δ4=Δ5(m2ΩrmrΩ2)u12(1γζjψ)u4,Δ5=(mrΩr)mr(m2Ω2)m2Γ(mr)Γ(m2)u_{3}=m_{r}-m_{2},u_{4}=0.5(m_{r}+m_{2}),\Delta_{4}=\Delta_{5}(\frac{m_{2}\Omega_{r}}{m_{r}\Omega_{2}})^{\frac{u_{1}}{2}}(\frac{1}{\gamma\zeta_{j}\psi})^{u_{4}},\Delta_{5}=\frac{(\frac{m{r}}{\Omega_{r}})^{m{r}}(\frac{m{2}}{\Omega_{2}})^{m_{2}}}{\Gamma(m_{r})\Gamma(m_{2})} and Δ6=2m2mrΩ2Ωrγζjψ\Delta_{6}=2\sqrt{\frac{m_{2}m_{r}}{\Omega_{2}\Omega_{r}\gamma\zeta_{j}\psi}}.

Finally, we obtain the e2e ABER expression of U1 as (16) by substituting (9) and (13) into (7). Also, by substituting (11) and (15) into (7), we get the e2e ABER expression of U2 as (17) (please see the top of the next page).

Pe2e,1(e)=1(112i=1212πΓ(mr+0.5)Γ(mr+1)×ϖγζiδ/2mr(1+ϖγζiδ/2mr)(mr+0.5)1(mr,ϖγζiδ))(112i=12Δ12πG3345(0,1u2,0.5u2,1u20.5u1,0.5u1,1u2,u2,0|Δ34))\begin{split}&P_{e2e,1}(e)=\\ &1-\left(1-\frac{1}{2}\sum_{{i}=1}^{2}\frac{1}{\sqrt{2\pi}}\frac{\Gamma({m_{r}+0.5})}{\Gamma(m_{r}+1)}\times\frac{\sqrt{\varpi\gamma\zeta_{i}\delta/2m_{r}}}{(1+\sqrt{\varpi\gamma\zeta_{i}\delta/2m_{r}})^{(m_{r}+0.5)}}\nabla_{1}(m_{r},\varpi\gamma\zeta_{i}\delta)\right)\left(1-\frac{1}{2}\sum_{{i}=1}^{2}\frac{\Delta_{1}}{2\sqrt{\pi}}G\begin{smallmatrix}3&3\\ 4&5\end{smallmatrix}\left(\begin{smallmatrix}0,1-u_{2},0.5-u_{2},1-u_{2}\\ 0.5u_{1},-0.5u_{1},1-u_{2},u_{2},0\end{smallmatrix}\middle|\frac{\Delta_{3}}{4}\right)\right)\\ \end{split} (16)
Pe2e,2(e)=1(112j=15νj12πΓ(mr+0.5)Γ(mr+1)×ϖγζjδ/2mr(1+ϖγζjδ/2mr)(mr+0.5)2(mr,ϖγζjδ))(112j=15νjΔ42πG3345(0,1u4,0.5u4,1u40.5u3,0.5u3,1u4,u4,0|Δ64))\begin{split}&P_{e2e,2}(e)=\\ &1-\left(1-\frac{1}{2}\sum_{{j}=1}^{5}\nu^{j}\frac{1}{\sqrt{2\pi}}\frac{\Gamma({m_{r}+0.5})}{\Gamma(m_{r}+1)}\times\frac{\sqrt{\varpi\gamma\zeta_{j}\delta/2m_{r}}}{(1+\sqrt{\varpi\gamma\zeta_{j}\delta/2m_{r}})^{(m_{r}+0.5)}}\nabla_{2}(m_{r},\varpi\gamma\zeta_{j}\delta)\right)\left(1-\frac{1}{2}\sum_{{j}=1}^{5}\nu^{j}\frac{\Delta_{4}}{2\sqrt{\pi}}G\begin{smallmatrix}3&3\\ 4&5\end{smallmatrix}\left(\begin{smallmatrix}0,1-u_{4},0.5-u_{4},1-u_{4}\\ 0.5u_{3},-0.5u_{3},1-u_{4},u_{4},0\end{smallmatrix}\middle|\frac{\Delta_{6}}{4}\right)\right)\\ \end{split} (17)

IV Numerical Results

In this section, we validate the theoretical analysis by Monte Carlo simulation. In all figures, the lines represent the simulation results, while the markers represent the analytical results. The channel parameters for all simulations are given in Table I and in all simulations, η=0.95\eta=0.95 and for a fair comparison, total transmit SNR is defined as PT/N0\nicefrac{{P_{T}}}{{N_{0}}}. In simulations, we use normalized i.e., T=1 transmission time.

In Fig. 1, we present ABER performances of both users for all protocols, where α1=0.9\alpha_{1}=0.9, and β=ρ=0.1\beta=\rho=0.1. Firstly, it is worthy noting that our analysis matches perfectly with simulations in all protocols and scenarios which proves the correctness of our closed-from expressions. Besides, the analysis refers to both integer and non-integer values of mm and it can be seen that this shape parameter mm drives the diversity order. In Fig. 1, we can observe that the proposed HEH-CNOMA outperforms all three benchmarks. By examining the comparisons in Fig. 1, we observe that the performance in Scenario I is better than Scenario III and in the same manner, the performance in Scenario IV is superior to Scenario II. These observations are a natural expectation of channel quality differences in various scenarios. The shape parameter in Scenario I is greater than Scenario III so that Scenario I has better ABER performance and improved diversity order since in both phases less erroneous detection is performed. Likewise, Scenario IV outperforms Scenario II since the spread parameters are higher which improves the performance. On the other hand, please note that in Scenario I and Scenario IV, the linear sum of spread parameters (i.e., Ωr+Ωk\Omega_{r}+\Omega_{k}) of links are fixed which can reflect the relay position. In other words, the higher Ωr\Omega_{r} means that the relay is closer to the source than the users. By comparing the performances of Scenario I and Scenario IV, we can observe that Scenario IV is superior. Thus, we can easily say that the CNOMA schemes have better performance when the relay is close to the middle of the source and users.

TABLE I: Channel parameters in simulations
Scenario I Scenario II Scenario III Scenario IV
[Ωr,Ω1,Ω2]\left[\Omega_{r},\ \Omega_{1},\ \Omega_{2}\right] [10, 2, 10][10,\ 2,\ 10] [2, 2, 10][2,\ 2,\ 10] [10, 2, 10][10,\ 2,\ 10] [8, 4, 12][8,\ 4,\ 12]
mr=m1=m2m_{r}=m_{1}=m_{2} 1.51.5 1.51.5 11 1.51.5
Refer to caption
(a)
Refer to caption
(b)
Figure 1: ABER performance comparisons between HEH-CNOMA and benchmarks w.r.t. SNR.
Refer to caption
Figure 2: ABER performance comparisons between HEH-CNOMA and benchmarks w.r.t. power allocation.
Refer to caption
(a)
Refer to caption
(b)
Figure 3: ABER performance comparisons between HEH-CNOMA and benchmarks w.r.t. EH parameters.

The power allocation has prevailing effects on the ABER performance. To this end, to discuss the optimal power allocation (PA), in Fig. 2, we present ABER performances of users in Scenario I for PT/N0=20\nicefrac{{P_{T}}}{{N_{0}}}=20 dB. Regardless of the PA changes, the HEH-CNOMA still outperforms all three benchmarks. Nevertheless, as expected, with the increase of α2\alpha_{2}, the performance of U1 is always decayed (i.e., ABER is increasing) since α1\alpha_{1} so the allocated power to U1 symbols is reduced. On the other hand, by increasing α2\alpha_{2}, we firstly have better performance for ‘U2; however, it becomes worse with so much increase in α2\alpha_{2}. This can be explained as follows. By increasing α2\alpha_{2}, we allocate more power to U2 symbols so that the performance is improved. However, U2 should implement SIC (i.e., detect U1 symbols firstly and subtract it) to detect its symbols. Therefore, if we increase α2\alpha_{2}, U1 symbols will be detected erroneously during SIC at the U2 and it will cause an error floor in detection of U2 symbols. Based on the above discussion and Fig. 2, by considering both users’ performances, we can say that the optimum value of power allocation should be 0.1α20.250.1\leq\alpha_{2}\leq 0.25.

In order to better illustrate the effects of the EH parameters (β\beta, ρ\rho) on the performance, in Fig. 3, we present the ABER versus β\beta and ρ\rho in Scenario I and Scenario IV. In both comparisons, we can say that an optimal β\beta-ρ\rho pair can be obtained to minimize the ABER in HEH-CNOMA. According to channel conditions, β\beta and ρ\rho can be chosen jointly. For instance, in Scenario I, since the relay is close to the source, the link between SS-RR has already a good quality and less erroneous detection is performed at the relay. Hence, we can increase β\beta and ρ\rho to increase the harvested energy to guarantee a reliable communication in the second phase by increasing the transmit power of the relay. On the other hand, when the relay is located nearby the middle (i.e., Scenario IV), we should select lower β\beta and ρ\rho values. If we choose β\beta and ρ\rho too much, the information processing in the first phase becomes unreliable since less power is devoted to information transfer. In addition to this, the channel quality of the first phase is relatively worse, so the relay makes more error in detection and this causes an error propagation from relay to users.

V Conclusion

This paper proposes a hybrid protocol-assisted CNOMA, which is a combination of existing EH protocols (i.e., PS and TS). The aim of the study is to improve the performance of the existing EH protocols with a more realistic protocol for practical assumptions. The e2e BER expressions of users are analyzed over the Nakagami-mm fading channel. We discuss the optimal hybrid protocol values with the power allocation factor. The results illustrate that proposed HEH-CNOMA outperforms all benchmarks. The optimum value of EH factors to achieve more satisfactory performance for both users jointly in HEH-CNOMA can be obtained according to spread and shape parameters of the channel.

References

  • [1] Y. Qi, X. Zhang, and M. Vaezi, “Over-the-air implementation of NOMA: New experiments and future directions,” IEEE Access, vol. 9, pp. 135 828–135 844, 2021.
  • [2] M. Vaezi, G. Amarasuriya, Y. Liu, A. Arafa, F. Fang, and Z. Ding, “Interplay between NOMA and other emerging technologies: A survey,” IEEE Trans. Cogn. Commun. Netw., vol. 5, no. 4, pp. 900–919, 2019.
  • [3] N. Ashraf, S. A. Sheikh, S. A. Khan, I. Shayea, and M. Jalal, “Simultaneous wireless information and power transfer with cooperative relaying for next-generation wireless networks: A review,” IEEE Access, vol. 9, pp. 71 482–71 504, 2021.
  • [4] M. Ju, K. M. Kang, K. S. Hwang, and C. Jeong, “Maximum transmission rate of PSR/TSR protocols in wireless energy harvesting DF-based relay networks,” IEEE J. Sel. Areas Commun., vol. 33, no. 12, pp. 2701–2717, 2015.
  • [5] X. Li, M. Liu, C. Deng, D. Zhang, X. C. Gao, K. M. Rabie, and R. Kharel, “Joint effects of residual hardware impairments and channel estimation errors on SWIPT assisted cooperative NOMA networks,” IEEE Access, vol. 7, pp. 135 499–135 513, 2019.
  • [6] D. T. Do and C. B. Le, “Impact of fixed power allocation in wireless energy harvesting NOMA networks,” Internat. J. Commun. Syst., vol. 32, no. 14, pp. 2–6, 2019.
  • [7] V. N. Nayak and K. K. Gurrala, “A novel resource allocation for SWIPT-NOMA enabled AF relay based cooperative network,” Wirel. Pers. Commun., no. 0123456789, 2021.
  • [8] D. T. Do, M. Vaezi, and T. L. Nguyen, “Wireless powered cooperative relaying using NOMA with imperfect CSI,” 2018 IEEE Globecom Workshops, pp. 1–6, 2019.
  • [9] S. Bisen, P. Shaik, and V. Bhatia, “On performance of energy harvested cooperative NOMA under imperfect CSI and imperfect SIC,” IEEE Trans. Veh. Technol., vol. 70, no. 9, pp. 8993–9005, 2021.
  • [10] L. Bariah, S. Muhaidat, and A. Al-Dweik, “Error probability analysis of NOMA-based relay networks with SWIPT,” IEEE Commun. Lett., vol. 23, no. 7, pp. 1223–1226, 2019.
  • [11] L. Mohjazi, S. Muhaidat, and M. Dianati, “Performance analysis of differential modulation in SWIPT cooperative networks,” IEEE Signal Proces. Lett., vol. 23, no. 5, pp. 620–624, 2016.
  • [12] S. Li, L. Bariah, S. Muhaidat, P. Sofotasios, J. Liang, and A. Wang, “Error analysis of NOMA-based user cooperation with SWIPT,” in Internat. Conf. Distr. Comput. Sens. Sys. (DCOSS), 2019, pp. 507–513.
  • [13] S. Li, L. Bariah, S. Muhaidat, P. C. Sofotasios, J. Liang, and A. Wang, “SWIPT-enabled cooperative NOMA with mth best relay selection,” IEEE Open J. Commun. Soc., vol. 1, pp. 1798–1807, 2020.
  • [14] F. Khennoufa, K. Abdellatif, and F. Kara, “Bit error rate evaluation of relay-aided cooperative NOMA with energy harvesting under imperfect SIC and CSI,” Phys. Commun., vol. 52, p. 101630, 2022.
  • [15] S. Atapattu and J. Evans, “Optimal energy harvesting protocols for wireless relay networks,” IEEE Trans. Wirel. Commun., vol. 15, no. 8, pp. 5789 – 5803, 2016.
  • [16] G. Li, D. Mishra, Y. Hu, and S. Atapattu, “Optimal designs for relay-assisted NOMA networks with hybrid SWIPT scheme,” IEEE Trans. Commun, vol. 68, no. 6, pp. 3588–3601, 2020.
  • [17] A. Amin, M. B. Uddin, and S. Y. Shin, “Performance enhancement of hybrid SWIPT protocol for cooperative NOMA downlink transmission,” arXiv preprint arXiv:1906.09159, 2019.
  • [18] A. A. Amin and S. Y. Shin, “Performance evaluation of cooperative NOMA-based improved hybrid SWIPT protocol,” arXiv preprint arXiv:2106.10799, 2021.
  • [19] M. Babaei, Ü. Aygölü, and E. Basar, “BER analysis of dual-hop relaying with energy harvesting in Nakagami-mm fading channel,” IEEE Trans. Wireless Commun, vol. 17, no. 7, pp. 4352–4361, 2018.
  • [20] F. Khennoufa, K. Abdellatif, and F. Kara, “Bit error rate and outage probability analysis for multi-hop decode-and-forward relay-aided NOMA with imperfect SIC and imperfect CSI,” AEU-Internat. J. Electron. and Commun., vol. 147, p. 154124, 2022.
  • [21] M. K. Simon and M.-S. Alouini, Digital Communication over Fading Channels.   Wiley-IEEE Press, 2005.
  • [22] M. Babaei, U. Aygolu, and E. Basar, “BER analysis of dual-hop relaying with energy harvesting in Nakagami-m fading channel,” IEEE Trans. Wirel. Commun., vol. 17, no. 7, pp. 4352–4361, 2018.
  • [23] I. Gradshteyn and I. Ryzhik, Table of Integrals, Series, and Products.   Academic Press.
  • [24] F. Kara, “Error performance of cooperative relaying systems empowered by SWIPT and NOMA,” Phys. Commun., vol. 49, p. 101450, 2021.