1

Outage Analysis of Energy Harvesting based
Opportunistic Cooperative Communication Systems

arXiv:1801.06854v1 [] 21 Jan 2018

Yanju Gu

Abstract
Wireless energy harvesting constitutes an effective way to prolong the lifetime of wireless networks.
In this paper, an opportunistic decode-and-forward cooperative communication system is investigated,
where the energy-constrained relays harvest energy from the received information signal and the cochannel interference (CCI) signals and then use that harvested energy to forward the correctly decoded
signal to the destination. Different from conventional relaying system with constant energy supplier,
the transmission power of the energy-constrained relay depends on the available energy that harvested,
which is a random process. Best relay selection that takes into account all the impacting ingredients
on the received signal quality at the destination is deployed. The exact closed-form expression of the
outage probability is derived, and the optimal value of the energy harvesting ratio for achieving minimum
outage is numerically investigated. In addition, the impacts of the CCI signals on the system’s outage
and diversity performances are analyzed. It is shown that the proposed relaying scheme can achieve full
diversity order equal to the number of relays without the need of fixed power supplier.

I. I NTRODUCTION
Energy harvesting has become a strong solution for providing green energy and is widely
adopted in green communications [1]–[4]. The energy is usually harvested from solar and wind
power and then transmitted to the power grid [5] with the knowledge of power states [6], [7].
In recent years, radio frequency (RF) based energy harvesting techniques [8], [9] have received
more and more attention as this RF energy is already pervasive in the wireless communication
network, such as those ambient RF from TV and cellular communications. In cooperative
communication networks, the intermediate relays are randomly deployed for forwarding signals
to the destination. These devices like sensors are typically equipped with batteries with limited
operation life. Replacing batteries for such devices is not an easy task as they may be located
in hostile environments. Therefore, equipping relays with the capability of RF energy harvesting
can prolong the lifetime of the network [10].

January 23, 2018

DRAFT

2

The relays can harvest energy from the received RF signal which is a superposition of the
desired signal and the CCI signals. The cooperative decode-and-forward (DF) relaying systems
with single energy-constrained relay have been studied in [3], [11], where co-channel interference
(CCI) signals as useful energy are taking into account. In this paper, we take a step further by
considering multi-relay cooperation, which can provide further spatial diversity. In conventional
DF multi-relaying systems with constant power supply for the relay [12], once the relay decodes
the information correctly, the outage performance at the destination only depends on the channel
quality between the relay and the destination. However, it is much more complicated for the
energy harvesting based relaying system, since the energy constrained relay is operating with
the random energy that harvested from the received information and interference signals. In this
work, the exact closed-form expression of the outage probability for the energy-harvesting-based
multi-relay DF cooperative system is derived, where best relay selection is deployed. The optimal
value of the energy harvesting ratio for achieving minimum outage is numerically investigated.
In addition, the impacts of the CCI signals on the system’s outage and diversity performances
are also analyzed.
II. E NERGY-H ARVESTING BASED R ELAYING
A. System and Channel Models
A cooperative DF relaying system is considered, where the source S communicates with
the destination D through the help of multiple energy-constrained intermediate relaying nodes
(R1 , R2 , · · · , RL ). Each node is equipped with a single antenna and operates in the half-duplex
mode in which the node cannot simultaneously transmit and receive signals in the same frequency
band. Both, the first hop (source-to-relay) and the second hop (relay-to-destination), experience
independent Rayleigh fading with the complex channel fading gains given by hi ∼ CN (0, Ωhi )
and gi ∼ CN (0, Ωgi ), respectively. The channels follow the block-fading model in which
the channel remains constant during the transmission of a block and varies independently
from one block to another. The channel state information is only available at the receiver. As
aforementioned, the system operates in the presence of external interferers. Specifically, we
assume that there is an aggregate CCI signal affecting the ith relay. The channel fading gain
between the interferer and the ith relay, denoted βi , is modeled as βi ∼ CN (0, Ωβi ). The desired
channels and the interference channel are supposed to be independent from each other.

January 23, 2018

DRAFT

3

B. Wireless Energy Harvesting at the Relay
The power-splitting based protocol, where a portion of the received power is utilized for
energy harvesting and the remaining power is used for information processing, is adopted at the
relay node. P is the received power of the signal at the relay and θ, with 0 ≤ θ ≤ 1, is the
fraction of power that the relay harvests from the received interference and information signal.
The remaining power is (1 − θ)P , for information transmission from the relay to the destination.
In the first-hop phase, the source S transmits signal s with power PS to the relays. We consider
a pessimistic case in which power splitting only reduces the signal power, but not to the noise
power, which can provide a lower bound for relaying networks in practice. Accordingly, the
received signal at the ith relay for information detection is given by
ySRi =

