Abstract

This paper aims to provide a comprehensive scheme with limited feedback for downlink millimeter wave (mmWave) multiuser multiple-input multiple-output (MIMO) nonorthogonal multiple access (NOMA) system. Based on the feedback of the best beam and the channel quality information (CQI) on this beam, the users are grouped into a cluster having the same or coherent best beam and the maximal CQI-difference. To further reduce the intercluster interference, only the candidate cluster can join the cluster set whose intercluster correlation with the existing clusters is lower than threshold. Based on the results of clustering, mmWave hybrid beamforming is designed. To improve the user experience, each cluster selects the best beam of the user with the higher guaranteed rate requirement as the analog beamforming vector. For digital beamforming, the weak user applies the block diagonalization algorithm based on the strong user’s effective channel to reduce its intracluster interference. Finally, an intracluster power allocation algorithm is developed to maximize the power difference in each cluster which is beneficial to improve the successive interference cancelation (SIC) performance of the strong user. Finally, simulation results show that the proposed MIMO-NOMA scheme offers a higher sum rate than the traditional orthogonal multiple access (OMA) scheme under practical conditions.

1. Introduction

In contrast to the conventional orthogonal multiple access (OMA), nonorthogonal multiple access (NOMA) allows different users to efficiently share the same resources (i.e., time, frequency, space, and code) at different power levels. Exploiting their respective channel gain differences, multiuser signals are separated at receivers by successive interference cancelation (SIC). Thus, it is an important technique for 5G systems which can significantly enhance system capacity and overall spectral efficiency [13].

NOMA can flexibly combine other 5G technologies, e.g., MIMO, cognitive radio, cooperative communication, and channel coding [3]. Because MIMO-NOMA systems can significantly increase the number of supportable users and in turn improve the system sum capacity, the research of MIMO-NOMA has been carried out widely [47]. Multiple users with distinct channel gains are grouped into MIMO-NOMA clusters. The users in the same cluster are scheduled on NOMA basis. The intracluster interference is mitigated by SIC based on the obvious channel gain difference. Multiple clusters may utilize the multiuser MIMO principles to cancel the intercluster interference. References [48] prove that NOMA combined with MIMO techniques could achieve great system performance gains over OMA.

Thanks to the abundant bandwidth resources, millimeter wave (mmWave) frequency becomes a natural choice to achieve Gigabit data rates [9, 10]. The number of supportable users can be more than the number of radio frequency (RF) chains at the same resources by exploiting MIMO-NOMA in mmWave system. On the other hand, the channel of mmWave is highly directional. For mmWave MIMO-NOMA system, the users in a same cluster use the same beam and the different clusters use the distinct different beams. For the users in a same cluster, the obvious power difference in a same beam could assure good SIC performance. For different clusters, the distinct different beams are beneficial to greatly reduce the intercluster inference [11, 12]. Therefore, mmWave MIMO-NOMA can effectively improve the spectral efficiency and cope with the demands of massive connectivity.

Considering mmWave MIMO-NOMA system, [11] analyzes the achievable sum rate of the proposed beamspace MIMO-NOMA in a typical mmWave channel mode, which shows an obvious performance gain compared with the existing beamspace MIMO. Zero-forcing beamforming (ZF-BF) is designed to reduce the interbeam interference in the beamspace. Furthermore, a dynamic power allocation is proposed to maximize the achievable sum rate which contains the intrabeam and interbeam power optimization. Reference [12] proposes a hybrid analog/digital beamforming scheme with a power allocation algorithm aiming to maximize the energy efficiency. Users having a high channel correlation and large channel gain difference are grouped as a cluster. Reference [13] investigates the maximization of the sum rate of the system with clusters of 2-user or more-user. It decomposes the original problem into power and beam gain allocation problems. Reference [14] proposes a hybrid analog/digital beamforming framework to maximize the sum rate and analyzes the effect of beam misalignment on rate performance.

Most of the existing works on mmWave MIMO-NOMA, such as [1114], have assumed perfect knowledge of channel state information (CSI) at the transmitter, which is difficult to realize in practice. Reference [15] proposes a hybrid precoding scheme based on one-stage feedback, which makes use of all feedback overhead to enable precise beam information and takes advantage of the second-order channel statistics to mitigate multiuser interference. Reference [16] proposes the random beamforming in which base station (BS) broadcasts the random beams and the users with the enough signal strength feed their channel quality information (CQI) back to the BS.

