A Comprehensive Analytical Framework under Practical Constraints for a Cooperative NOMA System Empowered by SWIPT IoT

: Over the past decade, there has been notable attention directed towards spectrum and energy efficiency in conjunction with simultaneous wireless information and power transfer (SWIPT), aimed at extending the operational lifespan of energy-constrained wireless devices within cooperative non-orthogonal multiple access (C-NOMA) systems. This article delves into a system model comprising a transmitter, a relay


Introduction 1.Information Background
NOMA represents a multiple access technique that integrates multiplexing and multiuser detection methods, enabling multiple users to share a single channel efficiently.In contrast to Frequency Division Multiple Access (FDMA) and Orthogonal Frequency Division Multiple Access (OFDM), which fall under the category of Orthogonal Multiple Access (OMA) techniques, NOMA presents a more effective means of leveraging the radio spectrum, exemplified by Orthogonal Frequency Division Multiple Access (OFDMA).NOMA facilitates more efficient resource utilization by allocating distinct portions of resources such as frequency, power, or codes to individual users.Moreover, it enables the provision of asymmetric data rates and diverse Quality of Service (QoS) standards, thus enhancing flexibility in service offerings [1,2].The Time-Splitting Relaying (TSR) protocol has demonstrated notable efficacy in enhancing the performance of NOMA-Simultaneous Wireless Information and Power Transfer (SWIPT) systems, evident in improvements in both energy efficiency and data rate.Consequently, TSR stands as a valuable asset for researchers and engineers engaged in the development of C-NOMA-SWIPT systems [3].In NOMA systems, the Decode-and-Forward (DF) approach serves as a cooperative relaying mechanism for SWIPT.Within NOMA setups, multiple users share the same frequency band, receiving identical signals from a singular base station (B).In SWIPT configurations, user equipment derives power from received wireless signals.Employing multiple relay nodes in the cooperative relaying process enhances communication link reliability and coverage area.In DF relaying, the relay node decodes the received signal, re-encodes it, and subsequently transmits it to the destination node.This method presents significant advantages in energy efficiency and communication performance over conventional single-hop communication systems.The fundamentals of DF relaying for SWIPT in NOMA systems and its advantages over other relaying strategies will be covered in this section [4][5][6].NOMA has recently emerged as a promising candidate for future wireless networks due to its capacity for significant improvements in spectral and energy efficiency, as well as fairness [7,8].Power domain NOMA employs techniques such as superimposed coding (SC) [9] and successive interference cancellation (SIC) [10,11] to combine transmission signals at the transmitter and decode reception signals at the receivers [12].Consequently, it enables the source to serve multiple users using the same frequency/time resources [5].Furthermore, in NOMA network systems, users are allocated different power levels to ensure fairness.This approach implies that users closer to the source receive less power compared to those farther away [13].In C-NOMA scenarios [14], a user located at a distance with weaker channel conditions may require assistance from another user to relay the signal from the source to itself.This signal-forwarding process can occur in systems with or without obstacles.The inclusion of a relay user in C-NOMA facilitates the uninterrupted propagation of information to the intended destination.However, in obstacle-free systems, referred to as C-NOMA with direct link [15][16][17], the involvement of a cooperative user ensures the integrity and reliability of the transmitted signal for the distant user, even though the signal from the source can reach this user directly.At user relay, amplify-and-decode (AF) [18], DF [19,20] schemes are utilized to harvest energy and forward information to the destination.As SWIPT [21] is exploited at the user relay in C-NOMA, it can combine with power-splitting relaying (PSR) and TSR [22] protocols for simultaneous EH and information decoding.The energy harvested at R is to maintain its lifetime due to battery limitations.In [23], a full-duplex (FD) cooperative NOMA model with one BS and two users was studied.Both direct link and no direct link between the BS and far user were considered.In [17], a C-NOMA with and without a direct link was investigated.The closed-form outage probability expression was derived for both no direct link and direct link.The OP for FD NOMA is better than that for half-duplex (HD) NOMA.Besides, throughput, ergodic rate and energy efficiency were evaluated for FD and HD NOMA.