p
p
(1 − θ)PS hi s + (1 − θ)Pi βi si + nRi ,

(1)

where si and Pi denote the signal and the corresponding power, respectively, from the interferer
at the ith relay, and nRi is the additive white Gaussian noise (AWGN) at the ith relay with zero
mean and variance σR2 .
According to (1), the received signal-to-interference-plus-noise ratio (SINR) at the relay is
given by
γSRi =
where γhi ,

(1−θ)PS
σ2
R

|hi |2 and Ii ,

γhi
(1 − θ)PS |hi |2
=
,
2
2
σR + (1 − θ)Pi |βi |
1 + Ii
(1−θ)Pi
|βi |2 .
2
σR

(2)

The relay harvests energy from the received

information signal and the interference signal for a duration of T /2 at each block, and thus, the
harvested energy at the ith relay is obtained as


2
2
EHi = ηθ PS |hi | + Pi |βi | T /2,

(3)

where η is the energy conversion efficiency of the relay with value varying from 0 to 1 depending
upon the harvesting circuitry. It is assumed that there is no additional device like an energy buffer
to store the harvested energy (Harvest-Use) which can decrease the complexity of the energy
harvesting nodes, so all the energy collected during the harvesting phase is consumed by the relay.
Since the processing power required by the transmit/receive circuitry at the relay is generally
negligible compared to the power used for signal transmission [13], [14], here we suppose that

January 23, 2018

DRAFT

4

all the energy harvested from the received signal is consumed by the relay for forwarding the
information to the destination node. The transmission power of the ith relay is then given by
P Ri =

EHi
ηθσR2
=
(γhi + Ii ) .
T /2
1−θ

(4)

C. Selective Transmission
Let L = {1, 2, . . . , L} denote the set of all the L relays and S represents the decoding subset
consisted of those relays that are able to decode the source message. That is,
o
n
S , i ∈ L : γSRi ≥ γth ,

(5)

where γth is a pre-defined threshold.
In the second-hop phase, only a single relay among the relays belonging to subset S is allowed
to transmit the signal. More specifically, from that decoding subset, the relay with the maximum
signal-to-noise ratio (SNR) retransmits the information. That is, î = arg maxi∈S {γRi D }. Such
technique is referred to as opportunistic relay selection. Note that the relay selection process
require centralized processing and how to extend it to distributed processing [15]–[17] is still
an open problem. Therefore, the received SNR at the destination, given set S, can be written as
γR D =
î

PR |gî |2
î

σD2

ηθ σR2
=
|gî |2 (γhî + Iî ),
2
(1 − θ) σD
{z
}
|

(6)

,Wî

where σD2 is the variance of AWGN at the destination. The defined random variable Wî follows
the same distribution as of |gî |2 .
Note that, in contrast to traditional DF relaying system with no rechargeable nodes, the
transmission power PRi at the relay in the energy harvesting system is not a constant but a
random variable, which depends on the replenished energy from the interference and information
signals. Therefore, the distribution of the received SNR at the destination is determined not only
by the distribution of the relay-to-destination channel power gain |gî |2 , but also by the distribution
of the information and interference signal power, i.e., γhî and Iî .
III. P ERFORMANCE E VALUATION
A. Outage Analysis
As an important performance measure of wireless systems, outage probability is defined
as the probability that the instantaneous output SNR falls below a pre-defined threshold γth .
January 23, 2018

DRAFT

5

Mathematically speaking, Pout (γth ) = Pr (γ < γth ). This SNR threshold guarantees the minimum
QoS requirement of the destination users. In the DF relaying system under study, the outage
probability at the destination, given set S, is expressed as
Pr{γR D < γth | |S| = l}
î
h
il
= Pr{γRi D < γth γSRi ≥ γth }
#l
"
Pr{γRi D < γth , γSRi ≥ γth }
,
=
Pr{γSRi ≥ γth }

(7)
(8)

where |S| denotes the number of relays in the decoding set S. Note that we assume a homogenous
scenario where the desired channels as well as the interference channels are independent and
identically distributed (i.i.d.), i.e., Ωhi = Ωh , Ωgi = Ωg , Ωβi = Ωβ , Pi = PI for i = 1, 2, . . . , L,
and hence, the outage probability is independent of the combination of relays forming the
decoding subset. Expression (7) is obtained by using the first principle of order statistics [18].
Since γRi D and γSRi are not known to be independent, this conditional probability is rewritten
as (8), in which, the probability in the denominator is given by
Pr{γSRi ≥ γth } =1 − Pr{γSRi < γth }
−1



