Abstract

In millimeter wave (mmWave) communication systems, beamforming-enabled directional transmission and network densification are usually used to overcome severe signal path loss problem and improve signal coverage quality. The combination of directional transmission and network densification poses a challenge to radio access resource management. The existing work presented an effective solution for dense mmWave wireless local area networks (WLANs). However, this scheme cannot adapt to network expansion when it is applied directly to dense mmWave cellular networks. In addition, there is still room for improvement in terms of energy efficiency and throughput. Therefore, we firstly propose an efficient hierarchical beamforming training (BFT) mechanism to establish directional links, which allows all the small cell base stations (SBSs) to participate in the merging of training frames to adapt to network expansion. Then, we design a BFT information-aided radio access resource allocation algorithm to improve the downlink energy efficiency of the entire mmWave cellular network by reasonably selecting beam directions and optimizing transmission powers and beam widths. Simulation results show that the proposed hierarchical BFT mechanism has the smaller overhead of BFT than the existing BFT mechanism, and the proposed BFT information-aided radio access resource allocation algorithm outperforms the existing corresponding algorithm in terms of average energy efficiency and throughput per link.

1. Introduction

The growing popularity of the mobile Internet and the Internet of things [1–3] as well as the increasing perceptual services [4–6] has led to a continued high demand for network capacity, in which more spectrum is needed to support [7]. In order to relieve the pressure of spectrum resource supply, various schemes for efficiently exploiting spectrum have been proposed [8–11], but they still cannot satisfy people’s thirst for high capacity and massive connectivity [12, 13]. Since spectrum resources in the traditional cellular frequency bands have dried up, millimeter wave (mmWave) band has become a new source of spectrum supply [14], which has rich spectrum resources and thus can support the wireless transmissions with high data rate. Now, mmWave communication technology has been applied to both the fifth-generation (5G) cellular systems and wireless local area networks (WLANs) (e.g., IEEE 802.11ad and IEEE 802.11ay) [15]. However, before making full use of its abundant spectral resources, it is necessary to overcome its inherent defects (e.g., severe path loss, high oxygen absorption, and blockages). Fortunately, the short mmWave wavelength facilitates the deployment of a large number of antennas on a small-sized terminal, and thus, it is also helpful to improve communication quality or extend communication distance by using antenna array technology and directional transmission.

In addition, by densely deploying base stations, the average communication distance can be reduced, and the signal coverage quality can also be improved because the probability of blockages goes down [16]. However, the combination of directional transmission and network densification will create new difficulties for radio access resource management, where improper mmWave beam allocation will lead to nonnegligible interferences between communication nodes. The authors of [17] believe that the beam management and interference coordination (BM-IC) in dense mmWave network is very challenging, while some researchers focus on the evolution of WLANs and the corresponding resource management, beam management, interference management, and mobility management [18–20].

Unlike the above studies that focus on the sum data rate improvement of the entire mmWave WLAN by proper mmWave beam allocation, we focus on the energy efficiency improvement of the entire mmWave cellular network, since the energy consumption problem becomes more prominent in high-capacity network applications. The main contributions of this paper mainly include two aspects. Firstly, in order to reduce the overhead of beamforming training (BFT) in a dense mmWave cellular network, inspired by the BFT mechanism in [19], we propose an efficient hierarchical BFT mechanism to establish directional links. Secondly, we model the radio access resource management for downlink communication in the dense mmWave cellular network as an optimization problem and design a hierarchical BFT information-aided radio access resource allocation algorithm to improve the downlink energy efficiency of the entire mmWave cellular network by reasonably selecting beam directions and optimizing transmission powers and beam widths.

The remainders of this paper are organized as follows. In Section 2, we review related work with respect to radio resource management (RRM). The network model and beamforming training are described in Section 3, while the problem formulation is given in Section 4. The design details regarding the proposed scheme are proposed in Section 5. The simulation results are evaluated in Section 6. Finally, the conclusion of this paper is summarized in Section 7.

2. Related Work

In this section, we mainly introduce the existing studies regarding RRM. For WLANs, some researchers think that distributed network architecture makes RRM very inefficient while centralized network architecture is helpful to achieve more efficient RRM [21–23]. Unlike these schemes that all work at low-frequency band and adopt the wired links between controller and access points (APs), in dense mmWave WLANs, a controller is used for managing and controlling multiple APs by directional mmWave links [24, 25]. The authors of [18–20] believed that cloud radio access network- (C-RAN-) based centralized control mode should be applied to WLANs. Therefore, efficient and smart algorithms are needed for mmWave beam management and interference management in WLANs. For example, a deep learning-based BM-IC method in dense mmWave WLANs was proposed in [19]. Based on the centralized architecture, the authors first proposed that an efficient BFT mechanism is to reduce the BFT overhead. Then, they formulated the BM-IC as an optimization problem and designed a BFT information-aided BM-IC algorithm to solve it. Finally, they trained a deep neural network (DNN) model to approximate the proposed algorithm by using the training data obtained from the proposed algorithm.

Since cellular networks also exploit mmWave frequencies to increase network capacity, some literatures reviewed the RRM studies on dense mmWave cellular networks. The authors of [26] reviewed the mmWave beam management procedures included in the 3GPP New Radio (NR) specifications based on both standalone and nonstandalone architectures and provided the performance comparison among different schemes and the design insights into the most important parameters associated with the beam management frameworks. The authors of [27] reviewed the existing measurement techniques for beam and mobility management in mmWave cellular networks and proposed the insights into the design of accurate, reactive, and robust control layer procedures (e.g., initial access, handover, and beam tracking). The authors of [28] reviewed the existing economic and pricing approaches for addressing resource management issues in the 5G wireless networks including user association, spectrum allocation, and interference and power management.

Some studies focused on BM-IC problem in dense mmWave cellular networks. The authors of [16] explored the power control and transmission duration allocation problem in self-backhauling dense mmWave cellular networks. The authors of [17] explored the potential gain of ultradensification in mmWave cellular networks and proposed a large-scale channel state information- (CSI-) based interference coordination approach to solve the cochannel interference. The authors of [29] explored the beam pair selection problem in heterogeneous multibeam cloud radio access network architecture, which is aimed at improving the mmWave cellular network performance. The authors of [30] studied the joint power control and user association problem in mmWave cellular communications with directional transmission and dense deployment and proposed a reinforcement learning framework to address this problem. The authors of [31] explored resource allocation problem in 5G mmWave small cell systems, presented a hybrid beamforming structure for small cell base station (SBS), and formulated the joint allocation of the beam and transmission power on the subcarriers into a mixed integer nonlinear programming problem and decomposed it into the beam selection and power allocation subproblems.

Exploiting the benefits from dense mmWave networks with multiple SBSs and many user equipment (UE), the authors of [32] discussed an approach for jointly optimizing SBS-UE association and power allocation, which can maximize the system energy efficiency (EE) while guaranteeing the quality of service (QoS) constraints for each UE. However, the mechanisms mentioned above do not consider the BFT problem in dense mmWave cellular networks. Since users and base stations need to align their beams during both initial access and data transmissions to ensure their end-to-end performance and quality of service, an effective BFT mechanism is needed to ensure the accuracy of beam selection and reduce the delay in performing initial access or updating beam pair in dense mmWave cellular networks.

3. Network Model and Beamforming Training

3.1. Network Model

As shown in Figure 1, there are SBSs overlapped in a micro cell and UEs randomly distributed in the same macro cell, where a micro base station (MBS) is located at the center of this micro cell. We denote the set of SBSs as and the set of UEs as , respectively. Also, there are APs overlapped in the small cells. We mark the set of APs as .

Figure 1 

A dense mmWave cellular network.

In this paper, a SBS is indirectly connected to the core network via the MBS in a wireless mode (e.g., mmWave communication) or directly connected to the core network in a wired mode (e.g., optical fiber), while an AP is only indirectly connected to the core network via a SBS in a mmWave communication mode.

Each AP can schedule a set of mmWave frequency resources and allocate them to the UEs associated by it in an orthogonal mode. The same mmWave frequency resources can be reused by the other APs. Also, these resources are shared among the MBS and SBSs in an underlay mode for the purpose of providing wireless access service.

Each UE selects one of the APs (or SBSs) around itself to perform channel access and data transmission through directional link (i.e., the gray beams in Figure 1), where the MBS is not selected for direct access by any UE, so it can focus on providing backhaul connection for SBSs and control plane information transmission. When the MBS and SBSs provide wireless backhaul service, they will use dedicated mmWave frequency bands to avoid interfering with each other. Based on the assumption of unrestricted backhaul capacity, we only focus on radio access resource management problem for downlink unicast communications in a dense mmWave cellular network.

To facilitate the evaluation of energy efficiency for dense mmWave cellular networks, we adopt the commonly used switch-based analog beam pattern [33]. The normalized beamforming gain is as follows. where is the beam width of the main lobe in radian, is the beam offset angle to the main lobe in radian, and is the gain of the side lobe and .

For analytical tractability, we assume that the number of beams of each of MBS, SBS, AP, and UE is limited, and each beam covers a unique direction in a nonoverlapping way. Therefore, the number of beams of each of MBS, SBS, AP, and UE is denoted as , which is estimated as follows. where is the maximum beam width of MBS, SBS, AP, and UE.

3.2. Efficient Beamforming Training Mechanism for Dense mmWave Cellular Network

As shown in Figure 1, when a UE is covered by multiple APs (e.g., UE is covered by AP and AP ), in order to allow the UE to select a suitable AP (e.g., AP ) among all APs for access in terms of smaller interference, the UE should perform beamforming training (BFT) with its surrounding APs in advance.

Since the BFT mechanism in IEEE 802.11ad/ay only considers the BFT between one AP and one wireless terminal or between one AP and multiple wireless terminals, it will generate serious collision problem if multiple APs are allowed to perform BFT simultaneously [18]. When this BFT mechanism is directly applied to dense mmWave networks, all APs have to perform BFT one by one to avoid collision, which makes BFT be inefficient. Therefore, the authors in [19] proposed an efficient BFT mechanism to coordinate the BFT processes of every AP and wireless terminals so that many duplicated BFT processes can be merged. However, the merge decision is made by a single node (e.g., a MBS), which is difficult to adapt to an increasing number of wireless terminals. It may be an effective solution to assign merge decision tasks to multiple nodes (e.g., SBSs). Therefore, the basic idea of our proposed efficient hierarchical BFT mechanism is shown as below.

Step 1. (BFT between the MBS and its surrounding SBSs). This step is initiated by the MBS, and the time slot structure of its training process is shown in Figure 2.

Figure 2 

The time slot structure of the training initiator.

In the transmitting training request phase, the MBS simultaneously transmits the training request beacon frames in all its beams to start the training process between the MBS and SBSs.

In the waiting for training response phase, each SBS that receives request beacon frame randomly selects a response slot to simultaneously transmit its training response beacon frames in all its beams.

If the response beacon frame conflict (that refers to the occurrence of more than one response beacon frame on the same beam at the MBS side) does not occur, the handling response conflicts phase can be ignored. Otherwise, the MBS simultaneously retransmits the training request beacon frames in the beams where there are conflicting events. In these retransmitted frames, the number of response slots is doubled to reduce the probability of conflict.

In the announcing training results phase, the MBS transmits the beacon frames for announcing training results to the SBSs that have successfully sent the training response beacon frames through the best beams estimated by it according to the information of the training response beacon frames.

The structure of beacon frames used in the above training process is shown in Figure 3. In the training request beacon frame, INI_ID represents the initiator’s identifier (e.g., the MBS); REQ_NUM denotes the serial number of this frame; REQ_BEAM_ID is the serial number of the beam in which this frame is transmitted; NARS records the number of allowed response slots; REQ_TX_POW is the power used to transmit this frame.

(a) The training request beacon frame

(b) The training response beacon frame

(c) The beacon frame of announcing training results

(a) The training request beacon frame
(b) The training response beacon frame
(c) The beacon frame of announcing training results

Figure 3 

The structure of beacon frames.

In the training response beacon frame, REP_ID represents the responder’s identifier (e.g., a SBS); INI_NUM denotes a pair of values, which comes from the INI_ID and REQ_NUM of the training request beacon frame received by this responder; REP_BEAM_ID is the serial number of the beam in which this frame is transmitted; PLV denotes the path loss value, which is estimated according to the REQ_TX_POW in the training request beacon frame received by this responder and the receiving signal strength perceived by this responder; REP_TX_POW is the power used to transmit this frame, which facilitates the initiator to estimate the path loss value.

In the beacon frame of announcing training results, SNFB is the serial number of frequency band, that is to say, it denotes the channel between the initiator (i.e., INI_ID) and the responder (i.e., REP_ID), where the initiator’s beam (i.e., REQ_BEAM_ID) can cover the responder and the responder’s beam (i.e., REP_BEAM_ID) can cover the initiator. This pair of beams should minimize the weighted sum of the two-way path loss values between the initiator and the responder. Otherwise, this pair of beams will not be placed in this frame. If different frequency bands are used on the same pair of beams, concurrent transmission can be achieved without conflict.

Step 2. (BFT between every SBS and its surrounding APs). This step is initiated by every SBS, and the time slot structure of its training process is similar to that of Figure 2. Since all SBSs must perform BFT one by one to avoid collision, there may be the case where an AP is associated with more than one SBS. Therefore, each SBS firstly transmits its built beacon frame of announcing training results to the MBS via mmWave link or optical fiber. Based on the results reported by all SBSs, the MBS performs comparative analysis and produces results and then feedbacks the APs with which each SBS should associate. After that, each SBS sends each expected associated AP the corresponding beacon frame.

Step 3. (BFT between every AP and its surrounding UEs). This step is initiated by every AP, and the time slot structure of its training process is also similar to that of Figure 2. Similar to Step 2, since all APs must perform BFT one by one to avoid collision, there may be the case where an UE is associated with more than one AP. Therefore, each AP firstly transmits its built beacon frame of announcing training results to its associated SBS (that is determined in Step 2). Based on the results reported by all the APs associated by a SBS, this SBS performs comparative analysis and produces results and then feedbacks the UEs with which each AP should associate. After that, each AP sends each expected associated UE the corresponding beacon frame.

By Steps 1–3, we only finish the transmitting BFT of SBSs, APs, and UEs instead of their receiving BFT. Due to the reciprocity of directional beams in mmWave band, both transmitting BFT and receiving BFT can be done in a similar way [15]. In order to finish the receiving BFT, the beam refinement protocol phase in IEEE 802.11ad/ay can be performed by applying the similar modifications of Steps 1–3. For more details (e.g., training processes and specific frame formats), please refer to [18].

3.3. Performance Analysis of Beamforming Training

Without loss of generality, we assume that each of SBSs, APs, and UEs can generate beams. Therefore, the total number of training frames required for our proposed efficient hierarchical BFT mechanism should be at least as follows.

If the BFT mechanism in [19] is used in this paper, the total number of training frames should be at least as follows.

Therefore, the total number of training frames saved by our proposed efficient hierarchical BFT mechanism is estimated by the following formula.

Compared to the BFT mechanism in [19], our proposed efficient hierarchical BFT mechanism can reduce BFT overhead. More importantly, the proposed BFT mechanism will be more efficient when the number of UEs (i.e., ) is much larger than that of APs (i.e., ).

The proof of equation (3) is given as follows. When the MBS performs the BFT with SBSs, where we assume that SBSs can be directly connected with the MBS via mmWave links and , it should transmit training frames in its beams to complete the transmission of training request phase in Step 1. In the waiting for training response phase, each of the SBSs should transmit training frames in a randomly selected response slot. If there is no conflict of response beacon frames, the MBS should transmit beacon frames for announcing training results to SBSs. Therefore, there are at least training frames required for completing Step 1 in Section 3.2.

For Step 2 in Section 3.2, SBSs should transmit training frames one by one to complete the transmission of training request phase. In the waiting for training response phase, each of the APs (e.g., AP , ) should transmit training frames in a randomly selected response slot. After the AP transmits training frames, SBSs that receive the training frames from the AP should transmit beacon frames for announcing training results to the MBS. Soon afterwards, the MBS transmits beacon frames for announcing training results to SBSs. Finally, the selected SBS (based on the signal-to-noise ratio (SNR) maximization) will transmit the beacon frames for announcing training results to APs. Therefore, APs require training frames to complete the transmission of response beacon frames in total. Therefore, there are at least training frames required for completing Step 2 in Section 3.2.

For Step 3 in Section 3.2, APs should transmit training frames one by one to complete the transmission of the training request phase. In the waiting for training response phase, each of the UEs (e.g., UE , ) should transmit training frames in a randomly selected response slot. After UE transmits training frames, APs that receive the training frames from UE should transmit beacon frames for announcing training results to their associated SBSs, where they are likely to be associated with several adjacent SBSs, even the same SBS. Soon afterwards, these SBSs transmit beacon frames for announcing training results to APs. Finally, the selected AP (based on the SNR maximization) will transmit the beacon frames for announcing training results to UEs. Therefore, UEs require training frames to complete the transmission of response beacon frames in total. Therefore, there are at least training frames required for completing Step 3 in Section 3.2. According to the results of Steps 1–3, equation (3) is true.

The proof of equation (4) is given as follows. In order to perform the BFT mechanism in [19], we assume that the MBS, SBSs, and APs can connect to each other via mmWave links, where the whole training process consists of only two steps. In the first step, the MBS performs the BFT with SBSs and APs. Thus, it should transmit training frames in its beams to complete the transmission of the training request phase. In the waiting for training response phase, each of the SBSs and APs should transmit training frames in a randomly selected response slot. If there is no conflict of response beacon frames, the MBS should transmit beacon frames for announcing training results. Therefore, there are at least training frames required for the first step.

In the second step, SBSs and APs should transmit training frames one by one to complete the transmission of the training request phase. In the waiting for training response phase, each of the UEs should transmit training frames in a randomly selected response slot. Then, SBSs and APs should transmit beacon frames for announcing training results to the MBS. Soon afterwards, the MBS sends back beacon frames for announcing training results. Finally, the selected AP or SBS (based on the SNR maximization) will transmit the beacon frames for announcing training results to UEs. So training frames are required in the waiting for the training response phase. Therefore, there are at least training frames required in the second step. According to the results of the first step and the second step, equation (4) is true.

Based on the above analysis, we know that the original two-stage BFT mechanism is based on the assumption that all the nodes can be directly connected to each other. Obviously, this assumption does not always hold in the dense mmWave cellular network scenario, where some nodes may be indirectly connected to each other. If the original two-stage BFT mechanism is applied to the dense mmWave cellular network scenario, the number of training frames would be also counted between the nodes that are not directly connected to each other. Therefore, there are more training frames from the perspective of theoretical statistics, which gives the illusion of more BFT overhead since the unnecessary training frames are counted.

Moreover, there are not only APs but also SBSs in the dense mmWave cellular network scenario, and also, there are the more UEs than the dense mmWave WLAN. Therefore, from formula (5), we know that the original two-stage BFT mechanism is really unsuitable for the dense mmWave cellular network scenario. As the foregoing, the original two-stage BFT mechanism relies on a single MBS to merge the duplicated BFT processes to reduce the number of BFT frames. If it is applied to the dense mmWave cellular network scenario, it would make the MBS become a performance bottleneck. By having all the SBS and the MBS work together to merge the duplicated BFT processes, our three-stage BFT mechanism can avoid this performance bottleneck.

4. Problem Formulation

In a dense mmWave cellular network, if the beams between UE and AP are aligned, according to equation (1), we can obtain the directional transmitting gain and directional receiving gain as where is the beam width of the transmitter and is the beam width of the receiver. According to [19, 34, 35], the channel between UE and AP is given by where PH is the number of paths and the superscript “(ph)” represents the ph-th path. and , as shown in equation (6), are the directional transmitting gain and directional receiving gain of the ph-th path between UE and AP is the amplitude of the ph-th path, and is the Dirac delta function. is the propagation delay of the ph-th path, which is estimated by the following formula. where is the distance of the ph-th path between UE and AP and is the speed of light. For analytical tractability, we assume the line-of-sight (LOS) path (i.e., ) always exists [22] and then the rest of the PH paths (i.e., from to ) are non-line-of-sight (NLOS). As for the LOS path, the amplitude is estimated by the following formula [35]. where is the wavelength and and is the carrier frequency. is the distance of the LOS path between UE and AP . As for the NLOS paths, the amplitude is related to both path loss and reflection coefficients, which is estimated by the following formula [35]. where is the reflection coefficient of the rph-th reflection for the ph-th path and RPH is the number of reflections of the path. Since the reflection loss is very high in the mmWave band [15], we only consider one reflection of a given path (i.e., ). Therefore, formula (10) is simplified as follows.

Finally, according to [36], the channel gain can be derived by

Let the directional beam’s transmission power from AP to UE be , and then, at the UE side, the received power from AP can be determined by

Similarly, the interference power received at UE from the other APs (e.g., AP ) and UEs (e.g., UE ) should be where and are the directional transmitting gain and directional receiving gain from the ph-th path between UE and AP , respectively. According to formula (1), the directional transmit-receive gain of each path can be derived by the following formula.

Let and be the beam offset angle from AP ’s (AP transmits to UE ) transmitting beam direction to the position of UE and that from the UE ’s (UE receives from AP ) receiving beam direction to the position of AP , respectively; Condition A is both and ; Condition B is both and ; Condition C is both and ; Condition D is both and .

We assume that there is no multiconnectivity capability for UEs, and thus, a UE can only connect to one AP at a time in a dense mmWave cellular network. We set as the binary association variable. If UE is connected to AP , then ; otherwise . When UE connects to AP , according to formulas (13) and (14), the signal-to-interference plus noise ratio (SINR) between AP and UE is denoted as , which can be derived by where and represent the bandwidth of mmWave band and the background noise power spectrum density, respectively. Then, the data rate (throughput) of the link from AP to UE is as follows.

The sum data rate of the entire mmWave cellular network (i.e., only the data rate from APs to UEs) is

In order to maximize the sum data rate in (18), the radio access resource management problem for downlink unicast communications can be transformed into the following optimization problem.

In (19), the constraints are specified as follows. where is a matrix with element , is a matrix with element , and is a matrix with element or . Constraint indicates that each UE can only connect to one AP. Constraint ensures that the number of UEs that an AP can serve is not more than . Constraint means that the power consumption constraint for each AP cannot exceed . Constraint indicates that the beam width is limited to []. Constraint stipulates that for an AP, the sum of beam widths used for serving the connected UEs cannot exceed , since the beams of a node are assumed to be nonoverlapping in Section 3.1. However, the optimization problem could not be solved directly because it is NP-hard and not convex. Accordingly, we design an algorithm to find a suboptimal solution for problem in Section 5.

5. Radio Access Resource Management in Dense mmWave Cellular Network

In this section, we propose a hierarchical BFT information-aided radio access resource management algorithm to solve problem so as to reduce the interference and improve the energy efficiency of dense mmWave cellular network.

Since the path loss of mmWave band is serious and the transmission power is limited in practical scenarios, a UE may not pick up signals from all the APs in the network. That is to say, a UE may only connect to its available surrounding APs. Each AP may only get a set of candidate UEs based on the hierarchical BFT results in Section 3.2. We judge whether there is a direct communication link between a UE (e.g., ) and an AP (e.g., ) by observing the receiving bit error rate (BER).

According to the formula of estimating the bit transmission success rate in [37] (i.e., ; here, is UE ’s bit transmission success rate, and is the corresponding SINR value), we know that the corresponding BER value can be estimated by . When the acceptable BER in this paper is determined and denoted as , which is the BER threshold of UE ’s receiving data from AP , the corresponding SINR threshold (i.e., ) can be estimated by

In order to ensure that the BER level of UE ’s receiving data from AP is not higher than , the transmission power of AP should not be lower than the power threshold , which is estimated as follows.

If is not less than , we believe that there is a direct communication link to UE from AP . For a set of candidate UEs of AP (), the directional link’s quality (e.g., ) between each UE () and AP should be greater than a predefined threshold (e.g., ), where the set of candidate UEs of AP is denoted as . Therefore, expressions (19) and (20) can be rewritten as

In (23), the constraints are specified as follows.

In the proposed BFT, when a UE concurrently transmits its training response beacon frames in all its beams, its surrounding APs that can receive these frames will firstly transmit their built beacon frames for announcing training results to their associated SBSs. If these SBSs independently perform the merge operation of training results and then feedback the results to the corresponding APs, some of these APs will reply to this UE with a response packet, respectively. That is to say, this UE may receive the multiple APs’ response packets, and thus, it should decide which AP to associate. Since these SBSs are adjacent to each other and also their numbers are usually small, we believe that they can exchange merge operation results with each other through the dedicated channels at a small cost, which can ensure that only the AP that has the best channel quality between itself and this UE will transmit one response packet to this UE.

For each SBS (e.g., ), it can obtain all the SNR values (i.e., ) in its coverage when the proposed BFT is done. However, using SNR directly to measure channel quality is not always an effective method. Therefore, considering the data rate per unit of power consumption (DRPC) in a link to measure channel quality can better reflect the goal of high energy efficiency pursued by the next-generation wireless network, which is formulated as follows.

In (25), is the energy efficiency of the link from AP to UE , and is the power consumed by a radio frequency (RF) chain, where is set to 0.0344 W in [38].

For the sake of simplicity of analysis, we store all the DRPC values into in descending order. Then, there is a one-to-one match between the elements of and those of , where and stands for the number of elements in a set and represents the set of all the UEs associated by all the APs. In addition, the number of elements in , , and is also .

The problem described in (23) and (24) is also NP-hard and not convex, so it cannot be solved directly. The proposed algorithm decomposes this problem into three subproblems, which can find the suboptimal solution. Firstly, we initialize and to obtain the optimized . Then, we use the optimized and initialized to get the optimized . Finally, we adopt the optimized and to find the optimized .

In the proposed algorithm description shown in Algorithm 1, to obtain the optimized , and are firstly initialized to the average power (i.e., ) and the maximum beam width (i.e., ), respectively. After that, all the elements in are initialized to 0. Then, the algorithm selects the AP and UE pair corresponding to and then calculates the sum data rate according to and sets to 1. Next, the algorithm selects the AP and UE pair corresponding to and calculates the sum data rate . If , the algorithm sets to 1; otherwise, it sets to 0.

By repeating the above steps, the algorithm selects the AP and UE pair corresponding to and calculates the sum data rate . If , is set to 1; otherwise, is set to 0, until all the are set. Therefore, the algorithm can get the optimized after mapping to . Since indicates which AP each UE should connect to, the beam directions are determined.

After this stage, the algorithm can use the optimized and the maximum beam width vector to obtain power allocation vector . The algorithm starts from the directional link with the highest DRPC value determined by , , and . Then, it decreases the transmission power from the average power in . If the sum data rate increases, the algorithm continues to reduce the transmission power until it reaches ; otherwise, it optimizes the transmission power of the link with the second highest DRPC value. Repeating the above steps until all the are optimized, we can get after mapping to .

Finally, the optimized can be obtained by applying and into the process similar to the one above. That is to say, the algorithm starts from the directional link with the highest DRPC value determined by , , and . Then, it decreases beam width from . If the sum data rate increases, the algorithm continues to reduce beam width until it reaches ; otherwise, it optimizes the beam width of the link with the second highest DRPC value. Repeating the above steps until all the are optimized, we can get after mapping to or . Therefore, , , and can be obtained.

It is noteworthy that is greater than in Algorithm 1, which means that the number of elements in is less than that of . When is equal to 0, is filled with 0 to make the number of elements in equal to that of . Similarly, the number of elements in should be equal to that of When is equal to 0, is also filled with 0 to reach the same number of elements (i.e., ).