Cooperative Non-Orthogonal Multiple Access with Energy Harvesting

In this paper, we analyze the outage probability of an amplify-and-forward (AF) cooperative non-orthogonal multiple access (NOMA) model in multi-relay multiuser networks. In contrast to conventional cooperative networks, relays in the considered network have no embedded energy supply; they need to rely on the energy harvested from the signals broadcasted by the source for their cooperative NOMA transmission. Based on this structure, a new relay selection scheme is proposed, considering both channel state information (CSI) and battery status of relays. Assuming each relay has infinite or finite energy storage for accumulating energy, we use the infinite or finite Markov chain to capture the evolution of relay batteries and certain simplified assumptions to reduce computational complexity of the Markov chain analysis. The approximate closed-form expressions for the outage probability of the proposed scheme are derived therefrom. All theoretical results are validated by numerical simulations. The impacts of the system parameters, such as relay number, energy harvesting threshold and battery size, on the performance are extensively investigated.


Introduction
With the aim of increasing system capacity and achieving higher spectral efficiency, non-orthogonal multiple access (NOMA) has been recognized as a promising multiple access technique for fifth generation (5G) networks due to its superior spectral efficiency [1].The key idea of NOMA is to serve multiple users in the same frequency band, but with different power levels, which is a domain fundamentally different from conventional orthogonal access technologies [2].In particular, power-domain NOMA allocates less transmitting power to users with better channel conditions and more transmit power to users with worse channel conditions in order to achieve a balanced tradeoff between system throughput and user fairness.Therefore, users can be separated by successive interference cancellation (SIC).Lan et al. [3] explored the impact of the error propagation of SIC and user velocity on the NOMA performance.Their results showed that, even in the worst error propagation scenario, NOMA outperforms conventional orthogonal multiple access and can yield performance gains for different user mobility.Chen et al. [4] studied NOMA for the downlink of a wireless system.Traditional minimum-mean-squared-error (MMSE) precoding matrices have been used.Considering that wireless relaying is an effective means to combat the effects of the fading and shadowing on transmitted signal, it is reasonable to combine NOMA with relaying networks.In [5], a new cooperative NOMA scheme was proposed and analyzed in terms of outage probability and diversity gain.The authors in [6] study the impact of relay selection on the performance of cooperative NOMA.In particular, a two-stage relay selection strategy is proposed; this strategy can achieve the minimal outage probability among all possible selection schemes, and realize the maximal diversity gain.
In addition to improving spectral efficiency, which is the motivation of NOMA, another key objective of future 5G networks is to maximize energy efficiency.Energy harvesting (EH), a technique to collect energy from the surrounding environment, has recently received considerable attention as a sustainable solution to overcome the bottleneck of energy-constrained wireless sensor networks.Apart from the conventional renewable energy sources, such as solar and wind, radio frequency (RF) signals radiated by ambient transmitters can be treated as a viable new source for EH.Such an approach can reduce the cost of a sensor network, as peripheral equipment can be avoided [7].However, a fundamental limitation of EH-based wireless communications lies in the restricted transmission range.Although longer ranges can be achieved through a stronger RF source, the available energy level to pick up at the distant receiver remains fairly small due to pathloss.In cellular and sensor networks, relays can be deployed to extend the coverage of base stations [8].Recently, relays with EH capabilities have received much attention as they use the energy harvested from the source signal to perform information forwarding [9].This can solve the problem of the energy supply of relays and expand the application of EH-based wireless communications.In [10], a wireless cooperative network is considered in which multiple source-destination pairs communicate with each other via an energy harvesting relay.In [11], both the source and the relay are EH devices and are charged by the destination serving as the power station.In [12], a novel best cooperative mechanism for wireless energy harvesting and spectrum sharing has been proposed, and this mechanism has been verified to be superior to the traditional schemes by simulation.The aforementioned literature concentrates on the achievable performance without considering energy storage at relays; i.e., the harvested energy within a transmission block is entirely consumed for forwarding information.Nevertheless, the energy harvested from RF radiation is often restricted, and thus it is desirable for relays to accumulate the harvested energy in energy storage such as super-capacitors or rechargeable batteries [13].In [14], a threshold-based "save-then-transmit" scheme is employed at relays; the stored energy level in each relay battery actually forms a Markov chain over time.By investigating the properties of this Markov chain, the asymptotical average throughput is derived.In [15], to support an efficient utilization of harvested energy to improve throughput, a harvest-use-store relaying strategy with distributed beamforming has been researched.
The aforementioned three communication concepts, NOMA, EH and cooperative communication can be combined naturally, which is the focus of this paper.In particular, we consider a NOMA-based downlink amplify-and-forward (AF) relaying network.The main contributions of this paper are summarized as follows: (1) We design a new NOMA-based relaying network and propose a new EH protocol for relays.
The relays have no other energy supplies, but they are equipped with a chargeable battery, and thus can harvest and store the wireless energy broadcasted by the source.Then, we model the capacity of the relay battery in infinite and finite cases, respectively.The outage behavior of the network is investigated, and the closed-form expressions for the outage probability are derived.(2) Simulation results are conducted to demonstrate our analytical results and the superiority of NOMA.In addition, the impacts of system parameters on the performance of the network are captured by simulation.(3) Finally, we compared the proposed scheme with a conventional multiple access (MA) cooperative scheme.Simulation results show that, although conventional MA obtains better outage performance than our scheme, our scheme can offer better fairness, since more users are served simultaneously.In addition, our scheme can guarantee acceptable system performance (the best outage probability is 10 −5 when transmission signal to noise ratio (SNR)is 20dB) even if the relays do not use their own batteries to power the relay transmission which demonstrates the superiority of our scheme compared to the common cooperative system or common NOMA system.