γ̄β
γth
exp −
,
= 1 + γth
γ̄h
γ̄h
where γ̄h =

(1−θ)PS
2
σR

Ωh and γ̄β =

(1−θ)PI
Ωβ .
2
σR

(9)

The calculation of the probability in the numerator

of (8) involves three random variables, γhi , γgi and Ii ; therefore, it is hard to get the result
directly. By taking into account the physical meaning of this joint probability, it can be divided
into two parts to simplify the derivation, that is,
Pr{γRi D < γth , γSRi ≥ γth }

= Pr γRi D 1C < γth − Pr{γSRi < γth },

(10)

where 1C is the indicator random variable for the set C = {γSRi ≥ γth }, i.e., 1C = 1 if γSRi ≥ γth ,
otherwise, 1C = 0.
Theorem 1: The first probability in (10) is given by

 γ

1
γ 
Pr γRi D 1C < γth =1 −
γ̄h Γ 1, th ; th − b−1
γ̄h − γ̄β
γ̄h γ̄h γ̄g

 

aγth
bγth 
× exp
Γ 1, bγth ;
,
1 + γth
γ̄g

January 23, 2018

(11)

DRAFT

6

σ2

ηθ R
and γ̄g = 1−θ
Ωg . Γ(a, x; b) is the generalized incomplete
σ2
R ∞ a−1 D
Gamma function defined by Γ(a, x; b) , x t exp(−t − bt−1 )dt.

where a ,

1
γ̄β

−

1
,
γ̄h

b,

1
γ̄h

+

a
1+γth

Proof: From (6), we get γRi D 1C = Wi (γhi + Ii )1C . Define Z , (γhi + Ii )1C , then the
cumulative distribution function (CDF) of Z is given by
Z Z
FZ (z) = Pr{Z < z} =
fγhi ,Ii (x, y)dxdy,

(12)

x,y∈A

where the set A =

n
x + y < z,

o
> γth , x ≥ 0, y ≥ 0 . After some set manipulations, we

x
1+y

have A =
6 ∅ if and only if z > γth . Since γhi and Ii are independent, we get the joint distribution
fγhi ,Ii (x, y) = fγhi (x)fIi (y). Then, after some straightforward algebraic derivations, we obtain
FZ (z) = 1Z

Z

z−γ
th
1+γ
th

Z

z−y

(1+y)γth

0

fγhi (x)fIi (y)dxdy.

(13)

Both γhi and Ii are of exponential distributions. Integrate with respect to x and y yielding the
CDF of Z, where [19, Eq.(3.351.1)] was used. Then the probability density function (PDF) of
Z, fZ (z), follows directly from differentiating FZ (z) with respect to z, and is given by



z
z − γth
1
exp(− ) 1−exp −a
,
fZ (z) = 1Z
γ̄h − γ̄β
γ̄h
1 + γth
1
γ̄β

where a ,

−

1
γ̄h

(14)

and 1Z is the indicator random variable for the set Z = {z > γth }, i.e.,

1Z = 1 if z > γth , otherwise, 1Z = 0. Finally, we have
Pr γRi D 1C < γth


where b ,

1
γ̄h

Γ(a, x; b) ,

January 23, 2018

+

R∞
x

a
1+γth





γth
=EZ 1 − exp −
γ̄g Z


Z ∞
γth
exp −
=1 −
fZ (z) dz
γ̄g z
0

 γ
1
γth 
th
=1 −
γ̄h Γ 1,
;
− b−1
γ̄h − γ̄β
γ̄h γ̄h γ̄g

 

aγth
bγth 
× exp
Γ 1, bγth ;
,
1 + γth
γ̄g

(15)

and Γ(a, x; b) is the generalized incomplete Gamma function defined by

ta−1 exp(−t − bt−1 )dt, which completes the proof.

DRAFT

7

According to the theorem of total probability, the unconditional outage probability at the
destination, Pout (γth ), is expanded as
Pout (γth ) =

L
X

Pr{|S| = l} Pr{γR D < γth | |S| = l}
î

l=0

h
iL
= Pr{γSRi < γth }
#l
" 
L  
X
Pr γRi D 1C < γth
L
−1 ,
×
Pr{γ
<
γ
}
l
SR
th
i
l=0

(16)

where the following expression is used to obtain (16),
Pr{|S| = l}
 h
