2. System Model
2.1. Distributed Antenna System
This paper takes a look at a DAS-assisted multi-cell system, as shown in
Figure 1, where each cell contains a BS and
S RRUs.
Figure 1 depicts a scenario with
. RRUs are located near the cell edge and have their own coverage regions that are subsets of the entire cell region. BS can provide service to all devices in the coverage region, whereas RRUs can only support users in their own coverage regions. We assume that
K cells in the vicinity of the target cell may interfere with the target cell. Here, the target cell index
k is set to 0, and we define
. Similarly, the index of the BS
s is set to 0, and we have
.
We look at two user scenarios: the cell-center user. who is only served by the BS, and the cell-edge user, who receives service not only from a BS but also from a nearby RRU. The goal of this research is to look into how power is allocated to these two types of users, as well as how to split power for ID and EH at these users. OMA is used to serve multiple users at the same time, and equal bandwidth is allocated to each user. Let and denote the BS and RRU transmit powers, respectively. Then, the total transmit power within the cell region is calculated as . The single-selection transmission scheme is used for the DAS in this study; in other words, only the RRU closest to the cell-edge user sends the cooperation signal. Herein, the cell-edge user is referred to as a user 1 or weak user, whereas the cell-center user is referred to as a user 2 or strong user.
2.2. Channel Model
In this study, we use the Rayleigh fading channel. The channel coefficient, , from the RRU i in region k to user j in the target region is expressed as , where and are slow and fast fading components, respectively, for , , and . is the path-loss exponent, and is the distance between user j of the target region and RRU i of cell k. For simplicity, we drop the index which indicates the target region, that is, . The BS uses the channel gain to determine user 1 (weak user) and user 2 (strong user) that satisfy .
2.3. DAS-Assisted Simultaneous Wireless Information and Power Transfer System
Let
be the symbol sent from RRU
i to user
j for
and
. The received signals of user 1 and user 2 can be derived as:
for
and
. Here,
and
indicate the inter-cell interference and the noise at user
j, respectively. The interference to user
j from
K interfering cells can be derived as:
where
and
represent the allocated powers of the BS in cell
k, and
represents the symbol sent from region
k. Because the single-selection scheme is adopted, in the interfering cell
k, only one RRU shares the frequency band with user 1 and user 2 of the target region, and the RRU’s index is denoted as
in cell
k. In addition, we suppose that
and
for all
and
.
Each user employs a PS receiver whose PS ratio is
. Here,
and
of the received signal are used for ID and EH, respectively. Because the portion that is in phase with the received signal is not very large, we assume that the noise and interference cannot be harvested. Furthermore, according to [
18], the power of the interference and noise are very difficult to be large enough to pass through the diode of rectifier.
The signal-to-interference-plus-noise ratio (SINR) of users 1 and 2 can be obtained by:
where
is the variance of the inter-cell interference which can be obtained as:
Then, the data rates of users 1 and 2 are derived as:
where
,
,
, and
is the bandwidth. We assume that the half of the total bandwidth
B is allocated to both users as in [
19].
Even if the PS receiver decides to harvest power from the the received signal, it is impossible for the receiver to harvest without losses. The RF-to-DC power conversion function, according to: [
20,
21,
22], can be modeled by a non-linear function as follows:
where
is the input power of the RF-to-DC rectifier and
is the maximum harvested power of the saturated EH circuit. In addition,
a and
b are constants determined by EH circuit properties, such as the resistance, capacitance, and diode turn-on voltage, as defined in [
20,
21,
22]. Then, the energies harvested by user 1 and user 2 are represented as:
The important notations describing the system model are summarized in
Table 1.
3. Optimal Sum-Rate Maximization Method
In this section, the optimal sum-rate maximization methods for SWIPT-OMA in DAS are presented. Let
be the optimal power allocation of user 1, then
is the optimal allocated power for user 2. Then, the sum-rate maximization problem can be derived as follows:
where
is the sum-rate and
is the minimum rate constraint.
is the harvested energy requirement of user
j. From now on, we propose an optimal closed-form solution of (
14). According to (
9) and (
10), the inequality (
14a) re-expressed
and
as follows:
where
The inequalities (
14b) are also re-expressed as:
where
is the solution of
for
, which can be obtained as:
To avoid an outage event, and should be satisfied. If not, the outage event occurs and . Now, for a given , the optimal PS ratios are obtained in Lemma 1.
Lemma 1. For , if is positive, the optimal maximizing the sum-rate is obtained as .
Proof. For
, the following derivatives are derived.
According to (
21)–(
23),
is non-decreasing for both
and
. As a result, the optimal
and
are the highest values in the feasible region. □
From now on, we find the joint optimal solution through Lemma 2 and Theorem 1.
Lemma 2. Let and . For , is increasing for and decreasing for , where .
Proof. Differentiating
with respect to
gives:
Thus,
for
and
for
. □
By substituting the obtained
and
into problem (
14), the problem (
14) becomes a single-variable problem for
; thus, the optimal allocated power
can be found as in Theorem 1.
Theorem 1. The optimal that solves the problem (
14)
can be obtained as:wherewhen . If is satisfied, an outage event occurs. Proof. Let
and
be the solutions of
and
, respectively. Then,
and
can be obtained as (
26) and (
27), respectively. Because
is decreasing for
and
is increasing for
, the constraints (
15) and (
18) can be converted into
. If
, the feasible region of
is empty, so an outage event occurs.
From Lemma 2, is increasing for and decreasing for . Thus, if , is a decreasing function of in the feasible region; therefore, . If , belongs to the feasible region of , so the optimal . Otherwise, if , is an increasing function of in the feasible region; therefore, . □
4. Optimal Fairness-Rate Maximization Method
In this section, we describe the optimal power allocation and PS control scheme for the max-min fairness problem. The optimization problem can be expressed as:
where
. Because constraints (
28a)–(
28c) of the max-min fairness problem are same as constraints (
14a)–(
14c), the feasible region of
becomes
. The optimal PS ratios are obtained in the following lemma, which is similar to
Section 3.
Lemma 3. When , the optimal is derived as for .
Proof. As shown in Lemma 1, since and are non-decreasing for both and , is also non-decreasing for both and . As a result, the optimal is the highest value in the feasible region. □
From now on, we find the joint optimal solution through Lemma 4 and Theorem 2.
Lemma 4. For , is increasing for and decreasing for , where .
Proof. By solving the inequality with respect to , we have . Thus, when , . Because is an increasing function for , is also increasing for . Otherwise, when , becomes . Because is decreasing for , is also decreasing for . □
Theorem 2. The power allocation that solves the problem (
28)
can be obtained as:when . If is satisfied, an outage event occurs. Proof. From Lemma 4, is increasing for and decreasing for . Thus, if , is a decreasing function of in the feasible region; therefore, . If , belongs to the feasible region of , so the optimal . Otherwise, if , is an increasing function of in the feasible region; therefore, . □
5. Energy-Efficient Power Allocation and Power Splitting Design
In this section, we jointly optimizes power allocation and PS for the cell-center user and the cell-edge user to maximize the energy efficiency of SWIPT-OMA in the DAS. In this paper, the energy efficiency
is defined by the sum-rate over the total power consumption [
2,
16] as follows:
where
is the overall power consumption, which is defined as
.
and
are the consumed powers of the amplifiers in BS and RRU, respectively.
is the power consumed by circuit parts, which includes the power used to run the digital signal processors, mixers, and so on. Then, the energy-efficiency maximization problem is formulated as follows:
where
is the minimum rate constraint and
and
are the harvested power requirement for user 1 and user 2, respectively. Furthermore,
and
are the maximum power budgets of the BS and RRU, respectively. The energy efficiency
is not convex with respect to
,
,
,
, and
; therefore, we propose an iterative algorithm for solving the problem of (
31)–(
31c).
5.1. Outage Event Condition
Before finding the optimal solution, we discuss the outage event of the problem (
31) in this subsection. The first constraint (
31a) of the problem (31) can be divided into
and
, which can be re-expressed with respect to
as follows:
The minimum harvested energy constraints
and
can also be represented as follows:
In order for
that satisfies (
32)–(
35) to exist, the following conditions should be satisfied:
where
Figure 2 shows a graphical representation of the feasible region of
. The boundary of the feasible region consists of (
36)–(
39). If the region is an empty set, an outage event occurs. Therefore, we can state the following lemma.
Lemma 5. If , an outage event occurs.
Proof. From
Figure 2, if
or
, the feasible region becomes empty and an outage event occurs. Those conditions are the same as
. □
5.2. Optimal Power Allocation of Base Station
In this subsection, for a given , , and , we obtain the optimal power allocation for the BS, that is, and .
First, we propose the optimal
for a given
. According to (
31a) and (
31c), the feasible region of
can be expressed as follows:
where
and
. Because the second-order derivative
is negative,
is a concave function for
. We define
, which can be derived as follows:
Because
is negative,
is a decreasing function for
. Therefore, we can find the optimal
for given the other optimization variables by applying the bisection method to
.
The optimal
for a given
is derived in a similar way to obtaining the optimal
. The constraint of
can be expressed as follows:
where
and
. Here, the second-order derivative
is also negative, so
is a concave function for
. With the definition
,
can be derived as follows:
The optimal
can also be obtained by the bisection method.
The method for finding the optimal
,
using the bisection method is summarized in Algorithm 1.
Algorithm 1: Power allocation of the base station based on the bisection method |
Set , and as the stopping criteria Initialize and , satisfying ( 32)–( 35) while doend |
5.3. Optimal Power Allocation of Remote Radio Units
In this subsection, we demonstrate the optimal power allocation of the RRU (
) for a given
,
,
, and
. The constraint of
can be summarized as follows:
where
. The second-order derivative of
with respect to
is also negative; therefore,
is a decreasing function that can be derived as follows:
Similar to Algorithm 1, the optimal
can be obtained by applying the bisection method to
.
5.4. Optimal Power Splitting Ratio
In this subsection, strategies are introduced for obtaining the optimal and for a given , and . We sequentially obtain the optimal and . Then, an iterative algorithm is used to find the optimal .
For a given
, the constraint of
can be expressed as follows:
where
and
. With the definition
,
can be derived as follows:
Because
is negative,
is a decreasing function for
.
For a given
, the constraint of
can be expressed as follows:
where
and
. Similar to before,
can be obtained as follows:
where
is also decreasing for
; therefore, the algorithm that finds the optimal
follows the same process as in Algorithm 1.
5.5. Joint Power Allocation and PS Ratio Assignments
Now, we introduce the joint optimal solution for the original optimization problem (31). In our design, we obtain the optimal allocated powers and PS ratios successively, while assuming that the other variables are given, and apply an alternating algorithm to obtain the joint optimum solution. To summarize, the complete algorithm for maximizing the energy efficiency of the SWIPT-OMA in DAS is presented in Algorithm 2. The proposed scheme does not guarantee a global optimum owing to the non-convexity of (31), but finds a local maximum. Thus,
random initial points are set.
Algorithm 2: Joint optimal power allocation and power splitting ratio assignment |
for do end Select , , , and which generate the maximum |
6. Numerical Results
This section provides numerical results to verify the DAS using SWIPT and OMA and the proposed jointly optimized power allocation and PS ratio assignments. For the simulations,
interfering cells and
RRUs in each cell are considered. The maximum transmit power budgets are assumed to be
. The minimum rate constraints of the sum-rate and fairness-rate are assumed to be
bps/Hz and
bps/Hz, respectively, and the required harvested energies are set as
mW. The pathloss exponent, cell radius, bandwidth, noise variance, and circuit power are set as
,
m,
MHz,
dBm/Hz, and
W, respectively. In addition, we use
mW,
, and
as in [
21]. The cell-center user is generated uniformly within a circle with a radius of
, whereas the cell-edge user is generated within outer and inner circle radii of
R and
, respectively.
For the sum-rate and the fairness-rate comparison, we compare the performance with SWIPT-OMA without DAS scheme named ‘SWIPT-OMA’. Note that the BS and the RRU always consume their full power, that is, and , respectively.
In
Figure 3, the averaged optimal PS ratios for the max-sum-rate and the max-min-fairness problems are shown. As the power budget increases, the harvested energy requirement constraint can be satisfied by performing EH using only a small portion of the received signal; therefore,
is also increased. As shown in
Figure 3,
has a smaller value than
. Because channel environment of user 1 is weaker than user 2, user 1 utilizes a larger portion of the received power for EH than user 2. In addition, the
values of ‘SWIPT-OMA-DAS’ is much larger than that of ‘SWIPT-OMA’, which implies that ’SWIPT-OMA’ has worse EH efficiency and requires more power for EH to satisfy the constraints.
In
Figure 4, the sum-rate and fairness-rate performances of the proposed technique and ‘SWIPT-OMA’ are compared. We assumed that the sum-rate and fairness-rate are zeros when the outage event occurs. As shown in
Figure 4, the proposed schemes show much higher data rates than ‘SWIPT-OMA’. This is because DAS can overcome the data rate and EH degradation for a user who is far from the BS.
To show the advantages of the proposed framework and the energy efficient joint power allocation and PS ratio assignments scheme, we compare the proposed technique with the following schemes:
SWIPT-OMA: DAS is not employed; therefore, users receive data only through the BS. For , the optimal solution is obtained. The total power budget of the BS is supposed as for fair comparison.
Rate-opt SWIPT-OMA in DAS: This scheme determines the power allocation and PS ratios to maximize the sum-rate rather than the energy efficiency, as in
Section 3.
Figure 5 presents a comparison of the average PS ratio for different power budgets. In all the methods, the PS ratio of user 2, which has a relatively large channel gain, has a larger value than that of user 1. A large
indicates that the PS receiver has decided to use a large portion of the received signal power for ID. Thus, ‘Rate-opt SWIPT-OMA in DAS’ selects a larger PS ratio than other schemes to boost the sum-rate without considering energy efficiency. The proposed scheme has a higher
than ‘SWIPT-OMA’, because the user of the proposed scheme is more likely to receive a signal with higher power due to cooperation from a nearby RRU. Furthermore, we can see that the
difference between users 1 and 2 is small in the proposed scheme, because the significant degradation of channel gain from the BS can be alleviated with the help of the RRU.
Figure 6 shows the multi-cluster impacts on data rates. For the multi-cluster scenario, we assume that the bandwidth and BS power budget are uniformly allocated to
L clusters as
and
, respectively. If the number of clusters increases, the bandwidth and BS power allocated to each cluster decrease, making it difficult to guarantee the minimum data rate and minimum harvested energy requirements, after which more outage events occur; therefore, the performance decreases.
Figure 7 and
Figure 8 show the energy efficiency performance of the proposed and comparative methods with respect to the maximum power budget and minimum harvest power constraints, respectively. The energy efficiency increases as the maximum power budget grows, except for ‘Rate-opt SWIPT-OMA in DAS’, which does not pursue high energy efficiency. Similarly, as the minimum harvest power constraint increases, the energy efficiency of every comparison method decreases. We see that the proposed method is much more energy efficient than conventional schemes, and we see the significant benefit of using the DAS for SWIPT-OMA in both figures. As shown in
Figure 5, ‘SWIPT-OMA’ has a smaller PS ratio
, because it uses a larger portion of the received signal power to satisfy the minimum harvested energy constraint, so the proposed method shows better performance than ‘SWIPT-OMA’. Furthermore, with a limited power budget and/or a tight energy harvest constraint, the relative performance gain of the proposed scheme over other comparison methods grows.