Next Article in Journal
Fast and Accurate Object Detection in Remote Sensing Images Based on Lightweight Deep Neural Network
Previous Article in Journal
An Ensembled Anomaly Detector for Wafer Fault Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time Switching Based Wireless Powered Relay Transmission with Uplink NOMA

1
Faculty of Information Technology, Macau University of Science and Technology, Avenida Wai Long, Taipa 999078, Macau, China
2
College of Electronic and Information Science, Fujian Jiangxia University, Fuzhou 350108, China
3
QianJiang College, Hangzhou Normal University, Hangzhou 310012, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(16), 5467; https://doi.org/10.3390/s21165467
Submission received: 6 July 2021 / Revised: 31 July 2021 / Accepted: 9 August 2021 / Published: 13 August 2021
(This article belongs to the Section Sensor Networks)

Abstract

:
Non-orthogonal multiple access (NOMA) utilizes power domain multiplexing to improve spectrum efficiency compared with orthogonal multiple access (OMA). In the Internet of Things (IoT) uplink NOMA networks, if the channel between the far-end node and the base station is in deep fading, allocating larger transmitting power for this node cannot achieve higher spectrum efficiency and overall system throughput. Relay cooperative communication reduces the transmitting power at the far-end node but leads to extra energy expenditure at the relay node. Fortunately, simultaneous wireless information and power transfer (SWIPT) is advocated in energy-constrained IoT networks to save energy consumption. However, early works all focus on energy harvesting (EH) from one source node or one dedicated power supply station. In this paper, we propose a time switching based wireless powered relay transmission model with uplink NOMA where our EH technique can harvest energy from two simultaneously transmitting nodes. More importantly, by optimizing relay position more energy is harvested from the near-end node at the relay and relay signal attenuation to the destination is reduced as well. Furthermore, the closed-form expressions of outage probability and overall system throughput are derived, and numerical results prove that NOMA in our EH scheme achieves better performance compared to the traditional EH scheme and OMA by optimizing the position of the relay node, time switching factor and so on.

1. Introduction

With tremendous deployment of machine-type nodes for the Internet of Things (IoT) applications in the coming years, wireless communication technology has faced severe scarcity of radio spectrum resources [1]. Conventional orthogonal multiple access (OMA) communication techniques, such as time division multiple access (TDMA), and orthogonal frequency division multiple access (OFDMA) require one orthogonal resource block to serve for one single node [2]. In a typical OFDMA scenario, a far-end node with poor channel condition has to be allocated enough subcarriers to meet its quality of service (QoS) requirement, while the near-end node with the perfect channel is forced to use the remaining narrow bandwidth, which may greatly reduce the spectrum utilization and overall system throughput. Unlike OMA, non-orthogonal multiple access (NOMA) utilizes power domain multiplexing to implement scarce spectrum resource sharing, where the transmitting signals of different nodes are assigned at different power levels but at the same time, frequency, and code [3]. For example, in a downlink NOMA scenario, a base station (BS) serving for multiple nodes simultaneously will assign a power ratio for each node in a superimposed transmitting signal. The key principle of NOMA is that more power is allocated for the far-end node, and less power for the near-end node in a given total power constraint since long distance transmission will cause the larger path loss. After receiving the superimposed signal, the near-end node first decodes the signal of the far-end node, and then decodes its own signal by successive interference cancellation (SIC), while the far-end node only decodes its high-power signal by treating other signals as interference. In an uplink NOMA scenario, all the messages from different nodes are transmitted simultaneously to the BS in the same frequency band. The strong signal, which may belong to the near-end node or the good channel condition, will be first decoded at the BS. The weak signal from the far-end node is decoded at last. Therefore, NOMA can achieve higher spectrum efficiency, connectivity and fairness, which is significant for 5G wireless networks to support IoT functionalities.
In terms of NOMA system performance, superior outage probability (OP) and sum rate are demonstrated by optimizing the power allocation scheme compared with OMA [3,4]. The work in [5] has further investigated NOMA performance in cognitive radio networks to guarantee QoS of the nodes with poor channel conditions. However, if the channel between the far-end node and the BS is in deep fading, larger energy allocation for this node at downlink NOMA is unable to achieve high spectrum efficiency and overall system throughput, but increase energy consumption. Moreover, it is unfeasible for a battery-operated IoT node to consume enormous transmitting power in uplink NOMA. Fortunately, relay cooperative communication has been widely employed to cope with wireless channel impairments, e.g., fading and interference [6,7]. For example, in a relay cooperative network, if no direct link exists between the source and the destination, with the aid of relay nodes, the transmission distance per hop is shortened. Due to lower path loss and lower interference, transmitting power is also reduced. Furthermore, cooperative communication may enhance the signal-to-noise ratio (SNR) due to two copies received at the destination if the destination can also receive the message from the source directly.
Relay cooperative communication techniques have been introduced to NOMA in [8,9,10,11]. Some suboptimal power allocation strategies are proposed and system performance, e.g., OP and achievable average rate, is analyzed in the cooperative relay NOMA system [8,9]. The work [10] adopts user pairing selection opportunistically to operate in direct NOMA mode or cooperative NOMA mode to guarantee transmission reliability and fairness. Additionally, an adaptive switching approach between OMA mode and cooperative NOMA mode is realized to improve the far-end node’s performance in [11].
Although cooperative relaying techniques increase overall system performance, forwarding also leads to extra energy expenditure at the relay node, which may prevent energy-constrained IoT nodes from participating in relaying operation [12]. Fortunately, energy harvesting (EH) techniques have emerged to reduce energy consumption in energy-constrained wireless IoT or sensor networks [13]. In particular, a new promising solution to harvest energy from ambient radio frequency (RF) signals is encouraged in wireless networks since RF signals can carry both power and information at the same time, thus those energy-constrained nodes can obtain energy and decode the information simultaneously. In a simultaneous wireless information and power transfer (SWIPT) system, an IoT relay node can scavenge energy and receive information from the wireless signals of the source node, and then use the harvested energy to realize relay forwarding. There are two main paradigms to implement wireless powered relay (WPR): time switching (TS) relay and power splitting (PS) relay [14]. In TS, one source node transmits RF signals to the relay node in β (TS factor) period for EH and the remaining 1- β period for information reception and forwarding. In PS, a portion δ (PS factor) of signal power received from one source node is used for EH, while the rest 1- δ is used for information decoding. The literature [15,16] investigates PS-based NOMA communication, while the achievable rate and OP in IoT relay NOMA system are analyzed under TS architecture [12,17].
Theoretically, OMA may be more suitable than NOMA in downlink communication, as it can gain higher system throughput [18]. Moreover, it is difficult for a tiny battery-operated IoT node to realize multi-user detection (MUD) in downlink NOMA decoding [19]. To the best of our knowledge, most of the literature only focuses on power allocation in downlink communication from BS to IoT nodes. However, in practice, it is more common in IoT applications, e.g., smart home, smart farming, and intelligent transportation that the IoT nodes, e.g., sensors or cameras, transmit their gathered data to BS or sink node via an uplink channel. Therefore, uplink NOMA is more worthy of investigation in wireless IoT networks. In the related works about uplink NOMA, the authors in [20] only plot energy efficiency in a direct downlink and uplink NOMA system, while the authors in [21] exploit relay uplink NOMA without taking wireless EH into account. Considering wireless EH in relay cooperative networks is a hot research topic, our paper focuses on a wireless powered relay system with uplink NOMA in IoT networks to improve spectrum efficiency and reduce energy consumption, which has never been investigated as far as we know. Different from the traditional EH scheme where the relay harvests energy only from one source or one dedicated power supply station, in our uplink NOMA model, the relay node can harvest RF power from two transmitting nodes due to power multiplexing transmission based on NOMA. Meanwhile, proper relay selection scheme can also reduce relay signal attenuation to improve the performance of the far-end node. In addition, in terms of selecting PS or TS scheme, PS is more complicated and expensive than TS for it uses a power splitter hardware component [22], which is not suitable for low-cost IoT nodes. Additionally, compared with OMA, the PS factor δ is more difficult to optimize due to complex superimposed signals in uplink NOMA. Therefore, TS scheme is employed in our wireless powered cooperative relay model with uplink NOMA.
In this paper, we aim to investigate system performance under different parameters, e.g., TS factor, power allocation factor, the position of relay node, and total transmitting power. Moreover, we also prove our proposed EH scheme achieves superior performance compared with the traditional EH scheme and OMA by numerical results.
The contributions of this paper are as follows:
  • We investigate a relay based uplink NOMA model where the relay can harvest energy from two simultaneously transmitting nodes, i.e., far-end node and near-end node, thanks to the nature of uplink NOMA. As far as we know, this model has never been exploited. Based on our proposed model, the relay close to near-end node not only harvests more energy from near-end node via power allocation factor but also forwards the signal via a shorter path between relay and destination, which reduces path attenuation greatly and achieves higher energy and spectrum efficiency. It is unfeasible in the traditional EH scheme because the relay closer to the source has to be chosen in order to harvest higher forwarding power regardless of long distance relay signal fading. Furthermore, it is impossible to deploy one dedicated power supply station to improve energy efficiency because the station is fixed, high-cost, and low energy-efficient for the far-end node in massive IoT networks;
  • We analyze EH from the superimposed signal and NOMA decoding at different cases in detail. Meanwhile, the closed-form expressions of outage probability and overall system throughput at different communication phases are derived for further numerical simulations;
  • Our numerical analysis provides practical insights into the impacts under different system parameters. Superior system performance is achieved by optimizing the position of relay, TS factor, power allocation factor, etc. Additionally, simulation results also prove that NOMA with our proposed EH scheme outperforms that with the traditional EH scheme and OMA.
The rest of this paper is organized as follows. Section 2 introduces the system model of uplink NOMA. Closed-form expressions about outage probability and system throughput are derived in Section 3. Section 4 conducts numerical simulations and results analysis. Finally, Section 5 concludes the paper.

2. System Model

2.1. Model Introduction

In our considered system model (see Figure 1), node A and node B are IoT nodes or sensors, which both intend to transmit their messages to the BS D with uplink NOMA to improve total spectrum efficiency. Due to long distance and deep fading between A and D, we assume that there is no direct link between them, and the relay R in the midst of A and B assists A to deliver the message to D. Considering that R is an energy-constrained node, and relay forwarding will consume additional energy expenditure, a TS relay protocol [14] is employed at R to scavenge energy from A and B.
We define d i and h i , i { 1 , 2 , 3 , 4 } , as distance and channel coefficient between two nodes, respectively, shown in Figure 1. We assume h i CN(0, 1) follows small scale Rayleigh fading. In addition, we also assume each node in our communication model is equipped with one single antenna.
There are three phases to perform TS relay uplink NOMA transmission. For simplicity, the whole transmission duration is set as a unity. Moreover, a total transmitting power for this uplink NOMA transmission is given as p, hence node A and node B at different communication phases are allocated different transmitting power under a total power constraint. It is noted that the total power constraint is a critical criterion. Firstly, OMA allocates all power for one user at a certain resource block, while NOMA requires a user pair to implement power domain multiplexing. As a whole of nodes A and B, it is rational to limit the total power sum. Secondly, judging system performance between NOMA and OMA requires a same total power constraint. Lastly, in a cluster where multiple users share the same frequency band, total power constraint is necessary to reduce inter-cluster interference.
Phase 1: During the beginning of β fraction transmission period, A sends a non-information signal x a to R for wireless EH. Meanwhile, B transmits its own information-bearing signal x b to D. Both of them use their respective power allocation factors ( w 1 for A, w 2 for B and w 1 + w 2 = 1 ). Thus, relay node R will receive a combined signal for EH. Additionally, due to no direct link between A and D, D can only receive and decode the signal from B.
Phase 2: A sends an information-bearing signal x a to R with w 3 ratio of the given power p during the ( 1 β ) / 2 fraction of a unity transmission period, and B concurrently transmits its own signal x b to D with w 4 ratio of power p. Similarly, we have w 3 + w 4 = 1 . After receiving the combined signal, the relay R uses the NOMA scheme to decode the signal of A, while D decodes the signal of B regularly.
Phase 3: At the rest ( 1 β ) / 2 period, R forwards the signal decoded from A with the energy harvested at phase 1, and B continues transmitting its own message to D. At this phase, B uses the total power p. After receiving the superimposed signal from B and R, D first decodes the stronger signal from B due to the shorter distance between B and D treating the signal from R as interference, and then decodes the signal from R after removing the signal of B in the superimposed signal.

2.2. NOMA Communication Process Analysis

2.2.1. Phase 1

At phase 1, the signal received at R is given by
y r 1 = h 1 w 1 p x a d 1 α + 1 + h 2 w 2 p x b d 2 α + 1 + n r
where α represents the path loss exponent, and n r denotes additive white Gaussian noise (AWGN) at R with variance σ 2 . Furthermore, for simplifying our analysis, we assume the variances of noise in all nodes are equal. It is worth noting that we use the bounded path loss model in (1) to ensure that the path loss is always larger than one for any given distance [23], which is suitable for short transmission communication or dense urban scenarios in our considered IoT network.
By defining c i = d i α + 1 , i { 1 , 2 , 3 , 4 } , (1) is rewritten by
y r 1 = h 1 w 1 p x a c 1 + h 2 w 2 p x b c 2 + n r .
The energy harvested from two signals at R is written by
E r = η β ( | h 1 | 2 w 1 p c 1 + | h 2 | 2 w 2 p c 2 )
where η is energy conversion coefficient, β is TS factor in TS relay protocol.
By defining X = | h 1 | 2 c 1 , Y = | h 2 | 2 c 2 , and V = w 1 X + w 2 Y , we have
E r = η β p V .
The destination D can only receive the signal x b at this phase, which is written by
y d 1 = h 4 w 2 p x b d 4 α + 1 + n d = h 4 w 2 p x b c 4 + n d
where n d denotes AWGN at D with variance σ 2 .
Therefore, the SNR of x b at D is
γ b d 1 = w 2 p W σ 2 = w 2 ρ W
where ρ = p σ 2 and W = | h 4 | 2 c 4 .

2.2.2. Phase 2

At phase 2, R receives two signals from A and B, thus the superimposed signal is given by
y r 2 = h 1 d 1 α + 1 w 3 p x a + h 2 d 2 α + 1 w 4 p x b + n r = h 1 c 1 w 3 p x a + h 2 c 2 w 4 p x b + n r .
In order to investigate NOMA decoding for the combined signal at R, there are two cases that need to be considered, case 1: X > Y , i.e., | h 1 | 2 c 1 > | h 2 | 2 c 2 and case 2: X < Y , i.e., | h 1 | 2 c 1 < | h 2 | 2 c 2 .
It is worth noting that since X and Y are two independent variables, it is difficult to distinguish the two cases. We expect X > Y if c 1 < c 2 or d 1 < d 2 , and X < Y if c 1 > c 2 or d 1 > d 2 . Though d 1 < d 2 cannot guarantee X > Y , and vice versa, it is a simple and efficient scheme to adopt the assumption under statistical channel state information [24].
Case 1: X > Y ( d 1 < d 2 ), R decodes x a due to its stronger signal by treating x b as interference, so the SNR of x a is given by
γ a r 2 = w 3 p X w 4 p Y + σ 2 = w 3 ρ X w 4 ρ Y + 1 .
Case 2: X < Y ( d 1 > d 2 ), R first decodes x b due to its stronger signal and treats x a as interference. The signal x a is obtained only after x b is successfully decoded and removed from the combined signal.
Hence, the signal to interference plus noise ratio (SINR) of x b is written by
γ b r 2 = w 4 ρ Y w 3 ρ X + 1 .
The SNR of x a after removing the strong signal x b can be expressed by
γ a r 2 = w 3 ρ X .
Additionally, the signal of B received at D at this phase is given by
y d 2 = h 4 c 4 w 4 p x b + n d .
Hence, the SNR of x b at D is
γ b d 2 = w 4 p W σ 2 = w 4 ρ W .

2.2.3. Phase 3

At phase 3, R uses the energy harvested at phase 1 to forward the signal x a to D. The transmitting power is given by
P r = E r ( 1 β ) / 2 = 2 η β p V 1 β .
At the same time, B continues to transmit its signal to D. The combined signal at D is
y d 3 = h 3 d 3 α + 1 P r x a + h 4 d 4 α + 1 p x b + n d = h 3 c 3 2 η β p V / ( 1 β ) x a + h 4 c 4 p x b + n d .
We assume the signal from B is stronger than from R due to the shorter distance between B and D, the SNR of x b at D is given by
γ b d 3 = p W P r Z + σ 2 = ρ ( 1 β ) W 2 η β ρ V Z + ( 1 β )
where Z = | h 3 | 2 d 3 α + 1 = | h 3 | 2 c 3 .
The SNR of the signal from R is
γ r d 3 = P r Z σ 2 = 2 η β ρ V Z 1 β .
Based on the expressions of SNR above, we can further analyze the system performance of our proposed model.

3. Performance Analysis

3.1. Outage Probability of Link B->D

We define r a and r b as the required target rate (RTR) of the link A->D via R and B->D, respectively, thus, if the effective transmission rate of A->D or B->D is less than the given r a or r b , link outage occurs. For instance, at phase 1, to avoid an outage from B to D, the following requirement should be satisfied,
log 2 ( 1 + γ b d 1 ) r b .
Namely,
γ b d 1 e b
where e b = 2 r b 1 .
Hence, we have the outage probability (OP) of B->D at phase 1
P b d 1 = Pr { γ b d 1 < e b } .
Substituting (6) into (19), we have
P b d 1 = Pr { W < e b w 2 ρ } .
Due to h 4 CN(0, 1) and W = | h 4 | 2 c 4 , the path gain W is exponentially distributed random variable, i.e., W E( 1 c 4 ). Similarly, we have X E( 1 c 1 ), Y E( 1 c 2 ) and Z E( 1 c 3 ). Accordingly, the probability density functions (PDFs) of X, Y, Z, and W are given by f X ( x ) = c 1 e c 1 x , f Y ( y ) = c 2 e c 2 y , f Z ( z ) = c 3 e c 3 z and f W ( w ) = c 4 e c 4 w respectively.
Based on the PDF of W, (20) can be further expressed by
P b d 1 = 0 e b w 2 ρ f W ( w ) d w = 1 e e b c 4 w 2 ρ .
At phase 2, the SNR of B should satisfy the requirement γ b d 2 e b , therefore, the OP of B->D is written by
P b d 2 = Pr { γ b d 2 < e b } .
By substituting (12) into (22), and some manipulations as (20) and (21), we have
P b d 2 = 1 e e b c 4 w 4 ρ .
At phase 3, the OP of B->D is given by
P b d 3 = Pr { γ b d 3 < e b } .
By substituting (15) into (24), and utilizing the PDFs of W and Z,
P b d 3 = Pr { ρ ( 1 β ) W 2 η β ρ V Z + ( 1 β ) < e b } = 1 c 3 e e b c 4 ρ 0 1 β 2 e b η β c 4 v + ( 1 β ) c 3 f V ( v ) d v
where f V ( v ) is the PDF of V = w 1 X + w 2 Y .
We have the following Theorem 1.
Theorem 1.
The PDF of V is expressed as follows.
f V ( v ) = c 1 c 2 c 2 w 1 c 1 w 2 ( e c 1 v w 1 e c 2 v w 2 ) i f c 2 w 1 c 1 w 2 0 c 1 c 2 v w 1 w 2 e c 2 v w 2 i f c 2 w 1 c 1 w 2 = 0
Proof. 
See Appendix A. □
According to Theorem 1, at c 2 w 1 c 1 w 2 0 , (25) can be further calculated as
P b d 3 = 1 c 1 c 2 ( c 2 w 1 c 1 w 2 ) e e b c 4 ρ × 0 1 2 e b η β c 4 v / ( 1 β ) c 3 + 1 ( e c 1 w 1 v e c 2 w 2 v ) d v [ I ] .
Define t = 2 e b η β c 4 v ( 1 β ) c 3 + 1 ,
[ I ] = ( 1 β ) c 3 2 e b η β c 4 × e ( 1 β ) c 1 c 3 2 e b η β w 1 c 4 1 1 t e ( 1 β ) c 1 c 3 2 e b η β w 1 c 4 t d t ( 1 β ) c 3 2 e b η β c 4 × e ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 1 1 t e ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 t d t = ( 1 β ) c 3 2 e b η β c 4 × e ( 1 β ) c 1 c 3 2 e b η β w 1 c 4 E 1 ( ( 1 β ) c 1 c 3 2 e b η β w 1 c 4 ) ( 1 β ) c 3 2 e b η β c 4 × e ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 E 1 ( ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 )
where E 1 ( x ) = 1 e x t t d t is the exponential integral.
Substituting (27) into (26), at c 2 w 1 c 1 w 2 0
P b d 3 = 1 ( 1 β ) c 1 c 2 c 3 e e b c 4 ρ 2 e b η β c 4 ( c 2 w 1 c 1 w 2 ) × [ e ( 1 β ) c 1 c 3 2 e b η β w 1 c 4 E 1 ( ( 1 β ) c 1 c 3 2 e b η β w 1 c 4 ) e ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 E 1 ( ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 ) ] .
Similarly, at c 2 w 1 c 1 w 2 = 0 ,
P b d 3 = 1 ( 1 β ) c 1 c 2 c 3 2 e b η β c 4 w 1 w 2 e e b c 4 ρ × 0 v v + ( 1 β ) c 3 / 2 e b η β c 4 e c 2 v w 2 d v .
According to 0 x x + β e μ x = β e μ β E i μ β   +   1 μ ,
P b d 3 = 1 ( 1 β ) c 1 c 2 c 3 e e b c 4 ρ 2 e b η β c 4 w 1 w 2 × [ ( 1 β ) c 3 e ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 2 e b η β c 4 E i ( ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 ) + w 2 c 2 ]
where E i · denotes the exponential integral function [25].
To summarize, the OP of B->D at phase 3 P b d 3 is given by (31).
P b d 3 = 1 ( 1 β ) c 1 c 2 c 3 e e b c 4 ρ 2 e b η β c 4 ( c 2 w 1 c 1 w 2 ) × [ e ( 1 β ) c 1 c 3 2 e b η β w 1 c 4 E 1 ( ( 1 β ) c 1 c 3 2 e b η β w 1 c 4 ) e ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 E 1 ( ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 ) ] i f c 2 w 1 c 1 w 2 0 1 ( 1 β ) c 1 c 2 c 3 2 e b η β c 4 w 1 w 2 e e b c 4 ρ [ ( 1 β ) c 3 2 e b η β c 4 e ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 E i ( ( 1 β ) c 2 c 3 2 e b η β w 2 c 4 ) + w 2 c 2 ] i f c 2 w 1 c 1 w 2 = 0

3.2. Outage Probability of Link A->D via R

In this subsection, we plot the OP of A->D via R. Due to relay communication, the OP of A->R and R->D should be derived first before we obtain the OP of the whole link A->D.

3.2.1. Outage Probability of Link A->R

In order to obtain the OP of A->R, we consider two cases according to the position of R, i.e., d 1 < d 2 and d 1 > d 2 .
Case 1: when X > Y   ( d 1 < d 2 ) , we have the OP of A->R
P a r 1 = Pr { γ a r 2 < e a }
where e a = 2 2 r a / ( 1 β ) 1 .
By substituting (8) into (32),
P a r 1 = Pr { X < e a w 4 w 3 Y + e a w 3 ρ } .
Considering X > Y and X < e a w 4 w 3 Y + e a w 3 ρ , if w 3 e a w 4 > 0 ,
P a r 1 = 0 e a ρ ( w 3 e a 1 w 4 ) y e a w 4 w 3 y + e a w 3 ρ f X ( x ) d x f Y ( y ) d y .
Apply the PDFs of X and Y, and some integral operations,
P a r 1 = 0 e a ρ ( w 3 e a w 4 ) ( e c 1 y e c 1 e a w 4 w 3 y c 1 e a w 3 ρ ) c 2 e c 2 y d y = c 2 c 1 + c 2 w 3 c 2 e c 1 e a w 3 ρ e a w 4 c 1 + w 3 c 2 + c 1 c 2 ( w 3 e a w 4 ) e e a ( c 1 + c 2 ) ρ ( w 3 e a w 4 ) ( e a w 4 c 1 + w 3 c 2 ) ( c 1 + c 2 ) .
If w 3 e a w 4 < 0 , with the similar approaches, we have
P a r 1 = c 2 c 1 + c 2 w 3 c 2 e a w 4 c 1 + w 3 c 2 e e a c 1 w 3 ρ
Case 2: when X < Y , x a can be decoded only after x b is decoded correctly, so we should derive the OP of x b first, and then obtain the OP of x a . The OP of x b is expressed by
P b r 2 = Pr { γ b r 2 < e b } = Pr { Y < e b w 3 X w 4 + e b w 4 ρ } .
Considering X < Y and Y < e b w 3 X w 4 + e b w 4 ρ , if w 4 e b w 3 > 0 , by some integral operations,
P b r 2 = c 2 c 1 + c 2 + ( e b w 3 w 4 ) c 1 c 2 ( e b w 3 c 2 + w 4 c 1 ) ( c 1 + c 2 ) e e b ( c 1 + c 2 ) ρ ( w 4 e b w 3 ) .
If w 4 e b w 3 < 0 ,
P b r 2 = c 2 c 1 + c 2 .
The OP of B->R at X < Y can be summarized in (40),
P b r 2 = c 2 c 1 + c 2 i f w 4 e b w 3 < 0 c 2 c 1 + c 2 + ( e b w 3 w 4 ) c 1 c 2 e e b ( c 1 + c 2 ) ρ ( w 4 e b w 3 ) ( e b w 3 c 2 + w 4 c 1 ) ( c 1 + c 2 ) i f w 4 e b w 3 > 0
Based on P b r 2 , the OP of x a is given by
P a r 2 = 1 ( 1 Pr { γ a r 2 < e a } ) × ( 1 P b r 2 )
Apply the result from (10), (41) is expressed as
P a r 2 = 1 e e a c 1 w 3 ρ × ( 1 P b r 2 )
where P b r 2 is shown in (40).
To summarize, P a r is shown in (43).
P a r = c 2 c 1 + c 2 w 3 c 2 e a w 4 c 1 + w 3 c 2 e c 1 e a w 3 ρ + c 1 c 2 ( w 3 e a w 4 ) ( e a w 4 c 1 + w 3 c 2 ) ( c 1 + c 2 ) e e a ( c 1 + c 2 ) ρ ( w 3 e a w 4 ) i f w 3 e a w 4 > 0 & X > Y c 2 c 1 + c 2 w 3 c 2 e a w 4 c 1 + w 3 c 2 e e a c 1 w 3 ρ i f w 3 e a w 4 < 0 & X > Y ( 1 e e a c 1 w 3 ρ ) ( c 1 c 1 + c 2 ( e b w 3 w 4 ) c 1 c 2 ( e b w 3 c 2 + w 4 c 1 ) ( c 1 + c 2 ) e e b ( c 1 + c 2 ) ρ ( w 4 e b w 3 ) ) i f w 4 e b w 3 > 0 & X < Y ( 1 e e a c 1 w 3 ρ ) × c 1 c 1 + c 2 i f w 4 e b w 3 < 0 & X < Y

3.2.2. Outage Probability of Link R->D

