Next Article in Journal
Robust Cylindrical Panorama Stitching for Low-Texture Scenes Based on Image Alignment Using Deep Learning and Iterative Optimization
Next Article in Special Issue
A Distributed Testbed for 5G Scenarios: An Experimental Study
Previous Article in Journal
Embodied Emotion Recognition Based on Life-Logging
Previous Article in Special Issue
Fog Computing Enabling Industrial Internet of Things: State-of-the-Art and Research Challenges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

User Association and Power Control for Energy Efficiency Maximization in M2M-Enabled Uplink Heterogeneous Networks with NOMA

1
Key Laboratory of Universal Wireless Communications, Ministry of Education, Beijing University of Posts and Telecommunications, Beijing 100876, China
2
Wuxi BUPT Sensory Technology and Industry Institute CO. LTD, Wuxi 214000, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(23), 5307; https://doi.org/10.3390/s19235307
Submission received: 8 October 2019 / Revised: 14 November 2019 / Accepted: 25 November 2019 / Published: 2 December 2019
(This article belongs to the Special Issue Internet of Things and Sensors Networks in 5G Wireless Communications)

Abstract

:
To support a vast number of devices with less energy consumption, we propose a new user association and power control scheme for machine to machine enabled heterogeneous networks with non-orthogonal multiple access (NOMA), where a mobile user (MU) acting as a machine-type communication gateway can decode and forward both the information of machine-type communication devices and its own data to the base station (BS) directly. MU association and power control are jointly considered in the formulated as optimization problem for energy efficiency (EE) maximization under the constraints of minimum data rate requirements of MUs. A many-to-one MU association matching algorithm is firstly proposed based on the theory of matching game. By taking swap matching operations among MUs, BSs, and sub-channels, the original problem can be solved by dealing with the EE maximization for each sub-channel. Then, two power control algorithms are proposed, where the tools of sequential optimization, fractional programming, and exhaustive search have been employed. Simulation results are provided to demonstrate the optimality properties of our algorithms under different parameter settings.

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 M , we always have the corresponding uppercase M to the denote the cardinality of M , i.e., M = M , [ · ] T denotes the transpose operator.

2. System Model and Problem Formation

2.1. System Model

As shown in Figure 1, we consider uplink HetNets with M2M communications, where all MUs are anchored to the control base station (CBS). The CBS performs the MU association algorithm to select the best serving BS for MUs and establishes a high BS–MU connection through backhaul links. Each MTCD selects the nearest MU as an MTCG. Since NOMA is adopted between MTCDs that select the same MU as their MTCGs, SIC is performed at the MU to gather the interference and channel gain will be obtained by channel estimation at the MU. The HetNets consist of a set F = 0 , , F of BSs and a set K = 1 , , K of MUs. Each MU is regarded as an MTCG, which can acts as a relay for some MTCDs. Denote U k as the specific set of MTCDs served by MU k ( MU k ). The index 0 denotes the MBS and other indexes stand for the SBSs in set F . Without special explanation, we always have f F , F = F . The system bandwidth shared by all BSs is divided into N orthogonal sub-channels, and each one is assigned with bandwidth B. For convenience, hereinafter we always have n N = { 1 , 2 , N } to denote the sub-channel. MUs are served by BSs according to the BSs’ coverage.

2.2. NOMA Strategy

