The goal of this paper is to minimize the network energy consumption with the aid of the DT-SWIPT scheme. This problem is termed D-ECm in the paper. In this section, we will construct a non-linear programming model for the D-ECm and use an SQP algorithm to solve this optimization problem.
3.1. Problem Formulation
According to the HTT protocol, at the beginning of a period, the HAP will transfer energy to all wireless stations. The power used by the HAP is denoted as
. There is an upper limit,
, for the transmission power. The constraint of
is shown in Equation (
1).
In the ET phase, the channel gain between the HAP (
) and
is denoted as
, which can be expressed in Equation (
2).
where
is the distance between the HAP (
) and
. Other parameters, such as antenna gain, antenna height of the HAP, etc., are neglected in the paper and termed
. For simplicity,
is assumed to be a constant in the paper [
23].
is the path loss exponent [
25].
With the channel gain in Equation (
2), when the HAP transfers energy in the ET phase, the received signal power of
is denoted as
, which can be expressed in Equation (
3).
By Equation (
3), the energy that
can obtain in the ET phase is denoted as
, which can be expressed in Equation (
4).
where
is the energy conversion efficiency of
when receiving energy during energy harvesting, and
.
is the duration of the ET phase.
In addition to the energy required to transmit information to the HAP, wireless stations also need the energy for basic operations, such as receiving beacons, modulating information, and maintaining basic circuits. Therefore, among the obtained energy
, the energy of
will have a proportion of
for transmission, and the remaining
energy will be used for basic operations, where
. Thus, through Equation (
4), the transmission power of
, denoted as
, can be expressed as Equation (
5).
According to the DT-SWIPT scheme,
completes the transmission of information during
and then shares the energy with other wireless stations during the remaining time, that is,
. Let
be the channel gain coefficient from
to
(assuming that
is transmitted earlier than
), which can be expressed in Equation (
6).
where
is the distance between
and
.
is a constant and stands for the other parameters, such as the antenna gain, antenna height of
, like
in Equation (
2).
Therefore, with the channel gain coefficient
, the additional energy shared from
and obtained by
is denoted as
, which can be expressed as Equation (
7).
By the TDMA protocol, wireless stations transmit information to the HAP sequentially. Therefore, each wireless station receives a different amount of energy shared by previously transmitted wireless stations. For this situation, let
be the transmission order of
, and
be the index of the wireless station whose transmission order is
k. That is, if
is the
kth transmitting wireless station, then
and
. If
is the
mth transmitting wireless station (
), then
will obtain the energy shared by
. Therefore, the total additional energy obtained by
, which is shared from the previously transmitted wireless stations, is denoted as
, which can be expressed in Equation (
8).
In addition to the energy obtained from the HAP,
also obtains additional energy shared by previously transmitted wireless stations. Therefore, the total energy obtained by
is denoted as
, which can be expressed in Equation (
9).
Since the first transmitting wireless station does not obtain additional energy, the wireless station with
does not need to calculate the additional energy. As a result, the transmission power of
will be replaced from Equation (
5) to Equation (
10).
Let
be the channel gain coefficient from
to the HAP (
), which can be expressed in Equation (
11).
While
starts to transmit information to the HAP, with the channel gain coefficient
, the received signal strength of the HAP from
is denoted as
, which can be expressed in Equation (
12).
To ensure that the HAP can decode the signal transmitted from
, the ratio of the received signal power
to the noise power
N must be larger than or equal to a certain threshold, which is denoted as
. The relationship is shown in Equation (
13).
According to the Shannon’s Capacity Theorem [
26], the channel capacity between HAP and
, which is denoted as
, can be expressed in Equation (
14).
Obviously,
will be affected by
and the bandwidth
B. According to Equation (
14), the throughput of
during the period
T is denoted as
, which can be expressed in Equation (
15).
Since the work performed by each wireless station may be different,
has its own throughput demand,
. The throughput of
must be larger than or equal to
. The relationship between the
and
is shown in Equation (
16).
This paper aims to minimize the energy consumption of the network under the conditions of Equations (
13) and (
16), for
i from 1 to
n. Since the energy source of wireless stations is harvested from the HAP, which can be considered as the energy consumption of the network. The energy consumption of the network is denoted as
, which can be expressed in Equation (
17).
From [
20], it shows that when
, the energy consumption of the network will be the lowest. Therefore, let
be
. By Equation (
17), the goal will be to minimize the WET time, that is,
.
Based on the goal and constraints, a non-linear programming model for the D-ECm problem can be formulated in Model 1.
| τ0 |
subject to | C1 :Equation (13), i = 1,…,n |
| C2 :Equation (16), i = 1,…,n |
| C3 :0 ≤ τ0, τi ≤ T, i = 1,…,n |
| C4 :τ0 + τi ≤ T |
| C5 :0 ≤ αi ≤ 1, i = 1,…,n − 1 |
In Model 1, denotes the SNR constraint. That is, the ratio of the received signal power to the noise power N must be greater than to ensure that the HAP can decode the signal transmitted from . denotes the throughput constraint, which requires that the throughput of must meet its throughput demand . denotes that the transmission time of the HAP and can not be longer than T and should be larger than 0. denotes that the sum of the transmission time of the HAP and wireless stations can not be longer than T. denotes that the DT-SWIPT ratio of during the transmission phase should be smaller than 1.
Since Model 1 is a nonlinear optimization problem, we use an SQP (Sequential Quadratic Programming) method [
27] to obtain the optimal solution.
3.2. SQP-Based Algorithm for Finding the Optimal Solution of the D-ECm Problem
SQP [
27] is an iterative method for solving constrained nonlinear optimization problems. SQP obtains the optimal solution by solving a sequence of QP subproblems. The standard form of the QP problem is formulated in Model 2. In Model 2,
and
represent the gradients of the function
and the transpose matrix of
, respectively.
M and
N represent the numbers of constraints of the equations and inequalities, respectively. The QP problem can get the corresponding optimal solution
w according to the input
. We can then take
into the original QP problem to form a new QP problem and obtain another optimal solution
. Therefore,
w can be regarded as the iterative step in SQP. When
approaches 0, the corresponding
is the optimal solution of the SQP problem.
|
subject to | |
|
|
To satisfy the standard form of the QP problem, we turn the objective function of Model 1 into minimize , and let , where and represent the set of non-negative real numbers. As a result, we rewrite Model 1 as Model 3.
| |
subject to | | |
C6 : | i = 1,…, n, |
C7 : | i = 1,…, n, |
C8 : | |
C9 : | i = 1,…, n, |
C10 : | |
C11 : | i = 1,…, n, n − 1 |
Model 3 can be converted into a QP problem, called QP-ECm, shown in Model 4. In Model 4, is the optimal solution corresponding to , . k is an iteration counter of the SQP.
|
subject to | | |
C12: | i = 1,…,n, |
C13: | i = 1,…,n, |
C14: | |
C15: | i = 1,…,n, |
C16: | |
C17: | i = 1,…,n − 1 |
To solve the D-ECm problem, an algorithm based on SQP is proposed, called DaTA (
DT-SWIPT
assisted
Time
Allocation), which is shown in Algorithm 1. In line 2, we convert the ECm problem into a QP-ECm problem. In line 3, set the iterative counter
k and select an initial point
that meets all constraints. In line 5, solve the QP-ECm problem according to
and get the optimal solution
. In lines 6 and 7, increase the iterative counter
k, and let
. If
approaches 0, return the optimal solution
; otherwise, return to line 5 to solve a new QP problem again.
Algorithm 1: DaTA: a SQP-based Algorithm for the D-ECm Problem |
|
In this section, we have formulated a nonlinear optimization problem, called D-ECm, and found the optimal solution through a SQP-based algorithm, named DaTA. However, it costs much time to find the optimal solution. However, to our best knowledge, there is no heuristic method proposed in the literature to find a plausible solution. Therefore, we propose a heuristic method, called DaTA-H, and describe it in the following section.