1

Outage Analysis of Ambient Backscatter
Communication Systems

arXiv:1711.09214v1 [] 25 Nov 2017

Wenjing Zhao, Gongpu Wang, Saman Atapattu, Chintha Tellambura, and Hao Guan

Abstract—This paper addresses the problem of outage characterization of an ambient backscatter communication system
with a pair of passive tag and reader. In particular, an exact
expression for the effective channel distribution is derived. Then,
the outage probability at the reader is analyzed rigorously.
Since the expression contains an infinite sum, a tight truncation
error bound has been derived to facilitate precise numerical
evaluations. Furthermore, an asymptotic expression is provided
for high signal-to-noise ratio (SNR) regime.
Index Terms—Ambient backscatter, Internet of Things (IoT),
outage probability, performance analysis.

I. I NTRODUCTION
Backscatter communication systems, such as Radio Frequency Identification (RFID) system, enable connecting massive small computing devices specially for applications in
Internet of Things (IoT) [1]. The traditional RFID system
typically consists of a tag and a reader. The reader first
generates and transmits an electromagnetic wave signal to the
tag, and then the tag receives and backscatters the signal to
the reader.
One disadvantage for the RFID system is that the reader
needs an oscillator to transmit a carrier signal, for which
dedicated encoding/decoding circuitry and power supply are
required [2]. While these are essential components for a
successful communication, such in-built technology may no
longer be promising for small-scale devices. To overcome
such overheads, ambient backscatter prototypes are proposed
in [3, 4].
The ambient backscatter technology utilizes environmental
wireless signals (e.g., digital TV broadcasting, cellular or
Wi-Fi) for both energy harvesting and information transmission, which avoids battery as well as manual maintenance.
Specifically, the tag indicates bit 1 or bit 0 through reflecting
or non-reflecting state, and the reader decodes the received
backscattered signal accordingly [5]. Ambient backscatter may
be widely used for future applications (e.g., many applications
in IoT with sensors located in dangerous spots filled with
poisonous gases/liquids, or inside building walls) that are
inconvenient and unsafe for wired communications [6].
W. Zhao and G. Wang are with School of Computer and Information
Technology, Beijing Jiaotong University, Beijing 100044, China (e-mail:
wenjing bjtu@163.com; gpwang@bjtu.edu.cn).
S. Atapattu are with the Department of Electrical and Electronic Engineering, The University of Melbourne, Melbourne, VIC 3010, Australia (e-mail:
saman.atapattu@unimelb.edu.au).
C. Tellambura is with the Department of Electrical and Computer Engineering, University of Alberta, Edmonton, AB T6G 2V4, Canada (e-mail:
chintha@ ece.ualberta.ca).
H. Guan is with Nokia Bell Labs, Beijing, China (e-mail:hao.guan@nokiasbell.com).

6

smartphone

laptop

RF Source

d st hst

hsr d sr
R

Reader

d tr
htr

T
Tag

Fig. 1. An ambient backscatter communication system.

The performance analysis of the ambient backscatter communication is considered over real Gaussian channels in [7],
and complex Gaussian channels in [8, 9]. The bit error rate
(BER) is derived and the BER-based outage probability is
obtained in [7] for ambient bakcscatter communication systems with energy detector. In addition, the outage capacity
optimization problem is investigated in [8] when successive
interference cancellation (SIC) method is applied. Besides, the
BER-based outage probability of a semi-coherent detection
scheme is calculated in [9] in the case of perfect and imperfect
channel state information (CSI), respectively.
To our best knowledge, effective channel distribution for
ambient backscatter communication systems has not be addressed and the outage performance based on signal-to-noise
ratio (SNR) remains an open problem, which is the focus of
this paper.
In this paper, we consider an ambient backscatter communication system over real Gaussian channels. We derive an
exact expression for the effective channel distribution in this
system. Particularly, we evaluate the outage performance and
analyse its asymptotic outage performance at high transmit
SNR. Moreover, since the derived outage probability is the
summation of infinite items, the corresponding truncation error
bound is calculated.
II. S YSTEM M ODEL
We consider an ambient backscatter communication system
comprised of an ambient RF source (S) and a pair of passive
tag (T ) and reader (R) (Fig. 1). While the RF source communicates with its legacy users (e.g., smartphones, laptops,
etc.), both tag and reader may also receive that source signal.
The tag first harvests energy from the source signal, and

2

then communicates with the reader via ambient backscatter.
Particularly, the tag can backscatter or consume the energy of
the received signal to represent two states “1” or “0” for the
reader, respectively [5].
The fading channels of S − R, S − T , and T − R links are
denoted as h̃sr , h̃st and h̃tr , respectively, which are real Gaus2
sian random variables (RVs) distributed as h̃sr ∼ N (0, σ̃sr
),
2
2
2
2
2
h̃st ∼ N (0, σ̃st
) and h̃tr ∼ N (0, σ̃tr
), where σ̃sr
, σ̃st
and σ̃tr
are channel variances. Further, the corresponding distances are
dsr , dst and dtr , respectively. Without loss of generality, we
consider time instance n. The signal received at the tag can
be given as
h̃st
(1)
yt (n) = p α s(n),
dst
where s(n) is the source signal with the average power P ,
and α is the path-loss exponent. The signal backscattered by
the tag can be written as
x(n) = ηB(n)yt (n),

(2)

where B(n) ∈ {0, 1} is a binary symbol and η ∈ [0, 1] is the
attenuation factor. Then, the received signal at the reader can
be given as [5]
h̃tr
h̃sr
yr (n) = p s(n) + p α x(n) + w(n)
dα
dtr
sr
= hs(n) + w(n),

(3)

where w(n) is the additive white Gaussian noise (AWGN)
2
at the reader with zero mean and σw
variance, and h is the
effective channel gain which can be given for two states as

hsr ,
B(n) = 0,
h=
(4)
hsr + ηhst htr , B(n) = 1,
2
2
),
), hst , √h̃stα ∼ N (0, σst
where hsr , √h̃srα ∼ N (0, σsr
dsr

and htr
2
σtr
=

2
2
) with σsr
=
, √h̃trα ∼ N (0, σtr

2
σ̃tr
.
dα
tr

dtr

dst
2
σ̃sr
2
, σst
dα
sr

=

2
σ̃st
dα
st

and

III. P ERFORMANCE A NALYSIS
A. Effective Channel and SNR Distributions
We first derive the probability density function (PDF) of the
effective channel, fh (x), which can be given as

2
− x2


√ 1 2 e 2σsr
,
B(n) = 0,


2πσsr


2


σsr
∞
x2 X 2k Γ2 k + 1 x2k
2 σ2
4η2 σst
− 2σ
tr
2
2
e√
fh (x) =
(5)
sr
e

2
2k
2π 3 σsr

(2k)!σsr

k=0 



2

σsr
 ×W−k,0
,
B(n) = 1,
2
2
2η 2 σst σtr

where Wa,b (·) and Γ(·) are the Whittaker function [10,
eq. (9.223)] and the Gamma function [10, eq. (8.310.1)],
respectively. The proof is in Appendix A.
2
The receive SNR at the reader is ρ = Pσh2 = ρ̄h2 where
w
ρ̄ = σP2 may be the average transmit SNR. With a variable
w
transformation for (5), we can derive the PDF of ρ, fρ (x), as
 q i
h q 
q 
x
x
x
−
fh
+
f
f
h
h
ρ̄
ρ̄
ρ̄
(6)
fρ (x) =
√
= √
,
2 ρ̄x
ρ̄x

where the second equality follows as the PDF fh (x) is an even
function, i.e., fh (x) = fh (−x).
B. Outage Probability
The outage probability is the probability that the SNR at the
reader falls below a certain predetermined threshold ρt . Thus,
it can be derived as
Z ρt
fρ (x)dx
Po = Pr[ρ ≤ ρt ] =
0
q

ρt


,
B(n) = 0,
erf
2

2ρ̄σsr





∞
 X
2
22k Γ2 (k + 12 )
σsr
√
W−k,0
(7)
=
2 σ2
2η 2 σst

π 3 (2k)!
tr

k=0


σ2




2 2
 ×e 4η2 σsr
st σtr γ
,
B(n) = 1,
k + 1 , ρt
2
2 2ρ̄σsr

where erf(·) and γ(·, ·) are the error function [10, eq. (8.250.1)]
and the incomplete gamma function [10, eq. (8.350.1)], respectively. The equation (7) can be obtained by following from [10,
eq. (3.321.2)], [10, eq. (8.250.1)] and [10, eq. (3.381.1)].

C. Truncation Error Bound
Since the outage probability expression for B(n) = 1 case
in (7) is with an infinite sum, it is a challenge for numerical
calculation. We therefore truncate it into a finite number
of terms T in order to ensure a given numerical accuracy
requirement. Then, we bound the truncation error as
 s


Ψ 12 , 0, ν
2ρt
ρt
|ǫT | ≤ √
γ
T
+
1,
2 ρ̄
2 ρ̄
πνT !
σsr
2σsr
(8)


3 ρt
− 2γ T + , 2
,
2 2σsr ρ̄
σ2

where ν = 2η2 σsr
2
2 , and Ψ(·, ·, ·) is the confluent hyst σtr
pergeometric function [10, eq. (9.211.4)]. The proof is in
Appendix B.

D. Asymptotic Analysis for High SNR
To further investigate the ambient backscatter system, we
approximate outage of the reader for large SNR as
 q 2ρ
1
t √

,
B(n) = 0,
2

πσsr
ρ̄

2
σsr


2
Po ≈
(9)
√
2 2 2
σsr
2ρt e 4η σst σtr W0,0

2 2 2


√ 2 2η σst σtr √1ρ̄ , B(n) = 1.
πσsr

The proof is in Appendix C.
Interestingly, when transmit SNR tends to infinity, we get
log Pout
1
that the diversity gain is − lim log
SNR = 2 for both cases
SNR→∞
B(n) = 0 and B(n) = 1. Accordingly, when power is
big enough, the outage probability of the reader is inversely
proportional to square root of the power.

3

100

10-1

101

Relative Error

10-2

100

10-3

B(n)=0
B(n)=0

10-4

B(n)=1

B(n)=1

10-1

-5

10

10-6

0

1

2

3

4

5

6

7

8

9

10

-20

-15

-10

-5

0

5

10

15

20

25

30

Number of Terms T

Fig. 2. The relative error versus the number of terms T for ρt = −3dB and
ρt = 7dB when ρ̄ = 3dB and B(n) = 1.

IV. N UMERICAL

AND

S IMULATION R ESULTS

This section provides simulation results based on the system
model in Section II and numerical results based on the
analytical results in Section III. Here, the attenuation factor η
2
2
2
is set as 0.7. Besides, we assume σsr
= σst
= 1 and σtr
= 3,
unless otherwise specified.
For the state B(n) = 1, we derive a truncation error bound
(8). Fig. 2 shows the relative error versus the number of
terms T when ρ̄ = 3dB for ρt = −3dB and ρt = 7dB. We
; the
calculate the relative error with truncation as |exact−truncated|
exact
bound in (8)
relative error with bound as errorexact
;
the
exact
value
with
value
numerical integration; and the truncated value for different T
by using (7). The relative error with truncation is less than
10−5 when T > 2 and T > 5 for ρt = −3dB and ρt = 7dB,
respectively. The relative error with bound is less than 10−5
when T > 3 and T > 8 for ρt = −3dB and ρt = 7dB,
respectively. This shows the tightness of the bound. Moreover,
by observation, we may say that a very accurate value can be
calculated using small T .
Fig. 3 illustrates the outage probability Po versus the
average transmit SNR ρ̄ when ρt = 2dB and ρt = 15dB.
The outage probability decreases with the increasing average
transmit SNR. The asymptotic expressions (9) also approach
the exact values asymptotically at high SNR. Since the diversity gain is 1/2, the slope of the asymptotic outage curves is
1/2.
Fig. 4 depicts the outage probabilities Po versus the distance
dtr between the tag and the reader. We consider both cases
of SNR ρ̄ = 10dB and ρ̄ = 20dB, respectively. We set ρt =
2 dB, α = 3, dsr = 20 meters and dst = 20 meters. Besides,
2
2
2
the channel variances σ̃sr
, σ̃st
and σ̃tr
are set as 1, 1 and 3,
respectively. The outage probability Po is a constant in the
case of B(n) = 0 due to no transmission between the tag and
the reader. However, when enlarging the distance dtr between
the tag and the reader, the outage probability Po witnesses a
upward trend in the case of B(n) = 1. For example, if we
expect Po < 10−4 , the distance dtr between the tag and the

Fig. 3. The outage probability Po versus the average transmit SNR ρ̄.

10-2

B(n)=0
10-3

10-4
B(n)=1

0

2

4

6

8

10

12

14

16

18

20

Fig. 4. The outage probabilities Po versus the distance dtr between the tag
and the reader.

reader should not exceed 3 meters when ρ̄ = 10dB and 7
meters when ρ̄ = 20dB.

V. C ONCLUSION

Ambient backscatter, a new form of wireless communication, has potential commercial value as well as a series of
open problems. In this paper, we first derived the effective
channel distribution for the ambient backscatter communication system. We next analyzed the outage probabilities,
its truncation error bound as well as the asymptotic outage
probabilities at high SNR. It was found that the asymptotic
outage probabilities could well approach the exact values,
and our truncation error bound could provide a reasonable
estimation of the truncation terms.

4

A PPENDIX

B. Proof of (8)
On the basis of (7), the truncation error ǫ(T ) with the
number of terms T can be bounded as

A. Proof of (5)

ǫ(T )

The distribution of any real Gaussian channel hab ∈
{hsr , hst , hrt } is
2

− x2
1
fhab (x) = p
e 2σab ,
2
2πσab

(10)

2
2
2
2
where σab
∈ {σsr
, σst
, σrt
}.

In the case of B(n) = 1, we have h = hsr + ηhst htr . The
distribution fξ (x) of ξ = ηhst htr can be shown as [11]

1
p
fξ (x) =
K0
2 σ2
πη σst
tr

|x|
p
2 σ2
η σst
tr

!

,

(11)

where Kα (·) is the modified Bessel function of the second
kind [10].
Since the two random variables hsr and ξ are independent,
the distribution fh (x) is the convolution of fhsr (x) and fξ (x).
Therefore, we can obtain

fh (x) =

Z

∞

−∞

(a)

=

e

fhsr (x − z)fξ (z) dz
2

x
− 2σ
2

sr

δϕ

=

2e

x
− 2σ
2

sr

δϕ
2

(c)

=

2e

∞

0

2

(b)

Z

x
− 2σ
2

sr

δϕ


z − z22  σxz
− xz
2
2
e sr
dz
+ e σsr
K0 ( )e 2σsr
ϕ

∞

xz
z − z22
cosh( 2 )dz
K0 ( )e 2σsr
ϕ
σsr
0
Z ∞
∞
X
z − 2σz22 2k
x2k
)e sr z dz
K
(
0
4k (2k)!
σsr
ϕ
0
Z

(12)

k=0

where

p
2 ,
2π 3 σsr

√ Z ∞ −νx 1
Z ρ2t
2σsr ρ̄ γ(T + 1, y)
ν
e
x2
<
dx
dy
√
πT ! 0 (1 + x) 12
y
0
s



Ψ( 12 , 0, ν)
ρt
2ρt
√
γ T + 1, 2
=
2 ρ̄
πνT !
σsr
2σsr ρ̄


3 ρt
− 2γ T + , 2
,
2 2σsr ρ̄

(d)

where ν is defined in (13), (a) is obtained from [10,
eq. (9.222.1)] and [10, eq. (8.350.1)], (b) follows by setting
j = k − T − 1 and leveraging the integral representation of
Hypergeometric function [10, eq. (9.211.4)], (c) utilizes the
following equation


−(T +1)

xy
xy
xy
= (T + 1)
e x+1
1 F1 1; T + 2;
x+1
x+1


xy
×γ T + 1,
,
x+1

x
(d) is based on x+1
< 1 for x > 0, and Ψ(·, ·, ·) is the
confluent hypergeometric function [10, eq. (9.211.4)].

k=0

ν
∞
e 2 X 2k Γ2 (k + 21 )W−k,0 (ν) − 2σx22 2k
e sr x ,
=
2k
δ
(2k)!σsr

(d)

δ=




ν
∞
X
22k Γ2 k + 12 W−k,0 (ν)
e2
1 ρt
=√
γ k+ , 2
(2k)!
2 2σsr ρ̄
ππ
k=T +1
ρ
√
Z 2t
Z ∞ −νx k+ 1
∞
2σsr ρ̄
1
e
x 2
(a) ν X 1
dx
=
e−y y k− 2 dy
1
k+
π
k! 0 (1 + x) 2
0
k=T +1
√
Z ∞ Z ρ2t
3
−νx T + 2
2σsr ρ̄ e
ν
x
(b)
−y T + 12
=
y
3 e
T
+
π(T + 1)! 0
2
(1 + x)
0


xy
× 1 F1 1; T + 2;
dydx
x+1
√ Z ∞ Z ρ2t
−νx 12
2σsr ρ̄ e
x −y − 1
ν
(c)
y 2
=
1 e
πT ! 0
(1 + x) 2
0


xy
xy
x+1
×e
dydx
γ T + 1,
x+1

ϕ=η

q
2 σ2 ,
σst
tr

ν=

2
σsr
2 σ2 ,
2
2η σst
tr

(13)

(a) follows by two integrals having identical bounds, (b)
follows by using the definition of hyperbolic cosine cosh x =
1 x
−x
by replacing cosh x with its series ex2 (e +e P), (c) follows
∞
x2k
pression k=0 (2k)!
, and (d) follows from [10, eq. (6.631.3)].

C. Proof of (9)
In the case of B(n) = 0 and using series representation of
the error function erf(·) [10, eq. (8.253.1)], we can expand (7)
as

r
∞
2 X
ρt
(−1)k ρt k+1/2
√
=
erf
.
2
2 )k+1/2
2ρ̄σsr
π
k!(2k + 1)(2ρ̄σsr
k=0

For ρ̄ → ∞, we consider the lowest exponent for 1/ρ̄, i.e.,
the index k = 0. Similarly, in the case of B(n) = 1, we can
expand γ (·, ·) in (7) as
γ (a, x) =

∞
X
(−1)j xa+j
j=0

j!(a + j)

.

We then consider the lowest exponent for 1/ρ̄, i.e., the indies
k = 0 and j = 0. Thereby, when average transmit SNR tends

5

to infinity, namely, ρ̄ → ∞, the asymptotic outage probability
can be simplified as (9).
R EFERENCES
[1] A. Al-Fuqaha, M. Guizani, M. Mohammadi, M. Aledhari, and
M. Ayyash, “Internet of things: A survey on enabling technologies,
protocols, and applications,” IEEE Commun. Surveys Tuts., vol. 17,
pp. 2347–2376, Jun. 2015.
[2] T. H. Lee, The design of CMOS Radio-Frequency integrated circuits.
Cambridge Univ. Press, 1998.
[3] V. Liu, A. Parks, V. Talla, S. Gollakota, D. Wetherall, and J. R. Smith,
“Ambient backscatter: wireless communication out of thin air,” in Proc.
ACM SIGCOMM, Aug. 2013, pp. 39–50.
[4] A. N. Parks, A. Liu, S. Gollakota, and J. R. Smith, “Turbocharging
ambient backscatter communication,” in Proc. ACM SIGCOMM, Aug.
2014, pp. 619–630.
[5] G. Wang, F. Gao, R. Fan, and C. Tellambura, “Ambient backscatter
communication systems: Detection and performance analysis,” IEEE
Trans. Commun., vol. 64, pp. 4836–4846, Nov. 2016.
[6] D. Kuester and Z. Popovic, “How good is your tag?: RFID backscatter
metrics and measurements,” IEEE Microwave Mag., vol. 14, pp. 47–55,
Jul. 2013.
[7] Y. Zhang, J. Qian, F. Gao, and G. Wang, “Ambient backscatter system
with real source,” in IEEE Workshop on Signal Processing Advances in
Wireless Commun., 2017, to be published.
[8] S. Xing, R. Long, Y. Liang, M. Ding, and Z. Lin, “Outage capacity analysis for ambient backscatter communication systems,” in IEEE Wireless
Commun. and Networking Conf. (WCNC), 2017, to be published.
[9] J. Qian, F. Gao, G. Wang, S. Jin, and H. Zhu, “Semi-coherent detection
and performance analysis for ambient backscatter system,” IEEE Trans.
Commun., 2017, to be published.
[10] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and
Products. Academic Press Inc, 7th revised edition ed., 2007.
[11] M. K. Simon, Probability distributions involving Gaussian Random
Variables: a handbook for engineers and scientists. Springer, 2006.