Multiple MTCDs can simultaneously transmit signals to the MU using NOMA. Since MUs and MTCDs use different transmission modes, we ignore the interference between MUs and MTCDs. The interference between MTCDs in different BSs is also not considered. According to the NOMA principle, the received signal of MU k is
Y k = j U k h j k k q j k s j k + n k ,
where h j k is the channel between MTCD j k and MU k ; q j k and s j k denote the transmit power and message of MTCD j k ; and n k represents the additive zero-mean Gaussian noise with variance σ 2 . U k represents the set of MTCDs which are served by MU k . Without loss of generality, the channels are sorted by h 1 k k 2 h 1 k k 2 σ k 2 σ k 2 > h 2 k k 2 h 2 k k 2 σ k 2 σ k 2 > h U k k 2 h U k k 2 σ k 2 σ k 2 > 0 . Applying SIC in NOMA [33], the achievable data throughput for MTCD j k at MU k is given by
R j k = log 2 1 + H j k k q j k 1 + I j k ,
where H j k k = h j k k 2 h j k k 2 σ k 2 σ k 2 , I j k = i U k H i k k < H j k k q i k k H i k k , and we define I j k = 0 for j k = U k . After MUs successfully decode the messages from MTCDs, all MUs simultaneously transmit data to the BS based on the NOMA principle. Denote h k f n = g k f n d k f α as the channel gain between MU k and BS f at sub-channel n ( SC n ). g k f n denotes the corresponding Rayleigh fading channel gain; α is the path loss factor; and d k f is the distance between MU k and BS f. In order to split the superimposed signals on SC n in BS f, SIC is carried out at BS f. Based on the uplink NOMA protocol [35], the signal of MU with the highest channel gain will be first decoded at BS f and experiences interference from other MUs having relatively weaker channel gains on SC n . Therefore, the channel gains of MUs over SC n in BS f are sorted as h 1 f n 2 h 1 f n 2 σ f n 2 σ f n 2 > h 2 f n 2 h 2 f n 2 σ f n 2 σ f n 2 > h k f n 2 h k f n 2 σ f n 2 σ f n 2 > h S f n f n 2 h S f n f n 2 σ f n 2 σ f n 2 , where S f n = S f n . Then the transmit data rate of MU k associated with BS f over SC n can be expressed as
R k f n = B log 2 1 + p k f n H k f n 1 + I k f n + ϕ k f n ,
where H k f n = h k f n 2 h k f n 2 σ f n 2 σ f n 2 . I k f n is the interference that MU k receives from other MUs whose channel gains are smaller than that on SC n of BS f, which can be given by
I k f n = i S f n H i f n < H k f n p i f n H i f n .
ϕ k f n = f F \ f i S f n p i f n H i f n is the interference from MUs associated to other BSs on SC n . Then the data rate of MUs at SC n is
R n = B f F k S f n log 2 1 + p k f n H k f n 1 + I k f n + ϕ k f n .

2.3. Problem Formation

In this paper, we focus on the EE maximization problem for all MUs considering the minimum data rate requirements of them. The MU association contains two parts: BS selection and sub-channel allocation. For a given MU association, the out-of-cell interference only come from the MUs associated with different BSs at the same sub-channel due to the orthogonality among the sub-channels. Then, each MU may not concern the whole EE, but the sub-channel EE it chooses. Therefore, the optimization problem is converted into solving the EE maximization of each sub-channel by appropriate MU association including BS selection and sub-channel allocation and power control.
From a physical standpoint, the efficiency with which a system uses a given resource is the ratio between the benefit obtained by using the resource and the corresponding incurred cost [17]. Applying this general definition to the uplink communication at SC n , then EE of SC n can be written as
E E n = R n f F k S f n p k f n + P c ,
where P c is the additional circuit power consumption over each sub-channel. Then the considered EE optimization problem can be formulated as
max P n N E E n
s . t . R k f n j U j k R j k + R r e q , k K , n N ,
p k f n > 0 , p k f n P m a x , k K , f F , n N ,
where P is the transmit power vector with elements p k f n ; P m a x is the maximum transmit power of each MU; and R r e q is the minimum data rate requirement of a MU. Since each MU is regarded as an MTCG for MTCDs, they should ensure the data rate that MTCDs can be uploaded to the SBS, therefore we have constraint (7b) as the data rate requirement of MU k associated to BS f at SC n [22]. Constraint (7c) is used to guarantee the feasible value ranges of P .

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 ( n , f ). 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, K = 1 , , K denotes the set for MUs, and M = ( 1 , 1 ) ( 1 , 2 ) , ( 2 , 1 ) , ( n , f ) , ( N , F ) represents the couples. A many-to-one matching Ψ is a mapping from the set K M into the set of all subsets of K M for f F , k K , n N satisfying
i) 
Ψ ( k ) M ;
ii) 
Ψ ( n , f ) K ;
iii) 
Ψ ( k ) = 1 , Ψ ( S C n , f ) = S f n ;
iv) 
( n , f ) = Ψ ( k ) k Ψ ( n , f ) .
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 S f n 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: L = K N F , K ^ = K
1: for l = 1 : L do
2:   M ^ = M , C o u n t = 1
3:  while ( C o u n t M ) do
4:    H k * f * n * = a r g m a x H k f n k K ^ , ( f , n ) M ^ . Assign k * to the couple ( f * , n * ) , K ^ = K ^ \ k * , M ^ = M ^ \ ( f * , n * ) , and set C o u n t = C o u n t + 1
5:  end while
6: end for
  Swap matching phase: I n d i c a t o r = 1