MmWave MIMO-NOMA system with limited feedback has been studied but with limited results. Reference [15] has not yet considered user clustering and the intracluster interference cancelation with limited feedback. Reference [16] reduces the feedback overhead, but the channel measurement is time-consuming and complicated. In this paper, we put forward a comprehensive scheme with imperfect CSI for mmWave MIMO-NOMA system. To reduce the hardware complexity, we apply low RF chains structure at the BS, where the hybrid analog/digital beamforming is considered. The major contributions of this paper can be listed as follows:

(1) The user clustering is designed based on the limited feedback scheme. References [1621] put forward some feedback designs for MIMO-NOMA, e.g., on bit feedback, the path loss information, and other second-order statistics. Considering mmWave user could acquire the big gain in the specific beam, one-bit feedback or the path loss information is not accurate enough. In this paper, the limited feedback uses the best beam and CQI value based on this beam. The users having the same or adjacent best beam are sorted based on their CQI values, and the users with the maximal CQI-difference are formed a cluster. To further reduce the intercluster interference, the user clustering algorithm applies the intercluster correlation threshold to assure only the cluster with the distinct different beam can join the cluster set. Numerical results prove that it acquires obvious sum rate improvement by adequately exploiting the high directional feature of mmWave.

(2) We propose a novel hybrid analog/digital beamforming algorithm. In general, the users in a same cluster will use the beam of the strong user as the beamforming vector [11, 12]. The beamforming gains are assured to achieve first for the strong users. However, the weak users would require higher data rate. To improve the quality of user experience and the user fairness, for each cluster, we select the beam of the user with the higher rate requirement as the analog beamforming vector. For digital beamforming, the prior work [12] applies ZF-BF to cancel the intercluster interference. For mmWave communication, the distinct different beams cause very little interference as to the sparsity of mmWave channel. In this paper, the digital beamforming is designed for the strong user and the weak user, respectively. The strong user utilizes the singular value decomposition (SVD) of its effective channel to acquire more gains from digital beamforming. Especially, the weak user exploits the block diagonalization (BD) algorithm to mitigate the intracluster interference from the strong user in the same cluster.

(3) As the power difference between the two users in a same cluster increases, not only is the SIC performance of the strong user improved, but also the intracluster interference of the weak user is reduced. So we derive an intracluster power allocation scheme that maximizes the power difference under the users’ quality of service (QoS) requirements. We derive a closed-form optimal power allocation for the two users in a cluster by deriving the exact bounds for the power allocation coefficient region.

(4) Finally, we evaluate the sum rate performance of the downlink MIMO-NOMA system using the proposed beamforming, user clustering and power allocation algorithms. The simulations are done with a wide range of beam correlation thresholds, including the interuser and intercluster thresholds in the user clustering algorithm. Both of the perfect CSI and random selection cases are considered in the simulations. Numerical results also compare the sum rate performances of MIMO-OMA and the proposed MIMO-NOMA and illustrate the significance of the proposed MIMO-NOMA scheme.

The rest of this paper is organized as follows. Section 2 presents the system and channel model of mmWave MIMO-NOMA system. Sections 3, 4, and 5 describe the designs of user clustering, hybrid beamforming and power allocation, respectively. The simulation configuration and results are presented and discussed in Section 6, and finally, Section 7 concludes this paper.

Notation. In the reminder of this paper, we use the following notations: and denote the transpose and Hermitian transpose, respectively, denotes the Frobenius norm, and denotes the distribution of circularly symmetric complex Gaussian random variable with mean and covariance .

2. System and Channel Model

Consider a downlink mmWave beamforming NOMA transmission system, in which one BS communicates with clusters. Each cluster consists of multiple single-antenna users. The BS is equipped with antennas and RF chains for reducing the hardware complexity. We consider the fully connected architecture that each RF chain is linked to all antennas by using multiple phase shifters. We assume that the number of clusters served by the BS is equal to that of RF chains . Although each cluster can contain more than two users, we suppose here that there are two users in each cluster for simplicity. This is consistent with the standard implementation of NOMA in long term evolution advanced (LTE-A) as well [22].

On the downlink, the BS applies an baseband beamforming followed by an RF beamforming . The sampled transmitted signal is thereforewhere , is the vector of transmitted symbols in which and denote the signals for the 1-th and 2-th in the -th cluster, respectively. Similarly, and stand for the power allocation coefficients for the two users in the -th cluster. is the power allocated to the -th cluster.

The received signal of the -th user in the -th cluster can be represented as follows:where is the vector that represents the mmWave channel of the -th user in the -th cluster. denotes additive white complex Gaussian noise. Moreover, the second term is the intercluster interference from other clusters. Without loss of generality, for each cluster, we assume that where . Following this, the 1-th and 2-th are defined as the strong user and weak user in each cluster, respectively.

Due to the limited scattering in mmWave channel, we adopt a well-established geometric channel model with scatterers [23, 24]. We assume that each scatterer contributes a single propagation path from the BS to user. Then, the channel can be expressed as where the first and second terms denote the line-of-sight (LoS) component and non-line-of-sight (NLoS) components of the -th user in the -th cluster, respectively. denotes the distance from the BS to the -th user in the -th cluster. and are the path loss exponents corresponding to the LoS and NLoS paths, respectively. where denotes the complex gain of the -th path between the BS and the i-th user in the l-th cluster. is the -th path’s normalized direction, whereas is the antenna array steering vector with respect to . If we assume a uniform linear array (ULA) is used at the BS, can be written as [25]where is related to the angle of departure as [26]. denotes the distance between antenna elements and denotes the signal wavelength satisfying at mmWave frequencies.

3. User Clustering

Appropriately selecting two users which are served in a cluster can help improve the performance of NOMA multiuser beamforming. On the one hand, the big difference of the channel gains can improve SIC performance of the users with high channel gains and reduce the intracluster interference of the users with low channel gains. On the other hand, the beamforming vector is shared by all users in the same cluster. If the channels of those users in a same cluster are highly correlated, the beamforming vectors can acquire the definite array gains and properly cancel the intercluster interference from other clusters. Reference [4] proposed a clustering algorithm selecting user-pair having a high channel correlation and large channel gain difference as set in the same cluster, which requires full CSI. However, the full CSI feedback to the BS is not feasible in practice due to prohibitively high feedback overhead. So we propose a user clustering algorithm with limited feedback information, which only requires the indexes of the best beams and the channel quality information (CQI) values.

The index of the best beam and the user’s CQI are measured and sent back as follows. In the first stage, the BS broadcasts pilot signals on the every vector in the beamforming codebook. In the second stage, each user chooses the best beam and measures the CQI with this beamforming vector. At last, the user sends the index of the best beam and CQI to the BS. The best beam of the i-th user can be selected according to the following criterion as [27]where represents the beamforming codebook with size (i.e., feedback bits) and consists of the steering vectors . They have the same form of the array response vector in (4). is the beamforming vector selected from .

As discussed in [24, 28], the effect of LoS component is dominant in mmWave channels, compared to those of NLoS components. Therefore, the above mmWave channel model in (3) can be simplified as where the subscript 0 of the LoS component is omitted for the simplification of notation. Then, the effective channel gain of the user on the beam in (5) can be expressed aswhere is the predefined normalized direction of . is the Fejér kernel, which goes to zero quickly when increases. Therefore, based on the criterion in (5), the normalized direction of the selected beam has the least difference from in the codebook. This means that the normalized direction of each user can be approximately the normalized direction of its selected beam. The correlation of the selected beam between user and user can be rewritten asBased on the above beam correlation expression in (8) and channel expression in (6), the channel correlation between users and user can be simplified aswhere denotes the channel correlation between user and user . (a) follows from that and are approximately and based on the above analysis, respectively. So the correlation between the users’ channels can be replaced with the correlation between their selected beams .

CQI is an indicator carrying information about the channel condition. In LTE, considering a single BS system, the received SNR for each user can be written as [29]where is the transmit power and is the power of received noise and interference. The measured SNR is then mapped to a discrete CQI value using a chipset vendor specific mapping table [30]. Because the positive correlation between channel gain and CQI value, we will just calculate the difference of CQI values between users instead of their channel gain difference.

It is proved in [8] that the more users are admitted to a cluster, the lower is the achieved sum rate, which illustrates the tradeoff between the sum rate and maximum number of admitted users. Considering the tradeoff, we assume that each cluster admits 2 users. Based on the above analysis, a limited feedback downlink MIMO-NOMA user clustering algorithm is described in Algorithm 1, in which the number of clusters and each cluster contains 2 users.

Input and Initialization:
Input: the number of clusters , feedback beam and CQI value , for .
Initialize: .
where is the candidate set of user-pair that can be selected as a cluster. is the number of
selected clusters. A is the set of the strong users of the selected clusters. is the set of the
selected users in the l-cluster. is the set of weak users that have been selected.
Output:the user set selected for each cluster
Step (1) Select all user-pairs having a beam correlation higher than the threshold:
The BS sorts users according to the ascending CQI value: .
for = 1  :  K-1
for j = K  :  -1: i + 1
if   AND
break
end if
end for
end for
where is the predefined beam correlation threshold .
Step (2) Pair selection for clustering:
while  l < L
(a) Select the user-pair that has the maximum CQI value difference in :
(b) Put the selected user-pair having correlations with the users in A below the threshold into a new cluster:
if  
remove the pairs in that include user or
else
end if
end while
where is the predefined inter-cluster beam correlation threshold .

Algorithm 1 can be summarized as follows. In the first step, the BS calculates the beam correlation of all possible user-pairs and puts the user-pair having a beam correlation higher than the threshold into set . In the second step, during each loop, the BS searches for a new user-pair having the maximum CQI value difference from T. Then the BS calculates the beam correlations between the user and the selected strong users in set A. If all the correlations are lower than the threshold , then the user-pair is added to a new cluster and the user is put into A. Otherwise, no users are added. The loop is repeated consecutively until L user clusters are selected.

4. Two-Stage Hybrid Beamforming

4.1. The Optimization of Analog Beamforming

For analog beamforming, one analog beamforming vector is designed for one user in traditional design [31], while two users in a same cluster should share one analog beamforming vector in MIMO-NOMA. Furthermore, considering the fairness of QoS, the rate requirement of the weak user may be greater than that of the strong user in a cluster. So we select the feedback beam of the user with the maximum guaranteed rate requirement as the analog beamforming vector in each cluster, i.e., , where is the guaranteed rate requirement of the i-user in the l-cluster.

Based on the above analysis, we can design the analog beamforming vector for the -th cluster according to the following criterion:Then, the analog beamforming matrix can be obtained as .

4.2. The Design of Digital Beamforming

In this subsection we design the digital beamforming. After the analog beamforming, the BS trains the effective channels with users asThe dimension of the effective channel is which is much less than that of its original channel . This is different from the algorithms proposed in [32, 33] in which the effective channel has a larger dimension. Then, each user quantizes its effective channel utilizing a random vector quantization (RVQ) codebook of size and feeds the index of the vector back to the BS with bits, where [34]. Finally, the BS designs its digital beamforming based on the quantized channels.

The user clustering and analog beamforming mitigate the intercluster interference. Then we can design the digital beamforming vector for the weak user to reduce its intracluster interference. We take the design of the digital beamforming vector for the weak user in the -th cluster as an example. First, we construct a complementary matrix that contains the other user in the same cluster, i.e., . Note that and . It is obviously . Then, we define the SVD of aswhere holds the last right singular vectors and forms an orthogonal basis for the null space of . We perform SVD on asFinally, we can design the digital beamforming vector aswhere denotes the first column of . For the strong user, we can design a digital beamforming vector to increase its capacity and in turn reduce its interuser interference. More exactly, we define the SVD of asThen we can design the digital beamforming vector for the strong user as

5. Power Allocation

In the proposed system, all the clusters contain two users (same cluster size), then the total transmit powers are equally allocated to the clusters, i.e., . This approach can be said as nearly optimal since each cluster contains users with nearly similar channel gain distinctness [6].

With the increasing of the power difference between two users in a same cluster, the SIC performance of the strong can be better, and the intracluster interference of the weak user can be also reduced because the power allocated to the signal of the strong user is decreased. So we propose a power allocation scheme for intracluster power allocation which maximizes the power difference to help perform SIC and decrease the intracluster interference for the weak user, while keeping the data rates of the two users are greater than or equal to their own minimum rate requirements. The formulation is as follows:where is the power coefficient difference between the two users. Here, and are defined as the strong and the weak users’ rates in the l-cluster if the two users would be supported by conventional TDMA hybrid beamforming system, respectively, i.e., . In conventional hybrid beamforming system with L RF chains, users need two time slots to be supported, while it only takes one time slot to support users in MIMO-NOMA system. So is multiplied by 1/2. is defined aswhere is the analog beamforming matrix containing all the feedback beams of the i-user in each cluster. is designed for the i-user in the l-th cluster by using ZF-BF based on the effective channels of the i-user in each cluster presented in [34].

For the strong user, the power allocation coefficient that ensures its rate is found by solving (19). Solving (19) for giveswhere and . The least value of power allocation coefficient to ensure that NOMA is fair to the strong user is given by the right side of (22), and any satisfying (22) will lead to (19).

Similarly, if the rate of the weak user is to be at least as good as conventional hybrid BF, then (20) leads towhere and . Therefore, the greatest value of the power allocation coefficient to ensure that NOMA is fair to the weak user is given by the right side of (23), and any satisfying (23) will lead to (20).

If is satisfied, the NOMA power allocation region is therefore defined as . Since the objective function is monotonically decreasing function of , the optimal power allocation coefficient for the strong user in the l-cluster, , can be given as In the case where , (19) and (20) cannot be satisfied simultaneously. We will set to allocate more power to the strong user to increase the system sum rate while keeping the rate of the weak user is equal to its minimum rate requirement, because the strong user has a larger channel gain and receives no intracluster interference by SIC.

6. Simulation Results

This section will present the computer simulation results to evaluate the performance of the proposed MIMO-NOMA scheme. We consider a multiuser downlink mmWave MIMO system with a BS employing a ULA of transmitting antennas and RF chains, i.e., 3 clusters in the MIMO-NOMA system. The mmWave bandwidth of the system is assumed to be 2 GHz, and its carrier frequency is assumed to be 60 GHz. The number of users . And the all users are randomly distributed in a cell with a radius of 20m. We assume that there is a LoS path and 2 NLoS paths for all users’ channels. We set the channel parameters of each user as follows: (1) , ; (2) the normalized direction for follows the uniform distribution over . The beamforming codebook and RVQ codebook are quantized with = 5 bits and = 13 bits, respectively.

Figure 1 shows the sum rate performances of the proposed MIMO-NOMA system and the conventional TDMA hybrid beamforming system mentioned in the Section 5. The key observation from Figure 1 is that the proposed MIMO-NOMA system outperforms conventional MIMO-OMA system when the beam correlation and intercluster correlation thresholds are set as and , respectively. Moreover, by comparing the three cases, i.e., = 0.7, 0.3, and 0.05, Figure 1 shows an improvement of sum rate of the MIMO-NOMA system as increases. This is because, with the increase of the beam correlation threshold , more highly correlated users are clustered through the proposed clustering algorithm. The MIMO-NOMA system with a higher channel correlation between the users can help reduce the intercluster interference for the users by user clustering and analog beamforming.

Figure 2 shows the sum rate performances of the proposed MIMO-NOMA system, where . Observe all the curves; we can see that the sum rate of the proposed system increases as the intercluster correlation threshold decreases. These results are expected because when is reduced, clusters having lower intercluster correlations can be selected per the proposed clustering algorithm. A lower correlation between the selected clusters can help reduce the intercluster interference for the users.

Figure 3 shows the sum rate performances in three cases: randomly selecting, perfect CSI, and the proposed scheme. Randomly selecting means selecting two users from all the users as a cluster randomly. Perfect CSI means that perfect CSI is assumed to be available at BS to generate beamforming matrices and cluster users with their channel correlation and channel gains instead of beam correlation and CQI values by Algorithm 1. Randomly selecting is the simplest of the three cases, but results in a degraded sum rate performance because this scheme randomly clusters users ignoring their intercluster interferences. On the other hand, the results show that the sum rate of the proposed scheme is lower than that of perfect CSI but requires less feedback overhead.

Next, we examine the advantage of the proposed hybrid beamforming for the MIMO-NOMA system with the conventional NOMA-ZF, which is performed by utilizing ZF-BF based on the effective channel information of the strong user of each cluster [12], cf. Figure 4. At very low SNR, the two beamforming schemes asymptotically approach the same sum rate. As SNR increases, the proposed beamforming scheme achieves a larger performance gain. This is because the proposed beamforming scheme reduces the interference of the weak users and the conventional NOMA-ZF does not.

Finally, we compare the proposed mmWave massive MIMO-NOMA scheme with limited feedback with the random beamforming scheme with one-bit feedback in [16]. The central angle of the sector is set as 2Δ = 0.2, and the power allocation coefficients of the two users in each cluster are set as 1/4 and 3/4. The number of the sectors in the random beamforming scheme is set as 3, which is the same as the number of RF chains in the proposed scheme. As can be observed from Figure 5, the performance of the proposed scheme degrades with decreasing the number of quantization bit . The figure also shows that the feedback overhead of the proposed scheme is larger than that of the random beamforming scheme with one-bit feedback, but it achieves better performance.

7. Conclusions

In this paper, we propose a novel mmWave MIMO-NOMA system with limited feedback. A low-complexity user clustering algorithm and hybrid beamforming scheme have been designed with the limited feedback CSI. On this basis, we propose a power allocation solution to maximize the power difference between the two users in each cluster under the rate constraints of the two users. Simulation results show that the proposed MIMO-NOMA scheme achieves a better sum rate performance than conventional MIMO-OMA, and the inter- and intracluster interferences and feedback overhead are reduced by applying the proposed clustering, beamforming, and effective power allocation algorithms. Therefore, the proposed scheme represents a feasible option to ensure that the mmWave MIMO-NOMA system provides an improved sum rate performance with limited feedback.

Data Availability

All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This paper was supported by National Natural Science Foundation of China under Grant 61871321, National Science and Technology Major Project under Grant 2016ZX03001016, and Innovation Team Project of Shaanxi Province under Grant 2017KCT-30-02.