Related Works and Motivation
C-NOMA can operate in both HD and FD modes [24,25].Due to the simultaneous reception and transmission of signals in the same time slot and frequency band, FD mode can improve the spectral efficiency of the NOMA system by up to two times compared to HD mode [14][15][16].Therefore, FD relay technology has been shown as a promising scheme for 5G-based wireless networks.In [25], an FD multi-antenna relay-assisted C-NOMA network was examined.The system comprised one BS, a cognitive near user, and a cognitive far user.An optimization problem regarding user rate was addressed, revealing that the impact of residual self-interference at the FD relay and inter-user interference in the near-user scenario could be significantly mitigated.Both analytical and simulation-based assessments were conducted to determine the outage performance in [23].The study in [23] also demonstrated that FD NOMA exhibits a lower OP than HD NOMA in the low Signal-to-Noise Ratio (SNR) region.Additionally, the system throughput for Delay-Limited Transmission (DLT) mode in FD NOMA surpassed that of HD NOMA.Furthermore, in [15], the outage performance of the proposed FD C-NOMA scheme, specifically an FD device-to-device (D2D) aided C-NOMA scheme, outperformed that of conventional NOMA and OMA.The expressions for outage probability were also derived in the same work.In [26], the authors investigated the confidentiality performance of cooperative NOMA systems integrating EH and FD relaying employing power splitting protocol over experience-independent, non-selective block Rayleigh fading.Their emphasis was on deriving closed-form expressions for EH-FD-NOMA, encompassing the lower bound on ergodic secrecy rates and the approximated secrecy OPs for the two users.In [27], the authors investigated a cooperative communication scheme that combines C-NOMA, FD relaying, and EH techniques through the utilization of a TS protocol over independent Rayleigh fading channels.They derived closed-form expressions to evaluate the performance of the proposed EH-FD-NOMA system, including OPs and ergodic rates.In [28], the authors introduced a Cooperative Spectrum-Sharing Transmission (CSST) scheme employing the DF relaying strategy.Their study incorporated realistic considerations such as FD-based loop self-interference, NOMA-based imperfect successive interference cancellation, and transceiver hardware impairments in Internet of Things (IoT) devices.Utilizing Nakagamim fading environments, a comprehensive performance analysis was conducted by deriving expressions for the OP for both primary and secondary networks within the FD-based CSST scheme under both TS and PS protocols.NOMA-SWIPT takes center stage as an indispensable requisite in the landscape of IoT networks.In a time marked by the explosive proliferation of IoT devices, where the optimization of both energy and spectral efficiency reigns supreme, NOMA-SWIPT emerges as a revolutionary remedy.By facilitating the concurrent sharing of scarce spectrum resources among multiple IoT devices, NOMA not only fine-tunes spectral utilization but also furnishes the critical capability of wireless power transfer.This dual-faceted functionality transcends mere convenience; it becomes imperative, particularly in scenarios where IoT devices are deployed in remote or inaccessible locales, challenging the provision of continuous power supply.Consequently, NOMA-SWIPT assumes the pivotal role of empowering IoT networks to unlock their full potential, cultivating sustainable, ever-connected ecosystems that catalyze innovation and drive efficiency across a spectrum of industries [29][30][31][32].In addition, Table 1 summarizes the related works.

•
An innovative C-NOMA scheme custom-tailored for SWIPT-enabled IoT systems takes center stage.This scheme orchestrates power allocation for two users through a pivotal relay node (R) in the DF of signals from the B. It is essential to emphasize that R also performs the critical task of EH from the power beacon (PB).This holistic integration of power management, signal manipulation, and EH within a unified relay node represents a groundbreaking strategy, aimed at optimizing overall system efficiency.The scheme utilizes a TS receiver architecture, enabling R to simultaneously perform both EH and Information Processing (IP) tasks.• EH Protocol Advancement: The paper introduces an advanced EH protocol with an HD-based TSR mechanism meticulously designed for SWIPT-based C-NOMA systems.This EH protocol assumes a pivotal role in the harnessing of energy from received signals, thereby bolstering the sustainability and autonomy of IoT devices interconnected within the network.

