1. Introduction
The increase of smartphones, laptops, and other mobile devices as well as data-hungry applications, need huge demands for ubiquitous coverage and very high data rates in cellular networks. However, homogeneous networks cannot satisfy these requirements [
1]. Then, two-fold efforts have been spent to meet the stringent requirements. On one hand, researchers have proposed heterogeneous networks (HetNets) where different types of base stations (BSs), e.g., macro BSs (MBSs) and small BSs (SBSs) are deployed in a multi-tier hierarchical structure. In this structure, all BSs have seamless coverage and reuse frequencies to achieve higher data rate [
2,
3]. On the other hand, the so-called non-orthogonal multiple access (NOMA) has been investigated as a potential technique to further improve the throughput of network [
4,
5,
6,
7]. Different from conventional orthogonal multiple access (OMA), NOMA serves multiple users at the same time/frequency/codes resource by allocating different powers for them, and the superposition coded signal can be decoded at receivers by successive interference cancellation (SIC). Therefore, the combination of HetNets and NOMA will exhibit great potential to satisfy the 1000-times increase of mobile broadband data for the upcoming fifth generation (5G) communication systems and beyond [
3].
However, the severe inter-tier and intra-tier interference make the NOMA-enabled HetNets challenging to achieve. Resource management plays an important role to alleviate these interference [
8]. For downlink communication, specifically, some work focuses on the sum rate maximization and shows higher spectral efficiency (SE) can be achieved by NOMA when considering the intercell interference [
9,
10,
11,
12]. Besides SE, energy efficiency (EE) is also a key performance metric investigated for resource allocation in NOMA-enabled HetNets [
8,
13,
14]. Moreover, EE is more important in uplink than in downlink NOMA-enabled HetNets since the devices in uplink communications are often battery-limited. It is a fact that the battery capacity has been improved at a very slow pace over the past decades [
15], and hence this increase cannot scale with the high energy consumption caused by the increasing traffic demands. Meanwhile, EE has emerged as a new prominent performance metric for wireless communication networks designs due to the economic, operational, and environmental concerns [
16,
17]. Therefore, it is a stringent work to improve EE for uplink transmission.
Machine-to-machine (M2M) communications, also known as machine-type communications (MTC), enable pervasive connections to support IoT. M2M communications are one of the potential applications of NOMA-enabled HetNets [
18], since NOMA-enabled HetNets provide a practical infrastructure to offer massive access opportunities for such a huge number of devices, especially for the cases in which each device only needs to send a small amount of data periodically in uplink. One of the challenges for HetNets with M2M communications is the access control, which can manage the engagement of massive MTC devices (MTCDs) to the core network. Among the existing access solutions, deploying MTC gateways (MTCGs) is an effective approach to connect M2M communication and cellular communication [
19,
20,
21]. When mobile user (MU) has more power and storage space than MTCDs (e.g., smart sensors), the MU can be configured as the MTCG, as proposed in [
22].
Since 5G will be HetNets including various network models (e.g., cellular networks, wireless networks (WSNs), and low power wireless area networks) to support high data rate and massive devices [
3], our work combing M2M communication and cellular network has a large significance for this heterogeneous scenario. The short distance communication in our system model can be realized with WSNs, which provide a new way to help the sink nodes in WSNs communicate to the core network. For example, the MTCDs can be the sensors in an environmental monitoring WSN, and they can transmit the collected data to the core network through a mobile device in cellular networks with NOMA. Therefore, our work also has a practical significance for sensors work.
Recently, there have been some studies addressing the aforementioned challenges of applying NOMA in HetNets for EE maximization. In [
23], a distributed user association algorithm based on inter-cell interference plus noise ratios of BS and a centralized user association based on the popular size of BS were both proposed. After user association was determined, a power control algorithm was proposed based on Lagrangian dual method, then a one-dimensional search algorithm was used to search Lagrangian multiplier, which added algorithm complexity. Two specific examples were provided to demonstrate the effectiveness of unified NOMA-enabled heterogeneous ultra-dense networks with user association and power control in [
18]. An alternated energy efficient resource allocation algorithm based on fixed power allocation was first proposed in [
13]. Then, two iterative energy-efficient resource allocation algorithms were proposed to update for better EE based on Lagrangian dual method. Joint base station association and power control optimization algorithms were proposed based on coalition formation games and interior-point method in [
24], but sub-channel allocation and fractional equation for EE maximization had not been considered. Moreover, the user association algorithms in the aforementioned work were all considered with fixed power allocation firstly, whereafter iterative algorithms were used to obtain the final optimal value.
There are also some studies on the usage of M2M communications in NOMA systems. For example, energy-efficient resource allocation with hybrid division multiple-access NOMA for cellular-enabled M2M communications was researched in [
25,
26]. With MTCDs cluster formation known beforehand, standard convex optimization and Lagrange duality methods were employed respectively for power control in [
25,
26]. User clustering in NOMA-aided cellular M2M communication systems was researched in [
27,
28] with millimeter-wave and narrow-band IoT separately. A joint power and sub-channel allocation for secrecy capacity algorithm was proposed in [
29] to obtain the suboptimal solution of the optimization problem. However, the aforementioned work deploys M2M user in single-cell networks. The trend of more and more intensive network deployment motivates us to deploy M2M-enabled NOMA in the scenarios with multi-tier HetNets and new resource allocation needs to be considered with the non-convexity caused by inter-cell interference in HetNets.
In this paper, we focus on the uplink EE maximization via user association and power control for M2M-enabled HetNets using NOMA. In this scenario, one macro base station (MBS) is located in the cell center. Each small cell has one small base station (SBS) located in the cell center. MUs are distributed randomly in the cell. An MU acting as an MTCG can decode and forward both the information of MTCDs and its own data to the BS. The EE (bits/Joule) maximization problem is formulated and solved to obtain the optimal MU association and power allocation. The main contributions of this paper are summarized as follows:
We propose a new framework of M2M-enabled HetNets with NOMA. In this framework, control data separation architecture, i.e., control information and data message are separated, which can reduce the signal overhead [
30]. NOMA is adopted by the MTCDs to transmit the information to MUs which is regarded as the relay. MUs decode the overlaid information and simultaneously transmit received data to the BSs based on the NOMA principle.
In order to solve the EE maximization optimization problem, a BS and a sub-channel are included in a couple, since a MU can only associate one BS at one sub-channel. Then, a many-to-one MU association algorithm is proposed based on matching game [
31]. Through swap operation among each couple, the EE maximization problem can be tackled by solving the power control problem at each sub-channel. Compared with the previous studies on the algorithms (user association and power allocation) [
13,
18,
23,
24], our algorithms are jointly optimized and fixed power allocation is not required for initialization.
Two power control algorithms are proposed based on sequential optimization [
32,
33]. The fractional programming [
34] and sequential optimization are combined to develop a novel sequential fractional power control algorithm (SFPCA), from which the original problem is transformed to be convex and requires less computational complexity. The other algorithm combines the exhaustive search method with sequential optimization, which can verify the correctness of SFPCA.
The rest of this paper is organized as follows. The system model and problem formulation are focused in
Section 2. The MU association matching algorithm is proposed in
Section 3. The power control problem is solved in
Section 4. Numerical results are provided in
Section 5, and concluding remarks are given in
Section 6.
Notations: Lowercase and uppercase boldface letters denote vectors and matrices, respectively. We use uppercase decorated letters to denote sets. For an arbitrary set , we always have the corresponding uppercase M to the denote the cardinality of , i.e., , denotes the transpose operator.
3. MU Association
The MTCDs associated to the corresponding MU are known beforehand. Since solving the optimization problem is equal to obtain the optimal EE of each sub-channel, the MU association will become the matching problem among BSs, sub-channels and MUs to achieve sub-channel EE maximization. Thus, we propose a MU association algorithm using matching game [
30] in the following parts.
3.1. Matching Problem Formulation
To develop a low-complexity MU association algorithm, we first regard a sub-channel and a BS as a couple, denoted as (). Then, the optimization problem is transformed to match the MUs to the couples and allocate power appropriately, such that the EE can be maximized. Finally, the matching problem is a many-to-one problem between MUs and couples based on matching game, which is described as follows.
Definition 1. Given two disjoint sets,denotes the set for MUs, andrepresents the couples. A many-to-one matching Ψ is a mapping from the setinto the set of all subsets offor, , satisfying
- i)
;
- ii)
;
- iii)
;
- iv)
.
Condition i indicates that each MU matches with a sub-channel-BS couple. On the other hand, each couple matches a subset of MUs, which is illustrated in condition ii. Condition iii states a MU can only associate one BS and choose one sub-channel while each couple matches MUs.
The aim of each couple is to maximize its own EE. To this end, we exploit the swap operation into our matching algorithm. A swap operation means two MUs matching with different couples exchange their matchings based on different cases, while other MUs remain their matchings. The EE of the exchanged couples will be recomputed by the power control algorithm. Note that how to allocate power to obtain the optimal EE for a given sub-channel will be presented in the next section, and we assume it is known in advance. A swap operation will be approved and the matching will be exchanged only when all EE of the sub-channels belonging to the exchanged couples increase if the swap is performed. The swap operation will be continued until no swap is further preferred. More details are described in Algorithm 1.
Algorithm 1 The MU association matching algorithm. |
Initialization phase:, |
1: for do |
2: , |
3: while () do |
4: . Assign to the couple , , , and set |
5: end while |
6: end for Swap matching phase: |
7: while () do |
8: for do |
9: for do |
10: if then |
11: continue; |
12: else if and are both in the coverage of the BSs of each other then |
13: switch () |
14: case and belonging to the same BS and different sub-channels: |
15: Calculate and compare the EE of the two sub-channels before and after the swap using the power control algorithm. If the EE of the two-subchannels both improve, exchange the sub-channel, form the new couple, and set . |
16: case and belonging to the different BSs and different sub-channels: |
17: Calculate and compare the EE of the two sub-channels before and after the swap using the power control algorithm. If the EE of the two sub-channels both improve, exchange the couple, form the new couple, and set . |
18: case and belonging to the different BSs and same sub-channels: |
19: Calculate the EE of the sub-channel before and after the swap using the power control algorithm. If the EE of the sub-channel has been improved, exchange the BS, form the new couple, and set . |
20: end switch |
21: end if |
22: end for |
23: end for |
24: end while |
3.2. Matching Algorithm
Algorithm 1 contains a initialization phase and a swap matching phase. Considering the user fairness, the number of MUs accommodated by one sub-channel in a given BS is at most . In the initialization step, the basic idea is to associate the MU to the couple providing the largest channel gain. This will lead to either a higher data rate for the MU, or a lower transmit power. Since the value of sub-channel gain between MU and the uncovered BSs is invalid and the maximized sub-channel gain is always chosen, there is no need to know whether the MUs are in the coverage of the exchanged BSs. However, in the swap matching phase, this judgment should be considered at first to avoid the invalid swap. Then, exchange will happen in the three cases. Iterations will continue until no swap operation can be approved in a new round.
3.3. Convergence and Complexity
Theorem 1. The proposed MU association and power control algorithm converges after a finite number of swap operations.
Proof of Theorem 1. For each swap operation, the matching changes from to . We have and to denote the corresponding EE of of and on . Based on the aim of swap operation, we have , that is, the EE of each sub-channel increases after each swap matching. Since each sub-channel is orthogonal to each other, the system EE will increase owing to the improved EE of each sub-channels. Moreover, the system EE has an upper bound due to the limited transmit power of each MU. Therefore, the MU association algorithm and power allocation converge after a finite number of swaps. □
4. Power Control
In this section, we will investigate the optimal power control design appearing in Algorithm 1 to obtain the maximum EE of
. Before we present the optimization problem for
maximization, we first deal with
, which can be rewritten as
and we can also obtain
Due to the multi-interference in the sum-rate function in (8),
in (6) is non-convex and cannot be directly solved by the generalized fractional programming approach. Then, we first transform the numerator into the difference of two non-negative functions and the
maximization can be rewritten as
with
where
denotes the transmit power vector for MUs on
. Moreover, we have
Note that
are concave functions regarding to
, then the numerator of (10a) and the constraint functions in (10b) are expressed as the difference of concave functions, which are not concave in general. Motivated by [
31,
32], where sequential optimization is used to solve the similar problem as (10), we adopt this method and combine it with fractional programming and exhaustive search to propose two power control algorithms. Before introducing the two algorithms, we first present the details of the sequential optimization theory in the next sub-section.
4.1. Sequential Optimization Theory
Sequential optimization is a powerful tool that can tackle a difficult optimization problem by solving a sequence of approximate problems in simple forms with affordable complexity. Specifically, we give a formal maximization problem
with a compact feasible set as [
31,
32], shown as
where
is the differentiable objective with constraints
. Let
be the problem solved in the
v-th iteration by the sequential method to tackle problem
, which can be written as
where
is the differentiable objective with the constraints
. Then, if
and
are suitable continuous functions and constraints, they must satisfy the following two properties:
- 1)
, ;
- 2)
,
.
is the optimal solution of the problem solved at iteration -th. This means the solution sequence of (14) monotonically increases the value of (13), i.e., for all v, which guarantees the convergence of the sequential method. Next, if the following third property is also satisfied:
- 3)
, .
then every limit point of of (14) fulfills the Karush–Kuhn–Tucker (KKT) conditions of problem in (13).
Therefore, if a maximization problem finds suitable approximate problems which can fulfill the above three properties, its optimal value can be approximated by solving the monotonically increased sequential problems. The critical issue is that the suitable approximate problems are solved easier than the original problem. In the rest of this section, we will first find the sequential approximate problems to the numerator in problem (10).
4.2. Sequential Fractional Power Control Algorithm
Based on sequential optimization, we should find a sequence problem to approximate optimization problem (10). To circumvent this issue, we obtain the following main result with the first-order Taylor expansion at of .
Proposition 1. For any given,
the sequence approximation problem of (10), denoted bycan be written aswith optimal solution,
where If, thenis monotonically increasing and converges to a value. Furthermore, any limit point of sequencethat achievesfulfills the KKT optimality conditions of (10a).
Proof of Proposition 1. As we know, any concave function is the upper-bounded of its first-order Taylor expansion at any point. Since
and
are concave functions, for any power vector
we have
Hence, (15a) and (15b) are lower bounds of (10a) and (11), respectively. Since the lower bounds in (16) are tight when evaluated by , it follows that (15a) and (15b) are equal to (10a) and (11), respectively, for . Similarly, the gradients of (15a) and (15b) are equal to those of (10a) and (11), for . Thus, (15) fulfills all the properties described in the above sub-section, which completes the proof of this proposition. □
For any
, problem (15) has a concave numerator and an affine denominator, while the constraint functions in (15b) and (15c) are both concave and affine. Therefore, (15) is a single-ratio problem, which can be solved by the generalized fractional programming. We adopt the widely used Dinkelbach’s algorithm to solve it. According to Dinkelbach’s method [
33], we first introduce the following auxiliary function
with
, and
.
Theorem 2. Let,
anddenote the optimal value, optimal solution and its elements of problem (15), respectively. Then, we haveif and only if Proof of Theorem 2. Theorem 2 was proved in [
33,
36], and we omit it due to the limited space. □
The optimal
can be obtained by Dinkelbach’s method, which is summarized in Algorithm 2. As shown in the algorithm, we need to solve the problem (23) for a given parameter
in each iteration. In Algorithm 2,
has been updated as
in each iteration until convergence. or reaching the maximum number of iterations.
denotes the optimal power of the following problem in the
c-th iteration, which can be obtained in Algorithm 3, as given by
Algorithm 2 The Dinkelbach’s algorithm. |
Initialization phase: |
Set iteration , , the maximum number of iterations , and error tolerance . |
1: repeat |
2: Solve the equivalent problem (23) for a given to obtain the solution . |
3: , |
4: . |
5: until or |
6: , . |
Algorithm 3 The algorithm for solving problem (23). |
Initialization phase: |
Set , iteration index , the maximum iterations and error tolerance . Calculate . |
1: repeat |
2: Solve the problem (23) to obtain the optimal solution for given and . |
3: . |
4: Set and cacluate . |
5: until |
6: . |
4.3. Computational Complexity Analysis
In above sub-section, we have proposed the SFPCA including two steps, i.e., Algorithms 2 and 3. The computational complexity of them are separately discussed. First, we use
C to denote the number of iterations for Algorithm 2, where
C is bounded by
. From Section V, we can see that Algorithm 2 will converge after a few number of iterations. Then we discuss the computational complexity of Algorithm 3, the complexity of this algorithm is mainly caused by (23), and denoted by X. The computational complexity of (23) is
[
37], where
is the number of MUs at
. The complexity of Algorithm 3 is
, where
V is the the number of iterations bounded by
. In summary, the computational complexity of the power control algorithm is
.
4.4. Sequential Exhaustive Algorithm
To evaluate the performance of the SFPCA, a sequential exhaustive algorithm (SEA) combined with sequential optimization and exhaustive search is proposed in this section. The detailed procedures of the compared algorithm is illustrated as follows. To solve the problem in an easier manner, we introduce the auxiliary variable
, as given by
if we fix
, the objective function (15a) can be recast as
Due to the multi-user interference, we cannot solve problem (25) by standard convex optimization tools. Similar to SFPCA, sequential optimization is applied and the approximate problem can be shown as
It can be observed that since
is fixed, (26) is equivalent to minimize the linear function
in the denominator, subject to convex constraints. Then, problem (26) can be solved by plain convex programming. To implement an efficient line search for
, the bound of
is given by
Then, the optimial can be obtained by searching an appropriate value of with stepsize .
5. Numerical Results and Discussions
In this section, the effectiveness of our proposed MU association and power control algorithms in M2M-enabled HetNets with NOMA was demonstrated by Monte Carlo simulations. The HetNets included one MBS and two SBSs, and the radius of the cells for them were 200 m and 80 m, respectively. MUs were randomly and uniformly distributed. The values of the simulation parameters are summarized in
Table 1.
We considered the EE performance obtained from EE maximization and sum-rate maximization with different
in
Figure 2. The latter could be obtained in the first iteration of Algorithm 2 due to
. In order to reflect the influence of
, we gave four schemes of different
. From the figure, we can see that all of the four schemes had a “green point”, where EE and sum rate could both achieve their optimal values. Different from sum rate, EE became gradually flat while sum rate decreased after “green point” as
grew. The reason is when the maximum EE is achieved, no more transmit power is needed. For sum rate maximization, larger sum rate requires more transmit power, and its ratio (EE) may decrease, since the numerator (sum rate) and denominator (sum transmit power) both grow. We can also see that the EE decreased as
increased, because the increase of
means the data rate requirement of MUs increases. It is worth noting that even though MTCDs have lower data rate and transmit power, they can also have a strong influence on the overall uplink EE with their massive number. Furthermore, Algorithm 1 with higher EE had similar tendency as Algorithm 4, which proves the correctness of our algorithms.
Figure 3 shows the EE performance with respect to different data rate requirements of MUs with the different transmit power of MTCDs. The four curves all decreased as
increased. This is due to the fact that higher data rate will narrow the feasible value regions of the transmit power. Note that the four curves decreased slightly first, when
, the EE of the four schemes all declined distinctly, since higher data rate requirement may require more transmit power, destroying the balance of sum transmit power and sum rate. As explained above, the increase of
leads to the increase of data rate requirement, and the variation of EE is in line with the reason as the figure shown.
Algorithm 4 Sequential exhaustive algorithm. |
Initialization phase: |
, |
1: for do |
2:
|
3: end for |
4: Obtain the optimal solution of . |
Figure 4 and
Figure 5 shows the convergence property of Algorithms 2 and 3. For simplicity, the numerical results in two figures are from a random chosen sub-channel, where
. In
Figure 4, we can see that the number of iterations are limited within four times. To show the influence of
, we give the different values of
in
Figure 5, where
. It is shown that the initial values have an effect on the number of iterations. Specifically, when
, less than 11 times is needed to reach the convergence. Although the initial values affect the number of iterations, it does not affect the final results.
To show the relationship between the different numbers of MUs and MTCDs and the EE, we have
Figure 6. It is not surprise to see that the EE performance of all these schemes increases as
grows. From the four schemes, we can find out that the EE of
is much larger than that of
, since the NOMA scheme can obtain much higher EE by supporting multiple MUs, and they can choose the suitable couples by swap operations for better EE. From the Algorithm 1, we know that power control algorithm needs to be executed after each swap operation, that is, the number of iterations increases with the increase of
K and more process time are required. From the
Figure 6, we can also see that the EE of
is larger than that of
under the same
K, i.e.,
K has much greater impact on EE than
, since the NOMA scheme can obtain much higher EE by supporting multiple MUs and the increasing
represents the increase data rate requirement of MUs.
Figure 7 presents the cumulative distribution function (CDF) of the number of swap operations of different scenarios when the matching algorithms reached convergence. From the figure we can see that more swap operations were needed for a larger number of MUs and sub-channels, such as,
needed more swap operations than that
and
. Especially, less than 70 swap operations were needed for
and
.