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
and
,
, as distance and channel coefficient between two nodes, respectively, shown in
Figure 1. We assume
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 to R for wireless EH. Meanwhile, B transmits its own information-bearing signal to D. Both of them use their respective power allocation factors ( for A, for B and ). 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 to R with ratio of the given power p during the fraction of a unity transmission period, and B concurrently transmits its own signal to D with ratio of power p. Similarly, we have . 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 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
where
represents the path loss exponent, and
denotes additive white Gaussian noise (AWGN) at R with variance
. 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
, (
1) is rewritten by
The energy harvested from two signals at R is written by
where
is energy conversion coefficient,
is TS factor in TS relay protocol.
By defining
,
, and
, we have
The destination D can only receive the signal
at this phase, which is written by
where
denotes AWGN at D with variance
.
Therefore, the SNR of
at D is
where
and
.
2.2.2. Phase 2
At phase 2, R receives two signals from A and B, thus the superimposed signal is given by
In order to investigate NOMA decoding for the combined signal at R, there are two cases that need to be considered, case 1: , i.e., and case 2: , i.e., .
It is worth noting that since
X and
Y are two independent variables, it is difficult to distinguish the two cases. We expect
if
or
, and
if
or
. Though
cannot guarantee
, and vice versa, it is a simple and efficient scheme to adopt the assumption under statistical channel state information [
24].
Case 1:
(
), R decodes
due to its stronger signal by treating
as interference, so the SNR of
is given by
Case 2: (), R first decodes due to its stronger signal and treats as interference. The signal is obtained only after is successfully decoded and removed from the combined signal.
Hence, the signal to interference plus noise ratio (SINR) of
is written by
The SNR of
after removing the strong signal
can be expressed by
Additionally, the signal of B received at D at this phase is given by
Hence, the SNR of
at D is
2.2.3. Phase 3
At phase 3, R uses the energy harvested at phase 1 to forward the signal
to D. The transmitting power is given by
At the same time, B continues to transmit its signal to D. The combined signal at D is
We assume the signal from B is stronger than from R due to the shorter distance between B and D, the SNR of
at D is given by
where
.
The SNR of the signal from R is
Based on the expressions of SNR above, we can further analyze the system performance of our proposed model.
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.,
m and
m) to investigate the system performance in two cases, i.e.,
and
, 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 and
In this subsection, we consider the variations on OPs and system throughput when or ranges from 0.1 to 0.9.
With the larger
, 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
= 1.5 shown in
Figure 2 and at
= 3 shown in
Figure 3. However, two curves
and
demonstrate the opposite results, where the curve
drops at
and rises at
, while
rises at
and drops at
. When the relay R is close to node A (
), R can harvest more power with increasing
, therefore, the higher power relay signal results in the lower
. On the contrary, when R is close to B (
= 3), with rising
, R harvests less power from the adjacent node B, which leads to the higher
. The reason for different results of
at two cases is that the stronger signal of A means the stronger interference with B and the higher value of
, and vice versa. As a result,
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
outperforms
greatly at different
and different positions of the relay node, illustrating the benefits from our EH scheme. Additionally, with our EH scheme the system throughput at
is much larger than that at
. When R is close to B, R still obtains plenty of energy if more energy is allocated to B at smaller
. The most important reason is that path fading of link R->D at
is much smaller than at
, leading to lower OP and higher throughput. Therefore, the maximum value on the curve
(
) occurs at
representing the 70% of total power for node B. Similarly, the maximum value on the curve
(
) occurs at
. We also observe that the lowest value on the curve
(
) occurs at
. 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 outperforms that at . 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
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
increasing, the OPs of B->D at phase 2 at two cases rise identically due to less power assigned for B and two curves
(
) and
(
) coincide with each other. In terms of OPs of links A->R and A->D at
, 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
(
) descends rapidly with larger
. However, the curves
(
) declines slowly, especially at bigger
. 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
, it is observed that the OP of A->R declines and obtains a minimum value at
= 0.2 but rise with larger
. When R is closer to B, R decodes the strong signal from B first, and then decodes the weak signal A. Hence, if
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
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
at two relay positions are basically stable and far below the curves of
due to lower energy harvested. Furthermore, The reason of performance degradation at the left side of the curves
(
) and the right side of the curves
(
) lies in higher OPs shown in
Figure 5 discussed above. In general, it is obvious that
outperforms
at different relay positions if a proper
value is given. Moreover, the best system throughput with our EH scheme at
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
. We do not show the case at
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
increases gradually. Additionally, the curve
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
, 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
and
decline gradually at the given reasonable range of
from 0.02 to 0.2. Additionally, the reason that
and
outperform
and
, respectively, lies in different EH schemes.
Figure 8 shows the relationship between system throughput and
. The curve
(
) 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
(
) and
(
) gain their maximum values at
= 0.04 and
= 0.08, respectively. In addition, the reason for a larger
at
than
is due to more energy required against long distance fading.
4.3. OP and System Throughput vs. Distance between A and R
Figure 9 plots the OP performance with different positions of relay. It is noted that we select
m instead of
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.,
, so we remove the points
at the curves. Since we have discussed two cases
<
and
>
in
Section 3.2, accordingly, we select two values
representing case 1 (
<
) and
representing case 2 (
>
) for investigation in this simulation.
It is observed that the curves of and are limited within the range of curves of and , 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 , 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., , 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
at
, 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
and
both decline rapidly (two curves are overlapped in the figure). However, at larger power given, the curves
and
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
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.