•
In an effort to provide a comprehensive understanding of system performance, the paper derives precise closed-form expressions for critical metrics, including but not limited to OP and System Throughput (ST).These analytical expressions serve as essential tools for rigorously evaluating the proposed C-NOMA and EH schemes.
Researchers and practitioners can use these insights to make informed decisions and assess the effectiveness of these systems in practical IoT deployments.

Organization and Notations
Organization: The paper is organized as follows: Section 2 introduces the system model of the relay-assisted C-NOMA along with its assumptions.Section 3 outlines the energy harvesting process at the relay node based on time-switching relaying.Section 4 delves into the analysis of performance evaluation.Section 5 focuses on determining the optimal time-switching factor to minimize the outage probability of users.Section 6 elaborates on the simulation results.Finally, Section 7 offers conclusions and summarizes the key findings of this study.The abbreviations and acronyms are presented in Table 2.
Notations: The notation E[.] stands for the operation of calculating expectation.f X (.) and F X (.) denote the probability density function (PDF) and the cumulative distribution function (CDF) of a random variable X, respectively.

System Model
The system model under investigation is illustrated in Figure 1.In this system, the signal is transmitted from source B to two users, namely U 1 and U 2 , via relay node R. R operates as a DF relay to forward the information from B to U 1 and U 2 .The relays are power-limited devices and rely on the support of a PB.We denote the parameters of the system model as follows: c r represents the channel coefficients from the source PB to R, while c 0 , c 1 , and c 2 represent the channel coefficients from the source B to R and from R to U 1 and U 2 , respectively.
Energy harvesting Information processing In this system model, we assume that there are no additional power sources embedded in the relay node to provide energy.Hence, in the initial phase, the relay node harvests energy from neighboring nodes, i.e., source B and the destinations U 1 and U 2 shown in the figure, to facilitate the forwarding of information.Specifically, this harvested energy can be stored in the batteries so that it can then be reused to forward the information from B to U 1 and U 2 .The channels have power The power gains are distributed according to exponential RVs with Ω i , i ∈ {0, 1, 2}, respectively.It is assumed that all wireless links are non-selection blocks and independent Rayleigh fading and are impacted by AWGN with average power N 0 .

The Energy Harvesting at the Relay Node
The EH scheme involving time and power is depicted in Figure 2. Within ϱT, the B transfers energy to R via a single antenna with consistent sending power (referred to as the energy-transfer phase).In this context, T and ϱ ∈ [0, 1], signify the coherence time block and the TS ratio, respectively.During the period referred to as the data transmission phase, which lasts for (1 − ϱ)T, R transmits message information signals to two users in the second hop, while B simultaneously broadcasts the aggregated signals from the users to the R node in the first hop.

R harvesting energy from PB B sends the signal to the relay node
Information is encoded and transmitted to U1 and U2 through the relay node

