Joint User Scheduling and MU-MIMO Hybrid Beamforming Algorithm for mmWave FDMA Massive MIMO System

,


Introduction
Millimeter wave (mmWave) communication will unleash the 30-300 GHz new spectrum and enabled gigabit-per-second data rates for the next generation mobile cellular systems [1][2][3].Massive MIMO is an essential part of mmWave communication to combat the stringent constraints imposed by the high propagation loss.A beneficial feature of millimeter wave is that large-scale antenna arrays can be packed into small dimensions thanks to the very small wavelength [4,5].However, the digital processing in traditional MIMO system requires dedicated baseband and RF hardware for every antenna element.With the large-scale antenna, the high cost and power consumption of mmWave RF hardware preclude such a transceiver architecture at present.In mmWave massive MIMO system, the trade-off between performance and simplicity drives the need to deploy beamforming at both the digital and analog domains, that is, hybrid beamforming.
Most current hybrid beamforming algorithms [6][7][8] assume user scheduling based on STDMA that a single user or multiple spatial multiplexing users are scheduled in the entire frequency band at a given time slot.Actually, the large bandwidth and multipath nature of mmWave channels in a cellular system assure the existence of frequency selective channels; it is necessary that mmWave communication remains with frequency division multiple access (FDMA) and user scheduling.OFDMA provides a natural multiple access method by assigning different users with orthogonal subcarriers, and multiuser diversity gain in frequency domain can be exploited by multiuser subcarrier scheduling [9,10].But for the hybrid beamforming system, the analog beamforming is implemented by the phase shifts which are constant in the entire frequency band.There are multiple frequency multiplexing users in the entire frequency band; they will experience the same analog beamforming processing.Thus it needs a joint optimization of user scheduling and the wideband processing of the analog beamforming.On the other hand, mmWave links are inherently directional, and the antenna array steers its beam towards any direction electronically and to achieve a high gain at this direction, 2 International Journal of Antennas and Propagation while offering a very low gain in all other directions.It is beneficial for MU-MIMO because the RF beams have sufficient degrees of freedom to be optimized for MU-MIMO.In this paper, we will design the joint user scheduling and MU-MIMO hybrid beamforming scheme for mmWave massive MIMO-OFDMA system.
References [11,12] have meaningful researches on OFDMA scheduling for the hybrid beamforming system.Reference [11] computes the analog beamforming matrix as the first   eigenvectors of the left singular value decomposition (SVD) of the combined digital precoding matrices of subcarriers having the highest sum rate (  is the amount of RF chains).Then, for fixed A, the digital beamforming matrix is computed and its corresponding users are scheduled such that the total sum rate is maximized for each subcarrier.Reference [12] enables users with high cochannel interference to be scheduled in different frequency channels in the same time slot while sharing the same RF chain and analog beam.
For MU-MIMO hybrid beamforming algorithms, [13] considers the zero forcing (ZF) hybrid beamforming which essentially applies phase-only control at the RF domain and then performs a low-dimensional baseband ZF precoding based on the effective channel seen from baseband.Reference [14] designs the hybrid beamforming by considering a weighed sum mean-square error (WSMSE) minimization problem incorporating the solution of the detected signals which is obtained from the block diagonalization technique.The resulting WSMSE problem is solved by applying the orthogonal matching pursuit algorithm.Reference [15] analyzes a low complexity hybrid precoding algorithm for downlink multiuser mmWave system, which configures the hybrid precoder at the transmitter and analog combiners at multiple receivers with a small training and feedback overhead.For this algorithm, a lower bound on the achievable rate for the case of single-path channels is derived.Reference [16] designs a hybrid block diagonalization scheme to approach the capacity performance of the traditional BD processing method, aiming to harvest the large array gain through the phase-only RF precoding and combining.
In [11], the analog beamforming matrix is acquired from SVD or QR decomposition, so it is hard to implement with the phase shift of the traditional analog beamforming.Reference [12] considers that the base station applies only one RF chain to transmit signals to users scheduled in different frequency channels.Actually there will be several RF chains in mmWave communication; the scheduled user in multiple subbands could transmit signals from different RF chains.Considering the contradiction between the optimization of the wideband analog beamforming and multiple users scheduled in different frequency subbands, we propose a joint user scheduling and MU-MIMO hybrid beamforming scheme for mmWave OFDMA system.The contribution of this paper can be summarized as follows: (i) User scheduling algorithm achieves frequency resource allocation and MU-MIMO user selection.Firstly the users with identical optimal beams are defined as an OFDMA user group.For an OFDMA user group, all members multiplex the entire frequency band and the one with the best channel gain is assigned in the corresponding frequency resources.
Then the frequency domain channel of each OFDMA user group is defined as an integrated channel, which is regarded as a virtual user in MU-MIMO user selection and analog beamforming design.Finally MU-MIMO users are selected to maximize the mmWave system throughput.
(ii) An OFDMA user group can be regarded as a virtual user that the RF beam of every member is the same; thus the hybrid beamforming could not only coordinate the contradiction of the wideband analog beamforming and multiple users scheduled in frequency resource, but also support arbitrary RF number and MU-MIMO algorithms.
(iii) For the proposed hybrid beamforming, the analog beamforming vectors adopt the optimal beam of each scheduled user, since the performance of each user is sensitive to beam direction.The digital beamforming algorithm is solved by weight MMSE, which not only achieves the optimal performance for the single user, but also mitigates the residual interuser interference.
(iv) Evaluate the sum rate of the existing and the proposed MU-MIMO hybrid beamforming under different number of BS antenna and scheduling users.In the simulation, the proposed user scheduling algorithm in this paper is based on FDMA, whereas other reference algorithms are based on TDMA.Simulation results show that hybrid beamforming together with user scheduling can greatly improve the performance of mmWave OFDMA massive MU-MIMO system.
The rest of this paper is organized as follows.Sections 2 and 3 introduce the system model and channel model.In Section 4, the proposed user scheduling algorithm is provided.In Section 5, the proposed hybrid beamforming algorithm is presented, and computer simulation results are shown in Section 6.Finally, conclusions are drawn in Section 7.
Notations.In this paper, upper-case/lower-case boldface letters denote matrices/column vectors.X  , X  , X −1 , and X(, ) denote transpose, conjugate transpose, inversion, and the (, )th element of X, respectively.[X]  denotes the th diagonal element of X. E(X) denotes the expected value of X. tr(X) denotes the trace of X.We define ‖X‖ 2  as tr(X  X) and ‖X‖ 2 as the square root of the maximum eigenvalue of X  X.I is an identity matrix with appropriate size; C × represent spaces of × matrices with complex entries.The acronyms s.t and i.i.d denote "subject to" and "independent and identically distributed," respectively.

System Model
In this paper, we will consider downlink OFDMA MU-MIMO mmWave system as shown in Figure 1  MU-MIMO spatial multiplexing users, which are selected out from all  all users; every user has   antennas.To enable multistream communication, the transmitter is equipped with  RF transmit chains such that   ≤  RF ≤   .We employ resource block based OFDMA transmission where each block occupies   adjacent subcarriers and   consecutive OFDM symbols.At each time slot of duration , the transmitter broadcasts   blocks and   is the amount of resource blocks [17,18].
The transmitted data streams at the BS are assumed to be processed by a digital beamformer W in the baseband, followed by an analog beamformer F before transmission.Notably, F can realize only phase changes (phase-only control), since it is implemented using analog phase shifters; each entry of F is constrained to satisfying ‖F , ‖ = 1/√  .Furthermore, the total power constraint is enforced by letting ‖FW‖ 2  =   .For simplicity, we will describe the system model in frequency domain.Analog beamforming vector f in time domain will be transformed to F in frequency domain just by employing fast Fourier transform (FFT) operation.At the MS, a digital combiner B is used to process the received signal.
Assume the transmitted signals of different users are independent from each other and from noises, the received signal Y , ∈ C   ×1 of the th user scheduled in the th block can be written as where s , ∈ C  ,  ×1 denotes the transmitted signal vector for th user in the th block, satisfying , n , ] = 0, and  ,  represents the number of the th user's data streams scheduled in the th block.W , ∈ C  RF × ,  denotes the transmitting digital beamformer for the th user in the th block, and F ∈ C   × RF denotes the transmitting analog beamformer in the frequency domain which is FFT transformed by time domain

Channel Model
Since mmWave channels are expected to have limited scattering, we adopt a geometric channel model with   rays for the channel of user  in block .Under this model, the channel H , can be expressed as [12,13] where   is the channel impulse response of the th path with E[|  | 2 ] = 1 and   is the path loss between BS and users.Considering the azimuth and elevation angles, the vectors a  (   ,    ) and a  (   ,    ) represent the normalized transmit and receive array response vectors at an azimuth (elevation) angle of departure    (   ) and that of arrival    (   ), respectively.For a uniform planar array (UPA) in the -plane with  row and  col elements on the and -axes, respectively, the array response vector at the BS is given by [19] = 2  cos  sin /, ] = 2  sin  sin /, where  is the wavelength and   and   are the distances between two adjacent antenna elements in the and -axes, respectively.vec(⋅) denotes the matrix of dimension  row ×  col stretches to the column vector of dimension   × 1.The antenna array at the user side utilizes a uniform linear array (ULA); thus array response vectors at the user side are given by where a  (   ) is the column vector of dimension   × 1.

User Scheduling Scheme
This section discusses the proposed user scheduling algorithm and MU-MIMO user selection to maximize the total sum rate of all subcarriers, which can be treated as a stage prior to the hybrid beamforming design.This section is divided into three subsections.In the first subsection, we define OFDMA user groups that the users in a same group have identical RF beam and multiplex the whole frequency resource.In the second subsection, we allocate the frequency resource for each member and define each OFDMA group as a virtual user.In the last subsection, we select several MU-MIMO users which can maximize the mmWave system throughput.
4.1.OFDMA Group Selection.Firstly, users transmit uplink sounding signals, and BS uses the sounding results to select the strongest beam for every downlink user.According to the indices of the selected beams, the users which have identical selected beam form an OFDMA user group.OFDMA group selection problem is mathematically formulated as follows: where where SNR , is the signal noise ratio for user  in the th block.
After the frequency resource scheduling, the users in this group are sorted as follows: For the th OFDAM user group Ω  OFDMA , the frequency channels of all members are merged into an integrated channel which represents the spatial characters of the users multiplexed in different frequency resources.The integrated channel of OFDAM user group Ω  OFDMA is defined as where H   represents the channel of the user scheduled in the th resource block for the th OFDMA user group and so on.

MU-MIMO User
Selection.Every integrated channel associated with an OFDMA user group is regarded as a virtual user.The virtual users and other users to be scheduled constitute the candidate user set Ω  .Specially, MU-MIMO channels of all spatial multiplexing users are defined as H 1 , . . ., H  .For OFDMA user group, H  is the integrated channel associated with this OFDMA user group.For other users which independently occupy the entire frequency resource, H  is the channel of the th MU-MIMO user.
This section discusses the proposed MU-MIMO user selection to maximize the total sum rate.The detailed solution can be given by two steps.
Firstly, to decrease the complexity of searching MU-MIMO user, the user with best channel gain is given a high priority in MU-MIMO user selection.Thus, the user with the maximum SNR is selected as the leader of MU-MIMO user selection.The leader is formulated as where H is the equivalent channel considering beamforming gains.As is known, an upper bound on the performance of hybrid beamforming approximates to that of full digital beamforming for any design criteria [14], and SVD is a typical beamforming method in numerous beamforming algorithms.Thus we utilize SVD as full digital beamforming weighting vectors when calculating SNR of the target user [20,21].The MIMO channel of the th user H  can be decomposed by SVD as H  = U  Λ  V   and the equivalent channel H is defined as H = U   H  V  .At the second step, the objective function of the other MU-MIMO user selection is defined as follows: where   sum is the sum rate when user  is scheduled as expressed in (11).
The MU-MIMO user selection process continues until iterating  − 1 times with employing exhaustive search and composes the selected user set Ω  , Ω  = {Ω 1  , . . ., Ω   },  ≤  RF .Note that the maximum number of multiplexed multiusers is equal to the number of RF links.
Outline of user scheduling algorithm is described as shown in Algorithm 1.
Initialization: OFDMA group set Ω OFDMA = Φ, the candidate user set Ω  = Φ, the selected user set Ω  = Φ Step 1. Select the users with same optimal beams to form OFDMA group set Ω OFDMA .
User Algorithm 1: User scheduling scheme.

Hybrid Beamforming Designs in Massive MIMO Systems
In this section, we design a hybrid beamforming algorithm for mmWave OFDMA massive MIMO system as illustrated in Figure 1.For the hybrid beamforming system, the analog beamforming is implemented by the phase shifts which are constant in the entire frequency domain.Considering the contradiction between the optimization of wideband analog beamforming and multiple users scheduled in OFDMA system, an effective solution is that every user or virtual user of MU-MIMO selected user set maps its own transmitting signals to a unique RF chain.That is, the analog beamforming vector for a RF chain will be optimized based on the channel character of only one user or the virtual user with same optimal beam.Because the integrated channel can be treated as a virtual user with similar optimal beam in the entire frequency band, the hybrid beamforming could not only reconcile the contradiction of the wideband analog beamforming and multiple users scheduled in frequency resource, but also support arbitrary user scheduling band and MU-MIMO algorithms.
In what follows, we split the proposed hybrid beamforming design into two steps: Firstly, for analog beamforming design, we adopt the optimal beam of each scheduled user.Then, we focus on digital beamforming design.The digital beamforming algorithm not only achieves the optimal beamforming gains for every user allocated in its own block, but also mitigates the interuser interference in the same block.Further, the digital beamforming is formulated by the weighted MMSE.

Analog Beamformer Design. Firstly, MU-MIMO channel matrix H in frequency domain is transformed to time domain h just by employing inverse fast Fourier transform (IFFT);
Consider the analog beamformer design, the achievable rate is  sum = log 2 (1+‖h  f  ‖ 2 / 2 ) for th MU-MIMO user and we seek to design the analog beamformer f  to maximize sum rate by scanning a codebook F  [22], which can be expressed as max Since every RF chain corresponds to one MU-MIMO user, then f = [f 1 , f 2 , . . ., f  RF ] is the selected time domain beam vectors from the predefined RF beamforming codebook F  .F  is specified in a quantized matrix a  (2  /  , 2  /  ), where each column is a weight vector corresponding to one beam pattern, for the variable   taking the values 0, 1, 2, . . .,   − 1 and   taking the values 0, 1, 2, . . .,   − 1, and   and   denote the quantized precision of azimuth and elevation angles, respectively.Note that the objective function of MU-MIMO user selection is similar to that of analog beamforming.The difference is that only MU-MIMO user selection considers the interuser interference.An enormous amount of simulation results indicates that the analog beamforming performance obviously decreases while the optimal beam direction of the target user is slightly changed.Thus for the analog beamforming, we still apply the optimal beam for each selected user.The interuser interference can be further mitigated by the digital beamforming.

Digital Beamformer Design.
In this stage, we design the digital beamformer by the weighted MMSE approach [23,24] to mitigate multiuser mutual interference.Assume g , ≥ 0 be a weight matrix for user , and the weighted sum-MSE minimization used to deal with the problem is formulated as where  denotes the power budget, and the mean-square estimation error matrix  , can be written as For fixed all W , , MMSE receive beamforming B , at user  is given as Let Then, the corresponding MSE error matrix for user  applying the receive beamforming B mmse , can be written as Because the object function of ( 13) is convex in each of the optimization variables g, W, B, the block coordinate descent method is adopted to solve (13).Specifically, the weighted sum-MSE object function is minimized by sequentially fixing two of the three variables g, W, B and updating the third variable.While the update of receiver beamforming B , is expressed by (16), the update of the weight variable g , is in closed form that can be written as The update of transmit digital beamforming W , can also be decoupled through transmitters, causing the following optimization problem: We can exploit standard convex optimization approaches to solve this convex quadratic optimization problem.Meanwhile, we can also apply the Lagrange multipliers method to get a closed form solution.In particular, assuming a which is equivalent to Notably the optimum   (denoted by  *  ) must be positive in this case and the left-hand side of ( 24) is a decreasing function in   for   > 0. Therefore ( 24) can be easily worked out by employing one-dimensional search techniques.Eventually, by plugging  *  in (24), we can achieve the solution for W , ( *  ),  = 1, 2, . . ., ,  = 1, 2, . . .,   .The digital beamforming algorithm for the mmWave massive MIMO is summarized in Algorithm 2.
The optimization problem (13) has a differentiable objective function and a constraint set that is separable in the variables g, W, B. The WMMSE algorithm is the block coordinate descent method applied to (13) and converges to a stationary point of (13), which is any limit point g * , W * , B * of the iterates generated by the WMMSE algorithm.

Simulation Results
In this section we present simulation results to characterize the performance of the proposed algorithm presented in Sections 4 and 5.The simulated channel is mmWave MIMO multiantennas channel model [25] extended from IEEE 802.11ad channel model [26].We consider a single-cell MIMO-OFDMA system consisting of BS and  all users, and the cell radius at the BS is the typical value for a microcellular system.The propagation environment is modeled as a   = 8 ray channel with uniformly random azimuth and elevation AoAs/AoDs distributed in [0, 2] and [0, ].We assume that the Channel State Information (CSI) is updated once per frame and group 5 OFDM symbols into a frame.One OFDM symbol has 512 subcarriers including 352 data subcarriers.The transmitter is assumed to employ a UPA antenna and the UE is assumed to have   = 2 antennas.The interelement spacing in both BS and UE antenna arrays is set to half a wavelength.The other simulation parameters are shown in Table 1.First, we compare the sum rates of different algorithms when  = 2 and   = 512 in Figure 2, where we consider the system model in Section 2. The rate achieved by the proposed hybrid beamforming algorithm is compared with the full digital beamforming, analog beamforming, and existing OMP + BD hybrid beamforming algorithm [14,27].To evaluate the proposed hybrid beamforming algorithm, we assume azimuth and elevation angles of the phase shifters at the BS to be quantized with   = 20 and   = 20.As can be seen from Figure 2, as expected the rate achieved by the proposed hybrid beamforming is higher than that of the analog beamforming, and superior performance is achieved by the digital beamforming approach.Meanwhile, the performance of hybrid beamforming is very close to the digital beamforming from the low to moderate SNR regions, and small performance gap is observed at SNR regions within −22∼−10 dB.When SNR increases to more than −10 dB, the proposed hybrid beamforming algorithm will achieve the same saturation value as digital beamforming.And the proposed hybrid beamforming yields a large improvement over the existing OMP + BD hybrid beamforming.

International Journal of Antennas and Propagation
Secondly, we examine and compare the performance of the proposed hybrid beamforming with different number of scheduled users per block.We set the same number of transmit antennas in Figure 2 as   = 512 and vary maximum number of scheduled users in different blocks as  = 2, 4. As can be seen from Figure 3, increasing the number of MU-MIMO users multiplexing in each block can improve the sum rate (for SNR values) due to canceling the residual multiuser interference.The sum rate with  = 4 has risen sharply at SNR = −30∼−14 dB.When the value of SNR approaches to −12 dB, it approximately tends to be maximum value and then starts to level off as SNR exceeds −12 dB, while the sum rate with  = 2 will also achieve the same maximum value as the value of SNR exceeds −10 dB.
Finally, we discuss the performance of hybrid beamforming with different number of transmit antennas.The sum rate is achieved by the hybrid beamforming with   = 128, 256, 512 antennas with  = 4.The SNR ranges from −30 dB to −10 dB.It can be found in Figure 4 that as the number of transmit antennas increases, the performance of system enhances significantly from the low to moderate SNR regions.When the SNR approaches to −12 dB, the sum rate with   = 512 approximately tends to be stable.And it saturates after a certain SNR which is around −12 dB.As the value of SNR continues to increase, the sum rate of other cases will also achieve the same maximum value.

Conclusions
In this paper, we propose a joint user scheduling and MU-MIMO hybrid beamforming algorithm for mmWave FDMA massive MU-MIMO system.The users with the same strongest beams direction form an OFDMA user group.For the same OFDMA user group, BS allocates the corresponding frequency resources to the member with the best channel gain.Then, to maximize the mmWave system throughput, MU-MIMO users are selected from OFDMA user groups and other users to be scheduled.For the proposed hybrid beamforming, the analog beamforming vectors apply the optimal beam of each MU-MIMO user.The digital beamforming algorithm is solved by weight MMSE, which not only achieves the optimal performance for each user, but also mitigates the residual interuser interference.Simulation results show that the performance of hybrid beamforming is very close to the full digital beamforming.Thus hybrid beamforming jointly designing with user scheduling can greatly improve the performance of mmWave OFDMA massive MU-MIMO system.Our further work will focus on reducing the complexity of the proposed hybrid beamforming algorithm and increasing the function of adaptive power allocation.

Figure 1 :
Figure 1: Block diagram of mmWave MU-MIMO system based on frequency domain user scheduling.
. The BS with   transmit antennas communicates   data streams to , ∈ C   ×  denotes the MIMO channel of the th user.n , is the vector of i.i.d additive complex Gaussian noise with zeromean and variance  2 .
,  denotes linear receive beamforming vectors to detect the transmit signals.H , User  ∈ Ω  Step 3. The candidate user set Ω  consists of virtual users and other users to be scheduled.Select the user with the maximum SNR as the leader from Ω  .