Next Article in Journal
Assessing Visual Exploratory Activity of Athletes in Virtual Reality Using Head Motion Characteristics
Next Article in Special Issue
Spectrum Occupancy Model Based on Empirical Data for FM Radio Broadcasting in Suburban Environments
Previous Article in Journal
Multi-Satellite Relative Navigation Scheme for Microsatellites Using Inter-Satellite Radio Frequency Measurements
Previous Article in Special Issue
Real-Time Implementation of Multiband Spectrum Sensing Using SDR Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Analysis of Underlay Cognitive Radio System with Self-Sustainable Relay and Statistical CSI

School of Electrical Engineering, University of Belgrade, 11000 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(11), 3727; https://doi.org/10.3390/s21113727
Submission received: 7 March 2021 / Revised: 13 May 2021 / Accepted: 24 May 2021 / Published: 27 May 2021
(This article belongs to the Special Issue Cognitive Radio Applications and Spectrum Management)

Abstract

:
The relentlessly increasing number of small-sized devices with limited powering and computational capabilities requires the adoption of new approaches to spectrum access. In this paper, we analyze an underlay cooperative cognitive wireless system based on available statistical channel state information (CSI) that is applicable to the cognitive system with limited computational resources due to its low complexity. We considered the scenario where the primary and the cognitive network coexist in the same spectrum band, under the constraints of interference threshold and maximal tolerable outage permitted by the primary user. The communication in the secondary decode-and-forward (DF) relaying system is established via a self-sustainable relay, which harvests energy from both cognitive and primary transmitters. The closed-form expressions for the outage probability of the cognitive network are derived, which are valid for both time-switching relaying (TSR) and power-splitting relaying (PSR) protocols. We analyze the influence of both cognitive and primary systems as well as the impact of channel parameters on the cognitive system outage performance. The derived analytical results are corroborated by an independent simulation method.

1. Introduction

Recent advances in telecommunication concepts have resulted in ubiquitous connectivity, as well as a variety of possible services of improved quality. However, due to increasing demands for dedicated spectrum resources, it is necessary to have extremely efficient and innovative approaches to the usage of limited spectrum resources. The traditional static spectrum allocation has led to its inefficient use (both in time and at various locations) and consequently to its scarcity for the deployment of new services. It is envisioned that forthcoming networks will satisfy a wide range of requirements regarding different types of services, so spectrum access should be adjusted with the aim to maximize overall spectrum utilization. This is especially important for the Internet of Things (IoT) systems that have an exceptionally important role in transforming the way communications are used in everyday life. IoT systems integrated with 5G offer new services and smart communications such as the smart grid, smart health, smart traffic, etc. Nevertheless, as IoT systems imply a huge number of spatially distributed connected devices and sensors with modest powering and computational capabilities, access to spectrum resources needs to be provided in a manner that is adapted to their characteristics and limitations [1,2].
It is well known that cognitive radio can significantly improve spectrum efficiency by adopting an underlay, overlay, or interweave transmission mode [3]. While the interweave approach involves energy consuming spectrum sensing, the underlay approach enables a coexistence of licensed and unlicensed networks. In the underlay concept, the cognitive, secondary user (SU) communicates concurrently in the same spectrum band dedicated to the existing primary user (PU) without causing harmful interference to the PU. SUs apply the real-time transmit power adaptation based on the available channel state information (CSI), such that the interference level caused at the licensed user’s receiver is within predetermined tolerable limits [4,5]. However, the fulfillment of this condition is possible only in the case when perfect CSI is available to the SU, which is not practically achievable. In the underlay cognitive system, the CSI can be provided to the SU through cooperation with the PU (if the channel is reciprocal) or by mediation of the band-manager. However, the acquired CSI is practically always imperfect, even in the case when there is cooperation with the PU [6,7,8,9,10]. Due to the time-varying nature of the channel and the finite backward channel latency, the acquired CSI is outdated. Even in the case when the propagation environment is slow varying, the backward channel quantization results in an imperfect CSI. In the absence of cooperation, only statistical CSI (i.e., long-term CSI) can be available [10,11,12,13]. Notably, in any case when perfect CSI is not available, the strict interference constraint cannot be satisfied, and the interference threshold will be exceeded during a certain percentage of time, resulting in degraded primary and secondary system performance. The adaptation of transmission power based on real-time updating of the CSI knowledge is not desirable in the case when there are a large number of cognitive nodes with limited computational capabilities. The underlay cognitive network based on statistical CSI enables a less complex solution for coexistence with the primary spectrum user’s network.
Apart from the constant demands for improved spectrum efficiency, the relentlessly increasing number of devices in a network also necessitates energy-efficient solutions. In this paper, we focus on the application of energy-harvesting techniques [14,15,16,17,18,19], which can improve the performance of an underlay cognitive radio system. In the rest of this section, we will provide an overview of related work, summary of the contributions, and explain the organization of the paper.

1.1. Related Work

The performances of underlay cognitive radio with available imperfect CSI are analyzed in [6,7,8]; these works introduce the concept of interference outage, as the strict interference threshold constraint at the PU receiver cannot be satisfied. The underlay cognitive network based on statistical CSI enables lower complexity solution for coexistence with the primary spectrum user’s network [9]. In the paper [10], the outage probability and ergodic capacity were obtained for the case when transmit power adaptation is done based on the CSI that is outdated, as well as for the case when only statistical properties of the channel from the SU transmitter to the PU receiver are known. The impact of CSI knowledge on the underlay system performances is analyzed in the context of relaying [11,12] and application of space-time codes under multiple interferers’ environment [13].
The energy-harvesting technology represents a promising way to enable a practically feasible and energy-efficient powering of the nodes [14]. Although energy can be harvested from various natural sources (solar, thermal, vibration energy, etc.) [15], harvesting radio frequency (RF) energy is very appealing, as it provides necessary reliability in the energy supply [16]. In the field of harvesting energy from RF sources, technology has achieved many advances both in industry [17] and academia [18] in recent years, and much more is expected in the future [19].
The performance analysis of the wireless network where a transmitter harvests energy from both interference and ambient RF sources is presented in [20]. The combined use of a cognitive network [21,22] with the cooperative energy-harvesting approach [23,24,25] can further enable the coexistence of two networks in the same spectrum band with the controlled mutual interference that enables energy sustainability of network relay nodes [26].
The concept of cooperative decode-and-forward (DF) relay that harvests energy from both source and interference signals in the Rayleigh fading environment is proposed in [24] and further analyzed in [27] for Nakagami-m environment and nonlinear energy harvester. Energy harvesting is analyzed in [28] for cooperative relaying cognitive network with perfect CSI, and in [29,30,31] for the case of imperfect CSI. However, to the best of authors’ knowledge, there is no available investigation concerning the impact of the statistical CSI on the performance of the cognitive cooperative relaying network that employs energy harvesting.

1.2. Summary and Organization

In this paper, we provide a comprehensive outage analysis of the cooperative cognitive system with a DF relaying scheme. The relay is energy-constrained and relies on the energy harvested from both primary and secondary transmitters, which coexist in the same spectrum band by using statistical CSI knowledge. Two important cases are analyzed: harvesting at the relay node is performed based on time-switching relaying (TSR) or power-splitting relaying (PSR) protocols, as proposed in [32]. It is assumed that the communication between nodes in the cooperative cognitive system is dominantly limited by the interference signals generated by the PU’s system.
We derive closed-form expressions for the outage probability of the cognitive energy-harvesting DF relaying system constrained by the interference threshold at the primary receiver and tolerable outage probability of the primary system. The derived results are valid for the Rayleigh fading environment. It is important to emphasize that the Rayleigh fading model encompasses the path-loss effect, while it does not provide shadowing effects. In the considered system model, we assume that there are no obstacles between nodes that would introduce a shadowing effect. The numerical results are obtained on the basis of derived expressions and corroborated by using an independent simulation method for various system parameters and propagation conditions.
The rest of the paper is organized as follows. System and channel models are described in Section 2, while outage probability analysis is provided in Section 3. Numerical results, based on both analytical and simulation approach, are presented in Section 4. The final conclusions are outlined in Section 5.

2. System and Channel Models

This paper considers a cognitive cooperative radio system in which secondary users coexist with the primary users of the spectrum in the same frequency band under the predefined conditions regarding permitted introduced interference. The model of the analyzed cognitive relay system is shown in Figure 1. Our further analysis is based on the following set of assumptions:
A1. The secondary system consists of the secondary source (SS), the energy constrained secondary relay node (SR), and the destination of the cognitive secondary system (SD). There is no direct link between the SS and the SD due to the occurrence of deep fading, so the SS and the SD can communicate only with the help of the SR [11,24,32,33,34].
A2. The secondary system shares the spectrum with the primary network based on the underlay paradigm, which means that secondary users can perform concurrent transmission as long as the interference generated at the primary receiver does not exceed a permissible threshold. It is assumed that only statistical CSI is available to the secondary network, i.e., that the mean channel power gain values are known (long-term CSI). The strict interference power constraint cannot be satisfied in this scenario, and the outage interference power constraint is applied as in [10,11,12,13].
A3. The primary network consists of N transmitters (PUTX) and M receivers (PURX). PU transmitters are closely positioned at one location, while all PURX are closely positioned at another location [12,33,34]. This corresponds to the case of cooperative distributed network where the corresponding PUTX form a transmit virtual antenna array [35], or the case where PU’s transmitters and PU’s receivers are grouped in two wireless sensor network node clusters [36].
A4. The fading envelopes in all channels follow the Rayleigh distribution. This model encompasses the path loss effect as the mean channel power gain equals d - β , where d denotes the distance between the transmitter and the receiver, while β is the path loss exponent.
A5. In the interference-limited case, the interference power caused by PU transmitters at the SR and the SD is dominant relative to the noise power [11,12,33,34].
A6. The SR does not have its own power supply, and it harvests energy from the secondary source and the N transmitters of the primary network. According to the applied protocol, all the energy harvested during a scheduled frame is used for information transmission [24,32,34].
As it is shown in Figure 1, the channel gain coefficients from the SS to the SR and from the SR to the SD are denoted by hSR and hRD, respectively. Furthermore, the channel gain coefficients from the SS and the SR to the PURX-j are denoted by g S j and g R j for j = 1, 2, …, M, respectively. Finally, the channel gain coefficients from the PUTX-i to the SS, the SR and the SD are denoted by f S i , f R i and f D i for i = 1, 2, …, N, respectively.
The secondary transmitter is allowed to transmit in the same frequency band under the constraint that the peak interference threshold, denoted by Qp, is exceeded with a maximum tolerable probability P o u t , P U . Therefore, the maximal transmitted power at the SS, PS, is determined by the permitted interference outage probability, i.e.,
Pr { P S g S j < Q p } = 1 P o u t , P U ,
where Pr{⋅} denotes probability.
The transmit power of the SS is adjusted by using the mean value of the power gain in the links from the SS to the PU’s receivers
P S , max = k S Q p Λ S ,
where kS is the coefficient of adjustment of the SS transmit power and Λ S is the mean value of all power gain coefficients | g S j | 2 (j = 1, 2, …, M) in the links from the SS to each of the PU’s receivers.
The transmit power of the SR is adjusted by using the mean value of the power gain in the links from the SR to the PU’s receivers
P R , m a x = k R Q p Λ R ,
where kR is the coefficient of the adjustment of the SR transmit power and Λ R is the mean value of all power gain coefficients | g R j | 2 (j = 1, 2, …, M) in the links from the SR to each of the PU’s receivers.

2.1. Time-Switching Relaying Protocol for Energy Harvesting

Figure 2 depicts the transmission block structure in the TSR scheme for energy harvesting and information processing at the SR. In the figure, T represents the time frame interval and α ∈ {0, 1} denotes the fraction of the block time in which the SR harvests energy from the RF signals from the SS and N PU’s transmitters. The remaining block time (1 − α)T is used for the information transmission in such a way that half of that time, (1 − α)T/2, is used for the information transmission from the SS to the SR, while the remaining half, (1 − α)T/2, is used for the information transmission from the SR to the SD.
The received baseband signal at the SR is given by
y R = P S , max h S R s + i = 1 N P P U T X , i f R i s P U T X , i + n R ,
where s is the information signal sent from the SS, PPU-TX,i is the power of the signal from the i -th PU transmitter, the interference signal from PUTX-i (with the unit power) is denoted by sPU-TX,i and the nR is additive white Gaussian noise component at the relay node.
We consider the case of interference-limited environment, where the interference power caused by N PU’s transmitters at the relay node SR is dominant relative to the noise power at the relay and therefore can be neglected in derivations. Then, the received SIR at the SR can be written in the form
γ R = P S , max | h S R | 2 i = 1 N P P U T X , i | f R i | 2 = γ S R I R ,
where γ S R = P S , max | h S R | 2 and I R = i = 1 N P P U T X , i | f R i | 2 .
The harvested energy at the relay node of the secondary system is
E H = η ( P S , max | h S R | 2 + i = 1 N P P U T X , i | f R i | 2 ) α T = η ( γ S R + I R ) α T ,
where η (0 < η < 1) is the energy conversion efficiency coefficient that depends on the energy harvester.
It is assumed that the total energy harvested at the SR during the dedicated time can be stored in a storage device (e.g., a supercapacitor or a short-term/high-efficiency battery) and be used afterwards for the power transmission in the same time frame. Therefore, for the harvested energy EH, the maximal available transmit power of the SR can be expressed as
P R , E H = E H ( 1 α ) T / 2 .
The transmit power of the energy constrained relay node is determined by both harvested energy and maximal power permitted due to interference outage constraint. Using Equations (3) and (7), the transmit power of the SR can be obtained as
P R = min ( P R , E H , k R Q p Λ R ) .
The transmit power of the secondary relay can be further expressed in the following form
P R = { 2 η α ( 1 α ) ( γ S R + I R ) , 2 η α ( 1 α ) ( γ S R + I R ) k R Q p Λ R , k R Q p Λ R , 2 η α ( 1 α ) ( γ S R + I R ) > k R Q p Λ R .
Then, the received information signal at the secondary destination is equal
y D = P R h R D s R + i = 1 N P P U T X , i f D i s P U T X , i + n D ,
where sR is the decoded information signal sent from the secondary relay, the sPU-TX,i is the interference signal from the i-th primary user transmitter, and nD is the additive white Gaussian noise component at the secondary destination. The interference power caused by N PU’s transmitters at the SD is dominant relative to the noise power; i.e., the communication is interference-limited. The received SIR at the SD can be written in the form
γ D = P R | h R D | 2 i = 1 N P P U T X , i | f D i | 2 = P R γ R D I D ,
where γ R D = | h R D | 2 and I D = i = 1 N P P U T X , i | f D i | 2 .
By using Equations (9) and (11), the expression for the received SIR at the SD becomes
γ D = { 2 η α ( 1 α ) γ R D I D ( γ S R + I R ) , 2 η α ( 1 α ) ( γ S R + I R ) k R Q p Λ R , k R Q p Λ R γ R D I D , 2 η α ( 1 α ) ( γ S R + I R ) > k R Q p Λ R .

2.2. Power-Splitting Relaying Protocol for Energy Harvesting

The transmission block structure for the information and power transmission according to the PSR protocol is presented in Figure 3. The total time T is split into two equal parts. In the first part of duration equal to T/2, the source transmits information to the relay, and in the second part of duration equal to T/2, the relay further transmits information to the destination. In addition to the information transmission, the first part is used for the energy harvesting. The fraction of the received signal power equal ρP (0 < ρ < 1) is used for that purpose, while the remaining received power equal (1 − ρ)P is used for the signal transmission from the source to the relay.
The received baseband signal at the SR is given by
y R = ( 1 ρ ) P S , max h S R s + i = 1 N ( 1 ρ ) P P U T X , i f R i s P U T X , i + n R ,
where s is the information signal sent from the SS, the sPU-TX,i is the interference signal from the i-th primary user transmitter and with the power equal PPU-TX,i, while nR is the additive white Gaussian noise component at the SR. As in the paper [24], the case when power splitting reduces the signal power and not the noise power is considered, which provides a lower-bound system performances.
The interference power caused by N primary user transmitters at the SR is dominant relative to the noise component power and can be neglected in the further calculations. The received SIR at the secondary relay node can be written in the following form
γ R = ( 1 ρ ) P S , max | h S R | 2 i = 1 N ( 1 ρ ) P P U T X , i | f R i | 2 = P S , max | h S R | 2 i = 1 N P P U T X , i | f R i | 2 = γ S R I R ,
where γ S R = P S , max | h S R | 2 and I R = i = 1 N P P U T X , i | f R i | 2 .
The energy harvested at the SR is equal to
E H = η ρ ( P S , max | h S R | 2 + i = 1 N P P U T X , i | f R i | 2 ) T / 2 = η ρ ( γ S R + I R ) T / 2 .
Similarly as in the case of TSR scheme, the maximal available transmit power of the SR for the harvested energy EH can be expressed as
P R , E H = E H T / 2 .
We assume that the SR transmits with the maximal allowed power, which is limited by both the quantity of the harvested energy and the interference constraint imposed by the primary network. By using Equations (3) and (16), the transmit power of SR can be expressed as
P R = min ( P R , E H , k R Q p Λ R ) ,
and the transmitted power of the secondary relay can be further presented in the following form
P R = { η ρ ( γ S R + I R ) , η ρ ( γ S R + I R ) k R Q p Λ R , k R Q p Λ R , η ρ ( γ S R + I R ) > k R Q p Λ R .
The received information signal at the SD is
y D = P R h R D s R + i = 1 N P P U T X , i f D i s P U T X , i + n D ,
where sR is the decoded information signal sent from the sR, sPU-TX,i is the interference signal from the i-th primary user transmitter, and nD is the AWGN component at the SD.
The interference power caused by N primary user transmitters at the SD is dominant relative to the noise power and the noise impact can be neglected in the further analysis. The received SIR at the SD can be written in the form
γ D = P R | h R D | 2 i = 1 N P P U T X , i | f D i | 2 = P R γ R D I D ,
where γ R D = | h R D | 2 and I D = i = 1 N P P U T X , i | f D i | 2 .
By using Equations (18) and (20), the expression for the received SIR at the SD becomes
γ D = { η ρ γ R D I D ( γ S R + I R ) , η ρ ( γ S R + I R ) k R Q p Λ R , k R Q p Λ R γ R D I D , η ρ ( γ S R + I R ) > k R Q p Λ R .

3. Outage Probability Analysis

The outage probability of the secondary network is defined as the probability that one of the links is below predefined threshold γ t h , and it is given by the following expression
P O U T ( γ t h ) = Pr { γ R γ t h } + Pr { γ D γ t h , γ R > γ t h }
It can be noticed that the SIR at the SD for both TSR and PSR protocols, given by Equations (12) and (21) respectively, can be expressed in the unique form and re-written as
γ D = { c γ R D I D ( γ S R + I R ) , c ( γ S R + I R ) q , q γ R D I D , c ( γ S R + I R ) > q ,
where q = k R Q p Λ R , c = c T S = 2 η α ( 1 α ) for the TSR protocol and c = c P S = η ρ for the PSR protocol.
By using Equations (22) and (23), the following expression is obtained
P O U T ( γ t h ) = Pr { γ S R I R γ t h } + Pr { c γ R D I D ( γ S R + I R ) γ t h , c ( γ S R + I R ) q , γ S R I R > γ t h } + Pr { q γ R D I D γ t h , c ( γ S R + I R ) > q , γ S R I R > γ t h } .
By substituting u = γ S R + I R , ν = γ S R I R and w = γ R D I D , the outage probability can be written as
P O U T ( γ t h ) = Pr { v γ t h } + Pr { u min ( γ t h c w , q c ) , v > γ t h } + Pr { w γ t h q , u > q c , v > γ t h } ,
and it can be further written in the form more convenient for the analysis that follows
P O U T ( γ t h ) = Pr { v γ t h } + Pr { u γ t h c w , w > γ t h q , v > γ t h } + Pr { w γ t h q , v > γ t h } = I 1 + I 2 + I 3 .
The propagation channels among primary and secondary transmitters and the receivers are subject to Rayleigh distributed fading. Therefore, the random variables (RVs) γ S R and γ R D are Gamma distributed, and the corresponding probability density functions (PDFs) are given by
p γ S R ( x ) = 1 Ω S R exp ( x Ω S R ) ,
p γ R D ( x ) = 1 Ω R D exp ( x Ω R D ) ,
with the mathematical expectations (mean values) equal to Ω S R = E [ P S , max | h S R | 2 ] and Ω R D = E [ | h R D | 2 ] , respectively. Furthermore, the sum of N independent Gamma distributed variables follows χ2-distribution with 2N degrees of freedom [37] and
p I R ( x ) = x N 1 Γ ( N ) Ω I R N exp ( x Ω I R ) ,   Ω I R = E [ P P U T X , i | f R i | 2 ] ,
p I D ( x ) = x N 1 Γ ( N ) Ω I D N exp ( x Ω I D ) ,   Ω I D = E [ P P U T X , i | f D i | 2 ] ,
where Γ(⋅) denotes Gamma function [38] (8.310.1).
Further calculations with appropriate mathematical operations presented in Appendix A lead to the final closed-form expression for the outage probability.
The probability defined by integral I1 can be calculated in the following way
I 1 = Pr { ν γ t h } = 1 1 ( 1 + Ω I R Ω S R γ t h ) N
The probability defined by the expression I2 can be further expressed as
I 2 = Pr { u γ t h c w , w > γ t h q , v > γ t h } = γ t h q ( 0 γ t h c w γ t h p u , v ( u , v ) d v d u ) p w ( w ) d w
By mathematical derivations described in Appendix B, Equation (32) is written in the following closed form
I 2 = Ω I R N ( 1 Ω I R 1 Ω S R ) N ( I 21 I 22 k = 0 N 1 ( 1 + γ t h ) Ω I R k ! ( γ t h Ω I R + Ω S R ) ( Ω S R Ω I R γ t h Ω I R + Ω S R ) k ( Γ ( 1 + k ) I 21 I 23 ) ) ,
where integrals I21 and I22 are given with
I 21 = γ t h q p w ( w ) d w = ( Ω I D Ω R D γ t h q + 1 ) N ,
and
I 22 = γ t h q e γ t h Ω S R c w p w ( w ) d w = N e γ t h c Ω S R Ω I D Ω R D i = 0 N 1 ( N 1 i ) ( 1 ) N 1 i ( Ω I D Ω R D ) N i × ( ( Ω I D Ω R D ) i N E i ( 1 i + N , γ t h c Ω S R Ω I D Ω R D ) ( Ω I D Ω R D + q γ t h ) i N E i ( 1 i + N , ( Ω I D Ω R D + q γ t h ) γ t h c Ω S R ) ) ,
while integral I23 is defined and solved as
I 23 = γ t h q Γ ( 1 + k , ξ w ) p w ( w ) d w = k ! N Ω I D Ω R D j = 0 k 1 j ! ( ξ ) j e ξ Ω I D Ω R D i = 0 N + j 1 ( N + j 1 i ) ( Ω I D Ω R D ) N + j 1 i × ( ( Ω I D Ω R D ) i N E i ( 1 i + N , Ω I D Ω R D ξ ) ( Ω I D Ω R D + q γ t h ) i N E i ( 1 i + N , ( Ω I D Ω R D + q γ t h ) ξ ) ) ,
where ξ = γ t h Ω I R + Ω S R Ω S R Ω I R ( 1 + γ t h ) γ t h c .
The detailed derivations of integrals I2i, i = 1, 2, 3 are provided in Appendix B.
The probability defined by integral I3 is equal to
I 3 = Pr { w γ t h q , v > γ t h } = Pr { w γ t h q } Pr { v > γ th } = ( 1 ( 1 + Ω I D Ω R D γ t h q ) N ) ( 1 + Ω I R Ω S R γ t h ) N .

4. Simulation Method and Numerical Results

In this section, we provide the numerical results for the outage probability and the throughput of the proposed cognitive EH relay system with available statistical CSI. The results are presented for various system parameters and propagation scenarios. The numerical results are obtained by applying the derived analytical closed-form expressions and verified by an independent simulation method. Simulation results are obtained by using a Monte Carlo method, based on waveform sequences with L = 107 samples. As it will be shown in the following part, the results obtained by both approaches are in excellent agreement.

4.1. Simulation Environment

The simulation environment is developed to verify the accuracy of the derived analytical expressions. It is based on Monte Carlo simulations, which is typically used to estimate performance of the wireless communication systems.
The channel gains are represented as mutually independent random processes with Rayleigh distribution. The mean channel gains from the group of primary users in the same cluster to the other node have the same value, i.e., the corresponding random variables are independent but identically distributed (i.i.d.). Waveform sequences are generated for all channel coefficients by using an improved Jakes fading simulator [39], with M0 = 20 oscillators and maximum Doppler frequency fD,max = 20 Hz (this value does not affect the simulation results, as temporal fading characteristics are not relevant in this system scenario; if the outdated CSI was analyzed instead of statistical SCI, it would be significant). Therefore, sequences hSR(k), hRD(k), gSj(k), gRj(k), fRi(k), and fDi(k) are generated for i = 1, 2, …, N, j = 1, 2, …, M, and k = 1, 2, …, L where L is the number of samples in each waveform and every sample corresponds to one time slot.
The system simulation is done based on the available statistical CSI, meaning that the transmitted power of the secondary user is adapted based on mean channel power gain values. Therefore, the maximum output power at the SS is determined by using Equation (2), where the coefficient of adjustment is calculated from (1) for given Pout,PU and Qp. Depending on the used protocol, further analysis differs, which is shown by using different equations for calculations. In both protocols, as described in Section 2, the SR harvests energy from the SS and from N transmitters of the primary network. As parameters PPU-TX-I, PS,max, α, η, and T are fixed and known and the waveforms hSR(k) and fR(k) are generated, the harvested energy in the k-th slot at the relay node of the secondary system is calculated using Equations (6) and (15) for TSR and PSR protocol, respectively. Depending on the harvested energy and maximal allowable power, the transmit power of the SR in the k-th slot is adjusted by using the mean value of the power gain in the channels from the SR to the PU’s receivers. It is obtained using Equations (8) and (17) for TSR and PSR protocol, respectively.
We consider the case of interference-limited environment, where the interference power caused by N PU’s transmitters at the SR or SD is dominant relative to the corresponding noise power. The signal-to-interference ratio on the SR in the k-th slot is obtained using Equations (5) and (14) for TSR and PSR protocol, respectively. In addition, the signal-to-interference ratio at the SD in k-th slot is obtained using Equations (11) and (20) for TSR and PSR protocol, respectively. The outage probability of the secondary system is estimated by comparing values of the signal-to-interference ratio at the SR and the SD nodes with the outage threshold in each particular slot.

4.2. Numerical Results

In this section, we present numerical results for outage probability and achievable throughput for both analyzed protocols. Numerical results that correspond to the analytical approach are obtained directly from Equations (26) and (31)–(37). We also wrote the scripts for estimating outage probability for a proposed cognitive system by using independent Monte Carlo simulation method, as described in previous subsection. The results are presented for various values of parameters Qp, Pout,PU, and for a few values of mean channel power gains.
The simulation parameters that are used for calculations of the results presented in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 are shown in Table 1. In all simulations, the energy conversion efficiency coefficient is η = 0.9. The number of primary transmitters is set to N = 1 in all simulations, except in the case of the results presented in Figure 6, where it is considered that the number of primary transmitters is N = 2. The simulation results are obtained based on waveform sequences with L = 107 samples. As the estimated outage probability values for various input parameters are not lower than 10−4 for any set of system and channel parameters, the obtained estimations can be considered accurate with high confidence level [40].
The dependence of outage probability on the interference threshold Qp is examined, and the results are presented in Figure 4 for the case of TSR protocol. Scenarios with different values of primary network outage probability Pout,PU and mean channel power gain ΛR from the secondary relay to the primary receiver are considered. The mean channel power gain ΛS from the secondary source to the primary receiver is set to 0.01, the energy harvesting ratio is 0.5, the mean channel power gains are ΩIR = 0 dB, ΩID = −10 dB, ΩSR = 10 dB, and ΩRD = 10 dB, while the primary user transmitters have the unit power and the outage threshold is γth = −5 dB. By increasing the interference threshold Qp, the outage probability decreases. This result can be explained by fact that the increase of the interference threshold Qp leads to higher permitted transmit power values at the secondary source, as well as more harvested energy and higher permitted power at the relay.
The outage probability values decrease for smaller values of Λ R , as the allowed transmit power of the secondary relay is higher. Moreover, for higher values of permitted primary network outage probability Pout,PU, the outage probability of the cognitive network is lower. The reason for this effect can be explained by the fact that with the higher maximal allowed value of outage probability (for the fixed interference threshold Qp), both the transmit power of the secondary source and the allowed power at the relay are increasing, and the received signal-to-interference ratios at both the secondary relay and the destination are higher. Similar results are obtained for PSR protocol using the same scenarios.
In the next figures, the achievable throughput of the energy-harvesting DF system with statistical CSI is analyzed. The achievable throughput depends on the applied protocol, so for TSR and PSR protocol, it can be expressed, respectively, as [24,41]
T O U T T S R = ( 1 α ) C O U T ,
T O U T P S R = C O U T ,
where COUT is outage capacity, which is defined as a maximum data rate that can be achieved in the channel with the outage probability P O U T ( γ t h ) = Pr ( γ < γ t h )
C O U T = 0.5 ( 1 P O U T ( γ t h ) ) l o g 2 ( 1 + γ t h ) .
More specifically, in the DF system, outage probability POUT is defined with Equation (22).
In Figure 5, the throughput is presented as a function of mean channel power gains ΩSR and ΩRD for both analyzed energy-harvesting protocols, where the outage probability of the primary link is set to Pout,PU = 0.05 and the permitted interference threshold is Qp = 10 dB. The impact of different values of the mean channel power gain Λ S from the secondary source to the primary receivers and the mean channel power gain Λ R from the secondary relay to the primary receivers is considered. The energy conversion efficiency coefficient is η = 0.9, while the time-splitting factor and power-splitting factor are equal to 0.5 for TSR and PSR protocols, respectively. The mean channel power gains in the interference channels from the primary transmitters to secondary relay and destination are ΩIR = 0 dB and ΩID =−10 dB, respectively. We analyzed the scenario where mean values of channel power gain in both source-to-relay and relay-to-destination links are equal and the outage threshold is γth = 8 dB. The obtained results demonstrate that in all considered cases, the throughput values increase with the raise of mean channel power gains ΩSR =RD when these values are below 10 dB, while the saturation effect occurs for the higher values of mean channel power gain.
Furthermore, the saturation occurs for higher ΩSR =RD when the mean channel power gain Λ S is higher, as the transmit power of the secondary relay is then dominantly determined by the interference constraint of the primary user. This can be explained by the fact that in this case, the maximal allowed transmit power of the secondary source is smaller, and the interference constraint at the relay becomes dominant only for higher mean channel power gain ΩSR. In accordance to the expectations, in all analyzed cases, greater throughput values are achieved in the case when PSR energy harvesting protocol is applied. It can be noticed that the throughput increases with the decrease of mean channel power gain Λ S , as in that case, the higher secondary source transmit power is permitted. Similarly, the throughput values increase for smaller values of Λ R , as in that case, the allowed transmit power of the secondary relay is higher. However, in the case when Λ S = 10 , the difference between cases of various Λ R values is diminishing, as outage performances are dominantly limited by the source-to-relay link.
The dependence of the cognitive system throughput on the time-splitting factor α and power-splitting factor ρ is presented in Figure 6, for different values of primary user’s interference threshold Qp and number of primary transmitters N. The obtainable throughput is analyzed for both TSR and PSR protocols, with the energy conversion efficiency coefficient η = 0.9. The maximal permitted probability of the primary link outage is set to Pout,PU = 0.05, the outage threshold of the cognitive system is γth = 8 dB. In this scenario, the mean channel power gains are equal to ΩIR = 0 dB, ΩID = −10 dB, ΩSR = 10 dB, and ΩRD = 10 dB. Moreover, in the case when TSR protocol is applied, one can identify the optimal value of time-splitting factor α when the throughput is maximized. This value is achieved for values of α < 0.05, while a higher increase of time-splitting factor leads to the throughput decrease. On the other hand, in the case of energy-harvesting PSR protocol, the throughput increases with the power-splitting factor for small values of ρ, as then, more power is dedicated to energy harvesting, and the transmit power of the relay is increasing. However, in the case of further increase, the saturation occurs, as the maximal relay transmit power is determined with the more strict restriction imposed by the primary network and the average SIR at the relay is independent of ρ.
Furthermore, better system performances are achieved with a smaller number of transmitters in the primary network for both energy-harvesting protocols. This can be explained by the fact that although the interference from the primary network is used for energy harvesting, the degrading effect of the interference is more dominant. Finally, from the obtained results, it can be concluded that the throughput always increases with the raise of the allowed interference threshold Qp. A bigger value of Qp increases both the total input power PS as well as the permitted transmit power at the relay, therefore resulting in better system throughput.
In Figure 7, the throughput is shown as the function of ΩSR and ΩRD for both analyzed energy harvesting protocols. We considered scenarios with different values of mean channel power gains ΩIR and ΩID in the interference channels from the primary transmitter to both relay and destination nodes, respectively. The outage probability of primary link is set to 0.05 and the interference threshold is Qp = 10 dB (the time-switching and power-splitting factors are α = ρ = 0.5 and the outage threshold is γth = 8 dB). The mean value of the mean channel power gain from both the secondary source and relay to the primary receivers are Λ S = Λ R = 1 . The best system performances are achieved for the smallest values of ΩID and ΩIR, for both analyzed protocols. Furthermore, the increase of ΩIR and ΩID has a negative influence on the system performances. In the case of PSR protocol, obtained results are better when ΩIR = 0 dB and ΩID = −10 dB than when ΩIR = −10 dB, ΩID = 0 dB, while in the case of TSR protocol, the difference is smaller, and these curves are almost overlapping. This effect can be explained by the fact that the interference at the destination has only a degrading effect, whereas at the relay interference, it also increases harvested energy. The throughput increases with the raise of ΩSR =RD, and in the region of high values, it is the same for all analyzed values of ΩIR and ΩID. That is because the transmit power of the energy-constrained relay node is determined by the harvested energy and the maximal power permitted due to the interference outage constraint. That effect is noticeable also in the case of PSR protocol for values of ΩIR and ΩID higher than 25 dB.
Throughput values are presented in Figure 8 as a function of the permitted interference threshold Qp for both energy-harvesting protocols (α = ρ = 0.2) and different values of primary network outage probability. The mean channel power gains are ΩIR = 0 dB, ΩID = −10 dB, ΩSR = 10 dB, ΩRD = 10 dB, and the threshold is γth = 8 dB. For low values of Qp, throughput values are very low, for both protocols and all analyzed values of primary network outage probability. For small values of Qp, the total transmit power is very low, the amount of harvested energy is limited, and the transmit power of the secondary relay is consequently low. By increasing the permitted interference threshold Qp, values of the system throughput are rising, as more energy is harvested and the obtainable transmit power of the secondary relay is higher. Moreover, the permitted values of the secondary relay transmit power is increasing for higher Qp as the interference constraint is relaxed. For higher values of Qp, better system performances are obtained with applied PSR protocol (than TSR). Moreover, better performances are achieved for higher values of permitted primary network outage probability Pout,PU. The reason for this effect can be explained by the fact that with higher maximal allowed value of outage probability (for the fixed interference threshold Qp), both the transmit power of the secondary source and the allowed power at the relay are increasing, and the received signal-to-interference ratios at the secondary relay and the destination are higher.
Although the used system model is widely accepted and utilized in various scenarios of cognitive and energy-harvesting systems, it has some weaknesses and limitations. As explained in assumption A3, our analysis is limited to the case when all PU transmitters are grouped in one cluster, and all PU receivers are grouped in the other cluster (as a special case, it is valid for the case of single PUTX and single PURX). In addition, we assumed in A4 that the line-of-sight propagation component is not present and the shadowing effect can be neglected, although this assumption is not reasonable in some propagation scenarios. More importantly, the energy-harvesting model is an idealization of the realistic process.
The realistic energy-harvesting model has the minimum amount of energy necessary for the harvesting process, as SR consumes energy in the process of energy harvesting [42]. Therefore, we will obtain simulation results for different values of minimum threshold constraints of the input power at the SR node, Po, which is needed for achieving energy harvesting gain.
Throughput values estimated by using Monte Carlo simulation method are presented in Figure 9, as a function of the permitted interference threshold Qp for both energy-harvesting protocols (α = ρ = 0.5) and different values of Po. The outage threshold is γth = 8 dB, while the mean channel power gain values are set to ΩIR = 0 dB, ΩID = −10 dB, ΩSR =10 dB, and ΩRD =10 dB. For low values of Qp, throughput values are very low, for both protocols and all analyzed values of Po. This can be explained by the fact that for small values of Qp, the total input power at the SR is very low, so it is difficult to reach a defined minimum threshold of input power needed for achieving energy-harvesting gain. By increasing the permitted interference threshold Qp, the total input power at the SR is also increasing, more energy is harvested, and the obtainable transmit power of the SR is higher, resulting in the improved system throughput values. Moreover, the permitted value of the SR transmit power is increasing for higher Qp as the interference constraint is relaxed. For higher values of Qp, better system performances are obtained with applied PSR protocol (than TSR).
It is obvious that better performances are achieved for lower values of minimum threshold Po, because the probability that the total input power of SR exceeds these values is higher. In the case when the total input power of the SR is below Po, there is no energy harvesting gain at the SR, and the output power of the SR is equal to zero. For values of minimum threshold constraints of the SR input power equal to 20 dBm, system performances are comparable with the ones obtained in the case when there is no threshold constraints (i.e., Po = 0 W). It is important to highlight that all considered values are significantly higher than values previously used in the literature (Po = −12 dBm or 64 μW in [42]).

5. Conclusions

In this paper, we analyzed a cognitive DF relaying system with multiple primary transmitters and primary receivers. The self-sustainable relay harvests energy from both cognitive and primary transmitters based on TSR and PSR energy-harvesting protocols. The transmission among nodes in the cognitive relaying network is constrained by the maximal tolerable outage probability and the interference threshold permitted by the primary user. The closed-form expressions for outage probability and achievable throughput have been derived for the interference-limited case and the Rayleigh propagation environment. All derived analytical results are confirmed by an independent simulation method.
It has been shown that a higher value of the interference threshold increases the permitted transmit power at the SS, enabling a larger amount of harvested energy and higher obtainable transmit power of the SR. As the increased interference threshold value also results in higher allowable transmit power at the SR (for the given amount of harvested energy), the resulting outage probability decreases, and system throughput improves. The performances of the proposed system are also improved for higher permitted primary network outage probability. The throughput values of the system increase with the raise of the mean channel power gains from the SS to the SR and from the SR to the SD when these values are below a certain value, while the saturation effect occurs for the higher values of the mean channel power gain. Furthermore, the decreased values of the mean channel power gain from the SS to the PU’s receivers and from the SR to the PU’s receivers result in better system performances. The opposite happens with the increase of the mean channel power gains from the PU’s transmitters to the SR and from the PU’s transmitters to the SD. In the case when TSR protocol is applied, one can identify the optimal value of the time-splitting factor α when the throughput is maximized. After that value, a further increase of the time-splitting factor leads to the throughput decrease. On the other hand, in the case of the energy-harvesting PSR protocol, the throughput increases with the power-splitting factor for small values of ρ, while further increase leads to saturation.
Finally, the provided analysis and derived closed-form expressions open up new research topics for further investigations. Due to limited computational capabilities of the proposed network and limited energy of the SR, it is a challenging task to provide secure communication in the secondary network. In such a scenario, physical layer security has some benefits when compared to complex cryptography techniques that ensure a high level of security. Therefore, the first task in the future work will be to determine the secrecy capacity of the secondary link in the scenario analyzed in this paper. Furthermore, the presented analysis can be extended to the case when the propagation environment is modeled with Nakagami-m distribution. This is a generalization of the considered scenario; as for the case when the fading parameter is m = 1, it corresponds to the analysis of the Rayleigh fading environment presented in this paper. Furthermore, as the analyzed statistical channel model encompasses the path-loss effect, the impact of network geometry and distances among nodes can be investigated based on the results provided in the paper. This type of analysis can offer some valuable guidelines for the optimal design of the cognitive energy-harvesting relaying system that maximizes both primary and secondary system performances. However, the considered model does not include the shadowing effect, so the analysis can be extended to the case when shadowing is also modeled (using log-normal or generalized Gamma fading distribution). Future work can also analyze some specific scenarios where channel power gains from the same group of users differ due to a non-uniform distribution of small-sized obstacles in the environment.

Author Contributions

Conceptualization, N.K. and V.B.; methodology, N.K., V.B. and P.I.; software, N.K.; validation, N.K., V.B. and P.I.; formal analysis, N.K. and V.B.; investigation, N.K.; resources, V.B.; data curation, N.K.; writing—original draft preparation, N.K.; writing—review and editing, V.B. and P.I.; visualization, N.K. and P.I.; supervision, P.I. and V.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this part, we derive the mathematical statistics for random variables v i w.
As random variable is defined as ν = γ S R / I R , using [43] (7.23), we obtain the corresponding PDF as
p v ( v ) = 0 x p γ S R ( v x ) p I R ( x ) d x .
Furthermore, by substituting corresponding PDF expressions, the following PDF is obtained
p ν ( ν ) = N Ω I R Ω S R 1 ( Ω I R ν Ω S R + 1 ) N + 1
On the basis of the previous expression, the CDF expression for the random variable v is obtained, and it is given by Equation (31).
Following the same approach, we obtain the PDF and the CDF expressions for the random variable w = γ R D / I D and they are given, respectively, as
p w ( w ) = N Ω I D Ω R D 1 ( Ω I D w Ω R D + 1 ) N + 1 ,
Pr { w γ t h } = 1 1 ( 1 + Ω I D Ω R D γ t h ) N

Appendix B

To solve the integral I2, the joint PDF for u = ( γ SR + I R ) and v = γ SR / I R is needed. With the help of p u , v ( u , v ) = | J | p γ SR ( u v v + 1 ) p I D ( u v + 1 ) [43] (pp. 7–37), where | J | = u ( v + 1 ) 2 is Jacobian determinant, and replacing appropriate PDFs, the joint PDF is found as
p u , v ( u , v ) = 1 Γ ( N ) Ω I R N Ω S R u N ( v + 1 ) N + 1 exp ( 1 Ω S R u v v + 1 ) exp ( 1 Ω I R u v + 1 ) .
We first solve integral over ν, γ t h p u , v ( u , v ) d v , by using the substitution t = ( v + 1 ) 1 and by using [37] (3.351.2) and following expression is obtained
γ t h p u , v ( u , v ) d v = ( 1 Ω I R 1 Ω S R ) N Γ ( N ) Ω I R N Ω S R exp ( u Ω S R ) ( Γ ( N ) Γ ( N , ( 1 Ω I R 1 Ω S R ) u 1 + γ t h ) ) ,
where Γ(⋅,⋅) is the incomplete Gamma function [38] (8.350.2).
Further, integration over u is performed by using serial representation of incomplete Gamma function [37] (3.351.2) and applying [38] (3.351.1)
0 γ t h c w γ t h p u , v ( u , v ) d v d u = 1 Ω I R N ( 1 Ω I R 1 Ω S R ) N ( 1 e γ t h Ω S R c w k = 0 N 1 ( 1 + γ t h ) Ω I R k ! ( γ t h Ω I R + Ω S R ) ( Ω S R Ω I R γ t h Ω I R + Ω S R ) k ( Γ ( 1 + k ) Γ ( 1 + k , γ t h Ω I R + Ω S R Ω S R Ω I R ( 1 + γ t h ) γ t h c w ) ) ) .
The last, the integral I2 is given by
I 2 = 1 Ω I R N ( 1 Ω I R 1 Ω S R ) N × γ t h q ( 1 e γ t h Ω S R c w k = 0 N 1 ( 1 + γ t h ) Ω I R k ! ( γ t h Ω I R + Ω S R ) ( Ω S R Ω I R γ t h Ω I R + Ω S R ) k × ( Γ ( 1 + k ) Γ ( 1 + k , γ t h Ω I R + Ω S R Ω S R Ω I R ( 1 + γ t h ) γ t h c w ) ) ) p w ( w ) d w .
Integral I2 can be written in the form as in Equation (33), where expressions I2 k, k = 1, 2, 3 are given as
I 21 = N Ω I D Ω R D γ t h q ( Ω I D w Ω R D + 1 ) N 1 d w = ( Ω I D Ω R D γ t h q + 1 ) N ,
I 22 = γ t h q e γ t h Ω S R c w p w ( w ) d w = N Ω I D Ω R D γ t h q e γ t h Ω S R c w ( Ω I D w Ω R D + 1 ) N 1 d w = = N Ω I D Ω R D e γ t h c Ω S R Ω I D Ω R D i = 0 N 1 ( N 1 i ) ( 1 ) N 1 i ( Ω I D Ω R D ) N 1 i × ( ( Ω I D Ω R D ) i N E i ( 1 i + N , γ t h c Ω S R Ω I D Ω R D ) ( Ω I D Ω R D + q γ t h ) i N E i ( 1 i + N , ( Ω I D Ω R D + q γ t h ) γ t h c Ω S R ) ) ,
I 23 = γ t h q Γ ( 1 + k , ξ w ) p w ( w ) d w = k ! N Ω I D Ω R D j = 0 k 1 j ! ( ξ ) j e ξ Ω I D Ω R D i = 0 N + j 1 ( N + j 1 i ) ( Ω I D Ω R D ) N + j 1 i × ( ( Ω I D Ω R D ) i N E i ( 1 i + N , Ω I D Ω R D ξ ) ( Ω I D Ω R D + q γ t h ) i N E i ( 1 i + N , ( Ω I D Ω R D + q γ t h ) ξ ) ) ,
where ξ = γ t h Ω I R + Ω S R Ω S R Ω I R ( 1 + γ t h ) γ t h c .

References

  1. Hu, F.; Chen, B.; Zhu, K. Full spectrum sharing in cognitive radio networks toward 5G: A survey. IEEE Access 2018, 6, 15754–15776. [Google Scholar] [CrossRef]
  2. Khan, A.A.; Rehmani, M.H.; Rachedi, A. Cognitive-radio-based internet of things: Applications, architectures, spectrum related functionalities, and future research directions. IEEE Wirel. Commun. 2017, 24, 17–25. [Google Scholar] [CrossRef]
  3. Goldsmith, A.; Jafar, S.A.; Maric, I.; Srinivasa, S. Breaking spectrum gridlock with cognitive radios: An information theoretic perspective. Proc. IEEE 2009, 97, 894–914. [Google Scholar] [CrossRef]
  4. Ghasemi, A.; Sousa, E.S. Fundamental limits of spectrum sharing in fading environments. IEEE Trans. Wirel. Commun. 2007, 6, 649–658. [Google Scholar] [CrossRef]
  5. Musavian, L.; Aissa, S. Capacity and power allocation for spectrum sharing communications in fading channels. IEEE Trans. Wirel. Commun. 2009, 8, 148–156. [Google Scholar] [CrossRef]
  6. Suraweera, H.A.; Smith, P.J.; Shafi, M. Capacity limits and performance analysis of cognitive radio with imperfect channel knowledge. IEEE Trans. Veh. Technol. 2010, 59, 1811–1822. [Google Scholar] [CrossRef] [Green Version]
  7. Jarrouj, J.; Blagojevic, V.; Ivanis, P. Outage probability of SINR for underlay cognitive radio systems in Nakagami fading. Frequenz 2014, 68, 563–572. [Google Scholar] [CrossRef]
  8. Kim, H.; Wang, H.; Lim, S.; Hong, D. On the impact of outdated channel information on the capacity of secondary user in spectrum sharing environments. IEEE Trans. Wirel. Commun. 2012, 11, 284–295. [Google Scholar] [CrossRef]
  9. Lim, S.; Wang, H.; Kim, H.; Hong, D. Mean Value-Based Power Allocation without Instantaneous CSI Feedback in Spectrum Sharing Systems. IEEE Trans. Wirel. Commun. 2012, 11, 874–879. [Google Scholar] [CrossRef]
  10. Jarrouj, J.; Blagojevic, V.; Ivanis, P. Outage Probability and Ergodic Capacity of Spectrum-Sharing Systems with MRC Diversity. Frequenz 2016, 70, 157–171. [Google Scholar] [CrossRef]
  11. Guan, X.; Yang, W.; Cai, Y. Outage Performance of Statistical CSI Assisted Cognitive Relay with Interference from Primary User. IEEE Commun. Lett. 2013, 17, 1416–1419. [Google Scholar] [CrossRef]
  12. Kabiri, C.; Zepernick, H.; Tran, H. Outage probability of a cognitive cooperative relay network with multiple primary users under primary outage constraint. In Proceedings of the 2016 International Conference on Advanced Technologies for Communications (ATC), Hanoi, Vietnam, 12–14 October 2016; pp. 38–42. [Google Scholar] [CrossRef]
  13. Lee, D. Performance analysis of scheduled STBC using statistical CSI under multiple interferers with different power in CR-MIMO systems. Phys. Commun. 2018, 27, 125–132. [Google Scholar] [CrossRef]
  14. Ulukus, S.; Yener, A.; Erkip, E.; Simeone, O.; Zorzi, M.; Grover, P.; Huang, K. Energy Harvesting Wireless Communications: A Review of Recent Advances. IEEE J. Sel. Areas Commun. 2015, 33, 360–381. [Google Scholar] [CrossRef] [Green Version]
  15. Yildiz, F. Potential Ambient Energy-Harvesting Sources and Techniques. J. Technol. Stud. 2009, 35, 40–48. [Google Scholar] [CrossRef]
  16. Paing, T.; Shin, J.; Zane, R.; Popovic, Z. Resistor emulation approach to low-power RF energy harvesting. IEEE Trans. Power Electron. 2008, 23, 1494–1501. [Google Scholar] [CrossRef]
  17. Powercast-Overview-rf-Energy-Harvesting-and-Wireless-Power-for-Micropower-Applications-1-728/. Available online: https://www.powercastco.com/company/about/ (accessed on 10 August 2020).
  18. Lee, K.; Ko, J. RF-Based Energy Transfer Through Packets: Still a Dream? or a Dream Come True? IEEE Access 2019, 7, 163840–163850. [Google Scholar] [CrossRef]
  19. Huang, J.; Zhou, Y.; Ning, Z.; Gharavi, H. Wireless power transfer and energy harvesting: Current status and future prospects. IEEE Wirel. Commun. 2019, 26, 163–169. [Google Scholar] [CrossRef]
  20. Doan, T.X.; Hoang, T.M.; Duon, T.Q.; Ngo, H.Q. Energy Harvesting-Based D2D Communications in the Presence of Interference and Ambient RF Sources. IEEE Access 2017, 5, 5224–5234. [Google Scholar] [CrossRef] [Green Version]
  21. Park, S.; Kim, H.; Hong, D. Cognitive radio networks with energy harvesting. IEEE Trans. Wirel. Commun. 2013, 12, 1386–1397. [Google Scholar] [CrossRef]
  22. Lee, S.; Zhang, R.; Huang, K. Opportunistic wireless energy harvesting in cognitive radio networks. IEEE Trans. Wirel. Commun. 2013, 12, 4788–4799. [Google Scholar] [CrossRef] [Green Version]
  23. Blagojevic, V.M.; Cvetkovic, A.M.; Ivanis, P. Performance analysis of energy harvesting DF relay system in generalized-K fading environment. Phys. Commun. 2018, 28, 190–200. [Google Scholar] [CrossRef]
  24. Gu, Y.; Aissa, S. RF-Based Energy Harvesting in Decode-and-Forward Relaying Systems: Ergodic and Outage Capacities. IEEE Trans. Wirel. Commun. 2015, 14, 6425–6434. [Google Scholar] [CrossRef]
  25. Yang, S.; Lu, G.; Ren, Y. Optimal power allocation and relay location for DF energy harvesting relaying sensor networks. Sensors 2019, 19, 2326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Xu, H.; Gao, H.; Zhou, C.; Duan, R.; Zhou, X. Resource allocation in cognitive radio wireless sensor networks with energy harvesting. Sensors 2019, 19, 5115. [Google Scholar] [CrossRef] [Green Version]
  27. Cvetkovic, A.; Blagojevic, V.; Ivanis, P. Performance Analysis of Nonlinear Energy-Harvesting DF Relay System in Interference-Limited Nakagami-m Fading Environment. ETRI J. 2017, 39, 803–812. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, F.; Jing, T.; Huo, Y.; Jiang, K. Outage probability minimization for energy harvesting cognitive radio sensor networks. Sensors 2017, 17, 224. [Google Scholar] [CrossRef]
  29. Toan, H.; Bao, V.; Le, K. Performance analysis of cognitive underlay two-way relay networks with interference and imperfect channel state information. EURASIP J. Wirel. Commun. Netw. 2018, 53. [Google Scholar] [CrossRef] [Green Version]
  30. Prasad, B.; Dhar Roy, S.; Kundu, S. Performance of cognitive relay network with energy harvesting relay under imperfect CSI. Int. J. Commun. Syst. 2018, 31, e3549. [Google Scholar] [CrossRef]
  31. Hoang, T.M.; Van, N.L.; Nguyen, B.C.; Dung, L.T. On the performance of energy harvesting non-orthogonal multiple access relaying system with imperfect channel state information over Rayleigh fading channels. Sensors 2019, 19, 3327. [Google Scholar] [CrossRef] [Green Version]
  32. Nasir, A.; Zhou, X.; Durrani, S.; Kennedy, R. Relaying protocols for wireless energy harvesting and information processing. IEEE Trans. Wirel. Commun. 2013, 12, 3622–3636. [Google Scholar] [CrossRef] [Green Version]
  33. Duong, T.; Yeoh, P.L.; Bao, V.N.Q.; Elkashlan, M.; Yang, N. Cognitive relay networks with multiple primary transceivers under spectrum-sharing. IEEE Signal Process. Lett. 2012, 19, 741–744. [Google Scholar] [CrossRef]
  34. Liu, Y.; Mousavifar, S.A.; Deng, Y.; Leung, C.; Elkashlan, M. Wireless Energy Harvesting in a Cognitive Relay Network. IEEE Trans. Wirel. Commun. 2016, 15, 2498–2508. [Google Scholar] [CrossRef] [Green Version]
  35. Ma, S.; Yang, Y.; Sharif, H. Distributed MIMO technologies in cooperative wireless networks. IEEE Commun. Mag. 2011, 49, 78–82. [Google Scholar] [CrossRef]
  36. Khalifeh, A.; Abid, H.; Darabkh, K.A. Optimal Cluster Head Positioning Algorithm for Wireless Sensor Networks. Sensors 2020, 20, 3719. [Google Scholar] [CrossRef]
  37. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; National Bureau of Standards: Washington, DC, USA, 1972; Volume 55, p. 319. [Google Scholar]
  38. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series and Products, 5th ed.; Academic Press Inc.: San Diego, CA, USA, 1994. [Google Scholar]
  39. Zheng, Y.R.; Xiao, C. Simulation models with correct statistical properties for Rayleigh fading channels. IEEE Trans. Commun. 2003, 51, 920–928. [Google Scholar] [CrossRef]
  40. Jeruchim, M.C.; Balaban, P.; Shanmugam, K.S. Simulation of Communication Systems; Plenum Press: New York, NY, USA, 1992. [Google Scholar]
  41. Goldsmith, A. Wireless Communications; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  42. Wang, S.; Xia, M.; Huang, K.; Wu, Y. Wirelessly Powered Two-Way Communication with Nonlinear Energy Harvesting Model: Rate Regions Under Fixed and Mobile Relay. IEEE Trans. Wirel. Commun. 2017, 16, 8190–8204. [Google Scholar] [CrossRef]
  43. Papoulis, A. Probability, Random Variables, and Stochastic Processes; McGraw-Hill: New York, NY, USA, 1991. [Google Scholar]
Figure 1. Model of the underlay cognitive relay network.
Figure 1. Model of the underlay cognitive relay network.
Sensors 21 03727 g001
Figure 2. Transmission block structure in the TSR scheme.
Figure 2. Transmission block structure in the TSR scheme.
Sensors 21 03727 g002
Figure 3. Transmission block structure in the PSR scheme.
Figure 3. Transmission block structure in the PSR scheme.
Sensors 21 03727 g003
Figure 4. Outage probability POUT vs. Qp for TSR protocol.
Figure 4. Outage probability POUT vs. Qp for TSR protocol.
Sensors 21 03727 g004
Figure 5. Throughput for different values of mean channel power gain in a cognitive relay system.
Figure 5. Throughput for different values of mean channel power gain in a cognitive relay system.
Sensors 21 03727 g005
Figure 6. Throughput vs. energy harvesting ratio.
Figure 6. Throughput vs. energy harvesting ratio.
Sensors 21 03727 g006
Figure 7. The impact of the interference power at the secondary source and the relay on the throughput values for various ΩSR = ΩRD.
Figure 7. The impact of the interference power at the secondary source and the relay on the throughput values for various ΩSR = ΩRD.
Sensors 21 03727 g007
Figure 8. The achievable throughput vs. interference threshold Qp.
Figure 8. The achievable throughput vs. interference threshold Qp.
Sensors 21 03727 g008
Figure 9. Simulation results for the achievable throughput vs. interference threshold Qp and minimum threshold constraints of the input power at the SR, Po.
Figure 9. Simulation results for the achievable throughput vs. interference threshold Qp and minimum threshold constraints of the input power at the SR, Po.
Sensors 21 03727 g009
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParameterFigure 4Figure 5Figure 6Figure 7Figure 8Figure 9
Pout,PU0.01–0.10.050.050.050.01–0.10.05
Qp [dB]−20–20105–1510−15–20−15–20
Λ R 0.1–20.1–11111
Λ S 0.010.1–101111
IR [dB]000−10–000
ID [dB]−10−10−10−10–0−10−10
SR [dB]10−10–2010−10–301010
RD [dB]10−10–2010−10–301010
γth [dB]−588888
η0.90.90.90.90.90.9
α0.50.50.01–0.990.50.20.5
ρ0.50.50.01–0.990.50.20.5
N111–2111
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kozić, N.; Blagojević, V.; Ivaniš, P. Performance Analysis of Underlay Cognitive Radio System with Self-Sustainable Relay and Statistical CSI. Sensors 2021, 21, 3727. https://doi.org/10.3390/s21113727

AMA Style

Kozić N, Blagojević V, Ivaniš P. Performance Analysis of Underlay Cognitive Radio System with Self-Sustainable Relay and Statistical CSI. Sensors. 2021; 21(11):3727. https://doi.org/10.3390/s21113727

Chicago/Turabian Style

Kozić, Nadica, Vesna Blagojević, and Predrag Ivaniš. 2021. "Performance Analysis of Underlay Cognitive Radio System with Self-Sustainable Relay and Statistical CSI" Sensors 21, no. 11: 3727. https://doi.org/10.3390/s21113727

APA Style

Kozić, N., Blagojević, V., & Ivaniš, P. (2021). Performance Analysis of Underlay Cognitive Radio System with Self-Sustainable Relay and Statistical CSI. Sensors, 21(11), 3727. https://doi.org/10.3390/s21113727

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