7: while ( I n d i c a t o r ) do
8:  for u = 1 : K do
9:   for k = 1 : K do
10:    if Ψ ( k ) = Ψ ( u ) then
11:     continue;
12:    else if MU k and MU u are both in the coverage of the BSs of each other then
13:     switch ( Ψ ( k ) , Ψ ( u ) )
14:     case MU k and MU u 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 I n d i c a t o r = 1 .
16:     case MU k and MU u 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 I n d i c a t o r = 1 .
18:     case MU k and MU u 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 I n d i c a t o r = 1 .
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 K F N . 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 Ψ e x to Ψ n o w . We have E E n , e x and E E n , n o w to denote the corresponding EE of of Ψ e x and Ψ n o w on SC n . Based on the aim of swap operation, we have E E n , n o w > E E n , e x , 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 S C n . Before we present the optimization problem for E E n maximization, we first deal with R n , which can be rewritten as
R n ( p f n ) = B f F log 2 1 + k S f n p k f n H k f n 1 + ϕ k f n ,
and we can also obtain
R t o t , k = j U k R j k = log 2 1 + j U k H j k k q j k .
Due to the multi-interference in the sum-rate function in (8), E E n 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 E E n maximization can be rewritten as
max p f n η ˜ n = F + ( p f n ) F ( p f n ) f F k S f n p k f n + P c
s . t . ( 7 c ) , ( 11 ) ,
with
C k f n + ( p f n ) C k f n ( p f n ) 0 , k K , f F , n N ,
where p f n = [ p 1 f n , p 2 f n , , p k f n , , p S f n F n ] T denotes the transmit power vector for MUs on S C n . Moreover, we have
F + ( p f n ) = B f F log 2 1 + k S f n p k f n H k f n + ϕ k f n , F ( p f n ) = B f F log 2 1 + ϕ k f n , C k f n + ( p f n ) = B log 2 1 + p k f n H k f n + I k f n + ϕ k f n R t o t , k R r e q , C k f n ( p f n ) = B log 2 1 + I k f n + ϕ k f n .
Note that F + , F , C + , and C are concave functions regarding to p f n , 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 F ¯ with a compact feasible set as [31,32], shown as
max x f 0 ( x )
s . t . f i ( x ) 0 , i { 1 , , I } ,
where f 0 ( x ) is the differentiable objective with constraints f i ( x ) 0 . Let G ( v ) be the problem solved in the v-th iteration by the sequential method to tackle problem F ¯ , which can be written as
max x g 0 ( v ) ( x )
s . t . g i ( v ) ( x ) 0 , i { 1 , , I } ,
where g 0 ( v ) ( x ) is the differentiable objective with the constraints g i ( v ) ( x ) . Then, if g 0 ( v ) ( x ) and g i ( v ) ( x ) are suitable continuous functions and constraints, they must satisfy the following two properties:
1)
g 0 ( v ) ( x ) f 0 ( x ) , g i ( v ) ( x ) f i ( x ) x ;
2)
g 0 ( v ) ( ( x * ) ( v 1 ) ) = f 0 ( ( x * ) ( v 1 ) ) , g i ( v ) ( ( x * ) ( v 1 ) )
f i ( ( x * ) ( v 1 ) ) .
( x * ) ( v 1 ) is the optimal solution of the problem solved at iteration ( v 1 ) -th. This means the solution sequence x * ( v ) of (14) monotonically increases the value of (13), i.e., f 0 ( ( x * ) ( v ) ) f 0 ( ( x * ) ( v 1 ) ) for all v, which guarantees the convergence of the sequential method. Next, if the following third property is also satisfied:
3)
g 0 ( v ) ( ( x * ) ( v 1 ) ) = f 0 ( ( x * ) ( v 1 ) ) , g i ( v ) ( ( x * ) ( v 1 ) ) = f i ( ( x * ) ( v 1 ) ) .
then every limit point of x ( v ) of (14) fulfills the Karush–Kuhn–Tucker (KKT) conditions of problem F ¯ 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 p f n ( v ) of F ( p f n ) .
Proposition 1.
For any given p f n ( v ) , the sequence approximation problem of (10), denoted by G ( v ) can be written as
max η n = F + ( p f n ) F ˜ ( p f n ) f F k S f n p k f n + P c
s . t . C k f n + ( p f n ) C ˜ ( p f n ) 0 ,
7 c ,
with optimal solution p f n * ( v ) , where
F ˜ ( p f n ) = F p f n ( v ) F p f n ( v ) T p f n p f n ( v )
C ˜ ( p f n ) = C k f n p f n ( v ) C p f n ( v ) T ( p f n p f n ( v ) )
If p f n ( v ) = p f n * ( v 1 ) , v 1 , then η n p f n * ( v ) ( v ) is monotonically increasing and converges to a value η ˜ n . Furthermore, any limit point of sequence η n p f n * ( v ) ( v ) that achieves η ˜ n fulfills 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 F ( p f n ) and C ( p f n ) are concave functions, for any power vector p f n ( v ) we have
F + ( p f n ) F ( p f n ) F + ( p f n ) F ˜ ( p f n ) = F + ( p f n ) F p f n ( v ) F p f n ( v ) T p f n p f n ( v ) ,
C + ( p f n ) C ( p f n ) C + ( p f n ) C ˜ ( p f n ) = C + ( p f n ) C k f n p f n ( v ) C p f n ( v ) T ( p f n p f n ( v ) ) .
Hence, (15a) and (15b) are lower bounds of (10a) and (11), respectively. Since the lower bounds in (16) are tight when evaluated by p f n ( v ) , it follows that (15a) and (15b) are equal to (10a) and (11), respectively, for p f n = p f n ( v ) . Similarly, the gradients of (15a) and (15b) are equal to those of (10a) and (11), for p f n = p f n ( v ) . Thus, (15) fulfills all the properties described in the above sub-section, which completes the proof of this proposition. □
For any p f n ( v ) , 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
T ( p f n , η n ) = f n ( p f n ) η n g n ( p f n ) ,
with f n ( p f n ) = F + ( p f n ) F ˜ ( p f n ) , and g n ( p f n ) = f F k S f n p k f n + P c .
Theorem 2.
Let η n * , p f n * and p k f n * denote the optimal value, optimal solution and its elements of problem (15), respectively. Then, we have
η n * = F + ( p f n * ) F ˜ ( p f n * ) f F k S f n p k f n * + P c = max F + ( p f n ) F ˜ ( p f n ) f F k S f n p k f n + P c ,
if and only if
max T ( p f n , η n * ) = f n ( p f n ) η n * g n ( p f n ) = f n ( p f n * ) η n * g n ( p f n * ) = 0 .
Proof of Theorem 2.
Theorem 2 was proved in [33,36], and we omit it due to the limited space. □
The optimal η n * 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 η n ( c ) in each iteration. In Algorithm 2, η n has been updated as η n ( c ) in each iteration until convergence. or reaching the maximum number of iterations. p f n ( c ) denotes the optimal power of the following problem in the c-th iteration, which can be obtained in Algorithm 3, as given by
max p f n T ( p f n , η n ( c ) ) = f n ( p f n ) η n ( c ) g n ( p f n ) s . t . C k f n + ( p f n ) C ˜ ( p f n ) 0 , 7 c ,
 Algorithm 2 The Dinkelbach’s algorithm.
Initialization phase:
 Set iteration c = 1 , η n ( c ) > 0 , the maximum number of iterations C max , and error tolerance τ > 0 .
1: repeat
2: Solve the equivalent problem (23) for a given η n ( c ) to obtain the solution p f n ( c ) .
3: η n ( c ) = F + ( p f n ( c ) ) F ˜ ( p f n ( c ) ) f F k S f n p k f n ( c ) + P c ,
4: c = c + 1 .
5: until T ( η n ( c 1 ) , p f n ( c 1 ) ) τ or c > C max
6: η n * = η n ( c 1 ) , p f n * = p f n ( c 1 ) .
 Algorithm 3 The algorithm for solving problem (23).
Initialization phase:
 Set p f n ( 0 ) , iteration index v = 0 , the maximum iterations V max and error tolerance μ . Calculate f n ( p f n ( 0 ) ) η n ( c ) g n ( p f n ( 0 ) ) .
1: repeat
2: Solve the problem (23) to obtain the optimal solution p f n * for given p f n ( v ) and η n ( c ) .
3: v = v + 1 .
4: Set p f n ( v ) = p f n * and cacluate f n ( p f n ( v ) ) η n ( c ) g n ( p f n ( v ) ) .
5: until f n ( p f n ( v ) ) η n ( c ) g n ( p f n ( v ) ) ( f n ( p f n ( v 1 ) ) η n ( c ) g n ( p f n ( v 1 ) ) ) μ
6: p f n * = p f n ( v ) .

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 C m a x . 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 O ( S n 3 ) [37], where S n is the number of MUs at SC n . The complexity of Algorithm 3 is X = V O ( S n 3 ) , where V is the the number of iterations bounded by V m a x . In summary, the computational complexity of the power control algorithm is O ( C X ) .

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 y n , as given by
y n = B f F log 2 1 + k S f n p k f n H k f n 1 + ϕ k f n .
if we fix y n , the objective function (15a) can be recast as
max p f n y n f F k S f n p k f n s . t . f F log 2 1 + k S f n p k f n H k f n 1 + ϕ k f n y n .
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
max p f n F + ( p f n ) F ˜ ( p f n ) f F k S f n p k f n + P c
s . t . F + ( p f n ) F ˜ ( p f n ) y n , n N ,
( 7 c ) , ( 11 ) .
It can be observed that since y n is fixed, (26) is equivalent to minimize the linear function f F k S f n p k f n + P c in the denominator, subject to convex constraints. Then, problem (26) can be solved by plain convex programming. To implement an efficient line search for y n , the bound of y n is given by
y n = F S f n R r e q f F log 2 1 + k S f n p k f n H k f n 1 + ϕ k f n < f F log 2 1 + k S f n P max H k f n = y n .
Then, the optimial p f n can be obtained by searching an appropriate value of y n 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 P m a x in Figure 2. The latter could be obtained in the first iteration of Algorithm 2 due to q ( 1 ) = 0 . In order to reflect the influence of q m a x , we gave four schemes of different q m a x . 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 P m a x 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 q m a x increased, because the increase of q m a x 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 R r e q 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 R r e q = 150 bps , 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 q m a x 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:
ε > 0 , ω = y n y n ε
1: for y f n y f n : ω : y f n do
2: p f n * = argmin f F k S f n p k f n + P c
3: end for
4: Obtain the optimal solution of p f n .
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 η n ( c ) = 1 . In Figure 4, we can see that the number of iterations are limited within four times. To show the influence of p f n ( 0 ) , we give the different values of p f n ( 0 ) in Figure 5, where f n ( p f n ( v ) ) η n ( c ) g n ( p f n ( v ) ) ( f n ( p f n ( v 1 ) ) η n ( c ) g n ( p f n ( v 1 ) ) ) = W ( p f n ( v ) ) . It is shown that the initial values have an effect on the number of iterations. Specifically, when p f n ( 0 ) = 0 × P max , 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 P m a x grows. From the four schemes, we can find out that the EE of K = 40 is much larger than that of K = 15 , 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 U k = 2 is larger than that of U k = 3 under the same K, i.e., K has much greater impact on EE than U k , since the NOMA scheme can obtain much higher EE by supporting multiple MUs and the increasing U k 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, K = 40 , N = 3 needed more swap operations than that K = 40 , N = 2 and K = 15 , N = 3 . Especially, less than 70 swap operations were needed for K = 40 and N = 3 .

6. Conclusions

This work investigated the uplink EE maximization problem in M2M-enabled HetNets with NOMA, where a MU acting as an MTCG can decode and forward both the information of MTCDs and its own data to the BS directly. Due to the limited spectrum resource, each BS shared the same sub-channels and NOMA was adopted between MUs in the same BS and sub-channel. The EE maximization problem was formulated, where MU association and power control were combined with each other. To solve it, a MU association matching algorithm was proposed based on the matching game. Under a given MU association, the uplink EE maximization was transformed into the EE maximization of each sub-channel. Two power control algorithms were provided to obtain the suboptimal power solutions based on sequential optimization. Simulation results showed that our proposed algorithms performed better than EE performance. It is known that cellular network is a key way to connect the M2M communications to the core network; our proposed scheme provided a new strategy for MTCDs to connect the cellular network with regard to MUs as their MTCGs based NOMA, and the power control of MTCDs was also considered as the constraints for the EE optimization problem. In fact, large scale devices are a more realistic scenario for 5G and next generation network, and since the number of MTCDs is considered on a small scale in this paper, the extension of our algorithms for large scale devices is one of the future works. Furthermore, the research of high computation complexity of the proposed algorithms for large scale devices is also a significant problem.

Author Contributions

The main contributions of S.Z. were to create the main ideas and execute performance evaluation by extensive simulation, while G.K. worked as the advisor to discuss, create, and advise the main ideas and performance evaluations together. All authors read and approved the final version of the paper.

Funding

This work was supported by the National Science and Technology Major Project of China (No. 2017ZX03001022)and National Natural Science Foundation of China (61471064).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ghosh, A.; Mangalvedhe, N.; Ratasuk, R.; Mondal, B.; Cudak, M.; Visotsky, E.; Thomas, T.A.; Andrews, J.G.; Xia, P.; Jo, H.S.; et al. Heterogeneous cellular networks: From theory to practice. IEEE Commun. Mag. 2012, 6, 54–64. [Google Scholar] [CrossRef]
  2. Damnjanovic, A.; Montojo, J.; Wei, Y.; Ji, T.; Luo, T.; Vajapeyam, M.; Yoo, T.; Song, O.; Malladi, D. A survey on 3GPP heterogeneous networks. IEEE Wirel. Commun. 2011, 3, 10–21. [Google Scholar] [CrossRef]
  3. Ericsson, A.B. 5G Radio Access, Research and Vision, Saito, Y., Benjebbour, A., Eds.; June 2013.
  4. Shin, W.; Vaezi, M.; Lee, B.; Love, D.J.; Lee, J.; Poor, H.V. Non-orthogonal multiple access in multi-cell networks: Theory, performance, and practical challenges. IEEE Commun. Mag. 2017, 10, 2–9. [Google Scholar] [CrossRef]
  5. Ding, Z.; Yang, Z.; Fan, P.; Poor, H.V. On the performance of non-orthogonal multiple access in 5G systems with randomly deployed users. IEEE Signal Process. Lett. 2014, 12, 1501–1505. [Google Scholar] [CrossRef]
  6. Ding, Z.; Peng, M.; Poor, H.V. Cooperative non-orthogonal multiple access in 5G systems. IEEE Commun. Lett. 2015, 8, 1462–1465. [Google Scholar] [CrossRef]
  7. Zhang, N.; Wang, J.; Kang, G.; Liu, Y. Uplink non-orthogonal multiple access in 5G systems. IEEE Commun. Lett. 2016, 3, 458–461. [Google Scholar] [CrossRef]
  8. Zhang, H.; Fang, F.; Cheng, J.; Long, K.; Wang, W.; Leung, V.C.M. Energy-efficient resource allocation in NOMA heterogeneous networks. IEEE Wirel. Commun. 2018, 2, 48–53. [Google Scholar] [CrossRef]
  9. Mokdad, A.; Azmi, P.; Mokari, N. Radio resource allocation for heterogeneous traffic in GFDM-NOMA heterogeneous cellular networks. IET Commun. 2016, 12, 1444–1455. [Google Scholar] [CrossRef]
  10. Zhao, J.; Liu, Y.; Chai, K.K.; Nallanathan, A.; Chen, Y.; Han, Z. Spectrum allocation and power control for non-orthogonal multiple access in HetNets. IEEE Trans. Wirel. Commun. 2017, 9, 5825–5837. [Google Scholar] [CrossRef]
  11. Ni, D.; Hao, L.; Tran, Q.T.; Qian, X. Power allocation for downlink NOMA heterogeneous networks. IEEE Access 2018, 6, 26742–26752. [Google Scholar] [CrossRef]
  12. Ni, D.; Hao, L.; Tran, Q.T.; Qian, X. Transmit power minimization for downlink multi-cell multi-carrier NOMA networks. IEEE Commun. Lett. 2018, 12, 2459–2462. [Google Scholar] [CrossRef]
  13. Fang, F.; Cheng, J.; Ding, Z. Joint energy efficient subchannel and power optimization for a downlink NOMA heterogeneous network. IEEE Trans. Veh. Technol. 2019, 68, 1351–1364. [Google Scholar] [CrossRef]
  14. Moltafet, M.; Azmi, P.; Mokari, N.; Javan, M.R.; Mokdad, A. Optimal and fair energy efficient resource allocation for energy harvesting-enabled-PD-NOMA-based HetNets. IEEE Trans. Wirel. Commun. 2018, 3, 2054–2067. [Google Scholar] [CrossRef]
  15. Pentikousis, K. In search of energy-efficient mobile networking. IEEE Commun. Mag. 2010, 1, 95–103. [Google Scholar] [CrossRef]
  16. IMT-2020(5G) Promotion Group. 5G Vision and Requirements, May 2014.
  17. Buzzi, S.; Chih-Lin, I.; Klein, T.E.; Poor, H.V.; Yang, C.; Zappone, A. A survey of energy-efficient techniques for 5G networks and challenges ahead. IEEE J. Sel. Areas Commun. 2016, 4, 697–709. [Google Scholar] [CrossRef]
  18. Qin, Z.; Yue, X.; Liu, Y.; Ding, Z.; Nallanathan, A. User association and resource allocation in unified NOMA enabled heterogeneous ultra dense networks. IEEE Commun. Mag. 2018, 6, 86–92. [Google Scholar] [CrossRef]
  19. Yang, Z.; Xu, W.; Xu, H.; Shi, J.; Chen, M. Energy efficient non-orthogonal multiple access for machine-to-machine communications. IEEE Commun. Lett. 2017, 4, 817–820. [Google Scholar] [CrossRef]
  20. Shirvanimoghaddam, M.M.; Dohler, M.; Johnson, S.J. Massive non-orthogonal multiple access for cellular IoT: Potentials and limitations. IEEE Commun. Mag. 2017, 9, 55–61. [Google Scholar] [CrossRef]
  21. Chen, S.; Ma, R.; Chen, H.; Zhang, H.; Meng, W.; Liu, J. Machine-to-machine communications in ultra-dense networks-a survey. IEEE Commun. Surv. Tutor. 2017, 3, 1478–1503. [Google Scholar] [CrossRef]
  22. Yang, Z.; Pan, Y.; Xu, W.; Guan, R.; Wang, Y.; Chen, M. Energy efficient resource allocation for machine-to-machine communications with NOMA and energy harvesting. In Proceedings of the 2017 IEEE Conference on Computer Communications Workshops (INFOCOM WKSHPS), Atlanta, GA, USA, 1–4 May 2017; pp. 1–6. [Google Scholar] [CrossRef]
  23. Xu, B.; Chen, Y.; Carrion, J.R.; Zhang, T. Resource allocation in energy-cooperation enabled two-tier NOMA HetNets towards green 5G. IEEE J. Sel. Areas Commun. 2017, 12, 2758–2770. [Google Scholar] [CrossRef]
  24. Qian, L.P.; Wu, Y.; Zhou, H.; Shen, X.S. Joint uplink base station sssociation and power control for small-cell networks with non-orthogonal multiple access. IEEE Trans. Wirel. Commun. 2017, 19, 5567–5582. [Google Scholar] [CrossRef]
  25. Yang, Z.; Xu, W.; Pan, Y.; Pan, C.; Chen, M. Energy efficient resource allocation in machine-to-machine communications with multiple access and energy harvesting for IoT. IEEE Internet Things J. 2018, 5, 229–245. [Google Scholar] [CrossRef]
  26. Li, Z.; Gui, J. Energy-efficient resource allocation with hybrid TDMA–NOMA for cellular-enabled machine-to-machine communications. IEEE Access 2019, 7, 105800–105815. [Google Scholar] [CrossRef]
  27. Lv, T.; Ma, Y.; Zeng, J.; Mathiopoulos, P.T. Millimeter-wave NOMA transmission in cellular M2M communications for Internet of Things. IEEE Internet Things J. 2018, 5, 1989–2000. [Google Scholar] [CrossRef]
  28. Shahini, A.; Ansari, N. NOMA aided narrowband IoT for machine type communications with user clustering. IEEE Internet Things J. 2019, 6, 7183–7191. [Google Scholar] [CrossRef]
  29. Han, S.; Xu, X.; Tao, X.; Zhang, P. Joint power and sub-channel allocation for secure transmission in NOMA-based mMTC networks. IEEE Syst. J. 2019, 13, 2476–2487. [Google Scholar] [CrossRef]
  30. Mohamed, A.; Onireti, O.; Imran, M.A.; Imran, A.; Tafazolli, R. Control-data separation architecture for cellular radio access networks: A survey and outlook. IEEE Commun. Surv. Tutor. 2016, 18, 446–465. [Google Scholar] [CrossRef]
  31. Bodine-Baron, E.; Lee, C.; Chong, A.; Hassibi, B.; Wierman, A. Peer effects and stability in matching markets. In Proceedings of the 4th Symposium on Algorithmic Game Theory (SAGT), Amalfi, Italy, 17–19 October 2011; pp. 117–129. [Google Scholar] [CrossRef]
  32. Aydin, O.; Jorswieck, E.A.; Aziz, D.; Zappone, A. Energy-spectral efficiency tradeoffs in 5G multi-operator networks with heterogeneous constraints. IEEE Trans. Wirel. Commun. 2017, 9, 5869–5881. [Google Scholar] [CrossRef]
  33. Zappone, A.; Bjornson, E.; Sanguinetti, L.; Jorswieck, E. Globally optimal energy-efficient power control and receiver design in wireless networks. IEEE Trans. Signal Process. 2017, 11, 2844–2859. [Google Scholar] [CrossRef]
  34. Dinkelbach, W. On nonlinear fractional programming. Manag. Sci. 1967, 7, 492–498. [Google Scholar] [CrossRef]
  35. Ali, M.S.; Tabassum, H.; Hossain, E. Dynamic user clustering and power allocation for uplink and downlink non-orthogonal multiple access (NOMA) systems. IEEE Access 2016, 4, 6325–6343. [Google Scholar] [CrossRef]
  36. Li, Y.; Sheng, M.; Yang, C.; Wang, X. Energy efficiency and spectral efficiency tradeoff in interference-limited wireless networks. IEEE Commun. Lett. 2013, 10, 1924–1927. [Google Scholar] [CrossRef]
  37. Kha, H.H.; Tuan, H.D.; Nguyen, H.H. Fast global optimal power allocation in wireless networks by local D.C. programming. IEEE Trans. Wirel. Commun. 2012, 2, 510–515. [Google Scholar] [CrossRef]
Figure 1. System Model.
Figure 1. System Model.
Sensors 19 05307 g001
Figure 2. EE versus the maximum transmit power for different schemes.
Figure 2. EE versus the maximum transmit power for different schemes.
Sensors 19 05307 g002
Figure 3. EE versus the date rate requirement of each MU for different schemes.
Figure 3. EE versus the date rate requirement of each MU for different schemes.
Sensors 19 05307 g003
Figure 4. Convergence property of Algorithm 2.
Figure 4. Convergence property of Algorithm 2.
Sensors 19 05307 g004
Figure 5. Convergence property of Algorithm 3.
Figure 5. Convergence property of Algorithm 3.
Sensors 19 05307 g005
Figure 6. EE versus the maximum transmit power for the different number of MUs and MTCDs.
Figure 6. EE versus the maximum transmit power for the different number of MUs and MTCDs.
Sensors 19 05307 g006
Figure 7. CDF of the number of swap operations for convergence.
Figure 7. CDF of the number of swap operations for convergence.
Sensors 19 05307 g007
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParametersMeaningsValues
FNumber of BSs3
BThe frequency bandwidth of each sub-channel 15 kHz
KNumber of MUs40
U k Number of MTCDs of each MU2
σ f n 2 Noise variance 2 dBm
μ , τ Error tolerance 10 3
P m a x The maximum of transmit power of MU 0.2 W
q m a x The maximum of transmit power of MTCD 0.08 W
α Path loss factor3
P c The circuit power at each sub-channel 0.1 W
R r e q The data rate requirement of each MU 100 bps

Share and Cite

MDPI and ACS Style

Zhang, S.; Kang, G. User Association and Power Control for Energy Efficiency Maximization in M2M-Enabled Uplink Heterogeneous Networks with NOMA. Sensors 2019, 19, 5307. https://doi.org/10.3390/s19235307

AMA Style

Zhang S, Kang G. User Association and Power Control for Energy Efficiency Maximization in M2M-Enabled Uplink Heterogeneous Networks with NOMA. Sensors. 2019; 19(23):5307. https://doi.org/10.3390/s19235307

Chicago/Turabian Style

Zhang, Shuang, and Guixia Kang. 2019. "User Association and Power Control for Energy Efficiency Maximization in M2M-Enabled Uplink Heterogeneous Networks with NOMA" Sensors 19, no. 23: 5307. https://doi.org/10.3390/s19235307

APA Style

Zhang, S., & Kang, G. (2019). User Association and Power Control for Energy Efficiency Maximization in M2M-Enabled Uplink Heterogeneous Networks with NOMA. Sensors, 19(23), 5307. https://doi.org/10.3390/s19235307

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop