1. Introduction
Energy harvesting (EH) is envisioned to be a promising approach for prolonging the lifetime of energy-constrained networks, typically the wireless sensor networks (WSNs) [
1,
2]. Traditionally, the energy is harvested from natural sources, e.g., solar, wind, heat, etc. [
3]. However, this energy harvesting approach relies on some uncontrollable factors such as weather conditions, thus is not suitable to wireless communication networks that require high stability in terms of the quality-of-service (QoS). Hence, to harvest energy from the radio-frequency (RF) signals, which is capable of providing controllable energy supplies, has gained ever-increasing attention among the wireless communities [
4]. An important application scenario of this method lies in wireless power transmission (WPT), which is enabled by the radio-frequency (RF) and is capable of providing convenient and continuous power supply for wireless powered devices [
5]. Specifically, classic
time-switching and
power-splitting receiving architectures to realize practical WPT are proposed and discussed extensively in the literature [
6,
7,
8]. To the state of the art, WPT technique has been investigated in numerous systems, e.g., cognitive radio networks [
9], multiple-input–multiple-output (MIMO) networks [
10], non-orthogonal multiple access (NOMA) networks [
11], etc.
While the above works all assume the perfect channel state information (CSI), some recent studies have focused on the effect of imperfect CSI on the system’s energy and information transmission. In [
12], the secrecy performance of a multiple-input–single-output (MISO) WPT system is studied, where the CSI used for transmit antenna selection (TAS) is outdated. In [
13], the power optimization to maximize the total capacity of small cell in wireless powered heterogeneous networks is investigated, where the imperfect CSI is handled by using a non-cooperative game approach. Meanwhile, Bi and Chen et al. [
14] and Liu et al. [
15] introduced the Markov chain theory to formulate the dynamic behaviors of energy storage at the wireless powered devices without taking into account the imperfect CSI.
In this paper, we analyze the impact of imperfect CSI on the energy and information transmissions in the wireless powered sensor networks using the Markov chain theory. The PB supplies energy to the wireless powered sensor using MRT protocol if the energy of the sensor is not enough to conduct the transmission operation. Otherwise, the sensor transmits information and the destination receives it using MRC protocol. The differences between this work and the work in [
16] are obvious. First, Tang et al. [
16] considered the transmission between multiple users and a single-antenna destination, while this work studies the transmission between a single user and a multi-antenna destination. Tang et al. [
16] adopted the user selection, while this work does not. Besides, this work adopts MRC at the multi-antenna destination while Tang et al. [
16] did not. Moreover, this work considers the effect of imperfect CSI, while Tang et al. [
16] only assumed that all the channels were perfect, which becomes a highlight of this work and makes this work much more practical.
The contributions of this paper are summarized as follows:
To the state of the art, the proposed MRT/MRC WPT system with imperfect CSI is firstly investigated in this paper.
The network transmit probability, and effective and overall ergodic capacities are derived under the condition of imperfect CSI during both MRT and MRC operations to assess the impact of imperfect CSI on the energy and information transmissions.
Our results demonstrate the detrimental effect of imperfect CSI on the network transmit probability, effective and overall ergodic capacities, and the validity to improve the overall ergodic capacity by optimizing the sensor transmit power.
2. System Model
We consider a wireless powered sensor network, as shown in
Figure 1, which consists of a PB node
B, a wireless powered sensor
S, and a destination node
D. It is assumed that
B and
D are equipped with
and
antennas, and
S is equipped with a single antenna. The sensor
S is equipped with an energy storage with a finite-capacity of
. We assume that all the channels experience quasi-static Rayleigh fading so that the channel coefficients keep constant during a block time
but change independently from one packet time to another [
7,
17]. Moreover, a standard path-loss model [
7] is adopted, namely the average channel power gain
, where
denotes the expectation operation,
is the path-loss factor, and
and
denote the channel coefficient and the distance between
a and
b, respectively.
To quantify the energy storage at the sensor, we adopt a discrete-level model similar to in [
14,
15]. Specifically, the energy capacity of the storage is discretized into
L units. As such, there will be
possible energy levels in total at
S with the
lth energy level defined as
where
is the single energy unit.
Obviously, an energy amount
is required to supply a transmission operation (note that the extra power consumption of the transmitting/receiving circuitry is neglected in this paper [
7,
8]), where
represents the sensor transmit power. We set
in this paper due to the energy discretization [
14], where
denotes the transmit energy level corresponding to the energy consumed at
S of each transmission with
. We highlight that the sensor energy will always transfer within
given in Equation (
1), and the transitions among energy levels form a Markov chain [
14]. For the notation convenience, we denote the stationary probability vector of the Markov chain as
, where
is the stationary probability of the
lth energy level
.
2.1. Information Transmitting for
When the index of the energy level at
S satisfies
, the sensor energy will be sufficient so that information transmission can occur at
S. As a result, the received signal-to-noise ratio (SNR) at
D could be given by
where
, and
is the variance of the additive white Gaussian noise (AWGN).
represents the channel coefficient vector between
S and
D with
,
is the normalized MRC weight vector applied at
D satisfying
due to imperfect CSI [
18], and
is the estimated channel coefficient vector of
that can be modeled as [
12,
19]
where
denotes the channel correlation factor between the actual channel coefficient vector
and its estimation
, and
is the Gaussian random estimating error vector with each element having the variance of
.
2.2. Information Transmitting for
When the index of the energy level at
S satisfies
, the energy harvesting will occur, and the harvested energy at
S would be
where
is the transmit power of
B and
denotes the energy conversion efficiency.
represents the channel coefficient vector between
B and
S with
,
is the normalized MRT weight vector applied at
B satisfying
due to imperfect CSI, and
is the estimated channel coefficient vector of
that can be modeled as [
12,
19]
where
denotes the channel correlation factor between the actual channel coefficient vector
and its estimation
, and
is the Gaussian random estimating error vector with each element having the variance of
. Then, the harvested energy that can be saved in the storage of
S after energy discretization is derived as [
14,
15]
3. Ergodic Capacity
The ergodic capacity is defined as the expected value of the instantaneous mutual information of the received SNR [
20]. By using the total probability theorem, the overall ergodic capacity of the proposed network can be calculated as
where
is the transmit probability of the network, and
denotes the effective ergodic capacity of the network on condition that the information transmission occurs, which is defined as [
19,
20]
Next, we elaborate on the derivation of and , and then focus on the optimization of to maximize .
3.1. Derivation of Effective Ergodic Capacity
Referring to ([
18], Equation (48)), the cumulative distribution function (CDF) of
is given by
where
, and
represents the CDF of random variable
. Differentiating Equation (
9) in the cases of
and
, we have
Substituting Equation (
10) into Equation (
8), we derive
as
Resorting to ([
21], eq. (3.352.4)) and ([
21], eq. (3.353.5)), we derive
and
as
Substituting Equations (
12) and (
13) into Equation (
11), we easily derive
.
3.2. Derivation of Transmit Probability
From the calculation expression of
given after Equation (
7), we need to derive the stationary probability vector of the Markov chain
. According to the Markov chain theory, to derive
, we need to figure out the transition probabilities among all the energy levels first. Without loss of generality, we examine the transition probability from
to
within one transition,
.
3.2.1. Transition for
As described in
Section 2, the sensor will transmit information when
so that the energy in its storage will decrease
. As a result, the transition probability from
to
will be
only when
, where
.
3.2.2. Transition for
For
, the sensor will be not able to transmit information and has to harvest energy from the PB. Therefore,
is not possible to occur because the energy saved in the storage is not possible to decrease if energy harvesting is occurred. For the case with
and
, we know that there would be an energy increment of
only when the harvested energy
satisfying
, which results in
with
given in Equation (
4). On the contrary, the event of
and
will occur as long as
, which leads to
.
As such, the transition probability from to within one transition is summarized as
We note that the CDF of
in Equation (
14) could be readily derived from Equation (
9) by making an appropriate replacement, i.e.,
,
,
,
. This can be easily concluded because the weight vector of MRC
has the similar form with the weight vector of MRT.
Denote
as the transition matrix with its
th element being
. It is easy to know that
is irreducible and row stochastic. Hence, the stationary probability vector can be derived as [
14,
15]
where
,
is the identity matrix, and
is an all-ones matrix. Hence,
can be then derived.
3.3. Sensor Transmit Power Optimization
From the overall ergodic capacity definition in Equation (
7), we see that there exists a trade-off between the value of
and the overall ergodic capacity of the network. On the one hand, increasing
leads to decreased
, and then degrades overall ergodic capacity. On the other hand, increasing
also results in increased
, and thus is beneficial to improve the overall ergodic capacity at the same time. As a result, there exists an optimum sensor transmit power
to maximize the overall ergodic capacity of the considered network. Mathematically, the optimization of
to maximize the overall ergodic capacity can be modeled as
Note that the closed-form expression of is intractable. As an alternative, it can be solved by applying the exhaustive method conveniently because it is an one-dimensional problem and the argument only takes finite values.
4. Numerical Results
In this section, we present the numerical results to illustrate the impacts of various system parameters on the performance of the proposed network. Without loss of generality, we set dBm, mJ, , , , , , m, m, s, and dBm, unless otherwise stated.
Figure 2 and
Figure 3 plot the transmit probability
of the proposed network versus the channel correlation factor between
B and
S,
, and the sensor transmit power
, respectively. We note that
is required to conduct MRT operation, and the case of
is presented as a benchmark when MRT is not applied. As can be expected, the energy transmission is greatly improved due to MRT with
increases, so that
is largely improved. However,
will be significantly degraded with decreasing
. Specifically, the MRT operation does not bring any benefit when
regardless of the value of
. This is because there is little correlation between
and
with
. Besides, we find from the results in
Figure 3 that
will severely decrease with the increase of
. The reason is that a larger transmit power is generally more difficult to be satisfied for the wireless powered sensor.
Figure 4 and
Figure 5 plot the effective ergodic capacity
of the proposed network versus the channel correlation factor between
S and
D,
, and the sensor transmit power
, respectively. Similarly, the effective ergodic capacity
can be enhanced by increasing the number of antennas at
D, which however degrades with the decrease of
. Specifically, the MRC becomes invalid when
regardless of the value of
. Moreover, it is observed in
Figure 5 that
can be improved by increasing the sensor transmit power, because a larger sensor transmit power generally leads to a greater received SNR at
D.
Figure 6 and
Figure 7 plot the overall ergodic capacity
of the proposed network versus the sensor transmit power
, the number of antennas at
B,
, and the channel correlation factor between
B and
S,
, respectively. In
Figure 6, a trade-off between the value of
and ergodic capacity is observed for various
and
. In addition,
Figure 7 presents the maximized overall ergodic capacity
with different sensor transmit power. Note that, in this figure,
and
correspond to
mW and
mW, respectively. As can be seen, the overall ergodic capacity can be greatly improved with the proposed sensor transmit power optimization.