iL−l h
il
L
Pr{γSRi < γth }
Pr{γSRi ≥ γth } .
=
l

(17)

B. Numerical Results and Discussion
Numerical examples are presented and corroborated by simulation results to examine the
outage of the DF cooperative communication system, where the energy-constrained relays harvest
energy from the received information signal and the CCI signals. Hereafter, and unless stated
otherwise, the threshold γth is set to 5dB and the energy conversion efficiency η is set to 1.
To better evaluate the effect of the interference on the system’s outage, we define
average signal-to-interference ratio (SIR) and

PI Ωβ
2
σR

PS Ωh
PI Ωβ

as the

as the average interference-to-noise ratio

(INR) at the first hop.
Figure 1 shows the outage probability versus the energy harvesting ratio θ for different number
of relays, where the first-hop average SNR is 20dB, and there is no interference affecting the
relay. It is observed that the analytical results of (16) match perfectly the simulation results.
As the number of available relays increases, the outage performance of the energy harvesting
based relaying system improves. The convex feature of the curves is due to the fact that the
energy harvested for the second-hop transmission increases with increasing θ, which effectively
decreases the outage of the second hop and, accordingly, improves the outage performance of
the system. Meanwhile, as θ increases, more power is harvested for information transmission
and less power is left for information decoding which increases the outage of the first hop and
reduces the number of reliable relays for relay selection, therefore, the outage probability first
drops down until reaching a minimum and then increases up.

January 23, 2018

DRAFT

8

0

10

L=1, 2, 3, 4
−1

Outage Probability

10

−2

10

−3

10

Simulation
Analytical

−4

10

0

0.2

0.4

0.6

0.8

1

θ

Fig. 1. Outage probability versus the energy harvesting ratio θ for different number of relays, where the first-hop average SNR
is 20dB and there is no interference at the relays.
0

10

−1

10

L=1

Outage Probability

−2

10

−3

L=2

10

−4

10

−5

10

L=3

−6

10

−7

10

5

Simulation
Analytical without Interference
Analytical with INR=10dB
10

15

20
25
SNR(dB)

30

35

40

Fig. 2. Outage probability versus the first-hop average SNR for different number of relays, where the first-hop average INR is
10dB and the energy harvesting ratio θ is set to 0.6.
0

10

L=1
−1

10

L=2
−2

Outage Probability

10

−3

L=3

10

−4

10

−5

10

−6

10

−7

10

5

Simulation
Analytical without Interference
Analytical with SIR=20dB
10

15

20
25
SNR(dB)

30

35

40

Fig. 3. Outage probability versus the first-hop average SNR for different number of relays, where the first-hop average SIR is
20dB and the energy harvesting ratio θ is set to 0.6.

January 23, 2018

DRAFT

9

In order to clearly demonstrate the diversity gain of the system and the impacts of CCI signals
on the system’s performance, Figs. 2 and 3 illustrate the outage probability of the system versus
the first-hop average SNR for different number of relays under two different situations where
the first-hop average INR is 10dB and the first-hop average SIR is 20dB, respectively. From
Fig. 1, it is clearly seen that the minimum outage is achieved when the energy harvesting ratio
θ is of the value around 0.6, so the ratio θ is set to 0.6 here for performance evaluation. It
is shown from each plot in Figs. 2 and 3 that, even though the CCI signals can be utilized
as useful energy during the energy harvesting phase, the existence of interference at the relay
during the information decoding phase still limits the outage performance. It is also shown that
when there is no interference affecting the relay, the spatial diversity order, which indicates the
decreasing speed of the outage with respect to SNR, increases as the number of available relays
increases. In Fig. 2, it is seen that the diversity order remains the same when the interference
keeps constant, while in Fig. 3, the diversity order is reduced to zero and error floors appear in
the high-SNR regime due to the constant interference to signal ratio.
IV. C ONCLUSIONS
In this paper, an opportunistic decode-and-forward (DF) cooperative communication system
was studied, where the relays are energy-constrained and need to replenish energy from the
received information signal and the co-channel interference (CCI) signals, and then use that
harvested energy to forward the correctly decoded signal to the destination. Different from
traditional DF relaying system with no rechargeable nodes, the transmission power of the energy
constrained relay is not a constant anymore but a random variable depending on the variation
of available energy harvested from the received information and CCI signals at the relay. In
order to better evaluate the system performance, the exact closed-form expression of the outage
probability was derived. The optimal value of the energy harvesting ratio for achieving minimum
outage was numerically investigated. The proposed relaying scheme can achieve full diversity
order equal to the number of relays without the need of fixed power supply. By applying some
interference cancelation schemes at the information decoder, the system performance can be
further improved.
R EFERENCES
[1] B. Gurakan, O. Ozel, J. Yang, and S. Ulukus, “Energy cooperation in energy harvesting communications,” IEEE Trans.
Commun., vol. 61, no. 12, pp. 4884–4898, Dec. 2013.