System Model
As shown in Figure 1, a source S and a number of potential users D j (j = 1, 2, ..., N) communicate over channels with flat fading.
Multiple potential relays are willing to amplify and forward the signal from S to the users.S has no direct link with users.The channels pertaining to the first hop and second hop undergo independent identical (i.i.d.) fading and the channel coefficients are denoted by h SR i and h R i D j , respectively.Assuming Rayleigh fading, h SR i is a circularly symmetric complex Gaussian random variable with zero mean and variance Likewise, . The channel power gain SR and σ 2 RD .Without loss of generality, we assume that the users' channel gains have been ordered as We also assume that the additive white Gaussian noise (AWGN) of all links has a zero mean and equal variance N 0 .

System Model
As shown in Figure 1, a source S and a number of potential users ( 1, 2,..., )  are willing to amplify and forward the signal from S to the users.S has no direct link with users.The channels pertaining to the first hop and second hop undergo independent identical (i.i.σ .Without loss of generality, we assume that the users' channel gains have been ordered as . We also assume that the additive white Gaussian noise (AWGN) of all links has a zero mean and equal variance In this paper, we use the conventional AF protocol for the NOMA downlink cooperative network.From Figure 1, we can see that each transmission block takes two time slots; for convenience but without loss of generality, we consider a normalized unit block time (i.e., is the transmission power and is sufficiently large such that the energy harvested from the noise is negligible.Thus, the amount of energy harvested from the source can be expressed as [16]  In this paper, we use the conventional AF protocol for the NOMA downlink cooperative network.From Figure 1, we can see that each transmission block takes two time slots; for convenience but without loss of generality, we consider a normalized unit block time (i.e., T = 1) hereafter, so each time slot takes T/2 and the energy harvesting time is also T/2.Before the transmission, each relay checks their battery at the beginning of a transmission block and sees if it has enough energy to forward the source information.If the relay does not have enough energy, it performs EH in this time block and stores the harvested energy into the individual battery.We assume z = y − θ * n is the transmission power and is sufficiently large such that the energy harvested from the noise is negligible.Thus, the amount of energy harvested from the source can be expressed as [16] where η is the energy harvesting efficiency.
For those relays with sufficient energy, they report their CSI to S for making the relay selection decision.Let e i denote the battery energy amount of relay R i ; we define those relays with enough energy as the eligible set where C is the energy harvesting threshold to activate the EH circuit and W represents the number of times transmit energy over C.Among this set of relays, R m is selected, which can be expressed as m = arg max i:R i ∈φ In Equation ( 3), the relay with the best channel power gain of the first hop in this set is selected for forwarding information.
The two-phase communication starts after relay selection.In the first slot, S transmits the unit-power superposition symbol x S = ∑ N j=1 a j P S x j to the selected R m , where x j contains the information required by D j , a j is the power allocation coefficient of D j which denotes the proportion of the transmit power allocated to x j .Following the principle of NOMA, we have The observation at R m can be expressed as where n R m is AWGN at R m .In this time slots, all other relays regard x S as energy signal and perform EH.
In the second slot, R m amplifies y R m with an amplifying coefficient β and then broadcasts it to all users.The signal received by the nth user can be expressed as where P R is the transmit power at relay and P R = WC/(T/2), n D n is AWGN at nth user, and β should be determined as follows: For each user, the desired signal is interfered by the other users' signals.SIC will be carried out at each user to mitigate the negative effect of the inter-user interference.The SIC decoding order is in increasing order of the effective users' channel gains (h . Therefore, at the nth user, the lth user's signal, l < n, will be detected and then cancelled out from the received signal of the nth user in a successive manner.The lth user's signal, l > n, will be treated as noise at the nth user.If we set γ 1 = P S /N 0 and γ 2 = P R /N 0 , then, we can compute the signal to interference and noise ratio (SINR) for the nth user to decode the lth user's signal, l < n, as follows: Accordingly, the instantaneous rate can be written as R l→n = 1 2 log 2 1 + γ R m D l→n .If the message x l can be decoded successfully, i.e., R l→n ≥ R l = 1 2 log 2 1 + γ l th , it will be removed from the nth user's observation, where R l and γ l th denote the target data rate and the target SINR for the lth user, respectively.This SIC will be implemented until n users' messages are all decoded, where the SINR for the nth user to decode its own signal is given by The nth user needs to decode all the other users' signals and the SNR for the nth user to decode its own signal can be expressed as Hence, the achievable data rate for the nth user is given by

Performance Analysis
In this section, the outage performance of the NOMA cooperative network will be characterized in terms of outage probability.According to the total probability law, the outage probability of the nth user can be written as where |•| denotes the number of elements in a set and P n,Ω out is the conditional outage probability when the number of φ is Ω.We will first study the P n,Ω out and then Pr(|φ| = Ω), finally, P n out will be derived.For ease of description, th , l ≤ n, is defined as the outage event where the nth user cannot decode the lth user's signal successfully, and Λ c n,l is defined as the complementary set of Λ n,l .As a result, the conditional outage probability for the nth user can be given by To proceed forward, we first rewrite Λ c N,N as Next, the other events Λ c n,l , l = N, can be attained as From the final step of Equation ( 13) we can obtain a necessary condition is: If the condition is not satisfied, the nth user can never decode the lth user's signal successfully irrespective of the channel SINR.
By defining θ * n = max{θ 1 , θ 2,••• θ n }, the outage probability for the nth user can be reformulated as Since h SR i follows the exponential distribution with mean σ 2 SR , the probability density function (PDF) of h SR m can be written as Similarly, the PDF of the ordered variable h R m D n is given by In what follows, J 1 and J 2 will be addressed.Firstly, according to Equation ( 16), J 1 can be calculated as On the other hand, by substituting Equations ( 16) and ( 17) into J 2 , it can be written as Defining = v, J 3 can be attained as Defining z = y − θ * n and with the aid of Equation (3.478.4) in [17], J 3 can be calculated as follows Until now, we have been able to obtain a conditional exact closed-form expression for the conditional outage probability P n,Ω out by substituting Equations ( 18) and ( 21) into (15).The next step is to get Pr(|φ| = Ω), which will be given in two cases based on energy storage capacity of relay.

Infinite Storage of Energy
In this case, we assume that each relay has infinite battery and define C as the energy harvesting threshold to activate the EH circuit.Although harvested energy can be more than C in one time slot, here, we consider each relay to only store C and use the rest for its own purpose.The eligible set is composed by relays whose amount of harvested energy are no less than WC.We model the amount of the stored energy at each relay using an infinite state Markov chain as Figure 2. When the relay does not accumulate enough energy, i.e., the stored energy is smaller than WC, which corresponds to states s = 0, 1, 2 . . . .W − 1 in Figure 2, the relay cannot forward the information and there is no transition back to the previous states.When the relay has enough stored energy, i.e., the relay is in state s ≥ W in Figure 2 and is selected as the best relay, its stored energy state transits to state s − W after the relaying.
When the relay does not accumulate enough energy, i.e., the stored energy is smaller than WC , which corresponds to states 0,1, 2.... 2, the relay cannot forward the information and there is no transition back to the previous states.When the relay has enough stored energy, i.e., the relay is in state s W ≥ in Figure 2 and is selected as the best relay, its stored energy state transits to state s W − after the relaying.When s < W, the transition from state s to state s + 1 happens only when the relay harvest no less than C amount of energy in one time slot.Hence, using Equation (1), we have where P a,b is the transition probability from state a to state b.Similarly, the transition from state s to state s happens when the relay performs EH but the harvested energy is not enough to increase its battery by one level, resulting in When s ≥ W, the transition from state s to state s + 1 happens when the relay harvest more than C amount of energy in one time slot and this relay is not selected at this time slot.We note that the i.i.d.fading assumption implies each relay in φ has an equal chance to be selected as the best relay, and the average number of relays in φ is approximately equal to M because of infinite storage capacity of relay.We denote Pr(R m = R i ) as the probability that R i is selected as the best relay, so Pr(R m = R i ) ≈ 1/M.This approximate expression may not always hold; we will discuss the effectiveness of the assumption through numerical results.Using this expression, the corresponding transition probability P s,s+1 is given by Similar to the above case, when s ≥ W, the transition from state s to state s happens when the relay is not selected and the harvested energy is less than C at this time slot.In this case, P s,s can be written as If the relay which is in state s ≥ W is selected as the best relay for transmitting, its harvested energy then transits from state s ≥ W to state s − W. Hence, the transition probability P s,s−W is given by With the transition probabilities derived in Equations ( 22)-( 26), we form an infinite dimensional transition matrix as P = [P a,b ] (a, b ∈ N), which can be written as Let v = (v 1 , v 2 , ...) denote the steady state probability vector; then, in the steady state, we have vP = v (28) By expanding Equation (28) and making correspondence on both sides, we notice that Based on the property of this transition matrix, we have Then, we have the following relation after some algebraic computation: Summing all the terms on both side of Equation (32), and using the fact that ∑ ∞ k=1 v k = 1, we can obtain v k , ∀k ∈ N.Then, the probability that an relay can be in eligible set, is equal to the probability that the relay is in state s ≥ W, which is given by With the above result, we can get Pr(|φ| = Ω), which follows the binomial distribution with the probability mass function given as With P n,Ω out and Equation (34), we can get the close form expression of outage probability.

Finite Storage of Energy
In this case, each relay accumulates the harvested energy using a finite energy storage with the size B. We also define C as the energy harvesting threshold and WC is the energy required for forwarding data, but relay can harvest UC amount of energy in one time slot and U is given by This assumption is closer to the practical scenario, and the evolution of the battery status of all relays can be modeled as a finite state Markov chain; using the transition probability matrix of this chain, we can get the steady state probability vector and Pr(|φ| = Ω).Once Pr(|φ| = Ω) is obtained, the outage probability can also be obtained.
Due to the lack of general form of steady state probability, the above analysis is computationally intense when M is large.To facilitate the computation, we propose an approximated approach based on two simplified assumptions.Firstly, we denote the relay energy amount at the selection epoch as a random variable Z.The exact distribution of Z is high computational complexity.To ease the computation, we approximate Z by a uniform random variable over [0, B].The adopted approximation is inspired by considering the amount of harvested energy in a transmission block follows the geometric distribution with parameter 1/2 [18].In general, these conditions may not always hold.We will discuss the effectiveness of the assumption through numerical results.Secondly, we found that an arbitrary relay may be either short of enough power to participate in relay selection or otherwise, so the evolution of relay energy amount is captured by using two states, either active or inactive.With this simplified two-state Markov chain, a relay is in s 0 if the relay lacked sufficient energy to transmit, or in s 1 when the relay has enough energy for transmission.Next, we explain how to obtain the transition probability matrix of the two-state Markov chain.
The transition from state s 0 to state s 0 happens when a relay has no enough energy to transmit (i.e., Z < WC) in the current block and the accumulated energy after harvesting remains below WC.The corresponding transition probability is given by where Ẑ is a truncated random variable defined as Since Z is approximated as uniformly distributed, the PDF of Ẑ can be obtained as where u(•) and δ(•) denote the unit step function and the Dirac delta function, respectively.Then Equation (36) can be solved as where E R i = E E R i = P S ησ 2 SR T/2.The transition from s 0 to s 1 happens when the relay enters the EH mode in the current bock and the accumulated energy exceeds WC.Hence, we have Similar to the derivation of Equation (36), we obtain If the relay which is in state s 1 is selected as the best relay for transmitting, its harvested energy then transits from state s 1 to state s 0 .Hence, the transition probability p 1,0 is given by where Z is a truncated random variable defined as Since Z is uniformly distributed, the PDF of Z can be obtained as Using this PDF, the first term in the right of Equation ( 42) can be obtained as As to the second term, we note that the i.i.d.fading assumption implies each relay in φ has an equal chance to be selected as the best relay.To simplify the analysis, we approximate the cardinality of φ by its mean such that where M e is the average number of relays in φ and v 1 is the steady-state probability of state s 1 .
The transition probability from s 1 to s 1 can be solved in the similar manner as the previous case, so we omit the derivation here.
When the two-state Markov chain formulated, the steady state probability vector can be easily obtained as We note that both v 0 and v 1 involve p 1,0 , which is a function of v 1 .By substituting p 1,0 and p 1,0 into Equation (47), v 1 can be solved explicitly.Take the condition 2W < B in Equation (45), for example; v 1 can be obtained in closed form as With v 1 at hand, we can get Pr(|φ| = Ω), which follows the binomial distribution with the probability mass function given as Submitting Pr(|φ| = Ω) into Equation ( 10), we can get outage probability of the system.

Simulation Results
In this section, computer simulations are performed to validate the mentioned theoretical analysis.In all simulations, we set the noise power to N 0 = 1, the energy harvesting efficiency to η = 0.5, and the fixed transmission rate of the source is 1 bit per channel use (bpcu).The battery size B is set to be a multiple of the energy harvesting threshold C, i.e., B = αC, where α ∈ N and α > W. We also set C is the multiple of the source transmission energy, i.e., C = δPT/2, where δ > 0 is the scaling factor and the length of a transmission block T = 1.To facilitate the analysis, σ 2 SR and σ 2 RD are set to be 1.Without loss of generality, we set N = 3, the power allocation coefficient are a 1 = 1/2, a 2 = 1/3 and a 3 = 1/6, the outage threshold are γ 1 th = 0.9 dB, γ 2 th = 1.5 dB and γ 3 th = 2 dB.Figures 3 and 4 illustrate the performance of the proposed protocol with infinite battery size and finite battery size versus the transmitting SNR of source for given values of C and W, respectively.We set W = 4 and δ = 0.5.In Figures 3 and 4, the approximated analysis of infinite and finite Markov chain is accurate when the number of relays is sufficient, but its results become deviated from thesimulated ones for M = 3.The outage probability reveals an error floor when M = 3.This is because, when M = 3, relays cannot maintain a enough energy state in the infinite storage case, and the energy distribution of the relay battery is no longer as uniform as assumed in the finite storage case, respectively.However, when M increases, the approximate cardinality of φ becomes more accurate.We can also see that the slope of the outage probability curve (namely, diversity order) increases with M. Since the theoretical analyses agree well with the simulations in medium and high SNR ranges, we will only plot the analytical results in the remaining figures. .In Figures 3 and 4, the approximated analysis of infinite and finite Markov chain is accurate when the number of relays is sufficient, but its results become deviated from the simulated ones for 3 M = .The outage probability reveals an error floor when 3 M = .This is because, when 3 M = , relays cannot maintain a enough energy state in the infinite storage case, and the energy distribution of the relay battery is no longer as uniform as assumed in the finite storage case, respectively.However, when M increases, the approximate cardinality of φ becomes more accurate.We can also see that the slope of the outage probability curve (namely, diversity order) increases with M .Since the theoretical analyses agree well with the simulations in medium and high SNR ranges, we will only plot the analytical results in the remaining figures.Next, we will investigate the impacts of the system parameters on the performance in medium SNR conditions (P S /N 0 = 20dB).In Figures 5 and 6, we illustrate the impacts of C on the performance of the proposed two battery scheme with different W. M = 6 is set in both two figures.For all the curves in Figures 5 and 6, the trends are the same for all schemes; namely, the probability first decreases then increases as δ varies from 0 to 1.This means that, when the other parameters are determined, there must be an optimal value of δ.However, the values of the inflection points are not always the same for different users and are inversely proportional to W. When W = 1, the eligible set has more relays than W = 3; the probability of choosing the best relay with more energy is larger than W = 3, and so the optimal value of δ can be greater than W = 3.The optimal value of δ can easily be obtained by a one-dimensional exhaustive search; with this optimal value of δ, the system can resist fading more effectively.Next, we will investigate the impacts of the system parameters on the performance in medium SNR conditions ( ).In Figures 5 and 6, we illustrate the impacts of C on the performance of the proposed two battery scheme with different W . 6 M = is set in both two figures.For all the curves in Figures 5 and 6, the trends are the same for all schemes; namely, the probability first decreases then increases as δ varies from 0 to 1.This means that, when the other parameters are determined, there must be an optimal value of δ .However, the values of the inflection points are not always the same for different users and are inversely proportional to W . 1 W = , the eligible set has more relays than 3 W = ; the probability of choosing the best relay with more energy is larger than 3 W = , and so the optimal value of δ can be greater than 3 W = .The optimal value of δ can easily be obtained by a one-dimensional exhaustive search; with this optimal value of δ , the system can resist fading more effectively.).In Figures 5 and 6, we illustrate the impacts of C on the performance of the proposed two battery scheme with different W . 6 M = is set in both two figures.For all the curves in Figures 5 and 6, the trends are the same for all schemes; namely, the probability first decreases then increases as δ varies from 0 to 1.This means that, when the other parameters are determined, there must be an optimal value of δ .However, the values of the inflection points are not always the same for different users and are inversely proportional to W .When 1 W = , the eligible set has more relays than 3 W = ; the probability of choosing the best relay with more energy is larger than 3 W = , and so the optimal value of δ can be greater than 3 W = .The optimal value of δ can easily be obtained by a one-dimensional exhaustive search; with this optimal value of δ , the system can resist fading more effectively.In Figure 7, we investigate the impact of battery size on the proposed finite battery scheme by varying the battery scaling factor with fixed number of relay (M = 6).Since three users have the same trend, we only analyze user 3 in this figure.From Figure 7, we can observe that the performance increases as α increases.However, the gain provided by a larger battery size does not increase when α exceeds a certain value and this value is decided by W. From this figure, we can also see that the order of performance for different W changes when the value of B changes.In Figure 7, we investigate the impact of battery size on the proposed finite battery scheme by varying the battery scaling factor with fixed number of relay ( 6 M = ).Since three users have the same trend, we only analyze user 3 in this figure.From Figure 7, we can observe that the performance increases as α increases.However, the gain provided by a larger battery size does not increase when α exceeds a certain value and this value is decided by W . From this figure, we can also see that the order of performance for different W changes when the value of B changes.From those two figures, we can see that, when M is small, there is a huge drop of outage probability with a slight increase of M ; this is because relays in the eligible set are insufficient when M is small, and the slight increase of M will increase the number of relays in the eligible set remarkably and a better relay can be selected.However, when M is huge ( 10 M =  In Figure 7, we investigate the impact of battery size on the proposed finite battery scheme by varying the battery scaling factor with fixed number of relay ( 6 M = ).Since three users have the same trend, we only analyze user 3 in this figure.From Figure 7, we can observe that the performance increases as α increases.However, the gain provided by a larger battery size does not increase when α exceeds a certain value and this value is decided by W . From this figure, we can also see that the order of performance for different W changes when the value of B changes.From those two figures, we can see that, when M is small, there is a huge drop of outage probability with a slight increase of M ; this is because relays in the eligible set are insufficient when M is small, and the slight increase of M will increase the number of relays in the eligible set remarkably and a better relay can be selected.However, when M is huge ( 10 M =  From those two figures, we can see that, when M is small, there is a huge drop of outage probability with a slight increase of M; this is because relays in the eligible set are insufficient when M is small, and the slight increase of M will increase the number of relays in the eligible set remarkably and a better relay can be selected.However, when M is huge (M = 10 or M = 20), relays in the eligible set are sufficient, and the increase of M enhances the performance slightly, thus our algorithm can be used when M is not more than 10.In addition, we compared the proposed scheme with conventional MA cooperative scheme.Here, we assume that an opportunistic MA approach is adopted for conventional MA, where the user with the best channel condition is scheduled.The target SINR γ th for conventional scheme satisfies 1  2 ∑ 3 i=1 log 2 1 + γ i th = 1 2 log 2 (1 + γ th ) Simulation results show that conventional MA can achieve the same diversity gain with user 3. Conventional MA obtains better outage performance than NOMA, but NOMA can offer better fairness since more users are served simultaneously.
performance slightly, thus our algorithm can be used when M is not more than 10.In addition, we compared the proposed scheme with conventional MA cooperative scheme.Here, we assume that an opportunistic MA approach is adopted for conventional MA, where the user with the best channel condition is scheduled.The target SINR th γ for conventional scheme satisfies

Conclusions
In this paper, we analyze the performance of the cooperative NOMA network.A new relay selection scheme is proposed considering both CSI and battery status of relays.We model the amount of harvested energy at each relay using an infinite or finite Markov chain and then derive the approximate closed-form expression of the outage probability in two cases, respectively.or 20 M = ), relays in the eligible set are sufficient, and the increase of M enhances the performance slightly, thus our algorithm can be used when M is not more than 10.In addition, we compared the proposed scheme with conventional MA cooperative scheme.Here, we assume that an opportunistic MA approach is adopted for conventional MA, where the user with the best channel condition is scheduled.The target SINR th γ for conventional scheme satisfies

Conclusions
In this paper, we analyze the performance of the cooperative NOMA network.A new relay selection scheme is proposed considering both CSI and battery status of relays.We model the amount of harvested energy at each relay using an infinite or finite Markov chain and then derive the approximate closed-form expression of the outage probability in two cases, respectively.

Conclusions
In this paper, we analyze the performance of the cooperative NOMA network.A new relay selection scheme is proposed considering both CSI and battery status of relays.We model the amount of harvested energy at each relay using an infinite or finite Markov chain and then derive the approximate closed-form expression of the outage probability in two cases, respectively.Simulations are carried out to verify the correctness of the theoretical analysis.We conclude the superiority of NOMA and find that, by carefully choosing the parameters of the network, (e.g., energy harvesting threshold or forwarding threshold), acceptable system performance can be guaranteed even if the relays do not use their own batteries to power the relay transmission.For future work, we will investigate the performance of cooperative NOMA in large-scale relaying networks using stochastic geometry.

2 T
and the energy harvesting time is also 2 T .Before the transmission, each relay checks their battery at the beginning of a transmission block and sees if it has enough energy to forward the source information.If the relay does not have enough energy, it performs EH in this time block and stores the harvested energy into the individual battery.We assume

Figure 1 .
Figure 1.A reference model of system.

Figure 2 .Figure 2 .
Figure 2. State transition diagram of the harvested energy amount.(The state s denotes the relay has sC amount of harvested energy at a relay.) Figure 2. State transition diagram of the harvested energy amount.(The state s denotes the relay has sC amount of harvested energy at a relay.)

5 , 2 SR σ and 2 RDσ 3 NFigures 3
Figures 3 and 4 illustrate the performance of the proposed protocol with infinite battery size and finite battery size versus the transmitting SNR of source for given values of C and W , respectively.We set 4 W = and

Figure 3 .
Figure 3. Outage probability with different M (infinite storage of energy).Figure 3. Outage probability with different M (infinite storage of energy).

Figure 3 .
Figure 3. Outage probability with different M (infinite storage of energy).Figure 3. Outage probability with different M (infinite storage of energy).

Figure5.
Figure5.Outage probability with different W (infinite storage of energy).Figure 5. Outage probability with different W (infinite storage of energy).

Figure 5 .
Figure5.Outage probability with different W (infinite storage of energy).Figure 5. Outage probability with different W (infinite storage of energy).

Figures 8 and 9
Figures 8 and 9 illustrate the performance of the proposed protocol for different relay number.We also only analyze user 3 in both two figures.
are set.In order to illustrate the scalability of our algorithm, we set 3 M = , 5 M = , 10 M = , 20 M = and 30 M = , respectively.

Figures 8 and 9
Figures 8 and 9 illustrate the performance of the proposed protocol for different relay number.We also only analyze user 3 in both two figures.W = 4 and δ = 0.5 are set.In order to illustrate the scalability of our algorithm, we set M = 3, M = 5, M = 10, M = 20 and M = 30, respectively.From those two figures, we can see that, when M is small, there is a huge drop of outage probability with a slight increase of M; this is because relays in the eligible set are insufficient when M is small, and the slight increase of M will increase the number of relays in the eligible set remarkably and a better relay can be selected.However, when M is huge (M = 10 or M = 20), relays in the eligible set are sufficient, and the increase of M enhances the performance slightly, thus our algorithm can be used when M is not more than 10.In addition, we compared the proposed scheme with conventional MA cooperative scheme.Here, we assume that an opportunistic MA approach is adopted for conventional MA, where the user with the best channel condition is scheduled.The target SINR γ th for conventional scheme satisfies1  2 ∑ 3 i=1 log 2 1 + γ i th = 1 2 log 2 (1 + γ th ) Simulation results show that conventional MA

.
Simulation results show that conventional MA can achieve the same diversity gain with user 3. Conventional MA obtains better outage performance than NOMA, but NOMA can offer better fairness since more users served simultaneously.

Figure 8 .
Figure 8.The impact of relay number on the outage probability (infinite storage of energy).

Figure 9 .
Figure 9.The impact of relay number on the outage probability (finite storage of energy, 15 α = ).

Figure 8 .
Figure 8.The impact of relay number on the outage probability (infinite storage of energy).

.
Simulation results show that conventional MA can achieve the same diversity gain with user 3. Conventional MA obtains better outage performance than NOMA, but NOMA can offer better fairness since more users are served simultaneously.

Figure 8 .
Figure 8.The impact of relay number on the outage probability (infinite storage of energy).

Figure 9 .
Figure 9.The impact of relay number on the outage probability (finite storage of energy, 15 α = ).

Figure 9 .
Figure 9.The impact of relay number on the outage probability (finite storage of energy, α = 15).