When successful decoding the signal from A at phase 2, R forwards x a with the power harvested from phase 1. Upon receiving the combined signal, D decodes it using NOMA. Since D should first decode and remove the strong signal x b correctly, therefore, the OP of x a at D is expressed by
P r d = 1 ( 1 Pr { γ r d 3 < e a } [ I I ] ) × ( 1 P b d 3 )
where P b d 3 is shown in (31).
Utilizing (16) and the PDF of Z, [II] in (44) is derived as
[ I I ] = 1 0 e e a ( 1 β ) c 3 2 η β ρ V f V ( v ) d v [ I I I ] .
If c 2 w 1 c 1 w 2 = 0 ,
[ I I I ] = c 1 c 2 w 1 w 2 0 v e e a ( 1 β ) c 3 2 η β ρ v c 2 v w 2 d v
Define s = c 2 v w 2 and b = e a ( 1 β ) c 2 c 3 2 η β ρ w 2 ,
[ I I I ] = w 2 c 1 w 1 c 2 0 s e s e b / s d s = w 2 c 1 w 1 c 2 Γ b ( 2 )
where Γ b ( 2 ) is generalized incomplete gamma function.
According to the equation Γ b ( a ) = 2 b a / 2 K a ( 2 b ) [26],
[ I I I ] = e a ( 1 β ) c 1 c 3 η β ρ w 1 K 2 ( 2 e a ( 1 β ) c 2 c 3 η β ρ w 2 )
where K 2 ( · ) is the second order modified Bessel function of the second kind.
Hence, at c 2 w 1 c 1 w 2 = 0 ,
P r d = 1 ( 1 P b d 3 ) × e a ( 1 β ) c 1 c 3 η β ρ w 1 K 2 ( 2 e a ( 1 β ) c 2 c 3 η β ρ w 2 )
If c 2 w 1 c 1 w 2 0 ,
[ I I I ] = c 1 c 2 c 2 w 1 c 1 w 2 × ( 0 e e a ( 1 β ) c 3 2 η β ρ v c 1 v w 1 d v [ I V ] 0 e e a ( 1 β ) c 3 2 η β ρ v c 2 v w 2 d v [ V ] )
According to the equation 0 e p 4 x q x d x = p q K 1 ( p q ) ,
[ I V ] = w 1 c 1 2 e a ( 1 β ) c 1 c 3 η β ρ w 1 K 1 ( 2 e a ( 1 β ) c 1 c 3 η β ρ w 1 )
and
[ V ] = w 2 c 2 2 e a ( 1 β ) c 2 c 3 η β ρ w 2 K 1 ( 2 e a ( 1 β ) c 2 c 3 η β ρ w 2 )
where K 1 ( · ) is the first order modified Bessel function of the second kind.
Hence, we have at c 2 w 1 c 1 w 2 0 ,
P r d = 1 ( 1 P b d 3 ) c 1 c 2 c 2 w 1 c 1 w 2 × [ w 1 c 1 2 e a ( 1 β ) c 1 c 3 η β ρ w 1 K 1 ( 2 e a ( 1 β ) c 1 c 3 η β ρ w 1 ) w 2 c 2 2 e a ( 1 β ) c 2 c 3 η β ρ w 2 K 1 ( 2 e a ( 1 β ) c 2 c 3 η β ρ w 2 ) ]

3.2.3. Outage Probability of the Whole Link A->D

Based on the OP of link A->R and R->D, it is easy to obtain the OP of link A->D,
P a d = 1 ( 1 P a r ) × ( 1 P r d )
where P a r is shown in (43), and P r d is shown in (49) at c 2 w 1 c 1 w 2 = 0 and (53) at c 2 w 1 c 1 w 2 0 .

3.3. System Throughput

During a unit transmission period, the system throughput from B to D is the sum of throughput at three phases, which is given by
S b d = β r b ( 1 P b d 1 ) + 1 β 2 r b ( 1 P b d 2 ) + 1 β 2 r b ( 1 P b d 3 )
Since the effective transmission period from A to D is ( 1 β ) / 2 , the system throughput of link A->D is given by
S a d = 1 β 2 r a ( 1 P a d )
As a result, the overall system throughput of our uplink NOMA model is expressed by
S N = S b d + S a d = β r b ( 1 P b d 1 ) + 1 β 2 r b ( 1 P b d 2 ) + 1 β 2 r b ( 1 P b d 3 ) + 1 β 2 r a ( 1 P a d )

4. Numerical Results

In this section, numerical results are presented to facilitate understanding the performance of our wireless powered relay model with uplink NOMA. Monte Carlo simulations are performed to verify our analysis. In the numerical simulation, we select two values of distance from A to R (i.e., d 1 = 1.5 m and d 1 = 3 m) to investigate the system performance in two cases, i.e., X > Y and X < Y , which have been discussed theoretically in Section 3. Since there are different OPs at uplink communication phases, in the following figures we only illustrate those OP curves impacted by the parameter discussed in this subsection. We also compare the performance of our proposed TS-based EH scheme, where relay harvests energy from two nodes, with the traditional TS-based EH scheme, where relay harvests energy just from the far-end source node. Furthermore, we also exhibit the differences in system throughput between NOMA and TDMA-based OMA. In order to distinguish OP curves in each phase and system throughput curves in different communication fashions, we make descriptions for each curve in Table 1. In addition, the default values of system parameters and their corresponding descriptions related to our simulation are all listed in Table 2.

4.1. OP and System Throughput vs. Power Allocation Factor w 1 and w 3

In this subsection, we consider the variations on OPs and system throughput when w 1 or w 3 ranges from 0.1 to 0.9.
With the larger w 1 , the near-end node B is allocated less transmitting power, which leads to greater OP of B->D at phase 1 at two cases, i.e., at d 1 = 1.5 shown in Figure 2 and at d 1 = 3 shown in Figure 3. However, two curves P a d and P b d 3 demonstrate the opposite results, where the curve P a d drops at d 1 = 1.5 and rises at d 1 = 3 , while P b d 3 rises at d 1 = 1.5 and drops at d 1 = 3 . When the relay R is close to node A ( d 1 = 1.5 ), R can harvest more power with increasing w 1 , therefore, the higher power relay signal results in the lower P a d . On the contrary, when R is close to B ( d 1 = 3), with rising w 1 , R harvests less power from the adjacent node B, which leads to the higher P a d . The reason for different results of P b d 3 at two cases is that the stronger signal of A means the stronger interference with B and the higher value of P b d 3 , and vice versa. As a result, w 1 should be carefully considered according to the position of relay node to achieve the required system OPs of two nodes.
Compared with the traditional EH scheme where relay harvests the energy only from a source, our EH scheme results in a higher-power relay signal for the far-end node, hence, the OP of A->D and total system throughput have been greatly improved at the expense of slight degradation for the OP of B->D, which is also illustrated in Figure 4.
From Figure 4, it is obvious that S N outperforms S N s greatly at different w 1 and different positions of the relay node, illustrating the benefits from our EH scheme. Additionally, with our EH scheme the system throughput at d 1 = 3 is much larger than that at d 1 = 1.5 . When R is close to B, R still obtains plenty of energy if more energy is allocated to B at smaller w 1 . The most important reason is that path fading of link R->D at d 1 = 3 is much smaller than at d 1 = 1.5 , leading to lower OP and higher throughput. Therefore, the maximum value on the curve S N ( d 1 = 3 ) occurs at w 1 = 0.3 representing the 70% of total power for node B. Similarly, the maximum value on the curve S N ( d 1 = 1.5 ) occurs at w 1 = 0.7 . We also observe that the lowest value on the curve S N ( d 1 = 1.5 ) occurs at w 1 = 0.5 . The reason is that R harvests the least energy from two nodes when total power is allocated equally for two nodes. Additionally, it is difficult for R to decode the signal A from the combined signal with NOMA when the SNRs of the signal A and B at R are close, which also reduces the system throughput.
Under the traditional EH scheme, it is observed that system throughput at d 1 = 1.5 outperforms that at d 1 = 3 . Due to energy harvesting only from one source node, relay R closer to the source node can obtain more power even though a longer distance from R to D causes deeper path loss. In general, our proposed EH scheme can achieve the highest throughput, where the relay node close to the near-end node not only harvest more energy but also leads to less path fading from R to D.
The parameter w 3 is the power allocation factor at phase 2, so it has no impact on the OPs of B->D at phase 1 and phase 3, which are not shown in Figure 5. With w 3 increasing, the OPs of B->D at phase 2 at two cases rise identically due to less power assigned for B and two curves P b d 2 ( d 1 = 1.5 ) and P b d 2 ( d 1 = 3 ) coincide with each other. In terms of OPs of links A->R and A->D at d 1 = 1.5 , the stronger transmitting signal from A leads to higher successful decoding rate at R and gains much lower OP of A->R. Hence, the curve P a r ( d 1 = 1.5 ) descends rapidly with larger w 3 . However, the curves P a d ( d 1 = 1.5 ) declines slowly, especially at bigger w 3 . The reason is that the OP of link R->D, which depends on the energy harvested at phase 1, is much larger than the OP of A->R, indicating that the OP of link A->D is determined by the OP of link R->D mainly. Additionally, at d 1 = 3 , it is observed that the OP of A->R declines and obtains a minimum value at w 3 = 0.2 but rise with larger w 3 . When R is closer to B, R decodes the strong signal from B first, and then decodes the weak signal A. Hence, if w 3 is bigger, the stronger transmitting signal from A will prevent the decoding the signal B at R, and finally result in higher OP of A->R. Moreover, if w 3 is smaller, weak signal of A also causes higher OP of A->R. Similarly, the OP of A->D presents the same result as the OP of A->R.
In Figure 6, the system throughput curves at different cases are shown. The performance with the traditional EH scheme is clearly poor and the curves of S N s at two relay positions are basically stable and far below the curves of S N due to lower energy harvested. Furthermore, The reason of performance degradation at the left side of the curves S N ( d 1 = 1.5 ) and the right side of the curves S N ( d 1 = 3 ) lies in higher OPs shown in Figure 5 discussed above. In general, it is obvious that S N outperforms S N s at different relay positions if a proper w 3 value is given. Moreover, the best system throughput with our EH scheme at d 1 = 3 is achieved mainly due to the lower path fading between R and D.

4.2. OP and System Throughput vs. Time Switching Factor β

Figure 7 shows the relationship between OP and β at d 1 = 3 . We do not show the case at d 1 = 1.5 because the two cases both reflect the approximate results. With larger β , the relay signal with more harvested energy from R will interfere with the transmission of B->D at phase 3. That is why the curve P b d 3 increases gradually. Additionally, the curve P a r is also observed to rise slightly. The reason is that larger β will reduce effective transmission time. Thus, at the same RTR and transmitting power for A by w 3 , shorter transmission time results in the higher OP of A->R. Moreover, although shorter effective transmission time may also increase the OP of R->D, more energy harvested at R decreases the OP of R->D. Therefore, the curves P r d and P a d decline gradually at the given reasonable range of β from 0.02 to 0.2. Additionally, the reason that P a d and P b d s outperform P a d s and P b d 3 , respectively, lies in different EH schemes.
Figure 8 shows the relationship between system throughput and β . The curve S N ( d 1 = 3 ) illustrates the best system performance among all curves. When β is smaller, the energy harvested by R is lower, and when β is greater, the effective transmission duration becomes shorter. These two cases both lead to the lower system throughput. Therefore, a maximum value in each curve exists at a certain β . The curves S N ( d 1 = 3 ) and S N ( d 1 = 1.5 ) gain their maximum values at β = 0.04 and β = 0.08, respectively. In addition, the reason for a larger β at d 1 = 1.5 than d 1 = 3 is due to more energy required against long distance fading.

4.3. OP and System Throughput vs. Distance d 1 between A and R

Figure 9 plots the OP performance with different positions of relay. It is noted that we select d 1 = 3.9 m instead of d 1 = 4 m because it is impossible for R to stay with B. Furthermore, NOMA cannot successfully decode two combined signals when they are similar path loss gain, i.e., d 1 = 2 , so we remove the points d 1 = 2 at the curves. Since we have discussed two cases d 1 < d 2 and d 1 > d 2 in Section 3.2, accordingly, we select two values d 1 = 1.9 representing case 1 ( d 1 < d 2 ) and d 1 = 2.1 representing case 2 ( d 1 > d 2 ) for investigation in this simulation.
It is observed that the curves of P b d 3 and P a d are limited within the range of curves of P b d s and P a d s , which explains that with our EH scheme the OP of the far-end node is improved at the cost of trivial performance degradation for the near-end node. When relay R is closer to A, higher OP of A->D is gained due to deeper path fading from R to D, while the OP of B->D is very lower due to minor interference from R. When R approaches B, the link A->D gains lower OP due to smaller fading of link R->D, as well as higher power harvested from B correspondingly. The curve OP of B->D becomes larger due to stronger interference from R. It is noted that when R is in the middle of two nodes A and B, i.e., at d 1 = 2.1 , due to approximate signal strength from A and B, NOMA decoding at R will cause higher OP. Furthermore, when R is getting closer to B, e.g., d 1 > 3.5 , the OP of A->D becomes bigger again. That is because the higher OP of B->D prevents the decoding of signal A.
Figure 10 shows the curves of system throughput at three cases. The first is NOMA with our proposed EH scheme, the second is NOMA with the traditional EH scheme, and the last is TDMA-based OMA with the traditional EH scheme. Obviously, NOMA at two cases achieves about twice system throughput compared with regular OMA, demonstrating a huge advantage of NOMA over OMA. In addition, our EH scheme outperforms the traditional EH scheme, and the maximum value is obtained at S N at d 1 = 3.5 , which indicates the position of relay plays a significant role on overall system throughput.

4.4. OP and System Throughput vs. Total Transmitting Power p

Figure 11 plots OPs under different transmitting power from 10 dB to 50 dB. All the OP curves definitely decrease given the larger transmitting power. At phase 1 and phase 2, the near-end node B can communicate with D regularly without being affected by the far-end node A, therefore, the curves P b d 1 and P b d 2 both decline rapidly (two curves are overlapped in the figure). However, at larger power given, the curves P b d 3 and P b d s tend to decline slowly. That is because greater power not only increases the SNR of B but also increases the SNR of A, which is interference with D to decode combined signal of A and B. One can infer that more power supply may not improve the OP of near-end node efficiently in NOMA. The conclusion can also be drawn from Figure 12. In this figure, the best performance of NOMA can be observed at d 1 = 3 under our proposed EH scheme, but with the larger power supply, all the curves except OMA converge gradually, showing that under different EH schemes and different relay positions NOMA achieves an approximate system throughput. The results also indicate NOMA plays the crucial role with our proposed EH scheme at the constrained power supply. Moreover, NOMA still demonstrates excellent system performance against OMA. One can infer that in energy-constrained IoT uplink networks, NOMA is a better selection than OMA.

5. Conclusions

In this paper, we investigate wireless powered relay with uplink NOMA in IoT networks to improve system performance, as well as reduce energy consumption. The closed-form expressions of OPs of uplink communication and overall system throughput are derived. Our numerical simulation demonstrates superior system performance by optimizing important parameters, e.g., TS factor, power allocation factor, the position of relay node at the constrained total transmitting power compared to the traditional EH scheme and OMA. The most important conclusion is that with our EH scheme, by optimizing relay position more energy is harvested from the near-end node at the relay and long distance relay signal attenuation from the relay to the destination is greatly reduced as well, which achieves the best system performance. Additionally, although TS relay communication reduces the effective transmission time compared with PS relay communication, adequate energy harvested from two transmitting nodes in our NOMA model can offset the effect and gain better performance.

Author Contributions

Z.L. wrote the manuscript. S.C. partially revised the paper. J.L. reviewed the paper prior to submission. Z.L. provided the concept and simulation design of the study. All authors discussed the results, analyzed the data and commented on the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Plan of Fujian Province, China in 2020 under Grant 2020J01933, and the Faculty Research Grants of Fujian Jiangxia University, Fuzhou, China under Grant JXZ2019008.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Theorem 1.
Here we derive the PDF of V = V 1 + V 2 where V 1 = w 1 X and V 2 = w 2 Y .
Due to X E   ( λ 1 = 1 c 1 ) and Y E   ( λ 2 = 1 c 2 ) , we have V 1 = w 1 X E   ( w 1 c 1 ) and V 2 = w 2 X E   ( w 2 c 2 ) .
Since V 1 and V 2 are independent, the PDF of V is calculated by
f V ( v ) = f V 1 ( v 1 ) f V 2 ( v v 1 ) d v 1 = c 1 c 2 w 1 w 2 e c 2 v w 2 0 v e c 2 w 1 c 1 w 2 w 1 w 2 v 1 d v 1
If c 2 w 1 c 1 w 2 = 0 , we have
f V ( v ) = c 1 c 2 v w 1 w 2 e c 2 v w 2
else
f V ( v ) = c 1 c 2 c 2 w 1 c 1 w 2 e c 2 v w 2 ( e c 2 w 1 v c 1 w 2 v w 1 w 2 1 ) = c 1 c 2 c 2 w 1 c 1 w 2 ( e c 1 v w 1 e c 2 v w 2 )
This ends the proof. □

References

  1. Gong, J.; Chen, X. Achievable rate region of non-orthogonal multiple access systems with wireless powered decoder. IEEE J. Sel. Areas Commun. 2017, 35, 2846–2859. [Google Scholar] [CrossRef] [Green Version]
  2. Ding, Z.; Lei, X.; Karagiannidis, G.K.; Schober, R.; Yuan, J.; Bhargava, V.K. A survey on non-orthogonal multiple access for 5g networks: Research challenges and future trends. IEEE J. Sel. Areas Commun. 2017, 35, 2181–2195. [Google Scholar] [CrossRef] [Green Version]
  3. Yang, Z.; Ding, Z.; Fan, P.; Al-Dhahir, N. A general power allocation scheme to guarantee quality of service in downlink and uplink noma systems. IEEE Trans. Wirel. Commun. 2016, 15, 7244–7257. [Google Scholar] [CrossRef]
  4. Ding, Z.; Yang, Z.; Fan, P.; Poor, H.V. On the performance of non-orthogonal multiple access in 5g systems with randomly deployed users. IEEE Signal Process. Lett. 2014, 21, 1501–1505. [Google Scholar] [CrossRef] [Green Version]
  5. Ding, Z.; Fan, P.; Poor, H.V. Impact of user pairing on 5g nonorthogonal multiple-access downlink transmissions. IEEE Trans. Veh. Technol. 2015, 65, 6010–6023. [Google Scholar] [CrossRef]
  6. Ibrahim, A.S.; Sadek, A.K.; Su, W.; Liu, K.J.R. Cooperative communications with relay-selection: When to cooperate and whom to cooperate with? IEEE Trans. Wirel. Commun. 2008, 7, 2814–2827. [Google Scholar] [CrossRef] [Green Version]
  7. Razeghi, B.; Hodtani, G.A.; Nikazad, T. Multiple criteria relay selection scheme in cooperative communication networks. Wirel. Pers. Commun. 2017, 96, 2539–2561. [Google Scholar] [CrossRef]
  8. Xu, M.; Ji, F.; Wen, M.; Duan, W. Novel receiver design for the cooperative relaying system with non-orthogonal multiple access. IEEE Commun. Lett. 2016, 20, 1679–1682. [Google Scholar] [CrossRef]
  9. Kim, J.-B.; Lee, I.-H. Capacity analysis of cooperative relaying systems using non-orthogonal multiple access. IEEE Commun. Lett. 2015, 19, 1949–1952. [Google Scholar] [CrossRef]
  10. Yang, L.; Jiang, H.; Ye, Q.; Ding, Z.; Fang, F.; Shi, J.; Chen, J.; Xue, X. Opportunistic adaptive non-orthogonal multiple access in multiuser wireless systems: Probabilistic user scheduling and performance analysis. IEEE Trans. Wirel. Commun. 2020, 19, 6065–6082. [Google Scholar] [CrossRef]
  11. Janghel, K.; Prakriya, S. Performance of adaptive oma/cooperative-noma scheme with user selection. IEEE Commun. Lett. 2018, 22, 2092–2095. [Google Scholar] [CrossRef]
  12. Rauniyar, A.; Engelstad, P.E.; Østerbø, O.N. Performance analysis of rf energy harvesting and information transmission based on noma with interfering signal for iot relay systems. IEEE Sens. J. 2019, 19, 7668–7682. [Google Scholar] [CrossRef] [Green Version]
  13. Lu, X.; Wang, P.; Niyato, D.; Kim, D.I.; Han, Z. Wireless networks with rf energy harvesting: A contemporary survey. IEEE Commun. Surv. Tutor. 2014, 17, 757–789. [Google Scholar] [CrossRef] [Green Version]
  14. Nasir, A.A.; Zhou, X.; Durrani, S.; Kennedy, R.A. Relaying protocols for wireless energy harvesting and information processing. IEEE Trans. Wirel. Commun. 2013, 12, 3622–3636. [Google Scholar] [CrossRef] [Green Version]
  15. Liu, Y.; Ding, Z.; Elkashlan, M.; Poor, H.V. Cooperative non-orthogonal multiple access with simultaneous wireless information and power transfer. IEEE J. Sel. Areas Commun. 2016, 34, 938–953. [Google Scholar] [CrossRef] [Green Version]
  16. Do, N.T.; Costa, D.B.D.; Duong, T.Q.; An, B. A bnbf user selection scheme for noma-based cooperative relaying systems with swipt. IEEE Commun. Lett. 2016, 21, 664–667. [Google Scholar] [CrossRef] [Green Version]
  17. Mondal, S.; Roy, S.D.; Kundu, S. Outage analysis for noma-based energy harvesting relay network with imperfect csi and transmit antenna selection. IET Commun. 2020, 14, 2240–2249. [Google Scholar] [CrossRef]
  18. Jang, J.; Lee, K.B. Transmit power adaptation for multiuser ofdm systems. IEEE J. Sel. Areas Commun. 2003, 21, 171–178. [Google Scholar] [CrossRef] [Green Version]
  19. Al-Imari, M.; Xiao, P.; Imran, M.A.; Tafazolli, R. Uplink non-orthogonal multiple access for 5g wireless networks. In Proceedings of the 2014 11th International Symposium on Wireless Communications Systems (ISWCS), Barcelona, Spain, 26–29 August 2014; pp. 781–785. [Google Scholar]
  20. Zamani, M.; Eslami, M.; Khorramizadeh, M.; Zamani, H.; Ding, Z. Optimizing weighted-sum energy efficiency in downlink and uplink noma systems. IEEE Trans. Veh. Technol. 2020, 69, 11112–11127. [Google Scholar] [CrossRef]
  21. Kader, M.F.; Shin, S.Y. Coordinated direct and relay transmission using uplink noma. IEEE Wirel. Commun. Lett. 2017, 7, 400–403. [Google Scholar] [CrossRef]
  22. Guo, W.; Zhou, S.; Chen, Y.; Wang, S.; Chu, X.; Niu, Z. Simultaneous information and energy flow for iot relay systems with crowd harvesting. IEEE Commun. Mag. 2016, 54, 143–149. [Google Scholar] [CrossRef]
  23. Ding, Z.; Krikidis, I.; Sharif, B.; Poor, H.V. Wireless information and power transfer in cooperative networks with spatially random relays. IEEE Trans. Wirel. Commun. 2014, 13, 4440–4453. [Google Scholar] [CrossRef] [Green Version]
  24. Wei, Z.; Dai, L.; Ng, D.W.K.; Yuan, J. Performance analysis of a hybrid downlink-uplink cooperative noma scheme. In Proceedings of the 2017 IEEE 85th Vehicular Technology Conference (VTC Spring), Sydney, NSW, Australia, 4–7 June 2017; pp. 1–7. [Google Scholar]
  25. Jeffrey, A.; Zwillinger, D. Table of Integrals, Series, and Products; Elsevier: Amsterdam, The Netherlands, 2007. [Google Scholar]
  26. Chaudhry, M.A.; Zubair, S.M. On A Class of Incomplete Gamma Functions with Applications; CRC Press: Boca Raton, FL, USA, 2001. [Google Scholar]
Figure 1. Wireless powered relay system with uplink NOMA.
Figure 1. Wireless powered relay system with uplink NOMA.
Sensors 21 05467 g001
Figure 2. OP vs. w 1 at d 1 = 1.5 .
Figure 2. OP vs. w 1 at d 1 = 1.5 .
Sensors 21 05467 g002
Figure 3. OP vs. w 1 at d 1 = 3 .
Figure 3. OP vs. w 1 at d 1 = 3 .
Sensors 21 05467 g003
Figure 4. System throughput vs. w 1 .
Figure 4. System throughput vs. w 1 .
Sensors 21 05467 g004
Figure 5. OP vs. w 3 .
Figure 5. OP vs. w 3 .
Sensors 21 05467 g005
Figure 6. System throughput vs. w 3 .
Figure 6. System throughput vs. w 3 .
Sensors 21 05467 g006
Figure 7. OP vs. β at d 1 = 3 .
Figure 7. OP vs. β at d 1 = 3 .
Sensors 21 05467 g007
Figure 8. System throughput vs. β at d 1 = 3 .
Figure 8. System throughput vs. β at d 1 = 3 .
Sensors 21 05467 g008
Figure 9. OP vs. d 1 .
Figure 9. OP vs. d 1 .
Sensors 21 05467 g009
Figure 10. System throughput vs. d 1 .
Figure 10. System throughput vs. d 1 .
Sensors 21 05467 g010
Figure 11. OP vs. total transmitting power.
Figure 11. OP vs. total transmitting power.
Sensors 21 05467 g011
Figure 12. System throughput vs. total transmitting power.
Figure 12. System throughput vs. total transmitting power.
Sensors 21 05467 g012
Table 1. Curves Description.
Table 1. Curves Description.
Curves NameDescription
P b d 1 OP of B->D at phase 1
P b d 2 OP of B->D at phase 2
P b d 3 OP of B->D at phase 3 with our proposed EH scheme
P b d s OP of B->D at phase 3 with the traditional EH scheme
P a r OP of A->R with our proposed EH scheme
P r d OP of R->D with our proposed EH scheme
P a d OP of A->D with our proposed EH scheme
P a d s OP of A->D with the traditional EH scheme
S N NOMA System throughput with our proposed EH scheme
S N s NOMA System throughput with the traditional EH scheme
S O s OMA System throughput with the traditional EH scheme
Table 2. Simulation Parameters by Default.
Table 2. Simulation Parameters by Default.
ParametersDescription
η = 1 energy conversion coefficient
α = 3 path loss exponent
d = 5 mdistance between A and D
d 1 = 1.5 or 3 mdistance between A and R
d 4 = 1 mdistance between B and D
d 2 = d d 1 d 4 distance between R and B
d 3 = d d 1 distance between R and D
β = 0.1 stime switching factor
r a = 0.2 bit/s/Hzrequired target rate of A->D
r b = 0.2 bit/s/Hzrequired target rate of B->D
p = 40 dBtotal transmitting power
w 1 = 0.5 power allocation factor for A at phase 1
w 2 = 1 w 1 power allocation factor for B at phase 1
w 3 = 0.5 power allocation factor for A at phase 2
w 4 = 1 w 3 power allocation factor for B at phase 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, Z.; Cao, S.; Li, J. Time Switching Based Wireless Powered Relay Transmission with Uplink NOMA. Sensors 2021, 21, 5467. https://doi.org/10.3390/s21165467

AMA Style

Lin Z, Cao S, Li J. Time Switching Based Wireless Powered Relay Transmission with Uplink NOMA. Sensors. 2021; 21(16):5467. https://doi.org/10.3390/s21165467

Chicago/Turabian Style

Lin, Zhihua, Shihua Cao, and Jianqing Li. 2021. "Time Switching Based Wireless Powered Relay Transmission with Uplink NOMA" Sensors 21, no. 16: 5467. https://doi.org/10.3390/s21165467

APA Style

Lin, Z., Cao, S., & Li, J. (2021). Time Switching Based Wireless Powered Relay Transmission with Uplink NOMA. Sensors, 21(16), 5467. https://doi.org/10.3390/s21165467

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop