1

A New Interweave Cognitive Radio Scheme for
Out-band Energy Harvesting Systems

arXiv:1801.10494v1 [] 31 Jan 2018

Nikolaos I. Miridakis and Theodoros A. Tsiftsis, Senior Member, IEEE

Abstract—In this paper, a new interweave cognitive radio
(CR) transmission scheme is analytically presented for energy
harvesting (EH)-enabled transmitting nodes. EH is performed
via a dedicated power beacon (PB) base station through wireless
power transfer, which harvests energy to both primary and
secondary systems. The out-band mode-of-operation between PB
and the two heterogeneous systems is considered (i.e., data communication and EH are realized at different frequency bands).
The performance of the new interweave CR system is studied
under independent and not necessarily identically distributed
κ − µ shadowed faded channels. Representative numerical and
simulation results corroborate the effectiveness of the proposed
scheme and presented analysis, while some useful engineering
insights are provided.
Index Terms—Cognitive radio (CR), energy harvesting (EH),
interweave secondary transmission, performance analysis.

I. I NTRODUCTION

E

NERGY harvesting (EH) via wireless power transfer
can be used for enabling perpetual operation of wireless
devices. Recently, EH has attracted increasing interests, mainly
due to its capability of converting received RF signals into
electricity, which in turn is able to provide stable and controllable power to prolong the lifetime of low-power energyconstrained autonomous networks [1]–[3]. Typical examples
are machine-type communications such as the Internet of
Things, Internet of Everything, Smart X, etc., where connectivity rather than high throughput seems to be of prime importance. On another front, cognitive radio (CR) has emerged as
one of the most promising technologies to resolve the issue of
spectrum scarcity. The key enabling technology of CR relies
on the provision of capability to share the spectrum in an
opportunistic manner, such that the secondary transmissions
do not cause any harmful interference to the primary communication [4].
The volume of the aforementioned machine-type communication devices (e.g., sensors) is expected to be massive
and widespread worldwide to support various upcoming 5Gand-beyond applications [5]. Thereby, spectrum- and energyefficiency should be some of the most prominent building
blocks for the design of such systems. Consequently, CR and
EH can be jointly utilized to achieve this goal due to their
complementary benefits [6], [7].
Essentially, systems in which EH and data transfer occur
at different frequency bands are known as out-band systems,
N. I. Miridakis is with the Department of Computer Systems Engineering,
Piraeus University of Applied Sciences, 12244, Aegaleo, Greece (e-mail:
nikozm@unipi.gr).
T. A. Tsiftsis is with the School of Engineering, Nazarbayev University,
Astana 010000, Kazakhstan (e-mail: theodoros.tsiftsis@nu.edu.kz).

while those working at the same frequencies are known as inband systems. Obviously, out-band systems are interferencefree between the power and data transfer segments [3]. To date,
most related research works are focused either on the underlay
CR transmission mode (i.e., CR and primary communication
overlaps to one another using the same resources) and/or the
in-band EH operation (i.e., see [8] and relevant references
therein).
In this paper, we propose a new interweave CR system,
which operates under the presence of a primary one. We
consider the case when both systems use EH-enabled transmission via a certain power beacon (PB) base station, which
transmits in a dedicated frequency band (i.e., an out-band EH
transmission). Notably, the out-band transmission mode is less
complex than its corresponding in-band counterpart [9]. The
main benefits of the proposed scheme are: (i) no co-channel
interference is produced between the two systems; (ii) the
CR system is purely EH-enabled since only an EH circuit
and a micro-supercapacitor (as an energy storage device) are
required, which is a much more cost-effective solution than
rechargeable batteries [10]. Moreover, some useful performance metrics are derived in closed-form when the signals
undergo κ − µ shadowed fading so as to approach realistic
channel conditions [11], [12] along with some impactful
engineering insights.
Notation: |x| takes the absolute value of x while kxk is the
Euclidean norm of vector x. E[·] is the expectation operator
d
and Pr[·] returns probability. The symbol = means equality in
distribution. The functions fX (·), FX (·) and F X (·) represent
probability density function (PDF), cumulative distribution
function (CDF) and complementary CDF of the random variable (RV) X, respectively. Complex-valued Gaussian RVs with
mean µ and variance σ 2 are denoted as CN (µ, σ 2 ). Finally,
Γ(·) denotes the gamma function [13, Eq. (8.310.1)], Γ(·, ·)
is the upper incomplete gamma function [13, Eq. (8.350.2)],
and ψ(·) is the digamma function [13, Eq. (8.360.1)].
II. S YSTEM AND C HANNEL M ODEL
Consider a point-to-point (secondary) CR system, which
operates under the presence of a primary system, as illustrated
in Fig. 1a, where PT and PR denote the single-antenna
primary Tx and Rx, respectively. Similarly, ST and SR denote
the single-antenna secondary Tx and Rx, respectively. It is
assumed that the Tx nodes of both systems harvest energy by
a dedicated PB base station, which is in their close vicinity.
Afterwards, based on this energy, they can transmit their
data streams to their corresponding Rx nodes. Independent
and not necessarily identically distributed κ − µ shadowed

2

T
(1 - τ)T

τT
ΡR

PT

PB

hp
ST

hs

SR

EH for
Primary System

DT for
Primary System

DT for
Secondary System

EH for
Secondary System

PB transmission in a dedicated frequency band

(a)

(b)

Fig. 1. (a) The considered system configuration. The parameters h{p,s}
denote the involved channel gains, which are explicitly defined hereinafter;
(b) The mode-of-operation for the proposed EH interweave CR scheme.

faded channels are assumed for the various involved links to
cope with realistic conditions. The reason why the κ − µ
shadowed distribution is adopted for the channel modeling
is threefold; generality, tractability and accuracy. In fact, this
model includes other popular distributions as special cases,
such as Rician, Nakagami-m and Rayleigh fading models.
Most importantly, it provides an improved fit to field measurements than the latter models and a relative mathematical
simplicity at the same time [12].
The communication activity of the single-antenna primary
system occurs in consecutive time frames, each with a fixed
duration T . The primary Tx uses a certain fraction of the
frame with duration τ T for EH and then utilizes the remaining
time (1 − τ )T for data transmission (in a similar basis as
in, e.g., [14], [15]). In out-band EH systems, PB operates
in a dedicated frequency band (different from the one used
for data communication).1 Capitalizing on this approach, the
secondary system operates in an interweave CR basis, such
that it causes no interference onto the primary system. To this
end, the secondary system uses the former time fraction for
data communication and the latter one for EH (see Fig. 1b). In
what follows, we focus on the analysis and evaluation of the
secondary system performance, which is the scope of current
work, given a fixed switching time τ specified by the primary
system.2
During the EH phase of the nth time instance, the received
signal at ST, yp , reads as
q
(1)
yp [n] = Pb d−α
PBST hp [n]sPB + nST [n],
where dPBST and α is the distance and path-loss factor of
the PB − ST channel link, respectively. Also, hp [n], sPB
and nST [n] denote, respectively, the channel fading, reference
signal and additive white Gaussian noise (AWGN) at ST,
regarding the PB − ST channel link.
In current work, a ‘one-shot’ policy is adopted; the data
transmission phase occurs in a given time instance only when
the energy buffer has been full (= M T Joules) in the previous
time instance(s), where M is the predetermined Tx power at
ST, while T M denotes the total capacity of the energy storage
1 As an illustrative example, typical frequencies that are used for prototype
development of far-field wireless power transfer systems are 2.45GHz and
5.8GHz [2].
2 The switchning time τ can be computed by the primary system without
taking into account any underlying secondary activity, e.g., as in [15].
Afterwards, τ is captured by the secondary system via feedback signaling,
which can essentially be used to inform the involved primary nodes.

circuit at ST. Then, it uses the required energy to transmit
during τ T fraction of time (i.e., using a predetermined Tx
power M to meet the appropriate transmission quality level)
in a single shot. Thus, the corresponding harvested energy at
the nth time instance is given by


