# Efficient Multiway Relaying for Data Sharing in Energy Harvesting Sensor Networks

Authors: Moslem Noori; Masoud Ardakani
Journal: Journal of Sensors
 (2015)

---

## Abstract

In a wireless sensor network (WSN), sensors often
need to share their measurements for applications like distributed
estimation and detection or data aggregation. Here, we suggest
using multiway relaying (MWR) for data sharing between
energy harvesting sensors that cannot directly communicate with
each other. We first start by studying the achievable data rate of
amplify-and-forward (AF) MWR for energy harvesting sensors.
Then, we show that, by backing off the transmit power at the
sensors, not only better energy efficiency and longer lifetime are
achieved, but also the data sharing rate increases. Based on this
result, we further improve the performance of AF MWR in the
assumed WSN by smartly adjusting the transmit power at the
sensors. Our power allocation is devised in a way to improve
the energy efficiency of MWR and increase the sum rate of data
sharing between the sensors over the network lifetime. Simulation
results are presented to verify the enhancement achieved by using
our proposed power allocation technique.


---

## Body

## 1. Introduction

Sensors often disseminate and share data with each other to accomplish a collaborative task in a wireless sensor network (WSN) [1]. Such collaboration could lead to achieving a common goal in the network or improving its general performance. For instance, by sharing their measurement data, sensors are able to track moving objects in their vicinity or improve their estimation from a parameter to make further decisions [2, 3]. Also, neighboring sensors may share their measurements to form an aggregated version of the data to be sent to the data sink.Sharing data between the sensors can take place in different forms. Flooding [1] is one of the data sharing/dissemination approaches where a sensor broadcasts its own message or a message received from another sensor. This continues until the data is shared with all intended sensors. Flooding, however, could be inefficient due to its high energy consumption, coming from repetitive packet transmission by the sensors and possible collision between the flooded packets [4]. Another strategy is to take a cooperative approach [5, 6] where a central (relay) node is used to coordinate data sharing between the sensors [7]. In this approach, the relay node asks the sensors to transmit their data to it and then it distributes the data between the sensors. This approach is more efficient in terms of energy consumption at the cost of extra hardware and more complex implementation. In addition, several other studies have considered developing new techniques or variations of the flooding and cooperative schemes with better energy efficiency, for example, [8–10].Multiway relaying (MWR) [11] is a new efficient data sharing paradigm for wireless networks. This approach is cooperative and a relay node is used to assist the users (sensors) to share their data. MWR can achieve better bandwidth/energy efficiency compared to conventional cooperative schemes. This is achieved through employingphysical layer network coding [12] that smartly leverages interference between different users’ transmitted signals. Despite many prior studies on the advantage of MWR including its achievable data rates and energy efficiency [11, 13–16], its application for WSNs has not been investigated.In this paper, we study the application of amplify-and-forward (AF) MWR for WSNs. More specifically, we consider a group of energy harvesting sensors that want to share their data with each other. The sensors do not have direct communication links, for example, due to surrounding obstacles; hence, employing a cooperative communication approach is essential. To this end, a relay node is used to assist the sensors with their data communication. For such a network, we first find the achievable rate of data sharing between sensors for a given time span using MWR. Then, we show that unlike conventional cooperative schemes for data sharing, by backing off the transmit power at the sensors, we can further improve the energy efficiency as well as the data rate of AF MWR.Considering the above, an important question is how to assign the transmit power at the sensors to get the best performance out of the network. To address this, we formulate an optimization problem to maximize the data sharing sum rate in the network while achieving a target lifetime is guaranteed. For this optimization problem, we show that the number of constraints grows exponentially with the number of sensors. As a result, solving it becomes drastically complex as the number of sensors increases. That said, we identify the conditions in the problem that are more likely to be active and formulate a relaxed version of the problem based on them. While the number of constraints in the relaxed problem grows linearly with the number of sensors, deriving a closed form solution for it is still quite involved. Therefore, we propose an algorithm with a reasonable complexity to find a suboptimal set of transmit powers.The performance of our proposed suboptimal solution is investigated via a set of numerical examples. The simulation results verify that our proposed power allocation technique results in a performance improvement over full power transmission. This improvement is more pronounced for smaller number of sensors and higher SNRs.The paper is organized as follows. First, we present the system model and some preliminary materials in Section2. The achievable data rates of AF MWR are studied in Section 3. The proposed power allocation schemes are discussed in Section 4 and numerical results for its performance are presented in Section 5. At the end, Section 6 concluded the paper.


## 2. System Model

Here, we assume a WSN withN
≥
2 sensors denoted by s
1
,
s
2
,
…
,
s
N. The sensors want to periodically share their data, with a period called atime slot, for at least T time slots to accomplish a common task. In-network processing at the sensors and network clustering are example situations where such setup may appear. Note that N does not represent the total number of sensors in the network.To enforce fairness between the sensors, all sensors share data with a common rateR
i
c during the ith time slot. The common rate is the data rate that all sensors can reliably transmit and receive data with. More formally, if R
i
,
j denotes the achievable transmission rate by s
j in the ith time slot, then, R
i
c is(1)
R
i
c
=
min
⁡
{
R
i
,
1
,
R
i
,
2
,
…
,
R
i
,
N
}
.The data sharing sum rate is defined using R
i
c and denotes the overall rate of data sharing between the sensors over a time span. More specifically, the data sharing sum rate until the ith time slot is defined as(2)
R
i
s
=
∑
k
=
1
i
‍
∑
j
=
1
N
‍
R
k
c
=
N
∑
k
=
1
i
‍
R
k
c
.


### 2.1. Data Communication Model

In this work, we assume that the sensors do not have a direct link for communication. This could be a result of a harsh environment with surrounding obstacles limiting the direct communication. For such scenario, we suggest using MWR to facilitate data sharing between the sensors. To this end, a relay (central) node is used to help the sensors with their data communication using AF MWR strategy. The relay node could be a sensor serving as a cluster head or a representative for a group of sensors or a special node with more capabilities than the ordinary sensors.To perform MWR, the time slot is equally divided into two phases, namely,uplink anddownlink. In the uplink phase, all sensors simultaneously transmit their messages to the relay node (Figure 1(a)). Hence, the relay receives a superposition of the sensors transmitted signals often referred to as physically network coded signal [12]. Then, the relay amplifies the received signal and broadcasts it to the sensors in the downlink phase (Figure 1(b)). After receiving the relay’s broadcast signal, each sensor removes self-interference from the broadcast message and decodes the data of all other sensors. For further information on MWR and how it can help improving the data rate compared to conventional one-way relaying, see [11, 15].Figure 1
Demonstration of the uplink and downlink phases.
(a)
Uplink phase
(b)
Downlink phaseThe receiver noise at the relay and the sensors is assumed to be additive white Gaussian with zero mean and unit variance. Further, the channels between the sensors and the relay are reciprocal and the channel gain is fixed for one uplink or downlink phase. Fors
j, h
i
,
j refers to the channel gain between s
j and the relay during the ith time slot.


### 2.2.  Energy Harvesting/Consumption Model

To provide the power needed for data communication, sensors are equipped with a battery with initial energyE
i Joules. In addition, during each time slot, sensors harvest energy to further supply power resource for data communication. More specifically, each sensor is capable of harvesting energy with rate P
h. Please note that the focus of this work is on the energy consumption for data transmission and the energy needed for sending and receiving by the sensors is out of our scope. However, they can be readily incorporated in our analysis.Here, the transmit power ofs
j at the ith round, 1
≤
i
≤
T, is denoted by P
i
,
j and the transmit power of the relay is assumed to be P
r. To guarantee network functionality until T, s
j takes a conservative approach and limits P
i
,
j such that P
i
,
j
≤
E
i
/
T
+
P
h. This guarantees that the sensor has enough energy for data sharing at time T even if the harvested energy is negligible.


## 2.1. Data Communication Model

In this work, we assume that the sensors do not have a direct link for communication. This could be a result of a harsh environment with surrounding obstacles limiting the direct communication. For such scenario, we suggest using MWR to facilitate data sharing between the sensors. To this end, a relay (central) node is used to help the sensors with their data communication using AF MWR strategy. The relay node could be a sensor serving as a cluster head or a representative for a group of sensors or a special node with more capabilities than the ordinary sensors.To perform MWR, the time slot is equally divided into two phases, namely,uplink anddownlink. In the uplink phase, all sensors simultaneously transmit their messages to the relay node (Figure 1(a)). Hence, the relay receives a superposition of the sensors transmitted signals often referred to as physically network coded signal [12]. Then, the relay amplifies the received signal and broadcasts it to the sensors in the downlink phase (Figure 1(b)). After receiving the relay’s broadcast signal, each sensor removes self-interference from the broadcast message and decodes the data of all other sensors. For further information on MWR and how it can help improving the data rate compared to conventional one-way relaying, see [11, 15].Figure 1
Demonstration of the uplink and downlink phases.
(a)
Uplink phase
(b)
Downlink phaseThe receiver noise at the relay and the sensors is assumed to be additive white Gaussian with zero mean and unit variance. Further, the channels between the sensors and the relay are reciprocal and the channel gain is fixed for one uplink or downlink phase. Fors
j, h
i
,
j refers to the channel gain between s
j and the relay during the ith time slot.


## 2.2.  Energy Harvesting/Consumption Model

To provide the power needed for data communication, sensors are equipped with a battery with initial energyE
i Joules. In addition, during each time slot, sensors harvest energy to further supply power resource for data communication. More specifically, each sensor is capable of harvesting energy with rate P
h. Please note that the focus of this work is on the energy consumption for data transmission and the energy needed for sending and receiving by the sensors is out of our scope. However, they can be readily incorporated in our analysis.Here, the transmit power ofs
j at the ith round, 1
≤
i
≤
T, is denoted by P
i
,
j and the transmit power of the relay is assumed to be P
r. To guarantee network functionality until T, s
j takes a conservative approach and limits P
i
,
j such that P
i
,
j
≤
E
i
/
T
+
P
h. This guarantees that the sensor has enough energy for data sharing at time T even if the harvested energy is negligible.


## 3. DerivingR
i
c for AF MWR

In this section, we first present the formal signaling model for AF MWR. Then, we study the achievable common rate of data sharing over an arbitrary time slot; say theith one. Using this analysis, we are able to determine the sum rate of data sharing. Without loss of generality, it is assumed that h
i
,
1
≤
h
i
,
2
≤
⋯
≤
|
h
i
,
N
|.After transmission by all sensors in the uplink phase ofith time slot, relay receives the following superimposed signal as depicted in Figure 1(a):(3)
y
i
r
=
∑
n
=
1
N
‍
h
i
,
n
x
i
,
n
+
n
r
,where x
i
,
n’s are sensors’ messages in the ith time slot and n
r is the Gaussian noise at the relay. Now, the relay forms its message, x
i
r, as x
i
r
=
α
i
y
i
r where(4)
α
i
≤
P
r
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1and broadcasts it to the sensors in the downlink phase as depicted in Figure 1(b). Each sensor then removes its self-interference from the broadcast signal. After doing such, the received signal at s
j is(5)
y
i
,
j
=
α
i
h
i
,
j
∑
n
≠
j
‍
h
i
,
n
x
i
,
n
+
α
i
h
i
,
j
n
r
+
n
j
,where n
j is the receiver noise at s
j. The signal model in (5) is basically similar to the one of a multiple-access (MAC) channel [17]. In this MAC channel, s
j wants to decode the message of all other sensors, where α
i
h
i
,
j
n
r
+
n
j is the effective noise. Achievable data rates for fading MAC channels have been studied in prior studies, for example, [18]. Please note that in conventional MAC channels, unlike MWR, there is only one communication hop between the sources and the destination. That said, the results in [18] can still be used to find achievable data rates in our WSN setup. More specifically, it can be shown that the data of all other sensors can be reliably decoded at s
j if their transmit data rates satisfy the following constraint:(6)
∑
n
∈
S
j
‍
R
i
,
n
≤
1
2
log
⁡
2
1
+
α
i
2
h
i
,
j
2
∑
n
∈
S
j
P
i
,
n
h
i
,
n
2
1
+
α
i
2
h
i
,
j
2
,for all nonempty S
j
⊆
1,2
,
…
,
N
∖
j. Note that for decoding other sensors’ messages, each sensor should have the knowledge of all channel gains. The relay node can obtain this information via pilot signaling and then pass it to the sensors. Looking at (6), one can identify another difference between conventional MAC channels and the considered WSN with MWR. Here, the achievable data rates at the receiver, s
j, depend not only on the other users’ transmit power but also on s
j transmit power through α
i. This in turn leads to significant interdependency between the allocation of the users’ transmit powers and the users’ achievable data rates.To find the achievable rate at other sensors, one can follow a model similar to (5). As a result, for a rate tuple R
=
(
R
i
,
1
,
R
i
,
2
,
…
,
R
i
,
N
) to be achievable at all sensors during time slot i, N constraints similar to (6) should be satisfied by R. Since (6) in fact refers to 2
N
-
1
-
1 constraints at s
j to cover all possible nonempty S
j, this means that the rate tuple R is achievable if N
(
2
N
-
1
-
1
) constraints are satisfied.Before finding the common rate for data sharing in theith time slot, let us define R
i
,
j
c as the achievable common rate at s
j referring to the rate that all other sensors can reliably transmit their data to s
j with. To find R
i
,
j
c, we set all R
i
,
n in (6) equal to R
i
,
j
c resulting in ∑
n
∈
S
j
R
i
,
n
=
|
S
j
|
R
i
,
j
c, where |
S
j
| is the cardinality of S
j. Thus,(7)
R
i
,
j
c
≤
min
⁡
S
j
R
i
,
j
S
j
,where(8)
R
i
,
j
S
j
=
1
2
S
j
log
⁡
2
1
+
α
i
2
h
i
,
j
2
∑
n
∈
S
j
P
i
,
n
h
i
,
n
2
1
+
α
i
2
h
i
,
j
2
.Now, using the common rate definition (1), we have(9)
R
i
c
≤
min
⁡
j
min
⁡
S
j
1
2
S
j
log
⁡
2
1
+
α
i
2
h
i
,
j
2
∑
n
∈
S
j
P
i
,
n
h
i
,
n
2
1
+
α
i
2
h
i
,
j
2
.Note that for a given set of power at the sensors, the common rate is maximum when the condition in (9) holds with equality. Furthermore, we have the following propositions based on our analysis.Proposition 1.
R
i
c is an increasing functions of α
i and as a consequence an increasing function of P
r.In the rest of the paper, we assume that the relay injects enough power into the broadcast signal such that (4) is satisfied with equality to increase the data sharing rate.Proposition 2.
For anyj, R
i
,
j
c is an increasing function of P
i
,
n for all n
≠
j and a decreasing function of P
i
,
j.As we discussed before, findingR
i
c originally requires checking N
(
2
N
-
1
-
1
) constraints. However, these conditions are not independent and all of them depend on h
i
,
1
,
h
i
,
2
,
…
,
h
i
,
N and the sensors transmit power. The following theorem indicates how their dependency can be exploited to reduce the number of conditions to be checked for finding R
i
c. This, as a result, decreases the complexity of calculating the sum rate.Theorem 3.
For any arbitrary sensors
j, R
i
,
j
c is bounded by those R
i
,
j
S
j’s where S
j
≥
j
-
1 and {
1,2
,
…
,
j
-
1
}
⊆
S
j.Proof.
See Appendix.Corollary 4.
R
i
c can be found by checking only 2
N
-
2 constraints rather than all possible N
(
2
N
-
1
-
1
) constraints.Proof.
See  Appendix.Having the achievable common rate of each time slot, the achievable data sharing sum rate is derived using (2).


## 4. Rate-Improving Power Allocation for AF MWR

As discussed earlier, the achievable common rate within each time slot and consequently the overall sum rate depend on the sensors’ transmit power. Thus, it is desired to seek for a power allocation at the sensors that give the best performance, meaning achieving the highest data sharing sum rate, while keeping the network operational for at leastT time slots. The following theorem helps us to understand the features of such power allocation.Theorem 5.
Full power transmission at the sensors, that is,P
i
,
j
=
E
i
/
T
+
P
h for all s
j and i, does not guarantee achieving maximum R
i
c.Proof.
See Appendix.The point stated in Theorem5 is different from what we see in conventional cooperative data sharing schemes based on one-way relaying and time sharing. That is, while full power transmission guarantees maximum data rate in conventional cooperative schemes, Theorem 5 implies that we can use an MWR strategy where sensors do not transmit with full power and gain a twofold benefit. That is to improve both data rate and energy efficiency compared to an MWR strategy sending with full power at the sensors. This gives a guideline to further improve MWR performance for WSN applications. Now, the question is if full power transmission does not give the best performance, what power allocation reaches the best performance, that is highest sum rate while achieving a lifetime of at least T time slots.To address the above question, we formulate an optimization problem. The goal of this optimization problem is to maximize the sum rate while guaranteeing that the network lifetime surpassesT time slots as follows:(10)
max
⁡
P
1,1
,
…
,
P
l
,
N
R
l
s
,
(11)
R
l
s
=
N
∑
i
=
1
l
‍
R
i
c
,
(12)
R
i
c
≤
min
⁡
j
min
⁡
S
j
1
2
S
j
log
⁡
2
1
+
α
i
2
h
i
,
j
2
∑
n
∈
S
j
P
i
,
n
h
i
,
n
2
1
+
α
i
2
h
i
,
j
2
,
(13)
T
≤
l
,
(14)
∑
k
=
1
i
‍
P
k
,
j
≤
E
i
+
i
P
h
,
∀
i
≤
l
,
∀
j
∈
{
1
,
…
,
N
}
,
(15)
0
≤
P
k
,
j
,
∀
k
≤
l
,
∀
j
∈
{
1
,
…
,
N
}
.In this optimization problem, the constraint in (13) ensures that the network lifetime is at least T time slots. Also, (14) guarantees that the sensors’ consumed energy for data transmission at any time slot is always less than the sum of the energy provided by the initial battery and the harvested energy.Based on the statement in Theorem3, the above optimization problem has at least 2
N
+
2
N
T constraints that grow exponentially with the number of sensors. That said, the optimization problem is quite difficult to solve especially for large N. In the following, we discuss how one can find a suboptimal solution to the optimization problem in (10) with significantly lower complexity.


### 4.1. Relaxing Problem

As we discussed in Section2, an arbitrary sensor s
j limits its transmit power such that P
i
,
j
≤
E
i
/
T
+
P
h. This guarantees that the sensor never runs out of power before T. Using this point, we relax the optimization problem in (10) and rather focus on maximizing R
i
c for any arbitrary i
≤
T assuming P
i
,
j
≤
E
i
/
T
+
P
h. Solving this relaxed optimization problem for all time slots, i
≤
T can give a suboptimal solution for the original sum rate maximization problem. To this end, we formulate an optimization problem to maximize the common rate at time slot i, R
i
c, as follows:(16)
max
⁡
P
i
,
1
,
…
,
P
i
,
N
min
⁡
R
i
,
j
c
,
P
i
,
j
≤
E
i
T
+
P
h
,
∀
j
,
0
≤
P
i
,
j
,
∀
j
.Using Fact 1, the above max-min problem is transformed into(17)
max
⁡
P
i
,
1
,
…
,
P
i
,
N
R
i
th
,
(18)
P
r
h
i
,
j
2
∑
n
∈
S
j
P
i
,
n
h
i
,
n
2
P
r
h
i
,
j
2
+
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1
≥
2
2
|
S
j
|
R
i
th
-
1
∀
j
,
S
j
,
P
i
,
j
≤
E
i
T
+
P
h
,
∀
j
,
0
≤
P
i
,
j
,
∀
j
.The number of constraints in (17) still grows exponentially with the number of sensors. To this end, we further simplify rate maximization problem by relaxing some of the rate constraints. More specifically, we consider the rate constraints, where |
S
j
|
=
N
-
1 for any j, that is, N
-
1-tuple constraints. These constraints have the smallest prelog factor, that is, 1
/
(
N
-
1
), among constraints presented in (6). Thus, in medium to high SNR range, they are more likely to limit R
i
c. Without loss of generality, one could also assume that during the ith time slot, h
i
,
1
≤
h
i
,
2
≤
⋯
≤
|
h
i
,
N
|. The relaxed problem also considers all single rate constraints at the worst user, that is, s
1, and the single rate constraint on the data rate of the worst user at the second worst user, that is, R
i
,
1 at s
2, as follows:(19)
max
⁡
P
i
,
1
,
…
,
P
i
,
N
R
i
th
,
(20)
P
r
h
i
,
j
2
∑
n
=
1
,
n
≠
j
N
P
i
,
n
h
i
,
n
2
P
r
h
i
,
j
2
+
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1
≥
2
2
(
N
-
1
)
R
i
th
-
1
∀
j
,
(21)
P
r
h
i
,
1
2
P
i
,
j
h
i
,
j
2
P
r
h
i
,
1
2
+
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1
≥
2
2
(
N
-
1
)
R
i
th
-
1
∀
j
≠
1
,
(22)
P
r
h
i
,
2
2
P
i
,
1
h
i
,
1
2
P
r
h
i
,
2
2
+
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1
≥
2
2
(
N
-
1
)
R
i
th
-
1
,
(23)
P
i
,
j
≤
E
i
T
+
P
h
,
∀
j
,
(24)
0
≤
P
i
,
j
,
∀
j
.To understand the intuition behind considering (21), recall that, for any k
≠
1, when S
k
=
{
j
} and j
∉
{
1
,
k
}, one can show that R
i
,
k
S
k
≥
R
i
,
1
S
1, where S
1
=
{
j
}. As a result, if (21) is satisfied, it is ensured that the smallest bound on s
j
’s single rate for any 2
≤
j has been taken into account. In a similar way, one can argue that (22) represents the smallest bound on the s
1
’s single rate. Thus, the relaxed optimization problem in (19) considers all N
-
1-tuple rate constraints as well as the worst (smallest) bounds on the single rates. This relaxation significantly decreases the number of constraints to 4
N which grows linearly with the number of sensors.


### 4.2. Proposed Suboptimal Solution for the Relaxed Problem

Although the relaxed optimization problem has significantly less number of constraints compared to the original problem, finding an analytical solution for it is quite involved. To address this challenge, we suggest a method to find a numerical suboptimal solution for the relaxed problem. The logic behind the proposed numerical solution is similar to the approach taken in the proof of Theorem5.The proposed algorithm for finding suboptimal power allocation is presented in Algorithm1. The proposed algorithm starts from an initial point for power allocations by assuming full power transmission for all sensors and then uses the point stated in Proposition 2 to further improve the rate. To this end, we find k
=
argmi
n
j
R
i
,
j
c, where R
i
,
j
c’s are calculated based on merely the relaxed constraints, that is, single rate constraints at s
1 and the constraint on s
1’s rate at s
2 plus all N
-
1-tuple constraints. We name s
k the bottleneck sensor since it limits the common rate R
i
c. Now, based on the statement in Proposition 2, R
i
,
k
c is a decreasing function of P
i
,
k. Thus, if we back off P
i
,
k, say with step ϵ
p, R
i
,
k
c sees an increase at the cost of a rate decrease for all other sensors. Backing off P
i
,
k continues to the point where s
k is not bottleneck anymore and another sensor, for example, s
l, limits the rate and becomes the bottleneck. Now, if s
l backs off its power to increase R
i
,
l
c, R
i
,
k
c sees a decline. However, if the transmit power at both s
k and s
l decreases with step ϵ
p, R
i
,
k
c and R
i
,
l
c could both see an increase. Backing off the power at s
k and s
l continues till another sensor becomes the bottleneck. The algorithm continues looking for bottleneck sensors and backing off their power up to the point where no further improvement at R
i
c can be achieved. Note that, in Algorithm 1, B is a set containing all sensors that have been a bottleneck sensor at some point. In our proposed algorithm, the power allocation at the sensors always moves in a direction where R
i
c, calculated based on the relaxed constraints, increases. As a result, the algorithms climb up toward local maxima.Algorithm 1:Algorithmic presentation of the proposed power allocation technique.
(1) Assign full power to all sensors by settingP
i
,
j
=
(
E
i
/
T
)
+
P
h
(2) Choose a power deduction stepϵ
p
(3) SetB
=
{
}
(4)repeat
(5)    CalculateR
i
,
j
c for all j based on the relaxed constraints
(6)    FindR
i
c
=
min
⁡
R
i
,
j
c and k
=
argmin
R
i
,
j
c
(7)B
=
B
∪
{
k
}
(8)    For allj
∈
B, P
i
,
j
=
P
i
,
j
-
ϵ
p
(9)until  R
i
c shows no improvement compared to the previous stepWhen the algorithm stops, it gives a set of power allocation denoted by(
P
i
,
1
∗
,
P
i
,
2
∗
,
…
,
P
i
,
N
∗
) giving a suboptimal value for R
i
c when only relaxed constraints are considered. This is not, however, the actual value for the common rate in time slot i. The actual R
i
c is found by considering all 2
N
-
2 rate constraints with (
P
i
,
1
∗
,
P
i
,
2
∗
,
…
,
P
i
,
N
∗
) as the set of transmit powers at the sensors. Please note that there is a performance gap between the data rate achieved through the power allocation given by the suboptimal algorithm and the optimal solution. However, finding the optimal solution has an exponential complexity compared to the linear complexity of the proposed suboptimal algorithm. Thus, finding the optimal solution is not feasible in a reasonable time.At time sloti and for an arbitrary s
j, a
i
,
j
=
P
i
,
j
∗
-
E
i
/
T
-
P
h represents the energy saving that is achieved through our proposed power allocation. If a
i
,
j
>
0, s
j stores the saved energy. At the end of T time slots, the total saved energy is A
j
=
∑
i
=
1
T
‍
a
i
,
j. Using this energy, s
j extends its lifetime beyond T and continues sharing data with other sensors in the same fashion as up to T. Let us denote the time slot when the first sensor runs out of energy by T
∗
≥
T. Now, the overall sum rate of AF MWR with our proposed power allocation is(25)
R
T
∗
s
=
N
∑
i
=
1
T
∗
‍
R
i
c
.Note that when full power transmission is in effect, all sensors deplete their batteries and run out of energy at the end of the Tth time slot.


## 4.1. Relaxing Problem

As we discussed in Section2, an arbitrary sensor s
j limits its transmit power such that P
i
,
j
≤
E
i
/
T
+
P
h. This guarantees that the sensor never runs out of power before T. Using this point, we relax the optimization problem in (10) and rather focus on maximizing R
i
c for any arbitrary i
≤
T assuming P
i
,
j
≤
E
i
/
T
+
P
h. Solving this relaxed optimization problem for all time slots, i
≤
T can give a suboptimal solution for the original sum rate maximization problem. To this end, we formulate an optimization problem to maximize the common rate at time slot i, R
i
c, as follows:(16)
max
⁡
P
i
,
1
,
…
,
P
i
,
N
min
⁡
R
i
,
j
c
,
P
i
,
j
≤
E
i
T
+
P
h
,
∀
j
,
0
≤
P
i
,
j
,
∀
j
.Using Fact 1, the above max-min problem is transformed into(17)
max
⁡
P
i
,
1
,
…
,
P
i
,
N
R
i
th
,
(18)
P
r
h
i
,
j
2
∑
n
∈
S
j
P
i
,
n
h
i
,
n
2
P
r
h
i
,
j
2
+
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1
≥
2
2
|
S
j
|
R
i
th
-
1
∀
j
,
S
j
,
P
i
,
j
≤
E
i
T
+
P
h
,
∀
j
,
0
≤
P
i
,
j
,
∀
j
.The number of constraints in (17) still grows exponentially with the number of sensors. To this end, we further simplify rate maximization problem by relaxing some of the rate constraints. More specifically, we consider the rate constraints, where |
S
j
|
=
N
-
1 for any j, that is, N
-
1-tuple constraints. These constraints have the smallest prelog factor, that is, 1
/
(
N
-
1
), among constraints presented in (6). Thus, in medium to high SNR range, they are more likely to limit R
i
c. Without loss of generality, one could also assume that during the ith time slot, h
i
,
1
≤
h
i
,
2
≤
⋯
≤
|
h
i
,
N
|. The relaxed problem also considers all single rate constraints at the worst user, that is, s
1, and the single rate constraint on the data rate of the worst user at the second worst user, that is, R
i
,
1 at s
2, as follows:(19)
max
⁡
P
i
,
1
,
…
,
P
i
,
N
R
i
th
,
(20)
P
r
h
i
,
j
2
∑
n
=
1
,
n
≠
j
N
P
i
,
n
h
i
,
n
2
P
r
h
i
,
j
2
+
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1
≥
2
2
(
N
-
1
)
R
i
th
-
1
∀
j
,
(21)
P
r
h
i
,
1
2
P
i
,
j
h
i
,
j
2
P
r
h
i
,
1
2
+
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1
≥
2
2
(
N
-
1
)
R
i
th
-
1
∀
j
≠
1
,
(22)
P
r
h
i
,
2
2
P
i
,
1
h
i
,
1
2
P
r
h
i
,
2
2
+
∑
n
=
1
N
P
i
,
n
h
i
,
n
2
+
1
≥
2
2
(
N
-
1
)
R
i
th
-
1
,
(23)
P
i
,
j
≤
E
i
T
+
P
h
,
∀
j
,
(24)
0
≤
P
i
,
j
,
∀
j
.To understand the intuition behind considering (21), recall that, for any k
≠
1, when S
k
=
{
j
} and j
∉
{
1
,
k
}, one can show that R
i
,
k
S
k
≥
R
i
,
1
S
1, where S
1
=
{
j
}. As a result, if (21) is satisfied, it is ensured that the smallest bound on s
j
’s single rate for any 2
≤
j has been taken into account. In a similar way, one can argue that (22) represents the smallest bound on the s
1
’s single rate. Thus, the relaxed optimization problem in (19) considers all N
-
1-tuple rate constraints as well as the worst (smallest) bounds on the single rates. This relaxation significantly decreases the number of constraints to 4
N which grows linearly with the number of sensors.


## 4.2. Proposed Suboptimal Solution for the Relaxed Problem

Although the relaxed optimization problem has significantly less number of constraints compared to the original problem, finding an analytical solution for it is quite involved. To address this challenge, we suggest a method to find a numerical suboptimal solution for the relaxed problem. The logic behind the proposed numerical solution is similar to the approach taken in the proof of Theorem5.The proposed algorithm for finding suboptimal power allocation is presented in Algorithm1. The proposed algorithm starts from an initial point for power allocations by assuming full power transmission for all sensors and then uses the point stated in Proposition 2 to further improve the rate. To this end, we find k
=
argmi
n
j
R
i
,
j
c, where R
i
,
j
c’s are calculated based on merely the relaxed constraints, that is, single rate constraints at s
1 and the constraint on s
1’s rate at s
2 plus all N
-
1-tuple constraints. We name s
k the bottleneck sensor since it limits the common rate R
i
c. Now, based on the statement in Proposition 2, R
i
,
k
c is a decreasing function of P
i
,
k. Thus, if we back off P
i
,
k, say with step ϵ
p, R
i
,
k
c sees an increase at the cost of a rate decrease for all other sensors. Backing off P
i
,
k continues to the point where s
k is not bottleneck anymore and another sensor, for example, s
l, limits the rate and becomes the bottleneck. Now, if s
l backs off its power to increase R
i
,
l
c, R
i
,
k
c sees a decline. However, if the transmit power at both s
k and s
l decreases with step ϵ
p, R
i
,
k
c and R
i
,
l
c could both see an increase. Backing off the power at s
k and s
l continues till another sensor becomes the bottleneck. The algorithm continues looking for bottleneck sensors and backing off their power up to the point where no further improvement at R
i
c can be achieved. Note that, in Algorithm 1, B is a set containing all sensors that have been a bottleneck sensor at some point. In our proposed algorithm, the power allocation at the sensors always moves in a direction where R
i
c, calculated based on the relaxed constraints, increases. As a result, the algorithms climb up toward local maxima.Algorithm 1:Algorithmic presentation of the proposed power allocation technique.
(1) Assign full power to all sensors by settingP
i
,
j
=
(
E
i
/
T
)
+
P
h
(2) Choose a power deduction stepϵ
p
(3) SetB
=
{
}
(4)repeat
(5)    CalculateR
i
,
j
c for all j based on the relaxed constraints
(6)    FindR
i
c
=
min
⁡
R
i
,
j
c and k
=
argmin
R
i
,
j
c
(7)B
=
B
∪
{
k
}
(8)    For allj
∈
B, P
i
,
j
=
P
i
,
j
-
ϵ
p
(9)until  R
i
c shows no improvement compared to the previous stepWhen the algorithm stops, it gives a set of power allocation denoted by(
P
i
,
1
∗
,
P
i
,
2
∗
,
…
,
P
i
,
N
∗
) giving a suboptimal value for R
i
c when only relaxed constraints are considered. This is not, however, the actual value for the common rate in time slot i. The actual R
i
c is found by considering all 2
N
-
2 rate constraints with (
P
i
,
1
∗
,
P
i
,
2
∗
,
…
,
P
i
,
N
∗
) as the set of transmit powers at the sensors. Please note that there is a performance gap between the data rate achieved through the power allocation given by the suboptimal algorithm and the optimal solution. However, finding the optimal solution has an exponential complexity compared to the linear complexity of the proposed suboptimal algorithm. Thus, finding the optimal solution is not feasible in a reasonable time.At time sloti and for an arbitrary s
j, a
i
,
j
=
P
i
,
j
∗
-
E
i
/
T
-
P
h represents the energy saving that is achieved through our proposed power allocation. If a
i
,
j
>
0, s
j stores the saved energy. At the end of T time slots, the total saved energy is A
j
=
∑
i
=
1
T
‍
a
i
,
j. Using this energy, s
j extends its lifetime beyond T and continues sharing data with other sensors in the same fashion as up to T. Let us denote the time slot when the first sensor runs out of energy by T
∗
≥
T. Now, the overall sum rate of AF MWR with our proposed power allocation is(25)
R
T
∗
s
=
N
∑
i
=
1
T
∗
‍
R
i
c
.Note that when full power transmission is in effect, all sensors deplete their batteries and run out of energy at the end of the Tth time slot.


## 5. Simulation Results

In this section, we present several numerical examples to evaluate the performance of our proposed power allocation. To this end, we compare the performance of AF MWR with the proposed power allocation with the performance of AF MWR with full power transmission. In our simulation results, it is assumed that the channel gains between the sensors and the relay follow a Rayleigh distribution with average power 1. To take into account the effect of channel gains randomness, the results are averaged over 100 realizations forT
=
50. For a given transmit SNR, we set E
i
=
1
0
SNR
/
10
T to ensure transmitting with the specified SNR over all time slots. In addition, P
r
=
1
0
SNR
/
10 and ϵ
p
=
0.001 SNR. Note that, here, we present actual rates by checking all 2
N
-
2 constraints for the set of transmit powers.Figure2 depicts the normalized sum rate of AF MWR with our proposed allocation for different values of N. In this figure, P
h
=
10dB and the normalized rate is derived by dividing the sum rate of our proposed technique to that of full power transmission. As seen, the improvement is more pronounced for larger SNRs and smaller N. It reaches around 18 percent for N
=
3 and SNR = 30 dB. The reason is that, for smaller SNRs, it is more likely that the constraints omitted in the relaxed problem actually limit the data rate.Figure 2
Normalized achievable common rate over SNR.In addition, in this figure, we compare our suboptimal results with the numerical solution for the original (nonrelaxed) optimization problem presented in (10)–(15). To find the numerical solution, we have used a Monte Carlo search method during each time slot. More specifically, for time slot i, we generate 1000 random N-tuples for (
P
i
,
1
,
P
i
,
2
,
…
,
P
i
,
N
). Then, after checking all rate conditions, the N-tuple that results in the maximum R
i
c is chosen as the optimal solution. As seen in Figure 2, there is a significant gap between the Monte Carlo numerical solution and our proposed suboptimal solution. However, one should keep in mind that this improvement comes at the price of significantly longer simulation run time. The performance gap between the suboptimal solution and the numerical solution decreases with increasing SNR since the constraints considered in the relaxed optimization problem are more likely to be active in the original optimization problem as well.The effect of the number of sensors on the performance of the proposed power allocation is captured with more detail in Figure3. As seen, the proposed power allocation results in a more noticeable improvement for smaller N. Note that when N increases, it is more likely to have a small channel gain for one of the links limiting the data rate through constraints not considered in the relaxed optimization problem.Figure 3
Normalized achievable common rate of the suboptimal solution for differentN.Figure4 presents the normalized rate of the proposed power allocation versus the energy harvest rate P
h. Higher P
h means having more power resources available for transmission that in turn increases the effective SNR, that is, provided by the initial battery and the harvested power. As a result, we observe a similar behavior to Figure 2 where rate improvement is more visible for higher SNRs.Figure 4
Normalized achievable common rate of the suboptimal solution for different energy harvesting rateP
h.


## 6. Conclusion

In this paper, we suggested the application of MWR for data sharing in an energy harvesting WSN. First, we studied the achievable common data rate for each round of data sharing and the overall data sharing sum rate between the sensors. Using our rate study, we then formulated an optimization problem to maximize the data sharing sum rate while a desired network lifetime is assured. Due to the complexity of this optimization problem, we then suggested a relaxed version of the problem to adjust the power allocation at the sensors for an improved data sharing sum rate. Then, a suboptimal solution for the relaxed optimization problem was proposed. Simulations results verified that our suboptimal solution indeed enhances the data sharing sum rate. The level of the rate improvement depends on the number of sensors as well as SNR. Alternatively, one can use MWR and focus on extending the lifetime rather than the data rate. To solve this problem, an approach similar to the proposed power allocation in this paper can be taken.



---

*Source: 285056-2015-04-15.xml*