January 23, 2018

DRAFT

10

[2] V. Raghunathan, S. Ganeriwal, and M. Srivastava, “Emerging techniques for long lived wireless sensor networks,” IEEE
Commun. Mag., vol. 44, no. 4, pp. 108–114, Apr. 2006.
[3] Y. Gu and S. Aissa, “RF-based energy harvesting in decode-and-forward relaying systems: Ergodic and outage capacities,”
IEEE Trans. Wireless Commun., vol. 14, no. 11, pp. 6425–6434, Nov. 2015.
[4] Z. Ding and H. V. Poor, “Multi-user swipt cooperative networks: Is the max-min criterion still diversity-optimal,” IEEE
Trans. Wireless Commun., vol. 15, no. 1, pp. 553–567, 2016.
[5] H. Farhangi, “The path of the smart grid,” IEEE Power and Energy Magazine, vol. 8, no. 1, pp. 18–28, 2010.
[6] J. Du, S. Ma, Y.-C. Wu, and H. V. Poor, “Distributed bayesian hybrid power state estimation with pmu synchronization
errors,” in Global Communications Conference (GLOBECOM), 2014 IEEE, 2014, pp. 3174–3179.
[7] ——, “Distributed hybrid power state estimation under pmu sampling phase errors,” IEEE Transactions on Signal
Processing, vol. 62, no. 16, pp. 4052–4063, 2014.
[8] J. S. Hoa, A. J. Yeha, E. Neofytoub, S. Kima, T. Tanabea, B. Patlollab, B. Beyguib, and A. Poon, “Wireless power transfer
to deep-tissue microimplants,” Proceedings of the National Academy of Sciences, vol. 111, no. 22, pp. 7974–7979, June
2014.
[9] L. Mateu and F. Moll, “Review of energy harvesting techniques and applications for microelectronics,” in Proc. SPIE
Circuits and Syst. II, 2005, pp. 359–373.
[10] Q. Wang, X. Liu, J. Du, and F. Kong, “Smart charging for electric vehicles: A survey from the algorithmic perspective,”
IEEE Communications Surveys & Tutorials, vol. 18, no. 2, pp. 1500–1517, 2016.
[11] Y. Gu and S. Aissa, “Interference aided energy harvesting in decode-and-forward relaying systems,” in Proc. 2014 IEEE
Int. Conf. Commun., June 2014, pp. 5378–5382.
[12] S. S. Ikki and M. H. Ahmed, “Performance analysis of cooperative diversity with incremental-best-relay technique over
rayleigh fading channels,” IEEE Transactions on Communications, vol. 59, no. 8, pp. 2152–2161, 2011.
[13] X. Zhou, R. Zhang, and C. K. Ho, “Wireless information and power transfer: Architecture design and rate-energy tradeoff,”
in Proc. 2012 IEEE Global Commun. Conf., Dec. 2012, pp. 3982–3987.
[14] B. Medepally and N. Mehta, “Voluntary energy harvesting relays and selection in cooperative wireless networks,” IEEE
Trans. Wireless Commun., vol. 9, no. 11, pp. 3543–3553, Nov. 2010.
[15] J. Du, S. Ma, Y.-C. Wu, S. Kar, and J. M. Moura, “Convergence analysis of distributed inference with vector-valued
Gaussian belief propagation,” arXiv preprint arXiv:1611.02010, 2016.
[16] ——, “Convergence analysis of belief propagation for pairwise linear Gaussian models,” IEEE Global Conference on
Signal and Information Processing (GlobalSIP), arXiv preprint arXiv:1706.04074, 2017.
[17] ——, “Convergence analysis of the information matrix in Gaussian belief propagation,” IEEE International Conference
on Acoustics, Speech and Signal Processing (ICASSP), arXiv preprint arXiv:1704.03969, 2017.
[18] J. Du and Y.-C. Wu, “Network-wide distributed carrier frequency offsets estimation and compensation via belief
propagation,” IEEE Transactions on Signal Processing, vol. 61, no. 23, pp. 5868–5877, 2013.
[19] I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, series and products, 7th ed.

January 23, 2018

Academic Press, 7th Ed., 2007.

DRAFT