η(1 − τ )T Pb |hp |2 [n]
EH[n] , min
, τT M ,
(2)
dα
PBST
where 0 < η ≤ 1 represents the power-to-energy conversion efficiency factor. Studying a finite-size energy buffer,
supercapacitors can be leveraged since they provide quite a
fast charging rate, long cycle lifetime, high storage efficiency,
while they can be integrated on chip, which makes them
suitable storage devices for EH nodes. [1], [16].
Afterwards, based on a full energy buffer at the nth time
instance, the secondary system enters the data transmission
phase in the (n+1)th time instance. Hence, the received signal
at SR, conditioned on the fact that EH[n] = M T , is expressed
as
√
M
ys [n + 1] = αs hs [n + 1]ss [n + 1] + nSR [n + 1], (3)
dSTSR
where dSTSR and αs is the distance and path-loss factor of the
ST − SR channel link, respectively. Also, hs [n + 1], ss [n + 1]
and nSR [n + 1] denote, respectively, the channel fading, unitenergy transmitted symbol and AWGN at SR, regarding the
ST − SR channel link. Otherwise, it transmits no data during
τ T [n + 1], whereas it enters to the EH phase for (at least) a
duration of T [n + 1] instead.3
d
Moreover, it is assumed that {nST , nSR } = CN (0, N0 ),
as time-invariant random processes. Also, the channel fading
gains are modeled as κ − µ shadowed RVs, such that
d

|hv |2 = S(N, µ, m),

v ∈ {p, s},

with a corresponding PDF defined as [11], [12]


N Cj xmj −1 exp − x
X
Ωj
,
f|hv |2 (x) =
mj
(m
−
1)!
Ω
j
j
j=0

(4)

(5)

while µ ≤ m are specific integer parameters, which are
explicitly defined in [11, Table I], N , m − µ, mj , m − j,
Ωj , (µK +m)/(mµ(1+K))
K denoting the Rician−K
 m with
µK
j
m−µ−j
factor and Cj , m−µ
(
)
(
. It is notej
µK+m
µK+m )
worthy that the µ parameter in (4) reflects the number of Tx
antenna elements. In particular, in the case of a single-antenna
d
PB (see Fig. 1a), µ = 1 and |hp |2 = S(N, 1, m), while in the
case of a L−antenna PB then µ = L and hence the effective
d
received channel gain becomes khp k2 = S(N, L, m).
Some other key statistic results directly arising from (5) and
are used in the rest of this paper are expressed as
 r
x


j −1 C
N m
X
X
j Ωj
x
, (6)
exp −
F|hv |2 (x) = 1 −
r!
Ωj
j=0 r=0


(7)
E |hv |2 = 1,
3 Hereinafter, for notational simplicity, we use a normalized time duration
by setting T = 1.

3


 
N
 X

1
Cj ψ(mj ) − ln
E ln |hv |2 =
.
Ωj
j=0

(8)

III. P ERFORMANCE M ETRICS
According to the proposed mode-of-operation, it is naturally
desirable to provide energy efficiency to ST, i.e., the EH phase
in each time instance should be enough to transmit in every
frame. Doing so, the total system throughput is also maximized
since ST never stays inactive at a given frame. To this end, it is
quite important to bound the effective range/distance between
the PB − ST channel link in order to satisfy this condition.
A. Effective EH Range
To derive the effective dPBST , according to (2), η(1 −
τ )Pb |hp |2 [n]d−α
PBST ≥ τ M should hold. Conditioned on the
latter constraint, the average channel capacity reads as C ,
η(1−τ )P |h |2 |h |2
E[τ log2 (1+ τ N0 dα b pdαs s )]. Further, it holds that C ≥ CL ,
PBST STSR
where CL denotes a lower bound of C. Using the fact that
log2 (1 + exp(z)) is convex with respect to z ∀z > 0, (8),
and based on Jensen’s inequality, we get
!
η(1 − τ )Pb exp E[ln(|hp |2 )]+E[ln(|hs |2 )]
CL = τ log2 1 +
αs
τ N0 dα
PBST dSTSR

h

i !
P
1
η(1 − τ )Pb exp 2 N
j=0 Cj ψ(mj )−ln Ωj
.
= τ log2 1 +
αs
τ N0 dα
PBST dSTSR
(9)
On the other hand, the actual average capacity experienced
A
at SR (using a fixed Tx power M at ST) is given by C ,
2
A
s|
E[τ log2 (1+ NM|h
)], where it holds that C ≥ CLA with CLA
αs
0 dSTSR
denoting the lower bound of the benchmark average channel
A
capacity C . Following a similar methodology as previously,
we get
h
 i !
P
N
1
ψ(m
)
−
ln
M
exp
C
j
j
j=0
Ωj
A
C = τ log2 1 +
.
s
N0 dα
STSR
(10)
Thereby, setting CL = CLA , we can directly interrelate the
lower achievable channel capacity of the secondary system
(given a fixed Tx power M ), with the effective distance of
the PB − ST channel link. Doing so, the effective range of
dPB−ST becomes dPB−ST ≤ d⋆ with
 α1



 
N
X
1
η(1
−
τ
)P
b
 .
exp 
Cj ψ(mj ) − ln
d⋆ , 
τM
Ω
j
j=0

(11)

B. Outage Probability and Average Throughput
Outage probability at the nth time instance is defined as
the probability that ST transmits yet the received SNR at SR
is below a prescribed outage threshold, γth , or the probability
that ST does not transmit and utilizes solely EH. It is expressed
as
(n)

Pout (γth ) , Ptr [n]FSNR (γth ) + (1 − Ptr [n]),

(12)

where Ptr [n] denotes the transmission probability at the nth
time instance, γth , 2R − 1 with R being the target
(normalized) data rate for the secondary service in bps/Hz and




γth N0
M |hs |2
≤ γth = F|hs |2
.
FSNR (γth ) = Pr αs
s
dSTSR N0
M d−α
STSR
(13)
The remaining issue in order to derive the outage probability
is Ptr [·]. The transmission probability at the nth time instance
is modeled as
( 

|hp |2 [n − 1]
+ (1 − τ )M ≥ M
Ptr [n] , Pr ηPb (1 − τ )
dα
PBST
) (
× Pr [dPBST ≤ d⋆ ]

×

+

Pr [dPBST > d⋆ ]

!)

|hp |2 [n − 1]
Pr ηPb
+ (1 − τ )M ≥ M + J (l)
,
dα
PBST
(14)


where
"

J (l) , Pr

Pn−1

|hp |2 [i]
ηPb i=n−1−l
dα
PBST
Pn−2
|hp |2 [i]
ηPb i=n−1−l
dα
PBST

+ (1 − τ )M ≥ M

!

+ (1 − τ )M < M

!#

&
,

(15)

and l stands for the number of previous consecutive time
instances where ST has been inactive (not transmitting) and
used solely EH instead. The first term of (14) in curly brackets
denotes the probability that ST is inside the effective PB− ST
range. Therefore, the EH activity during the remaining fraction
at the (n−1)th time instance (along with the remaining energy
(1 − τ )M after transmission) is enough to meet the prescribed
standards of the energy buffer and one-shot policy.
On the other hand, the second term of (14) in curly brackets
denotes the case when the PB−ST link distance is outside the
effective range. In this case, transmission in every consecutive
time instance is not satisfied because the received harvested
energy at ST does not meet the aforementioned condition due
to severe channel fading and/or strong propagation attenuation.
There are two possible scenarios: (i) ST does not transmit at
the previous (n− 1)th time instance, whereas it uses the entire
duration for EH, and the received energy is enough (i.e., ≥ M )
to transmit at the nth time instance; (ii) The harvested energy
was not enough for l − 1 consecutive time instances, while it
satisfies the aforementioned energy standards for l ones, up to
the (n − 1)th time instance.
The latter scenario (ii) is modeled by J (l), which includes
and
correlated
Pn−2RVs, such2 that
Pl mutually
Pn−1 two ordered
2
2
|hp | [i] =
|h
|
[i]
and
|h
|
[i]
=
p
p
i=1
Pl i=n−1−l
Pl−1
Pi=n−1−l
l−1
2
4
2
2
|h
|
[i].
Hence,
|h
|
[i]
≤
|h
|
[i],
while
p
p
p
i=1
i=1
i=1
a similar approach as in [17, Eq. (14) and Fig. 3a] can be used
to show that
J (l) =
4 It is assumed that small-scale channel fading coefficients are mutually
independent across different time instances.

4

FPl−1 |hp |2 [i]
i=1



τM
ηPb d−α
PBST



− FPl

i=1

|hp

|2 [i]




τM
.
ηPb d−α
PBST
(16)

Pl−1
d
2
=
S(N, l − 1, m) and
that
i=1 |hp | [i]
d
2
i=1 |hp | [i] = S(N, l, m). Thus, (16) can be easily
computed via (6), given a certain value of l ≥ 2.
In addition, assuming that the random placement of ST
within a given range follows the uniform (spatial) distribution,
the PDF of dPBST is given by

2x
 d2max −d2min , when dmin ≤ x ≤ dmax ,
fdPBST (x) =

0,
otherwise,
(17)

Recall
Pl

where dmin denotes a minimum allowable distance from PB
(mainly for human safety reasons due to the high radiated
power of PB [3, Remark 4]; otherwise, dmin = 0) and dmax
is the radius of coverage (circular) area of PB.
Capitalizing on the above results, the first term in curly
brackets of (14) can be rewritten as
 R ⋆


α
d

dmin ≤ d⋆ ,
 dmin F |hp |2 ηPτbMx
(1−τ ) fdPBST (x)dx,



0,

otherwise.
(18)

Similarly, the second term in curly brackets reads as
 R
h
i


dmax
⋆
M xα

+J (l) fdPBST (x)dx, d ≤ dmax ,
F |hp |2 τ ηP
 d⋆
b


≈

0,

otherwise
(19)

 R dmax


⋆
xα
fdPBST (x)dx, d ≤ dmax ,
 d⋆ F |hp |2 τ M
ηPb


0,

(20)

otherwise,

where (20) follows by neglecting J (l) from (19). For a
numerical setting of practical interest (provided in the next
section), J (l) takes values in the order of 10−6 by setting
l = 2, while it takes values ≪ 10−8 for l = 3. On the
other hand, evaluating (20) yields values ≫ 10−3 . Hence, the
approximation of (20) is very tight.
Using (6) and (17), (18) and (20) are, respectively, evaluated
with the aid of [13, Eq. (3.381.9)] as
Φ1 ,

− 2

α
τM

PN Pmj −1 2Cj ηPb (1−τ
)Ωj



2
r=0
⋆
2
j=0

r!α((d ) −dmin )



i
h 

α

τ M(d⋆ )α
2

min
 × Γ r + α2 , ηPτ Md
−
Γ
r
+
,
,
α ηPb (1−τ )Ωj
b (1−τ )Ωj











dmin ≤ d⋆ ,

0,

and
Φ2 ,

otherwise,
(21)



− 2

PN Pmj −1 2Cj ηPτbMΩj α





−(d⋆ )2 )
i
hj=0 r=0 r!α(d⋆2max

α

τ Mdα
2 τ M(d )
max

− Γ r + α2 , ηPb Ω
,
 × Γ r + α , ηPb Ωj
j











d⋆ ≤ dmax ,

0,

otherwise.
(22)

Finally, by inserting (21) and (22) into (14), we get Ptr [n] ≈
Φ1 + Φ2 . Also, inserting (13) into (12), the outage probability
is derived in a closed formula.
Furthermore, the average effective throughput (in bps/Hz)
is defined as


(23)
T , τ R 1 − Pout 2R − 1 ,
which can be directly computed via (12).

C. Practical Energy Harvesting Imperfections
In practice, the EH and power amplifier circuits may not
perform ideally. In this case, for a given Tx power at ST, i.e.,
M , the corresponding power amplifier consumes a total power
ρM , where ρ ≥ 1 denotes the power amplifier inefficiency
factor. In addition, the EH circuit consumes a constant power
Pc during transmission [10]. Hence, according to the proposed
mode-of-operation, ST transmits using a constant power M
whenever its stored energy reaches M , ρM + Pc ; otherwise
it does not transmit and performs EH instead. Thus, in the
case of EH imperfections, the transmission probability in (14)
becomes
( 

|hp |2 [n − 1]
+ (1 − τ )M ≥ M
Ptr [n] ≈ Pr ηPb (1 − τ )
dα
PBST
) (
× Pr [dPBST ≤ d⋆ ]

+

Pr [dPBST > d⋆ ]

|hp |2 [n − 1]
× Pr ηPb
+ (1 − τ )M ≥ M
dα
PBST


)

.

(24)

To this end, after some simple manipulations, we can easily
show that the previously presented analytical results for the
considered performance metrics are valid in this case, by
substituting M with M in (21) and (22).
IV. N UMERICAL R ESULTS AND C ONCLUDING R EMARKS
In this section, numerical results are presented and compared with Monte-Carlo simulation results, which are illustrated as line-curves and circle-marks, respectively. In the
ensuing simulation experiments, the Tx power at PB is set
to be Pb = 33dBm, the Tx power at ST is M = 20dBm,
{α = 2.4, m = 20, K = 7}, which reflect line-of-sight
propagation conditions for the PB−ST channel link [11], [12],
the power-to-energy conversion efficiency is η = 0.85 [3], and
N0 = −101dBm. For the ST − SR signal path loss, as = 3 is
assumed. Also, for the scenario of EH imperfection, ρ = 1.2
and Pc = −30dBm (i.e., M ≈ 0.12mJ for T = 1ms). Finally,
{dmin = 1m and dmax = 15m}, which correspond to the
coverage of the PB − ST channel link, while dSTSR = 30m.

5

2

L = 1 (solid line)
L = 16 (dashed line)

Average Throughput (bps/Hz)

Outage Probability

0.2

Non-ideal

0.15

0.1

Ideal

0.05

0
0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

Switching Time τ

Ideal
(solid line)
Non-ideal (dashed line)
1.5

R=4
1

R=2
0.5

0
0.1

R=1
0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

Switching Time τ

Fig. 2. Outage probability as per (12) vs. various values of the switching
time τ for a single-antenna PB (i.e. L = 1) and for a multi-antenna PB (i.e.,
L = 16), when R = 1bps/Hz.

Fig. 3. Average effective throughput as per (23) vs. various values of the
switching time τ for a single-antenna PB and for three different scenarios of
target data rate R (in bps/Hz).

In Fig. 2, the outage performance of the considered secondary system is presented for two different setups; the ideal
EH case (where ρ = 1 and Pc = 0) and the non-ideal case
reflecting the practical imperfection of EH circuit. Obviously,
the outage performance is worse for a higher switching time τ .
This occurs since ST transmits a longer duration in each time
instance for higher τ , which corresponds to more energy losses
of the energy buffer at ST. In turn, the probability that ST
stays inactive (i.e., non-transmitting) in a given time instance
is increased, which also results to an increased probability
of an outage event. It is also noteworthy that the case when
a multi-antenna PB is adopted drastically impacts on the
enhancement of the outage performance since the amount of
collected energy at ST is proportionally enhanced; hence, the
probability that ST transmits using a constant Tx power M ,
given a full energy buffer, is correspondingly increased.
In Fig. 3, the average effective throughput of the considered
system is presented, which is enhanced for higher τ values.
As expected, this is a reasonable outcome since for a longer
transmit duration, in each time instance, more data streams
can potentially be transferred from ST to SR. As noticed
from Fig. 2 and also verified in Fig. 3, the efficiency of EH
circuit (i.e., the quality of its hardware gear) plays a key role
to the overall performance of the considered system. To this
end, network practitioners working into this field should pay
attention to the efficient design of EH equipment to maintain
sufficient communication quality in EH-enabled transmission
systems.

[5] C. Bockelmann, N. Pratas, H. Nikopour, K. Au, T. Svensson, C. Stefanovic, P. Popovski, and A. Dekorsy, “Massive machine-type communications in 5G: Physical and mac-layer solutions,” IEEE Commun. Mag.,
vol. 54, no. 9, pp. 59–65, Sep. 2016.
[6] X. Huang, T. Han, and N. Ansari, “On green-energy-powered cognitive
radio networks,” IEEE Commun. Surveys Tuts., vol. 17, no. 2, pp. 827–
842, Second Quarter 2015.
[7] N. I. Miridakis, T. A. Tsiftsis, G. C. Alexandropoulos, and M. Debbah,
“Green cognitive relaying: Opportunistically switching between data
transmission and energy harvesting,” IEEE J. Sel. Areas Commun.,
vol. 34, no. 12, pp. 3725–3738, Dec. 2016.
[8] M. L. Ku, W. Li, Y. Chen, and K. J. R. Liu, “Advances in energy
harvesting communications: Past, present, and future challenges,” IEEE
Commun. Surveys Tuts., vol. 18, no. 2, pp. 1384–1412, Second Quarter
2016.
[9] S. Lee and R. Zhang, “Distributed wireless power transfer with energy
feedback,” IEEE Trans. Signal Process., vol. 65, no. 7, pp. 1685–1699,
Apr. 2017.
[10] R. Morsi, D. Michalopoulos, and R. Schober, “Performance analysis
of near-optimal energy buffer aided wireless powered communication,”
IEEE Trans. Wireless Commun., To appear, 2018.
[11] F. J. Lopez-Martinez, J. F. Paris, and J. M. Romero-Jerez, “The κ-µ
shadowed fading model with integer fading parameters,” IEEE Trans.
Veh. Technol., vol. 66, no. 9, pp. 7653–7662, Sep. 2017.
[12] U. Fernández-Plazaola, L. Moreno-Pozas, F. J. López-Martı́nez, J. F.
Paris, E. Martos-Naya, and J. M. Romero-Jerez, “A tractable product
channel model for line-of-sight scenarios,” Under review, 2017.
[Online]. Available: http://arxiv.org/abs/1711.05650
[13] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and
Products. Academic Press, 2007.
[14] C. Zhang and Y. Chen, “Wireless power transfer strategies for cooperative relay system to maximize information throughput,” IEEE Access,
vol. 5, pp. 2573–2582, 2017.
[15] F. Zhao, H. Lin, C. Zhong, Z. Hadzi-Velkov, G. K. Karagiannidis, and
Z. Zhang, “On the capacity of wireless powered communication systems
over Rician fading channels,” IEEE Trans. Commun., vol. 66, no. 1, pp.
404–417, Jan. 2018.
[16] K. Grigoras, J. Keskinen, L. Grnberg, E. Yli-Rantala, S. Laakso,
H. Vlimki, P. Kauranen, J. Ahopelto, and M. Prunnila, “Conformal
titanium nitride in a porous silicon matrix: A nanomaterial for in-chip
supercapacitors,” Nano Energy, vol. 26, pp. 340 – 345, 2016.
[17] H.-C. Yang, “New results on ordered statistics and analysis of minimumselection generalized selection combining (GSC),” IEEE Trans. Wireless
Commun., vol. 5, no. 7, pp. 1876–1885, Jul. 2006.

R EFERENCES
[1] A. Kansal, J. Hsu, S. Zahedi, and M. B. Srivastava, “Power management
in energy harvesting sensor networks,” ACM Trans. Embed. Comput.
Syst., vol. 6, no. 4, Sep. 2007.
[2] E. Falkenstein, D. Costinett, R. Zane, and Z. Popovic, “Far-field RFpowered variable duty cycle wireless sensor platform,” IEEE Trans.
Circuits and Systems II: Express Briefs, vol. 58, no. 12, pp. 822–826,
Dec. 2011.
[3] M. Xia and S. Aı̈ssa, “On the efficiency of far-field wireless power
transfer,” IEEE Trans. Signal Process., vol. 63, no. 11, pp. 2835–2847,
Jun. 2015.
[4] S. Haykin, “Cognitive radio: brain-empowered wireless communications,” IEEE J. Sel. Areas Commun., vol. 23, no. 2, pp. 201–220, Feb.
2005.