Phase: Energy Transfer
During this phase, the receiver uses its resources to collect radio frequency energy from the power beacon during a time slot labeled as ϱT.Here, ϱ represents the time switching ratio, and T denotes the duration of the block time period.The energy harvested through this process can be precisely described as follows: (1) In this case, η ∈ (0, 1) indicates the energy conversion efficiency of the B, while P B stands for the transmit power of the B. Let us denote |c 0 | 2 as the channel coefficient from B to R. Following the nonlinear energy-harvesting approach outlined in [38], the transmit power of R acquired over the duration of (1 − ϱ)T 2 can be expressed as: where ι is the power saturation threshold.During the initial time slot, B transmits the overlaid mixture, √ α 1 P B s 1 + √ α 2 P B s 2 , where s i and α i represent the signal and power allocation coefficient of user i, respectively.Notably, α 1 + α 2 = 1, and it is assumed that α 1 < α 2 .The baseband received signal at R can thus be expressed as follows: where n r ∼ CN (0, N 0 ) is the AWGN component at R. At first, R decodes the U 2 's message using SIC and based on the NOMA principle it decodes its own message.To decode the messages s 2 and s 1 at R, the corresponding SINR and SNR can be expressed as follows: and Here, γr N 0 depicts the SNR level at the transmitter side and the residual noise interference is represented by δ.In real-world scenarios of NOMA systems, achieving pSIC might not always be possible.Therefore, if δ = 0 , U 1 employs pSIC to decode the received signal.However, if δ ranges between 0 < δ ≤ 1, U 1 resorts to ipSIC.

Stage. 2: Communication
During the second time slot, relays transmit this information to U 1 and U 2 .Thus, the received signal at U i , i ∈ {1, 2} can be expressed as follows: Here, c i represents the channel gain between the relay and U i , and n U i ∼ CN (0, N 0 ) represents the AWGN noise at U i .Following the SIC principle, U 2 decodes symbol s 2 while treating s 1 as noise.Based on Equation ( 6), SINR at U 2 is given by: Replace Equation ( 2) with Equation ( 7) we have: where γ0,2 It is important to note that at U 1 , both s 1 and s 2 are present simultaneously.Consequently, using the SIC method, U 1 needs to decode s 1 .By considering the low-power signal s 1 as interference and utilizing SIC to eliminate it, U 1 can then extract the symbol s 1 by decoding the high-power symbol s 2 .Consequently, the received SINR for s 2 at U 1 can be formulated as follows: Here, γ0,1 Once canceled, the signal s 2 is only partially removed and trans- forms into interference.The SINRs of symbol s 1 at U 1 are given as follows: Similarly, at R, when δ = 0 and 0 < δ ≤ 1, it corresponds to the cases of pSIC and ipSIC, respectively.

Performance Assessment 4.1. Channel Characteristics
Let us examine two distinct random variables following the κ − µ distribution, each characterized by the PDF as shown in Equation ( 11).
f T e (r) = 2µ e (1+κ e ) µe +1 In the given context, where 'e' belongs to the set of real numbers, including 0, 1, and 2, κ e , which is greater than 0, represents the ratio between the combined power of the primary components and that of the scattered waves.µ e , also greater than 0, is influenced by the quantity of multipath clusters.The function I v (s) refers to the modified Bessel function of the first kind and order 'v', while re represents the root-mean-square (RMS) value of the received signal envelope, denoted as T e .As described in reference [23], we have: √ µ e (1+κ e )r e Γ(µ e ) 1 F 1 (µ e +0.5;µ e ;κ e µ e )e −κe µe Γ(µ e +0.5) , E T l e = 1 F 1 (µ e +0.5l;µ e ;κ e µ e )e −κe µe rl e Γ(µ e +0.5l) [µ e (1+κ e )] In this context, re = E[T i ], where E[s] denotes the expectation operator, and re = E[T l e ] signifies the lth moment of T e .Additionally, 1 F 1 (s; y; z) refers to the confluent hypergeometric function, and Γ(.) represents the Gamma function.
If we let γe denote the instantaneous SNR of κ − µ fading channel, we can derive the PDF of its instantaneous SNR γe from the envelope PDF provided in Equation ( 11) by employing a variable transformation r = xr 2 e Ω e , which converts γe into a fading channel, as described in Equation ( 14). where The modified Bessel function of the first kind from [39] (Eq. (9.6.16)) , ( 14 is substituted as and defining ψ e = µ e (1 + κ e ) Ω e , ( 14) can be re-written as follows: where Γ(Ξ e , ∆ e ) refer to the PDF of the gamma distribution as in Equation ( 16) The form parameter is denoted by Ξ e , while the rate parameter, represented by ∆ e , serves as an inverse scale parameter.Table 3 summarizes the notations employed in this document.

R i
The target rate is set at U i , where i belongs to the set {1, 2}.
The AWGN term was succeeded by CN (0, N 0 ).

P B
Transmit power at BS.

P r
Transmit power at R T The total time used for EH and IP.
e Subindex for defining the e th RV.Here, e ∈ {R, 0, 1, 2} The ratio of the dispersed waves to the total power of the dominant components.
µ e An extension utilizing real values that are linked to the quantity of multipath clusters. re The root-mean-square (RMS) value of the received signal envelope. re The average value of the received signal envelope.

γ e
The instantaneous signal-to-noise ratio (SNR).

Ω e
The average SNR.

Ξ e
The shaping parameter refers to a crucial aspect of a Nakagami-m RV.

∆ e
The rate parameter (also known as the inverse scale parameter) represents an optimistic measure of a Gamma RV.
The CDF of γe is given by: e µ e , 2ψ e s , s > 0, (17) where Q c (e, b) = e 1−c ∞ b s c e −(s 2 +e 2 ) 2 I c−1 (es)ds specifies the generalized Marcum Q- function for integer values of µ, and with the help of [25] (Eq.(4.63)) and [39] (Eq.(8.352.6)),F γe (s) can be simplified as follows: In this section, following a methodology similar to that outlined in [40,41], we determine the product distribution of X 0,l = γ0 γl , where l ∈ {1, 2}, for the double κ − µ fading channel.Given that the PDFs of γ0 , γ1 , and γ2 are expressed as linear combinations of the gamma distribution, it becomes straightforward to compute the PDF of the product term X 0,l = γ0 γl using [42] (Lemma 2), as depicted below: where I n,m 0,l = (µ 0 κ 0 ) m (µ l κ l ) n n!m!Γ(m+µ 0 )Γ(n+µ l ) , Φ 0,l = e µ 0 κ 0 +µ l κ l and G m,n p,q [.] is the Meijer G-function.After performing some mathematical manipulations on Equation (20), the PDF of the product distribution simplifies to Equation (21).Here, we set . ., e p ; b 1 , b 2 , . . ., b q ; z represents the regularized Hypergeometric function [43], and K j (.) stands for the jth-order modified Bessel function of the second kind.Subsequently, f X 0,l (s) is computed as: The proof of ( 21) is shown in Appendix A.
Similar to this, according to [41] (Lemma 2), the CDF of the product of two κ − µ Random Variables (RVs) can be calculated as shown in Equation ( 22).
where Φ 0,l and I n,m 0,l were defined previously.

The OP
OP is a critical metric in NOMA systems, impacting user fairness, QoS, system capacity, and overall network performance.By conducting thorough analysis and optimization, NOMA systems can ensure reliable and efficient communication, supporting various applications and services in modern wireless networks.
We define R 1 and R 2 as the target data rates for s 1 and s 2 , respectively.We then introduce ε 1 and ε 2 as defined by Subsequently, we provide expressions for the output performance measures O out U 1 and O out U 2 .

The OP in the Scenario of LEH
OP of U 1 : The OP of U 1 encompasses two scenarios: firstly, U 1 fails to detect the signal s 1 at the receiving end R. Secondly, even if U 1 successfully detects s 1 , it still might not be able to correctly identify its own transmitted message s 1 .Considering these situations, the OP of U 1 is articulated as follows: where Fγ r (.) = 1 − F γr (.) and Fγ 0,1 (.) = 1 − F γ0,1 (.) are the complementary cumula- tive distribution function, Substituting Equations ( 19) and ( 22) into Equation ( 23), we obtain the OP of U 1 as calculated in Equation (24).
When both the relay and U 2 manage to decode the signal s 2 , we denote this outcome as O out U 2 .The expression for O out U 2 can be derived from Equations ( 4) and ( 8) and can be expressed as, where Ξ 2 and Θ 2 are as defined previously.Substituting ( 19) and ( 22) into Equation ( 25), we obtain the OP of U 2 is given by

The OP in the Scenario of NLEH
In the scenario of NLEH (ι < ∞), the OP at U 1 is expressed as follows: The OP at U 1 can be rewritten as follows: where Θmax = max . It is noted that we can rewrite Equation ( 28) as follows: Combining Equations ( 17) and ( 22), the OP at U 1 with NLEH can be calculated as follows: Finally, the OP at U 2 with NLEH can be given by: We have that Equation ( 31) is eventually simplified as follows: Similarly, by solving for can be obtained as follows: (33)

System Throughput Analysis (STA)
In the delay-limited transmission mode, data are sent from the transmitter R to users at a fixed data rate.During this mode, the transmitter B maintains a steady transmission rate denoted as R, which can be impacted by the occurrence probability of outages owing to fluctuating wireless channel conditions.The relative throughput of the NOMA system is expressed by Equation [22] (Eq.( 14)) (refer to (34)).
Remark 1.Because the cumulative throughput depends on the OPs of U i and the desired rate R i , it is upper-bounded by the ceiling 2 ∑ i=1 R i .Moreover, power allocation at both the transmitter B and the relay R also contributes to a decrease in throughput due to their detrimental effects on outage performance.

Pseudo Code
To offer a better understanding of the power beacon-assisted scheme, a flowchart outlining the methodology is also provided in Figure 3. From expression (12), the study presents an algorithm for the Monte Carlo simulation with two users (U i , i ∈ {1, 2}), detailed in Algorithm 1.

Optimization of the TS Factor to Minimize the OP for Users
This section focuses on optimizing the TS factor ϱ, which is a crucial parameter for effectively balancing the EH and data transmission phases, aiming to minimize the OP of users.Below is the formula for calculating the OP of users: The OP function of U i versus ϱ exhibits an optimal ϱ * on the interval τ, 1 − 2R 2 log 2 (1/α 2 ) − τ τ ≃ 0 with τ ≃ 0, similar to references [44,45].The credibility of these convex function characteristics was validated in [46].However, no supporting evidence is presented.It is reasonable to reformulate the problem outlined in Equation ( 35) as follows: Given the computational challenges associated with employing gradient methods to solve the formula expression P2, we utilize the Golden Section Search technique, as elaborated in Table 2. Here, ν represents the search accuracy, while O[log 2 (1/ν)] indicates the com- plexity of the method, as stated in the work by Nguyen in [21].This technique is widely acknowledged for its combination of high precision and low computational complexity, as highlighted by Nguyen.Furthermore, Nguyen demonstrated that for unimodal functions, this method consistently achieves convergence to the global optimum.

Energy Efficiency
To enhance throughput in delay-constrained transmission scenarios, further research is crucial to explore the energy efficiency of C-NOMA systems.The necessity for additional studies is highlighted by the reference in [17] (Eq.( 52)), which indicates that the formula for the energy efficiency coefficient can yield valuable insights into the performance metrics of these systems.Therefore, examining the details of this coefficient can provide a deeper understanding of how C-NOMA setups can optimize energy consumption while maintaining high throughput levels.The energy efficiency coefficient can be expressed as follows: Total data rate Total energy consumption .
Consequently, the energy efficiency metrics of C-NOMA systems can be redefined as follows: In this context, the variable T represents the duration of the transmission procedure, covering all aspects of the communication process.

Simulation Results
In this section, we evaluate the OP and throughput of the proposed system through analytical analysis, supplemented by simulation results obtained using the Monte Carlo method to validate our analytical findings.The parameters of the κ − µ random variables are assumed to be identical, namely Unless otherwise stated, all simulation configurations are provided in Table 4.Moreover, we employed the κ − µ fading model, and a variety of examples can be found in Table 5.For scenarios involving NLEH, we set ι = 0 dBm and N 0 = −90 dBm as specified in [47].Table 5.The κ − µ distribution yields common fading distributions [48].
The two κ − µ random variables associated with the EH-NOMA system discussed in Section 3 are clearly depicted in Figure 4.This visualization illustrates that the simulation outcomes closely align with the numerical analysis we conducted, which was derived from simulation data.Figure 5 shows the representation of OP according to ρ B in a cooperative NOMA system with different κ and µ values.Specifically, Figure 5a presents the OP of the Rice distribution (µ = 1) versus two values κ = 1 and κ = 4. From the figure, it can be observed that the larger the κ, the better the OP of the system.A similar explanation is provided for Figure 5b for the Nakagami-m distribution (κ → 0) versus two values µ = 1 and µ = 3.From the figure, we can also see that the larger the value of µ, the lower the OP.These results can be explained based on Equations ( 23)- (26).Generally, the alignment between simulation results and analytical results in NOMA systems can be attributed to several key factors: accurate user pairing, optimized power allocation, and precise interference management through SIC.These elements collectively enhance the predictability and stability of the channel environment, ensuring that theoretical models closely match simulations.The the OP values for U 1 and U 2 as a function of ρ B are presented in Figure 7a and Figure 7b respectively.The figure reveals that the Nakagami-m distribution has the lowest OP value, while the One-sided Gaussian distribution has the highest OP value.This suggests that the transmission line with the Nakagami-m distribution is the most stable.This observation can also be explained using Equations ( 23)- (26).Upon applying NLEH, an enhancement in outage performance for the downlink scenario is apparent in Figure 8.Here, we set δ = 0.05.It is notable that as the SNR ρ B during transmission increases, a considerable gap emerges between the outage probabilities attained by the linear and non-linear EH models.This disparity arises due to the extended duration required for EH to accumulate the minimum power influenced by ι, consequently diminishing the OP. Figure 9 illustrates OP versus the power allocation factor a 2 , where a larger value of ρ B indicates better optimization.Based on the figure, we observe that the OP occurs at U 1 when the ipSIC ceases to function entirely within the α 2 range from 0.86 to 1. Similarly, the OP occurs at U 1 in the case of pSIC complete outage within the α 2 range from 0.97 to 1.These observations can be explained similarly to Figures 3 and 4, as well as based on Equations ( 23)- (26).Figure 10 plots the OP versus ϱ for various target rate values, with parameters δ = 0.02, α 1 = 0.1, α 2 = 0.9, κ = 0, and µ = 3.Specifically, as the target rate values change from R 1 = R 2 = 0.2 to R 1 = R 2 = 0.5, the OP also increases.We can explain these observations based on Equations ( 23)-( 26). Figure 11 describes the OP at U 1 and U 2 versus target rate values (R 1 = R 2 ).Specifically, the figure illustrates that OP increases with increasing target rate values.
Figure 12 illustrates the ST at the destination node versus ϱ, where R 1 = 0.5, R 2 = 0.25, α 1 = 0.3, α 2 = 0.7, κ = 0, and µ = 2.The exact theoretical curves of the system throughput are plotted according to Equation (21).From the graph, it is evident that the system throughput increases as the value of ρ B increases.When ρ B changes from 10 [dB] to 15 [dB], the maximum system throughput nearly doubles.Additionally, the system throughput is dependent on δ, with δ = 0 corresponding to the dashed line and δ = 0.05 corresponding to the solid line.As δ increases at the same value of ρ B , the system throughput decreases rapidly, and vice versa.This implies that the quality of the system is dependent on the quality of the SIC.In Figure 13, we adjust the simulation settings to match those [36], and the outage probability of an energy harvesting-based NOMA system falls dramatically in the high SNR zone.Such patterns in outage performance are observed for both time and power splitting systems, which are consistent with our findings.However, target rates remain a restriction of our system's outage performance, as indicated in [36].It is concluded that our relay-based NOMA system outperforms that of [36] in terms of outage performance in the case of time-switching energy harvesting at SNRs greater than 20 dB.We show the outage probability at the users versus ρ for the LEH and NLEH cases in Figure 14.It can be noted that when ρ increases, the NOMA system's outage performance improves.This is because a greater ρ corresponds to a higher diversity order for each user, resulting in a reduced outage probability.Furthermore, the study in Figure 8 shows that the users in the NLEH scenario outperform those in the LEH case. Figure 15 illustrates the relationship between system energy efficiency and SNR when operating in delay-limited transmission mode for user relaying within NOMA systems.The plotted graph includes dashed curves representing the performance of two users relaying NOMA with the LEH scenario, derived from Equations ( 24) and (26) based on throughput in delay-limited transmission mode.Additionally, solid curves on the graph depict the results for two users relaying NOMA with the NLEH setup, calculated using Equations (A1) and (A4) considering throughput in transmission mode.Upon observation, it becomes apparent that the energy efficiency of user relaying in NLEH NOMA during delay-limited transmission mode exhibits a higher throughput compared to that of LEH NOMA under similar transmission conditions.Through the comparison of these scenarios, it is evident that the NLEH configuration offers improved energy efficiency for user relaying within NOMA systems, particularly in delay-limited transmission mode.The findings from the plotted curves emphasize the advantages of utilizing NLEH over LEH in terms of energy efficiency and throughput performance when implementing user relaying in NOMA systems.This insight underscores the significance of considering different energy harvesting mechanisms in NOMA setups to enhance system performance and efficiency.

Conclusions
This paper investigates the DF relay-based protocol for EH and IP in wireless C-NOMA relaying networks.We derived closed-form expressions for the OP, ST, and the optimal TS factor to minimize OP for two users.Additionally, expressions for achievable OP, ST, and the optimal TS factor for users of the TSR protocol for SWIPT were obtained.Analytical and simulation results consistently demonstrate the superiority of the Nakagami-m fading channel over other fading channels in terms of energy efficiency performance.For instance, simulations with parameters such as κ = 0, µ = 3, and SNR ranging from 10 dB to 20 dB showed up to 20% lower outage probabilities compared to Rice and One-sided Gaussian channels.Despite these promising findings, our study is not without limitations.One limitation is the assumption of identical parameters (κ = κ R = κ 0 = κ 1 = κ 2 and µ = µ R = µ 0 = µ 1 = µ 2 ) across different scenarios, which may not fully capture real-world variations in channel conditions.Additionally, while our simulations provide valuable insights, they are based on specific configurations and may not generalize to all practical deployment scenarios.Moving forward, we propose the development of a C-NOMA network system utilizing multiple-input-multiple-output (MIMO) NOMA to further enhance system performance and address these limitations.This avenue presents promising opportunities for future research and practical implementation, aiming to optimize wireless communication networks through advanced NOMA techniques.

Figure 3 .
Figure3.The flowchart for calculating the primary system metrics for the power beacon-assisted scheme.

Figure 4 .
Figure 4. Visualization of the derived closed-form expression for the PDF and CDF, alongside the graphical representation of the predicted actual distributions.

4 Figure 5 .
Figure 5. OP ρ B with LEH for the Cooperative NOMA system with varying κ and µ values.

Figure 6
Figure 6 findings are identical to Figure 5 in terms of general channel characteristics;we can see that for the situation, NLEH will be better than LEH, as shown in Figure5.

4 Figure 6 .
Figure 6.OP versus ρ B with NLEH for the Cooperative NOMA system with different κ and µ values.

Figure 7 .
Figure 7. OP of U 1 and U 2 versus ρ B for different fading channels, with δ = 0.1.

Figure 14 .
Figure 14.The outage probabilities at the other users versus ρ B .

Table 1 .
Comparison between the novelty of our work and previous papers (X: Not Discussed; ✓: Discussed.

Table 4 .
Definition of system parameters.